-
2500 vibronic eigenstates of the NO$_3$ radical
Authors:
Henrik R. Larsson,
Alexandra Viel
Abstract:
The nitrate radical NO$_3$ plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic $\tilde X ^2A_2'$ state. Using recently developed tensor network state…
▽ More
The nitrate radical NO$_3$ plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic $\tilde X ^2A_2'$ state. Using recently developed tensor network state methods, we are able to compute more than 2500 vibronic states, thereby increasing the number of computed full-dimensional states by a factor of 50, compared to previous work. While we obtain good agreement with experiment for most of the assigned vibronic levels, for several others, we observe striking disagreement. Further, for the antisymmetric bending motion we find remarkably large symmetry-induced level splittings that are larger than the zero-order reference. We discuss non-negligible nonadiabatic effects and show that the Born-Oppenheimer approximation leads to significant errors in the spectrum.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
Rovibrational (de-)excitation of H$_{2}$ by He revisited
Authors:
Hubert Jóźwiak,
Franck Thibault,
Alexandra Viel,
Piotr Wcisło,
François Lique
Abstract:
Collisional (de-)excitation of H$_{2}$ by helium plays an important role in the thermal balance and chemistry of various astrophysical environments, making accurate rate coefficients essential for the interpretation of observations of the interstellar medium. Our goal is to utilize a state-of-the-art potential energy surface (PES) to provide comprehensive state-to-state rate coefficients for He-in…
▽ More
Collisional (de-)excitation of H$_{2}$ by helium plays an important role in the thermal balance and chemistry of various astrophysical environments, making accurate rate coefficients essential for the interpretation of observations of the interstellar medium. Our goal is to utilize a state-of-the-art potential energy surface (PES) to provide comprehensive state-to-state rate coefficients for He-induced transitions among rovibrational levels of H$_{2}$. We perform quantum scattering calculations for the H$_{2}$-He system and provide state-to-state rate coefficients for 1 089 transitions between rovibrational levels of H$_{2}$ with internal energies up to 15 000 cm$^{-1}$ for temperatures ranging from 20 to 8 000 K. Our results show good agreement with previous calculations for pure rotational transitions between low-lying rotational levels, but we find significant discrepancies for rovibrational processes involving highly-excited rotational and vibrational states. We attribute these differences to two key factors: the broader range of intramolecular distances covered by ab initio points, and the superior accuracy of the PES, resulting from the utilization of the state-of-the-art quantum chemistry methods, compared to the previous lower-level calculations. Radiative transfer calculations performed with the new collisional data indicate that the population of rotational levels in excited vibrational states experiences significant modifications, highlighting the critical need for this updated dataset in models of high-temperature astrophysical environments.
△ Less
Submitted 16 November, 2023;
originally announced November 2023.
-
The principle underlying antiaromaticity
Authors:
Raphael J. F. Berger,
Alexandre Viel
Abstract:
Aromaticity is one of the most widely used chemical concepts. Current definitions are purely phenomenological and relate symmetry, reactive stability and the occurence of molecular diamagnetic response currents. The antithetical concept of antiaromaticity provides a connection between the contrary properties of structural instability or distortion out of higher symmetry, a small HOMO-LUMO gap, and…
▽ More
Aromaticity is one of the most widely used chemical concepts. Current definitions are purely phenomenological and relate symmetry, reactive stability and the occurence of molecular diamagnetic response currents. The antithetical concept of antiaromaticity provides a connection between the contrary properties of structural instability or distortion out of higher symmetry, a small HOMO-LUMO gap, and paramagnetic response currents. We reveal the principle that is underlying antiaromaticity in showing an intimate and strict symmetry induced relation between these properties. This principle can be proven and is formulated like: First order (and related) Jahn-Teller distorted molecules out of non-isometric point groups are prone to paramagnetic current susceptibility parallel to the main axis of symmetry. We show by the exemplary cases of cyclobutadiene, cylcooctatetraene, pentalene and manganese trifluoride how this principle works and discuss this new perspective on antiaromaticity.
△ Less
Submitted 22 February, 2019; v1 submitted 21 November, 2018;
originally announced November 2018.
-
Application of the spectral element method to the solution of the multichannel Schrödinger equation
Authors:
Andrea Simoni,
Alexandra Viel,
Jean-Michel Launay
Abstract:
We apply the spectral element method to the determination of scattering and bound states of the multichannel Schrödinger equation. In our approach the reaction coordinate is discretized on a grid of points whereas the internal coordinates are described by either purely diabatic or locally diabatic (diabatic-by-sector) bases. Bound levels and scattering matrix elements are determined with spectral…
▽ More
We apply the spectral element method to the determination of scattering and bound states of the multichannel Schrödinger equation. In our approach the reaction coordinate is discretized on a grid of points whereas the internal coordinates are described by either purely diabatic or locally diabatic (diabatic-by-sector) bases. Bound levels and scattering matrix elements are determined with spectral accuracy using relatively small numbers of points. The scattering problem is cast as a linear system solved using state-of-the-art sparse matrix non iterative packages. Boundary conditions can be imposed so to compute a single column of the matrix solution. A comparison with log-derivative propagators customarily used in molecular physics is performed. The same discretization scheme can also be applied to bound levels that are computed using direct scalable sparse-matrix solvers.
△ Less
Submitted 11 May, 2017;
originally announced May 2017.
-
Structure and spectroscopy of doped helium clusters using quantum Monte Carlo techniques
Authors:
Alexandra Viel,
K. Birgitta Whaley
Abstract:
We present a comparative study of the rotational characteristics of various molecule-doped 4He clusters using quantum Monte Carlo techniques. The theoretical conclusions obtained from both zero and finite temperature Monte Carlo studies confirm the presence of two different dynamical regimes that correlate with the magnitude of the rotational constant of the molecule, i.e., fast or slow rotors.…
▽ More
We present a comparative study of the rotational characteristics of various molecule-doped 4He clusters using quantum Monte Carlo techniques. The theoretical conclusions obtained from both zero and finite temperature Monte Carlo studies confirm the presence of two different dynamical regimes that correlate with the magnitude of the rotational constant of the molecule, i.e., fast or slow rotors. For a slow rotor, the effective rotational constant for the molecule inside the helium droplet can be determined by a microscopic two-fluid model in which helium densities computed by path integral Monte Carlo are used as input, as well as by direct computation of excited energy levels. For a faster rotor, the conditions for application of the two-fluid model for dynamical analysis are usually not fulfilled and the direct determination of excitation energies is then mandatory. Quantitative studies for three molecules are summarized, showing in each case excellent agreement with experimental results.
△ Less
Submitted 10 June, 2002;
originally announced June 2002.
-
Efficient implementation of the Projection Operator Imaginary Time Spectral Evolution (POITSE) method for excited states
Authors:
Patrick Huang,
Alexandra Viel,
K. Birgitta Whaley
Abstract:
We describe and systematically analyze new implementations of the Projection Operator Imaginary Time Spectral Evolution (POITSE) method for the Monte Carlo evaluation of excited state energies. The POITSE method involves the computation of a correlation function in imaginary time. Decay of this function contains information about excitation energies, which can be extracted by a spectral transfor…
▽ More
We describe and systematically analyze new implementations of the Projection Operator Imaginary Time Spectral Evolution (POITSE) method for the Monte Carlo evaluation of excited state energies. The POITSE method involves the computation of a correlation function in imaginary time. Decay of this function contains information about excitation energies, which can be extracted by a spectral transform. By incorporating branching processes in the Monte Carlo propagation, we compute these correlation functions with significantly reduced statistical noise. Our approach allows for the stable evaluation of small energy differences in situations where the previous POITSE implementation was limited by this noise.
△ Less
Submitted 5 March, 2002;
originally announced March 2002.
-
Importance Sampling in Rigid Body Diffusion Monte Carlo
Authors:
Alexandra Viel,
Mehul V. Patel,
Parhat Niyaz,
K. Birgitta Whaley
Abstract:
We present an algorithm for rigid body diffusion Monte Carlo with importance sampling, which is based on a rigorous short-time expansion of the Green's function for rotational motion in three dimensions. We show that this short-time approximation provides correct sampling of the angular degrees of freedom, and provides a general way to incorporate importance sampling for all degrees of freedom.…
▽ More
We present an algorithm for rigid body diffusion Monte Carlo with importance sampling, which is based on a rigorous short-time expansion of the Green's function for rotational motion in three dimensions. We show that this short-time approximation provides correct sampling of the angular degrees of freedom, and provides a general way to incorporate importance sampling for all degrees of freedom. The full importance sampling algorithm significantly improves both calculational efficiency and accuracy of ground state properties, and allows rotational and bending excitations in molecular van der Waals clusters to be studied directly.
△ Less
Submitted 4 September, 2001;
originally announced September 2001.
