Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation
218
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time-evolution of the nuclear and electronic phase space variables as well as the Monodromy matrix, under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve-crossing in model two-level systems and show that in the quantum limit MQC-IVR is in good agreement with the exact quantum results, whereas in the classical limit the method yields results in keeping with mean-field approaches like the Linearized Semiclassical IVR. Finally, exploiting the ability of MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy.
rate research
Read More
A method for carrying out semiclassical initial value representation calculations using first-principles molecular dynamics (FP-SC-IVR) is presented. This method can extract the full vibrational power spectrum of carbon dioxide from a single trajectory providing numerical results that agree with experiment even for Fermi resonant states. The computational demands of the method are comparable to those of classical single-trajectory calculations, while describing uniquely quantum features such as the zero-point energy and Fermi resonances. By propagating the nuclear degrees of freedom using first-principles Born-Oppenheimer molecular dynamics, the stability of the method presented is improved considerably when compared to dynamics carried out using fitted potential energy surfaces and numerical derivatives.
The mapping approach addresses the mismatch between the continuous nuclear phase space and discrete electronic states by creating an extended, fully continuous phase space using a set of harmonic oscillators to encode the populations and coherences of the electronic states. Existing quasiclassical dynamics methods based on mapping, such as the linearised semiclassical initial value representation (LSC-IVR) and Poisson bracket mapping equation (PBME) approaches, have been shown to fail in predicting the correct relaxation of electronic-state populations following an initial excitation. Here we generalise our recently published modification to the standard quasiclassical approximation for simulating quantum correlation functions. We show that the electronic-state population operator in any system can be exactly rewritten as a sum of a traceless operator and the identity operator. We show that by treating the latter at a quantum level instead of using the mapping approach, the accuracy of traditional quasiclassical dynamics methods can be drastically improved, without changes to their underlying equations of motion. We demonstrate this approach for the seven-state Frenkel-Exciton model of the Fenna-Matthews-Olson light harvesting complex, showing that our modification significantly improves the accuracy of traditional mapping approaches when compared to numerically exact quantum results.
We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in protonated glycine molecule [C. Aieta et. al. Nat. Commun. 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, and to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond-lengths, angles and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching the typical harmonic nodal pattern is absent in the anharmonic distribution.
Simulating the nonadiabatic dynamics of condensed-phase systems continues to pose a significant challenge for quantum dynamics methods. Approaches based on sampling classical trajectories within the mapping formalism, such as the linearized semiclassical initial value representation (LSC-IVR), can be used to approximate quantum correlation functions in dissipative environments. Such semiclassical methods however commonly fail in quantitatively predicting the electronic-state populations in the long-time limit. Here we present a suggestion to minimize this difficulty by splitting the problem into two parts, one of which involves the identity, and treating this operator by quantum-mechanical principles rather than with classical approximations. This strategy is applied to numerical simulations of spin-boson model systems, showing its potential to drastically improve the performance of LSC-IVR and related methods with no change to the equations of motion or the algorithm in general, but rather by simply using different functional forms of the observables.
We derive an exact quantum propagator for nonadiabatic dynamics in multi-state systems using the mapping variable representation, where classical-like Cartesian variables are used to represent both continuous nuclear degrees of freedom and discrete electronic states. The resulting expression is a Moyal series that, when suitably approximated, can allow for the use of classical dynamics to efficiently model large systems. We demonstrate that different truncations of the exact propagator lead to existing approximate semiclassical and mixed quantum-classical methods and we derive an associated error term for each method. Furthermore, by combining the imaginary-time path-integral representation of the Boltzmann operator with the exact propagator, we obtain an analytic expression for thermal quantum real-time correlation functions. These results provide a rigorous theoretical foundation for the development of accurate and efficient classical-like dynamics to compute observables such as electron transfer reaction rates in complex quantized systems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
