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: Benchmarking the semi-stochastic CC( P ; Q ) approach for singlet–triplet gaps in biradicals
We recently proposed a semi-stochastic approach to converging high-level coupled-cluster (CC) energetics, such as those obtained in the CC calculations with singles, doubles, and triples (CCSDT), in which the deterministic CC(P;Q) framework is merged with the stochastic configuration interaction Quantum Monte Carlo propagations [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017)]. In this work, we investigate the ability of the semi-stochastic CC(P;Q) methodology to recover the CCSDT energies of the lowest singlet and triplet states and the corresponding singlet–triplet gaps of biradical systems using methylene, (HFH)−, cyclobutadiene, cyclopentadienyl cation, and trimethylenemethane as examples.  more » « less
Award ID(s):
1763371
PAR ID:
10440304
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
157
Issue:
13
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The energies and geometries of the lowest lying singlet and triplet states of the four diradicals formed by removing two H atoms from thiophene have been characterized. We utilized the highly correlated, multireference methods configuration interaction with single and double excitations with and without the Pople correction for size‐extensivity (MR‐CISD+Q and MR‐CISD) and averaged quadratic coupled cluster theory (MR‐AQCC). CAS (8,7) and CAS (10,8) active spaces involving σ, σ*, π, and π* orbitals were employed along with the cc‐pVDZ and cc‐pVTZ basis sets. The larger active space included the two electrons in the nonbonding sp2hybrid orbital on sulfur. We find that all didehydro isomers exist as planar, stable ground state singlets. The singlet‐triplet (S‐T) adiabatic gaps range from 15 to 25 kcal/mol while the vertical splittings are 21–35 kcal/mol. The 2,3 isomer has the lowest absolute ground state singlet energy and the largest adiabatic and vertical S‐T splitting. The ground states of the 2,3‐, and 2,5‐didehydrothiophene isomers are predicted to exhibit the smallest and largest diradical character, respectively, based on their electronic structures, spin densities and bonding analysis. To our knowledge, no experimental excitation energies of any of the didehydrothiophene isomers are available, and our computed MR‐AQCC/cc‐pVTZ data are believed to be among the most accurate computed results. This extensive study shows a competitive performance between MR‐AQCC and MR‐CISD+Q. 
    more » « less
  2. Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as 𝒪(𝑁9) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the 𝒪(𝑁7) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the 𝒪(𝑁9) CCSD(TQ). Further, we examine the performance of Goodson’s continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a “new” method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol–1 and an L2D of 0.52 kcal mol–1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably “ballpark” how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math. 
    more » « less
  3. null (Ed.)
    A novel pyranine derivative, Et HPTA-OH, was synthesized via the substitution of the anionic sulfonate groups with neutral diethylsulfonamide groups. The photophysical and photochemical properties of Et HPTA-OH were studied using photoluminescence quenching and transient absorption spectroscopy. The singlet state of Et HPTA-OH was found to be highly photoacidic (p K a * = 8.74 in acetonitrile). A series of aniline and pyridine bases were used to investigate excited-state proton transfer (ESPT) from singlet Et HPTA-OH, and rate constants for singlet quenching via ESPT were determined ( k q = 5.18 × 10 9 to 1.05 × 10 10 M −1 s −1 ). Et HPTA-OH was also found to exhibit a long-lived triplet state which reacts through a triplet–triplet annihilation (TTA) process to reform singlet Et HPTA-OH on timescales of up to 80 μs. Detection of ESPT photoproducts on timescales comparable to that of TTA singlet regeneration provides strong evidence for photoacidic behavior stemming from the regenerated singlet Et HPTA-OH. 
    more » « less
  4. We report results from experiments with the quinoline-O2 complex, which was photodissociated using light near 312 nm. Photodissociation resulted in formation of the lowest excited state of oxygen, O2 a 1Δg, which we detected using resonance enhanced multiphoton ionization and velocity map ion imaging. The O2+ ion image allowed for a determination of the center-of-mass translational energy distribution, P(ET), following complex dissociation. We also report results of electronic structure calculations for the quinoline singlet ground state and lowest energy triplet state. From the CCSD/aug-cc-pVDZ//(U)MP2/cc-pVDZ calculations, we determined the lowest energy triplet state to have ππ* electronic character and to be 2.69 eV above the ground state. We also used electronic structure calculations to determine the geometry and binding energy for several quinoline-O2 complexes. The calculations indicated that the most strongly bound complex has a well depth of about 0.11 eV and places the O2 moiety above and approximately parallel to the quinoline ring system. By comparing the experimental P(ET) with the energy for the singlet ground state and the lowest energy triplet state, we concluded that the quinoline product was formed in the lowest energy triplet state. Finally, we found the experimental P(ET) to be in agreement with a Prior translational energy distribution, which suggests a statistical dissociation for the complex. 
    more » « less
  5. Abstract The chemical stability and the low‐lying singlet and triplet excited states of BN‐n‐acenes (n = 1–7) were studied using single reference and multireference methodologies. From the calculations, descriptors such as the singlet‐triplet splitting, the natural orbital (NO) occupations and aromaticity indexes are used to provide structural and energetic analysis. The boron and nitrogen atoms form an isoelectronic pair of two carbon atoms, which was used for the complete substitution of these units in the acene series. The structural analysis confirms the effects originated from the insertion of a uniform pattern of electronegativity difference within the molecular systems. The covalent bonds tend to be strongly polarized which does not happen in the case of a carbon‐only framework. This effect leads to a charge transfer between neighbor atoms resulting in a more strengthened structure, keeping the aromaticity roughly constant along the chain. The singlet‐triplet splitting also agrees with this stability trend, maintaining a consistent gap value for all molecules. The BN‐n‐acenes molecules possess a ground state with monoconfigurational character indicating their electronic stability. The low‐lying singlet excited states have charge transfer character, which proceeds from nitrogen to boron. 
    more » « less