skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A study of core-excited states of organic molecules computed with the generalized active space driven similarity renormalization group
This work examines the accuracy and precision of x-ray absorption spectra computed with a multireference approach that combines generalized active space (GAS) references with the driven similarity renormalization group (DSRG). We employ the x-ray absorption benchmark of organic molecule (XABOOM) set, consisting of 116 transitions from mostly organic molecules [Fransson et al., J. Chem. Theory Comput. 17, 1618 (2021)]. Several approximations to a full-valence active space are examined and benchmarked. Absolute excitation energies and intensities computed with the GAS-DSRG truncated to second-order in perturbation theory are found to systematically underestimate experimental and reference theoretical values. Third-order perturbative corrections significantly improve the accuracy of GAS-DSRG absolute excitation energies, bringing the mean absolute deviation from experimental values down to 0.32 eV. The ozone molecule and glyoxylic acid are particularly challenging for second-order perturbation theory and are examined in detail to assess the importance of active space truncation and intruder states.  more » « less
Award ID(s):
1900532
PAR ID:
10407070
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
158
Issue:
12
ISSN:
0021-9606
Page Range / eLocation ID:
124112
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The molecules 1,4-cyclohexadiene (unconjugated 1,4-CHD) and 1,3-cyclohexadiene (conjugated 1,3-CHD) both have two double bonds, but these bonds interact in different ways. These molecules have long served as examples of through-bond and through-space interactions, respectively, and their electronic structures have been studied in detail both experimentally and theoretically, with the experimental assignments being especially complete. The existence of Rydberg states interspersed with the valence states makes the quantum mechanical calculation of their spectra a challenging task. In this work, we explore the electronic excitation energies of 1,4-CHD and 1,3-CHD for both valence and Rydberg states by means of complete active space second-order perturbation theory (CASPT2), extended multi-state CASPT2 (XMS-CASPT2), and multiconfiguration pair-density functional theory (MC-PDFT); it is shown by comparison to experiment that MC-PDFT yields the most accurate results. We found that the inclusion of Rydberg orbitals in the active space not only enables the calculation of Rydberg excitation energies but also improves the accuracy of the valence ones. A special characteristic of the present analysis is the calculation of the second moments of the excited-state orbitals. Because we find that the CASPT2 densities agree well with the CASSCF ones and since the MC-PDFT methods gets accurate excitation energies based on the CASSCF densities, we believe that we can trust these moments as far as giving a more accurate picture of the diffuseness of the excited-state orbitals in these prototype molecules than has previously been available. 
    more » « less
  2. A practical ab initio composite method for modeling x-ray absorption and non-resonant x-ray emission is presented. Vertical K-edge excitation and emission energies are obtained from core-electron binding energies calculated with spin-projected ΔHF/ΔMP and outer-core ionization potentials/electron affinities calculated with electron propagator theory. An assessment of the combined methodologies against experiment is performed for a set of small molecules containing second-row elements. 
    more » « less
  3. Accurate modeling of conformational energies is key to the crystal structure prediction of conformational polymorphs. Focusing on molecules XXXI and XXXII from the seventh blind test of crystal structure prediction, this study employs various electronic structure methods up to the level of domain-local pair natural orbital coupled cluster singles and doubles with perturbative triples [DLPNO-CCSD(T1)] to benchmark the conformational energies and to assess their impact on the crystal energy landscapes. Molecule XXXI proves to be a relatively straightforward case, with the conformational energies from generalized gradient approximation (GGA) functional B86bPBE-XDM changing only modestly when using more advanced density functionals such as PBE0-D4, ωB97M-V, and revDSD-PBEP86-D4, dispersion-corrected second-order Møller–Plesset perturbation theory (SCS-MP2D), or DLPNO-CCSD(T1). In contrast, the conformational energies of molecule XXXII prove difficult to determine reliably, and variations in the computed conformational energies appreciably impact the crystal energy landscape. Even high-level methods such as revDSD-PBEP86-D4 and SCS-MP2D exhibit significant disagreements with the DLPNO-CCSD(T1) benchmarks for molecule XXXII, highlighting the difficulty of predicting conformational energies for complex, drug-like molecules. The best-converged predicted crystal energy landscape obtained here for molecule XXXII disagrees significantly with what has been inferred about the solid-form landscape experimentally. The identified limitations of the calculations are probably insufficient to account for the discrepancies between theory and experiment on molecule XXXII, and further investigation of the experimental solid-form landscape would be valuable. Finally, assessment of several semi-empirical methods findsr2SCAN-3c to be the most promising, with conformational energy accuracy intermediate between the GGA and hybrid functionals and a low computational cost. 
    more » « less
  4. X-ray absorption spectroscopy (XAS) is a powerful experimental technique to probe the local order in materials with core electron excitations. Experimental interpretation requires supporting theoretical calculations. For water, these calculations are very demanding and, to date, could only be done with major approximations that limited the accuracy of the calculated spectra. This prompted an intense debate on whether a substantial revision of the standard picture of tetrahedrally bonded water was necessary to improve the agreement of theory and experiment. Here, we report a first-principles calculation of the XAS of water that avoids the approximations of prior work, thanks to recent advances in electron excitation theory. The calculated XAS spectra, and their variation with changes of temperature and/or with isotope substitution, are in good quantitative agreement with experiments. The approach requires accurate quasiparticle wave functions beyond density functional theory approximations, accounts for the dynamics of quasiparticles, and includes dynamic screening as well as renormalization effects due to the continuum of valence-level excitations. The three features observed in the experimental spectra are unambiguously attributed to excitonic effects. The preedge feature is associated with a bound intramolecular exciton, the main-edge feature is associated with an exciton localized within the coordination shell of the excited molecule, and the postedge feature is delocalized over more distant neighbors, as expected for a resonant state. The three features probe the local order at short, intermediate, and longer range relative to the excited molecule. The calculated spectra are fully consistent with a standard tetrahedral picture of water. 
    more » « less
  5. The focal-point approximation can be used to estimate a high-accuracy, slow quantum chemistry computation by combining several lower-accuracy, faster computations. We examine the performance of focal-point methods by combining second-order Møller–Plesset perturbation theory (MP2) with coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] for the calculation of harmonic frequencies and that of fundamental frequencies using second-order vibrational perturbation theory (VPT2). In contrast to standard CCSD(T), the focal-point CCSD(T) method approaches the complete basis set (CBS) limit with only triple-ζ basis sets for the coupled-cluster portion of the computation. The predicted harmonic and fundamental frequencies were compared with the experimental values for a set of 20 molecules containing up to six atoms. The focal-point method combining CCSD(T)/aug-cc-pV(T + d)Z with CBS-extrapolated MP2 has mean absolute errors vs experiment of only 7.3 cm−1 for the fundamental frequencies, which are essentially the same as the mean absolute error for CCSD(T) extrapolated to the CBS limit using the aug-cc-pV(Q + d)Z and aug-cc-pV(5 + d)Z basis sets. However, for H2O, the focal-point procedure requires only 3% of the computation time as the extrapolated CCSD(T) result, and the cost savings will grow for larger molecules. 
    more » « less