Complete Basis Set Limit Extrapolation Research Articles

ExoMol line lists – LIII: empirical rovibronic spectra of yttrium oxide

ABSTRACT Empirical line lists BRYTS for the open shell molecule 89Y16O (yttrium oxide) and its isotopologues are presented. The line lists cover the six lowest electronic states: $X\, {}^{2}\Sigma ^{+}$, $A\, {}^{2}\Pi$, $A^{\prime }\, {}^{2}\Delta$, $B\, {}^{2}\Sigma ^{+}$, $C\, {}^{2}\Pi$, and $D\, {}^{2}\Sigma ^{+}$ up to 60 000 cm−1 (<0.167 $\mu$m) for rotational excitation up to J = 400.5. An ab initio spectroscopic model consisting of potential energy curves (PECs), spin–orbit, and electronic angular momentum couplings is refined by fitting to experimentally determined energies of YO, derived from published YO experimental transition frequency data. The model is complemented by empirical spin-rotation and Λ-doubling curves and ab initio dipole moment and transition dipole moment curves computed using multireference configuration interaction (MRCI). The ab initio PECs computed using the complete basis set limit extrapolation and the coupled-cluster CCSD(T) method with its higher quality provide an excellent initial approximation for the refinement. Non-adiabatic coupling curves for two pairs of states of the same symmetry A/C and B/D are computed using a state-averaged complete active space self-consistant field theory (CASSCF) and used to build diabatic representations for the $A\, {}^{2}\Pi$, $C\, {}^{2}\Pi$, $B\, {}^{2}\Sigma ^{+}$, and $D\, {}^{2}\Sigma ^{+}$ curves. The experimentally derived energies of 89Y16O utilized in the fit are used to replace the corresponding calculated energy values in the BRYTS line list. Simulated spectra of YO show excellent agreement with the experiment, where it is available. Calculated lifetimes of YO are tuned to agree well with the experiment, where available. The BRYTS YO line lists are available from the ExoMol database (www.exomol.com).

A theoretical study of the photoelectron spectra of dichloroketene with accurate computation of ionization energies via complete basis set limit extrapolation

AbstractWe calculated the equilibrium geometries and harmonic vibrational frequencies of the ground state and five cationic states of dichloroketene using (TD‐)B3LYP, PBE0, and M06/M06‐2X approaches. The photoelectron spectra of dichloroketene were simulated by computing Franck‐Condon factors. The ionization energies were computed using the CCSD(T) approach with extrapolation to the complete basis set (CBS) limit. We propose two new CBS energy formulas (E = ECBS + Aexp(‐x) + B/(x−1) n, n = 2 or 3) and compare the performance of different CBS approaches. A new ionic state of dichloroketene belonging to the Cs point group is reported. This state is identified as the first excited state of Cl2CCO+ having a double‐well potential‐energy curve along the CCO bending mode with a barrier height of 1.335 eV. The simulated photoelectron spectra are in agreement with the experiment. The vertical ionization energies calculated via spectral simulation are more accurate compared with those obtained at the ground‐state structure. Among the CBS formulas used, the proposed ansatz with n = 2 performs best, with a mean absolute error of 0.021 and 0.012 eV for the adiabatic and vertical ionization energies, respectively.

Anion Photoelectron Spectroscopy and High Level Ab Initio Calculations of the Halide–Nitric Oxide Dimer Complexes

Anion photoelectron spectra are presented for gas phase complexes formed between halide anions and nitric oxide, X−⋯NO where X− = Cl−, Br−, and I−. Electron binding energies are experimentally determined to be 3.82, 3.51, and 3.17 eV. Results from CCSD(T)/aug-cc-pVTZ calculations are presented for the anion species, whereby a single minimum of Cs symmetry is predicted. Binding energies (D0) of 15.3, 13.3, and 11.7 kJ mol−1 are predicted from complete basis set limit extrapolation, and are found to be in line with previous experimental studies.

Towards a quantum chemical protocol for the prediction of rovibrational spectroscopic data for transition metal molecules: Exploration of CuCN, CuOH, and CuCCH.

High accuracy electronic structure computations for small transition metal-containing molecules have been a long term challenge. Due to coupling between electronic and nuclear wave functions, even experimental/theoretical identification of the ground electronic state requires tremendous efforts. Quartic force fields (QFFs) are effective ab initio tools for obtaining reliable anharmonic spectroscopic properties. However, the method that employs complete basis set limit extrapolation ("C"), consideration of core electron correlation ("cC"), and inclusion of scalar relativity ("R") to produce the energy points on the QFF, the composite CcCR methodology, has not yet been utilized to study inorganic spectroscopy. This work takes the CcCR methodology and adapts it to test whether such an approach is conducive for the closed-shell, copper-containing molecules CuCN, CuOH, and CuCCH. Gas phase rovibrational data are provided for all three species in their ground electronic states. Equilibrium geometries and many higher-order rovibrational properties show good agreement with earlier studies. However, there are notable differences, especially in computation of fundamental vibrational frequencies. Even with further additive corrections for the inner core electron correlation and coupled cluster with full single, double, and triple substitutions (CCSDT), the differences are still larger than expected indicating that more work should follow for predicting rovibrational properties of transition metal molecules.

Ab initio potential energy surface for methane and carbon dioxide and application to vapor-liquid coexistence.

A six-dimensional intermolecular potential energy surface for a rigid methane (CH4) and carbon dioxide (CO2) dimer was developed from the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory. A total of 466 grid points distributed to 46 orientations were calculated from the complete basis set limit extrapolation based on up to aug-cc-pVQZ basis set. A modified site-site pair potential function was proposed for rapid representation of the high level ab initio calculations. A nonadditive three-body interaction was represented by the Axilrod-Teller-Muto expression for mixtures with the polarizability and the London dispersion constant of each molecule. Second to fourth virial coefficients of CH4 and CO2 mixtures were calculated using both the Mayer sampling Monte Carlo method and the present potential functions. The virial equation of state derived from these coefficients was used to predict the pVT values and showed good agreement with experimental data below 200 bar at 300 K. The vapor-liquid coexistence curves of pure CH4, CO2 and their mixtures were presented with the aid of Gibbs ensemble Monte Carlo simulations. The predicted tie lines agreed with the experimental data within the uncertainties up to near the critical point.

Complete basis set limit extrapolation calculations with PNOF3

Complete basis set (CBS) extrapolated total energies have been calculated for the ground state and lowest-lying excited state with different spin of the second-row atoms (Li–Ne) using the recently developed spin-conserving PNOF3 natural orbital functional. CBS results are obtained by extrapolating power and exponential functions to the results obtained using the correlation-consistent aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets. The extrapolated total energy values are very accurate as compared to the benchmark data, and reproduce near 100% of the total energy. Interestingly, the benchmark data are lower bounds for the PNOF3 results. Finally, we have found that the calculated PNOF3 excitation energies are very accurate as comparing with both their corresponding CCSD(T, full) values and their experimental marks.

H2O, H2, HF, F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin–spin coupling constants (SSCCs) in the BHandH/pcJ‐n and BHandH/XZP Kohn–Sham limits

Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H(2)O, H(2), HF, F(2) and F(2)O nuclear isotropic shielding constants in the BHandH Kohn-Sham basis set limit was shown. The results of two- and three-parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization-consistent basis sets pcS-n and pcJ-n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ-n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems.