A Phase-Space Electronic Hamiltonian For Vibrational Circular Dichroism.

We show empirically that a phase-space non-Born-Oppenheimer electronic Hamiltonian approach to quantum chemistry (where the electronic Hamiltonian is parametrized by both nuclear position and momentum, (,)) is both a practical and accurate means to recover vibrational circular dichroism spectra. We further hypothesize that such a phase-space approach may lead to very new dynamical physics beyond spectroscopic circular dichroism, with potential implications for understanding chiral induced spin selectivity (CISS), noting that classical phase-space approaches conserve the total nuclear plus electronic momentum, whereas classical Born-Oppenheimer approaches do not (they conserve only the nuclear momentum).

Practical phase-space electronic Hamiltonians for ab initio dynamics.

Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ') and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(X,P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie).

A simple one-electron expression for electron rotational factors.

Within the context of fewest-switch surface hopping (FSSH) dynamics, one often wishes to remove the angular component of the derivative coupling between states J and K. In a previous set of papers, Shu et al. [J.

Total angular momentum conservation in Ehrenfest dynamics with a truncated basis of adiabatic states.

We show that standard Ehrenfest dynamics does not conserve linear and angular momentum when using a basis of truncated adiabatic states. However, we also show that previously proposed effective Ehrenfest equations of motion [M. Amano and K.

Linear and angular momentum conservation in surface hopping methods.

We demonstrate that, for systems with spin-orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phase-space electronic Hamiltonians H(R, P) (i.

Surface hopping, electron translation factors, electron rotation factors, momentum conservation, and size consistency.

For a system without spin-orbit coupling, the (i) nuclear plus electronic linear momentum and (ii) nuclear plus orbital electronic angular momentum are good quantum numbers. Thus, when a molecular system undergoes a nonadiabatic transition, there should be no change in the total linear or angular momentum. Now, the standard surface hopping algorithm ignores the electronic momentum and indirectly equates the momentum of the nuclear degrees of freedom to the total momentum.

Modeling Spin-Dependent Nonadiabatic Dynamics with Electronic Degeneracy: A Phase-Space Surface-Hopping Method.

Nuclear Berry curvature effects emerge from electronic spin degeneracy and can lead to nontrivial spin-dependent (nonadiabatic) nuclear dynamics. However, such effects are not captured fully by any current mixed quantum-classical method such as fewest-switches surface hopping. In this work, we present a phase-space surface-hopping (PSSH) approach to simulate singlet-triplet intersystem crossing dynamics.

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.

Uniform semiclassical approximation for the Wigner 6j-symbol in terms of rotation matrices.

A new uniform asymptotic approximation for the Wigner 6j-symbol is given in terms of Wigner rotation matrices (d-matrices). The approximation is uniform in the sense that it applies for all values of the quantum numbers, even those near caustics. The derivation of the new approximation is not given, but the geometrical ideas supporting it are discussed and numerical tests are presented, including comparisons with the exact 6j-symbol and with the Ponzano-Regge approximation.

Phase space deformation and basis set optimization.

By deforming a given region of phase space---occupied by some unknown eigenfunctions one wishes to find---into a standard, integrable region, effective reductions in basis set size can be achieved. In one-dimensional problems we are able to achieve B/C=1+O(Planck' s constant), where B is the basis set size and C is the number of "converged" eigenfunctions. This result is confirmed by numerical examples, which also indicate exponential convergence as the basis set size is increased.

Phase integral theory, coupled wave equations, and mode conversion.

Phase integral or WKB theory is applied to multicomponent wave equations, i.e., wave equations in which the wave field is a vector, spinor, or tensor of some kind.