Computing quantum dynamics in condensed matter systems is an open challenge due to the exponential scaling of exact algorithms with the number of degrees of freedom. Current methods try to reduce the cost of the calculation using classical dynamics as the key ingredient of approximations of the quantum time evolution. Two main approaches exist, quantum classical and semi-classical, but they suffer from various difficulties, in particular when trying to go beyond the classical approximation. It may then be useful to reconsider the problem focusing on statistical time-dependent averages rather than directly on the dynamics. In this paper, we discuss a recently developed scheme for calculating symmetrized correlation functions. In this scheme, the full (complex time) evolution is broken into segments alternating thermal and real-time propagation, and the latter is reduced to classical dynamics via a linearization approximation. Increasing the number of segments systematically improves the result with respect to full classical dynamics, but at a cost which is still prohibitive. If only one segment is considered, a cumulant expansion can be used to obtain a computationally efficient algorithm, which has proven accurate for condensed phase systems in moderately quantum regimes. This scheme is summarized in the second part of the paper. We conclude by outlining how the cumulant expansion formally provides a way to improve convergence also for more than one segment. Future work will focus on testing the numerical performance of this extension and, more importantly, on investigating the limit for the number of segments that goes to infinity of the approximate expression for the symmetrized correlation function to assess formally its convergence to the exact result.
Computing quantum dynamics in condensed matter systems is an open challenge due to the exponential scaling of exact algorithms with the number of degrees of freedom. Current methods try to reduce the cost of the calculation using classical dynamics as the key ingredient of approximations of the quantum time evolution. Two main approaches exist, quantum classical and semi-classical, but they suffer from various difficulties, in particular when trying to go beyond the classical approximation. It may then be useful to reconsider the problem focusing on statistical time-dependent averages rather than directly on the dynamics. In this paper, we discuss a recently developed scheme for calculating symmetrized correlation functions. In this scheme, the full (complex time) evolution is broken into segments alternating thermal and real-time propagation, and the latter is reduced to classical dynamics via a linearization approximation. Increasing the number of segments systematically improves the result with respect to full classical dynamics, but at a cost which is still prohibitive. If only one segment is considered, a cumulant expansion can be used to obtain a computationally efficient algorithm, which has proven accurate for condensed phase systems in moderately quantum regimes. This scheme is summarized in the second part of the paper. We conclude by outlining how the cumulant expansion formally provides a way to improve convergence also for more than one segment. Future work will focus on testing the numerical performance of this extension and, more importantly, on investigating the limit for the number of segments that goes to infinity of the approximate expression for the symmetrized correlation function to assess formally its convergence to the exact result.
The ability to efficiently and accurately calculate equilibrium time correlation functions of many-body condensed phase quantum systems is one of the outstanding problems in theoretical chemistry. The Nakajima-Zwanzig-Mori formalism coupled to the self-consistent solution of the memory kernel has recently proven to be highly successful for the computation of nonequilibrium dynamical averages. Here, we extend this formalism to treat symmetrized equilibrium time correlation functions for the spin-boson model. Following the first paper in this series [A. Montoya-Castillo and D. R. Reichman, J. Chem. Phys. 144, 184104 (2016)], we use a Dyson-type expansion of the projected propagator to obtain a self-consistent solution for the memory kernel that requires only the calculation of normally evolved auxiliary kernels. We employ the approximate mean-field Ehrenfest method to demonstrate the feasibility of this approach. Via comparison with numerically exact results for the correlation function Czz(t)=Re⟨σz(0)σz(t)⟩, we show that the current scheme affords remarkable boosts in accuracy and efficiency over bare Ehrenfest dynamics. We further explore the sensitivity of the resulting dynamics to the choice of kernel closures and the accuracy of the initial canonical density operator.
We present a way of combining real-time path integral Monte Carlo simulations with a maximum entropy numerical analytic continuation scheme in a new approach for calculating time correlation functions for finite temperature many body quantum systems. The real-time dynamics is expressed in the form of the symmetrized time correlation function, which is suitable for Monte Carlo methods, and several simulation techniques are presented for evaluating this function accurately up to moderate values of time. The symmetrized time correlation function is then analytically continued in combination with imaginary time data to obtain the real-time correlation function. We test this approach on several exactly solvable problems, including two one-dimensional systems, as well two cases of vibrational relaxation of a system coupled to a dissipative environment. The computed time correlation functions are in good agreement with exact results over several multiples of the thermal time βℏ, and exhibit a significant improvement over analytic continuation of imaginary time correlation functions. Moreover, we show how the method can be systematically improved.
An efficient method, the self-consistent hybrid method, is proposed for accurately simulating time-dependent quantum dynamics in complex systems. The method is based on an iterative convergence procedure for a dynamical hybrid approach. In this approach, the overall system is first partitioned into a “core” and a “reservoir” (an initial guess). The former is treated via an accurate quantum mechanical method, namely, the time-dependent multiconfiguration self-consistent field or multiconfiguration time-dependent Hartree approach, and the latter is treated via a more approximate method, e.g., classical mechanics, semiclassical initial value representations, quantum perturbation theories, etc. Next, the number of “core” degrees of freedom, as well as other variational parameters, is systematically increased to achieve numerical convergence for the overall quantum dynamics. The method is applied to two examples of quantum dissipative dynamics in the condensed phase: the spin-boson problem and the electronic resonance decay in the presence of a vibrational bath. It is demonstrated that the method provides a practical way of obtaining accurate quantum dynamical results for complex systems.
This thesis is focused on classical and semiclassical approximations of a specific quantum time correlation function, the “quantum fidelity.” Namely, we rigorously study the efficiency of a continuous class of algorithms for the evaluation of its classical limit and both the cost and the error of a particular semiclassical approximation of it. In the physics domain, quantum fidelity can be used to study the cross section of inelastic neutron scattering and phenomena like decoherence and irreversibility in NMR experiments. In chemical physics, quantum fidelity can be employed to compute the spectrum of a quantum state, to evaluate the accuracy of quantum dynamics on an approximate potential energy surface, and to verify the importance of nonadiabatic electronic transitions in molecules. Unfortunately, as for any numerical simulation of a quantum dynamical object, many-dimensional calculations of quantum fidelity are extremely expensive computationally. Classical methods, not including any quantum effect, are generally believed to be efficient approximations for the numerical evaluation of large systems’ properties when quantum effects are not important. We prove that, previous to this research, the computational costs of available algorithms for the simulation of classical fidelity increased exponentially with the number of degrees of freedom and we single out, within a continuous family of algorithms, the only one for which the number of trajectories needed for convergence is independent of the system’s dimensionality. Recently a semiclassical method, called “dephasing representation” (DR), has been proposed to approximate quantum fidelity accurately and efficiently. Although semiclassical methods, approximately including all types of quantum effects, are generally more accurate but much less efficient than classical methods, we prove that for the semiclassical DR there exists at least one algorithm for which the number of trajectories needed for convergence is independent of dimensionality and also that this semiclassical algorithm is even faster than the most efficient classical fidelity algorithm previously singled out. Encouraged by the surprising efficiency of the specific semiclassical method, we investigate the error of the DR compared to the quantum benchmark.
We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.
An effective and flexible numerical scheme is proposed to calculate the quantum and statistical mechanics of multidimensional systems in a simple economic way. The basic idea is to split the Hamiltonian operator into a reference separable part, whose solution can be obtained by a combination of analytic and numerical techniques, and a term containing nonseparable interactions and then to employ a symmetric decomposition of the time evolution operator, which is exact up to a high order in the time step. The method is applicable to a wide range of coupling potentials and requires numerical effort that scales only linearly with the number of degrees of freedom involved. To verify the utility of the present approach, two model systems with strongly anharmonic mode coupling are considered. The applications show that the method accurately describes the dynamics for fairly long times with moderate coupling strengths and is still much less arduous than a general numerically exact calculation.