Thermal Gaussian molecular dynamics for quantum dynamics simulations of many-body systems: application to liquid para-hydrogen

A new method, here called thermal Gaussian molecular dynamics (TGMD), for simulating the dynamics of quantum many-body systems has recently been introduced [I. Georgescu and V. A. Mandelshtam, Phys. Rev. B 82, 094305 (2010)]. As in the centroid molecular dynamics (CMD), in TGMD the N-body quantum system is mapped to an N-body classical system. The associated both effective Hamiltonian and effective force are computed within the variational Gaussian wave-packet approximation. The TGMD is exact for the high-temperature limit, accurate for short times, and preserves the quantum canonical distribution. For a harmonic potential and any form of operator A?, it provides exact time correlation functions C(AB)(t) at least for the case of B, a linear combination of the position, x, and momentum, p, operators. While conceptually similar to CMD and other quantum molecular dynamics approaches, the great advantage of TGMD is its computational efficiency. We introduce the many-body implementation and demonstrate it on the benchmark problem of calculating the velocity time auto-correlation function for liquid para-hydrogen, using a system of up to N = 2592 particles.     

Relativistic, quantum electrodynamic and many-body effects in the water molecule

Ab initio calculations are reported of the electronic structure of the water molecule, based on the Dirac theory of the electron using correlation consistent basis sets. We calculate electron correlation corrections by second-order many-body perturbation theory, NMR shielding constants in the four-component relativistic interaction Hamiltonian formulation of QED, and report the first calculation of the electronic structure of a molecular system in which the Breit interaction is included both perturbatively and variationally in the calculation of the total energy.     

Density dynamics in some quantum systems

The quantum domain behavior of classical nonintegrable systems is well‐understood by the implementation of quantum fluid dynamics and quantum theory of motion. These approaches properly explain the quantum analogs of the classical Kolmogorov–Arnold–Moser type transitions from regular to chaotic domain in different anharmonic oscillators. Field‐induced tunneling and chaotic ionization in Rydberg atoms are also analyzed with the help of these theories. Quantum fluid density functional theory may be used to understand different time‐dependent processes like ion‐atom/molecule collisions, atom‐field interactions, and so forth. Regioselectivity as well as confined atomic/molecular systems and their reactivity dynamics have also been explained. © 2013 Wiley Periodicals, Inc.     

Few-body quantum and many-body classical hyperspherical approaches to reactions and to cluster dynamics

The hyperspherical method is a widely used and successful approach for the quantum treatment of elementary chemical processes. It has been mostly applied to three-atomic systems, and current progress is here outlined concerning the basic theoretical framework for the extension to four-body bound state and reactive scattering problems. Although most applications only exploit the advantages of the hyperspherical coordinate systems for the formulation of the few-body problem, the full power of the technique implies representations explicitly involving quantum hyperangular momentum operators as dynamical quantities and hyperspherical harmonics as basis functions. In terms of discrete analogues of these harmonics one has a universal representation for the kinetic energy and a diagonal representation for the potential (hyperquantization algorithm). Very recently, advances have been made on the use of the approach in classical dynamics, provided that a hyperspherical formulation is given based on “classical” definitions of the hyperangular momenta and related quantities. The aim of the present paper is to offer a retrospective and prospective view of the hyperspherical methods both in quantum and classical dynamics. Specifically, regarding the general quantum hyperspherical approaches for three- and four-body systems, we first focus on the basis set issue, and then we present developments on the classical formulation that has led to applications involving the implementations of hyperspherical techniques for classical molecular dynamics simulations of simple nanoaggregates.     

Model study of coherent quantum dynamics of hole states in functionalized semiconductor nanostructures

Functionalization of semiconductor nanocrystals can be achieved by anchoring organic ligands to the surface dangling bonds. The resulting surface complexes often introduce electronic states in the semiconductor band gap. These interband states sensitize the host material for photoabsorption at frequencies characteristic of the molecular adsorbates, leading to the well-known process of photoexcitation and subsequent femtosecond interfacial electron transfer. This paper investigates the relaxation dynamics of hole states, energetically localized deep in the semiconductor band gap, after the ultrafast electron-hole pair separation due to interfacial electron transfer. Mixed quantum-classical methods, based on mean-field nuclear dynamics approximated by ab initio density functional theory molecular dynamics simulations, reveal superexchange hole tunneling between adjacent adsorbate molecules in a model study of functionalized TiO2-anatase nanostructures. It is shown that electronic coherences can persist for hundreds of picoseconds under cryogenic and vacuum conditions, despite the partial intrinsic decoherence induced by thermal ionic motion, providing results of broad theoretical and experimental interest.     

Analysis of long-range effects in many-body correlation approaches for one-dimensional periodic systems

We analyze the divergent contributions to the Hamiltonian for extended, nonmetallic systems in one dimension, to both the ground-state correlation energy and to the correlated band structure. It is shown that the contribution from the long-range divergent part of the Hamiltonian tends to zero as 1/M⁻², where M is the extent of the troublesome lattice summations. Therefore, it is well justified to neglect such contributions in higher order many-body perturbation theory or coupled cluster treatments of the electronic structure for polymers. © 1997 John Wiley & Sons, Inc.     

Non-adiabatic molecular dynamics with quantum solvent effects

Three novel approaches extending quantum-classical non-adiabatic (NA) molecular dynamics (MD) to include quantum effects of solvent environments are described. In a standard NA-MD the solute subsystem is treated quantum mechanically, while the larger solvent part of a system is treated classically. The three novel approaches presented here are based on the Bohmian formulation of quantum mechanics, the stochastic Schrödinger equation for the evolution of open quantum systems and the quantized Hamilton dynamics generalization of classical mechanics. The approaches extend the standard NA-MD to incorporate the following quantum effects of the solvent. (1) Branching, i.e. the ability of solvent quantum wave packets to split and follow asymptotically diverging trajectories correlated with different quantum states of the solute. (2) Decoherence, i.e. loss of quantum interference within the solute subsystem induced by the diverging solvent trajectories. (3) Zero point energy that contributes to NA coupling and must be preserved during the energy exchange between solvent and solute degrees of freedom. The Bohmian quantum-classical mechanics, stochastic mean-field and quantized mean-field approximations incorporate the quantum solvent effects into the standard quantum-classical NA-MD in a straightforward and efficient way that can be easily applied to quantum dynamics of condensed phase chemical systems.     

Including quantum effects in the dynamics of complex (i.e., large) molecular systems

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.     

Relativistic many-body and QED calculations on atomic systems

A review is given of many-body perturbation methods, particularly in the all-order and coupled-cluster forms. Relativistic many-body schemes are analyzed in terms of one- and two-photon potentials, derived by means of QED. A complete second-order (nonradiative) calculation for He-like ions is presented, including repeated Breit interactions as well as the effects of retardation and of negative-energy states, but omitting the Lamb shift. Numerical results of some Lamb-shift calculations are also given. From the analysis, conclusions can be drawn concerning the accuracy of certain relativistic many-body approaches. © 1996 John Wiley & Sons, Inc.     

Communication: Conditions for one-photon coherent phase control in isolated and open quantum systems

Coherent control of observables using the phase properties of weak light that induces one-photon transitions is considered. Measurable properties are shown to be categorizable as either class A, where control is not possible, or class B, where control is possible. Using formal arguments, we show that phase control in open systems can be environmentally assisted.     

Angular momentum decoupling approximations in the quantum dynamics of reactive systems

Two methods for implementing angular momentum decoupling approximations in quantum mechanical reactive scattering examined. Applications of both reactive and nonreactive H + H₂ collisions indicate that for the most intense individual reactive transitions and for all degeneracy-averaged ones, these decoupling methods (especially the “proper” decoupling method) yield results in good agreement with those of fully-coupled calculations. However, for the less intense reactive transitions and all individual nonreactive transitions, very large errors can result from use of these approximate methods.     

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.     

Understanding quantum interference in coherent molecular conduction

Theory and experiment examining electron transfer through molecules bound to electrodes are increasingly focused on quantities that are conceptually far removed from current chemical understanding. This presents challenges both for the design of interesting molecules for these devices and for the interpretation of experimental data by traditional chemical mechanisms. Here, the concept of electronic coupling from theories of intramolecular electron transfer is extended and applied in the scattering theory (Landauer) formalism. This yields a simple sum over independent channels, that is then used to interpret and explain the unusual features of junction transport through cross-conjugated molecules and the differences among benzene rings substituted at the ortho, meta, or para positions.     

Simulation of silicon clusters from “quantum” Langevin molecular dynamics

We present a simulation method to determine from first principles the structure of low symmetry atomic systems. Our method is based on Langevin molecular dynamics and quantum mechanical interactions derived fromab initio pseudopotential calculations. The molecular dynamics time step with this approach can be one to two orders of magnitude larger than in the Car-Parrinello method, compensating for the time required for self-consistency at each step. Moreover, because the simulation is constrained to reside on the Born-Oppenheimer surface, this method can be used for insulating as well as metallic and charged systems. Application will be made to small silicon clusters.     

Polaronic quantum master equation theory of inelastic and coherent resonance energy transfer for soft systems

This work extends the theory of coherent resonance energy transfer [S. Jang, J. Chem. Phys. 131, 164101 (2009)] by including quantum mechanical inelastic effects due to modulation of donor-acceptor electronic coupling. Within the approach of the second order time local quantum master equation (QME) in the polaron picture and under the assumption that the bath degrees of freedom modulating the electronic coupling are independent of other modes, a general time evolution equation for the reduced system density operator is derived. Detailed expressions for the relaxation operators and inhomogeneous terms of the QME are then derived for three specific models of modulation in distance, axial angle, and dihedral angle, which are all approximated by harmonic oscillators. Numerical tests are conducted for a set of model parameters. Model calculation shows that the torsional modulation can make significant contribution to the relaxation and dephasing mechanisms.     

Influence of environment induced correlated fluctuations in electronic coupling on coherent excitation energy transfer dynamics in model photosynthetic systems

Two-dimensional photon-echo experiments indicate that excitation energy transfer between chromophores near the reaction center of the photosynthetic purple bacterium Rhodobacter sphaeroides occurs coherently with decoherence times of hundreds of femtoseconds, comparable to the energy transfer time scale in these systems. The original explanation of this observation suggested that correlated fluctuations in chromophore excitation energies, driven by large scale protein motions could result in long lived coherent energy transfer dynamics. However, no significant site energy correlation has been found in recent molecular dynamics simulations of several model light harvesting systems. Instead, there is evidence of correlated fluctuations in site energy-electronic coupling and electronic coupling-electronic coupling. The roles of these different types of correlations in excitation energy transfer dynamics are not yet thoroughly understood, though the effects of site energy correlations have been well studied. In this paper, we introduce several general models that can realistically describe the effects of various types of correlated fluctuations in chromophore properties and systematically study the behavior of these models using general methods for treating dissipative quantum dynamics in complex multi-chromophore systems. The effects of correlation between site energy and inter-site electronic couplings are explored in a two state model of excitation energy transfer between the accessory bacteriochlorophyll and bacteriopheophytin in a reaction center system and we find that these types of correlated fluctuations can enhance or suppress coherence and transfer rate simultaneously. In contrast, models for correlated fluctuations in chromophore excitation energies show enhanced coherent dynamics but necessarily show decrease in excitation energy transfer rate accompanying such coherence enhancement. Finally, for a three state model of the Fenna-Matthews-Olsen light harvesting complex, we explore the influence of including correlations in inter-chromophore couplings between different chromophore dimers that share a common chromophore. We find that the relative sign of the different correlations can have profound influence on decoherence time and energy transfer rate and can provide sensitive control of relaxation in these complex quantum dynamical open systems.     

Time-correlation function approach to many-body scattering in molecular systems

Energy transfer cross sections for collisions between atoms and molecular systems are related to time-correlation functions (TCFs) of the transition operators for scattering. The frequency Fourier transforms of the TCFs are studied within a Liouville-space formalism, and are shown to be lorentzian functions with frequency-dependent widths and shifts.