Nonequilibrium processes from generalized Langevin equations: Realistic nanoscale systems connected to two thermal baths

We extend the generalized Langevin equation (GLE) method [L. Stella, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nosé-Hoover thermostats). We concentrate on the steady-state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e., ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Expansion tube operating conditions for studying nonequilibrium radiation relevant to Titan Aerocapture

The predictions of nonequilibrium radiation in the shock layer for a Titan aerocapture aeroshell vary significantly amongst Computational Fluid Dynamics (CFD) analyses and are limited by the physical models of the nonequilibrium flow processes. Of particular interest are nonequilibrium processes associated with the CN molecule which is a strong radiator. It is necessary to have experimental data for these radiating shock layers which will allow for validation of the CFD models. This paper describes the development of a test flow condition for subscale aeroshell models in a superorbital expansion tunnel. We discuss the need for a Titan gas condition that closely simulates the atmospheric composition and present experimental data of the free stream test flow conditions. Furthermore, we present finite-rate CFD calculations of the facility to estimate the remaining free stream conditions, which cannot be directly measured during experiments.

A kinetic model for nanocrystal morphology evolution

A kinetic model is developed with the goal of understanding and predicting the morphology evolution of nonocrystals in nonequilibrium growth conditions. The model is based on the assumption that under such conditions, different crystal planes have different kinetic parameters. This model focuses on the morphology-developing stage and is successfully related to the nucleation process and other crystal evolution mechanisms. It is believed to be a universal model and is applied to discuss the morphology evolution of CdSe nanocrystals, including the aspect ratio, injection I schemes, ligands effect and morphology distribution.

Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Implications of macromolecular crowding for signal transduction and metabolite channeling

The effect of different total enzyme concentrations on the flux through the bacterial phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) in vitro was determined by measuring PTS-mediated carbohydrate phosphorylation at different dilutions of cell-free extract of Escherichia coli. The dependence of the flux on the protein concentration was more than linear but less than quadratic. The combined flux–response coefficient of the four enzymes constituting the glucose PTS decreased slightly from values of ≈1.8 with increasing protein concentrations in the assay. Addition of the macromolecular crowding agents polyethylene glycol (PEG) 6000 and PEG 35000 led to a sharper decrease in the combined flux–response coefficient, in one case to values of ≈1. PEG 6000 stimulated the PTS flux at lower protein concentrations and inhibited the flux at higher protein concentrations, with the transition depending on the PEG 6000 concentration. This suggests that macromolecular crowding decreases the dissociation rate constants of enzyme complexes. High concentrations of the microsolute glycerol did not affect the combined flux–response coefficient. The data could be explained with a kinetic model of macromolecular crowding in a two-enzyme group-transfer pathway. Our results suggest that, because of the crowded environment in the cell, the different PTS enzymes form complexes that live long on the time-scale of their turnover. The implications for the metabolic behavior and control properties of the PTS, and for the effect of macromolecular crowding on nonequilibrium processes, are discussed.

Entropy, chaos, and excited-state quantum phase transitions in the Dicke model

We study nonequilibrium processes in an isolated quantum system-the Dicke model-focusing on the role played by the transition from integrability to chaos and the presence of excited-state quantum phase transitions. We show that both diagonal and entanglement entropies are abruptly increased by the onset of chaos. Also, this increase ends in both cases just after the system crosses the critical energy of the excited-state quantum phase transition. The link between entropy production, the development of chaos, and the excited-state quantum phase transition is more clear for the entanglement entropy.

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.

Nonequilibrium dynamics of photoexcited electrons in graphene: Collinear scattering, Auger processes, and the impact of screening

We present a combined analytical and numerical study of the early stages (sub-100-fs) of the nonequilibrium dynamics of photoexcited electrons in graphene. We employ the semiclassical Boltzmann equation with a collision integral that includes contributions from electron-electron (e-e) and electron-optical phonon interactions. Taking advantage of circular symmetry and employing the massless Dirac fermion (MDF) Hamiltonian, we are able to perform an essentially analytical study of the e-e contribution to the collision integral. This allows us to take particular care of subtle collinear scattering processes - processes in which incoming and outgoing momenta of the scattering particles lie on the same line - including carrier multiplication (CM) and Auger recombination (AR). These processes have a vanishing phase space for two-dimensional MDF bare bands. However, we argue that electron-lifetime effects, seen in experiments based on angle-resolved photoemission spectroscopy, provide a natural pathway to regularize this pathology, yielding a finite contribution due to CM and AR to the Coulomb collision integral. Finally, we discuss in detail the role of physics beyond the Fermi golden rule by including screening in the matrix element of the Coulomb interaction at the level of the random phase approximation (RPA), focusing in particular on the consequences of various approximations including static RPA screening, which maximizes the impact of CM and AR processes, and dynamical RPA screening, which completely suppresses them. © 2013 American Physical Society.

Cross-relaxation dynamics on anomalous saturation processes in low pressure supersonic cw HF chemical laser amplifier

It is assumed that both translational and rotational nonequilibrium cross-relaxations play a role simultaneoulsy in low pressure supersonic cw HF chemical laser amplifier. For two-type models of gas flow medium with laminar and turbulent flow diffusion mixing, the expressions of saturated gain spectrum are derived respectively, and the numerical calculations are performed as well. The numerical results show that turbulent flow diffusion mixing model is in the best agreement with the experimental result.

Assessing the Nonequilibrium Thermodynamics in a Quenched Quantum Many-Body System via Single Projective Measurements

We analyze the nature of the statistics of the work done on or by a quantum many-body system brought out of equilibrium. We show that, for the sudden quench and for an initial state that commutes with the initial Hamiltonian, it is possible to retrieve the whole nonequilibrium thermodynamics via single projective measurements of observables. We highlight, in a physically clear way, the qualitative implications for the statistics of work coming from considering processes described by operators that either commute or do not commute with the unperturbed Hamiltonian of a given system. We consider a quantum many-body system and derive an expression that allows us to give a physical interpretation, for a thermal initial state, to all of the cumulants of the work in the case of quenched operators commuting with the unperturbed Hamiltonian. In the commuting case, the observables that we need to measure have an intuitive physical meaning. Conversely, in the noncommuting case, we show that, although it is possible to operate fully within the single-measurement framework irrespectively of the size of the quench, some difficulties are faced in providing a clear-cut physical interpretation to the cumulants. This circumstance makes the study of the physics of the system nontrivial and highlights the nonintuitive phenomenology of the emergence of thermodynamics from the fully quantum microscopic description. We illustrate our ideas with the example of the Ising model in a transverse field showing the interesting behavior of the high-order statistical moments of the work distribution for a generic thermal state and linking them to the critical nature of the model itself.

Real-time simulation of nonequilibrium transport of magnetization in large open quantum spin systems driven by dissipation

Using quantum Monte Carlo, we study the nonequilibrium transport of magnetization in large open strongly correlated quantum spin-12 systems driven by purely dissipative processes that conserve the uniform or staggered magnetization, disregarding unitary Hamiltonian dynamics. We prepare both a low-temperature Heisenberg ferromagnet and an antiferromagnet in two parts of the system that are initially isolated from each other. We then bring the two subsystems in contact and study their real-time dissipative dynamics for different geometries. The flow of the uniform or staggered magnetization from one part of the system to the other is described by a diffusion equation that can be derived analytically.