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.     

Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III

Several thermodynamic properties of ice Ih, II, and III are studied by a quasi-harmonic approximation and compared to results of quantum path integral and classical simulations. This approximation allows to obtain thermodynamic information at a fraction of the computational cost of standard simulation methods, and at the same time permits studying quantum effects related to zero-point vibrations of the atoms. Specifically, we have studied the crystal volume, bulk modulus, kinetic energy, enthalpy, and heat capacity of the three ice phases as a function of temperature and pressure. The flexible q-TIP4P/F model of water was employed for this study, although the results concerning the capability of the quasi-harmonic approximation are expected to be valid independently of the employed water model. The quasi-harmonic approximation reproduces with reasonable accuracy the results of quantum and classical simulations showing an improved agreement at low temperatures (T< 100 K). This agreement does not deteriorate as a function of pressure as long as it is not too close to the limit of mechanical stability of the ice phases.     

Including quantum subsystem character within classical equilibrium simulations

A mixed quantum/classical density matrix approximation is derived. The density matrix makes use of quantum subsystem vibrational wave functions. The diagonal of the density matrix can be used as an equilibrium distribution in Monte Carlo simulations. The approximate distribution compares well with the path integral distribution for a model system. Since it includes quantum subsystem information, it performs much better than the quadratic Feynman-Hibbs distribution. These types of distributions can aid in including quantum vibrational information in otherwise classical simulations.     

Phase control in the vibrational qubit

In order to use molecular vibrations for quantum information processing one should be able to shape infrared laser pulses so that they can play the role of accurate quantum gates and drive the required vibrational transitions. In this paper we studied theoretically how the relative phase of the optimized transitions affects accuracy of the quantum gates in such a system. Optimal control theory and numerical propagation of laser-driven vibrational wave packets were employed. The dependencies observed for one-qubit gates NOT, pi-rotation, and Hadamard transform are qualitatively similar to each other. The results of the numerical tests agree well with the analytical predictions.     

The r-centroid approximation and molecular transitions

The implications of the r-centroid approximation on vibrational transition integrals are explored with reference to simple harmonic oscillator wavefunctions. It is quantitatively established for molecular bands of low vibrational quantum number that the existence of a localised well-developed uncancelled peak in the wavefunction product is necessary for the approximation to hold well.     

Quantum dot with N interacting electrons confined in a power-law external potential

The ground-state physical properties, such as electron density, chemical potential, and total energy, of a two-dimensional quantum dot with N interacting electrons confined in a power-law external potential are numerically determined by the Thomas-Fermi approximation. The effect of the confining potential on properties such as electron density and chemical potential is examined for both interacting and non-interacting systems. It is shown that the results of the calculations are in excellent agreement with those given in the literature. The results indicate that interactions and the shape of the confinement affect the density and thus the ground-state properties of the electrons significantly.     

Calculation of autocorrelation functions for electronic spectra

The autocorrelation function C(t) associated with the vibrational motion in the A²II_u state of N₂⁺ is calculated by three different methods. The classical method is based on classical trajectory calculations and provides good results except for secondary maxima appearing in the experimental autocorrelation function and resulting from interference effects. In the semiclassical method, the wavefunction is expanded in a basis of frozen gaussian functions. The secondary maxima are reproduced and the result is quite satisfactory. The superposition procedure consists in an expansion of the wavefunction in the eigenfunctions of the A state of N₂⁺ and leads to a very simple formula for the autocorrelation function C(t). The secondary maxima are accounted for and the calculated function agrees very well with experimental results.     

Quantum versus classical electron transfer energy as reaction coordinate for the aqueous Ru2+/Ru3+ redox reaction

Applying density functional theory (DFT)-based molecular dynamics simulation methods we investigate the effect of explicit treatment of electronic structure on the solvation free energy of aqueous Ru²⁺ and Ru³⁺.Our approach is based on the Marcus theory of redox half reactions, focussing on the vertical energy gap for reduction or oxidation of a single aqua ion. We compare the fluctuations of the quantum and classical energy gap along the same equilibrium ab initio molecular dynamics trajectory for each oxidation state. The classical gap is evaluated using a standard point charge model for the charge distribution of the solvent molecules (water). The quantum gap is computed from the full DFT electronic ground state energies of reduced and oxidized species, thereby accounting for the delocalization of the electron in the donor orbital and reorganization of the electron cloud after electron transfer (ET). The fluctuations of the quantum ET energy are well approximated by gaussian statistics giving rise to parabolic free energy profiles. The curvature is found to be independent of the oxidation state in agreement with the linear response assumption underlying Marcus theory. By contrast, the diabatic free energy curves evaluated using the classical gap as order parameter, while also quadratic, are asymmetric reflecting the difference in oxidation state. The response of these two order parameters is further analysed by a comparison of the spectral density of the fluctuations and the corresponding reorganization free energies.     

A semiclassical initial-value representation for quantum propagator and boltzmann operator

Starting from the position-momentum integral representation, we apply the correction operator method to the derivation of a uniform semiclassical approximation for the quantum propagator and then extend it to approximate the Boltzmann operator. In this approach, the involved classical dynamics is determined by the method itself instead of given beforehand. For the approximate Boltzmann operator, the corresponding classical dynamics is governed by a complex Hamiltonian, which can be described as a pair of real Hamiltonian systems. It is demonstrated that the semiclassical Boltzmann operator is exact for linear systems. A quantum propagator in the complex time is thus proposed and preliminary numerical results show that it is a reasonable approximation for calculating thermal correlation functions of general systems. © 2018 Wiley Periodicals, Inc.     

Complex trajectory method in time-dependent WKB

We present a significant improvement to a complex time-dependent WKB (CWKB) formulation developed by Boiron and Lombardi [J. Chem. Phys. 108, 3431 (1998)] in which the time-dependent WKB equations are solved along classical trajectories that propagate in complex space. Boiron and Lombardi showed that the method gives very good agreement with the exact quantum mechanical result as long as the wavefunction does not exhibit interference effects such as oscillations and nodes. In this paper, we show that this limitation can be overcome by superposing the contributions of crossing trajectories. Secondly, we demonstrate that the approximation improves when incorporating higher order terms in the expansion. Thirdly, equations of motion for caustics and Stokes lines are implemented to help overcome Stokes discontinuities. These improvements could make the CWKB formulation a competitive alternative to current time-dependent semiclassical methods.     

Simulating quantum dynamics in classical environments

Methods for simulating the dynamics of composite systems, where part of the system is treated quantum mechanically and its environment is treated classically, are discussed. Such quantum–classical systems arise in many physical contexts where certain degrees of freedom have an essential quantum character while the other degrees of freedom to which they are coupled may be treated classically to a good approximation. The dynamics of these composite systems are governed by a quantum–classical Liouville equation for either the density matrix or the dynamical variables which are operators in the Hilbert space of the quantum subsystem and functions of the classical phase space variables of the classical environment. Solutions of the evolution equations may be formulated in terms of surface-hopping dynamics involving ensembles of trajectory segments interspersed with quantum transitions. The surface-hopping schemes incorporate quantum coherence and account for energy exchanges between the quantum and classical degrees of freedom. Various simulation algorithms are discussed and illustrated with calculations on simple spin-boson models but the methods described here are applicable to realistic many-body environments.     

Surface-hopping dynamics and decoherence with quantum equilibrium structure

In open quantum systems, decoherence occurs through interaction of a quantum subsystem with its environment. The computation of expectation values requires a knowledge of the quantum dynamics of operators and sampling from initial states of the density matrix describing the subsystem and bath. We consider situations where the quantum evolution can be approximated by quantum-classical Liouville dynamics and examine the circumstances under which the evolution can be reduced to surface-hopping dynamics, where the evolution consists of trajectory segments exclusively evolving on single adiabatic surfaces, with probabilistic hops between these surfaces. The justification for the reduction depends on the validity of a Markovian approximation on a bath averaged memory kernel that accounts for quantum coherence in the system. We show that such a reduction is often possible when initial sampling is from either the quantum or classical bath initial distributions. If the average is taken only over the quantum dispersion that broadens the classical distribution, then such a reduction is not always possible.     

On the classical behavior, coherent and incoherent contributions to scattering of atoms from a surface

In this work we calculate the coherent and incoherent contributions to the diffraction probability of atoms scattered by a surface. We are interested in exploring the relative importance of each of these contributions, and compare them with results obtained from molecular dynamics calculations. To achieve this goal, we employed a method developed by Heller that consists of writing the incident plane wave as a sum of Gaussian wave packets, propagating them independently by using the time dependent Schrödinger equation, and constructing the scattered wave function by adding coherently the scattered packets. For the system studied, the molecular dynamics results show the largest intensity in the specular region and also display a classical rainbow structure. On the other hand, the quantum results exhibit diffraction features, with the coherent contribution accounting for most of the total intensity probability.     

Grid‐based density functional calculations of many‐electron systems

Exploratory variational pseudopotential density functional calculations are performed for the electronic properties of many‐electron systems in the 3D cartesian coordinate grid (CCG). The atom‐centered localized gaussian basis set, electronic density, and the two‐body potentials are set up in the 3D cubic box. The classical Hartree potential is calculated accurately and efficiently through a Fourier convolution technique. As a first step, simple local density functionals of homogeneous electron gas are used for the exchange‐correlation potential, while Hay‐Wadt‐type effective core potentials are employed to eliminate the core electrons. No auxiliary basis set is invoked. Preliminary illustrative calculations on total energies, individual energy components, eigenvalues, potential energy curves, ionization energies, and atomization energies of a set of 12 molecules show excellent agreement with the corresponding reference values of atom‐centered grid as well as the grid‐free calculation. Results for three atoms are also given. Combination of CCG and the convolution procedure used for classical Coulomb potential can provide reasonably accurate and reliable results for many‐electron systems. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Path integral based calculations of symmetrized time correlation functions. I

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.