Search

The direct random-phase approximation (dRPA) is used to calculate and compare atomization energies for the HEAT set and ten selected molecules of the G2-1 set using both plane waves and Gaussian-type orbitals. We describe detailed procedures to obtain highly accurate and well converged results for the projector augmented-wave method as implemented in the Vienna Ab initio Simulation Package as well as the explicitly correlated dRPA-F12 method as implemented in the TURBOMOLE package. The two approaches agree within chemical accuracy (1 kcal/mol) for the atomization energies of all considered molecules, both for the exact exchange as well as for the RPA. The root mean-square deviation is 0.41 kcal/mol for the exact exchange (evaluated using density functional theory orbitals) and 0.33 kcal/mol for exact exchange plus correlation from the RPA. [ABSTRACT FROM AUTHOR]

With a He– H 3 + interaction potential obtained from advanced electronic structure calculations, we computed the vibration-rotation-tunneling (VRT) states of this complex for total angular momenta J from 0 to 9, both for the vibrational ground state and for the twofold degenerate v2 = 1 excited state of H 3 + . The potential has three equivalent global minima with depth De = 455.3 cm−1 for He in the plane of H 3 + , three equatorial saddle points that separate these minima with barriers of 159.5 cm−1, and two axial saddle points with energies of 243.1 cm−1 above the minima. The dissociation energies calculated for the complexes of He with ortho- H 3 + (o H 3 + ) and para- H 3 + (p H 3 + ) are D0 = 234.5 and 236.3 cm−1, respectively. Wave function plots of the VRT states show that they may be characterized as weakly hindered internal rotor states, delocalized over the three minima in the potential and with considerable amplitude at the barriers. Most of them are dominated by the jk = 10 and 11 rotational ground states of o H 3 + and p H 3 + , with the intermolecular stretching mode excited up to v = 4 inclusive. However, we also found excited internal rotor states: 33 in He–o H 3 + , and 22 and 21 in He–p H 3 + . The VRT levels and wave functions were used to calculate the frequencies and line strengths of all allowed v2 = 0 → 1 rovibrational transitions in the complex. Theoretical spectra generated with these results are compared with the experimental spectra in Paper II [Salomon et al., J. Chem. Phys. 156, 144308 (2022)] and are extremely helpful in assigning these spectra. This comparison shows that the theoretical energy levels and spectra agree very well with the measured ones, which confirms the high accuracy of our ab initio He– H 3 + interaction potential and of the ensuing calculations of the VRT states. [ABSTRACT FROM AUTHOR]

The rotationally resolved infrared (IR) spectrum of the He– H 3 + complex has been measured in a cryogenic ion trap experiment at a nominal temperature of 4 K. Predissociation of the stored complex has been invoked by excitation of the degenerate ν2 mode of the H 3 + sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He– H 3 + complex in its vibrational ground state, v = 0. Experimental and theoretical term energies agree within experimental accuracy for the rotational fine structure associated with the total angular momentum quantum number J and the parity e/f as well as for the coarse spacing of the lowest K states of the complex. This favorable comparison shows that the potential energy surface (PES) calculated in Paper I is accurate. The barriers between the three equivalent global minima in this PES are relatively low and the He– H 3 + complex is extremely floppy, with nearly unhindered internal rotation of the H 3 + sub-unit. The resulting Coriolis interactions couple the internal and end-over-end rotation of the complex and contribute significantly to the energy terms. They are observed both in experiment and theory and are, e.g., the origin of different rotational constants for states of e and f parity. Also in this respect, experiment and theory agree very well. Despite the assignment and analysis of many bands of the extremely rich IR spectrum of He– H 3 + , higher levels of excitation, including the complex stretching mode, need further attention. [ABSTRACT FROM AUTHOR]

An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Finally, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined. [ABSTRACT FROM AUTHOR]

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted. [ABSTRACT FROM AUTHOR]

The site-specific first microsolvation step of furan and some of its derivatives with methanol is explored to benchmark the ability of quantum-chemical methods to describe the structure, energetics, and vibrational spectrum at low temperature. Infrared and microwave spectra in supersonic jet expansions are used to quantify the docking preference and some relevant quantum states of the model complexes. Microwave spectroscopy strictly rules out in-plane docking of methanol as opposed to the top coordination of the aromatic ring. Contrasting comparison strategies, which emphasize either the experimental or the theoretical input, are explored. Within the harmonic approximation, only a few composite computational approaches are able to achieve a satisfactory performance. Deuteration experiments suggest that the harmonic treatment itself is largely justified for the zero-point energy, likely and by design due to the systematic cancellation of important anharmonic contributions between the docking variants. Therefore, discrepancies between experiment and theory for the isomer abundance are tentatively assigned to electronic structure deficiencies, but uncertainties remain on the nuclear dynamics side. Attempts to include anharmonic contributions indicate that for systems of this size, a uniform treatment of anharmonicity with systematically improved performance is not yet in sight. [ABSTRACT FROM AUTHOR]

A variety of approaches are presented for the computation of atomic and molecular correlation energies based on the Bethe–Salpeter equation in the framework of the adiabatic-connection fluctuation–dissipation theorem. The performance of the approaches is assessed by computing the total energies of the atoms H—Ne and the atomization energies of the high-accuracy extrapolated ab initio thermochemistry set of small molecules as well as by determining the bond lengths and harmonic vibrational frequencies of the metal monoxides MO with M=Ca—Zn. [ABSTRACT FROM AUTHOR]

The rotational spectrum of the parent isotopic species of HCCNCS, along with those of three 13 C singly substituted variants and one 34 S minor isotopologue, has been observed with high resolution using Fourier transform microwave spectroscopy from 4 to 26 GHz. Based on the observed spectral pattern and assignment of the 14N nuclear quadrupole hyperfine structure, the geometry of HCCNCS is effectively linear, which is further supported by the observation of ϱ-type doubling for rotational transitions in an excited bending state of the parent HCCNCS. The experimental observations are supported by high level quantum chemical calculations using the coupled-cluster singles and doubles model augmented by a perturbative correction for triple excitations, CCSD(T), level of theory that reveal an extremely shallow potential energy well along the CNC bending angle with a surprising basis set dependence. [ABSTRACT FROM AUTHOR]

Herein we present the results of a blind challenge to quantum chemical methods in the calculation of dimerization preferences in the low temperature gas phase. The target of study was the first step of the microsolvation of furan, 2-methylfuran and 2,5-dimethylfuran with methanol. The dimers were investigated through IR spectroscopy of a supersonic jet expansion. From the measured bands, it was possible to identify a persistent hydrogen bonding OH-O motif in the predominant species. From the presence of another band, which can be attributed to an OH-л interaction, we were able to assert that the energy gap between the two types of dimers should be less than or close to 1 kJ/mol across the series. These values served as a first evaluation ruler for the 12 entries featured in the challenge. A tentative stricter evaluation of the challenge results is also carried out, combining theoretical and experimental results in order to define a smaller error bar. The process was carried out in a double-blind fashion, with both theory and experimental groups unaware of the results on the other side, with the exception of the 2,5-dimethylfuran system which was featured in an earlier publication. [ABSTRACT FROM AUTHOR]

Guided by theoretical predictions, the hyperfine structures of the rotational spectra of mono- and bideuterated-water containing 17O have been experimentally investigated. To reach sub-Doppler resolution, required to resolve the hyperfine structure due to deuterium quadrupole coupling as well as to spin-rotation (SR) and dipolar spin-spin couplings, the Lamb-dip technique has been employed. The experimental investigation and in particular, the spectral analysis have been supported by high-level quantum-chemical computations employing coupled-cluster techniques and, for the first time, a complete experimental determination of the hyperfine parameters involved was possible. The experimentally determined 17O spin-rotation constants of D2 17O and HD17O were used to derive the paramagnetic part of the corresponding nuclear magnetic shielding constants. Together with the computed diamagnetic contributions as well as the vibrational and temperature corrections, the latter constants have been employed to confirm the oxygen nuclear magnetic shielding scale, recently established on the basis of spin-rotation data for H2 17O [ABSTRACT FROM AUTHOR]

The high accuracy extrapolated ab initio thermochemistry (HEAT) protocol is applied to compute the total atomization energy (TAE) and the heat of formation of benzene. Large-scale coupled-cluster calculations with more than 1500 basis functions and 42 correlated electrons as well as zero-point energies based on full cubic and (semi)diagonal quartic force fields obtained with the coupled-cluster singles and doubles with perturbative treatment of the triples method and atomic natural orbital (ANO) triple- and quadruple-zeta basis sets are presented. The performance of modifications to the HEAT scheme and the scaling properties of its contributions with respect to the system size are investigated. A purely quantum-chemical TAE and associated conservative error bar of 5463.0 ± 3.1 kJ mol-1 are obtained, while the corresponding 95% confidence interval, based on a statistical analysis of HEAT results for other and related molecules, is ± 1.8 kJ mol-1. The heat of formation of benzene is determined to be 101.5 ± 2.0 kJ mol-1 and 83.9 ± 2.1 kJ mol-1 at 0 K and 298.15 K, respectively. [ABSTRACT FROM AUTHOR]

The hyperfine structure in the rotational spectrum of water containing 17O has been investigated experimentally and by means of quantum-chemical calculations. The Lamb-dip technique has been used to resolve the hyperfine structure due to spin-rotation as well as spin-spin interactions and allowed the determination of the corresponding hyperfine parameters with high accuracy. The experimental investigation and, in particular, the analysis of the spectra have been supported by quantum-chemical computations at the coupled-cluster level. The experimental 17O isotropic spin-rotation constant of H217O has been used in a further step for the determination of the paramagnetic part of the corresponding nuclear magnetic shielding constant, whereas the diamagnetic contribution as well as vibrational and temperature corrections have been obtained from quantum-chemical calculations. This joint procedure leads to a value of 325.3(3) ppm for the oxygen shielding in H217O at 300 K, in good agreement with pure theoretical predictions, and in this way provides the basis for a new absolute oxygen shielding scale. [ABSTRACT FROM AUTHOR]

Results of large-scale coupled-cluster calculations of selected C2H2Si and CNHSi structural isomers are reported. Equilibrium molecular structures of a total of 12 molecules in their singlet electronic states have been calculated systematically employing the coupled-cluster singles and doubles model augmented by a perturbative correction for triple excitations (CCSD(T)) in combination with Dunning’s hierarchy of correlation consistent basis sets. In addition, anharmonic force fields were calculated to yield fundamental vibrational frequencies and rotation-vibration interaction constants αiA,B,C. The latter were used to determine empirical equilibrium structures reemp of two molecules — silacyclopropenylidene, c-C2H2Si, and silapropadienylidene, H2CCSi — for which sufficient isotopic data are available from literature. Very good agreement with theoretical equilibrium structures from CCSD(T) calculations employing core-valence basis sets of quadruple and quintuple-zeta quality — i.e., cc-pwCVQZ (337 basis functions), cc-pCV5Z, and cc-pwCV5Z (581 basis functions) is found — to within 0.001 Å for bond lengths and 0.1° for bond angles. Theoretical ground state rotational constants of HSiCN and HSiNC compare very favorably with experimental microwave data from literature, to within 0.15% (HSiCN) and 0.1% (HSiNC) for the B0 and C0 rotational constants. In the case of c-C2H2Si and H2CCSi this agreement is even better than 0.1%. For the latter two molecules effects of higher-level electron-correlation and relativity to the equilibrium geometry as well as the electronic contributions to the rotational constants are investigated. For eight molecules not yet studied at high spectral resolution in the gas-phase theoretical molecular parameters are provided to support future laboratory investigations. Theoretical vibrational fundamentals compare well with data of eight species studied previously with infrared matrix isolation spectroscopy. [ABSTRACT FROM AUTHOR]

Benchmark calculations of 19F nuclear magnetic shielding constants are presented for a set of 28 molecules. Near-quantitative accuracy (ca. 2 ppm deviation from experiment) is achieved if (1) electron correlation is adequately treated by employing the coupled-cluster singles and doubles (CCSD) model augmented by a perturbative correction for triple excitations [CCSD(T)], (2) large (uncontracted) basis sets are used, (3) gauge-including atomic orbitals are used to ensure gauge-origin independence, (4) calculations are performed at accurate equilibrium geometries [obtained from CCSD(T)/cc-pVTZ calculations correlating all electrons], and (5) vibrational averaging and temperature corrections via second-order vibrational perturbation theory (VPT2) are included. For the CCSD(T)/13s9p4d3f calculations corrected for vibrational effects, mean and standard deviation from experiment are -1.9 and 1.6 ppm, respectively. Less elaborate theoretical treatments result in larger errors. Consideration of relative shifts can reduce the mean deviation (through an appropriately chosen reference compound), but does not change the standard deviation. Density-functional theory calculations of absolute and relative 19F nuclear magnetic shielding constants are found to be, at best, as accurate as the corresponding Hartree–Fock self-consistent-field calculations and are not improved by consideration of vibrational effects. Molecular systems containing fluorine-oxygen, fluorine-nitrogen, and fluorine-fluorine bonds are found to be more challenging than the other investigated molecules for the considered theoretical methods. [ABSTRACT FROM AUTHOR]

Effects of increased basis-set size as well as a correlated treatment of the diagonal Born-Oppenheimer approximation are studied within the context of the high-accuracy extrapolated ab initio thermochemistry (HEAT) theoretical model chemistry. It is found that the addition of these ostensible improvements does little to increase the overall accuracy of HEAT for the determination of molecular atomization energies. Fortuitous cancellation of high-level effects is shown to give the overall HEAT strategy an accuracy that is, in fact, higher than most of its individual components. In addition, the issue of core-valence electron correlation separation is explored; it is found that approximate additive treatments of the two effects have limitations that are significant in the realm of <1 kJ mol-1 theoretical thermochemistry. [ABSTRACT FROM AUTHOR]

Electronic excitation energy transfer (EET) rates in rylene diimide dyads are calculated using second-order approximate coupled-cluster theory and time-dependent density functional theory. We investigate the dependence of the EET rates on the interchromophoric distance and the relative orientation and show that Förster theory works quantitatively only for donor-acceptor separations larger than roughly 5 nm. For smaller distances the EET rates are over- or underestimated by Förster theory depending on the respective orientation of the transition dipole moments of the chromophores. In addition to the direct transfer rates we consider bridge-mediated transfer originating from oligophenylene units placed between the chromophores. We find that the polarizability of the bridge significantly enhances the effective interaction. We compare our calculations to single molecule experiments on two types of dyads and find reasonable agreement between theory and experiment. [ABSTRACT FROM AUTHOR]

The rotationally resolved infrared (IR) spectrum of the He-H 3 + complex has been measured in a cryogenic ion trap experiment at a nominal temperature of 4 K. Predissociation of the stored complex has been invoked by excitation of the degenerate ν 2 sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He-H 3 + sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He-H 3 + complex in its vibrational ground state, v = 0. Experimental and theoretical term energies agree within experimental accuracy for the rotational fine structure associated with the total angular momentum quantum number J and the parity e/f as well as for the coarse spacing of the lowest K states of the complex. This favorable comparison shows that the potential energy surface (PES) calculated in Paper I is accurate. The barriers between the three equivalent global minima in this PES are relatively low and the He-H 3 + complex is extremely floppy, with nearly unhindered internal rotation of the H 3 + sub-unit. The resulting Coriolis interactions couple the internal and end-over-end rotation of the complex and contribute significantly to the energy terms. They are observed both in experiment and theory and are, e.g., the origin of different rotational constants for states of e and f parity. Also in this respect, experiment and theory agree very well. Despite the assignment and analysis of many bands of the extremely rich IR spectrum of He-H 3 + , higher levels of excitation, including the complex stretching mode, need further attention.

The rotational spectrum of the parent isotopic species of HCCNCS, along with those of three 13 C singly substituted variants and one 34 S minor isotopologue, has been observed with high resolution using Fourier transform microwave spectroscopy from 4 to 26 GHz. Based on the observed spectral pattern and assignment of the 14 N nuclear quadrupole hyperfine structure, the geometry of HCCNCS is effectively linear, which is further supported by the observation of ℓ -type doubling for rotational transitions in an excited bending state of the parent HCCNCS. The experimental observations are supported by high level quantum chemical calculations using the coupled-cluster singles and doubles model augmented by a perturbative correction for triple excitations, CCSD(T), level of theory that reveal an extremely shallow potential energy well along the CNC bending angle with a surprising basis set dependence.

Guided by theoretical predictions, the hyperfine structures of the rotational spectra of mono- and bideuterated-water containing (17)O have been experimentally investigated. To reach sub-Doppler resolution, required to resolve the hyperfine structure due to deuterium quadrupole coupling as well as to spin-rotation (SR) and dipolar spin-spin couplings, the Lamb-dip technique has been employed. The experimental investigation and in particular, the spectral analysis have been supported by high-level quantum-chemical computations employing coupled-cluster techniques and, for the first time, a complete experimental determination of the hyperfine parameters involved was possible. The experimentally determined (17)O spin-rotation constants of D2 (17)O and HD(17)O were used to derive the paramagnetic part of the corresponding nuclear magnetic shielding constants. Together with the computed diamagnetic contributions as well as the vibrational and temperature corrections, the latter constants have been employed to confirm the oxygen nuclear magnetic shielding scale, recently established on the basis of spin-rotation data for H2 (17)O [Puzzarini et al., J. Chem. Phys. 131, 234304 (2009)].

The hyperfine structure in the rotational spectrum of water containing (17)O has been investigated experimentally and by means of quantum-chemical calculations. The Lamb-dip technique has been used to resolve the hyperfine structure due to spin-rotation as well as spin-spin interactions and allowed the determination of the corresponding hyperfine parameters with high accuracy. The experimental investigation and, in particular, the analysis of the spectra have been supported by quantum-chemical computations at the coupled-cluster level. The experimental (17)O isotropic spin-rotation constant of H(2)(17)O has been used in a further step for the determination of the paramagnetic part of the corresponding nuclear magnetic shielding constant, whereas the diamagnetic contribution as well as vibrational and temperature corrections have been obtained from quantum-chemical calculations. This joint procedure leads to a value of 325.3(3) ppm for the oxygen shielding in H(2)(17)O at 300 K, in good agreement with pure theoretical predictions, and in this way provides the basis for a new absolute oxygen shielding scale.

Results of large-scale coupled-cluster calculations of selected C(2)H(2)Si and CNHSi structural isomers are reported. Equilibrium molecular structures of a total of 12 molecules in their singlet electronic states have been calculated systematically employing the coupled-cluster singles and doubles model augmented by a perturbative correction for triple excitations (CCSD(T)) in combination with Dunning's hierarchy of correlation consistent basis sets. In addition, anharmonic force fields were calculated to yield fundamental vibrational frequencies and rotation-vibration interaction constants alpha(i) (A,B,C). The latter were used to determine empirical equilibrium structures r(e) (emp) of two molecules - silacyclopropenylidene, c-C(2)H(2)Si, and silapropadienylidene, H(2)CCSi - for which sufficient isotopic data are available from literature. Very good agreement with theoretical equilibrium structures from CCSD(T) calculations employing core-valence basis sets of quadruple and quintuple-zeta quality - i.e., cc-pwCVQZ (337 basis functions), cc-pCV5Z, and cc-pwCV5Z (581 basis functions) is found - to within 0.001 A for bond lengths and 0.1 degrees for bond angles. Theoretical ground state rotational constants of HSiCN and HSiNC compare very favorably with experimental microwave data from literature, to within 0.15% (HSiCN) and 0.1% (HSiNC) for the B(0) and C(0) rotational constants. In the case of c-C(2)H(2)Si and H(2)CCSi this agreement is even better than 0.1%. For the latter two molecules effects of higher-level electron-correlation and relativity to the equilibrium geometry as well as the electronic contributions to the rotational constants are investigated. For eight molecules not yet studied at high spectral resolution in the gas-phase theoretical molecular parameters are provided to support future laboratory investigations. Theoretical vibrational fundamentals compare well with data of eight species studied previously with infrared matrix isolation spectroscopy.

Benchmark calculations of (19)F nuclear magnetic shielding constants are presented for a set of 28 molecules. Near-quantitative accuracy (ca. 2 ppm deviation from experiment) is achieved if (1) electron correlation is adequately treated by employing the coupled-cluster singles and doubles (CCSD) model augmented by a perturbative correction for triple excitations [CCSD(T)], (2) large (uncontracted) basis sets are used, (3) gauge-including atomic orbitals are used to ensure gauge-origin independence, (4) calculations are performed at accurate equilibrium geometries [obtained from CCSD(T)/cc-pVTZ calculations correlating all electrons], and (5) vibrational averaging and temperature corrections via second-order vibrational perturbation theory (VPT2) are included. For the CCSD(T)/13s9p4d3f calculations corrected for vibrational effects, mean and standard deviation from experiment are -1.9 and 1.6 ppm, respectively. Less elaborate theoretical treatments result in larger errors. Consideration of relative shifts can reduce the mean deviation (through an appropriately chosen reference compound), but does not change the standard deviation. Density-functional theory calculations of absolute and relative (19)F nuclear magnetic shielding constants are found to be, at best, as accurate as the corresponding Hartree-Fock self-consistent-field calculations and are not improved by consideration of vibrational effects. Molecular systems containing fluorine-oxygen, fluorine-nitrogen, and fluorine-fluorine bonds are found to be more challenging than the other investigated molecules for the considered theoretical methods.