Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “excitations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Exploring Excited State Proton Transfer Dynamics upon Ultraviolet Excitation

Excited state intramolecular proton transfer (ESIPT) has been investigated in two prototypical systems─salicylaldehyde azine (SAA) and 1,5-dihydroxyanthraquinone (DHAQ)─using transient absorption spectroscopy upon ultraviolet excitation into the less studied higher excited (S n ) manifold. Excitation with sub-30 fs pulses and broadband visible probing has allowed for direct measurement of the ESIPT rate. In conjunction with steady-state measurements and TD-DFT calculations, a complete delineation of the ultrafast photophysics has been carried out. In SAA, ESIPT remains ultrafast (∼30 fs), consistent with previous S 1 excitation studies. Coherent vibrational beats maps reveal significant wavelength dependence, however. Theoretical analysis suggests that the observed modes and their intensities in coherent vibrational spectra are modulated by the nature of the electronically excited state. In DHAQ, the first direct observation of ESIPT presents a time-constant of ∼85 fs, and a slower component of 9 ps, akin to previous reports on double-proton transfer systems. Collectively, the results suggest that while the ultrafast ESIPT rate remains largely invariant vis-à-vis the excitation energy, the product yield, as well as accompanying coherent oscillations, may be substantively altered, owing to the existence of alternative decay pathways.

excited states↗
DOE OSTI2025

Determination of proton and neutron contributions to the 0$^+_{𝑔⁡𝑠}$ → 2$^+_{1}$ excitations in 42 Si and 44 S using inelastic proton scattering in inverse kinematics and intermediate-energy Coulomb excitation

Here, we have measured the 0$^+_{𝑔⁡𝑠}$ → 2$^+_1$ transition in the neutron-rich 𝑁 = 28 isotope 42 Si using the probes of intermediate-energy Coulomb excitation and inelastic proton scattering in inverse kinematics at the Facility for Rare Isotope Beams with beam particle rates of ≈ 5 particles/s. The results of these two measurements allowed us to determine 𝑀 𝑛 /𝑀 𝑝 , the ratio of the neutron and proton transition matrix elements for the 0$^+_{𝑔⁡𝑠}$ → 2$^+_1$ transition. In addition, we have measured the 0$^+_{𝑔⁡𝑠}$ → 2$^+_1$ transition in the isotone 44 S using inverse kinematics inelastic proton scattering. By comparing the 44 S proton-scattering result with a recent intermediate-energy Coulomb excitation result on the same transition, we were able to determine 𝑀 𝑛 /𝑀 𝑝 for the 0$^+_{𝑔⁡𝑠}$ → 2$^+_1$ transition in this nucleus as well. This work strengthens the evidence that 42 Si has a stable quadrupole deformation in its ground state and that 44 S does not. Both conclusions are further supported by shell-model calculations carried out with the FSU interaction.

39 ≤ A ≤ 58↗
DOE OSTI2025

Rotational excitation and de-excitation of magnesium mono-sulphide (MgS) by He collisions

ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.

Hendaoui, Hamza (ORCID:0000000218641872)↗
DOE OSTI2024

Exciting DeePMD: Learning excited-state energies, forces, and non-adiabatic couplings

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. Furthermore, the efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH 2 NH 2 + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Excited State Covalency, Dynamics, and Photochemistry of Square Planar Ni-Thiolate Complexes Revealed by Ultrafast X-ray Absorption

Highly covalent Ni bis(dithiolene) and related complexes provide an ideal platform for investigating the effects of metal–ligand orbital hybridization on excited state character and dynamics. In particular, we focus on the ligand field excited states that dominate the photophysics of first-row transition metal complexes. Here, we investigate if they can be significantly delocalized off the metal center, possibly yielding photochemical reactivity more similar to charge transfer excited states than metal-centered ligand field excited states. Here [Ni(mpo) 2 ] (mpo = 2-mercaptopyridine- N-oxide) provides a representative example for the larger chemical class and is an active electro- and photo-catalyst for proton reduction. A detailed characterization of the excited state electronic structure, dynamics, and photochemistry of [Ni(mpo) 2 ] is presented based on ultrafast transient X- ray absorption spectroscopy at the Ni and S 1s core absorption K-edges. By comparing the ultrafast Ni K-edge absorption to ab initio calculations, we identify an excited state relaxation mechanism where an initial ligand-to-metal charge transfer excitation results in both excited state electron transfer (generating a catalytically relevant reduced photoproduct [Ni(mpo) 2 ] − ) and relaxation to a pseudo-tetrahedral triplet ligand field excited state. From the ultrafast S K-edge absorption, the ligand field excited state is found to be highly delocalized onto the thiolate ligands and a tetrahedral structural distortion is shown to substantially influence the degree of delocalization. The results identify a significant structural coordinate to target when aiming to control the excited state covalency in square planar complexes.

Lim, Hyeongtaek [SLAC National Accelerator Laborat↗
DOE OSTI2026

Highly-destabilized ligand field excited states of iron carbene complexes and their relation to charge transfer state lifetimes

Lifetimes of photoexcited charge transfer (CT) states in transition metal chromophores are influenced by low-lying ligand field (LF) excited states, especially for 3d metal complexes. To manipulate interactions between LF and CT states, it is important to be able to control LF excited state energies using tunable synthetic variables. In this report, we use Fe 2p3d L 3 -edge resonant inelastic X-ray scattering (RIXS) to measure LF excited state energies of three homoleptic iron chromophores coordinated by strong-field N-heterocyclic carbenes (NHCs). We investigate the effect of oxidation state and ligand scaffold on LF energies and covalency parameters. A cyclometalated bis(NHC) ligand affords both high LF excited state energies (and thus high 10 Dq) as well as high metal–ligand covalency compared to other iron complexes with very strong-field ligands. However, for the set of complexes investigated, we do not observe meaningful correlation between the LF excited state energies and the CT excited state lifetimes. These results illustrate that targeting long-lived CT excited states necessitates control of multiple molecular excited state properties, with destabilization of the LF excited state energies proving necessary, but insufficient, to control the CT excited state lifetime in Fe carbene complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Analytic Gradients for Equation-of-Motion Coupled Cluster with Single, Double, and Perturbative Triple Excitations

Understanding the process of molecular photoexcitation is crucial in various fields, including drug development, materials science, photovoltaics, and more. The electronic vertical excitation energy is a critical property, for example in determining the singlet-triplet gap of chromophores. However, a full understanding of excited-state processes requires additional explorations of the excited-state potential energy surface and electronic properties, which is greatly aided by the availability of analytic energy gradients. Owing to its robust high accuracy over a wide range of chemical problems, equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) is a powerful method for predicting excited state properties, and the implementation of analytic gradients of many EOM-CCSD (excitation energies, ionization potentials, electron attachment energies, etc.) along with numerous successful applications high- lights the flexibility of the method. In specific cases where a higher level of accuracy is needed or in more complex electronic structures, the inclusion of triple excitations becomes essential, for example, in the EOM-CCSD* approach of Saeh and Stanton. In this work, we derive and implement for the first time the analytic gradients of EOMEE-CCSD*, which also provides a template for analytic gradients of related ex- cited state methods with perturbative triple excitations. Here, the capabilities of analytic EOMEE-CCSD* gradients are illustrated by several representative examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Benchmarking third-order cluster perturbation theory for electronically excited states

In this study, we investigate the reliability of cluster perturbation (CP) theory applied to the calculation of electronically excited states through a comprehensive benchmark. In CP theory, perturbative corrections are added to the properties of a parent excitation space, which converge toward the properties of a target excitation space. For the CPS(D-n) model, perturbative corrections through order n are added to the coupled cluster singles (CCS) excitation energies to target the coupled cluster singles and doubles (CCSD) excitation energies. Through a comparative analysis of excitation energy calculations across a diverse set of molecules and wavefunction methods, we present a comprehensive evaluation of the accuracy of the third-order CPS(D) model, CPS(D-3), in calculating excitation energies. Further, our findings demonstrate that CPS(D-3) is a reliable alternative to established methods, particularly CCSD, while systematically overestimating the excitation energies compared to high-level coupled cluster methods such as CC3. These results highlight the strengths and limitations of CPS(D-3), as well as the promising directions for its future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S 1 and S 2 . Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S 1 –S 6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S 0 –S 6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. In conclusion, we compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Rishi, Varun↗
DOE OSTI2026

Luminescent Zr Complexes with Long-lived Intraligand Charge Transfer Excited States

Understanding and establishing design principles to tune the excited states of earth-abundant sensitizers is crucial for identifying new photosensitizers for sustainable technological advancements. Here, we report a series of three bis Zr (IV) complexes of tridentate, dianionic ligands incorporating two phenoxide donors and aza arene acceptors of increasing electron affinity (pyridine < pyrimidine < pyrazine ) that are air and water-stable. These complexes emit via thermally assisted delayed fluorescence from an intraligand charge transfer excited state with varying metal contributions. The electronic structural changes from different acceptors vary the excited state character and the photophysical properties significantly between the complexes. The variation in photophysics that includes emission lifetime (76 μs to 265 ns), intersystem crossing (ISC) lifetime ( 390 – 290 ps), and the energy difference between the singlet and triplet excited states (∆EST. 170 -100 meV). We observe solvent-independent ISC rates when the excited state has metal contributions (acceptor = pyridine), and the ISC rates vary significantly with solvent polarity when the excited state is an intra-ligand charge transfer state (acceptor = pyrimidine and pyrazine). Transient absorption measurements provided the basis spectra and verified the metal contribution in the excited state. This work provides a basis for developing new sensitizers based on Zr (IV) by providing design principles that can be used to modulate the character of the excited state and its photophysics.

Nattikallungal, Thabassum A. [University of Southe↗
DOE OSTI2025

Excitation Laser Energy Dependence of the Gap-Mode TERS Spectra of WS 2 and MoS 2 on Silver

In this work we present a systematic study of the dependence of the gap mode tip-enhanced Raman scattering (TERS) response of the mono- and bi-layer WS 2 and MoS 2 on silver as a function of the excitation laser energy in a broad spectral range from 473nm to 830 nm. For this purpose, we collected consecutive TERS maps of the same area in the sample containing mono-and bi-layer regions with the same TERS probe with 6 different excitation lasers. To decrease the number of collected TERS maps, we used for the first time to the best of our knowledge, concurrent excitation and collection with two lasers simultaneously. We found that the E 2g /A 1g peak intensity ratio for the bilayer WS 2 @Ag and the ratio of the A’/A 1g peak intensity of the out-of-plane mode for the mono- and the bilayer, change in a significantly non-monotonous way as the excitation laser energy is swept from 1.58 eV to 2.62 eV. The former ratio increases at energies corresponding to A and B excitons (~2.0 eV and 2.4 eV correspondingly) in bilayer WS 2 . The absolute intensity of the A’ peak in the monolayer, and correspondingly the A’/A 1g ratio, is surprisingly high at lower excitation energies, but dips dramatically at the energy corresponding to the A exciton, being restored partially in between A and B excitons, but still showing the descending trend as the excitation laser energy increases. A somewhat similar picture was observed in mono- and bi-layers of MoS 2 @Ag, though the existing set of excitation lasers did not match the excitonic profile of this material as nicely as for the case of WS 2 . We attribute the observed behavior to the presence of intermediate (Fano resonance) or strong (Rabi splitting) coupling between the excitons in transition metal dichalcogenides (TMDs) and the plasmons in the tip-substrate nanocavity. This is akin to the so-called Fano (Rabi) transparency experimentally observed in far field scattering from transition metal dichalcogenides between two plasmonic metals. The possibility of the formation of intermediate/strong coupling between the excitonic resonances in TMDs and the nanocavity re-evaluates the role of various resonances in the gap-mode TERS and should become an important factor to be considered by TERS practitioners during planning the experiments. Lastly, based on observed phenomena and its explanation, we propose the “ideal” substrate for efficient TERS and tip enhanced photoluminescence (TEPL) measurements.

Raman spectroscopy↗
DOE OSTI2025

Excited-State Densities from Time-Dependent Density Functional Response Theory

While the variational principle for excited-state energies leads to a route to obtaining excited-state densities from time-dependent density functional theory, relatively little attention has been paid to the quality of the resulting densities in real space obtained with different exchange-correlation functional approximations or how nonadiabatic approximations developed for energies of states of double-excitation character perform for their densities. Here we derive an expression directly in real space for the excited-state density, which includes the case of nonadiabatic kernels and consequently is able, for the first time, to yield densities of states of double-excitation character. Under some well-defined simplifications, we compare the performance of the local-density approximation and exact-exchange approximation, which are in a sense at the opposite extremes of the fundamental functional approximations, on local and charge-transfer excitations in one-dimensional model systems and show that the dressed Time-Dependent Density Functional Theory (TDDFT) approach gives good densities of double excitations.

approximation↗
DOE OSTI2026

Excitation and tunneling spectra of a fractional quantum Hall system in the thin-cylinder limit

Here, the excitations of fractional quantum Hall effect (FQHE) states have been largely inaccessible to experimental probes until recently. New electron scanning tunneling microscopy (STM) results from Hu et al. [Nat. Phys. 21, 716 (2025)] show promise in detecting and identifying these excited states via the local density of states (LDOS) spectrum. On a torus, there exists a mapping from the lowest Landau level states to a 1D lattice with a Hamiltonian that features dipole moment conservation. In this work, we apply perturbation theory starting from the thin-cylinder limit (𝐿 𝑥 → ∞, 𝐿 𝑦 < 𝑙 𝐵 for torus dimensions 𝐿 𝑥 and 𝐿 𝑦 and magnetic length 𝑙 𝐵 ) to obtain an analytical approach to the low-lying neutral and charged excitations of the 𝜈 = 1/3 FQHE state. Notably, in the thin cylinder, we can systematically enumerate all the low-lying excitations by the patterns of “dipoles” formed by the electron occupation pattern on the 1D lattice. We find that the thin-cylinder limit predicts a significant dispersion of the low-lying neutral excitations but sharpness of the LDOS spectra, which measure charged excitations. We also discuss connections between our work and several different approaches to the FQHE STM spectra, including those using the composite fermion theory. Numerical exact diagonalization beyond the thin-cylinder limit suggests that the energies of charged excitations remain largely confined to a narrow range of energies, which in experiments might appear as a single peak.

Adhidewata, Jyesta M. [University of California, B↗
DOE OSTI2025

Nonpigmented PsbR is involved in the integrity of excitation landscape in higher plant photosystem II, a case study in Arabidopsis thaliana and a mutant

PsbR is a nonpigmented 10 kDa protein in Photosystem II (PSII) in algae and plants. A recent structural study clarified its enigmatic structural location in a Photosystem II megacomplex that has baffled the community for more than four decades. Our current study interrogates whether absence of PsbR affects the overall dynamics of excitation energy migration within light harvesting complexes (LHC) and PSII super assemblies using highly-active PSII membrane particles, so-called BBY particles, isolated from a PsbR deletion mutant (ΔPsbR) of Arabidopsis thaliana. A femto-second (fs)-time-resolved transient absorption experimentation recorded at 77 K with selective excitation of Chl b which is exclusively present in LHCs enabled us to resolve the temporal differences in LHC→LHC and LHC→PSII excitation energy transfer steps. By applying specific target spectro-kinetic models to the transient absorption datasets, we demonstrated that the time constants of Chl a LHC → Chl a LHC excitation transfer significantly elongates in the ΔPsbR LHC-PSII particles, suggestive of the decreased aggregation level of photosynthetic proteins in the mutant. These findings highlight excitation energy transfer integrity in LHC-PSII assembly is not only determined by the pigmented light-harvesting complexes, but also synergistically by the nonpigmented PSII components. Furthermore, the disturbed integrity in dynamics of excitation energy transfer pathway within LHC-PSII supercomplex is discussed in the context of the altered LHC-PSII megacomplexes type I and II architectures which result from the absence of the PsbR protein in higher plant PSII.

Chlorophyll↗
DOE OSTI2025

Leveraging Intramolecular π-Stacking to Access an Exceptionally Long-Lived 3 MC Excited State in an Fe(II) Carbene Complex

The ability to manipulate excited-state decay cascades using molecular structure is essential to the application of abundant-metal photosensitizers and chromophores. Ligand design has yielded some spectacular results elongating charge-transfer excited state lifetimes of Fe(II) coordination complexes, but triplet metal-centered ( 3 MC) excited states - recently demonstrated to be critical to the photoactivity of isoelectronic Co(III) polypyridyls - have to date remained elusive, with temporally isolable examples limited to the picosecond regime. Here, with this report, we show how strong-field donors and intramolecular π-stacking can conspire to stabilize a long-lived 3 MC excited state for a remarkable 4.1 ± 0.3 ns in fluid solution at ambient temperature. Analysis of variable-temperature time-resolved absorption data using theoretical models ranging from Arrhenius to semiclassical Marcus theory, combined with computational modeling and X-ray crystallography, reveal a Jahn−Teller stabilized excited state with a high activation barrier for ground-state recovery. The net result is a chromophore with a 3 MC excited-state lifetime that is orders of magnitude longer than anything yet observed for an Fe(II) complex.

carbene compounds↗
DOE OSTI2025

Search for excited tau leptons in the ττγ final state in proton-proton collisions at $\sqrt{\text{s}}$ = 13 TeV

Results are presented for a test of the compositeness of the heaviest charged lepton, τ, using data collected by the CMS experiment in proton-proton collisions at a center-of-mass energy of 13 TeV at the CERN LHC. The data were collected in 2016–2018 and correspond to an integrated luminosity of 138 fb −1 . This analysis searches for tau lepton pair production in which one of the tau leptons is produced in an excited state and decays to a ground state tau lepton and a photon. The event selection consists of two isolated tau lepton decay candidates and a high-energy photon. The mass of the excited tau lepton is reconstructed using the missing transverse momentum in the event, assuming the momentum of the neutrinos from each tau lepton decay are aligned with the visible decay products. No excess of events above the standard model background prediction is observed. This null result is used to set lower bounds on the excited tau lepton mass. For a compositeness scale Λ equal to the excited tau lepton mass, excited tau leptons with masses below 4700 GeV are excluded at 95% confidence level; for Λ = 10 TeV this exclusion is set at 2800 GeV. This is the first experimental result covering this production and decay process in the excited tau mass range above 175 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Here, our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗