An efficient, fully quantum mechanical real-time path integral method for including the effects of static disorder in the dynamics of systems coupled to common or local harmonic baths is presented. Rather than performing a large number of demanding calculations for different realizations of the system Hamiltonian, the influence of the bath is captured through a single evaluation of the path sum by grouping the system paths into equivalence classes of fixed system amplitudes. The method is illustrated with several analytical and numerical examples that show a variety of nontrivial effects arising from the interplay among coherence, dissipation, thermal fluctuations and geometric phases.
more »
« less
PathSum : A C++ and Fortran suite of fully quantum mechanical real-time path integral methods for (multi-)system + bath dynamics
This paper reports the release of PathSum, a new software suite of state-of-the-art path integral methods for studying the dynamics of single or extended systems coupled to harmonic environments. The package includes two modules, suitable for system–bath problems and extended systems comprising many coupled system–bath units, and is offered in C++ and Fortran implementations. The system–bath module offers the recently developed small matrix path integral (SMatPI) and the well-established iterative quasi-adiabatic propagator path integral (i-QuAPI) method for iteration of the reduced density matrix of the system. In the SMatPI module, the dynamics within the entanglement interval can be computed using QuAPI, the blip sum, time evolving matrix product operators, or the quantum–classical path integral method. These methods have distinct convergence characteristics and their combination allows a user to access a variety of regimes. The extended system module provides the user with two algorithms of the modular path integral method, applicable to quantum spin chains or excitonic molecular aggregates. An overview of the methods and code structure is provided, along with guidance on method selection and representative examples.
more »
« less
- Award ID(s):
- 1955302
- PAR ID:
- 10443266
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 22
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)We investigate the time evolution of the reduced density matrix (RDM) and its purity in the dynamics of a two-level system coupled to a dissipative harmonic bath, when the system is initially placed in one of its eigenstates. We point out that the symmetry of the initial condition confines the motion of the RDM elements to a one-dimensional subspace and show that the purity always goes through its maximally mixed value at some time during relaxation, but subsequently recovers and (under low-temperature, weakly dissipative conditions) can rise to values that approach unity. These behaviors are quantified through accurate path integral calculations. Under low-temperature, weakly dissipative conditions, we observe unusual, nonmonotonic population dynamics when the two-level system is initially placed in its ground state. We also analyze the origin of the system-bath interactions responsible for the nonmonotonic behavior of purity during relaxation. Our results show that classical dephasing processes arising from site level fluctuations lead to a monotonic decay of purity, and that the quantum mechanical decoherence events associated with spontaneous phonon emission are responsible for the subsequent recovery of purity. Last, we show that coupling with a low-temperature bath can purify a mixed two-level system. In the case of the maximally mixed initial RDM, the purity increases monotonically even during short time.more » « less
-
null (Ed.)The time evolution of the purity (the trace of the square of the reduced density matrix) and von Neumann entropy in a symmetric two-level system coupled to a dissipative harmonic bath is investigated through analytical arguments and accurate path integral calculations on simple models and the singly excited bacteriochlorophyll dimer. A simple theoretical analysis establishes bounds and limiting behaviors. The contributions to purity from a purely incoherent term obtained from the diagonal elements of the reduced density matrix, a term associated with the difference of the two eigenstate populations, and a third term related to the square of the time derivative of a site population, are discussed in various regimes. In the case of tunneling dynamics from a localized initial condition, the complex interplay among these contributions leads to the recovery of purity under low-temperature, weakly dissipative conditions. Memory effects from the bath are found to play a critical role to the dynamics of purity. It is shown that the strictly quantum mechanical decoherence process associated with spontaneous phonon emission is responsible for the long-time recovery of purity. These analytical and numerical results show clearly that the loss of quantum coherence during the evolution toward equilibrium does not necessarily imply the decay of purity, and that the time scales relevant to these two processes may be entirely different.more » « less
-
The influence functional (IF) encodes all the information required for calculating dynamical properties of a system in contact with its environment. A direct and simple procedure is introduced for extracting from a few numerical evaluations of the IF, without computing time correlation functions or evaluating integrals, the parameters required for path integral calculations, either within or beyond the harmonic mapping, and for assessing the accuracy of the harmonic bath approximation. Further, the small matrix decomposition of the path integral (SMatPI) is extended to anharmonic environments and the required matrices are constructed directly from the IF.more » « less
-
DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. The flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.more » « less
An official website of the United States government

