Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse 10(d) x 10(d) matrix, with d > 6. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in nonequilibrium setups without mass transport.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Nonequilibrium noise in electrophoresis: The microion wind

A colloid supported against gravitational settling by means of an imposed electric field behaves, on average, as if it is at equilibrium in a confining potential T. M. Squires, J. Fluid Mech. 443, 403 (2001)]. We show, however, that the effective Langevin equation for the colloid contains a nonequilibrium noise source, proportional to the field, arising from the thermal motion of dissolved ions. The position fluctuations of the colloid show strong, experimentally testable signatures of nonequilibrium behavior, including a highly anisotropic, frequency-dependent ``effective temperature'' obtained from the fluctuation-dissipation ratio.

Nonequilibrium ``Melting'' of a Charge Density Wave Insulator via an Ultrafast Laser Pulse

We employ an exact solution of the simplest model for pump-probe time-resolved photoemission spectroscopy in charge-density-wave systems to show how, in nonequilibrium, the gap in the density of states disappears while the charge density remains modulated, and then the gap reforms after the pulse has passed. This nonequilibrium scenario qualitatively describes the common short-time experimental features in TaS2 and TbTe3, indicating a quasiuniversality for nonequilibrium ``melting'' with qualitative features that can be easily understood within a simple picture.

Experimental signatures of a nonequilibrium phase transition governing the yielding of a soft glass

We present direct experimental signatures of a nonequilibrium phase transition associated with the yield point of a prototypical soft solid-a binary colloidal glass. By simultaneously quantifying single-particle dynamics and bulk mechanical response, we identified the threshold for the onset of irreversibility with the yield strain. We extracted the relaxation time from the transient behavior of the loss modulus and found that it diverges in the vicinity of the yield strain. This critical slowing down is accompanied by a growing correlation length associated with the size of regions of high Debye-Waller factor, which are precursors to yield events in glasses. Our results affirm that the paradigm of nonequilibrium critical phenomena is instrumental in achieving a holistic understanding of yielding in soft solids.

Study of motion around a static black hole in Einstein and Lovelock gravity

We study motion around a static Einstein and pure Lovelock black hole in higher dimensions. It is known that in higher dimensions bound orbits exist only for a pure Lovelock black hole in all even dimensions, D = 2N + 2, where N is the degree of Lovelock polynomial action. In particular, we compute periastron shift and light bending, and the latter is given by one of the transverse spatial components of the Riemann curvature tensor. We also consider the pseudo-Newtonian potentials and Kruskal coordinates.

Nonequilibrium Microstructures for Ag-Ni Nanowires

This work illustrates that a variety of nanowire microstructures can be obtained either by controlling the nanowire formation kinetics or by suitable thermal processing of as-deposited nanowires with nonequilibrium metastable microstructure. In the present work, 200-nm diameter Ag-Ni nanowires with similar compositions, but with significantly different microstructures, were electrodeposited. A 15 mA deposition current produced nanowires in which Ag-rich crystalline nanoparticles were embedded in a Ni-rich amorphous matrix. A 3 mA deposition current produced nanowires in which an Ag-rich crystalline phase formed a backbone-like configuration in the axial region of the nanowire, whereas the peripheral region contained Ni-rich nanocrystalline and amorphous phases. Isothermal annealing of the nanowires illustrated a phase evolution pathway that was extremely sensitive to the initial nanowire microstructure.

Modified Einstein's gravity as a possible missing link between sub- and super-Chandrasekhar type Ia supernovae

We explore the effect of modification to Einstein's gravity in white dwarfs for the first time in the literature, to the best of our knowledge. This leads to significantly sub- and super-Chandrasekhar limiting masses of white dwarfs, determined by a single model parameter. On the other hand, type Ia supernovae (SNeIa), a key to unravel the evolutionary history of the universe, are believed to be triggered in white dwarfs having mass close to the Chandrasekhar limit. However, observations of several peculiar, under- and over-luminous SNeIa argue for exploding masses widely different from this limit. We argue that explosions of the modified gravity induced sub- and super-Chandrasekhar limiting mass white dwarfs result in under- and over-luminous SNeIa respectively, thus unifying these two apparently disjoint sub-classes and, hence, serving as a missing link. Our discovery raises two fundamental questions. Is the Chandrasekhar limit unique? Is Einstein's gravity the ultimate theory for understanding astronomical phenomena? Both the answers appear to be no!

Imprint of modified Einstein's gravity on white dwarfs: Unifying Type Ia supernovae

We establish the importance of modified Einstein's gravity (MG) in white dwarfs (WDs) for the first time in the literature. We show that MG leads to significantly sub- and super-Chandrasekhar limiting mass WDs, depending on a single model parameter. However, conventional WDs on approaching Chandrasekhar's limit are expected to trigger Type Ia supernovae (SNeIa), a key to unravel the evolutionary history of the universe. Nevertheless, observations of several peculiar, under-and over-luminous SNeIa argue for the limiting mass widely different from Chandrasekhar's limit. Explosions of MG induced sub-and super-Chandrasekhar limiting mass WDs explain under-and over-luminous SNeIa respectively, thus unifying these two apparently disjoint sub-classes. Our discovery questions both the global validity of Einstein's gravity and the uniqueness of Chandrasekhar's limit.

On nonequilibrium of pyrolysis process in the manufacture of ethene

Nonequilibrium process for cracking ethane and n-buthane in the manufacture of ethene has been analytically and numerically investigated in a Heavi-side function temperature field and through a normal shock wave. The results demonstrate that, while the reaction temperature increases, the maximum value of ethene yield is increased, and the optimal reaction duration is sharply shortened. For the identical initial reaction temperature, the maximum value of ethene yield through a stationary normal shock wave is less than that in a Heavi-side function temperature field. However, the ethene consumption after the maximum value in the former case is less than that in the latter. Higher ethene yield will be obtained by using the gasdynamic heating method than by using the current methods.

Nonequilibrium microstructures and their evolution in a Fe-Cr-W-Ni-C laser clad coating

The laser-solidified microstructural and compositional characterization and phase evolution during tempering at 963 K were investigated using an analytical transmission electron microscope with energy dispersive X-ray analysis. The cladded alloy, a powder mixture of Fe, Cr, W, Ni, and C with a weight ratio of 10:5:1:1:1, was processed with a 3 kW continuous wave CO2 laser. The processing parameters were 16 mm/s beam scanning speed, 3 mm beam diameter. 2 kW laser power, and 0.3 g/s feed rate. The coating was metallurgically bonded to the substrate, with a maximum thickness of 730 mu m, a microhardness of about 860 Hv and a volumetric dilution ratio of about 6%. Microanalyses revealed that the cladded coating possessed the hypoeutectic microstructure comprising the primary dendritic gamma-austenite and interdendritic eutectic consisted of gamma-austenite and M7C3 carbide. The gamma-austenite was a non-equilibrium phase with extended solid solution of alloying elements and a great deal of defect structures, i.e. a high density of dislocations, twins, and stacking faults existed in gamma phase. During high temperature aging, in situ carbide transformation occurred of M7C3 to M23C6 and M6C. The precipitation of M23C6, MC and M2C carbides from austenite was also observed.

Partial-Dislocation-Mediated Processes in Nanocrystalline Ni with Nonequilibrium Grain Boundaries

The partial-dislocation-mediated processes have so far eluded high-resolution transmission electron microscopy studies in nanocrystalline nc Ni with nonequilibrium grain boundaries. It is revealed that the nc Ni deformed largely by twinning instead of extended partials. The underlying mechanisms including dissociated dislocations, high residual stresses, and stress concentrations near stacking faults are demonstrated and discussed.

Implication of scaling hierarchy associated with nonequilibrium: field and particulate

The advent of nanotechnology has necessitated a better understanding of how material microstructure changes at the atomic level would affect the macroscopic properties that control the performance. Such a challenge has uncovered many phenomena that were not previously understood and taken for granted. Among them are the basic foundation of dislocation theories which are now known to be inadequate. Simplifying assumptions invoked at the macroscale may not be applicable at the micro- and/or nanoscale. There are implications of scaling hierrachy associated with in-homegeneity and nonequilibrium. of physical systems. What is taken to be homogeneous and equilibrium at the macroscale may not be so when the physical size of the material is reduced to microns. These fundamental issues cannot be dispensed at will for the sake of convenience because they could alter the outcome of predictions. Even more unsatisfying is the lack of consistency in modeling physical systems. This could translate to the inability for identifying the relevant manufacturing parameters and rendering the end product unpractical because of high cost. Advanced composite and ceramic materials are cases in point. Discussed are potential pitfalls for applying models at both the atomic and continuum levels. No encouragement is made to unravel the truth of nature. Let it be partiuclates, a smooth continuum or a combination of both. The present trend of development in scaling tends to seek for different characteristic lengths of material microstructures with or without the influence of time effects. Much will be learned from atomistic simulation models to show how results could differ as boundary conditions and scales are changed. Quantum mechanics, continuum and cosmological models provide evidence that no general approach is in sight. Of immediate interest is perhaps the establishment of greater precision in terminology so as to better communicate results involving multiscale physical events.

Effects of Translational Nonequilibrium on Performance of Flowing Chemical Oxygen-Iodine Laser

The effect of the translational nonequilibrium on performance modeling of flowing chemical oxygen-iodine lasers (COIL) is emphasized in this paper. The spectral line broadening (SLB) model is a basic factor for predicting the performances of flowing COIL. The Voigt profile function is a well-known SLB model and is usually utilized. In the case of gas pressure in laser cavity less than 5 torr, a low pressure limit expression of the Voigt profile function is used. These two SLB models imply that ail lasing particles can interact with monochromatic laser radiation. Basically, the inhomogeneous broadening effects are not considered in these two SLB models and they cannot predict the spectral content. The latter requires consideration of finite translational relaxation rate. Unfortunately, it is rather difficult to solve simultaneously the Navier-Stokes (NS) equations and the conservation equations of the number of lasing particles per unit volume and per unit frequency interval. In the operating condition of flowing COIL, it is possible to obtain a perturbational solution of the conservational equations for lasing particles and deduce a new relation between the gain and the optical intensity, i.e., a new gain-saturation relation. By coupling the gain-saturation relation with other governing equations (such as the NS equations, chemical reaction equations and the optical model of gain-equal-loss), We have numerically calculated the performances of flowing COIL. The present results are compared with those obtained by the common rate-equation (RE) model, in which the Voigt profile function and its low pressure limit expression are used. The difference of different model's results is great. For instance, in the case of lasing frequency coinciding with the central frequency of line profile and very low gas pressure, the gain-saturation relation of the present model is quite different with that of the RE model.