Electronic energy levels with the help of trajectory-guided random grid of coupled wave packets. I. Six-dimensional simulation of H2

Mixed quantum-classical approach to model non-adiabatic electron-nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron-nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional "virtual" nuclear wave packets whose interference with the "real" (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

Static Coherent States Method: One- and Two-Electron Laser-Induced Systems with Classical Nuclear Dynamics

In this report, we introduce the static coherent states (SCS) method for investigating quantum electron dynamics in a one- or two-electron laser-induced system. The SCS method solves the time-dependent Schrödinger equation (TDSE) both in imaginary and real times on the basis of a static grid of coherent states (CSs). Moreover, we consider classical dynamics for the nuclei by solving their Newtonian equations of motion. By implementing classical nuclear dynamics, we compute the electronic-state potential energy curves of H2+ in the absence and presence of an ultra-short intense laser field. We used this method to investigate charge migration in H2+. In particular, we found that the charge migration time increased exponentially with inter-nuclear distance. We also observed substantial charge localization for sufficiently long molecular bonds.

Complementary version of fermion coupled coherent states method and gram-schmidt algorithm: Theory and applications for electronic states of H2 and H2

In our previous report, we introduced a new version of Fermion coupled coherent states method (FCCS) which was especially suited for simulating the first symmetric spatial electronic state of two-electron systems. In this manuscript, we report a complementary version for FCCS method to simulate both of the first symmetric and antisymmetric spatial electronic states of two-electron systems. Moreover, the Gram-Schmidt orthogonalization process is employed to reach the excited states of the system. We apply this FCCS method and the original coupled coherent state method to simulate the energy of different electronic states of H2 and H2+, respectively. The results for the energy of computed electronic states of H2 and H2+ show a pretty good consistency with the exact values. © 2017 Wiley Periodicals, Inc.

Frozen Gaussian Wavepacket Study of the Ground State of the He Atom

The Rayleigh-Ritz functional is used in conjunction with an approximate time evolution to improve ab initio estimates of ground-state energies. The improvement is due in part to the introduction of a novel variational "normalization function" for the approximate propagator. An additional variational parameter was introduced in the form of a constant shift energy of the Hamiltonian. The approximate propagator used was the frozen Gaussian propagator; however, the trajectories evolved on the coherent-state averaged Hamiltonian (Q representation). For Coulombic forces, this removes the singularity, easing the computation. An additional variational parameter was the width parameter used for the coherent states appearing in the frozen Gaussian propagator. Using an initial combination of nine Gaussian functions for He, with an initial energy of -2.5115 au, the variational method, with a very short time interval of integration, led to an improved energy of -2.81 ± 0.04 au.

A local coherent-state approximation to system-bath quantum dynamics

A novel quantum method to deal with typical system-bath dynamical problems is introduced. Subsystem discrete variable representation and bath coherent-state sets are used to write down a multiconfigurational expansion of the wave function of the whole system. With the help of the Dirac-Frenkel variational principle, simple equations of motion--a kind of Schrodinger-Langevin equation for the subsystem coupled to (pseudo) classical equations for the bath--are derived. True dissipative dynamics at all times is obtained by coupling the bath to a secondary, classical Ohmic bath, which is modeled by adding a friction coefficient in the derived pseudoclassical bath equations. The resulting equations are then solved for a number of model problems, ranging from tunneling to vibrational relaxation dynamics. Comparison of the results with those of exact, multiconfiguration time-dependent Hartree calculations in systems with up to 80 bath oscillators shows that the proposed method can be very accurate and might be of help in studying realistic problems with very large baths. To this end, its linear scaling behavior with respect to the number of bath degrees of freedom is shown in practice with model calculations using tens of thousands of bath oscillators.

A version of diffusion Monte Carlo method based on random grids of coherent states. II. Six-dimensional simulation of electronic states of H2

We report a new version of the diffusion Monte Carlo (DMC) method, based on coherent-state quantum mechanics. Randomly selected grids of coherent states in phase space are used to obtain numerical imaginary time solutions of the Schrodinger equation, with an iterative refinement technique to improve the quality of the Monte Carlo grid. Accurate results were obtained, for the appropriately symmetrized two lowest states of the hydrogen molecule, by Monte Carlo sampling and six-dimensional propagation in the full phase space.