Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ultrafast Dynamics of Excited Electronic States in Nitrobenzene Measured by Ultrafast Transient Polarization Spectroscopy

We investigate ultrafast dynamics of the lowest singlet excited electronic state in liquid nitrobenzene using ultrafast transient polarization spectroscopy, extending the well-known technique of optical Kerr effect spectroscopy to excited electronic states. The third-order nonlinear response of the excited molecular ensemble is measured using a pair of femtosecond pulses following a third femtosecond pulse that populates the S 1 excited state. By measuring this response, which is highly sensitive to details of the excited state character and structure, as a function of time delays between the three pulses involved, we extract the dephasing time of the wave packet on the excited state. The dephasing time, measured as a function of time delay after pump excitation, shows oscillations indicating oscillatory wave packet dynamics on the excited state. Finally, from the experimental measurements and supporting theoretical calculations, we deduce that the wave packet completely leaves the S 1 state potential energy surface after three traversals of the intersystem crossing between the singlet S 1 and triplet T 2 states.

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

Formation of Phenanthrene via Recombination of Indenyl and Cyclopentadienyl Radicals: A Theoretical Study

This work presents quantum chemical G3(MP2,CC)//B2PLYPD3/6-311G(d,p) calculations of the potential energy surface for the indenyl (C 9 H 7 ) + cyclopentadienyl (C 5 H 5 ) reaction followed by unimolecular decomposition of the C 14 H 11 radicals formed as the primary products, as well as the Rice–Ramsperger–Kassel–Marcus master equation (RRKM-ME) calculations to predict temperature- and pressure-dependent reaction rate constants and product branching ratios. The reaction begins with the barrierless recombination of indenyl and cyclopentadienyl forming a C 14 H 12 molecule with a new C–C bond connecting two five-membered rings, which subsequently dissociates to C 14 H 11 radicals by H losses. The primary products of the C 9 H 7 + C 5 H 5 → C 14 H 11 + H reaction can directly decompose by another H loss to benzofulvalene, and this pathway is most favorable in terms of the entropy factor and hence is preferable at higher temperatures. Otherwise, the initial C 14 H 11 isomers can undergo significant structural rearrangements before eliminating an H atom and producing phenanthrene, anthracene, or benzoazulenes, among which the formation of phenanthrene via the “spiran” pathway is clearly preferred. The calculated barriers along the computed favorable dissociation pathways are relatively low, in the ~30–40 kcal/mol range, making the C 14 H 11 radicals unstable at temperatures above 1000–1250 K at 1 atm. The results of RRKM-ME calculations show that, under typical combustion conditions, the decomposition of the C 14 H 11 radicals predominantly leads to benzofulvalene. However, the latter can be rapidly converted to phenanthrene via H-assisted isomerization with the rate constant for the benzofulvalene + H → phenanthrene + H reaction being close to 10 –11 cm 3 molecule –1 s –1 at 1000–1500 K and 1 atm. Here, the results provide further support for the hypothesis that recombination of two π radicals containing five-membered rings can lead to a growth of PAH with the formation of two fused six-membered rings, but the reaction mechanism may not be direct and is likely to involve two consecutive H atom losses leading to a fulvalene-like product, with subsequent H-assisted isomerization of the latter to a benzenoid PAH.

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

Electrostatics, Charge Transfer, and the Nature of the Halide–Water Hydrogen Bond

Binary halide–water complexes X – (H 2 O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X – (H 2 O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C 2v -symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F – (H 2 O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X – (H 2 O) and provides a driving force for the formation of quasi-linear X··· H–O hydrogen bonds. Here, charge-transfer energies correlate well with measured O–H vibrational redshifts for the halide–water complexes and also for OH – (H 2 O) and NO 2 – (H 2 O), providing some indication of a general mechanism.

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

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

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

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

Theoretical Study of the Phenoxy Radical Recombination with the O( 3 P) Atom, Phenyl plus Molecular Oxygen Revisited

Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.

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

Reactions of U + with H 2 , D 2 , and HD Studied by Guided Ion Beam Tandem Mass Spectrometry and Theory

The kinetic energy-dependent reactions of the atomic actinide uranium cation (U + ) with H 2 , D 2 , and HD were examined by guided ion beam tandem mass spectrometry. An average 0 K bond dissociation energy of D 0 (U + – H) = 2.48 ± 0.06 eV is obtained by analysis of the endothermic product ion cross sections. Quantum chemistry calculations were performed for comparison with experimental thermochemistry, including high-level CASSCF–CASPT2–RASSI calculations of the spin–orbit corrections. CCSD(T) and the CASSCF levels show excellent agreement with experiment, whereas B3LYP and PBE0 slightly overestimate and the M06 approach badly underestimates the bond energy for UH + . Theory was also used to investigate the electronic structures of the reaction intermediates and potential energy surfaces. The experimental product branching ratio for the reaction of U + with HD indicates that these reactions occur primarily via a direct reaction mechanism, despite the presence of a deep-well for UH 2 + formation according to theory. The reactivity and hydride bond energy for U + are compared with those for transition metal, lanthanide, and actinide cations, and periodic trends are discussed. Furthermore, these comparisons suggest that the 5f electrons on uranium are largely core and uninvolved in the reactive chemistry.

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

Theoretical Study of the Reaction of the Methylidyne Radical (CH; X 2 Π) with 1-Butyne (CH 3 CH 2 CCH; X 1 A')

In this study, ab initio CCSD(T)-F12/cc-pVTZ-f12//ωB97X-D/6-311G(d,p) + ZPE(ωB97X-D/6-311G(d,p)) calculations were carried out to unravel the area of the C 5 H 7 potential energy surface accessed by the reaction of the methylidyne radical with 1-butyne. The results were utilized in Rice-Ramsperger-Kassel-Marcus calculations of the product branching ratios at the zero pressure limit. The preferable reaction mechanism has been shown to involve (nearly) instantaneous decomposition of the initial reaction adducts, whose structures are controlled by the isomeric form of the C 4 H 6 reactant. If CH adds to the triple CΞC bond in the entrance reaction channel, the reaction is predicted to predominantly form the methylenecyclopropene + methyl (CH 3 ) and cyclopropenylidene + ethyl (C 2 H 5 ) products roughly in a 2:1 ratio. CH insertion into a C-H bond in the methyl group of 1-butyne is anticipated to preferentially form ethylene + propargyl (C 3 H 3 ) by the C-C bond β-scission in the initial complex, whereas CH insertion into C-H of the CH 2 group would predominantly produce vinylacetylene + methyl (CH 3 ) also by the C-C bond β-scission in the adduct. The barrierless and highly exoergic CH + 1-butyne reaction, facile in cold molecular clouds, is not likely to lead to the carbon skeleton molecular growth but generates C 4 H 4 isomers methylenecyclopropene, vinylacetylene, and 1,2,3-butatriene and smaller C 2 and C 3 hydrocarbons such as methyl, ethyl, and propargyl radicals, ethylene, and cyclopropenylidene.

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

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

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

Recent Advances toward Efficient Calculation of Higher Nuclear Derivatives in Quantum Chemistry

In this article, we provide an overview of state-of-the-art techniques that are being developed for efficient calculation of second and higher nuclear derivatives of quantum mechanical (QM) energy. Calculations of nuclear Hessians and anharmonic terms incur high costs and memory and scale poorly with system size. Three emerging classes of methods—machine learning (ML), automatic differentiation (AD), and matrix completion (MC)—have demonstrated promise in overcoming these challenges. We illustrate studies that employ unsupervised ML methods to reduce the need for multiple Hessian calculations in dynamics simulations and those that utilize supervised ML to construct approximate potential energy surfaces and estimate Hessians and anharmonic terms at reduced cost. By extension, if electronic structure operations could be written in a manner similar to functions underlying ML methods, rapid differentiation or AD routines can be employed to inexpensively calculate higher arbitrary-order derivatives. While ML approaches are typically black-box, we describe methods such as compressed sensing (CS) and MC, which explicitly leverage problem-specific mathematical properties of higher derivatives such as sparsity and low-rank, to complete higher derivative information using only a small, incomplete sample. The three classes of methods facilitate reliable predictions of observables ranging from infrared spectra to thermal conductivity and constitute a promising way forward in accurately capturing otherwise intractable higher-order responses of QM energy to nuclear perturbations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Radical–Radical Reactions in Molecular Weight Growth: The Phenyl + Propargyl Reaction

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C 6 H 5 ) with propargyl radical (H 2 CCCH) provides an important prototype for radical–radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300–1000 K and pressure range of 4–10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C 9 H 8 and C 9 H 7 + H product channels and report experimental isomer-resolved product branching fractions for the C 9 H 8 product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. Here, these ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical–radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C 9 H 7 + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C 9 H 7 ) recombination to indene and H-assisted isomerization that converts less stable C 9 H 8 isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

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

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. Here, in this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S 1 /S 0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

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

Attosecond Probing of Coherent Vibrational Dynamics in CBr 4

A coherent vibrational wavepacket is launched and manipulated in the symmetric stretch (a1) mode of CBr 4 , by impulsive stimulated Raman scattering (ISRS) from nonresonant 400 nm laser pump pulses with various peak intensities on the order of tens of 10 12 W/cm 2 . Extreme ultraviolet (XUV) attosecond transient absorption spectroscopy (ATAS) records the wavepacket dynamics as temporal oscillations in XUV absorption energy at the bromine M 4,5 3d 3/2,5/2 edges around 70 eV. The results are augmented by nuclear timedependent Schrödinger equation simulations. Slopes of the (Br 3d 3/2,5/2 ) –1 10a 1 * core-excited state potential energy surface (PES) along the a 1 mode are calculated to be –9.4 eV/Å from restricted open-shell Kohn–Sham calculations. Using analytical relations derived for the small-displacement limit and the calculated slopes of the core-excited state PES, a deeper insight into the vibrational dynamics is obtained by retrieving the experimental excursion amplitude of the vibrational wavepacket and the amount of population transferred to the vibrational first-excited state as a function of pump-pulse peak intensity. Experimentally, the results show that XUV ATAS is capable of resolving oscillations in the XUV absorption energy on the order of a few to tens of meV with tens of femtosecond time precision. This corresponds to change in C–Br bond length on the order of 10 –4 to 10 –3 Å. The results and the analytic relationships offer a clear physical picture, on multiple levels of understanding, of how the pump-pulse peak intensity controls the vibrational dynamics launched by nonresonant ISRS in the small-displacement limit.

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

Identification of the Elusive Methyl-Loss Channel in the Crossed Molecular Beam Study of Gas-Phase Reaction of Dicarbon Molecules (C 2 ; X 1 Σ g + /a 3 Π u ) with 2-Methyl-1,3-butadiene (C 5 H 8 ; X 1 A')

The crossed molecular beams technique was utilized to explore the reaction of dicarbon C 2 (X 1 Σ g + /a 3 Π u ) with 2-methyl-1,3-butadiene (isoprene, CH 2 C(CH 3 )CHCH 2 ; X 1 A') at a collision energy of 28 ± 1 kJ mol⁻¹ using a supersonic dicarbon beam generated via photolysis (248 nm) of helium-seeded tetrachloroethylene (C 2 Cl 4 ). Here, experimental data combined with previous ab initio calculations provide evidence of the detection of the hitherto elusive methyl elimination channels leading to acyclic resonantly stabilized hexatetraenyl radicals: 1,2,4,5-hexatetraen-3-yl (CH 2 CC•CHCCH 2 ) and/or 1,3,4,5-hexatetraen-3-yl (CH 2 CHC•CCCH 2 ). These pathways are exclusive to the singlet potential energy surface, with the reaction initiated by the barrierless addition of dicarbon to one of the carbon-carbon double bonds in the diene. In combustion systems, both hexatetraenyl radicals can isomerize to the phenyl radical (C 6 H 5 ) through a hydrogen atom assisted isomerization – the crucial reaction intermediate and molecular mass growth species step toward the formation of polycyclic aromatic hydrocarbons (PAHs) and soot.

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

A Combined Crossed Molecular Beam and Theoretical Investigation of the Elementary Reaction of Tricarbon (C 3 (X 1 Σ g + )) with Diacetylene (C 4 H 2 (X 1 Σ g + )): Gas Phase Formation of the Heptatriynylidyne Radical ( l -C 7 H(X 2 Π))

An elucidation of the underlying formation pathways to acyclic hydrocarbons such as polyynes (C n H 2 ), cumulenes (C n H 2 ), and linear resonantly stabilized linear radicals (l-C n H) is indispensable to understand the hydrocarbon chemistry in extreme low and high temperature environments. In this study, we exploited the crossed molecular beam technique to investigate the reaction of tricarbon C 3 (X 1 Σ g + ) with diacetylene (butadiyne; HCCCCH; X 1 Σ g + ) at a collision energy of 47 ± 1 kJ mol ⁻1 . The experimental data were merged with ab initio calculations of the singlet C 7 H 2 potential energy surface (PES) revealing that the reaction is initiated via the formation of an initial van der Waals reactant complex in the entrance channel. Subsequent rearrangements lead to various carbene-type and cyclic intermediates via ring-opening, ring-closure, and hydrogen migration processes eventually forming acyclic C 7 H 2 isomers prior to their barrierless unimolecular decomposition to the most stable linear isomer, heptatriynylidyne (C 7 H, X 2 Π) in an overall endoergic reaction (+57 kJ mol ⁻1 ). The reaction exhibits strong similarities to the tricarbon – acetylene (C 3 – C 2 H 2 ). Furthermore, the significant energy threshold suggests that the tricarbon reaction with (poly)acetylenes forming resonantly stabilized linear radicals are open in high-temperature environments such as combustion flames and circumstellar envelopes of carbon stars and planetary nebulae as their descendants; however, these reactions are closed in low-temperature environments as in cold molecular clouds and hydrocarbon-rich atmospheres of planets and their moons such as in Titan.

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

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
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↗