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.

Intense-field dissociation dynamics of D2+ molecular ions using ultrafast laser pulses

The dynamics of dissociation of pre-ionized D2+ molecules using intense (10^12–10^15 W cm-2), ultrashort (50 fs), infrared (? = 790 nm) laser pulses are examined. Use of an intensity selective scan technique has allowed the deuterium energy spectrum to be measured over a broad range of intensity. It is found that the dominant emission shifts to lower energies as intensity is increased, in good agreement with corresponding wavepacket simulations. The results are consistent with an interpretation in terms of bond softening, which at high intensity (approximately >3 × 10^14 W cm-2) becomes dominated by dissociative ionization. Angular distribution measurements reveal the presence of slow molecular dissociation, an indication that vibrational trapping mechanisms occur in this molecule.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Dynamics of nonlinearly coupled magnetic-field-aligned electromagnetic electron-cyclotron waves near the zero-group-dispersion point in magnetized plasmas

The nonlinear coupling between two magnetic-field-aligned electromagnetic electron-cyclotron (EMEC) waves in plasmas is considered. Evaluating the ponderomotive coupling between the EMEC waves and quasistationary plasma density perturbations, a pair of coupled nonlinear Schrodinger equations (CNLSEs) is obtained. The CNLSEs are then used to investigate the occurrence of modulational instability in magnetized plasmas. Waves in the vicinity of the zero-group-dispersion point are considered, so that the group dispersion terms may either bear the same or different signs. It is found that a stable EMEC wave can be destabilized due to its nonlinear interactions with an unstable one, while a pair of unstable EMEC waves yields an increased instability growth rate. Individually stable waves remain stable while interacting with one another. Stationary nonlinear solutions of the coupled equations are presented. The relevance of our investigation to nonlinear phenomena in space plasmas is discussed. (c) 2005 American Institute of Physics.

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.

Fragility of the spin-glass-like collective state to a magnetic field in an interacting Fe-C nanoparticle system

The effect of applied magnetic fields on the collective nonequilibrium dynamics of a strongly interacting Fe-C nanoparticle system has been investigated. It is experimentally shown that the magnetic aging diminishes to finally disappear for fields of moderate strength. The field needed to remove the observable aging behavior increases with decreasing temperature. The same qualitative behavior is observed in an amorphous metallic spin glass (Fe0.15Ni0.85)(75)P16B6Al3.

An Investigation of the Flow Field through a Variable Geometry Turbine Stator with Vane Endwall Clearance

Variable geometry turbines provide an extra degree of flexibility in air management in turbocharged engines. The pivoting stator vanes used to achieve the variable turbine geometry necessitate the inclusion of stator vane endwall clearances. The consequent leakage flow through the endwall clearances impacts the flow in the stator vane passages and an understanding of the impact of the leakage flow on stator loss is required. A numerical model of a typical variable geometry turbine was developed using the commercial CFX-10 computational fluid dynamics software, and validated using laser doppler velocimetry and static pressure measurements from a variable geometry turbine with stator vane endwall clearance. Two different stator vane positions were investigated, each at three different operating conditions representing different vane loadings. The vane endwall leakage was found to have a significant impact on the stator loss and on the uniformity of flow entering the turbine rotor. The leakage flow changed considerably at different vane positions and flow incidence at vane inlet was found to have a significant impact.

Reconstructing the dynamics of a movable mirror in a detuned optical cavity

We consider the dynamics of a movable mirror in a Fabry-Perot cavity coupled through radiation pressure to the cavity field and in contact with a thermal bath at finite temperature. In contrast to previous approaches, we consider arbitrary values of the effective detuning between the cavity and an external input field. We analyse the radiation-pressure effect on the Brownian motion of the mirror and its significance in the density noise spectrum of the output cavity field. Important properties of the mirror dynamics can be gathered directly from this noise spectrum. The presented reconstruction provides an experimentally useful tool in the characterization of the energy and rigidity of the mirror as modified by the coupling with light. We also give a quantitative analysis of the recent experimental observation of self-cooling of a micromechanical oscillator.

A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.