Context. Despite the fact that many sulfur-bearing molecules, ranging from simple diatomic species up to astronomical complex molecules, have been detected in the interstellar medium, the sulfur chemistry in space is largely unknown and a depletion in the abundance of S-containing species has been observed in the cold, dense interstellar medium. The chemical form of the missing sulfur has yet to be identified. Aims. For these reasons, in view of the fact that there is a large abundance of triatomic species harbouring sulfur, oxygen, and hydrogen, we decided to investigate the HSO radical in the laboratory to try its astronomical detection. Methods. High-resolution measurements of the rotational spectrum of the HSO radical were carried out within a frequency range well up into the THz region. Subsequently, a rigorous search for HSO in the two most studied high-mass star-forming regions, Orion KL and Sagittarius (Sgr) B2, and in the cold dark cloud Barnard 1 (B1-b) was performed. Results. The frequency coverage and the spectral resolution of our measurements allowed us to improve and extend the existing dataset of spectroscopic parameters, thus enabling accurate frequency predictions up to the THz range. These were used to derive the synthetic spectrum of HSO, by means of the MADEX code, according to the physical parameters of the astronomical source under consideration. For all sources investigated, the lack of HSO lines above the confusion limit of the data is evident. Conclusions. The derived upper limit to the abundance of HSO clearly indicates that this molecule does not achieve significant abundances in either the gas phase or in the ice mantles of dust grains.

We present an interface of the wavefunction-based quantum chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantum mechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/molecular mechanical (MM) operations are treated according to the previous MiMiC-based QM/MM implementation. Long-range electrostatic interactions are treated by a multipole expansion of the potential from the QM electron density to reduce the computational cost without loss of accuracy. Testing on model water/water systems, we verified that the CFOUR interface to MiMiC is robust, guaranteeing fast convergence of the self-consistent field cycles and optimal conservation of the energy during the integration of the equations of motion. Finally, we verified that the CFOUR interface to MiMiC is compatible with the use of a QM/QM multiple time-step algorithm, which effectively reduces the cost of ab initio MD (AIMD) or QM/MM-MD simulations using higher level wavefunction-based approaches compared to cheaper density functional theory-based ones. The new wavefunction-based AIMD and QM/MM-MD implementations were tested and validated for a large number of wavefunction approaches, including Hartree–Fock and post-Hartree–Fock methods like Møller–Plesset, coupled-cluster, and complete active space self-consistent field.

We combine internally contracted multireference coupled cluster theory with a four-component treatment of scalar-relativistic effects based on the spin-free Dirac-Coulomb Hamiltonian. This strategy allows for a rigorous treatment of static and dynamic correlation as well as scalar-relativistic effects, which makes it viable to describe molecules containing heavy transition elements. The use of a spin-free formalism limits the impact of the four-component treatment on the computational cost to the non-rate-determining steps of the calculations. We apply the newly developed method to the lowest singlet and triplet states of the monoxides of titanium, zirconium, and hafnium and show how the interplay between electronic correlation and relativistic effects explains the electronic structure of such molecules.

We present a scheme for the calculation of energies and analytic energy gradients within spin-free exact two-component (SFX2C) theory in its mean-field variant, which we refer to as SFX2C-mf. In the presented scheme, the Foldy-Wouthuysen transformation is carried out after the spin-free four-component Hartree-Fock treatment such that in electron-correlated calculations only the non-mean-field part of the two-electron interactions is handled in an untransformed manner. The formulation of analytic gradients requires some adjustments in comparison with the nonrelativistic case, i.e., the additional solution of the spin-free Dirac Coulomb coupled-perturbed Hartee-Fock equations together with a simplified treatment of orbital relaxation at the SFX2C-mf level. The improved accuracy of SFX2C-mf in comparison with SFX2C-1e is demonstrated in the calculation of energies, dipole moments, and electric-field gradients for the hydrogen halides HX, X = F-At. It is shown that the main contribution to the improvement stems from the elimination of the error at the Hartree-Fock (HF) level; however, the corresponding correlation contribution is also improved such that SFX2C-mf can be considered a suitable scheme for the treatment of heavy-element compounds for which the error of SFX2C-1e is rather substantial.

No abstract