SEARCH · Engineering Papers

Results for “high mass loading”

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.

773 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2023

FAST Discovery of a Fast Neutral Hydrogen Outflow

Abstract In this letter, we report the discovery of a fast neutral hydrogen outflow in SDSS J145239.38+062738.0, a merging radio galaxy containing an optical type I active galactic nucleus (AGN). This discovery was made through observations conducted by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) using redshifted 21 cm absorption. The outflow exhibits a blueshifted velocity likely up to ∼−1000 km s −1 with respect to the systemic velocity of the host galaxy with an absorption strength of ∼−0.6 mJy beam −1 corresponding to an optical depth of 0.002 atv= −500 km s −1 . The mass outflow rate ranges between 2.8 × 10 −2 and 3.6M ⊙ yr −1 , implying an energy outflow rate ranging between 4.2 × 10 39 and 9.7 × 10 40 erg s −1 , assuming 100 K s< 1000 K. Plausible drivers of the outflow include the starbursts, AGN radiation, and radio jet, the last of which is considered the most likely culprit according to the kinematics. By analyzing the properties of the outflow, AGN, and jet, we find that if the Hioutflow is driven by the AGN radiation, the AGN radiation does not seem powerful enough to provide negative feedback, whereas the radio jet shows the potential to provide negative feedback. Our observations contribute another example of a fast outflow detected in neutral hydrogen and demonstrate the capability of FAST in detecting such outflows.

Astronomy & Astrophysics↗

DOE OSTI2025

Detection of visible-wavelength aurora on Mars

Mars hosts various auroral processes despite the planet’s tenuous atmosphere and lack of a global magnetic field. To date, all aurora observations have been at ultraviolet wavelengths from orbit. We describe the discovery of green visible-wavelength aurora, originating from the atomic oxygen line at 557.7 nanometers, detected with the SuperCam and Mastcam-Z instruments on the Mars 2020 Perseverance rover. Near–real-time simulations of a Mars-directed coronal mass ejection (CME) provided sufficient lead-time to schedule an observation with the rover. The emission was observed 3 days after the CME eruption, suggesting that the aurora was induced by particles accelerated by the moving shock front. To our knowledge, detection of aurora from a planetary surface other than Earth has never been reported, nor has visible aurora been observed at Mars. This detection demonstrates that auroral forecasting at Mars is possible, and that during events with higher particle precipitation, or under less dusty atmospheric conditions, aurorae will be visible to future astronauts.

Science & Technology - Other Topics↗

DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗

DOE OSTI2026

Interfacial Adhesion Mechanism of Anionic Polyelectrolyte Brushes Induced by Oppositely Charged Macromolecular Counterions

Abstract Polyelectrolyte brushes are widely used as model systems for investigating electrostatic interactions at soft interfaces and to achieve exceptional lubrication properties. Previous studies have primarily focused on their behavior in the presence of multivalent counterions, which induce brush collapse, ionic crosslinking, and pronounced changes in interfacial structures. However, interactions between polyelectrolyte brushes and oppositely charged macromolecular counterions, such as polycations, remain poorly understood. Here, we prepare well-defined polystyrene sulfonate (PSS) brushes fabricated via surface-initiated grafting and employ surface forces apparatus measurements to investigate their interactions with oppositely charged polycations. In contrast to multivalent ions, polycations do not induce noticeable brush collapse, but instead generate significant adhesion between symmetric PSS brush layers. This adhesion increases with both contact time and applied load, eventually reaching a steady plateau, indicating that the interaction is governed not simply by electrostatic screening but by the gradual formation of polycation-mediated bridging under confinement. Furthermore, the introduction of monovalent Na+ ions disrupts the adhesive interaction even at very low concentrations, suggesting that the bridging function of the adsorbed polycation is highly sensitive to competitive ionic screening. In comparison, measurements with trivalent counterions reveal the expected brush collapse behavior but minimal dependence of adhesion on contact time, highlighting a clear mechanistic distinction from the polymeric counterion case. Collectively, these results demonstrate that molecular size, configurational restriction, and confinement-induced rearrangement, rather than charge valency alone, govern the interaction behavior of macromolecular counterions at brush interfaces. This work provides new insight into the molecular origins of adhesion and the regulation of interfacial interactions in charged polymer brush systems.

Park, Jinwoo [Argonne National Laboratory , , , ,;↗

DOE OSTI2023

Enhanced Optical Contrast and Switching in Near‐Infrared Electrochromic Devices by Optimizing Conjugated Polymer Oligo(Ethylene Glycol) Sidechain Content and Gel Electrolyte Composition

Abstract A detailed investigation addressing the effects of functionalizing conjugated polymers with oligo(ethylene glycol) (EG n ) sidechains on the performance and polymer‐electrolyte compatibility of electrochromic devices (ECDs) is reported. The electrochemistry for a series of donor‐acceptor copolymers having near‐infrared (NIR)‐optical absorption, where the donor fragment is 3,4‐ethylenedioxythiophene (EDOT) or an EG n functionalized bithiophene (g2T) and the acceptor fragment is diketopyrrolopyrrole (DPP) functionalized with branched alkyl or EG n sidechains, is extensively probed. ECDs are next fabricated and it is found that EG n sidechain incorporation must be finely balanced to promote polymer‐electrolyte compatibility and provide efficient ion exchange. Proper electrolyte‐cation pairing and polymer structural tuning affords a 2x increase in optical contrast (from 12% to 24%) and >60x reduction in switching time (from 20 to 0.3 s). Atomic force microscopy (AFM)/grazing incidence wide‐angle X‐ray scattering (GIWAXS) characterization of the polymer film morphology/microstructure reveals that an over‐abundance of EG n sidechains generates large polymer crystallites, which can suppress ion exchange. Lastly, time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) indicates sidechain/electrolyte identity does not influence the electrolyte penetration depth into the films, and EG n sidechain inclusion increases electrolyte cation uptake. The material structural design insight and guidelines regarding the polymer‐electrolyte ion insertion/expulsion dynamics reported here should be of significant utility for developing next‐generation mixed ionic‐electronic conducting materials.

Chemistry↗

DOE OSTI2027

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗

DOE OSTI2026

Dynamics of Radiation Damage Buildup in Ultrathin Hexagonal Boron Nitride Films under Ion Bombardment

Two-dimensional hexagonal boron nitride (hBN) is attractive for several emerging applications. Ion bombardment can be used to modify the hBN properties. However, the understanding of radiation damage buildup in hBN remains limited. Here, we investigate the effects of the dose rate and ion mass on radiation damage buildup by studying 40 nm-thick hBN films bombarded at room temperature with 500 keV 4 He, 15 N, 40 Ar, and 129 Xe ions and comparing with results for ion bombardment of polycrystalline hBN ceramics. Raman spectroscopy is used to quantify damage buildup, and transmission electron microscopy is used for microstructural analysis. Experiments are complemented by molecular dynamics simulations of the formation and evolution of point defects. Lighter ions are found to be more efficient at disordering hBN than heavier ions. This observation points to a critical role of intracascade defect processes. In contrast, a negligible dose rate effect observed suggests limited intercascade defect dynamic annealing processes for these irradiation conditions. These findings provide a fundamental basis for hBN defect engineering.

2D materials↗

DOE OSTI2026

Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES)

The C2SAGES project evaluated the feasibility of a community geothermal system for the planned Windy Ridge affordable housing development in Hinesburg, Vermont. Led by GTI Energy with Vermont Gas Systems, LN Consulting, NREL, and Frontier Energy, the work assessed technical design, energy performance, costs, business models, community engagement, maintenance, workforce development, and permitting. The proposed system was designed to serve 100% of the development’s heating, cooling, and domestic hot water loads. Compared with a baseline using air-source heat pumps and natural gas water heating, the geothermal system was estimated to reduce HVAC and domestic hot water energy use by about 45% to 48%, lower operating and maintenance costs, and reduce 30-year life-cycle costs by 37% for Phase 1 and 10% for Phase 2. Technical testing and modeling indicated that the Windy Ridge site is suitable for a community-scale geothermal system. The project also developed borehole field layouts, piping concepts, pump house designs, controls, maintenance plans, and supporting engineering drawings. The business model analysis found that first cost, ownership structure, and customer affordability remain major deployment challenges. Utility-led maintenance and operation were viewed favorably, but traditional utility cost-recovery models may require subsidy or revised financing structures to be practical for affordable housing. Community engagement highlighted the need for clear public education, transparent financing, reliable long-term maintenance, trained technicians, and the potential to pair geothermal systems with weatherization. Overall, the report concludes that community geothermal is technically feasible and offers meaningful energy, emissions, and life-cycle cost benefits, but broader deployment will depend on workable financing models and workforce readiness.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Specific Bacterial Taxa and Their Metabolite, DHPS, May Be Linked to Gut Dyshomeostasis in Patients with Alzheimer’s Disease, Parkinson’s Disease, and Amyotrophic Lateral Sclerosis

Background: Neurodegenerative diseases (NDDs) are multifactorial disorders frequently associated with gut dysbiosis, oxidative stress, and inflammation; however, the pathophysiological mechanisms remain poorly understood. Methods: Using untargeted mass spectrometry-based metabolomics and 16S sequencing of human stool, we investigated bacterial and metabolic dyshomeostasis in the gut microbiome associated with early disease stages across three NDDs—amyotrophic lateral sclerosis (ALS), Alzheimer’s disease (AD), Parkinson’s disease (PD)—and healthy controls (HC). Results: We discovered a previously unrecognized link between a microbial-derived metabolite with an unknown role in human physiology, 2,3-dihydroxypropane-1-sulfonate (DHPS), and gut dysbiosis in NDDs. DHPS was downregulated in AD, ALS, and PD, while bacteria involved in DHPS metabolism, Eubacterium and Desulfovibrio, were increased in all disease cohorts. Additionally, select taxa within the Clostridia class had strong negative correlations to DHPS, suggesting a potential role in DHPS metabolism. A catabolic product of DHPS is hydrogen sulfide, and when in excess, it is known to promote inflammation, oxidative stress, mitochondrial damage, and gut dysbiosis, known hallmarks of NDDs. Conclusions: These findings suggest that cryptic sulfur metabolism via DHPS is a potential missing link in our current understanding of gut dysbiosis associated with NDD onset and progression. As this was a hypothesis generating study, more work is needed to elucidate the role of DHPS in gut dysbiosis and neurodegenerative diseases.

Nutrition & Dietetics↗

DOE OSTI2026

Expanding the horizon of bio-naphtha beyond gasoline blend: property characterization and conversion opportunity assessment through technoeconomic and life-cycle analyses

Bio-naphtha, a common by-product of biorefineries, is expected to experience substantial growth in supply due to increasing demands for renewable diesel and synthetic aviation fuel (SAF). However, demand for bio-naphtha itself as a gasoline blendstock is limited because of the electrification of light-duty vehicles. This work investigated valorization opportunities for bio-naphtha from catalytic fast pyrolysis, hydrothermal liquefaction, Fischer–Tropsch synthesis, and hydrotreated esters and fatty acids pathways. These opportunities include producing polymer-grade olefin via steam cracking, SAF via steam cracking followed by olefin oligomerization, and renewable aromatics benzene, toluene, and xylene (BTX) and hydrogen via catalytic reforming. Process models were developed in Aspen Plus V14 and Aspen HYSYS V14 to calculate the mass and energy balances for each conversion step. Technoeconomic assessment and life-cycle analysis were conducted to evaluate the minimum fuel/product selling price, conversion cost, and life-cycle CO2 equivalent (CO2e) emission reduction. Technoeconomic assessment results suggest a minimum fuel/product selling price as low as $1.9/kg of olefins, $6.30/gal of gasoline-equivalent SAF, and $1.2/kg of BTX without any incentives. For all pathways, these prices are dominated by bio-naphtha feedstock costs, which account for at least 76% of the total cost. Compared with petroleum baselines, bio-naphtha-derived SAF, olefins, and BTX can achieve significant CO2e emission reductions from the use of renewable carbon resources. The results of life-cycle analysis and subsequent technoeconomic assessment, incorporating carbon credits, indicate the economic viability of using bio-naphtha for polymer-grade olefin and BTX production, with product costs comparable to market prices.

Xu, Shuang↗

DOE OSTI2026

SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization

Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.

08 HYDROGEN↗

DOE OSTI2026

Envelope-driven comfort risk in residential demand response

Residential demand response (DR) is a valuable resource for grid reliability, but remains challenging because the highly heterogeneous residential building stock leads to widely varying and hard-to-predict load and comfort responses during DR events. Although prior research has estimated the technical potential of DR-capable technologies for achieving energy demand savings, little is known about how they affect thermal comfort. In particular, it remains unclear how indoor thermal conditions due to DR depend on the thermal envelope characteristics of the housing stock. To address this gap, this study provides a systematic, location-specific assessment of indoor thermal performance during DR-events across the US housing stock using both typical DR weather data and detailed building metadata. We evaluate how envelope characteristics influence indoor temperatures during realistic simulated summer and winter DR events across 37 US locations, applying both temperature threshold and rate of temperature change criteria to estimate region-level probabilities of discomfort. Additionally, we show the impact of distinct weather patterns that intensify or abate thermal stress on comfort outcomes. Results show a near-universal overheating risk in summer DR events, where comfort outcomes are strongly influenced by rapid risk of comfort violations. In contrast, overall winter DR discomfort risk is lower, risk escalation is more gradual and shows greater sensitivity to event duration. These findings offer a data-driven quantification of comfort risk across diverse climates and building envelopes, demonstrating the need for region-specific DR scheduling and discomfort mitigation strategies tailored to local weather patterns and the performance of existing residential buildings.

Demand response↗

DOE OSTI2023

Host galaxies of ultra-strong Mg ii absorbers at z ∼ 0.7

ABSTRACT We report spectroscopic identification of the host galaxies of 18 ultra-strong Mg ii systems (USMg ii) at 0.6 ≤ z ≤ 0.8. We created the largest sample by merging these with 20 host galaxies from our previous survey within 0.4 ≤ z ≤ 0.6. Using this sample, we confirm that the measured impact parameters ($\rm 6.3\leqslant D[kpc] \leqslant 120$ with a median of 19 kpc) are much larger than expected, and the USMg ii host galaxies do not follow the canonical $\rm {\it W}_{2796}-{\it D}$ anticorrelation. We show that the presence and significance of this anticorrelation may depend on the sample selection. The $\rm {\it W}_{2796}-{\it D}$ anticorrelation seen for the general Mg ii absorbers show a mild evolution at low $\rm W_{2796}$ end over the redshift range 0.4 ≤ z ≤ 1.5 with an increase of the impact parameters. Compared to the host galaxies of normal Mg ii absorbers, USMg ii host galaxies are brighter and more massive for a given impact parameter. While the USMg ii systems preferentially pick star-forming galaxies, they exhibit slightly lower ongoing star-forming rates compared to main sequence galaxies with the same stellar mass, suggesting a transition from star-forming to quiescent states. For a limiting magnitude of mr < 23.6, at least 29 per cent of the USMg ii host galaxies are isolated, and the width of the Mg ii absorption in these cases may originate from gas flows (infall/outflow) in isolated haloes of massive star forming but not starbursting galaxies. We associate more than one galaxy with the absorber in $\ge 21~{{\ \rm per\ cent}}$ cases, where interactions may cause wide velocity spread.

Astronomy & Astrophysics↗

DOE OSTI2023

Constraints on Nuclear Symmetry Energy Parameters

A review is made of constraints on the nuclear symmetry energy parameters arising from nuclear binding energy measurements, theoretical chiral effective field predictions of neutron matter properties, the unitary gas conjecture, and measurements of neutron skin thicknesses and dipole polarizabilities. While most studies have been confined to the parameters 𝑆 𝑉 and L, the important roles played by, and constraints on 𝐾 sym , or, equivalently, the neutron matter incompressibility 𝐾 𝑁 , are discussed. Strong correlations among 𝑆 𝑉 ,𝐿, and 𝐾 𝑁 are found from both nuclear binding energies and neutron matter theory. However, these correlations somewhat differ in the two cases, and those from neutron matter theory have smaller uncertainties. To 68% confidence, it is found from neutron matter theory that 𝑆 𝑉 = 32.0 ± 1.1 MeV, 𝐿 = 51.9 ± 7.9 MeV and 𝐾 𝑁 = 152.2 ± 38.1 MeV. Theoretical predictions for neutron skin thickness and dipole polarizability measurements of the neutron-rich nuclei 48 Ca, 120 Sn, and 208 Pb are compared to recent experimental measurements, most notably the CREX and PREX neutron skin experiments from Jefferson Laboratory. By themselves, PREX I+II measurements of 208 Pb and CREX measurement of 48 Ca suggest 𝐿 = 121 ± 47 MeV and 𝐿 = −5 ± 40 MeV, respectively, to 68% confidence. However, we show that nuclear interactions optimally satisfying both measurements imply 𝐿 = 53 ± 13 MeV, nearly the range suggested by either nuclear mass measurements or neutron matter theory, and is also consistent with nuclear dipole polarizability measurements. This small parameter range implies 𝑅 1.4 = 11.6 ± 1.0 km and Λ 1.4 = 228$^{+148}_{−90}$, which are consistent with NICER X-ray and LIGO/Virgo gravitational wave observations of neutron stars.

neutron skins↗

DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗

DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2027

Low thermal budget dopant activation in shallow junctions of implanted Si via Tunable plasma enhanced annealing

Low thermal budget (LTB) annealing has become increasingly critical for shallow junctions as semiconductor devices are scaled down. Plasma Enhanced Annealing (PEA) has emerged as a promising LTB annealing process, however its mechanisms remain unclear. In this study, arsenic and boron implanted silicon wafers were subjected to helium or argon ion bombardment generated by annealing plasmas under controlled ion energies and doses in a home-built annealing reactor. The structural and electrical changes were characterized by four-point probe measurements, secondary ion mass spectrometry, Raman spectroscopy and spectroscopic ellipsometry. A pronounced reduction in sheet resistance, accompanied by a decrease in the damage layer thickness demonstrates the surface selective annealing capability of PEA process. Comparative studies using He and Ar plasmas and related ion bombardment reveal a dependence of activation efficacy on both ion species and the dopant implantation profile. Furthermore, a multi-step annealing mechanism is proposed, consisting of an initial (low dose) structural recovery stage followed by the generation of “extended” defects and concluded at large doses by an enhanced dopant activation. The proposed annealing kinetics are thought to be controlled by ion flux, dose and energy. Finally, these results highlight PEA as an effective, surface-selective, LTB approach for shallow-dopant activation and near surface annealing, and provide mechanistic insights into PEA processes.

dopant activation↗

DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗