-
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 16 September, 2024; v1 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.
-
Feshbach resonances and weakly bound molecular states of boson-boson and boson-fermion NaK pairs
Authors:
Alexandra Viel,
Andrea Simoni
Abstract:
We study theoretically magnetically induced Feshbach resonances and near-threshold bound states in isotopic NaK pairs. Our calculations accurately reproduce Feshbach spectroscopy data on Na$^{40}$K and explain the origin of the observed multiplets in the p-wave [Phys. Rev. A 85, 051602(R) (2012)]. We apply the model to predict scattering and bound state threshold properties of the boson-boson Na…
▽ More
We study theoretically magnetically induced Feshbach resonances and near-threshold bound states in isotopic NaK pairs. Our calculations accurately reproduce Feshbach spectroscopy data on Na$^{40}$K and explain the origin of the observed multiplets in the p-wave [Phys. Rev. A 85, 051602(R) (2012)]. We apply the model to predict scattering and bound state threshold properties of the boson-boson Na$^{39}$K and Na$^{41}$K systems. We find that the Na$^{39}$K isotopic pair presents broad magnetic Feshbach resonances and favorable ground-state features for producing non-reactive polar molecules by two-photon association. Broad s-wave resonances are also predicted for Na$^{41}$K collisions.
△ Less
Submitted 1 April, 2016;
originally announced April 2016.
-
Transition from molecular complex to quantum solvation in OCS(He)_N
Authors:
F. Paesani,
A. Viel,
F. A. Gianturco,
K. B. Whaley
Abstract:
We present quantum calculations of the rotational energy levels and spectroscopic rotational constants of the linear OCS molecule in variable size clusters of 4He, using spectral evolution quantum Monte Carlo methods that allow excited states to be accessed without nodal constraints. The rotational constants of OCS are found to decrease monotonically from the gas phase value as the number of hel…
▽ More
We present quantum calculations of the rotational energy levels and spectroscopic rotational constants of the linear OCS molecule in variable size clusters of 4He, using spectral evolution quantum Monte Carlo methods that allow excited states to be accessed without nodal constraints. The rotational constants of OCS are found to decrease monotonically from the gas phase value as the number of helium atoms increases to N=6, after which the average constant increases to saturation at a value in excellent agreement with experimental measurements made on significantly larger clusters (N>1000). The minimum is shown to indicate a transition from a molecular complex to a quantum solvated molecule, with the former characterized by floppy but near rigid behavior, while the latter is characterized by non-zero permutation exchanges and a smaller extent of rigid coupling.
△ Less
Submitted 12 December, 2002;
originally announced December 2002.
-
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.