Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

How generic is iodide-tagging photoelectron spectroscopy: An extended investigation on the Gly· X - (Gly = glycine, X = Cl or Br) complexes

We report a joint negative ion photoelectron spectroscopy (NIPES) and quantum chemical computational study on glycine-chloride/bromide complexes (denoted Gly·X - , X = Cl/Br) in close comparison to the previously studied Gly·I - cluster ion. Combining experimental NIPE spectra and theoretical calculations, various Gly·X - complexes were found to adopt the same types of low-lying isomers, albeit with different relative energies. Despite more congested spectral profiles for Gly·Cl - and Gly·Br - , spectral assignments were accomplished with the guidance of the knowledge learned from Gly·I - , where a larger spin–orbit splitting of iodine afforded well-resolved, recognizable spectral peaks. Three canonical plus one zwitterionic isomer for Gly·Cl - and four canonical conformers for Gly·Br - were experimentally identified and characterized in contrast to the five canonical ones observed for Gly·I - under similar experimental conditions. Taken together, this study investigates both genericity and variations in binding patterns for the complexes composed of glycine and various halides, demonstrating that iodide-tagging is an effective spectroscopic means to unravel diverse ion-molecule binding motifs for cluster anions with congested spectral bands by substituting the respective ion with iodide.

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

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

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

Learning the stable and metastable phase diagram to accelerate the discovery of metastable phases of boron

Boron, an element of captivating chemical intricacy, has been surrounded by controversies ever since its discovery in 1808. The complexities of boron stem from its unique position between metals and insulators in the Periodic Table. Recent computational studies have shed light on some of the stable boron allotropes. However, the demand for multifunctionality necessitates the need to go beyond the stable phases into the realm of metastability and explore the potentially vast but elusive metastable phases of boron. Traditional search for stable phases of materials has focused on identifying materials with the lowest enthalpy. Here, we introduce a workflow that uses reinforcement learning coupled with decision trees, such as Monte Carlo tree search, to search for stable and metastable boron phases, with enthalpy as the objective. We discover new boron metastable phases and construct a phase diagram that locates their phase space (T, P) at different levels of metastability (ΔG) from the ground state and provides useful information on the domains of relative stability of the various stable and metastable boron phases.

42 ENGINEERING↗
DOE OSTI2024

Deprotonated sulfamic acid and its homodimers: Does sulfamic acid adopt zwitterion during cluster growth?

Here we present a joint experimental and computational study on the geometric and electronic structures of deprotonated sulfamic acid (SA) clusters [(SA) n –H] – (n = 1, 2) employing negative ion photoelectron spectroscopy and high-level ab initio calculations. The photoelectron spectra provide the vertical/adiabatic detachment energy (VDE/ADE) of the sulfamate anion (SM – ) H 2 N•SO 3 – at 4.85 ± 0.05 and 4.58 ± 0.08 eV, respectively, and the VDE and ADE of the SM – •SA dimer at 6.41 ± 0.05 and 5.87 ± 0.08 eV, respectively. The significantly increased electron binding energies of the dimer confirm the enhanced electronic stability upon the addition of one SA molecule. The CCSD(T)-predicted VDEs/ADEs agree excellently with the experimental data, confirming the identified structures as the most stable ones. Two types of dimer isomers possessing different hydrogen bonding (HB) motifs are identified, corresponding to SM – binding to a zwitterionic SA (SM – •SA z ) and a canonical SA (SM – •SA c ), respectively. Two N–H$\cdots$O HBs and one superior O–H$\cdots$O HB are formed in the lowest-lying SM – •SA c , while SM – •SA z has three moderate N–H$\cdots$O HBs, with the former being 4.71 kcal/mol more stable. Further theoretical analyses reveal that the binding strength advantage of SM – •SA c over SM – •SA z arises from its significant contributions of orbital interactions between fragments, illustrating that sulfamate strongly interacts with its parent SA acid and preferably chooses the canonical SA in the subsequent cluster formations. Given the prominent presence of SA, this study provides the first evidence that the canonical dimer model of sulfamic acid should exist as a superior configuration during cluster growth.

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

Compositionally complex carbide ceramics: A perspective on irradiation damage

Extensive experimental and computational studies have demonstrated outstanding physical and chemical properties of the novel materials of compositionally complex carbides (CCCs), enabling their promising applications in advanced fission and fusion energy systems. This perspective provides a comprehensive overview of radiation damage behavior reported in the literature to understand the fundamental mechanisms related to the impact of multi-principal metal components on phase stability, irradiation-induced defect clusters, irradiation hardening, and thermal conductivity of compositionally complex carbides. Several future research directions are recommended to critically evaluate the feasibility of designing and developing new ceramic materials for extreme environments using the transformative “multi-principal component” concept. Compared to the existing materials for nuclear applications including stainless steels, nickel alloys, ZrC, SiC, and potentially high-entropy alloys, as well as certain other compositionally complex ceramic families. CCCs appear to be more resistant to amorphization, growth of irradiation defect clusters, and void swelling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2025

Effect of nozzle translation on a bounded vortex flow

A bounded vortex flow is generated by a nozzle that combines azimuthally tilted downward jets in a circular array and a central suction port. When the nozzle is directed toward a bottom surface (called the impingement surface), the flow develops a strong intake vortex spanning from the impingement surface to the suction port, which acts to generate strong shear stress on the impingement surface. This strong shear stress can aid in removing particles attached to the impingement surface or in enhancement of heat or mass transfer, where in typical applications, the nozzle would be translated along the impingement surface. This paper reports on a computational study of the effect of nozzle translation on the bounded vortex flow field. The nozzle translation induces a bending of the bounded vortex in the downstream flow direction, which becomes larger with increase in the nozzle translation speed. We identify an interesting oscillatory phenomenon by which the bounded vortex can resist the bending induced by the translational flow and rotate about a mean position downstream of the suction port. At sufficiently high translation speeds, the oscillatory process can no longer resist the downwash effect from nozzle translation and the vortex eventually detaches from the impingement surface, limiting the effectiveness of this flow device.

42 ENGINEERING↗
DOE OSTI2021

Cryo-EM structures of Escherichia coli cytochrome bo 3 reveal bound phospholipids and ubiquinone-8 in a dynamic substrate binding site

Two independent structures of the proton-pumping, respiratory cytochrome bo3 ubiquinol oxidase (cyt bo 3 ) have been determined by cryogenic electron microscopy (cryo-EM) in styrene–maleic acid (SMA) copolymer nanodiscs and in membrane scaffold protein (MSP) nanodiscs to 2.55- and 2.19-Å resolution, respectively. The structures include the metal redox centers (heme b, heme o 3 , and Cu B ), the redox-active cross-linked histidine–tyrosine cofactor, and the internal water molecules in the proton-conducting D channel. Each structure also contains one equivalent of ubiquinone-8 (UQ8) in the substrate binding site as well as several phospholipid molecules. The isoprene side chain of UQ8 is clamped within a hydrophobic groove in subunit I by transmembrane helix TM0, which is only present in quinol oxidases and not in the closely related cytochrome c oxidases. Both structures show carbonyl O1 of the UQ8 headgroup hydrogen bonded to D75I and R71 I . In both structures, residue H98 I occupies two conformations. In conformation 1, H98 I forms a hydrogen bond with carbonyl O4 of the UQ8 headgroup, but in conformation 2, the imidazole side chain of H98 I has flipped to form a hydrogen bond with E14 I at the N-terminal end of TM0. We propose that H98 I dynamics facilitate proton transfer from ubiquinol to the periplasmic aqueous phase during oxidation of the substrate. Computational studies show that TM0 creates a channel, allowing access of water to the ubiquinol headgroup and to H98 I .

bioenergetics↗
DOE OSTI2021

Combining pressure and electrochemistry to synthesize superhydrides

Significance Superhydrides are a materials system where near–room-temperature superconductivity has been achieved but only at very high (megabar) pressures. This work proposes an approach that combines pressure and electrochemistry to stabilize superhydrides at moderate pressures. Through a computational study of the palladium–hydrogen system, we construct electrochemical phase diagrams and show that electrochemically synthesizing superhydrides may be possible when combined with moderate pressures. We generalize this to other binary metal superhydrides of interest for superconductivity, including La, Y, and Mg hydrides.

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

Hydrogen ejection from hydrocarbons: Characterization and relevance in soot formation and interstellar chemistry

Polycyclic aromatic hydrocarbons (PAHs) play a major role in the chemistry of combustion, pyrolysis, and the interstellar medium. Production (or activation) of radical PAHs and propagation of their resulting reactions require efficient dehydrogenation, but the preferred method of hydrogen loss is not well understood. Unimolecular hydrogen ejection (i.e., direct C─H bond fission) and bimolecular radical abstraction are two main candidate pathways. We performed a computational study to characterize the role of H ejection, particularly as a driver for radical-centric hydrocarbon-growth mechanisms and particle formation. Electronic structure calculations establish that C─H bond strengths span a broad range of energies, which can be weaker than 30 kcal/mol in some C 9 and C 13 PAH radicals. At T > 1200 K, calculated thermal rates for hydrogen ejection from weak C─H bonds at zigzag sites on PAH radicals are significantly larger than typical H-abstraction rates. These results are highly relevant in the context of chain reactions of radical species and soot inception under fuel-rich combustion conditions. Furthermore, calculated microcanonical rates that include the additional internal energy released by bond formation (e.g., ring closure to yield C 9 H 9 ) yield significantly higher rates than those associated with full thermalization. These microcanonical considerations are relevant to the astrochemical processes associated with hydrocarbon growth and processing in the low-density interstellar environment.

08 HYDROGEN↗
DOE OSTI2026

Photoenzymatic Csp 3 –Csp 3 bond formation via enzyme-templated radical–radical coupling

Cross-couplings are essential reactions in modern chemical synthesis, enabling the rapid construction of complex molecules from simple precursors. Transition metal catalysts are prized for these transformations because their reactivity and selectivity can be tuned via judicious selection of the metal and ligand. Although enzymes offer analogous opportunities for tuning via protein engineering, their application to cross-coupling remains limited, as nature relies on alternative paradigms for building molecular complexity. Here, we report the cross-coupling of alkyl halides and benzylic carboxylic acids using an engineered flavin-dependent lactate monooxygenase—a photoenzyme. The enzyme achieves this feat by exploiting the redox versatility of the flavin cofactor. Stoichiometric experiments, ultrafast spectroscopy, and computational studies support a mechanism in which photoexcited flavin quinone initiates the reaction via oxidative decarboxylation to generate a benzylic radical. The resulting flavin semiquinone can reduce the alkyl halide to form a second organic radical within the protein active site, which rapidly engages in C(sp 3 )–C(sp 3 ) bond formation. A variant was engineered to control the stereochemical outcome of this radical–radical coupling event, highlighting the ability of the protein to alter the energetic barrier for a mechanistic step that is traditionally understood to be near barrierless. This work demonstrates that the scope for nonnative reaction mechanisms in biocatalysis far exceeds previously established bounds and has potential to solve a variety of reactivity challenges in cross-coupling chemistry.

biocatalysis↗
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

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

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A novel measurement of marginal Alfvén eigenmode stability during high power auxiliary heating in JET

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2026

Evidence of topological charge polarization at active-passive interfaces in acoustically powered active liquid crystals

Theoretical and computational studies predict topological charge separation at active-passive interfaces in active nematics, but reliable experimental validation has been lacking. We utilize a synthetic acoustically energized liquid crystal to experimentally investigate the dynamics of topological defects across an active-passive interface. The interference pattern induced by the acoustic wave inside the experimental cell creates a spatial distribution of activity, resulting in the formation of active-passive interfaces with the liquid-crystalline director field aligned parallel to those interfaces. The activity gradient drives the reorientation of positively charged topological defects along its direction, causing them to migrate into the passive zone and form a layer with net topological charge. That layer “discharges” upon cessation of activity through annihilation of positive and negative topological charges, resembling the discharge of a capacitor.

Active defects↗