Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete Element Method”

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.

700 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Novel Protection Scheme for Unbalanced Faults in Inverter Dominated Networks: A Computationally Efficient Algorithm for Entry-Level Relays

Microgrids are now a common practice in distribution systems to increase resilience and reliability. However, microgrid protection remains a critical challenge, considering its requirement to operate in both grid connected and islanded, and the variability in fault characteristics under each mode of operation. This paper presents unbalanced power (S unb ) based fault detection algorithm, which considers local voltage and current unbalances to determine faults in the system. S unb is a computationally efficient fault detection algorithm that is suitable for implementation in the programmable logic of entry level protective relays. In addition, the difference in current and voltage unbalance (D n ) is used to determine the fault type. The proposed method demonstrates high sensitivity and selectivity for line-to-ground (LG), line-to-line (LL), and double line-to-ground (LLG) faults, representing the most common faults in distribution systems. It also allows relay coordination with upstream and downstream protection devices in both island and grid connected operation, while preserving grading margins. The same pickup and time multiplier settings of a particular relay for both modes of operation eliminates the need for adaptive settings, which rely on communication networks. Validation was performed with a hardware-in-the-loop (HIL) setup using Typhoon HIL real time simulator interfaced with three entry-level, SEL 751 relays. Results confirmed the algorithm’s ability to discriminate fault conditions, and determine the fault type under both operating modes, maintain fast detection times, and ensure proper protection coordination.

fault classification↗
DOE OSTI2026

Progress and Perspectives: Zirconium Electrodeposition from Different Electrolytes

Zirconium (Zr) possesses outstanding properties, including exceptional chemical resistance and a high melting point, and is therefore desirable for use in a wide variety of challenging environments, such as in nuclear reactors and the chemical processing industry. To minimize the amount of Zr required for any given application, developing methods for generating metallic Zr coatings is highly advantageous. Owing to Zr’s highly negative reduction potential, electrodeposition in traditional solvents at near ambient temperatures remains challenging and thus not well understood. Due to the extreme conditions required to deposit this metal, extensive work has been conducted in molten salt electrolytes, however the broad applicability of this methodology is limited due to its corrosivity. The primary focus of this article is to present an overview of Zr’s electrochemical behavior and to consolidate the efforts of researchers in exploring electrodeposition techniques for Zr involving aqueous, organic, ionic liquid, deep eutectic, and molten salt solvents. With this information, we highlight trends across solvent systems and opportunities for future research.

42 ENGINEERING↗
DOE OSTI2026

Improved Protein Semi-Synthesis Enables Biophysical Studies of Thioamide Destabilization of β-Sheet Interactions

Abstract Thioamides are natural post-translational modifications of the peptide backbone and can be introduced synthetically to probe protein folding or functionalize peptides for translational applications. In this work, we demonstrate that thioamide-containing peptides with C-terminal thioesters can be efficiently generated using Knorr pyrazole activation and used in subsequent native chemical ligation reactions to generate thioamide-containing proteins. We compare this method to acyl azide activation and find that both routes provide similar yields. We also investigate ultrasound-mediated desulfurization of the ligation site cysteine for potential advantages over chemical radical initiators. Scaling up our syntheses allows us to study thioamide perturbations to the β-sheet region of the B1 domain of protein G (GB1) as well as β-strand interactions in amyloid fibrils of the Parkinson’s disease protein α-synuclein. In both contexts, we observe dramatic destabilization of the β-sheet networks, manifested in decreased GB1 thermal stability and altered folding and slowed aggregation of α-synuclein. These findings illustrate the impact that a single atom substitution can have on cooperative hydrogen-bonding networks and prompt future study of both systems.

Yanagawa, Evan S. K. [University of Pennsylvania ,↗
DOE OSTI2025

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗
DOE OSTI2025

Photon-Counting Interferometry to Detect Geontropic Space-Time Fluctuations with GQuEST

The gravity from the quantum entanglement of space-time (GQuEST) experiment uses tabletop-scale Michelson laser interferometers to probe for fluctuations in space-time. We present a practicable interferometer design featuring a novel photon-counting readout method that provides unprecedented sensitivity, as it is not subject to the interferometric standard quantum limit. We evaluate the potential of this design to measure space-time fluctuations motivated by recent “geontropic” quantum gravity models. The accelerated accrual of Fisher information offered by the photon-counting readout enables GQuEST to detect the predicted quantum gravity phenomena within measurement times at least 100 times shorter than equivalent conventional interferometers. The GQuEST design, thus, enables a fast and sensitive search for signatures of quantum gravity in a laboratory-scale experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

DEPRECATED AI-Batt-OS (Autonomous Identification of Battery Life Models - Open Source) [SWR 21-17]

DEPRECATED. This repository was archived by the owner on Jun 30, 2026. It is now read-only. Open source implementation of some of the methods utilized by AI-Batt, a battery lifetime modeling and analysis toolkit provided by the National Laboratory of the Rockies (NLR). This software demonstrates the use of bi-level optimization and symbolic regression techniques to semi-autonomously identify algebraic models predicting the capacity fade of lithium-ion batteries during calendar aging. Modeling the degradation of batteries is a complex task, due to the difficulty in separating the time-dependent and time-independent factors impacting cell level degradation, across multiple data series with different numbers of measurements and/or data quality. Bi-level optimization enables model parameters to be optimized to either the entire data set or to individual data series, allowing statistical disambiguation of global behaviors (data series independent) and local behaviors (data series dependent). Symbolic regression is used to automatically search for optimal low-dimesional models predicting the variation of locally optimized parameters versus time-independent experimental variables from millions of possible models, resulting in a more accurate and repeatable model identification process than is possible by a manual search. The provided tools also implement cross-validation and bootstrap resampling schemes, empowering statistical model comparison/selection and quantification of model uncertainties. An example script replicates the results from the manuscript "Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning", submitted to ECS. All code is written in MATLAB. Requires the Statistics and Machine Learning Toolbox. Contact Dr. Paul Gasper at Paul.Gasper@nlr.gov for any questions.

Gasper, Paul [National Renewable Energy Lab. (NREL↗
DOE OSTI2025

Crystal Growth Scale-Up and Stability of RbSr 2 X 5 :Eu Scintillators

The discovery and development of new scintillation materials support national security needs. In these applications, large volumes of scintillator crystals are needed to achieve efficient screening for contraband. Therefore, one important step in the discovery of new scintillator compositions is testing their feasibility for scale-up and their stability. In this work, high-quality Ø22 mm crystals of two new scintillators RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu were grown via the Vertical Bridgman method, and their scintillation properties were characterized. Here, the Ø22 mm RbSr 2 I 5 :Eu crystals could be grown with fast translation rates up to 3.5 mm/h. Both RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu had high scintillation performance, including light yields of 46,000 and 61,000 ph/MeV, respectively, for Ø22 × 35 mm crystals. Additionally, properties related to physical stability were investigated, including the coefficients of thermal expansion via high-temperature X-ray diffraction (HTXRD) as well as moisture sensitivity. HTXRD confirmed the absence of solid–solid phase transitions and showed that RbSr 2 Br 5 had minimal thermal expansion anisotropy compared to RbSr 2 I 5 and some other inorganic metal halide scintillators, which favors the growth of large-sized crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Cross‐Linked Flexible Metaferroelectrolyte Regulated by 2D/2D Perovskite Heterostructures for High‐Performance Compact Solid‐State Sodium Batteries

Abstract To address the issues of limited ionic conductivity and poor interface stability at room and low temperatures in solid‐state electrolytes, a robust intrinsic ferroelectrolyte or nanoferroelectrolyte strategy for engineering solid‐state flexible ferroelectric composite electrolytes utilizing strongly coupled intrinsic ion conducting 2D/2D sodium‐rich anti‐perovskite (NaRAP)/ferroelectric perovskite heterostructures is introduced. Herein, highly scalable PVDF‐based metaferroelectrolytes with Na 2.99 Ba 0.005 OCl/Ca 2 Na 2 Nb 5 O 16 − (CNNO − ) nanosheets into a ferroelectric poly(vinylidene fluoride‐co‐hexafluoropropylene) (PVDF‐HFP) matrix, through an in situ cross‐linking and spontaneous bridging method, for compact solid‐state sodium batteries (SSBs), are reported. Benefiting from unique well‐dispersed 3D ferroelectric coupled network and the Na 2.99 Ba 0.005 OCl/CNNO − ‐induced PVDF‐HFP ferroelectric β phase, the Na + flux is regulated, thereby inhibiting Na dendrite growth at the interface. Notably, the optimized PH‐5% NC metaferroelectrolyte exhibits rapid ion transport (1.11 × 10 −4 S cm −1 at 25 °C), a wide electrochemical window (> 4.8V), superior conformal mechanical compatibility, improved flexibility, good elasticity and flame retardancy. The solid‐state Na 3 V 2 (PO 4 ) 3 /PH‐5% NC/Na batteries present a stable cycling performance (remaining 56.4 mAh g −1 after 500 cycles at 1 C) even at 0 °C, potential for cost‐effective, safe, stable and compact SSB energy storage over 600 Wh L −1 , vastly surpassing 365 Wh L −1 of the current commercial sodium‐ion liquid‐electrolyte batteries.

Chemistry↗
DOE OSTI2026

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗
NASA NTRS2003

Concept to Reality: Contributions of the Langley Research Center to US Civil Aircraft of the 1990s

This document is intended to be a companion to NASA/SP-2000-4519, Partners in Freedom: Contributions of the Langley Research Center to U.S. Military Aircraft of the 1990s . Material included in the combined set of volumes provides informative and significant examples of the impact of Langley's research on U.S. civil and military aircraft of the 1990s. This volume, 'Concept to Reality: Contributions of the NASA Langley Research Center to U.S. Civil Aircraft of the 1990s', highlights significant Langley contributions to safety, cruise performance, takeoff and landing capabilities, structural integrity, crashworthiness, flight deck technologies, pilot-vehicle interfaces, flight characteristics, stall and spin behavior, computational design methods, and other challenging technical areas for civil aviation. The contents of this volume include descriptions of some of the more important applications of Langley research to current civil fixed-wing aircraft (rotary-wing aircraft are not included), including commercial airliners, business aircraft, and small personal-owner aircraft. In addition to discussions of specific aircraft applications, the document also covers contributions of Langley research to the operation of civil aircraft, which includes operating problems. This document is organized according to disciplinary technologies, for example, aerodynamics, structures, materials, and flight systems. Within each discussion, examples are cited where industry applied Langley technologies to specific aircraft that were in operational service during the 1990s and the early years of the new millennium. This document is intended to serve as a key reference for national policy makers, internal NASA policy makers, Congressional committees, the media, and the general public. Therefore, it has been written for a broad general audience and does not presume any significant technical expertise. An extensive bibliography is provided for technical specialists and others who desire a more in-depth discussion of the contributions.

NASA History Series↗
NASA NTRS2026

Swamp Works Regolith Compaction Technologies

While the level of compaction below the lunar surface increases quickly after only a few cm of depth, in many cases during a construction mission there will be a need to excavate and transport regolith to a new location for cut-and-fill or horizontal construction of structures such as berms. In these cases to achieve high levels of bulk density, compaction must be per-formed. Additionally, in some cases surface technolo-gies such as systems that sinter/melt the surface may desire the maximum possible compaction at the sur-face to improve melting/heating performance and the final material strength properties. Kennedy Space Center’s (KSC) Swamp Works has developed two means of compaction, lunar and mar-tian compaction. Planetary Autonomous Compaction Technology (PACT) which is part of the Multifunction End Effector for Regolith Compaction Acquisition and Transfer (MEERCAT [1]) robotic arm end effector system’s capabilities and the Site Preparation Tooling for Operations on Mobility Platforms (STOMP [2]) vibratory roller compactor. PACT on MEERCAT has been demonstrated to a TRL 5 and STOMP to a TRL 4 in ambient testing. The results of PACT on MEERCAT and STOMP testing will be shared with results for various simulants including BP-1, ICN-LHT-1G (aka CSM-LHT-1G), RDW-LHT-1GH (a simulant developed for the Mason Tipping Point to match characteristics of ICN-LHT-1G), and Exolith LHS-1E. This will also include discus-sions on methods used to verify relative density before and after compaction and means to verify density effects below depth. To calculate relative density, maximum and minimum densities for simulants were taken from literature [3] [4] [5] and additional lab test-ing (publication in work).

vibration↗
DOE OSTI2026

ClassNMSW- a real-time classification approach for non-recycled municipal solid waste using hyperspectral imaging

Real-time classification of non-recycled municipal solid waste (NMSW) is essential for efficient valorization. This study introduces ClassNMSW, a comprehensive framework for classifying 22 NMSW subclasses under industrial constraints by using hyperspectral imaging (HSI). A primary innovation of this work is the development of a variance-controlled spectral extraction algorithm. Unlike traditional methods that rely on simple averaging, this approach systematically investigates the extent of pixel extraction to minimize the loss of critical chemical information while maximizing data reduction thus ensuring high spectral fidelity with low computational cost. The approach developed in this work integrates automated, computer-vision-based background removal, eliminating the need for the manual thresholding common in current literature. To resolve ambiguities among chemically similar subclasses, a tiered classification and multi-camera fusion strategy (NIR17 and NIR22) is implemented. Results demonstrate that ClassNMSW achieves an object-wise weighted accuracy of 98.70% for single-sensor configurations and 100% under sensor fusion. A novel rolling-window strategy satisfies desired end-to-end latency of <2 s, satisfying the strict deterministic requirements of high-speed industrial sorting environments. The ClassNMSW framework provides a scalable foundation for advancing circularity and resource recovery in large-scale waste valorization operations.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

High-Intensity UV Exposure for the Rapid Screening of Silicon Photovoltaic Architectures

Advanced Si photovoltaic architectures incorporate different materials and processing pathways that influence degradation modes. Ultraviolet-induced degradation (UVID) is an understudied degradation mode for advanced cell architectures and is of increasing concern to industry due to growing adoption of UV-transparent encapsulation and bifacial technologies. In order to adopt new and evolving technologies confidently, novel component materials and processing techniques must be evaluated and designed for long-term stability, in addition to the conventional design focus on efficiency. In this work, a study protocol framework is presented for the rapid screening of unencapsulated devices against UVID. Unencapsulated passivated emitter rear contact (PERC) and tunnel oxide passivated contact (TOPCon) devices were aged under different UV irradiance intensities and measured via conventional nondestructive electrical characterization methods to assess performance degradation. Based on the results, protocol efficacy and recommendations for further study are discussed. As a result, this work is part of a broader effort to develop rapid screening processes that cut across architectures and exposure conditions to aid module manufacturers in vetting new materials choices for long-term stability.

Accelerated exposure↗
DOE OSTI2026

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion↗
DOE OSTI2024

Coercivity enhancement in hematite/permalloy heterostructures across the Morin transition

Interfacial effects between antiferromagnetic (AFM) and ferromagnetic (FM) materials have long been a center of magnetism studies. Aside from the exchange bias occurring at the AFM/FM interface, controlling the coercivity is another significant topic in magnetic recordings. The coercivity of FM materials is often determined through varying grain size, alloy composition, density of defects, etc., which is set during material growth and offers limited room for modification after growth. Hematite (α - Fe 2 O 3 ) is an AFM material that undergoes a temperature-controlled spin-flip transition, the so-called Morin transition. This transition gives an extra degree of freedom making hematite an intriguing component to study the exchange coupling when interfaced with an FM material. Here, in this work, changes in the magnetic properties of soft magnetic permalloy (Ni 81 Fe 19 , or Py) thin films grown on hematite were studied across the Morin transition. Surprisingly, these samples showed a remarkable change in coercivity during the Morin transition. We attribute this effect to the magnetic domain mixture of hematite during the Morin transition. Our findings present a novel method of controlling the coercivity of plain ferromagnetic thin films.

Antiferromagnet↗
DOE OSTI2026

Rinse-Free, Sodium-Efficient Synthesis of O3-Type Layered Oxide Materials Enabled by Acetate Precursors

Sodium-ion batteries (SIBs) are a sustainable alternative to lithium-ion systems for global electrification, with O3-type layered oxide cathodes offering high specific capacity and feasibility of scalable synthesis. Industrial co-precipitation synthesis of these cathodes typically uses transition metal sulfates, requiring extensive water rinsing to remove Na 2 SO 4 impurities, a process that consumes significant water and risks residual inactive phases if incomplete. Here, this work introduces a rinse-free, resource-efficient approach using metal acetate precursors. Residual sodium acetate in non-rinsed precursors decomposes during sintering to generate Na 2 CO 3 in situ, partially substituting an external sodium resource (e.g., NaOH and Na 2 CO 3 ) and reducing its consumption by ∼18–20%. Phase-pure O3-Na 1.0 Ni 1/3 Fe 1/3 Mn 1/3 O 2 (NFM111) cathodes synthesized via this method exhibit microstructure and electrochemical performance comparable to rinsed sulfate-derived counterparts, with initial capacities of 141 mAh g −1 at C/20 (7.5 mA g −1 ). By eliminating rinsing and minimizing sodium reagent use, this acetate-based route enhances sustainability and scalability of layered oxide production for SIBs.

acetate vs. sulfate↗
DOE OSTI2026

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗
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↗