SEARCH · Search NASA

Results for “Cold Atoms”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Directed Gas-Phase Formation of the 1-Cyanovinyl Radical (H 2 CCCN, X 2 A′) in the Interstellar Medium

The formation pathways to nitrogen-containing molecules and radicals are crucial to the understanding of the carbon–nitrogen chemistry in interstellar and atmospheric environments. While over 65 nitrogen-containing neutral species have been observed in deep space to date, their formation mechanisms─in particular, those of radical species─remain largely speculative. The crossed molecular beam technique in conjunction with electronic structure and statistical calculations was utilized to offer a detailed overview of the fundamental pathways in the gas-phase bimolecular reaction of ground-state atomic carbon (C, 3 P) with acetonitrile-d 3 (CD 3 CN, X 1 A 1 ) under single-collision conditions leading to the formation of the 1-cyanovinyl radical (D 2 CCCN, X 2 A′) coupled with deuterium atom loss. Here, the indirect reaction was initiated by barrierless carbon-atom addition, with the most probable route involving carbon addition across the carbon–nitrogen nitrile triple bond of acetonitrile, forming a three-membered ring intermediate followed by ring-opening and unimolecular decomposition via atomic deuterium loss from the C3 carbon atom. The reaction was overall exoergic, and intermediates and transition states lie lower in energy than the separated reactants, unlocking the reaction of carbon with acetonitrile in low-temperature environments such as cold molecular clouds, e.g., Taurus Molecular Cloud (TMC-1), and planetary atmospheres, e.g., Saturn’s moon Titan. In these environments, the 1-cyanovinyl radical may act as a building block for cyano-substituted polycyclic aromatic hydrocarbons and N-heterocycles, thus furthering our understanding of the complex carbon–nitrogen chemistry in deep space.

Chemical reactions↗

DOE OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Design of a shipping fixture for a compact cryomodule hermetic assembly

Two conduction-cooled 915 MHz superconducting radio frequency hermetic assemblies must be safely transported from the Jefferson Lab in Newport News, VA to General Atomics in San Diego, CA for perfor-mance testing in a custom horizontal test cryostat. One hermetic assembly consists of a 2-cell 915 MHz cavity, a coaxial fundamental power coupler, and the warm-to-cold transition beam tubes. The second hermetic assembly consists of a 2-cell 915 MHz cavity only. The assemblies will be transported on a flatbed air-ride trailer over the approximate 4000 km distance. Design requirements included adequate attenuation of 4g vertical axis, 5g beamline axis, and 1.5g lateral axis shock events. The isolation system was designed using helical wire-rope isolators with modal and transient finite element analysis performed in Ansys. Results show shock attenuation of a 10 ms half-sine pulse input to < 1g in the vertical axis, < 1.5g in the beam-line axis, and < 0.5g in the lateral axis for both assem-blies at the specified design loads and all structural stresses are kept below the material yield limits. Addi-tionally, the natural frequencies of both isolation sys-tems adequately attenuate the fundamental modes of the critical structures.

Accelerator Physics↗

DOE OSTI2024

A near X-point charge exchange neutral spectroscopy (CENS) system for DIII-D

A 16 channel spectroscopy system has been installed on DIII-D to provide information about the energy distribution of the atomic neutrals using the Doppler shift and broadening of passive Balmer-α emission. Here, the lines of sight are dominantly in the toroidal direction, with tangencies near the lowest point on closed magnetic flux surfaces moving from the lower divertor upward into the confined plasma. This allows the system to provide details of the neutrals as they undergo various atomic physics processes while traveling upward from the lower divertor. The spectrally resolved measurements provide several advantages that complement typical optical filter based measurements of hydrogenic spectral lines. These include more direct measurements of the neutrals that fuel deeper in the plasma and capturing neutral velocity distribution information via the Doppler broadened and shifted line. In addition, the spectral separation of higher energy “thermal” neutrals and bright cold emission from the scrape-off layer allows for a more straightforward calculation of the underlying neutral densities based on the emission because of the dependence of the thermal emission on confined plasma properties that are approximately flux functions and well measured.

Atomic and molecular physics↗

DOE OSTI2025

Thermodynamically leveraged solventless aerobic deconstruction of polyethylene-terephthalate plastics over a single-site molybdenum-dioxo catalyst

Here, we describe the solventless catalytic deconstruction of polyethylene-terephthalate (PET) under an aerobic atmosphere, mediated by an earth-abundant, low-cost activated carbon (AC)-supported single-site molybdenum-dioxo catalyst (AC/MoO 2 ). Catalytic amounts of AC/MoO 2 selectively convert waste PET into its monomer, terephthalic acid (TPA), within 4 h at 265 °C with yields as high as 94% under 1 atm air. Pure crystalline TPA product sublimes from the reaction hot zone, crystallizing on the reactor cold zone, thus avoiding the need for separation and purification steps. This process does not employ any hazardous/toxic reducing agents or solvents, and the catalyst can be recycled multiple times without loss of activity, rendering this process highly atom-efficient. According to computational and experimental mechanistic studies, the AC/MoO 2 catalyst mediates a thermoneutral metal-catalyzed β-scission step, followed by an exothermic step that converts the vinyl benzoate intermediate to TPA and acetaldehyde using trace amounts of moisture in the air. The formation of gaseous acetaldehyde makes the isolation of TPA from the reaction mixture facile and industrially favorable, especially since solvents are unnecessary. The present methodology is also extended to the deconstruction of other frequently used polyester plastics, polybutylene terephthalate (PBT), polyethylene naphthalate (PEN), and polyethylene furanoate (PEF), and operates equally well with post-consumer waste products. Notably, this process is also compatible with plastic mixtures of polyesters with polyolefins, polyamides, and polycarbonates, leading to the selective conversion of each polyester to the corresponding monomer, leaving the residual polymer unchanged and polyester-free.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Effects of HCP/BCC element ratios on the room-temperature tensile properties of Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys

Equiatomic and non-equiatomic Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys (RHEAs) were arc melted, homogenized, cold rolled, and recrystallized to produce single-phase, body-centered cubic (BCC), microstructures with weak texture and equiaxed grains 76–199 μm in size. Here, the non-equiatomic alloys had either a 60:40 or 80:20 atomic ratio of hexagonal close-packed (HCP) elements (Ti + Zr + Hf) to BCC elements (Nb + Ta). Alloy compositions were measured after thermomechanical processing to determine the concentrations of the major (substitutional) and minor (interstitial) elements. We investigated how elastic constants and uniaxial tensile properties were affected by changes in the relative concentrations of the constituent elements at fixed HCP:BCC ratios. Yield strengths ranged from 801 to 922 MPa and ultimate tensile strengths from 815 to 933 MPa. Good agreement is obtained between the experimental yield strengths and those predicted by a strength theory based on edge dislocations indicating that the observed compositional effects are due to their effects on shear modulus and volume misfit. Fracture occurred by dimpled rupture with fracture strains of 19.4%–25.7%, but uniform strains were an order of magnitude lower at 1.1%–3.2%, calling into question the useable ductility (prior to necking) of RHEAs considered to be ductile based on their fracture strain. Contrary to predictions in the literature that HCP elements promote ductility, our present results show that increasing the HCP:BCC ratio decreases both the total strain and the uniform strain. Similar trends were not evident in the yield or ultimate strengths; rather, strengths were affected mainly by shear modulus and volume misfit.

BCC high-entropy alloys↗

DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

The atomic gas sequence and mass–metallicity relation from dwarfs to massive galaxies

Galaxy scaling relations provide insights into the processes that drive galaxy evolution. The extension of these scaling relations into the dwarf galaxy regime is of particular interest. This is because dwarf galaxies represent a crucial stage in galaxy evolution, and understanding them could also shed light on their role in reionizing the early Universe. There is currently no consensus on the processes that dominate the evolution of dwarfs. In this work, we constrain the atomic gas sequence (stellar mass versus atomic gas fraction) and mass–metallicity relation (stellar mass versus gas-phase metallicity) from dwarf (⁠10 6.5 M ⊙ ⁠) to massive (⁠10 11.5 M ⊙ ) galaxies in the local Universe. The combined optical and 21-cm spectroscopic observations of the Dark Energy Spectroscopic Instrument and Arecibo Legacy Fast ALFA surveys allow us to constrain both scaling relations simultaneously. We find a slope change of the atomic gas sequence at a stellar mass of ~10 9 M ⊙ ⁠. We also find that the shape and scatter of the atomic gas sequence and mass–metallicity relation are strongly linked for both dwarfs and more massive galaxies. Consequently, the low-mass slope change of the atomic gas sequence is imprinted onto the mass–metallicity relation of dwarf galaxies. The mass scale of the measured slope change is consistent with a predicted escape velocity threshold below which low-mass galaxies experience significant supernova-driven gas loss, as well as with a reduction in cold gas accretion onto more massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Energy shift of Fe-K fluorescence lines due to low ionization demonstrated with XRISM in Centaurus X-3

The Fe K$\alpha$ fluorescence line at 6.4 keV is a powerful probe of cold matter surrounding X-ray sources and has been widely used in various astrophysical contexts. The X-ray microcalorimeter spectrometer onboard XRISM can measure line shifts with unprecedented precision of ${\sim }$0.2 eV, equivalent to a line-of-sight velocity of ${\sim }$10 km s$^{-1}$. At this level of accuracy, however, several factors that influence the line energy must be carefully considered prior to astrophysical interpretation. One such important factor is the ionization degree, Fe$^{q+}$. The K$\alpha$ line shifts redward by ${\sim }$4 eV as q increases from 0 (neutral) to 8 (Ar-like). Additionally, the accompanying Fe K$\beta$ line at 7.06 keV shifts blueward by ${\sim }$30 eV from $q=0$ to 8. We demonstrate that this effect is actually observable in the XRISM data of the high-mass X-ray binary Centaurus X-3 (Cen X-3). We advocate that the differential energy shift between the K$\alpha$ and K$\beta$ line provides a robust estimate of q by decoupling from other effects that shift the two lines in the same direction. We derived $q \sim 5$ (Sc-like) for the fluorescing matter by comparing the observation with atomic structure calculations of our own and in the literature. By accounting for the derived charge state and the corresponding shift in the rest-frame line energy, we made corrections for this effect and reached a consistent residual shift among the K$\alpha$, K$\beta$, and the optical measurement attributable to the systemic velocity of the system. Consequently, we obtained a new constraint on the location of the cold matter. This ionization effect needs to be assessed in all use cases of the Fe K$\alpha$ line shift beyond Cen X-3, and the proposed metric is generally applicable to all of them.

Astronomy and AstroPhysics↗

DOE OSTI2025

Microstructural features underpinning the mechanical behavior of powder metallurgy Cr-based alloys

Refractory materials such as Cr-based alloys offer the potential of enhanced elevated temperature performance but have been limited by their poor formability. However, powder metallurgy has been shown to be a viable pathway to fabricate these alloys. Nanophase separation sintering (NPSS) in particular has been used in the literature to accelerate the densification of powder metallurgy Cr- and W-based alloys. Here, we explore microstructure evolution during NPSS in a binary Cr 85 Ni 15 alloy consolidated via (i) cold pressing & pressureless sintering and (ii) hot isostatic pressing followed by hot extrusion & hot upsetting, and the role of these different processing routes on resulting material properties. The alloy lacked room temperature tensile ductility regardless of consolidation process, with multi-length scale characterization, including scanning electron microscopy, transmission electron microscopy, atom probe tomography, X-ray diffraction, and uniaxial tensile testing, revealing that brittleness was due to intrinsically poor Cr grain boundary cohesion. Tensile testing conducted at 760 °C showed marked strength reduction for the extruded & upset (94 %) condition compared to the pressed & sintered (18 %). The formation of orthorhombic CrNi 2 intermetallics functioned as precipitate strengtheners and prevented elevated temperature softening in the pressed and sintered condition. The findings offer foundational insights into aiding the future development of Cr-based alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Crystal Growth of Cesium Hafnium Chloride (Cs2HfCl6) Scintillating Crystal

Crystal growth is the process where crystalline substance forms and increases in size through the addition of atoms, ions, or molecules. In this research we used the Bridgman technique to grow Cesium Hafnium Chloride (C2HfCl6). This technique is ideal for growing large volume crystal. Using a vertical furnace with two heating zones for heating and cooling. the material solidifies into a crystal as it slowly transitions from hot to cold zone.

Tribble, Jalen↗

DOE OSTI2024

Cold Seeded Epitaxy and Flexomagnetism in GdAuGe Membranes Exfoliated From Graphene/Ge(111)

Remote and van der Waals epitaxy are promising approaches for synthesizing single crystalline membranes for flexible electronics and discovery of new properties via extreme strain; however, a fundamental challenge is that most materials do not wet the graphene surface. We develop a cold seed approach for synthesizing smooth intermetallic films on graphene that can be exfoliated to form few nanometer thick single crystalline membranes. Our seeded GdAuGe films have narrow X-ray rocking curve widths of 9–24 arc seconds, which is 2 orders of magnitude lower than their counterparts grown by typical high temperature methods, and have atomically sharp interfaces observed by transmission electron microscopy. Upon exfoliation and rippling, strain gradients in GdAuGe membranes induce an antiferromagnetic to ferri/ferromagnetic transition. Finally, our smooth, ultrathin membranes provide a clean platform for discovering new flexomagnetic effects in quantum materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Tailoring Cu-Zr gradient nanoglass structures: Influence of nanoparticle size and cooling rates on glass-glass interfaces

The study of gradient nanoglasses (GNGs) has gained attention due to their unique mechanical properties and potential applications in advanced materials. This study employs molecular dynamics simulations to synthesize a GNG using Cu-Zr metallic glass nanoparticles (NPs) sized from 3 to 15 nm. The NPs were produced by melting and quenching metallic clusters at a relatively slow quench rate of 10 9 K/s. The synthesis of GNG is elucidated along with the characterization of its heterogeneous metallic glass nanostructure. A seamless GNG structure is formed through cold compression of Cu 64 Zr 36 amorphous NPs of varying sizes. The influence of NP size on the GNG structure is investigated, utilizing deeply relaxed NPs, which exhibit a characteristic Cu segregation pattern on their surfaces. The results highlight an increase in structural heterogeneity due to heterogeneous mass transport and the development of local composition and density variations caused by Cu segregation at glass-glass interfaces (GGIs). A reduction in NP size is correlated with decreased Cu atomic displacements and local density at GGIs, suggesting that larger NPs may produce stronger GGIs. This research presents a novel methodology for synthesizing heterogeneous metallic glasses, demonstrating the capacity to control and customize nanostructure heterogeneity through the manipulation of NP sizes and cooling rates. Furthermore, these findings enhance our understanding of structural evolution during nanoglass synthesis and lay the foundation for further exploration in nanomaterial synthesis and characterization.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Investigations of atomic and molecular processes of NBI-heated discharges in the MAST Upgrade Super-X divertor with implications for reactors

This experimental study presents an in-depth investigation of the performance of the MAST-U Super-X divertor during NBI-heated operation (up to 2.5 MW) focussing on volumetric ion sources and sinks as well as power losses during detachment. The particle balance and power loss analysis revealed the crucial role of Molecular Activated Recombination and Dissociation (MAR and MAD) ion sinks in divertor particle and power balance, which remain pronounced in the change from ohmic to higher power (NBI heated) L-mode conditions. The importance of MAR and MAD remains with double the absorbed NBI heating. MAD results in significant power dissipation (up to ${\sim} 20\%$ of $P_\textrm{SOL}$), mostly in the cold ($T_e \lt 5$ eV) detached region. Theoretical and experimental evidence is found for the potential contribution of $D^-$ to MAR and MAD, which warrants further study. These results suggest that MAR and MAD can be relevant in higher power conditions than the ohmic conditions studied previously. Post-processing reactor-scale simulations suggests that MAR and MAD can play a significant role in divertor physics and synthetic diagnostic signals of reactor-scale devices, which are currently underestimated in exhaust simulations. This raises implications for the accuracy of reactor-scale divertor simulations of particularly tightly baffled (alternative) divertor configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Tunable pairing with local spin-dependent Rydberg molecule potentials in an atomic Fermi superfluid

We explore the energy spectrum and eigenstates of two-component atomic Fermi superfluids with tunable pairing interactions in the presence of spin-dependent ultralong-range Rydberg molecule (ULRM) potentials, within the Bogoliubov–de Gennes formalism. The attractive ULRM potentials lead to local-density accumulation, while their difference results in a local polarization potential and induces the in-gap Yu-Shiba-Rusinov (YSR) states whose energies lie below the bulk energy gap. A transition from equal population to population imbalance occurs as the pairing strength falls below a critical value, accompanied by the emergence of local Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)–like states characterized by out-of-phase wave functions and lower energies compared to the YSR states. The negative contribution emanating from the FFLO-like states also causes a sign change in the gap function within the ULRM potentials. Depending on the Rydberg excitation, the transition towards population imbalance can be on either the BCS or the Bose-Einstein condensation side of the Fermi superfluid. Additionally, spin-polarized bound states arise along with oscillatory “clumpy states” to compensate for the local-density difference. Here, we discuss possible experimental realizations of the composite Rydberg-atom–Fermi-superfluid system.

Cold and ultracold molecules↗

DOE OSTI2024

MHONGOOSE: A MeerKAT nearby galaxy H I survey

The MHONGOOSE (MeerKAT H IObservations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (H I) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low H Icolumn densities. The H Icolumn density sensitivity (3σover 16 km s −1 ) ranges from ∼5 × 10 17 cm −2 at 90″ resolution to ∼4 × 10 19 cm −2 at the highest resolution of 7″. The H Imass sensitivity (3σover 50 km s −1 ) is ∼5.5 × 10 5 M ⊙ at a distance of 10 Mpc (the median distance of the sample galaxies). The velocity resolution of the data is 1.4 km s −1 . One of the main science goals of the survey is the detection of cold accreting gas in the outskirts of the sample galaxies. The sample was selected to cover a range in H Imasses from 10 7 M ⊙ to almost 10 11 M ⊙ in order to optimally sample possible accretion scenarios and environments. The distance to the sample galaxies ranges from 3 to 23 Mpc. In this paper, we present the sample selection, survey design, and observation and reduction procedures. We compared the integrated H Ifluxes based on the MeerKAT data with those derived from single-dish measurement and find good agreement, indicating that our MeerKAT observations are recovering all flux. We present H Imoment maps of the entire sample based on the first ten percent of the survey data, and find that a comparison of the zeroth- and second-moment values shows a clear separation in the physical properties of the H Ibetween areas with star formation and areas without related to the formation of a cold neutral medium. Finally, we give an overview of the H I-detected companion and satellite galaxies in the 30 fields, five of which have not previously been cataloged. We find a clear relation between the number of companion galaxies and the mass of the main target galaxy.

Astronomy & Astrophysics↗

DOE OSTI2025

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗