Search

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

An accurate system-bath model to investigate the quantum dynamics of hydrogen atoms chemisorbed on graphene is presented. The system comprises a hydrogen atom and the carbon atom from graphene that forms the covalent bond, and it is described by a previously developed 4D potential energy surface based on density functional theory ab initio data. The bath describes the rest of the carbon lattice and is obtained from an empirical force field through inversion of a classical equilibrium correlation function describing the hydrogen motion. By construction, model building easily accommodates improvements coming from the use of higher level electronic structure theory for the system. Further, it is well suited to a determination of the system-environment coupling by means of ab initio molecular dynamics. This paper details the system-bath modeling and shows its application to the quantum dynamics of vibrational relaxation of a chemisorbed hydrogen atom, which is here investigated at T = 0 K with the help of the multi-configuration time-dependent Hartree method. Paper II deals with the sticking dynamics. [ABSTRACT FROM AUTHOR]

In this paper, quantum fluctuations of a carbonyl sulfide molecule in helium-4 clusters are studied as a function of cluster size N in a small-to-large size regime (2≤=N≤=64). The molecular rotation of the dopant shows nonmonotonic size dependence in the range of 10≤=N≤=20, reflecting the density distribution of the helium atoms around the molecule. The size dependence on the rotational constant shows a plateau for N>=20, which is larger than the experimental nanodroplet value. Superfluid response of the doped cluster is found to show remarkable anisotropy especially for N≤=20. The superfluid fraction regarding the axis perpendicular to the molecular axis shows a steep increase at N=10, giving the significant enhancement of the rotational fluctuation of the molecule. On the other hand, the superfluid fraction regarding the axis parallel to the molecular axis reaches 0.9 at N=5, arising from the bosonic exchange cycles of the helium atoms around the molecular axis. The anisotropy in the superfluid response is found to be the direct consequence of the configurations of the bosonic exchange cycles. [ABSTRACT FROM AUTHOR]

We report quantum-mechanical calculations on the formation of aligned rotational wave packets in the ground state of I[sub 2], following the excitation of a rotationally cold sample of molecules with an intense picosecond laser pulse. Using these calculations, we have studied (1) the dynamic alignment during or shortly after the intense pump pulse, (2) the alignment at characteristic revival times following the laser excitation, and (3) the alignment between the revivals, as a function of the laser-pulse duration and peak intensity and the initial rotational temperature of the molecular beam. We conclude that the alignment at a revival is maximum for intermediate pulse durations (∼3 ps for I[sub 2]) that are long enough to get efficient population transfer out of the initially populated state, yet short enough for the interaction to remain partially diabatic and partially adiabatic. The alignment at the revivals increases with laser intensity, although the time-dependent structure of the revivals becomes increasingly complicated as the laser intensity is increased, and the maximum degree of alignment that is obtained is eventually limited by the initial rotational temperature of the molecular beam. The findings in this paper qualitatively agree with experimental results on short-pulse laser-induced alignment that were presented in our experimental article (preceding paper). © 2002 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Advances in experimental methodology aligned with technological developments, such as 3rd generation light sources, X-ray Free Electron Lasers, and High Harmonic Generation, have led to a paradigm shift in the capability of X-ray spectroscopy to deliver high temporal and spectral resolution on an extremely broad range of samples in a wide array of different environments. Importantly, the complex nature and high information content of this class of techniques mean that detailed theoretical studies are often essential to provide a firm link between the spectroscopic observables and the underlying molecular structure and dynamics. In this paper, we present approaches for simulating dynamical processes in X-ray spectroscopy based upon on-the-fly quantum dynamics with a Gaussian basis set. We show that it is possible to provide a fully quantum description of X-ray spectra without the need of precomputing highly multidimensional potential energy surfaces. It is applied to study two different dynamical situations, namely, the core-hole lifetime dynamics of the water monomer and the dissociation of C F 4 + recently studied using pump-probe X-ray spectroscopy. Our results compare favourably to previous experiments, while reducing the computational effort, providing the scope to apply them to larger systems. [ABSTRACT FROM AUTHOR]

In this paper, we examine how and when quantum evolution can be approximated in terms of (generalized) classical dynamics in calculations of correlation functions, with a focus on the symmetrized time correlation function introduced by Schofield. To that end, this function is expressed as a path integral in complex time and written in terms of sum and difference path variables. Taylor series expansion of the path integral's exponent to first and second order in the difference variables leads to two original developments. The first order expansion is used to obtain a simple, path integral based, derivation of the so-called Schofield's quantum correction factor. The second order result is employed to show how quantum mechanical delocalization manifests itself in the approximation of the correlation function and hinders, even in the semiclassical limit, the interpretation of the propagators in terms of sets of guiding classical trajectories dressed with appropriate weights. [ABSTRACT FROM AUTHOR]

Gas-phase atomic anions lack bound electronic excited states, yet in solution many of these anions exhibit intense absorption bands due to the presence of excited states, referred to as charge-transfer-to-solvent (CTTS) states that are bound only by the presence of the solvent. CTTS spectra thus serve as delicate probes of solute-solvent interactions, but the fact that they are created by the interactions of a solute with many solvent molecules makes them a challenge to describe theoretically. In this paper, we use mixed quantum/classical molecular dynamics with the two-electron Fourier-grid (2EFG) electronic structure method presented in the previous paper [W. J. Glover, R. E. Larsen, and B. J. Schwartz, J. Chem. Phys. 132, 144101 (2010)] to simulate the CTTS states of a sodium anion in liquid tetrahydrofuran, Na-/THF. Since our 2EFG method is based on configuration interaction with single and double excitations in a grid basis, it allows for an exact treatment of the two valence electrons of the sodium anion. To simulate Na-/THF, we first develop a new electron–THF pseudopotential, and we verify the accuracy of this potential by reproducing the experimental absorption spectrum of an excess electron in liquid THF with near quantitative accuracy. We also are able to reproduce the CTTS spectrum of Na-/THF and find that the CTTS states of Na- exhibit a Rydberg-like progression due to the pre-existing long-range solvent polarization around the anion. We also find that the CTTS states are highly mixed with the disjoint electronic states supported by naturally occurring solvent cavities that exist in liquid THF. This mixing explains why the solvated electrons that are ejected following CTTS excitation appear with their equilibrium absorption spectrum. The mixing of the CTTS and solvent-cavity states also explains why the recombination of the electron and its geminate Na0 partner occurs on slower time scales when photoexciting in the blue compared to in the red of the CTTS band: blue excitation accesses CTTS states whose charge densities lies further from the position of the anion, whereas red excitation accesses CTTS states that lie primarily within the anion’s first solvation shell. Finally, we see that the radial character of the CTTS states near the Na+ core matches that of Na0, explaining why the spectrum of this species appears instantly after photoexciting Na-. [ABSTRACT FROM AUTHOR]

A numerical solution for the quantum dynamics of the spin boson problem is obtained using the semiclassical initial value series representation approach to the quantum dynamics. The zeroth order term of the series is computed using the new forward-backward representation for correlation functions presented in the preceding adjacent paper. This leads to a rapid convergence of the Monte Carlo sampling, as compared to previous attempts. The zeroth order results are already quite accurate. The first order term of the series is small, demonstrating the rapid convergence of the semiclassical initial value representation series. This is the first time that the first order term in the semiclassical initial value representation series has been converged for systems with the order of 50 degrees of freedom. [ABSTRACT FROM AUTHOR]

In this paper, we present a path integral hybrid Monte Carlo (PIHMC) method for rotating molecules in quantum fluids. This is an extension of our PIHMC for correlated Bose fluids [S. Miura and J. Tanaka, J. Chem. Phys. 120, 2160 (2004)] to handle the molecular rotation quantum mechanically. A novel technique referred to be an effective potential of quantum rotation is introduced to incorporate the rotational degree of freedom in the path integral molecular dynamics or hybrid Monte Carlo algorithm. For a permutation move to satisfy Bose statistics, we devise a multilevel Metropolis method combined with a configurational-bias technique for efficiently sampling the permutation and the associated atomic coordinates. Then, we have applied the PIHMC to a helium-4 cluster doped with a carbonyl sulfide molecule. The effects of the quantum rotation on the solvation structure and energetics were examined. Translational and rotational fluctuations of the dopant in the superfluid cluster were also analyzed. [ABSTRACT FROM AUTHOR]

In the hidden variable theory, Bohm proved a connection between the Schrödinger and Hamilton–Jacobi equations and showed the existence of classical paths, for which the generalized Bohr quantization condition is valid. In this paper we prove similar properties, starting from the equivalence between the Schrödinger and wave equations in the case of the conservative bound systems. Our approach is based on the equations and postulates of quantum mechanics without using any additional postulate. Like in the hidden variable theory, the above properties are proven without using the approximation of geometrical optics or the semiclassical approximation. Since the classical paths have only a mathematical significance in our analysis, our approach is consistent with the postulates of quantum mechanics. [ABSTRACT FROM AUTHOR]

Expanding the set of stable, accurate, and scalable methods for simulating molecular quantum dynamics is important for accelerating the computational exploration of molecular processes. In this paper, we adapt the signed particles Monte Carlo algorithm for solving the transient Wigner equation to scenarios of chemical interest. This approach was used in the past to study electronic processes in semi-conductors, but to the best of our knowledge, it had never been applied to molecular modeling. We present the algorithm and demonstrate its excellent performance on harmonic and double well potentials for electronic and nuclear systems. We explore the stability of the algorithm, discuss the choice of hyper-parameters, and cautiously speculate that it may be used in quantum molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

The formulation of quantum statistical mechanics based on the path centroid variable in Feynman path integration is generalized to a phase space perspective, thereby including the momentum as an independent dynamical variable. By virtue of this approach, operator averages and imaginary time correlation functions can be expressed in terms of an averaging over the multidimensional phase space centroid density. The imaginary time centroid-constrained correlation function matrix for the phase space variables is then found to define the effective thermal width of the phase space centroid variable. These developments also make it possible to rigorously analyze the centroid molecular dynamics method for computing quantum dynamical time correlation functions. As a result, the centroid time correlation function as calculated from centroid molecular dynamics is shown to be a well-defined approximation to the exact Kubo transformed position correlation function. This analysis thereby clarifies the underlying role of the equilibrium path centroid variable in the quantum dynamical position correlation function and provides a sound theoretical basis for the centroid molecular dynamics method. [ABSTRACT FROM AUTHOR]

In this work, we examine the efficiency of the recently developed [P. Piecuch et al., J. Chem. Phys. 99, 6732 (1993)] state-selective (SS) multi-reference (MR) coupled-cluster (CC) method for calculation of molecular properties. In our earlier papers, we demonstrated that the SSMRCC method with inclusion of single, double, and internal and semi-internal triple excitations [SSCCSD(T) approach] is capable of providing an accurate description of the ground-state potential energy surfaces. In this paper, we present the dipole moment and polarizability values of the HF molecule at equilibrium and stretched geometries calculated using finite field technique and SSCCSD(T) ansatz. The calculations use double zeta quality basis sets with and without polarization functions. Molecular orbital basis sets include both relaxed and nonrelaxed orbitals. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Classical transition state theory (TST) provides the rigorous basis for the application of molecular dynamics (MD) to infrequent events, i.e., reactions that are slow due to a high energy barrier. The TST rate is simply the equilibrium flux through a surface that divides reactants from products. In order to apply MD to infrequent events, corrections to the TST rate that account for recrossings of the dividing surface are computed by starting trajectories at the dividing surface and integrating them backward and forward in time. Both classical TST and conventional MD invoke the adiabatic approximation, i.e., the assumption that nuclear motion evolves on a single potential energy surface. Many chemical rate processes involve multiple potential energy surfaces, however, and a number of ‘‘surface-hopping’’ MD methods have been developed in order to incorporate nonadiabatic transitions among the potential energy surfaces. In this paper we generalize TST to processes involving multiple potential energy surfaces. This provides the framework for a new method for MD simulation of infrequent events for reactions that evolve on multiple potential energy surfaces. We show how this method can be applied rigorously even in conjunction with phase-coherent surface-hopping methods, where the probability of switching potential energy surfaces depends on the history of the trajectory, so integrating trajectories backward to calculate the recrossing correction is problematic. We illustrate this new method by applying it in conjunction with the ‘‘molecular dynamics with quantum transitions’’ (MDQT) surface-hopping method to a one-dimensional two-state barrier crossing problem. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Trajectory-based Gaussian basis sets have been tremendously successful in describing highdimensional quantum molecular dynamics. In this paper, we introduce a pseudospectral Gaussianbased method that achieves accurate quantum dynamics using efficient, real-space sampling of the time-dependent basis set. As in other Gaussian basis methods, we begin with a basis set expansion using time-dependent Gaussian basis functions guided by classical mechanics. Unlike other Gaussian methods but characteristic of the pseudospectral and collocation methods, the basis set is tested with N Dirac delta functions, where N is the number of basis functions, rather than using the basis function as test functions. As a result, the integration for matrix elements is reduced to function evaluation. Pseudospectral Gaussian dynamics only requires O(N) potential energy calculations, in contrast to O(N²) evaluations in a variational calculation. The classical trajectories allow small basis sets to sample high-dimensional potentials. Applications are made to diatomic oscillations in a Morse potential and a generalized version of the Henon-Heiles potential in two, four, and six dimensions. Comparisons are drawn to full analytical evaluation of potential energy integrals (variational) and the bra-ket averaged Taylor (BAT) expansion, an O(N) approximation used in Gaussian-based dynamics. In all cases, the pseudospectral Gaussian method is competitive with full variational calculations that require a global, analytical, and integrable potential energy surface. Additionally, the BAT breaks down when quantum mechanical coherence is particularly strong (i.e., barrier reflection in the Morse oscillator). The ability to obtain variational accuracy using only the potential energy at discrete points makes the pseudospectral Gaussian method a promising avenue for on-the-fly dynamics, where electronic structure calculations become computationally significant. [ABSTRACT FROM AUTHOR]

Basing on the earlier works on the hierarchical equations of motion for quantum transport, we present in this paper a first principles scheme for time-dependent quantum transport by combining time-dependent density functional theory (TDDFT) and Keldysh's non-equilibrium Green's function formalism. This scheme is beyond the wide band limit approximation and is directly applicable to the case of non-orthogonal basis without the need of basis transformation. The overlap between the basis in the lead and the device region is treated properly by including it in the self-energy and it can be shown that this approach is equivalent to a lead-device orthogonalization. This scheme has been implemented at both TDDFT and density functional tight-binding level. Simulation results are presented to demonstrate our method and comparison with wide band limit approximation is made. Finally, the sparsity of the matrices and computational complexity of this method are analyzed. [ABSTRACT FROM AUTHOR]

Most of today's molecular-dynamics simulations of materials are based on the Born-Oppenheimer approximation. There are many cases, however, in which the coupling of the electrons and nuclei is important and it is necessary to go beyond the Born-Oppenheimer approximation. In these methods, the non-adiabatic coupling vectors are fundamental since they represent the link between the classical atomic motion of the nuclei and the time evolution of the quantum electronic state. In this paper we analyze the calculation of non-adiabatic coupling vectors in a basis set of local orbitals and derive an expression to calculate them in a practical and computationally efficient way. Some examples of the application of this expression using a local-orbital density functional theory approach are presented for a few simple molecules: H3, formaldimine, and azobenzene. These results show that the approach presented here, using the Slater transition-state density, is a very promising way for the practical calculation of non-adiabatic coupling vectors for large systems. [ABSTRACT FROM AUTHOR]

This paper shows how a compact finite difference Hessian approximation scheme can be proficiently implemented into semiclassical initial value representation molecular dynamics. Effects of the approximation on the monodromy matrix calculation are tested by propagating initial sampling distributions to determine power spectra for analytic potential energy surfaces and for 'on the fly' carbon dioxide direct dynamics. With the approximation scheme the computational cost is significantly reduced, making ab initio direct semiclassical dynamics computationally more feasible and, at the same time, properly reproducing important quantum effects inherent in the monodromy matrix and the pre-exponential factor of the semiclassical propagator. [ABSTRACT FROM AUTHOR]

An efficient time-dependent importance sampling method is developed for the Monte Carlo calculation of time correlation functions via the initial value representation (IVR) of semiclassical (SC) theory. A prefactor-free time-dependent sampling function weights the importance of a trajectory based on the magnitude of its contribution to the time correlation function, and global trial moves are used to facilitate the efficient sampling the phase space of initial conditions. The method can be generally applied to sampling rare events efficiently while avoiding being trapped in a local region of the phase space. Results presented in the paper for two system-bath models demonstrate the efficiency of this new importance sampling method for full SC-IVR calculations. [ABSTRACT FROM AUTHOR]

Internal conversion is an inherently quantum mechanical process. To date, 'on the fly' computation of internal conversion rates is limited to harmonic approximations, which would seem to be especially unsuitable, given that the typical transition to the ground electronic state occurs at energies which are far from the harmonic limit. It is thus of interest to study the applicability of the semiclassial initial value representation (SCIVR) approach which is in principle amenable to on the fly studies even with 'many' degrees of freedom. In this paper we study the applicability of the Herman-Kluk (HK) SCIVR to a model system with two coupled and anharmonic degrees of freedom. We find that (a) the HK SCIVR is a good approximation to the exact quantum dynamics; (b) computation of the first order correction to the HK-SCIVR approximation corroborates the accuracy; (c) by studying a large parameter range, we find that the harmonic approximation is mostly unsatisfactory; and (d) for the specific model used, the coupling between the modes was found to be relatively unimportant. These results imply that the HK-SCIVR methodology is a good candidate for on the fly studies of internal conversion processes of 'large' molecules. [ABSTRACT FROM AUTHOR]

Intrapulse quantum cascade laser spectrometers are able to produce both saturation and molecular alignment of the gas sample. This is due to the rapid sweep of the radiation through the absorption features. The intrapulse time domain spectra closely resemble those recorded in coherent optical nutation experiments. In the present paper, the frequency down-chirped technique is employed to investigate the nitrous oxide–foreign gas collisions. We have demonstrated that the measurements may be characterized by the induced polarization dominated and collision dominated measurement limits. The first of these is directly related to the time dependence of the long range collision cross sections. Among the collisional partners considered, carbon dioxide shows a very unusual behavior of rapid polarization damping, resulting in the production of symmetrical line shapes at very low gas buffer pressures. In the collision dominated regime, the pressure broadening parameters, which we have derived, are comparable at slow chirp rates, with those derived from other experimental methods. By comparing the pressure broadening coefficients of Ar, N2, and CO2 with those of He, making use of the chirp rate independence of the pressure broadening by helium, we have shown that at higher chirp rates there is clear evidence of the chirp-rate dependence of the pressure broadening parameters of N2 and CO2. [ABSTRACT FROM AUTHOR]

Initial value representations (IVRs) of semiclassical (SC) theory provide a general approach for adding quantum mechanical effects to classical molecular dynamics simulations of large molecular systems. Of the various versions of SC-IVR methodology for evaluating time correlation functions, the Fourier transform forward-backward (FB) approach is the simplest one that is able to describe true quantum coherence effects, so it is of considerable importance to find efficient and systematic ways for implementing it. It is shown in this paper that a Gaussian approximation for the “structure function”—the dependence of the correlation function on the (typically) momentum jump parameter—provides an efficient and accurate way for doing so. The approach is illustrated by an application to the time-dependent radial distribution function of I2 (after photoexcitation) in a cluster of (up to 16) argon atoms. [ABSTRACT FROM AUTHOR]

The important role of liquid water in many areas of science from chemistry, physics, biology, geology to climate research, etc., has motivated numerous theoretical studies of its structure and dynamics. The significance of quantum effects on the properties of water, however, has not yet been fully resolved. In this paper we focus on quantum dynamical effects in liquid water based on the linearized semiclassical initial value representation (LSC-IVR) with a quantum version of the simple point charge/flexible (q-SPC/fw) model [Paesani et al., J. Chem. Phys. 125, 184507 (2006)] for the potential energy function. The infrared (IR) absorption spectrum and the translational diffusion constants have been obtained from the corresponding thermal correlation functions, and the effects of intermolecular and intramolecular correlations have been studied. The LSC-IVR simulation results are compared with those predicted by the centroid molecular dynamics (CMD) approach. Although the LSC-IVR and CMD results agree well for the broadband for hindered motions in liquid water, the intramolecular bending and O–H stretching peaks predicted by the LSC-IVR are blueshifted from those given by CMD; reasons for this are discussed. We also suggest that the broadband in the IR spectrum corresponding to restricted translation and libration gives more information than the diffusion constant on the nature of quantum effects on translational and rotational motions and should thus receive more attention in this regard. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations are performed on two sodium chloride solutions in TIP4P water with concentrations c=1.36 mol/kg and c=2.10 mol/kg upon supercooling. The isotherms and isochores planes are calculated. The temperature of maximum density line and the limit of mechanical stability line are obtained from the analysis of the thermodynamic planes. The comparison of the results shows that for densities well above the limit of mechanical stability, the isotherms and isochores of the sodium chloride aqueous solution shift to lower pressures upon increasing concentration while the limit of mechanical stability is very similar to that of bulk water for both concentrations. We also find that the temperature of maximum density line shifts to lower pressures and temperatures upon increasing concentration. Indications of the presence of a liquid-liquid coexistence are found for both concentrations. [ABSTRACT FROM AUTHOR]

A new method called adaptive force matching (AFM) has been developed that is capable of producing high quality force fields for condensed phase simulations. This procedure involves the parametrization of force fields to reproduce ab initio forces obtained from condensed phase quantum-mechanics/molecular-mechanics (QM/MM) calculations. During the procedure, the MM part of the QM/MM is iteratively improved so as to approach ab initio quality. In this work, the AFM method has been tested to parametrize force fields for liquid water so that the resulting force fields reproduce forces calculated using the ab initio MP2 and the Kohn–Sham density functional theory with the Becke–Lee–Yang–Parr (BLYP) and Becke three-parameter LYP (B3LYP) exchange correlation functionals. The AFM force fields generated in this work are very simple to evaluate and are supported by most molecular dynamics (MD) codes. At the same time, the quality of the forces predicted by the AFM force fields rivals that of very expensive ab initio calculations and are found to successfully reproduce many experimental properties. The site-site radial distribution functions (RDFs) obtained from MD simulations using the force field generated from the BLYP functional through AFM compare favorably with the previously published RDFs from Car–Parrinello MD simulations with the same functional. Technical aspects of AFM such as the optimal QM cluster size, optimal basis set, and optimal QM method to be used with the AFM procedure are discussed in this paper. [ABSTRACT FROM AUTHOR]

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed α-helix and extended β and PII conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the PII conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is αR. Polarizable CHARMM force field gives almost exclusively PII conformation and other force fields predict that both α-helical and extended (β and PII) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the PII conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper. [ABSTRACT FROM AUTHOR]

This paper describes the application of a relatively simple, but computationally tractable, “particle-on-a-sphere” (POS) model for quantum-mechanical calculation of large-amplitude, H atom dynamics in polyatomic hydrides (XHn), based on radially relaxed, two-dimensional angular motion of H atoms on the surface of a sphere. This work focuses on systems with many degrees of freedom, i.e., XH4 (eight dimensional) and XH5 (ten dimensional), with corresponding molecular analogs of CH4 and CH5+ and is applicable to rovibrationally excited states with J>=0. A pairwise-additive potential fit for CH5+, which yields remarkable agreement with geometries, energies, and barrier heights on the full-dimensional surface of Brown et al. [J. Chem. Phys. 121, 4105 (2004)] is presented. Comparisons with experimental data and diffusion quantum Monte Carlo (DMC) methods test convergence for the POS model and provide insight into multidimensional quantum rovibrational dynamics. In particular, POS energy-level patterns for a series of scaled CH5+ potentials indicate an absence of strong tunneling behavior, consistent with the highly delocalized wave functions, large zero-point energies, and small interconversion barriers noted in previous DMC studies of Brown et al. [ABSTRACT FROM AUTHOR]

Many naturally occurring fluids, such as crude oils, consist of a very large number of components. It is often of interest to determine the composition of the fluids in situ. Diffusion coefficients and nuclear magnetic resonance (NMR) relaxation times can be measured in situ and depend on the size of the molecules. It has been shown [D. E. Freed et al., Phys. Rev. Lett. 94, 067602 (2005)] that the diffusion coefficient of each component in a mixture of alkanes follows a scaling law in the chain length of that molecule and in the mean chain length of the mixture, and these relations were used to determine the chain length distribution of crude oils from NMR diffusion measurements. In this paper, the behavior of NMR relaxation times in mixtures of chain molecules is addressed. The author explains why one would expect scaling laws for the transverse and longitudinal relaxation times of mixtures of short chain molecules and mixtures of alkanes, in particular. It is shown how the power law dependence on the chain length can be calculated from the scaling laws for the translational diffusion coefficients. The author fits the literature data for NMR relaxation in binary mixtures of alkanes and finds that its dependence on chain length agrees with the theory. Lastly, it is shown how the scaling laws in the chain length and the mean chain length can be used to determine the chain length distribution in crude oils that are high in saturates. A good fit is obtained between the NMR-derived chain length distributions and the ones from gas chromatography. [ABSTRACT FROM AUTHOR]

An approach is developed in the replica exchange framework to enhance conformational sampling for the quantum mechanical (QM) potential based molecular dynamics simulations. Importantly, with our enhanced sampling treatment, a decent convergence for electronic structure self-consistent-field calculation is robustly guaranteed, which is made possible in our replica exchange design by avoiding direct structure exchanges between the QM-related replicas and the activated (scaled by low scaling parameters or treated with high “effective temperatures”) molecular mechanical (MM) replicas. Although the present approach represents one of the early efforts in the enhanced sampling developments specifically for quantum mechanical potentials, the QM-based simulations treated with the present technique can possess the similar sampling efficiency to the MM based simulations treated with the Hamiltonian replica exchange method (HREM). In the present paper, by combining this sampling method with one of our recent developments (the dual-topology alchemical HREM approach), we also introduce a method for the sampling enhanced QM-based free energy calculations. [ABSTRACT FROM AUTHOR]

Trajectories of duration τ joining two points q0 and q1 in the configuration space of a classical system satisfy Hamilton’s principle: they are stationary points of the classical action. The second variation (fluctuations) of the action around the stationary point signals whether the latter is or not a minimum and delivers the density in trajectory space around the points q0 and q1. This concept of density of paths is of great importance in semiclassical quantum theory, where it weights the contribution to the propagator from the single classical trajectories. In this paper, two algorithms based on the concepts of molecular dynamics simulation are introduced for computing the density of paths, also called van Vleck [Proc. Natl. Acad. Sci. U.S.A. 14, 178 (1928)] determinant. Examples for realistic systems are presented, together with a suggestion about possible applications in the field of rare events in physics and chemistry. [ABSTRACT FROM AUTHOR]

This paper presents a general theoretical formulation for proton-coupled electron transfer (PCET) reactions. The solute is represented by a multistate valence bond model, and the active electrons and transferring proton(s) are treated quantum mechanically. This formulation enables the classical or quantum mechanical treatment of the proton donor-acceptor vibrational mode, as well as the dynamical treatment of the proton donor-acceptor mode and the solvent. Nonadiabatic rate expressions are presented for PCET reactions in a number of well-defined limits for both dielectric continuum and molecular representations of the environment. The dynamical rate expressions account for correlations between the fluctuations of the proton donor-acceptor distance and the nonadiabatic PCET coupling. The quantities in the rate expressions can be calculated with a dielectric continuum model or a molecular dynamics simulation of the full system. The significance of the quantum and dynamical effects of the proton donor-acceptor mode is illustrated with applications to model PCET systems.© 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The three-dimensional hydrogen cyanide/isocyanide isomerization problem is taken as an example to present a general theory for computing the phase space structures which govern classical reaction dynamics in systems with an arbitrary (finite) number of degrees of freedom. The theory, which is algorithmic in nature, comprises the construction of a dividing surface of minimal flux which is locally a ''surface of no return.'' The theory also allows for the computation of the global phase space transition pathways that trajectories must follow in order to react. The latter are enclosed by the stable and unstable manifolds of a so-called normally hyperbolic invariant manifold (NHIM). A detailed description of the geometrical structures and the resulting constraints on reaction dynamics is given, with particular emphasis on the three degrees of freedom case. A procedure is given which uses these structures to compute orbits homoclinic to, and heteroclinic between, NHIMs. The role of homoclinic and heteroclinic orbits in global recrossings of dividing surfaces and transport in complex systems is explained. The complete description provided here is inherently one within phase space; it cannot be inferred from a configuration space picture. A complexification of the classical phase space structures to incorporate quantum effects is also discussed. The results presented here call into question certain assumptions routinely made on the global dynamics; this paper provides methods that enable one to understand and quantify the phase space dynamics of reactions without making such assumptions. [ABSTRACT FROM AUTHOR]

In this paper we describe an application of the trajectory-based semiclassical Liouville method for modeling coherent molecular dynamics on multiple electronic surfaces to the treatment of the evolution and decay of quantum electronic coherence in many-body systems. We consider a model representing the coherent evolution of quantum wave packets on two excited electronic surfaces of a diatomic molecule in the gas phase and in rare gas solvent environments, ranging from small clusters to a cryogenic solid. For the gas phase system, the semiclassical trajectory method is shown to reproduce the evolution of the electronic–nuclear coherence nearly quantitatively. The dynamics of decoherence are then investigated for the solvated systems using the semiclassical approach. It is found that, although solvation in general leads to more rapid and extensive loss of quantum coherence, the details of the coupled system–bath dynamics are important, and in some cases the environment can preserve or even enhance quantum coherence beyond that seen in the isolated system. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In the de Broglie–Bohm formulation of quantum mechanics the time-dependent Schrödinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is accurate for semiclassical systems, we use approximate quantum potentials. Recently, we have shown that optimization of the nonclassical component of the momentum operator in terms of fitting functions leads to the energy-conserving approximate quantum potential. In particular, linear fitting functions give the exact time evolution of a Gaussian wave packet in a locally quadratic potential and can describe the dominant quantum-mechanical effects in the semiclassical scattering problems of nuclear dynamics. In this paper we formulate the Bohmian dynamics on subspaces and define the energy-conserving approximate quantum potential in terms of optimized nonclassical momentum, extended to include the domain boundary functions. This generalization allows a better description of the non-Gaussian wave packets and general potentials in terms of simple fitting functions. The optimization is performed independently for each domain and each dimension. For linear fitting functions optimal parameters are expressed in terms of the first and second moments of the trajectory distribution. Examples are given for one-dimensional anharmonic systems and for the collinear hydrogen exchange reaction. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper, we describe a new method for simulating quantum processes using classical-like molecular dynamics. The approach is based on solving the quantum Liouville equation in the Wigner representation using ensembles of classical trajectories in phase space. The nonlocality of quantum mechanics is incorporated in the trajectory representation as nonclassical interactions between the members of the ensemble, leading to an entanglement of their evolution. The statistical independence of the individual trajectories making up an ensemble in the classical limit is lost when quantum effects are included, and the entire state of the system must be propagated as a unified whole. We develop the formalism and its numerical implementation, and illustrate its application on two model problems of quantum mechanical tunneling: escape from a metastable well and wave packet penetration of the Eckhart barrier. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Low-frequency intermolecular dynamics in liquids is studied by ultrafast four- and six-wave mixing. The theory of these nonlinear optical processes is given for electronically nonresonant optical interactions up to fifth order in the electric field. The Born–Oppenheimer approximation is used to separate the motional part of the response functions from coordinate independent electronic hyperpolarizabilities. A large variety of experiments, involving far-infrared absorption, ordinary Rayleigh–Raman or hyper Rayleigh–Raman scattering is covered by this theory. The response in nonresonant six-wave mixing comprises four dynamically different processes. It is shown that one of the terms contains information on the time scale(s) of intermolecular dynamics, that is not available from lower-order nonresonant experiments. For instance, homogenous and inhomogeneous contributions to line broadening can be distinguished. The optical response of harmonic nuclear motion is calculated for nonlinear coordinate dependence of the polarizabilities. Results for level-dependent and level-independent damping of the motion are compared. It is shown that level-dependent damping destroys the interference between different quantum mechanical pathways, yielding an extra contribution to the fifth-order response that has not been discussed before. When two or more nuclear modes determine the optical response, their relative contributions to the four- and six-wave mixing signals are in general different. These contributions are determined by the coordinate dependence of the electronic polarizability, which is usually not fully known. Model calculations are presented for the dynamic parameters of liquid CS2. The theory of this paper will be employed in Part II, to analyze experimental results on femtosecond four- and six-wave mixing. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

This paper presents a mathematical derivation of a model for quantum-classical molecular dynamics (QCMD) as a partial classical limit of the full Schrödinger equation. This limit is achieved in two steps: separation of the full wave function and short wave asymptotics for its ‘‘classical’’ part. Both steps can be rigorously justified under the same smallness assumptions. This throws some light on the time-dependent self-consistent-field method and on mixed quantum-semiclassical models, which also depend on the separation step. On the other hand, the theory leads to a characterization of the critical situations in which the QCMD model is in danger of largely deviating from the solution of full Schrödinger equation. These critical situations are exemplified in an illustrative numerical simulation: the collinear collision of a classical particle with a harmonic quantum oscillator. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

This paper presents a way of improving second-order perturbation theory calculations by summing contributions of uncoupled excitations to infinite order. For problems involving molecular vibrations, the new theory is shown to give similar results to conventional second-order perturbation theory when the system treated has no near resonances but also to give accurate and stable results even very close to resonance. The new theory is tested by comparison to converged variational calculations for vibrational energy levels of formaldehyde, formaldehyde-d2, and two two-dimensional model subsystems based on formaldehyde. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The routine location of the global minimum energy structures of a molecular system remains a challenging problem in chemical physics. Recently Straub and co-workers [J. Ma et al., J. Chem. Phys. 99, 4024 (1993) and 101, 533 (1994)] have introduced a number of novel dynamic simulated annealing algorithms for global optimization. Instead of using Newton’s equations, they are based on Liouville’s equation and the particles are described with phase space distributions. The methods were shown to have significantly enhanced effectiveness over traditional methods when used to locate the global minima of Lennard-Jones clusters. In this paper the extensions necessary to use these methods with ab initio quantum mechanical potentials are presented and their application is illustrated by the simulated annealing of a small lithium atom cluster. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. The proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms. [ABSTRACT FROM AUTHOR]

We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich–Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer–Miller–Stock–Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately. [ABSTRACT FROM AUTHOR]

We study theoretically the quantum dynamics and spectroscopy of rovibrational polaritons formed in a model system composed of a single rovibrating diatomic molecule, which interacts with two degenerate, orthogonally polarized modes of an optical Fabry–Pérot cavity. We employ an effective rovibrational Pauli–Fierz Hamiltonian in length gauge representation and identify three-state vibro-polaritonic conical intersections (VPCIs) between singly excited vibro-polaritonic states in a two-dimensional angular coordinate branching space. The lower and upper vibrational polaritons are of mixed light–matter hybrid character, whereas the intermediate state is purely photonic in nature. The VPCIs provide effective population transfer channels between singly excited vibrational polaritons, which manifest in rich interference patterns in rotational densities. Spectroscopically, three bright singly excited states are identified when an external infrared laser field couples to both a molecular and a cavity mode. The non-trivial VPCI topology manifests as pronounced multi-peak progression in the spectral region of the upper vibrational polariton, which is traced back to the emergence of rovibro-polaritonic light–matter hybrid states. Experimentally, ubiquitous spontaneous emission from cavity modes induces a dissipative reduction of intensity and peak broadening, which mainly influences the purely photonic intermediate state peak as well as the rovibro-polaritonic progression. [ABSTRACT FROM AUTHOR]

We perform path integral molecular dynamics (PIMD) simulations of a monatomic liquid that exhibits a liquid–liquid phase transition and liquid–liquid critical point. PIMD simulations are performed using different values of Planck's constant h, allowing us to study the behavior of the liquid as nuclear quantum effects (NQE, i.e., atoms delocalization) are introduced, from the classical liquid (h = 0) to increasingly quantum liquids (h > 0). By combining the PIMD simulations with the ring-polymer molecular dynamics method, we also explore the dynamics of the classical and quantum liquids. We find that (i) the glass transition temperature of the low-density liquid (LDL) is anomalous, i.e., T g LDL (P) decreases upon compression. Instead, (ii) the glass transition temperature of the high-density liquid (HDL) is normal, i.e., T g HDL (P) increases upon compression. (iii) NQE shift both T g LDL (P) and T g HDL (P) toward lower temperatures, but NQE are more pronounced on HDL. We also study the glass behavior of the ring-polymer systems associated with the quantum liquids studied (via the path-integral formulation of statistical mechanics). There are two glass states in all the systems studied, low-density amorphous ice (LDA) and high-density amorphous ice (HDA), which are the glass counterparts of LDL and HDL. In all cases, the pressure-induced LDA–HDA transformation is sharp, reminiscent of a first-order phase transition. In the low-quantum regime, the LDA–HDA transformation is reversible, with identical LDA forms before compression and after decompression. However, in the high-quantum regime, the atoms become more delocalized in the final LDA than in the initial LDA, raising questions on the reversibility of the LDA–HDA transformation. [ABSTRACT FROM AUTHOR]

Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum–classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum–classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems. [ABSTRACT FROM AUTHOR]

In this paper we present coupled quantum mechanical/molecular dynamics simulations on the solvation of naphthalene in acetonitrile. A semiempirical treatment is introduced into classical molecular dynamics simulations to deal with the electronic polarizability of the solute. We study the structure of the solvation shell, the features of the electrostatic intermolecular interactions, the dynamics of solute and solvent molecules and the spectral shifts of the solute. For this particular system it is found that the influence of induced electrostatic interaction on the structure of the solvation shell is negligible. On the other hand, net induced dipolar contributions arise in the electrostatic interactions, in spite of the fact that the electric field is found to be inhomogeneous at the scale of the solute molecule. These induced dipolar interactions enhance the calculated absorption and emission shifts and contribute up to 40% to the electrostatic spectral shifts. They also modify or even dominate the fluctuations of the solute-solvent interaction and thus the solvation response. Both shifts and solvation dynamics are discussed in terms of dipole and quadrupole jumps. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of nonadiabatic electronic transitions and test it against exact integration of the time-dependent Schrödinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the nonadiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed. [ABSTRACT FROM AUTHOR]

The development in the 1950s and 1960s of crossed molecular beam methods for studying chemical reactions at the single-collision molecular level stimulated the need and desire for theoretical methods to describe these and other dynamical processes in molecular systems. Chemical dynamics theory has made great strides in the ensuing decades, so that methods are now available for treating the quantum dynamics of small molecular systems essentially completely. For the large molecular systems that are of so much interest nowadays (e.g., chemical reactions in solution, in clusters, in nanostructures, in biological systems, etc.), however, the only generally available theoretical approach is classical molecular dynamics (MD) simulations. Much effort is currently being devoted to the development of approaches for describing the quantum dynamics of these complex systems. This paper reviews some of these approaches, especially the use of semiclassical approximations for adding quantum effects to classical MD simulations, also showing some new versions that should make these semiclassical approaches even more practical and accurate. [ABSTRACT FROM AUTHOR]

A combined molecular dynamics (MD)+quantum mechanics (QM) method for studying processes on ionic surfaces is presented. Through the combination of classical MD and ab initio embedded-cluster calculations, this method allows the modeling of surface processes involving both the structural and dynamic features of the substrate, even for large-scale systems. The embedding approach used to link the information from the MD simulation to the cluster calculation is presented, and rigorous tests have been carried out to ensure the feasibility of the method. The electrostatic potential and electron density resulting from our embedded-cluster model have been compared with periodic slab results, and confirm the satisfying quality of our embedding scheme as well as the importance of applying embedding in our combined MD+QM approach. We show that a highly accurate representation of the Madelung potential becomes a prerequisite when the embedded-cluster approach is applied to temperature-distorted surface snapshots from the MD simulation. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su–Schrieffer–Heeger (Holstein–Peierls) model Hamiltonian. [ABSTRACT FROM AUTHOR]

Compares surface hopping and mean field approaches for simulating proton transfer reactions. Molecular dynamics with quantum transitions method used for calculations in the surface hopping method; Consistency between quantum probabilities and the fraction of trajectories in each adiabatic state.