Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerator Simulation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Particle Acceleration in Two Converging Shocks

Observations by spacecraft such as ACE , STEREO , and others show that there are proton spectral “breaks” with energy E {sub br} at 1–10 MeV in some large CME-driven shocks. Generally, a single shock with the diffusive acceleration mechanism would not predict the “broken” energy spectrum. The present paper focuses on two converging shocks to identify this energy spectral feature. In this case, the converging shocks comprise one forward CME-driven shock on 2006 December 13 and another backward Earth bow shock. We simulate the detailed particle acceleration processes in the region of the converging shocks using the Monte Carlo method. As a result, we not only obtain an extended energy spectrum with an energy “tail” up to a few 10 MeV higher than that in previous single shock model, but also we find an energy spectral “break” occurring on ∼5.5 MeV. The predicted energy spectral shape is consistent with observations from multiple spacecraft. The spectral “break,” then, in this case is caused by the interaction between the CME shock and Earth’s bow shock, and otherwise would not be present if Earth were not in the path of the CME.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Orchestrating Spontaneous Emission With Metasurfaces: Recent Advances in Engineering Thermal, Luminescent, and Quantum Emissions

Metasurfaces have emerged as powerful tools for controlling spontaneous emission, offering unprecedented control over light-matter interactions at sub-wavelength scales. While metasurfaces are traditionally utilized for shaping coherent electromagnetic waves, they have recently extended their capabilities to control incoherent or spontaneous emission. This examines review how metasurfaces can enhance and precisely control properties of thermal, luminescent, and quantum emission. In thermal emission, metasurfaces enable control over spatial, temporal, and spin coherence, offering new possibilities for applications such as energy harvesting, radiative cooling and heat assisted ranging and detection. For luminescent emission, metasurfaces significantly improve emission rates, quantum efficiency, and directionality, driving innovations in lighting and display technologies. For controlling quantized spontaneous emission, metasurfaces are instrumental in enhancing single-photon sources and enabling novel functionalities in quantum states through photon-pair generation, which is vital for quantum communication, meteorology, and computing. Here, despite these advancements several challenges to increase the operational bandwidths, accelerate and develop simulation strategies, and fabrication complexities persist. Emerging trends are also discussed, such as dynamic metasurfaces and their integration with nanophotonic platforms, which could further expand the capabilities of light-emitting metasurfaces.

Luminescence↗
DOE OSTI2020

Thermodynamic Modeling of the Al-Ce-Cu-Mg-Si System and Its Application to Aluminum-Cerium Alloy Design

Recently discovered AlCe alloys have shown promise in a number of applications, but the propensity of Ce to react with Al and other alloying elements can complicate the phase equilibria and design approach. To solve this, the CALPHAD method is used to explore an alloy within the quinary Al-Ce-Cu-Mg-Si system by developing a thermodynamic database with self-consistent parameters. The database includes a description of all 10 binary systems and 8 ternary systems consisting of: (i) 6 Al-containing ternaries (Al-Ce-Cu, Al-Ce-Mg, Al-Ce-Si, Al-Cu-Mg, Al-Cu-Si and Al-Mg-Si); and (ii) 2 additional ternaries that include Mg and Si (i.e., Ce-Mg-Si and Cu-Mg-Si). The thermodynamic description for the Al-Ce-Mg and Al-Mg-Si systems were reassessed to ensure consistency with the binary systems and the Ce-Mg-Si system is presented for the first time and compared to theoretical data from DFT (Density Functional Theory). In addition to the ternary interactions, the quaternary compound Al3Cu2Mg9Si7 and solid solution extending from the ternary Al2CuMg phase (Al,Si)2CuMg are incorporated. The CALPHAD method is employed and leveraged through the use of a Materials Design Simulator (MDS) to accelerate the design of novel aluminum-cerium-based alloys. The combination of a CALPHAD-based framework with experimental efforts and industrial insight permits the development of three new Al-Ce alloys: Al-3.5Ce-0.4Mg-7Si (Ce-modified A356), Al-5Ce-1Cu-0.5Mg-10Si and Al-19Ce-0.9Mg-1.1Si.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Next-generation anion exchange membrane water electrolyzers operating for commercially relevant lifetimes

Alkaline anion exchange membrane (AEM) water electrolysis has gained increasing attention due to its potential to achieve low-cost, high performance hydrogen production. However, most existing membranes are not durable in industrial settings. In this work, we demonstrate good performance relative to industrial parameters using Sustainion® anion exchange membranes. Long-duration tests showed stable performance of 1 A/cm 2 at 1.85 V with a degradation rate of less than 1 μV/h over 10,000 h. The projected lifetime is thus over 20 years. Daily on/off cycling performance over the course of 30 years was simulated experimentally through accelerated voltage shock tests, resulting in a performance loss of only 0.15 μV/cycle over 11,000 cycles. As shown through impact and crossover testing, an improvement in performance is achieved by the addition of zirconia to the polymer matrix and mechanically reinforcing the membrane.

08 HYDROGEN↗
DOE OSTI2015

Meso-scale modeling of irradiated concrete in test reactor

In this paper, we detail a numerical model accounting for the effects of neutron irradiation on concrete at the mesoscale. Irradiation experiments in test reactor (Elleuch et al.,1972), i.e., in accelerated conditions, are simulated. Concrete is considered as a two-phase material made of elastic inclusions (aggregate) subjected to thermal and irradiation-induced swelling and embedded in a cementitious matrix subjected to shrinkage and thermal expansion. The role of the hardened cement paste in the post-peak regime (brittle-ductile transition with decreasing loading rate), and creep effects are investigated. Radiation-induced volumetric expansion (RIVE) of the aggregate cause the development and propagation of damage around the aggregate which further develops in bridging cracks across the hardened cement paste between the individual aggregate particles. The development of damage is aggravated when shrinkage occurs simultaneously with RIVE during the irradiation experiment. The post-irradiation expansion derived from the simulation is well correlated with the experimental data and, the obtained damage levels are fully consistent with previous estimations based on a micromechanical interpretation of the experimental post-irradiation elastic properties (Le Pape et al.,2015). In conclusion, the proposed modeling opens new perspectives for the interpretation of test reactor experiments in regards to the actual operation of light water reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

A pathway to synthesizing single-crystal Fe and FeCr films

Nuclear reactor environments provide a unique scientific and engineering challenge wherein materials must tolerate prolonged exposure to concurrent irradiation, elevated temperatures, and corrosive media. However, uncontrolled variability in material composition and structure often prohibits truly single-variable experiments that can reveal basic aspects of environmental damage. Here, Magnetron sputtering is used to provide a more controlled model system for these fundamental studies, yielding reproducible single-crystal Fe and FeCr thin films containing 8 and 18 at.% Cr. Electron microscopy is used to determine the systematic correlations between growth conditions and the resulting film microstructure and surface morphology. It is found that the substrate temperature and applied radio frequency (RF) bias can be tuned to obtain consistent homogeneous and single crystal films with a minimal amount of Ar impurities from the RF bias process. Epitaxial, single-crystal Fe films are obtained on MgO substrates at 500 °C with 10 Watt (W) RF bias deposition. However, when Cr is alloyed with Fe, higher substrate temperatures (600 °C) and applied RF biases (15 W) are required to achieve a similar epitaxial single-crystal FeCr film. Accelerated molecular dynamics simulations reveal that Cr impedes surface transport, explaining the need for higher temperature and bias during the growth of the Cr-bearing films.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Dynamics of rapidly depressurized multiphase shock tubes

Rapid depressurization is a fluid phenomenon that occurs in many industrial and natural applications. Its behaviour is often complicated by the formation, propagation and interaction of waves. In this work, we perform computer simulations of the rapid depressurization of a gas–solid mixture in a shock tube. Our problem set-up mimics previously performed experiments, which have been historically used as a laboratory surrogate for volcanic eruptions. The simulations are carried out with a discontinuous Galerkin compressible fluid solver with four-way coupled Lagrangian particle tracking capabilities. The results give an unprecedented look into the complex multiphase physics at work in this problem. Different regimes have been characterized in a regime map that highlights the key observations. While the mean flow behaviour is in good agreement with experiments, the simulations show unexpected accelerations of the particle front as it expands. Additionally, a new lifting mechanism for gas bubble (void) growth inside the gas–solid mixture is detailed.

Mechanics↗
DOE OSTI2024

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

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

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

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

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2019

Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration

In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ~6–8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. Finally, the combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.

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

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

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

Evidence for percolation diffusion of cations and reordering in disordered pyrochlore from accelerated molecular dynamics

Diffusion in complex oxides is critical to ionic transport, radiation damage evolution, sintering, and aging. In complex oxides such as pyrochlores, anionic diffusion is dramatically affected by cation disorder. However, little is known about how disorder influences cation transport. Here, we report results from classical and accelerated molecular dynamics simulations of vacancy-mediated cation diffusion in Gd 2 Ti 2 O 7 pyrochlore, on the microsecond timescale. We find that diffusion is slow at low levels of disorder, while higher disorder allows for fast diffusion, which is then accompanied by antisite annihilation and reordering, and thus a slowing of cation transport. Cation diffusivity is therefore not constant, but decreases as the material reorders. We also show that fast cation diffusion is triggered by the formation of a percolation network of antisites. This is in contrast with observations from other complex oxides and disordered media models, suggesting a fundamentally different relation between disorder and mass transport.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

The role of charge in microdroplet redox chemistry

In charged water microdroplets, which occur in nature or in the lab upon ultrasonication or in electrospray processes, the thermodynamics for reactive chemistry can be dramatically altered relative to the bulk phase. Here, we provide a theoretical basis for the observation of accelerated chemistry by simulating water droplets of increasing charge imbalance to create redox agents such as hydroxyl and hydrogen radicals and solvated electrons. We compute the hydration enthalpy of OH - and H + that controls the electron transfer process, and the corresponding changes in vertical ionization energy and vertical electron affinity of the ions, to create OH • and H • reactive species. We find that at ~ 20 - 50% of the Rayleigh limit of droplet charge the hydration enthalpy of both OH - and H + have decreased by >50 kcal/mol such that electron transfer becomes thermodynamically favorable, in correspondence with the more favorable vertical electron affinity of H + and the lowered vertical ionization energy of OH - . We provide scaling arguments that show that the nanoscale calculations and conclusions extend to the experimental microdroplet length scale. The relevance of the droplet charge for chemical reactivity is illustrated for the formation of H 2 O 2 , and has clear implications for other redox reactions observed to occur with enhanced rates in microdroplets.

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

Laboratory realization of relativistic pair-plasma beams

Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black-hole and neutron-star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with electron-positron pairs. Their role in the dynamics of such environments is in many cases believed to be fundamental, but their behavior differs significantly from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies, which are rather limited. We present the first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN’s Super Proton Synchrotron (SPS) accelerator. Monte Carlo simulations agree well with the experimental data and show that the characteristic scales necessary for collective plasma behavior, such as the Debye length and the collisionless skin depth, are exceeded by the measured size of the produced pair beams. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Near-Field Corrosion Interactions between Glass and Corrosion Resistant Alloys

This study explores the corrosion interactions between model nuclear waste glass materials and corrosion resistant alloys, under conditions that simulate the near field of a nuclear waste repository. The interactions between the corrosion of stainless steel (SS) 316, alloy G30, or alloy 625, and international simple glass or soda-lime silica glass are systematically studied. The dissimilar materials were exposed in close proximity to each other in different electrolytes at 90 oC. After exposure, the glass surface exposed near metals showed different regimes of corrosion, with distinct surface morphologies and chemistries that were likely affected by the local environment created by the localized corrosion of metals. Surface and solution analyses showed that the corrosion rate of glass was enhanced by the presence of metals. Infrared spectroscopy data suggested the local build-up of stresses in the contact area of glass, which may lead to the mechanical instability of the glass alteration layer. On the other hand, the effect of glass on metal corrosion is strongly dependent on the leaching solution. In electrolytes containing abundant aggressive anions such as Cl-, glass seems to suppress the localized corrosion of SS by the precipitation of a Si-rich surface film that protects the SS substrate from solutions. However, in less aggressive electrolytes, the corrosion rate of SS was increased by the presence of glass corrosion products. Overall, our study showed that the hidden and localized damage on glass in contact with metals may enhance the release of incorporated radionuclides compared to typical uniform glass corrosion. Therefore, more attention should be paid to the corrosion interactions between metals and glass under repository conditions.

Guo, Xiaolei↗