Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled cluster”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

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

Study of self-interaction errors in density functional predictions of dipole polarizabilities and ionization energies of water clusters using Perdew-Zunger and locally scaled self-interaction corrected methods

Here, we studied the effects of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modelled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz. the local spin density approximation (LDA), the Perdew-Burke-Ernzherof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Lowdin orbital SIC (FLO-SIC) framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.58 Bohr 3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but over-corrects, resulting in underestimated polarizabilities in SIC-LDA, -PBE, and -SCAN. Finally, we applied a recently proposed local-scaling SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute error of 0.08 Bohr 3 for LSIC-LDA and 0.06 Bohr 3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as an absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAE of 0.4 eV for LSIC-LDA and 0.04 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies while the LSIC-PBE ionization energies are more accurate than reported GW results.

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

Vertical detachment energies of ammonia cluster anions using self-interaction-corrected methods

Systems with weakly bound extra electrons impose great challenges to semilocal density functional approximations (DFAs), which suffer from self-interaction errors. Small ammonia clusters are one such example of weakly bound anions where the extra electron is weakly bound. We applied two self-interaction correction (SIC) schemes, viz., the well-known Perdew–Zunger and the recently developed locally scaled SIC (LSIC) with the local spin density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation (GGA), and the SCAN meta-GGA functionals to calculate the vertical detachment energies (VDEs) of small ammonia cluster anions (NH3)n–. Our results show that the LSIC significantly reduces the errors in calculations of VDE with LSDA and PBE-GGA functionals leading to better agreement with the reference values calculated with coupled cluster singles and doubles with perturbative triples [CCSD(T)]. Accurate prediction of VDE as an absolute of the highest occupied molecular orbital (HOMO) is challenging for DFAs. Our results show that VDEs estimated from the negative of HOMO eigenvalues with the LSIC-LSDA and Perdew–Zunger SIC-PBE are within 11 meV of the reference CCSD(T) results. Here, the LSIC method performs consistently well for the VDE estimates, from both the total energy differences and the absolute HOMO eigenvalues.

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

Towards complete assignment of the infrared spectrum of the protonated water cluster H + (H 2 O) 21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H + (H 2 O) n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H + (H 2 O) 21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

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

Quantitative Account of the Bonding Properties of a Rubredoxin Model Complex [Fe(SCH3)4]q, q = -2, -1, +2, +3

Iron-sulfur clusters play important roles in biology as parts of electron transfer chains and catalytic cofactors. Here, we report a detailed computational analysis of a structural model of the simplest natural iron-sulfur cluster of rubredoxin and its cationic counterparts. Specifically, we report results for the ground and low-lying electronically excited states of the complex [Fe(SCH3)4]2-/1-/2+/3+, using Multi-Reference (CASSCF, MRCISD), and Coupled Cluster [CCSD(T)] methodology in order to provide accurate adiabatic reduction energies, dissociation energies and insights into the bonding analysis. The nature of the Fe-S chemical bond and the magnitude of the ionization potentials in the anionic and cationic [Fe(SCH3)4] complexes offer a physical rationale for the relative stabilization, structure and speciation of these complexes. Anionic and cationic complexes present different types of chemical bonds: prevalently ionic in [Fe(SCH3)4]2-/1- complexes and covalent in [Fe(SCH3)4]2+/3+ complexes. The ionic bonds result in an energy gain for the transition [Fe(SCH3)4]2-®[Fe(SCH3)4]- (i.e., FeII®FeIII) of 1.5 eV, while the covalent bonds result in an energy loss for the transition [Fe(SCH3)4]2+®[Fe(SCH3)4]3+ of 16.6 eV, almost half of the IP of Fe2+. The ionic vs covalent bond character influences the Fe-S bond strength and length, i.e., ionic Fe-S bonds are longer than covalent ones by about 0.2 Å (for FeII) and 0.04 Å (for FeII). Finally, the average Fe-S heterolytic bond strength is 6.7 eV (FeII) and 14.6 (FeIII) eV at the RCCSD(T) level of theory.

Tzeli, Demeter↗
DOE OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Interplay of Open-Shell Spin-Coupling and Jahn–Teller Distortion in Benzene Radical Cation Probed by X-ray Spectroscopy

Here we report a theoretical investigation and elucidation of the X-ray absorption spectra of neutral benzene and of the benzene cation. The generation of the cation by multiphoton ultraviolet (UV) ionization and the measurement of the carbon K-edge spectra of both species using a table-top high-harmonic generation source are described in the companion experimental paper. We show that the 1sC → π transition serves as a sensitive signature of the transient cation formation, as it occurs outside of the spectral window of the parent neutral species. Moreover, the presence of the unpaired (spectator) electron in the π-subshell of the cation and the high symmetry of the system result in significant differences relative to neutral benzene in the spectral features associated with the 1sC → π* transitions. High-level calculations using equation-of-motion coupled-cluster theory provide the interpretation of the experimental spectra and insight into the electronic structure of benzene and its cation. The prominent split structure of the 1sC → π* band of the cation is attributed to the interplay between the coupling of the core → π* excitation with the unpaired electron in the π-subshell and the Jahn–Teller distortion. The calculations attribute most of the splitting (~1–1.2 eV) to the spin coupling, which is visible already at the Franck–Condon structure, and we estimate the additional splitting due to structural relaxation to be around ~0.1–0.2 eV. These results suggest that X-ray absorption with increased resolution might be able to disentangle electronic and structural aspects of the Jahn–Teller effect in the benzene cation.

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

An Ab Initio Electronic Structure Investigation of the Ground and Excited States of ScH + , YH + , and LaH +

Multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), coupled-cluster singles, doubles, and perturbative triples [CCSD(T)], and frozen-core full configuration interaction (fcFCI) calculations were carried out using large, correlation-consistent basis sets to investigate the excited states of the Sc atom and the spin–free and spin–orbit coupled potential energy profiles, energetics, spectroscopic constants, and electron populations of low-lying states of MH + (M = Sc, Y, La). The core electron correlation effects, complete basis set effects, and spin–orbit coupling effects were also evaluated. The first four electronic states of all MH + are 1 2 Δ, 1 2 Σ + , 1 2 Π, and 2 2 Σ + with 1σ 2 1δ 1 , 1σ 2 2σ 1 , 1σ 2 1π 1 , and 1σ 2 3σ 1 single-reference electron configurations, respectively. These states of MH + can be represented by the M 2+ H – ionic structure. The ground states of ScH + , YH + , and LaH + are 1 2 Δ 3/2 , 1 2 Σ + 1/2 , and 1 2 Δ 3/2 with 55.45, 60.54, and 62.34 kcal/mol bond energies, respectively. The core electron correlation was found to be vital for gaining accurate predictions on the ground and excited state properties of MH + . The spin–orbit coupling effects are minor for ScH + but become substantial moving to YH + and LaH + . Overall, the results of this work are in good agreement with the limited set of experimental findings of MH + available in the literature and will be of use for future investigations. Furthermore, the theoretical approaches, findings, and trends reported here are expected to aid studies of similar species.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2023

SPT-CL J2215−3537: A Massive Starburst at the Center of the Most Distant Relaxed Galaxy Cluster

Abstract We present the discovery of the most distant, dynamically relaxed cool core cluster, SPT-CL J2215−3537 (SPT2215), and its central brightest cluster galaxy (BCG) at z = 1.16. Using new X-ray observations, we demonstrate that SPT2215 harbors a strong cool core with a central cooling time of 200 Myr (at 10 kpc) and a maximal intracluster medium cooling rate of 1900 ± 400 M ⊙ yr −1 . This prodigious cooling may be responsible for fueling the extended, star-forming filaments observed in Hubble Space Telescope imaging. Based on new spectrophotometric data, we detect bright [O ii ] emission in the BCG, implying an unobscured star formation rate (SFR) of 320 − 140 + 230 M ⊙ yr −1 . The detection of a weak radio source (2.0 ± 0.8 mJy at 0.8 GHz) suggests ongoing feedback from an active galactic nucleus (AGN), though the implied jet power is less than half the cooling luminosity of the hot gas, consistent with cooling overpowering heating. The extreme cooling and SFR of SPT2215 are rare among known cool core clusters, and it is even more remarkable that we observe these at such high redshift, when most clusters are still dynamically disturbed. The high mass of this cluster, coupled with the fact that it is dynamically relaxed with a highly isolated BCG, suggests that it is an exceptionally rare system that must have formed very rapidly in the early universe. Combined with the high SFR, SPT2215 may be a high- z analog of the Phoenix cluster, potentially providing insight into the limits of AGN feedback and star formation in the most massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

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

Emulating ab initio computations of infinite nucleonic matter

We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Threshold Ionization Spectroscopy and Theoretical Calculations of LnO (Ln = La and Ce)

Mass-analyzed threshold ionization (MATI) spectroscopy was used to measure the vibronic spectra of LnO (Ln = La and Ce). Single-reference coupled cluster and relativistic multireference configuration calculations were carried out to compare with the measured vibronic energies. The spectrum of LaO displays a single vibronic band system, while that of CeO shows multiple ones. The ionization energies of LaO and CeO are measured as 5.2446 (6) and 5.3332(6) eV, respectively, which are a hundredfold improvement over the literature values. The vibrational energies of the neutral molecule and corresponding ion reveal the charge effect on the metal-oxygen bond of both species. The single band system in the spectrum of LaO arises from the transition of the ground state of the neutral molecule with the La(6s 1 )O(2p 6 ) valence configuration to the ground state of the singly charged ion with the La(6s o )O(2p 6 ) configuration. Here, the multiple band systems in the spectrum of CeO are attributed to the spin-orbit coupling for the Ce(4f 1 6s 1 )O(2p 6 ) configuration of the neutral molecule and an excited state for the Ce(4f 1 )O(2p 6 ) configuration of the ion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

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

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗