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 1,099 records · Page 61

Thin‐Film Ferroelectrics

Abstract Over the last 30 years, the study of ferroelectric oxides has been revolutionized by the implementation of epitaxial‐thin‐film‐based studies, which have driven many advances in the understanding of ferroelectric physics and the realization of novel polar structures and functionalities. New questions have motivated the development of advanced synthesis, characterization, and simulations of epitaxial thin films and, in turn, have provided new insights and applications across the micro‐, meso‐, and macroscopic length scales. This review traces the evolution of ferroelectric thin‐film research through the early days developing understanding of the roles of size and strain on ferroelectrics to the present day, where such understanding is used to create complex hierarchical domain structures, novel polar topologies, and controlled chemical and defect profiles. The extension of epitaxial techniques, coupled with advances in high‐throughput simulations, now stands to accelerate the discovery and study of new ferroelectric materials. Coming hand‐in‐hand with these new materials is new understanding and control of ferroelectric functionalities. Today, researchers are actively working to apply these lessons in a number of applications, including novel memory and logic architectures, as well as a host of energy conversion devices.

Fernandez, Abel↗

Probing the Kinetic Origin of Varying Oxidative Stability of Ethyl- vs. Propyl-spaced Amines for Direct Air Capture

Amine-based adsorbents are promising for direct air capture of CO 2 , yet oxidative degradation remains a key unmitigated risk hindering wide-scale deployment. Borrowing wisdom from the basic auto-oxidation scheme, insights are gained into the underlying degradation mechanisms of polyamines by quantum chemical, advanced sampling simulations, adsorbent synthesis, and accelerated degradation experiments. The reaction kinetics of polyamines are contrasted with that of typical aliphatic polymers and they elucidate for the first time the critical role of aminoalkyl hydroperoxide decomposition in the oxidative degradation of amino-oligomers. The experimentally observed variation in oxidative stability of polyamines with different backbone structures is explained by the relationship between the local chemical structure and the free energy barrier of aminoalkyl hydroperoxide decomposition, suggesting that its energetics can be used as a descriptor to screen and design new polyamines with improved stability. In conclusion, the developed computational capability sheds light on radical-induced degradation chemistry of other organic functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanism of Nonequilibrium Formation of Novel Interface Structures in Explosive Welding

Novel structures are commonly observed in explosive welding products. This paper presents the nonequilibrium mechanism of the welding interface. A theoretical model was established first to deduce the motion equation of the rapid accelerated flyer plate. Numerical simulations of the welding interface were carried out with the smooth particle hydrodynamic method. Both the theoretical and numerical results were compared with experimental phenomena and showed that under the very rapid application of an impact load usually associated with high-strength plastic deformation, a significant strain rate, and temperature rise, the material volumes of the welding interface are in the viscoelastic state and prone to rotation, detrusion, and dilatation. In addition, the extreme values from the nonequilibrium development of the interface structure are a consequence of the competition between the loading and unloading stress waves during the explosive welding process.

36 MATERIALS SCIENCE↗

FlowPM: Distributed TensorFlow implementation of the FastPM cosmological N-body solver

Here, we present FlowPM, a Particle-Mesh (PM) cosmological N-body code implemented in Mesh-TensorFlow for GPU-accelerated, distributed, and differentiable simulations. We implement and validate the accuracy of a novel multi-grid scheme based on multiresolution pyramids to compute large-scale forces efficiently on distributed platforms. We explore the scaling of the simulation on large-scale supercomputers and compare it with corresponding Python based PM code, finding on an average 10x speed-up in terms of wallclock time. We also demonstrate how this novel tool can be used for efficiently solving large scale cosmological inference problems, in particular reconstruction of cosmological fields in a forward model Bayesian framework with hybrid PM and neural network forward model. We provide skeleton code for these examples and the entire code is publicly available at https://github.com/modichirag/flowpm[Formula presented].

79 ASTRONOMY AND ASTROPHYSICS↗

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗

Machine learning accelerated discrete element modeling of granular flows

Granular flows are widely encountered in many industrial processes and natural phenomena. Discrete Element Modeling (DEM) is a useful tool for understanding and troubleshooting devices, which handle granular materials. However, its applicability is significantly limited by the huge computational cost associated with detecting and computing collisions. In this research, the computation speed of DEM was accelerated by orders of magnitude using a convolutional neural network to replace the direct calculation of particle-particle and particle-boundary collisions. The MFiX software was used to generate the training and testing dataset. Additionally, a GPU accelerated TensorFlow model was used to train the neural network and test the results. The model fluctuations caused by different training steps were reduced with a multi-scale loss function. The accuracy was improved with more frames within one training step. The modeling of a rotating drum and a hopper demonstrated the accuracy and efficiency of this machine learning accelerated DEM in the simulation of granular flows.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗

Re-additization of commercial biodiesel blends during long-term storage

Commercial biodiesel blends were aged at 43 degrees C while monitoring stability. The oxidation stability -- or oxidation reserve expressed as Rancimat induction period (IP) -- gradually decreased from its initial value. At a predetermined IP threshold, an antioxidant was used to restore IP to the ASTM D7467 specification minimum of 6 h, referred to as re-additization. At lower IP values, the amount of antioxidant required increased significantly, and the effectiveness tended to be reduced. Once IP fell to essentially zero, the acid content increased to above the allowable limit and insoluble material was also detected. Storage life was increased relative to the as-received fuels as evidenced by longer time to produce acids. Experience in the field may vary based on storage conditions; however, these results indicate re-additization can significantly increase storage life of biodiesel blends when used with regular monitoring of IP and acid number. An assessment of the storage stability of the as-received fuels showed that the initial IP did not predict storage behavior, although fuels above the specification minimum remained stable for >12 weeks accelerated aging (1 year simulated).

25 ENERGY STORAGE↗

Interpolation of computed gamma-ray detector response functions

Gamma-ray spectra measured by traditional detectors contain features that result from a combination of the effects of detector materials/geometry, the incident gamma-ray energy, and the angle of entry. The features, such as the full-energy photopeak, Compton continuum, annihilation peak, and escape peaks, are governed by simple relationships depending on incident energy and have been known for a long time. Monte Carlo computer simulations of gamma rays interacting with a detector will show these features, and with a resolution function applied, the results should look similar to real measurements. The traditional approach to creating a detector response function requires many separate simulations of monoenergetic gamma rays striking the detector. This paper presents a new approach to developing computed detector response functions. The new approach involves a much smaller number of monoenergetic gamma-ray simulations and uses interpolation to quickly generate the responses of gamma rays that were not simulated. During the interpolation process, the underlying physics equations are used to accurately compute the response of a given energy gamma ray from the small set of simulations. Such work enables accelerated generation of synthetic radiation detector data.

Detector response↗

Zirconium nanoparticle coating development for FCCI diffusion barrier in nuclear cladding

Fuel clad chemical interaction (FCCI) is a pervasive issue for metallic nuclear fuels. FCCI can damage cladding, reduce thermal efficiency, and compromise the structural integrity of the fuel rod. In this study, we evaluated the FCCI mitigation capabilities of electrophoretically-deposited Zr nanoparticle coatings on ferritic-martensitic stainless-steel cladding (HT9). Zr was deposited on the internal surface of 6 mm outer diameter cladding using electrophoretic deposition (EPD). Thickness varied among samples with similar EPD conditions. A uniformly coated sample was selected and sent to the Argonne Tandem Linac Accelerator System (ATLAS) to simulate the coating's behavior under irradiation conditions. Furthermore, the coating prevented FCCI in diffusion couple studies. Results of this study confirm that Zr could potentially make an effective FCCI barrier in metallic fuel rods but the deposition process needs improvement.

36 MATERIALS SCIENCE↗

Assessment of Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

We report the availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H 2 ) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 data set consisting of equilibrium geometries and assessed the performance of 55 density functionals over this data set (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput. 2020, 16, 4963-4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we extended this data set with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H 2 binding, forming the H2Bind78 × 7 data set. By assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78 × 7 data set, we identified the best performing density functionals for H 2 binding applications: PBE0-DH, ωB97X-V, ωB97M-V, and DSD-PBEPBE-D3(BJ). The addition of Hartree-Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of ωB97X-V and ωB97M-V density functionals as they offer good performance for both geometries and energies. In addition, we also identified B97M-V and B97M-rV as the best semilocal density functionals for predicting H 2 binding energy at its equilibrium geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Dimensional Atomistic Neural Network Potentials for Molecule–Surface Interactions: HCl Scattering from Au(111)

Ab initio molecular dynamics (AIMD) simulations of molecule–surface scattering allow first-principles characterization of the dynamics. However, the large number of density functional theory calculations along the trajectories is very costly, limiting simulations of long-time events and giving rise to poor statistics. To avoid this computational bottleneck, we report here the development of a high-dimensional molecule–surface interaction potential energy surface (PES) with movable surface atoms, using a machine learning approach. With 60 degrees of freedom, this PES allows energy transfer between the energetic impinging molecule and thermal surface atoms. Classical trajectory calculations for the scattering of DCl from Au(111) on this PES are found to agree well with AIMD simulations, with ~105-fold acceleration. Scattering of HCl from Au(111) is further investigated and compared with available experimental results.

Kolb, Brian↗

Role of Greenland Freshwater Anomaly in the Recent Freshening of the Subpolar North Atlantic

The cumulative Greenland freshwater flux anomaly has exceeded 5,000 km 3 since the 1990s. The volume of this surplus freshwater is expected to cause substantial freshening in the North Atlantic. Analysis of hydrographic observations in the subpolar seas reveals freshening signals in the 2010s. The sources of this freshening are yet to be determined. In this study, the relationship between the surplus Greenland freshwater flux and this freshening is tested by analyzing the propagation of the Greenland freshwater anomaly and its impact on salinity in the subpolar North Atlantic based on observational data and numerical experiments with and without the Greenland runoff. A passive tracer is continuously released during the simulations at freshwater sources along the coast of Greenland to track the Greenland freshwater anomaly. Tracer budget analysis shows that 44% of the volume of the Greenland freshwater anomaly is retained in the subpolar North Atlantic by the end of the simulation. This volume is sufficient to cause strong freshening in the subpolar seas if it stays in the upper 50–100 m. However, in the model the anomaly is mixed down to several hundred meters of the water column resulting in smaller magnitudes of freshening compared to the observations. Therefore, the simulations suggest that the accelerated Greenland melting would not be sufficient to cause the observed freshening in the subpolar seas and other sources of freshwater have contributed to the freshening. Impacts on salinity in the subpolar seas of the freshwater transport through Fram Strait and precipitation are discussed.

54 ENVIRONMENTAL SCIENCES↗

Deep Learning Emulation of Subgrid-Scale Processes in Turbulent Shear Flows

Deep neural networks (DNNs) are developed from a data set obtained from the dynamic Smagorinsky model to emulate the subgrid-scale (SGS) viscosity (ν sgs ) and diffusivity (κ sgs ) for turbulent stratified shear flows encountered in the oceans and the atmosphere. These DNNs predict ν sgs and κ sgs from velocities, strain rates, and density gradients such that the evolution of the kinetic energy budget and density variance budget terms is similar to the corresponding values obtained from the original dynamic Smagorinsky model. These DNNs also compute ν sgs and κ sgs ~2–4 times quicker than the dynamic Smagorinsky model resulting in a ~2–2.5 times acceleration of the entire simulation. This study demonstrates the feasibility of deep learning in emulating the subgrid-scale (SGS) phenomenon in geophysical flows accurately in a cost-effective manner. In a broader perspective, deep learning-based surrogate models can present a promising alternative to the traditional parameterizations of the subgrid-scale processes in climate models.

58 GEOSCIENCES↗

Effects of presenilin-1 familial Alzheimer’s disease mutations on γ-secretase activation for cleavage of amyloid precursor protein

Presenilin-1 (PS1) is the catalytic subunit of γ-secretase which cleaves within the transmembrane domain of over 150 peptide substrates. Dominant missense mutations in PS1 cause early-onset familial Alzheimer’s disease (FAD); however, the exact pathogenic mechanism remains unknown. Here we combined Gaussian accelerated molecular dynamics (GaMD) simulations and biochemical experiments to determine the effects of six representative PS1 FAD mutations (P117L, I143T, L166P, G384A, L435F, and L286V) on the enzyme-substrate interactions between γ-secretase and amyloid precursor protein (APP). Biochemical experiments showed that all six PS1 FAD mutations rendered γ-secretase less active for the endoproteolytic (ε) cleavage of APP. Distinct low-energy conformational states were identified from the free energy profiles of wildtype and PS1 FAD-mutant γ-secretase. The P117L and L286V FAD mutants could still sample the “Active” state for substrate cleavage, but with noticeably reduced conformational space compared with the wildtype. The other mutants hardly visited the “Active” state. The PS1 FAD mutants were found to reduce γ-secretase proteolytic activity by hindering APP residue L49 from proper orientation in the active site and/or disrupting the distance between the catalytic aspartates. Therefore, our findings provide mechanistic insights into how PS1 FAD mutations affect structural dynamics and enzyme-substrate interactions of γ-secretase and APP.

60 APPLIED LIFE SCIENCES↗

Three-dimensional magnetic optimization of accelerator magnets using an analytic strip model

The end design is a critical step in the design of superconducting accelerator magnets. First, the strain energy of the conductors must be minimized, which can be achieved using differential geometry. The end design also requires an optimization of the magnetic field homogeneity. A mechanical and magnetic model for the conductors, using developable strips, is described in this paper. This model can be applied to superconducting Rutherford cables, and it is particularly suitable for High Temperature Superconducting tapes. The great advantage of this approach is analytic simplifications in the field computation, allowing for very fast and accurate computations, which save a considerable computational time during the optimization process. Some 3D designs for dipoles are finally proposed, and it is shown that the harmonic integrals can be easily optimized using this model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An accurate optical design method for synchrotron radiation beamlines with wave-front aberration theory

An aberration calculation method which was developed by Lu [1] can treat individual aberration term precisely. Spectral aberration is the linear sum of these aberration terms, and the aberrations of multi-element systems also can be calculated correctly when the stretching ratio, defined herein, is unity. Evaluation of focusing mirror-grating systems which are optimized according to Lu’s method, along with the Light Path Function (LPF) and the Spot Diagram method (SD) are discussed to confirm the advantage of Lu’s methodology. Lu’s aberration terms are derived from a precise wave-front treatment, whereas the terms of the power series expansion of the light path function do not yield an accurate sum of the aberrations. Moreover, Lu’s aberration terms can be individually optimized. This is not possible with the analytical spot diagram formulae.

43 PARTICLE ACCELERATORS↗

Dephasing of ion beams as magnetic vortex acceleration regime transitions into a bubble-like field structure

The interaction of an ultra-intense laser pulse with a near critical density target results in the formation of a plasma channel, a strong azimuthal magnetic field and moving vortices. An application of this is the generation of energetic and collimated ion beams via magnetic vortex acceleration. The optimized regime of magnetic vortex acceleration is becoming experimentally accessible with new high intensity laser beamlines coming online and advances made in near critical density target fabrication. The robustness of the acceleration mechanism with realistic experimental conditions is examined with three-dimensional simulations. Of particular interest is the acceleration performance with different laser temporal contrast conditions, in some cases leading to pre-expanded target profiles prior to the arrival of the main pulse. Preplasma effects on the structure of the accelerating fields are explored, including a detailed analysis of the ion beam properties and the efficiency of the process. Improved scaling laws for the magnetic vortex acceleration mechanism, including the laser focal spot size effects, are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗