Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHASE SHIFT”

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

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Flash flourishing of Northern Hemisphere vegetation and its drivers

Rapid surges in vegetation growth—defined by thresholds in growth rate and duration—are critical yet understudied indicators of ecosystem responses to environmental change. Here, we investigate spatiotemporal patterns of such abrupt, short-lived flash flourishing events across the northern extratropical latitudes (NEL) from 2003 to 2022. We find more frequent occurrence of flash flourishing events at high latitudes (≥45° N), where their incidence is 1.6 times higher than at mid-latitudes. Moreover, there is an increasing tendency in frequency, duration, and intensity of flash flourishing events over the past two decades, alongside consistent rises in vegetation indices across onset, post-onset, and entire phases. Model simulations attribute these multiyear increases primarily to elevated atmospheric CO 2 , while temperature and radiation predominantly control phase-specific variability, with onset traits strongly predicting subsequent phenological responses. Together, these findings identify the patterns and drivers of NEL flash flourishing and highlight their large-scale impacts on ecosystem dynamics, offering critical insights for model improvement and the assessment of ecological shifts.

Kong, Xiangxu [Nanjing University of Information S↗
DOE OSTI2024

Origin of spin-driven ferroelectricity and effect of external pressure on the complex magnetism of the 6⁢ H perovskite Ba 3 ⁢Ho⁢Ru 2⁢ O 9

The compound Ba 3 HoRu 2 O 9 magnetically orders at 50 K (T N1 ), followed by another complex magnetic ordering at 10.2 K (T N2 ). The second magnetic phase transition was characterized by the coexistence of two competing magnetic ground states associated with two different magnetic wave vectors (K 1 = 0.5 0 0 and K 2 = 0.25 0.25 0). The multiferroicity and magnetoelectric coupling were predicted below T N2 in these 4d-based materials. Here, in this work, we have discussed the origin of spin-driven ferroelectricity, which is not known yet, and the nature of magnetoelectric domains. We have investigated the compound through time-of-flight neutron diffraction, synchrotron x-ray diffraction (XRD), ac susceptibility, frequency-dependent complex dielectric spectroscopy, and dc magnetization under external pressure. We have demonstrated that the noncollinear structure involving two different magnetic ions, Ru (4d) and Ho (4f), breaks the spatial inversion symmetry via inverse Dzyaloshinskii-Moriya (DM) interaction through strong 4d-4f magnetic correlation, which shifts the oxygen atoms and results in nonzero polarization. Such an observation of inverse DM interaction from two different magnetic ions which cause ferroelectricity is rarely observed. The stronger spin-orbit coupling of 4d orbital might play a major role in creating DM interaction of noncollinear spins. We have systematically studied the spin and dipolar dynamics, which exhibit intriguing behavior with shorter coherence lengths of second magnetic phase associated with the k 2 wave vector. The results manifest the development of finite-size magnetoelectric domains instead of true long-range ordering, which justifies the experimentally obtained low value of ferroelectric polarization. The lattice parameters and volume show a sharp anomaly at T N2 obtained by analyzing the temperature-dependence XRD, which is consistent with the ferroelectric transition, predicting a noncentrosymmetric space group, $P\bar{6}2c$, for this compound. Furthermore, we have investigated the effect of external pressure on this complex magnetism. The result reveals an enhancement of ordering temperature by the application of external pressure (~1.6 K/GPa). The external pressure might favor stabilizing the magnetic ground state associated with second magnetic phase. Our study shows an unconventional mechanism of spin-driven ferroelectricity involving inverse DM interaction between Ru (4d) and Ho (4f) magnetic ions due to strong 4d-4f cross coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Elucidating the Gas-Phase Behavior of Nitazene Analog Protomers Using Structures for Lossless Ion Manipulations Ion Mobility-Orbitrap Mass Spectrometry

2-benzylbenzimidazoles, or “nitazenes”, are a class of novel synthetic opioids (NSOs) that are increasingly being detected alongside fentanyl analogs and other opioids in drug overdose cases. Nitazenes can be 20x more potent than fentanyl but are not routinely tested for during postmortem or clinical toxicology drug screens; thus, their prevalence in drug overdose cases may be under-reported. Traditional analytical workflows utilizing liquid chromatography-tandem mass spectrometry (LC-MS/MS) often require additional confirmation with authentic reference standards to identify a novel nitazene. However, additional analytical measurements with ion mobility spectrometry (IMS) may provide a path towards reference-free identification, which would greatly accelerate NSO identification rates in toxicology labs. Presented here are the first IMS and collision cross section (CCS) measurements on a set of fourteen nitazene analogs using a Structures for Lossless Ion Manipulations (SLIM)-Orbitrap MS. All nitazenes exhibited two high intensity baseline-separated IMS distributions, which fentanyls and other drug and drug-like compounds also exhibit. Incorporating water into the electrospray ionization (ESI) solution caused the intensities of the higher mobility IMS distributions to increase the intensities of the lower mobility IMS distributions to decrease. Nitazenes lacking a nitro group at the R1 position exhibited the greatest shifts in signal intensities due to water. Furthermore, IMS-MS/MS experiments showed that the higher mobility IMS distributions of all nitazenes produced fragment ions with m/z 72, 100, and other low intensity fragments while the lower mobility IMS distributions only produced fragment ions with m/z 72 and 100. The IMS, solvent, and fragmentation studies provide experimental evidence that nitazenes potentially exhibit three gas-phase protomers. In conclusion, the cyclic IMS capability of SLIM was also employed to partially resolve four sets of structurally similar nitazene isomers (e.g., protonitazene/isotonitazene, butonitazene/isobutonitazene/secbutonitazene), showcasing the potential of using high-resolution IMS separations in MS-based workflows for reference-free identification of emerging nitazenes and other NSOs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Significant acceleration of solid-state NMR simulations via three-angle powder averaging

The anisotropic frequency shifts imparted onto the NMR resonance frequency depend on the spherical angular coordinates that describe the orientations of the NMR interaction tensors with respect to the applied magnetic field direction. Experiments performed using magic-angle spinning, however, gain a dependence on a third angle: the rotor phase γ. Traditionally, a carousel average is performed to integrate over γ, which leads to a slow convergence of intensities without contributing to the underlying powder patterns. Herein, we show an order of magnitude acceleration in computation time may be obtained by including the γ-averaging into the main powder average to eliminate redundant calculation of resonance frequencies.

magic-angle spinning↗
DOE OSTI2025

Formation of Carbon–Carbon Interlinkage Bonds under High Pressure

The formation of carbon–carbon interlinkage bonds (CCIBs) via the chemical binding of interlayer carbon atoms of many sp 2 -bonded carbon precursors is an essential step for synthesizing various diamond and diamond-like materials. Although the existence of CCIBs may be reasonably assumed under high-pressure conditions, direct experimental evidence has been scarce. Micro-Raman spectroscopy is here employed to track in situ the evolution of C–C bonds in a pressure range from ambient to 54 GPa. A pressure-induced two-stage (polynomial and linear) shift of the G peak and new generation of the CCIB peak at about 1550 cm –1 are observed in multiple types of layer-structured carbon precursors, including glassy carbon, natural graphite, and carbon nanotubes. In conclusion, the experimental discovery of CCIBs holds significance in comprehending phase transitions of sp 2 -bonded carbon materials and has implications for the advancement of novel carbon structures.

Anode materials↗
DOE OSTI2025

Selective reduction in epitaxial SrFe 0.5 Co 0.5 O 2.5 and its reversibility

Oxygen-vacancy engineering in transition metal oxides enables programmable functionalities by modulating the valence states and local coordination of constituents. Here, we report the selective reduction of cobalt ions in epitaxial SrFe 0.5 Co 0.5 O 2.5 thin films under reducing gas environments, while iron ions remain unchanged. X-ray absorption spectroscopy reveals an absorption edge shift of 1.65 eV in the Co L-edge upon reduction, and multiplet simulations estimate a decrease in the average Co valence from Co 2.91+ to Co 2.00+ . This site- and element-specific reduction leads to the formation of a structurally distinct oxygen-deficient phase stabilized by oxygen vacancies at tetrahedral sites, as confirmed by density functional theory. Optical spectroscopy reveals an increase in the bandgap from 2.47 eV to 3.04 eV, accompanied by enhanced transparency. Furthermore, simultaneous in situ diffraction and transport measurements confirm fully reversible redox-driven transitions among three phases: reduced defective perovskite, brownmillerite, and oxygen-rich perovskite phases. These findings demonstrate that selective redox control in multi-cation oxides enables the realization of chemically and functionally distinct oxygen-deficient phases.

Fuel cells↗
DOE OSTI2026

Preliminary Assessment for the Electric Load Shifting Potential of Integrating Thermal Energy Storage with Heat Pumps in Residential Buildings in Texas of United States

The widespread adoption of electric-driven heat pumps for heating and cooling is expected to significantly increase electric demand. Cooling electric demand will result in an increase in peak hours electric loads, placing additional strain on the grid during peak hours. This challenge, combined with the current rapidly increasing demand from data centers, exacerbates the electric demand duck curve problem. Integrating thermal energy storage (TES) with heat pump can shift electric use for heating and cooling from peak to off-peak hours of the electric grid, which can help flatten the daily electric demand profile of a building. This paper presents a novel design that integrates heat pump with TES (HP-TES), which uses phase change materials (PCM). Heat pump charges TES by melting or freezing PCM during off-peak hours when there is no thermal demand from the building. TES is then discharged (i.e., by freezing or melting PCM) during peak hours to provide a more favorable heat source or heat sink for the heat pump to meet the thermal demand of the building with lower electricity use than conventional air-source heat pumps. Computer simulations were developed to predict the performance of HP-TES applied to a typical single-family house in the US. The building-level simulation results were scaled up to preliminarily assess the aggregated impacts of deploying HP-TES across all single-family houses in Texas of the United States, including reduction of peak demand of the electric grid.

Anees, Fady [ORNL]↗
DOE OSTI2026

Half-ice, half-fire-driven ultranarrow phase crossover in one-dimensional decorated 𝑞-state Potts ferrimagnets: An AI-co-led exploration

OpenAI’s reasoning model o3-mini-high was used to carry out an exact analytic study of one-dimensional ferrimagnetic site- and bond-decorated 𝑞-state Potts models. We demonstrate that the finite-temperature ultranarrow phase crossover (UNPC), driven by a hidden “half-ice, half-fire” state recently discovered in the 𝑞=2 case (Ising model), persists for 𝑞>2. Moreover, we identify unique features for 𝑞>2, including the dome structure in the field-temperature phase diagram, and for large 𝑞 a secondary high-temperature UNPC to the fully disordered paramagnetic state. As the UNPC quickly approaches a genuine transition by enhancing 𝐽, the interaction between the backbone spins, two distinct behaviors emerge: In the site-decorated Potts model, 𝑇 0 is independent of 𝐽 and thus remains unchanged (Type-I UNPC), and in the bond-decorated Potts model with 𝑞>2, 𝑇 0 depends on 𝐽 and quickly shifts toward a finite temperature as 𝐽 increases (Type-II UNPC). These results establish a versatile framework for engineering controlled fast state-flipping switches in low-dimensional systems. Our nine-dan artificial intelligence (AI)-contribution framework assigns AI the meritorious status of AI-co-led discovery in this work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A machine learning approach to quantify degradation of nuclear fuels and the effects of fission products

Nuclear fuel performance is critically dependent on understanding the evolution of fuel properties under operational conditions, a complex challenge driven by chemical changes and substantial radiation damage during fission. Traditionally, property evolution has been determined via empirical data collected following irradiation. However, these empirical correlations are limited in their applicability beyond the specific conditions in which they were obtained. This study explores a novel approach to address this challenge by applying materials informatics to develop a machine learning random forest (ML-RF) model that captures the effects of fission products on fuel compounds. The model predicts formation enthalpy (ΔH f ) by leveraging extensive quantum materials property data and correlating it with material descriptors such as composition, atomic and site features, and crystal lattice properties. This ML-RF model enables rapid interpolation across the compositional and structural spaces covered by the training data, thus supporting high-throughput screening and energetic ranking of candidate phases. The model demonstrates the ability to predict ΔH f with a mean absolute error (MAE) of approximately 0.1 to 0.2 eV/atom across a wide range of compounds, including key nuclear fuel systems (U-O, U-N, U-C, U-Si, and U-Mo). For example, it was used to assess shifts in stoichiometry for UO 2 (O/M) and UN (N/M) fuels, revealing their distinct tendencies in chemical potential variation and enabling preliminary convex hull analyses. Furthermore, the model provides insights into how individual fission products affect fuel properties. Results indicate that larger fission products (e.g., Nd, Pu, Ce) have a more pronounced impact on UO 2 , while lighter ones (e.g., Zr) strongly influence UN. Here, the model developed in this work can be used to support the Accelerated Fuel Qualification approach by facilitating preliminary evaluations prior to extensive materials modeling and experimentation. To this end, the trained model has been made available to the fuel community to support ongoing fuel development efforts.

Accelerated fuel qualification↗
DOE OSTI2026

Ultrafast Dynamics of Excitonic Complexes in Electrolyte-Gated Monolayer MoS 2

Understanding and controlling ultrafast excitons and trion dynamics in monolayer transition metal dichalcogenides (TMDs) is critical for optoelectronic and photonic applications. These dynamics depend strongly on carrier density, but most studies use fluence-dependence to modulate electron-hole populations rather than direct electrostatic gating that is more relevant to optoelectronic devices. We utilize electrolyte gating to tune the ground-state carrier density of monolayer MoS 2 and probe excitonic dynamics using transient absorption spectroscopy, thereby modifying absorption through bandgap renormalization, screening, phase-space filling, and exciton-trion crossover. Spectral deconvolution reveals distinct but coupled exciton and trion dynamics. Excitons form independently of voltage but relax faster with increasing carrier density, consistent with Auger-like processes, while trions decay more rapidly through multiple channels. At higher voltages, trion relaxation shifts from a sub-picosecond cooling to slower trapping-assisted processes. Furthermore, our results provide missing insights into the interplay between excitons and trions in monolayer MoS 2 , relevant for TMD-based optoelectronics.

chalcogenides↗
DOE OSTI2026

Unraveling In-Situ Formation of Surface Nickel Nitride Structures in Plasma-Assisted Catalytic Ammonia Synthesis

We report the in situ formation of Ni nitride for plasma-assisted ammonia synthesis. Both the surface nitrogen concentration and the ammonia formation rate exhibit dependence on the N 2 :H 2 feed ratio. The maximum surface nitrogen concentration occurs at a N 2 :H 2 ratio of 4:1, and the maximum catalytic activity occurs at 2:1. In contrast, the formation of gas phase radicals is less sensitive to feed composition, indicating that Ni nitride is more kinetically relevant to ammonia production than gas-phase radicals. The plasma-induced formation of Ni nitride is therefore proposed to be a critical contributor to the synergistic effects in plasma-assisted catalytic ammonia synthesis. Additionally, Ni nitride alters the surface reaction mechanism of plasma-assisted ammonia synthesis, with the rate-determining-step (RDS) shifting to surface-bound NH 3 formation rather than N 2 activation at temperatures below 373 K. These findings provide mechanistic insight that opens opportunities for optimizing the performance of plasma-assisted catalytic ammonia synthesis

Ammonia↗
DOE OSTI2024

Nitrogen availability and summer drought, but not N:P imbalance, drive carbon use efficiency of a Mediterranean tree-grass ecosystem

All ecosystems contain both sources and sinks for atmospheric carbon (C). A change in their balance of net and gross ecosystem carbon uptake, ecosystem-scale carbon use efficiency (CUE ECO ), is a change in their ability to buffer climate change. However, anthropogenic nitrogen (N) deposition is increasing N availability, potentially shifting terrestrial ecosystem stoichiometry towards phosphorus (P) limitation. Depending on how gross primary production (GPP, plants alone) and ecosystem respiration (R ECO , plants and heterotrophs) are limited by N, P or associated changes in other biogeochemical cycles, CUE ECO may change. Seasonally, CUE ECO also varies as the multiple processes that control GPP and respiration and their limitations shift in time. We worked in a Mediterranean tree-grass ecosystem (locally called ‘dehesa’) characterized by mild, wet winters and summer droughts. We examined CUE ECO from eddy covariance fluxes over 6 years under control, +N and + NP fertilized treatments on three timescales: annual, seasonal (determined by vegetation phenological phases) and 14-day aggregations. Finer aggregation allowed consideration of responses to specific patterns in vegetation activity and meteorological conditions. We predicted that CUE ECO should be increased by wetter conditions, and successively by N and NP fertilization. Milder and wetter years with proportionally longer growing seasons increased CUE ECO , as did N fertilization, regardless of whether P was added. Using a generalized additive model, whole ecosystem phenological status and water deficit indicators, which both varied with treatment, were the main determinants of 14-day differences in CUE ECO . The direction of water effects depended on the timescale considered and occurred alongside treatment-dependent water depletion. Overall, future regional trends of longer dry summers may push these systems towards lower CUE ECO .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Significantly Promoting the Thermal Conductivity and Machinability of Negative Thermal Expansion Alloy via In Situ Precipitation of Copper Networks

Rapid advancements in electronic devices yield an urgent demand for high-performance electronic packaging materials with high thermal conductivity, low thermal expansion, and great mechanical properties. However, it is a great challenge for current design philosophies to fulfill all the requirements simultaneously. Here, an effective strategy is proposed for significantly promoting the thermal conductivity and machinability of negative thermal expansion alloy (Zr,Nb)Fe 2 through eutectic precipitation of copper networks. The eutectic dual-phase alloy exhibits an isotropic chips-matched thermal expansion coefficient and a thermal conductivity enhancement exceeding 200% compared with (Zr,Nb)Fe 2 , along with an ultimate compressive strength of 550 MPa. The addition of copper reorganizes the composition of (Zr,Nb)Fe 2 , which smooths the magnetic transition and shifts it toward higher temperature, resulting in linear low thermal expansion in a wide temperature range. The highly fine eutectic copper lamellae construct high thermal conductivity networks within (Zr,Nb)Fe 2 , serving as highways for heat transfer electrons and phonons. The in situ forming of eutectic copper lamellae also facilitates the mechanical properties by enhancing interfacial bonding and bearing additional stress after yielding of (Zr,Nb)Fe 2 . This work provides a novel strategy for promoting thermal conductivity and mechanical properties of negative thermal expansion alloys via eutectic precipitation of copper networks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying and visualizing the microscopic degrees of freedom of grain boundaries in the Wigner–Seitz cell of the displacement-shift-complete lattice

Here we introduce a grain boundary (GB) translation vector, t ws , to describe and quantify the domain of the microscopic degrees of freedom of GBs. It has long been recognized that for fixed macroscopic degrees of freedom of a GB there exists a large multiplicity of states characterized by different relative grain translations. More recently another degree of freedom, [n], the number of GB atoms, has emerged and is now recognized as an equally important component of GB structural multiplicity. In this work, we show that all GB microstates can be uniquely characterized by their value of t ws , which is located within the Wigner–Seitz (WS) cell of the Displacement-Shift-Complete lattice (DSCL) of the GB. The GB translation vector captures information about both the translation state and the number of GB atoms. We show that the density of GB microstates inside the WS cell of the DSCL is not uniform and can form clusters that correspond to different GB phases. Further, the vectors connecting the centers of the clusters correspond to the Burgers vectors of GB phase junctions, which can be predicted without building the junctions. Using t ws , we quantify GB excess shear and argue that it is defined up to a DSCL vector, which has implications for thermodynamic equilibrium conditions. Additionally, this work generalizes the definition of the number of GB atoms [n] to asymmetric boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

I/O in Machine Learning Applications on HPC Systems: A 360-degree Survey

Growing interest in Artificial Intelligence (AI) has resulted in a surge in demand for faster methods of Machine Learning (ML) model training and inference. This demand for speed has prompted the use of high performance computing (HPC) systems that excel in managing distributed workloads. Because data is the main fuel for AI applications, the performance of the storage and I/O subsystem of HPC systems is critical. In the past, HPC applications accessed large portions of data written by simulations or experiments or ingested data for visualizations or analysis tasks. ML workloads perform small reads spread across a large number of random files. This shift of I/O access patterns poses several challenges to modern parallel storage systems. In this paper, we survey I/O in ML applications on HPC systems, and target literature within a 6-year time window from 2019 to 2024. We define the scope of the survey, provide an overview of the common phases of ML, review available profilers and benchmarks, examine the I/O patterns encountered during offline data preparation, training, and inference, and explore I/O optimizations utilized in modern ML frameworks and proposed in recent literature. Lastly, we seek to expose research gaps that could spawn further R&D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Imaging Bragg Edge Analysis TooLs for Engineering Structures (iBeatles)

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) provides pulsed neutrons with energies varying from epithermal to cold. In preparation for VENUS, the neutron imaging beamline to be located at beam port 10, we have performed a series of experiments focused on wavelength-dependent radiography and computed tomography for a broad range of applications, from materials science to biological tissues.One of the time-of-flight (TOF) techniques that is of interest to the scientific community is the 2-dimensional mapping of phases and average crystalline plane orientation in samples both ex-situ and during applied stresses such as tensile loading and heating. This technique is known as Bragg edgeimaging and relies on the identification of changes of transmission values, fitting of the edge to measure its displacement, and thus identify the shift in lattice parameter due to stresses. One of the challenges of TOF imaging measurements is the amount of data and the inability to observe Bragg edge shifts in real time during an experiment. Thus, we have been focusing on creating a Python-based interface that allows fast data processing and instantaneous mapping and fitting of the Bragg edges, and their evolution through time. Python libraries and Jupyter notebooks have been implemented to facilitate decision making during an experiment. The advantage of the notebooks is the possibility to guide an experiment as they can quickly process and display Bragg edge data. These notebooks can be used independently, or can be combined in a Python Graphical User Interface (GUI) tool called iBeatles. This interface permits visualization and fitting of the Bragg edges, and ultimately back-projects the fitting results onto the radiographs to display a strain map. Assuming data collection has sufficient statistics, the strain mapping analysis can be performed on a pixel-by-pixel basis. This development is a step forward toward a better user experience at the future VENUS beamline in terms of live feedback and productivity. Analysis that used to take days of switching between different applications can now be done in minutes within the

Bilheux, JeanChristophe [Oak Ridge National Labora↗
DOE OSTI2026

Directed Gas-Phase Formation of The Propargyl Family of Resonance-Stabilized Radicals in The Reactions of Ground-State Carbon Atoms (C; 3 P j ) with Butene Isomers (C 4 H 8 ): Dimethylpropargyl and Ethylpropargyl

The propargyl radical (C 3 H 3 ) is the simplest resonance-stabilized free radical (RSFR), but how does stepwise methyl substitution in the alkene reactant affect its dynamics of their formation? We report a crossed molecular beam study of the reactions of atomic carbon (C, 3 P j ) with four butene isomers (C 4 H 8 ) under single collision conditions at a collision energy of 28 ± 2 kJ mol −1 . Barrierless addition of atomic carbon to the alkene C=C bond triggers ring opening to triplet substituted allenes—a de facto insertion mechanism—followed by unimolecular decomposition via atomic hydrogen (H), methyl (CH 3 ) or ethyl (C 2 H 5 ) loss, yielding a family of propargyl‑type RSFRs. RRKM calculations reveal that the branching ratios are highly sensitive to the alkene structure. While the methyl loss channel, affording 1‑methylpropargyl, dominates for 2 butenes (80–90%), the predicted hydrogen atom loss channel (≈10%) leading to 1,3 dimethylpropargyl is identified in the experiment by comparison with theoretical energetics. For isobutene, a near‑equal competition is seen, with the reaction producing 3‑methylpropargyl (≈50%) and 1,1‑dimethylpropargyl (≈40%), along with 2‑vinylallyl (≈5%), whose formation is supported by the experimental data. Most notably, the reaction with 1 butene uniquely favors an enthalpically driven hydrogen shift, eventually producing 1 vinylallyl (≈38%), which is assigned based on the excellent agreement between the measured and calculated reaction exothermicity. Rapid entropically favored fragmentation channels yield ≈40% of propargyl type species (propargyl, 1- and 3-ethylpropargyls), slightly outcompeting the allyl type product. Furthermore, these results establish a systematic progression from C 2 H 4 to C 4 H 8 , where the increasing alkyl substitution unlocks new fragmentation channels, providing a versatile gas phase route to alkylated RSFRs—key intermediates in the growth of methylated and ethylated PAHs and aliphatic chains in combustion and cold interstellar environments (molecular clouds).

Alkyls↗