Single Spin Asymmetry in Charmonium Production

We present estimates of single spin asymmetry (SSA) in the electroproduction of taking into account the transverse momentum dependent (TMD) evolution of the gluon Sivers function and using Color Evaporation Model of charmonium production. We estimate SSA for JLab, HERMES, COMPASS and eRHIC energies using recent parameters for the quark Sivers functions which are fitted using an evolution kernel in which the perturbative part is resummed up to next-to-leading logarithms accuracy. We find that these SSAs are much smaller as compared to our first estimates obtained using DGLAP evolution but are comparable to our estimates obtained using TMD evolution where we had used approximate analytical solution of the TMD evolution equation for the purpose.

QUANTITATIVE OH PLANAR LASER INDUCED FLUORESCENCE DIAGNOSTICS OF SYNGAS AND METHANE COMBUSTION IN A CAVITY COMBUSTOR

The current work reports quantitative OH species concentration in the cavity of a trapped vortex combustor (TVC) in the context of mixing and flame stabilization studies using both syngas and methane fuels. Planar laser induced fluorescence (PLIF) measurements of OH radical obtained using a Nd: YAG pumped dye laser are quantified using a flat flame McKenna burner. The momentum flux ratio (MFR), defined as the ratio of the cavity fuel jet momentum to that of the guide vane air stream, is observed to be a key governing parameter. At high MFRs similar to 4.5, the flame front is observed to form at the interface of the fuel jet and the air jet stream. This is substantiated by velocity vector field measurements. For syngas, as the MFR is lowered to similar to 0.3, the fuel-air mixing increases and a flame front is formed at the bottom and downstream edge of the cavity where a stratified charge is present. This trend is observed for different velocities at similar equivalence ratios. In case of methane combustion in the cavity, where the MFRs employed are extremely low at similar to 0.01, a different mechanism is observed. A fuel-rich mixture is now observed at the center of the cavity and this mixture undergoes combustion. On further increase of the cavity equivalence ratio, the rich mixture exceeds the flammability limit and forms a thin reaction zone at the interface with air stream. As a consequence, a shear layer flame at the top of the cavity interface with the mainstream is also observed. The equivalence ratio in the cavity also determines the combustion characteristics in the case of fuel-air mixtures that are formed as a result of the mixing. Overall, flame stabilization mechanisms have been proposed, which account for the wide range of MFRs and premixing in the mainstream as well.

Observation of correlated spin-orbit order in a strongly anisotropic quantum wire system

Quantum wires with spin-orbit coupling provide a unique opportunity to simultaneously control the coupling strength and the screened Coulomb interactions where new exotic phases of matter can be explored. Here we report on the observation of an exotic spin-orbit density wave in Pb-atomic wires on Si(557) surfaces by mapping out the evolution of the modulated spin-texture at various conditions with spin-and angle-resolved photoelectron spectroscopy. The results are independently quantified by surface transport measurements. The spin polarization, coherence length, spin dephasing rate and the associated quasiparticle gap decrease simultaneously as the screened Coulomb interaction decreases with increasing excess coverage, providing a new mechanism for generating and manipulating a spin-orbit entanglement effect via electronic interaction. Despite clear evidence of spontaneous spin-rotation symmetry breaking and modulation of spin-momentum structure as a function of excess coverage, the average spin polarization over the Brillouin zone vanishes, indicating that time-reversal symmetry is intact as theoretically predicted.

Spectral intensity distribution of trapped

To calculate static response properties of a many-body system, local density approximation (LDA) can be safely applied. But, to obtain dynamical response functions, the applicability of LDA is limited bacause dynamics of the system needs to be considered as well. To examine this in the context of cold atoms, we consider a system of non-interacting spin4 fermions confined by a harmonic trapping potential. We have calculated a very important response function, the spectral intensity distribution function (SIDF), both exactly and using LDA at zero temperature and compared with each other for different dimensions, trap frequencies and momenta. The behaviour of the SIDF at a particular momentum can be explained by noting the behaviour of the density of states (DoS) of the free system (without trap) in that particular dimension. The agreement between exact and LDA SIDFs becomes better with increase in dimensions and number of particles.

Thermalization of Green functions and quasinormal modes

We develop a new method to study the thermalization of time dependent retarded Green function in conformal field theories holographically dual to thin shell AdS Vaidya space times. The method relies on using the information of all time derivatives of the Green function at the shell and then evolving it for later times. The time derivatives of the Green function at the shell is given in terms of a recursion formula. Using this method we obtain analytic results for short time thermalization of the Green function. We show that the late time behaviour of the Green function is determined by the first quasinormal mode. We then implement the method numerically. As applications of this method we study the thermalization of the retarded time dependent Green function corresponding to a minimally coupled scalar in the AdS 3 and AdS 5 thin Vaidya shells. We see that as expected the late time behaviour is determined by the first quasinormal mode. We apply the method to study the late time behaviour of the shear vector mode in AdS 5 Vaidya shell. At small momentum the corresponding time dependent Green function is expected to relax to equilibrium by the shear hydrodynamic mode. Using this we obtain the universal ratio of the shear viscosity to entropy density from a time dependent process.

Quark and gluon distribution functions in a viscous quark-gluon plasma medium and dilepton production via q(q)over-bar annihilation

Viscous modifications to the thermal distributions of quark-antiquarks and gluons have been studied in a quasiparticle description of the quark-gluon-plasma medium created in relativistic heavy-ion collision experiments. The model is described in terms of quasipartons that encode the hot QCD medium effects in their respective effective fugacities. Both shear and bulk viscosities have been taken in to account in the analysis, and the modifications to thermal distributions have been obtained by modifying the energy-momentum tensor in view of the nontrivial dispersion relations for the gluons and quarks. The interactions encoded in the equation of state induce significant modifications to the thermal distributions. As an implication, the dilepton production rate in the q (q) over bar annihilation process has been investigated. The equation of state is found to have a significant impact on the dilepton production rate along with the viscosities.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

Scaling analysis and numerical studies on water vapour adsorption in a columnar porous silica gel bed

This article presents a theoretical analysis of heat and mass transfer in a silica gel + water adsorption process using scaling principles. A two-dimensional columnar packed adsorber domain is chosen for the study, with side and bottom walls cooled and vapour inlet from the top. The adsorption process is initiated from the cold walls with a temperature jump of 15 K, whereas the water vapour supply is maintained at a constant inlet pressure of 1 kPa. The first part of the study is dedicated to deriving relevant scales for the adsorption process by an order of magnitude analysis of energy, continuity and momentum equations. In the latter part, the derived scales are compared with the outcome of numerical studies performed for various domain widths and aspect ratio of bed. A good correlation between scaling and simulation results is observed, thereby validating the scaling approach. (C) 2015 Elsevier Ltd. All rights reserved.

The legacy of G. N. Ramachandran and the development of structural biology in India

G. N. Ramachandran is among the founding fathers of structural molecular biology. He made pioneering contributions in computational biology, modelling and what we now call bioinformatics. The triple helical coiled coil structure of collagen proposed by him forms the basis of much of collagen research at the molecular level. The Ramachandran map remains the simplest descriptor and tool for validation of protein structures. He has left his imprint on almost all aspects of biomolecular conformation. His contributions in the area of theoretical crystallography have been outstanding. His legacy has provided inspiration for the further development of structural biology in India. After a pause, computational biology and bioinformatics are in a resurgent phase. One of the two schools established by Ramachandran pioneered the development of macromolecular crystallography, which has now grown into an important component of modern biological research in India. Macromolecular NMR studies in the country are presently gathering momentum. Structural biology in India is now poised to again approach heights of the kind that Ramachandran conquered more than a generation ago.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

Sound radiation from a perforated panel set in a baffle with a different perforation ratio

A finite flexible perforated panel set in a differently perforated rigid baffle is considered. The radiation efficiency from such a panel is derived using a 2-D wavenumber domain formulation. This generalization is later used to represent a more practical case of a perforated panel fixed in an unperforated baffle. The perforations are in the form of an array of uniformly distributed circular holes. A complex impedance model for the holes available in the literature is used. An averaged fluid particle velocity is derived using the continuity equation and the surface pressure is derived using an appropriate momentum equation. The discontinuity in the perforate impedance (due to different hole dimensions or perforation ratio) at the panel-baffle interface is carefully taken into account. It is found that there exists a `coupling' of different wavenumbers of the spatially mean fluid particle velocity field. The change in the resonance frequencies and the modeshapes of the panel due to the perforations is taken into account using the Receptance method. Analytical expressions for the radiated power and radiation efficiency are derived in an integral form and numerical results are presented. Several comparisons are made to understand the radiation efficiency curves. Since both the resistive and reactive components of the hole impedance are taken into account, the model is directly applicable to micro-perforated panels also. (C) 2016 Elsevier Ltd. All rights reserved.

Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Laboratory Study on Sediment Diffusion and Deposition into Blind Channels

Any waterway with one end closed and the other open is generally called a blind channel. The main flow tends to expand, separate, and cause circulation at the mouth of blind channels. The main flow continuously transfers momentum and sediment into the circulation region through the turbulent mixing region (TMR) between them, thus leading to a large amount of sediment deposition in the blind channels. This paper experimentally investigated the properties of the water flow and sediment diffusion in TMR, demonstrating that both water flow and sediment motion in TMR approximately coincide with a similar structure as in the free mixing layer induced by a jet. The similarity functions of flow velocity and sediment concentration are then assumed, based on observation, and the resulting calculation of these functions is substantially facilitated. For the kind of low velocity flow system of blind channels with a finite width, a simple formula for the sediment deposition rate in blind channels is established by analyzing the gradient of crosswise velocity and sediment concentration in TMR.