Non-Born-Oppenheimer path in anti-Hermitian dynamics for nonadiabatic transitions

A quantum-classical Liouville formalism in a preconditioned basis and its connection with phase-space surface hopping

We revisit a recent proposal to model nonadiabatic problems with a complex-valued Hamiltonian through a phase-space surface hopping (PSSH) algorithm employing a pseudo-diabatic basis. Here, we show that such a pseudo-diabatic PSSH (PD-PSSH) ansatz is consistent with a quantum-classical Liouville equation (QCLE) that can be derived following a preconditioning process, and we demonstrate that a proper PD-PSSH algorithm is able to capture some geometric magnetic effects (whereas the standard fewest switches surface hopping approach cannot capture such effects). We also find that a preconditioned QCLE can outperform the standard QCLE in certain cases, highlighting the fact that there is no unique QCLE. Finally, we also point out that one can construct a mean-field Ehrenfest algorithm using a phase-space representation similar to what is done for PSSH. These findings would appear extremely helpful as far as understanding and simulating nonadiabatic dynamics with complex-valued Hamiltonians and/or spin degeneracy.

Quantum Chaos in the Dynamics of Molecules

Quantum chaos is reviewed from the viewpoint of "what is molecule?", particularly placing emphasis on their dynamics. Molecules are composed of heavy nuclei and light electrons, and thereby the very basic molecular theory due to Born and Oppenheimer gives a view that quantum electronic states provide potential functions working on nuclei, which in turn are often treated classically or semiclassically. Therefore, the classic study of chaos in molecular science began with those nuclear dynamics particularly about the vibrational energy randomization within a molecule. Statistical laws in probabilities and rates of chemical reactions even for small molecules of several atoms are among the chemical phenomena requiring the notion of chaos. Particularly the dynamics behind unimolecular decomposition are referred to as Intra-molecular Vibrational energy Redistribution (IVR). Semiclassical mechanics is also one of the main research fields of quantum chaos. We herein demonstrate chaos that appears only in semiclassical and full quantum dynamics. A fundamental phenomenon possibly giving birth to quantum chaos is "bifurcation and merging" of quantum wavepackets, rather than "stretching and folding" of the baker's transformation and the horseshoe map as a geometrical foundation of classical chaos. Such wavepacket bifurcation and merging are indeed experimentally measurable as we showed before in the series of studies on real-time probing of nonadiabatic chemical reactions. After tracking these aspects of molecular chaos, we will explore quantum chaos found in nonadiabatic electron wavepacket dynamics, which emerges in the realm far beyond the Born-Oppenheimer paradigm. In this class of chaos, we propose a notion of Intra-molecular Nonadiabatic Electronic Energy Redistribution (INEER), which is a consequence of the chaotic fluxes of electrons and energy within a molecule.

A phase-space semiclassical approach for modeling nonadiabatic nuclear dynamics with electronic spin

Chemical relaxation phenomena, including photochemistry and electron transfer processes, form a vigorous area of research in which nonadiabatic dynamics plays a fundamental role. However, for electronic systems with spin degrees of freedom, there are few if any applicable and practical quasiclassical methods. Here, we show that for nonadiabatic dynamics with two electronic states and a complex-valued Hamiltonian that does not obey time-reversal symmetry (as relevant to many coupled nuclear-electronic-spin systems), the optimal semiclassical approach is to generalize Tully’s surface hopping dynamics from coordinate space to phase space. In order to generate the relevant phase-space adiabatic surfaces, one isolates a proper set of diabats, applies a phase gauge transformation, and then diagonalizes the total Hamiltonian (which is now parameterized by both R and P). The resulting algorithm is simple and valid in both the adiabatic and nonadiabatic limits, incorporating all Berry curvature effects. Most importantly, the resulting algorithm allows for the study of semiclassical nonadiabatic dynamics in the presence of spin–orbit coupling and/or external magnetic fields. One expects many simulations to follow as far as modeling cutting-edge experiments with entangled nuclear, electronic, and spin degrees of freedom, e.g., experiments displaying chiral-induced spin selectivity.

Exploring dynamical electron theory beyond the Born-Oppenheimer framework: from chemical reactivity to non-adiabatically coupled electronic and nuclear wavepackets on-the-fly under laser field

Chemical theory and its application to dynamical electrons in molecules under intense electromagnetic fields is explored, in which we take an explicit account of nuclear nonadiabatic (kinematic) interactions along with simultaneous coupling with intense optical interactions. All the electronic wavefunctions studied here are necessarily time-dependent, and thereby beyond stationary state quantum chemistry based on the Born-Oppenheimer framework. As a general and tractable alternative framework with which to track the electronic and nuclear simultaneous dynamics, we propose an on-the-fly method to calculate the electron and nuclear wavepackets coupled along the branching non-Born-Oppenheimer paths, through which their bifurcations, strong quantum entanglement between nuclear electronic motions, and coherence and decoherence among the phases associated with them are properly represented. Some illustrative numerical examples are also reported, which are aimed at our final goals; real time tracking of nonadiabatic electronic states, chemical dynamics in densely degenerate electronic states coupled with nuclear motions and manipulation and/or creation of new electronic states in terms of intense lasers, and so on. Other examples are also presented as to how the electron wavepacket dynamics can be used to analyze chemical reactions, shedding a new light on some typical and conventional chemical reactions such as proton transfer followed by tautomerization.

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.