Solution of the master equation for Wigner's quasiprobability distribution in phase space for the Brownian motion of a particle in a double well potential

Quantum effects in the Brownian motion of a particle in the symmetric double well potential V(x)=ax(2)2+bx(4)4 are treated using the semiclassical master equation for the time evolution of the Wigner distribution function W(x,p,t) in phase space (x,p). The equilibrium position autocorrelation function, dynamic susceptibility, and escape rate are evaluated via matrix continued fractions in the manner customarily used for the classical Fokker-Planck equation. The escape rate so yielded has a quantum correction depending strongly on the barrier height and is compared with that given analytically by the quantum mechanical reaction rate solution of the Kramers turnover problem. The matrix continued fraction solution substantially agrees with the analytic solution. Moreover, the low-frequency part of the spectrum associated with noise assisted Kramers transitions across the potential barrier may be accurately described by a single Lorentzian with characteristic frequency given by the quantum mechanical reaction rate.     

Solution of quantum Langevin equation: approximations, theoretical and numerical aspects

Based on a coherent state representation of noise operator and an ensemble averaging procedure using Wigner canonical thermal distribution for harmonic oscillators, a generalized quantum Langevin equation has been recently developed [Phys. Rev. E 65, 021109 (2002); 66, 051106 (2002)] to derive the equations of motion for probability distribution functions in c-number phase-space. We extend the treatment to explore several systematic approximation schemes for the solutions of the Langevin equation for nonlinear potentials for a wide range of noise correlation, strength and temperature down to the vacuum limit. The method is exemplified by an analytic application to harmonic oscillator for arbitrary memory kernel and with the help of a numerical calculation of barrier crossing, in a cubic potential to demonstrate the quantum Kramers' turnover and the quantum Arrhenius plot.     

Non-Markovian theory of open systems in classical limit

A fully classical limit of the recently published quantum-classical approximation [A. A. Neufeld, J. Chem. Phys. 119, 2488 (2003)] is obtained and analyzed. The resulting kinetic equations are capable of describing the evolution of an open system on the entire time axis, including the short-time non-Markovian stage, and are valid beyond linear response regime. We have shown, that proceeding to the classical mechanics limit we restrict the class of allowed correlations between an open system and a canonical bath, so that the initial conditions and the relaxation operator has to be appropriately modified (projected). Disregard of the projection may lead to unphysical behavior, since mechanism of the decay of some correlations is essentially of quantum-mechanical nature, and is not correctly described by classical mechanics. The projection (quantum correction to the kinetics) is particularly important for the non-Markovian regime of relaxation towards canonical equilibrium. The conformity of the developed method to the conventional approaches is demonstrated using a model of Brownian motion (heavy particle in the bath of light ones), for which the obtained non-Markovian equations are reduced to the standard Fokker-Planck equation in phase space.     

Real time correlation function in a single phase space integral beyond the linearized semiclassical initial value representation

It is shown how quantum mechanical time correlation functions [defined, e.g., in Eq. (1.1)] can be expressed, without approximation, in the same form as the linearized approximation of the semiclassical initial value representation (LSC-IVR), or classical Wigner model, for the correlation function [cf. Eq. (2.1)], i.e., as a phase space average (over initial conditions for trajectories) of the Wigner functions corresponding to the two operators. The difference is that the trajectories involved in the LSC-IVR evolve classically, i.e., according to the classical equations of motion, while in the exact theory they evolve according to generalized equations of motion that are derived here. Approximations to the exact equations of motion are then introduced to achieve practical methods that are applicable to complex (i.e., large) molecular systems. Four such methods are proposed in the paper--the full Wigner dynamics (full WD) and the second order WD based on "Wigner trajectories" [H. W. Lee and M. D. Scully, J. Chem. Phys. 77, 4604 (1982)] and the full Donoso-Martens dynamics (full DMD) and the second order DMD based on "Donoso-Martens trajectories" [A. Donoso and C. C. Martens, Phys. Rev. Lett. 8722, 223202 (2001)]--all of which can be viewed as generalizations of the original LSC-IVR method. Numerical tests of the four versions of this new approach are made for two anharmonic model problems, and for each the momentum autocorrelation function (i.e., operators linear in coordinate or momentum operators) and the force autocorrelation function (nonlinear operators) have been calculated. These four new approximate treatments are indeed seen to be significant improvements to the original LSC-IVR approximation.     

Two more approaches for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics

We show two more approaches for generating trajectory-based dynamics in the phase space formulation of quantum mechanics: "equilibrium continuity dynamics" (ECD) in the spirit of the phase space continuity equation in classical mechanics, and "equilibrium Hamiltonian dynamics" (EHD) in the spirit of the Hamilton equations of motion in classical mechanics. Both ECD and EHD can recover exact thermal correlation functions (of even nonlinear operators, i.e., nonlinear functions of position or momentum operators) in the classical, high temperature, and harmonic limits. Both ECD and EHD conserve the quasi-probability within the infinitesimal volume dx(t)dp(t) around the phase point (x(t), p(t)) along the trajectory. Numerical tests of both approaches in the Wigner phase space have been made for two strongly anharmonic model problems and a double well system, for each potential auto-correlation functions of both linear and nonlinear operators have been calculated. The results suggest EHD and ECD are two additional potential useful approaches for describing quantum effects for complex systems in condense phase.     

Fewest switches adiabatic surface hopping as applied to vibrational energy relaxation

In this contribution quantum/classical surface hopping methodology is applied to vibrational energy relaxation of a quantum oscillator in a classical heat bath. The model of a linearly damped (harmonic) oscillator is chosen which can be mapped onto the Brownian motion (Caldeira-Leggett) Hamiltonian. In the simulations Tully's fewest switches surface hopping scheme is adopted with inclusion of dephasing in the adiabatic basis using a simple decoherence algorithm. The results are compared to the predictions of a Redfield-type quantum master equation modeling using the classical heat bath force correlation function as input. Thereby a link is established between both types of quantum/classical approaches. Viewed from the latter perspective, surface hopping with dephasing may be interpreted as "on-the-fly" stochastic realization of a quantum/classical Pauli master equation.     

Projection of two-dimensional diffusion in a curved midline and narrow varying width channel onto the longitudinal dimension

This study focuses on the derivation of a general effective diffusion coefficient to describe the two-dimensional (2D) diffusion in a narrow and smoothly asymmetric channel of varying width, in the simple diffusional motion of noninteracting pointlike particles under no external field. We present a generalization to the case of an asymmetric channel using the projection method introduced earlier by Kalinay and Percus [J. Chem. Phys. 122, 204701 (2005); and Phys. Rev. E 74, 041203 (2006)] to project the 2D diffusion equation into an effective one-dimensional generalized Fick-Jacobs equation. The expression for the diffusion coefficient given in Eq. (23) is our main result. This expression is a more general effective diffusion coefficient for narrow channels in 2D, which contains the well-known previous results as special cases, namely, those obtained by Bradley [Phys. Rev. E 80, 061142 (2009)], and more recently by Berezhkovskii and Szabo [J. Chem. Phys. 135, 074108 (2011)]. Finally, we study some specific 2D asymmetric channel configurations to test and show the broader applicability of Eq. (23).     

Superposition of quantum and classical rotational motions in Sc2C2@C84 fullerite

The superposition of the quantum rotational motion (tunneling) of the encapsulated Sc(2)C(2) complex with the classical rotational motion of the surrounding C(84) molecule in a powder crystal of Sc(2)C(2)@C(84) fullerite is investigated by theory. Since the quantum rotor is dragged along by the C(84) molecule, any detection method which couples to the quantum rotor (in casu the C(2) bond of the Sc(2)C(2) complex) also probes the thermally excited classical motion (uniaxial rotational diffusion and stochastic meroaxial jumps) of the surrounding fullerene. The dynamic rotation-rotation response functions in frequency space are obtained as convolutions of quantum and classical dynamic correlation functions. The corresponding Raman scattering laws are derived, and the overall shape of the spectra and the width of the resonance lines are studied as functions of temperature. The results of the theory are confronted with experimental low-frequency Raman spectra on powder crystals of Sc(2)C(2)@C(84) [M. Krause et al., Phys. Rev. Lett. 93, 137403 (2004)]. The agreement of theory with experiment is very satisfactory in a broad temperature range.     

Direct calculation of (1)H(2)O T(1) NMRD profiles and EPR lineshapes for the electron spin quantum numbers S = 1, 3/2, 2, 5/2, 3, 7/2, based on the stochastic Liouville equation combined with Brownian dynamics simulation

Direct calculation of electron spin relaxation and EPR lineshapes, based on Brownian dynamics simulation techniques and the stochastic Liouville equation approach (SLE-L) [Mol. Phys., 2004, 102, 1085-1093], is here generalized to high spin systems with spin quantum number S = 3/2, 2, 5/2, 3 and 7/2. A direct calculation method is demonstrated for electron spin-spin and spin-lattice relaxation, S-, X- and Q-band EPR-lineshapes and paramagnetic enhanced water proton T(1)- NMRD profiles. The main relaxation mechanism for the electron spin system is a stochastic second rank zero field splitting (ZFS). Brownian dynamics simulation techniques are used in describing a fluctuating ZFS interaction which comprises two parts namely the "permanent" part which is modulated by isotropic reorientation diffusion, and the transient part which is modulated by fast local distortion, which is also modelled by the isotropic rotation diffusion model. The SLE-L approach present is applicable both in the perturbation (Redfield) regime as well as outside the perturbation regime, in the so called slow motion regime.     

Does ℏ play a role in multidimensional spectroscopy? Reduced hierarchy equations of motion approach to molecular vibrations

To investigate the role of quantum effects in vibrational spectroscopies, we have carried out numerically exact calculations of linear and nonlinear response functions for an anharmonic potential system nonlinearly coupled to a harmonic oscillator bath. Although one cannot carry out the quantum calculations of the response functions with full molecular dynamics (MD) simulations for a realistic system which consists of many molecules, it is possible to grasp the essence of the quantum effects on the vibrational spectra by employing a model Hamiltonian that describes an intra- or intermolecular vibrational motion in a condensed phase. The present model fully includes vibrational relaxation, while the stochastic model often used to simulate infrared spectra does not. We have employed the reduced quantum hierarchy equations of motion approach in the Wigner space representation to deal with nonperturbative, non-Markovian, and nonsecular system-bath interactions. Taking the classical limit of the hierarchy equations of motion, we have obtained the classical equations of motion that describe the classical dynamics under the same physical conditions as in the quantum case. By comparing the classical and quantum mechanically calculated linear and multidimensional spectra, we found that the profiles of spectra for a fast modulation case were similar, but different for a slow modulation case. In both the classical and quantum cases, we identified the resonant oscillation peak in the spectra, but the quantum peak shifted to the red compared with the classical one if the potential is anharmonic. The prominent quantum effect is the 1-2 transition peak, which appears only in the quantum mechanically calculated spectra as a result of anharmonicity in the potential or nonlinearity of the system-bath coupling. While the contribution of the 1-2 transition is negligible in the fast modulation case, it becomes important in the slow modulation case as long as the amplitude of the frequency fluctuation is small. Thus, we observed a distinct difference between the classical and quantum mechanically calculated multidimensional spectra in the slow modulation case where spectral diffusion plays a role. This fact indicates that one may not reproduce the experimentally obtained multidimensional spectrum for high-frequency vibrational modes based on classical molecular dynamics simulations if the modulation that arises from surrounding molecules is weak and slow. A practical way to overcome the difference between the classical and quantum simulations was discussed.     

Stochastic quantum molecular dynamics for finite and extended systems

We present a detailed account of the technical aspects of stochastic quantum molecular dynamics, an approach introduced recently by the authors [H. Appel, M. Di Ventra, Phys. Rev. B 80 (2009) 212303] to describe coupled electron-ion dynamics in open quantum systems. As example applications of the method we consider both finite systems with and without ionic motion, as well as describe its applicability to extended systems in the limit of classical ions. The latter formulation allows the study of important phenomena such as decoherence and energy relaxation in bulk systems and surfaces in the presence of time-dependent fields.     

Limitations of differential electrophoresis for measuring colloidal forces: a Brownian dynamics study

Differential electrophoresis experiments are often used to measure subpiconewton forces between two spheres of a heterodoublet. The experiments have been interpreted by solving the electrokinetic equations to obtain a simple Stokes law-type equation. However, for nanocolloids, the effects of Brownian motion alter the interpretation: (1) Brownian translation changes the rate of axial separation. (2) Brownian rotation reduces the alignment of the doublet with the applied electric field. (3) Particles can reaggregate by Brownian motion after they break, forming either heterodoublets or homodoublets, and because homodoublets cannot be broken by differential electrophoresis, this effectively terminates the experiment. We tackle points 1 and 2 using Brownian dynamics simulations (BDS) with electrophoresis as an external force, accounting for convective translation and rotation as well as Brownian translation and rotation. Our simulations identify the lower particle size limit of differential electrophoresis to be about 1 microm for desired statistical accuracy. Furthermore, our simulations predict that particles around 10 nm in size and at ambient conditions will break primarily by Brownian motion, with a negligible effect due to the electric field.     

Chemical reaction dynamics within anisotropic solvents in time-dependent fields

The dynamics of low-dimensional Brownian particles coupled to time-dependent driven anisotropic heavy particles (mesogens) in a uniform bath (solvent) have been described through the use of a variant of the stochastic Langevin equation. The rotational motion of the mesogens is assumed to follow the motion of an external driving field in the linear response limit. Reaction dynamics have also been probed using a two-state model for the Brownian particles. Analytical expressions for diffusion and reaction rates have been developed and are found to be in good agreement with numerical calculations. When the external field driving the mesogens is held at constant rotational frequency, the model for reaction dynamics predicts that the applied field frequency can be used to control the product composition.     

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.     

Quantum generalized Langevin equation approach to gas/solid collisions

The classical generalized Langevin equation (GLE) approach to gas/solid collisions is generalized to quantum scattering. Using Feynman's method of partial path integration, the full gas/solid propagator is reduced to a form in which only the dynamics of the incident atom and the surface oscillator(s) directly struck appear explicitly. Solving this effective dynamical problem in the semiclassical limit yields a stationary phase equation of motion identical in form to the classical GLE. The noise, however, is distributed according to quantum rather than classical statistics. From the GLE a quantum phase can be constructed and an S-matrix computed. The resulting theory is capable of describing inelastic-diffractive scattering which has been seen experimentally by Williams.     

Anomalous diffusion of particles with inertia in external potentials

Recently a new type of Kramers-Fokker-Planck equation has been proposed [Friedrich, R.; Jenko, F.; Baule, A.; Eule, S. Phys. Rev. Lett. 2006, 96, 230601] describing anomalous diffusion in external potentials. In the present paper, the explicit cases of a harmonic potential and a velocity-dependent damping are incorporated. Exact relations for moments for these cases are presented, and the asymptotic behavior for long times is discussed. Interestingly, the bounding potential and the additional damping by itself lead to a subdiffusive behavior. While acting together, the particle becomes localized for long periods of time.     

Incorporation of quantum effects for selected degrees of freedom into the trajectory-based dynamics using spatial domains

The approach of defining quantum corrections on nuclear dynamics of molecular systems incorporated approximately into selected degrees of freedom, is described. The approach is based on the Madelung-de-Broglie-Bohm formulation of time-dependent quantum mechanics which represents a wavefunction in terms of an ensemble of trajectories. The trajectories follow classical laws of motion except that the quantum potential, dependent on the wavefunction amplitude and its derivatives, is added to the external, classical potential. In this framework the quantum potential, determined approximately for practical reasons, is included only into the "quantum" degrees of freedom describing light particles such as protons, while neglecting with the quantum force for the heavy, nearly classical nuclei. The entire system comprised of light and heavy particles is described by a single wavefunction of full dimensionality. The coordinate space of heavy particles is divided into spatial domains or subspaces. The quantum force acting on the light particles is determined for each domain of similar configurations of the heavy nuclei. This approach effectively introduces parametric dependence of the reduced dimensionality quantum force, on classical degrees of freedom. This strategy improves accuracy of the quantum force and does not restrict interaction between the domains. The concept is illustrated for two-dimensional scattering systems, where the quantum force is required to reproduce vibrational energy of the quantum degree of freedom.     

Transient response of a Brownian particle with general damping

We study the transient response of a Brownian particle with general damping in a system of metastable potential well. The escape rate is evaluated as a function of time after an infinite wall is removed from the potential barrier. It takes a relaxation time for the rate to reach its limit value and this rate relaxation time differs from the relaxation time of the majority of the probability around the bottom of the potential well. The rate relaxation time is found to depend on the temperature as well as the damping constant. It involves the diffusion time and the instanton time, in general agreement with recent studies of the overdamped case by Bier et al. [Phys. Rev. E 59, 6422 (1999)].     

Equation of state of nitrogen (N2) at high pressures and high temperatures: molecular dynamics simulation

Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures above 800 K was studied by molecular dynamics (MD) simulations. The dynamics of the N(2) molecules is treated in hard rotor approximation, i.e., it accounts both translational and rotational degrees of freedom. The rotational motion of the N(2) molecule is treated assuming constant moment of inertia of the nitrogen molecule. The new MD program fully accounts anisotropic molecular nitrogen interaction. The N(2)-N(2) interaction potential has been derived by van der Avoird et al. [J. Chem. Phys. 84, 1629 (1986)] using the results of high precision Hartree-Fock ab initio quantum mechanical calculations. The potential, fully accounts rotational symmetry of the N(2)-N(2) system, by employing 6-j Wigner symbols, i.e., preserving full rotational symmetry of the system. Various numerical algorithms were tested, in order to achieve the energy preservation during the simulation. It has been demonstrated that the standard Verlet algorithm was not preserving the energy for the standard MD time step, equal to 5x10(-16) s. Runge-Kutta fourth order method was able to preserve the energy within 10(-4) relative error, but it requires calculation of the force four times for each time step and therefore it is highly inefficient. A predictor-corrector method of the fifth order (PC5) was found to be efficient and precise and was therefore adopted for the simulation of the molecular nitrogen properties at high pressure. Singer and Fincham algorithms were tested and were found to be as precise as PC5 algorithm and they were also used in the simulation of the equation of state. Results of MD simulations are in very good agreement with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars). For higher pressures, up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 85, 1262 (1984)], the obtained numerical results provide new data of the experimentally unexplored region. These data were formulated in the analytical form of pressure-density-temperature equation of state.