More Like this

1. Calculation of rovibrational eigenstates of H3+ using ScalIT

   https://doi.org/10.1063/5.0047823  

   Sarka, János; Das, Debojyoti; Poirier, Bill (April 2021, AIP Advances)

   H 3 + is a key player in molecular astrophysics, appearing in the interstellar medium and in the atmospheres of gas giants. It also plays an important role in star formation, and it has also been detected in supernova remnants. In theoretical chemistry, H3+ has long been a benchmark polyatomic system for high-level electronic-structure computations, as well as for quantum dynamics studies. In this work, exact quantum dynamical calculations are carried out for H3+, using the ScalIT suite of parallel codes, applied to two spectroscopically accurate potential energy surfaces. Specifically, rovibrational energy levels and wavefunctions are computed and labeled. Sixty vibrational states (for J = 0) are first determined, and then, rotational excitations for each of these “vibrational parent” states are computed up to total angular momentum J = 46, which is the highest value for which bound states of this molecule exist (D0 ∼ 35 000 cm−1). For these calculations, a very tight basis set convergence of a few 10−4 cm−1 (or less) has been achieved for almost all the computed energy levels. Where comparisons can be made, our results are found to agree well with earlier calculations and experimental data. 

   more » « less
2. Classical and nonclassical effects in surface hopping methodology for simulating coupled electronic-nuclear dynamics

   https://doi.org/10.1039/C9FD00042A  

   Martens, Craig C. (January 2019, Faraday Discussions)

   In this paper, we analyze the detailed quantum-classical behavior of two alternative approaches to simulating molecular dynamics with electronic transitions: the popular fewest switches surface hopping (FSSH) method, introduced by Tully in 1990 [Tully, \textit{J.~Chem.~Phys.}, 1990, \textbf{93}, 1061] and our recently developed quantum trajectory surface hopping (QTSH) method [Martens, \textit{J.~Phys.~Chem.~A}, 2019 \textbf{123}, 1110]. Both approaches employ an independent ensemble of trajectories that undergo stochastic transitions between electronic surfaces. The methods differ in their treatment of energy conservation, with FSSH imposing conservation of the classical kinetic plus potential energy by rescaling of the classical momenta when a surface hop occurs, while QTSH incorporates quantum forces throughout the dynamics which lead naturally to the conservation of the full quantum-classical energy. We investigate the population transfer and energy budget of the surface hopping methods for several simple model systems and compare with exact quantum result. In addition, the detailed dynamics of the trajectory ensembles in phase space are compared with the quantum evolution in the Wigner representation. Conclusions are drawn. 

   more » « less
3. Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

   https://doi.org/10.1039/d1cp04090a  

   Talbot, Justin J.; Head-Gordon, Martin; Miller, William H.; Cotton, Stephen J. (February 2022, Physical Chemistry Chemical Physics)

   Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology. 

   more » « less
4. Surface crossing and energy flow in many-dimensional quantum systems

   https://doi.org/10.1073/pnas.2221690120  

   Zhang, Chenghao; Gruebele, Martin; Logan, David E.; Wolynes, Peter G. (February 2023, Proceedings of the National Academy of Sciences)

   Energy flow in molecules, like the dynamics of other many-dimensional finite systems, involves quantum transport across a dense network of near-resonant states. For molecules in their electronic ground state, the network is ordinarily provided by anharmonic vibrational Fermi resonances. Surface crossing between different electronic states provides another route to chaotic motion and energy redistribution. We show that nonadiabatic coupling between electronic energy surfaces facilitates vibrational energy flow and, conversely, anharmonic vibrational couplings facilitate nonadiabatic electronic state mixing. A generalization of the Logan–Wolynes theory of quantum energy flow in many-dimensional Fermi resonance systems to the two-surface case gives a phase diagram describing the boundary between localized quantum dynamics and global energy flow. We explore these predictions and test them using a model inspired by the problem of electronic excitation energy transfer in the photosynthetic reaction center. Using an explicit numerical solution of the time-dependent Schrödinger equation for this ten-dimensional model, we find quite good agreement with the expectations from the approximate analytical theory. 

   more » « less
5. Probing the Dependence of Long-Range, Four-Atom Interactions on Intermolecular Orientation. 4. The Dissociation Dynamics of H2/D2···ICl(B,v′=3) and the Observation of Efficient Vibrational–Rotational Energy Transfer

   https://doi.org/10.1021/acs.jpca.2c05817  

   Darr, Joshua P.; Loomis, Richard A. (November 2022, The Journal of Physical Chemistry A)

   The vibrational predissociation dynamics of H2/D2···I35Cl(B,v′=3) complexes containing both para- and ortho-hydrogen prepared in different intermolecular vibrational levels were investigated. The Δv = −1 I35Cl(B,v = 2,j) rotational product-state distributions measured for excitation to the lowest-energy T-shaped levels of these complexes are mostly bimodal. The rotational distributions measured for excitation of the H2···I35Cl(B,v′=3) complexes are colder than those of the D2···I35Cl(B,v′=3) complexes, and there are only slight differences between those measured for the para- and ortho-hydrogen containing complexes. Excitation of the delocalized bending levels results in slightly colder rotational product-state distributions. The distributions suggest the dynamics result from more than impulsive dissociation off of the inner repulsive wall of the lower-energy H2/D2 + I35Cl(B,v = 2) potential surfaces of the products. The depths of these potentials and the energies available to these products also contribute to the dynamics. The formation of the Δv = −2, I35Cl(B,v = 1) product channel was only identified for excitation of levels within the ortho(j = 0)-D2 + I35Cl(B,v′=3) potential. The formation of this channel occurs via I35Cl(B,v′=3) vibrational to D2 rotational energy transfer forming the ortho(j = 2)-D2 + I35Cl(B,v = 1,j) products. 

   more » « less