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.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

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.

Berry phase in open quantum systems: a quantum Langevin equation approach

The evolution of a two level system with a slowly varying Hamiltonian, modeled as a spin 1/2 in a slowly varying magnetic field, and interacting with a quantum environment, modeled as a bath of harmonic oscillators is analyzed using a quantum Langevin approach. This allows to easily obtain the dissipation time and the correction to the Berry phase in the case of an adiabatic cyclic evolution.

CRITICAL PROCESSES, LANGEVIN EQUATION AND UNIVERSALITY

In nature there are ubiquitous systems that can naturally approach critical states, The Langevin equation in the discrete version can be used to describe a class of critical processes, which are characterized by power-law behaviors and scaling relations. As an example, we present a simple model for a clinical thermometer, whose reading cannot fall even when its temperature decreases. The fibers bundle model and the spring-block model are also shown to belong to such a class.

Current-induced atomic dynamics, instabilities, and Raman signals: Quasiclassical Langevin equation approach

We derive and employ a semiclassical Langevin equation obtained from path integrals to describe the ionic dynamics of a molecular junction in the presence of electrical current. The electronic environment serves as an effective nonequilibrium bath. The bath results in random forces describing Joule heating, current-induced forces including the nonconservative wind force, dissipative frictional forces, and an effective Lorentz-type force due to the Berry phase of the nonequilibrium electrons. Using a generic two-level molecular model, we highlight the importance of both current-induced forces and Joule heating for the stability of the system. We compare the impact of the different forces, and the wide-band approximation for the electronic structure on our result. We examine the current-induced instabilities (excitation of runaway "waterwheel" modes) and investigate the signature of these in the Raman signals.

Fokker-Planck equation for a test-particle weakly coupled to a magnetized plasma

We consider the derivation of a kinetic equation for a charged test particle weakly interacting with an electrostatic plasma in thermal equilibrium, subject to a uniform external magnetic field. The Liouville equation leads to a generalized master equation to second order in the `weak' interaction; a Fokker-Planck-type equation then follows as a `Markovian' approximation. It is shown that such an equation does not preserve the positivity of the distribution function f(x,v;t). By applying techniques developed in the theory of open systems, a correct Fokker-Planck equation is derived. Explicit expressions for the diffusion and drift coefficients, depending on the magnetic field, are obtained.

Statistical-mechanical description of classical test-particle dynamics in the presence of an external force field: modelling noise and damping from first principles

Aiming to establish a rigorous link between macroscopic random motion (described e.g. by Langevin-type theories) and microscopic dynamics, we have undertaken a kinetic-theoretical study of the dynamics of a classical test-particle weakly coupled to a large heat-bath in thermal equilibrium. Both subsystems are subject to an external force field. From the (time-non-local) generalized master equation a Fokker-Planck-type equation follows as a "quasi-Markovian" approximation. The kinetic operator thus defined is shown to be ill-defined; in specific, it does not preserve the positivity of the test-particle distribution function f(x, v; t). Adopting an alternative approach, previously introduced for quantum open systems, is proposed to lead to a correct kinetic operator, which yields all the expected properties. A set of explicit expressions for the diffusion and drift coefficients are obtained, allowing for modelling macroscopic diffusion and dynamical friction phenomena, in terms of an external field and intrinsic physical parameters.

Symmetries and exact solutions of KP equation with an arbitrary nonlinear term

The generalized KP (GKP) equations with an arbitrary nonlinear term model and characterize many nonlinear physical phenomena. The symmetries of GKP equation with an arbitrary nonlinear term are obtained. The condition that must satisfy for existence the symmetries group of GKP is derived and also the obtained symmetries are classified according to different forms of the nonlinear term. The resulting similarity reductions are studied by performing the bifurcation and the phase portrait of GKP and also the corresponding solitary wave solutions of GKP
equation are constructed.

Score for neonatal acute physiology-II and neonatal pain predict corticospinal tract development in premature newborns

Premature infants are at risk for adverse motor outcomes, including cerebral palsy and developmental coordination disorder. The purpose of this study was to examine the relationship of antenatal, perinatal, and postnatal risk factors for abnormal development of the corticospinal tract, the major voluntary motor pathway, during the neonatal period. In a prospective cohort study, 126 premature neonates (24-32 weeks' gestational age) underwent serial brain imaging near birth and at term-equivalent age. With diffusion tensor tractography, mean diffusivity and fractional anisotropy of the corticospinal tract were measured to reflect microstructural development. Generalized estimating equation models examined associations of risk factors on corticospinal tract development. The perinatal risk factor of greater early illness severity (as measured by the Score for Neonatal Acute Physiology-II [SNAP-II]) was associated with a slower rise in fractional anisotropy of the corticospinal tract (P = 0.02), even after correcting for gestational age at birth and postnatal risk factors (P = 0.009). Consistent with previous findings, neonatal pain adjusted for morphine and postnatal infection were also associated with a slower rise in fractional anisotropy of the corticospinal tract (P = 0.03 and 0.02, respectively). Lessening illness severity in the first hours of life might offer potential to improve motor pathway development in premature newborns.

Procedural pain and brain development in premature newborns

Objective: Preterm infants are exposed to multiple painful procedures in the neonatal intensive care unit (NICU) during a period of rapid brain development. Our aim was to examine relationships between procedural pain in the NICU and early brain development in very preterm infants.

Methods: Infants born very preterm (N ¼ 86; 24–32 weeks gestational age) were followed prospectively from birth, and studied with magnetic resonance imaging, 3-dimensional magnetic resonance spectroscopic imaging, and diffusion tensor imaging: scan 1 early in life (median, 32.1 weeks) and scan 2 at term-equivalent age (median, 40 weeks). We calculated N-acetylaspartate to choline ratios (NAA/choline), lactate to choline ratios, average diffusivity, and white matter fractional anisotropy (FA) from up to 7 white and 4 subcortical gray matter regions of interest. Procedural pain was quantified as the number of skin-breaking events from birth to term or scan 2. Data were
analyzed using generalized estimating equation modeling adjusting for clinical confounders such as illness severity, morphine exposure, brain injury, and surgery.

Results: After comprehensively adjusting for multiple clinical factors, greater neonatal procedural pain was associated with reduced white matter FA (b ¼ 0.0002, p ¼ 0.028) and reduced subcortical gray matter NAA/choline (b ¼ 0.0006, p ¼ 0.004). Reduced FA was predicted by early pain (before scan 1), whereas lower NAA/choline was predicted by pain exposure throughout the neonatal course, suggesting a primary and early effect on subcortical structures with secondary white matter changes.

Interpretation: Early procedural pain in very preterm infants may contribute to impaired brain development.

Changes in epilithic biomasses and invertebrate community structure over a deposit metal concentration gradient in upland headwater streams

Stream bed metal deposits affect the taxon richness, density and taxonomic diversity of primary and secondary producers by a variety of direct or indirect abiotic and biotic processes but little is known about the relative importance of these processes over a deposit metal concentration gradient. Inorganic matter (IM), algal and non-photosynthetic detrital (NPD) dry biomasses were estimated for 10 monthly samples, between 2007 and 2008, from eight sites differing in deposit density. Invertebrate abundance, taxon richness and composition were also determined. Relations between these variables were investigated by canonical correspondence analysis (CCA), generalized estimating equation models and path analysis. The first CCA axis correlates with deposit density and invertebrate abundance, with lumbriculids and chironomids increasing in abundance with deposit density and all other taxa declining. Community structure changes significantly above a deposit density of approximately 8 mg cm, when algal biomass, invertebrate richness and diversity decline. Invertebrate richness and diversity were determined by direct effects of NPD biomass and indirect effects of IM. Algal biomass only had an effect on invertebrate abundance. Possible pH, oxygen, food and ecotoxicological effects of NPD biomass on the biota are discussed.

Central corneal thickness and corneal hysteresis associated with glaucoma damage.

PURPOSE:

We sought to measure the impact of central corneal thickness (CCT), a possible risk factor for glaucoma damage, and corneal hysteresis, a proposed measure of corneal resistance to deformation, on various indicators of glaucoma damage.

DESIGN:

Observational study.

METHODS:

Adult patients of the Wilmer Glaucoma Service underwent measurement of hysteresis on the Reichert Ocular Response Analyzer and measurement of CCT by ultrasonic pachymetry. Two glaucoma specialists (H.A.Q., N.G.C.) reviewed the chart to determine highest known intraocular pressure (IOP), target IOP, diagnosis, years with glaucoma, cup-to-disk ratio (CDR), mean defect (MD), pattern standard deviation (PSD), glaucoma hemifield test (GHT), and presence or absence of visual field progression.

RESULTS:

Among 230 subjects, the mean age was 65 +/- 14 years, 127 (55%) were female, 161 (70%) were white, and 194 (85%) had a diagnosis of primary open-angle glaucoma (POAG) or suspected POAG. In multivariate generalized estimating equation models, lower corneal hysteresis value (P = .03), but not CCT, was associated with visual field progression. When axial length was included in the model, hysteresis was not a significant risk factor (P = .09). A thinner CCT (P = .02), but not hysteresis, was associated with a higher CDR at the most recent examination. Neither CCT nor hysteresis was associated with MD, PSD, or GHT "outside normal limits."

CONCLUSIONS:

Thinner CCT was associated with the state of glaucoma damage as indicated by CDR. Axial length and corneal hysteresis were associated with progressive field worsening.

Corneal biomechanics, refractive error, and axial length in Chinese primary school children.

PURPOSE: Low corneal hysteresis is associated with longer axial length in Chinese secondary school children. The authors sought to explore this association in primary school children. METHODS: LogMAR presenting visual acuity, cycloplegic refractive error, ocular biometry, central corneal thickness (CCT), and corneal hysteresis (CH) was assessed for children in grades 1 to 3 at an academically competitive urban school in Shantou, China. RESULTS: Among 872 eligible children (mean age, 8.6 ± 2.1 years), 651 (74.7%) completed the examination. Among 1299 examined eyes, 111 (8.5%) had uncorrected vision ≤6/12. Mean spherical equivalent refractive error for all eyes was +0.26 ± 1.41 D, and axial length (AL) was 22.7 ± 0.90 mm. CH for the lowest (mean AL, 21.7 ± 0.39 mm), two middle (mean AL, 22.4 ± 0.15 and 22.9 ± 0.15 mm), and highest quartiles (mean AL, 23.7 ± 0.74 mm) of AL were 10.6 ± 2.1 mm Hg, 10.4 ± 2.1 mm Hg, 10.3 ± 2.3 mm Hg, and 10.2 ± 2.3 mm Hg respectively (age- and gender-adjusted Pearson's correlation coefficient r = -0.052; P = 0.001). In generalized estimating equation models adjusting for age, gender, and CCT, lower CH was significantly associated with longer AL (P < 0.001) and more myopic refractive error (P = 0.001). CONCLUSIONS: CH measurement is practical in young children because this is when myopia undergoes its most rapid progression. Prospective follow-up of this cohort at high risk for myopia is under way to determine whether low CH is predictive, or a consequence, of long AL.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.