No Arabic abstract
Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in nonadiabatic molecular dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, we introduce a regularized nuclear Bohm potential admitting solutions comprising a train of $delta$-functions which provide a finite-dimensional sampling of the hydrodynamic flow paths. The bohmion method inherits all the basic conservation laws from its underlying variational structure and captures electronic decoherence. After reviewing the general theory, the method is applied to the well-known Tully models, which are used here as benchmark problems. In the present case of study, we show that the new method accurately reproduces both electronic decoherence and nuclear population dynamics.
The Hamiltonian action of a Lie group on a symplectic manifold induces a momentum map generalizing Noethers conserved quantity occurring in the case of a symmetry group. Then, when a Hamiltonian function can be written in terms of this momentum map, the Hamiltonian is called `collective. Here, we derive collective Hamiltonians for a series of models in quantum molecular dynamics for which the Lie group is the composition of smooth invertible maps and unitary transformations. In this process, different fluid descriptions emerge from different factorization schemes for either the wavefunction or the density operator. After deriving this series of quantum fluid models, we regularize their Hamiltonians for finite $hbar$ by introducing local spatial smoothing. In the case of standard quantum hydrodynamics, the $hbar e0$ dynamics of the Lagrangian path can be derived as a finite-dimensional canonical Hamiltonian system for the evolution of singular solutions called `Bohmions, which follow Bohmian trajectories in configuration space. For molecular dynamics models, application of the smoothing process to a new factorization of the density operator leads to a finite-dimensional Hamiltonian system for the interaction of multiple (nuclear) Bohmions and a sequence of electronic quantum states.
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.
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.
Motivated by some models arising in quantum plasma dynamics, in this paper we study the Maxwell-Schrodinger system with a power-type nonlinearity. We show the local well-posedness in $H^2(mathbb{R}^3)times H^{3/2}(mathbb{R}^3)$ and the global existence of finite energy weak solutions, these results are then applied to the analysis of finite energy weak solutions for Quantum Magnetohydrodynamic systems.
We show that a novel, general phase space mapping Hamiltonian for nonadiabatic systems, which is reminiscent of the renowned Meyer-Miller mapping Hamiltonian, involves a commutator variable matrix rather than the conventional zero-point-energy parameter. In the exact mapping formulation on constraint space for phase space approaches for nonadiabatic dynamics, the general mapping Hamiltonian with commutator variables can be employed to generate approximate trajectory-based dynamics. Various benchmark model tests, which range from gas phase to condensed phase systems, suggest that the overall performance of the general mapping Hamiltonian is better than that of the conventional Meyer-Miller Hamiltonian.