Multireference quantum Monte Carlo study of the O4 molecule

S2O2q+ (q = 0, 1, and 2) Molecular Systems: Characterization and Atmospheric Planetary Implications

We use accurate ab initio methodologies at the coupled cluster level ((R)CCSD(T)) and its explicitly correlated version ((R)CCSD(T)-F12) to investigate the electronic structure, relative stability, and spectroscopy of the stable isomers of the [S₂O₂] system and of some of its cations and dications, with a special focus on the most relevant isomers that could be involved in terrestrial and planetary atmospheres. This work identifies several stable isomers (10 neutral, 8 cationic, and 5 dicationic), including trigonal-OSSO, cis-OSSO, and cyc-OSSO. For all these isomers, we calculated geometric parameters, fragmentation energies, and simple and double ionization energies of the neutral species. Several structures are identified for the first time, especially for the ionic species. Computations show that in addition to cis-OSSO and trans-OSSO proposed for the absorption in the near-UV spectrum of the Venusian atmosphere other S₂O₂, S₂O₂⁺, and S₂O₂²⁺ species may contribute. Moreover, the characterization of the stability of singly and doubly charged S₂O₂ entities can also be used for their identification by mass spectrometry and UV spectroscopy in the laboratory or in planetary atmospheres. In sum, the quest for the main UV absorber in Venus' atmosphere is not over, since the physical chemistry of sulfur oxides in Venus' atmosphere is far from being understood.

Global potential energy surface of ground state singlet spin O4

A new global potential energy for the singlet spin state O₄ system is reported using CASPT2/aug-cc-pVTZ ab initio calculations. The geometries for the six-dimensional surface are constructed using a novel point generation scheme that employs randomly generated configurations based on the beta distribution. The advantage of this scheme is apparent in the reduction of the number of required geometries for a reasonably accurate potential energy surface (PES) and the consequent decrease in the overall computational effort. The reported surface matches well with the recently published singlet surface by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. In addition to the O₄ PES, the ground state N₄ PES is also constructed using the point generation scheme and compared with the existing PES [Y. Paukku et al., J. Chem. Phys. 139, 044309 (2013)]. The singlet surface is constructed with the aim of studying high energy O₂-O₂ collisions and predicting collision induced dissociation cross section to be used in simulating non-equilibrium aerothermodynamic flows.

Predicting crystal structures and properties of matter under extreme conditions via quantum mechanics: the pressure is on

The role of quantum mechanical calculations in understanding and predicting the behavior of matter at extreme pressures is discussed in this feature contribution.

Using perturbatively selected configuration interaction in quantum Monte Carlo calculations

Defining accurate and compact trial wavefunctions leading to small statistical and fixed-node errors in quantum Monte Carlo (QMC) calculations is still a challenging problem. Here we propose to make use of selected configuration interaction (CI) expansions obtained by selecting the most important determinants through a perturbative criterion. A major advantage with respect to truncated CASSCF wavefunctions or CI expansions limited to a maximum number of excitations (e.g, CISD) is that much smaller expansions can be considered (many unessential determinants are avoided), an important practical point for efficient QMC calculations. The most important determinants entering first during the selection process (hierarchical construction) the main features of the nodal structure of the wavefunction can be expected to be obtained with a moderate number of determinants. Thanks to this property, the delicate problem of optimizing in a Monte Carlo framework the numerous linear and (or) nonlinear parameters of the determinantal part of the trial wavefunction could be avoided. As a first numerical example, the calculation of the ground-state energy of the oxygen atom is presented. The best DMC value reported so far is obtained.

Ab initio molecular dynamics studies of tetrasulfur. Dynamics coupling the C2v open and D2h closed forms of S4

Ab initio molecular dynamics (AIMD) simulations were performed on the closed D(2h) and open C(2v) isomers of tetrasulfur. After a careful calibration of the electronic structure method, the calculations were done using the BPW91/aug-cc-pVTZ method. This combination of method/basis set adequately reproduces the relative benchmark CCSD(T) energy difference [Matus, M.; Dixon, D.; Peterson, K. A.; Harkless, J. A. W.; Francisco, J. S. J. Chem. Phys. 2007, 127, 174305] between these two isomers and, crucially, the fact that the D(2h) structure is a transition state linking two equivalent (mirror images) C(2v) isomers. The trajectories show that the symmetric open C(2v) isomers interconvert when passing through the D(2h) closed transition state structure and that, unlike tetraoxygen, no three-dimensional structures arise. The dynamic vibrational analysis yields peaks in good agreement with the static CCSD(T) harmonic frequencies and explains higher peaks as overtones, thus showing that unlike previous AIMD DFT-based approaches, carefully calibrated exchange-correlation functionals can produce reliable molecular dynamics results for complex PESs as the one corresponding to the lowest singlet of S(4).