Search NASA⌕ Search

SEARCH · Search NASA

Results for “PARTICLE INTERACTIONS”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Essential barrier height and a probabilistic approach in characterizing potential landscape

In this work we propose a probabilistic approach to investigate the shape of landscapes of multi-dimensional potential functions. Under a suitable coupling scheme, two copies of the overdamped Langevin dynamics associated with the potential function are coupled, and the coupling times are collected. Assuming a set of intuitive yet technically challenging conditions on the coupling scheme, it is shown that the tail distributions of the coupling times exhibit qualitatively different dependencies on the noise magnitude for single-well versus multi-well potential functions. More specifically, for convex single-well potentials, the negative tail exponent of the coupling time distribution is uniformly bounded away from zero by the convexity parameter and is independent of the noise magnitude. In contrast, for multi-well potentials, the negative tail exponent decreases exponentially as the noise vanishes, with the decay rate governed by the essential barrier height, a quantity introduced in this paper to characterize the non-convex nature of the potential function. Numerical investigations are conducted for a variety of examples, including the Rosenbrock function, interacting particle systems, and loss functions arising in artificial neural networks. These examples not only illustrate the theoretical results in various contexts but also provide crucial numerical validation of the conjectured assumptions, which are essential to the theoretical analysis yet lie beyond the reach of standard technical tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Solvent and A-Site Cation Control Preferred Crystallographic Orientation in Bromine-Based Perovskite Thin Films

Preferred crystallographic orientation in polycrystalline films is desirable for efficient charge carrier transport in metal halide perovskites and semiconductors. However, the mechanisms that determine the preferred orientation of halide perovskites are still not well understood. In this work, we investigate crystallographic orientation in lead bromide perovskites. We show that the solvent of the precursor solution and organic A-site cation strongly affect the preferred orientation of the deposited perovskite thin films. Specifically, we show that the solvent, dimethylsulfoxide, influences the early stages of crystallization and induces preferred orientation in the deposited films by preventing colloidal particle interactions. Additionally, the methylammonium A-site cation induces a higher degree of preferred orientation than the formamidinium counterpart. We use density functional theory to show that the lower surface energy of the (100) plane facets in methylammonium-based perovskites, compared to the (110) planes, is the reason for the higher degree of preferred orientation. In contrast, the surface energy of the (100) and (110) facets is similar for formamidinium-based perovskites, leading to lower degree of preferred orientation. Furthermore, we show that different A-site cations do not significantly affect ion diffusion in bromine-based perovskite solar cells but impact ion density and accumulation, leading to increased hysteresis. Our work highlights the interplay between the solvent and organic A-site cation which determine crystallographic orientation and plays a critical role in the electronic properties and ionic migration of solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Including Physics-Informed Atomization Constraints in Neural Networks for Reactive Chemistry

Machine learning interatomic potentials (MLIPs) have emerged as powerful tools for investigating atomistic systems with high accuracy and a relatively low computational cost. However, a common and unaddressed challenge with many current neural network (NN) MLIP models is their limited ability to accurately predict the relative energies of systems containing isolated or nearly isolated atoms, which appear in various reactive processes. To address this limitation, we present a mathematical technique for modifying any existing atom-centered NN architecture to account for the energies of isolated atoms. The result produces a consistent prediction of the atomization energy (AE) of a system using minimal constraints on the model. Using this technique, we build a model architecture that we call hierarchically interacting particle neural network (HIP-NN)-AE, an AE-constrained version of the HIP-NN, as well as ANI-AE, the AE-constrained version of the accurate NN engine for molecular energies (ANI). Our results demonstrate AE consistency of AE-constrained models, which drastically improves the AE predictions for the models. We compare the AE-constrained approach to unconstrained models as well as models from the literature in other scenarios, such as bond dissociation energies, bond dissociation pathways, and extensibility tests. These results show that the constraints improve the model performance in some of these tasks and do not negatively affect the performance on any tasks. The AE constraint approach thus offers a robust solution to the challenges posed by isolated atoms in energy prediction tasks.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗
DOE OSTI2018

Transferable Dynamic Molecular Charge Assignment Using Deep Neural Networks

We report that the ability to accurately and efficiently compute quantum-mechanical partial atomistic charges has many practical applications, such as calculations of IR spectra, analysis of chemical bonding, and classical force field parametrization. Machine learning (ML) techniques provide a possible avenue for the efficient prediction of atomic partial charges. Modern ML advances in the prediction of molecular energies [i.e., the hierarchical interacting particle neural network (HIP-NN)] has provided the necessary model framework and architecture to predict transferable, extensible, and conformationally dynamic atomic partial charges based on reference density functional theory (DFT) simulations. Utilizing HIP-NN, we show that ML charge prediction can be highly accurate over a wide range of molecules (both small and large) across a variety of charge partitioning schemes such as the Hirshfeld, CM5, MSK, and NBO methods. To demonstrate transferability and size extensibility, we compare ML results with reference DFT calculations on the COMP6 benchmark, achieving errors of 0.004e – (elementary charge). This is remarkable since this benchmark contains two proteins that are multiple times larger than the largest molecules in the training set. An application of our atomic charge predictions on nonequilibrium geometries is the generation of IR spectra for organic molecules from dynamical trajectories on a variety of organic molecules, which show good agreement with calculated IR spectra with reference method. Critically, HIP-NN charge predictions are many orders of magnitude faster than direct DFT calculations. Lastly, these combined results provide further evidence that ML (specifically HIP-NN) provides a pathway to greatly increase the range of feasible simulations while retaining quantum-level accuracy.

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

Simulating Crystallization in a Colloidal System Using State Predictive Information Bottleneck Based Enhanced Sampling

Here, we investigate crystal nucleation in supersaturated colloid suspensions using enhanced molecular dynamics simulations augmented with machine learning techniques. The simulations reveal that crystallization in the model colloidal system studied here, with particles interacting through a repulsive screened Coulomb Yukawa potential, proceeds from vapor to dense liquid droplet to crystalline phases across multiple high barriers. Employing a one-dimensional reaction coordinate derived from the State Predictive Information Bottleneck framework, our simulations capture back-and-forth phase transitions across multiple barriers effectively in biased metadynamics simulations. We obtain relative free energy differences between different phases and also quantify the roles of different molecular level features in driving the phase changes.

Chemistry↗
DOE OSTI2018

Electrolyte-Assisted Hydrogen Storage Reactions

Use of electrolytes, in the form of LiBH 4 /KBH 4 and LiI/KI/CsI eutectics, is shown to significantly improve (by more than a factor of 10) both the dehydrogenation and full rehydrogenation of the MgH 2 /Sn destabilized hydride system and the hydrogenation of MgB 2 to Mg(BH 4 ) 2 . The improvement revealed that interparticle transport of atoms heavier than hydrogen can be an important rate-limiting step during hydrogen cycling in hydrogen storage materials consisting of multiple phases in powder form. Electrolytes enable solubilizing heavy ions into a liquid environment and thereby facilitate the reaction over full surface areas of interacting particles. In conclusion, the examples presented suggest that use of electrolytes in the form of eutectics, ionic liquids, or solvents containing dissolved salts may be generally applicable for increasing reaction rates in complex and destabilized hydride materials.

08 HYDROGEN↗
DOE OSTI2018

Energization of the Ring Current by Substorms

The substorm process releases large amounts of energy into the magnetospheric system, although where the energy is transferred to and how it is partitioned remains an open question. In this study, we address whether the substorm process contributes a significant amount of energy to the ring current. The ring current is a highly variable region, and understanding the energization processes provides valuable insight into how substorm-ring current coupling may contribute to the generation of storm conditions and provide a source of energy for wave driving. In order to quantify the energy input into the ring current during the substorm process, we analyze Radiation Belt Storm Probes Ion Composition Experiment and Helium Oxygen Proton Electron ion flux measurements for H + , O + , and He + . The energy content of the ring current is estimated and binned spatially for L and magnetic local time. The results are combined with an independently derived substorm event list to perform a statistical analysis of variations in the ring current energy content with substorm phase. We show that the ring current energy is significantly higher in the expansion phase compared to the growth phase, with the energy enhancement persisting into the substorm recovery phase. The characteristics of the energy enhancement suggest the injection of energized ions from the tail plasma sheet following substorm onset. The local time variations indicate a loss of energetic H + ions in the afternoon sector, likely due to wave-particle interactions. Overall, we find that the average energy input into the ring current is ~9% of the previously reported energy released during substorms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Pitch Angle Dependence of Electron and Ion Flux Changes During Local Magnetic Dipolarization Inside Geosynchronous Orbit

In this study we have statistically examined flux changes of 1–1,000-keV electrons, and hydrogen, helium, and oxygen ions during local magnetic dipolarization inside geosynchronous orbit, with a focus on their pitch angle dependence. Here, using 144 dipolarization events that were selected in a previous study with Van Allen Probes observations in 2012–2016, we have performed a superposed epoch analysis of differential flux changes after the dipolarization onset for each species. On average the electron flux increases primarily around pitch angle (α) = 90° at >80 keV and almost isotropically at 10–50 keV. The electron flux at <5 keV increases at α = 0° and 180° but slightly decreases (or remains unchanged) at α = 90°. On the other hand, the ion flux at >80 keV increases around α = 90°, while at <30 keV it decreases nearly independent of pitch angle. Only the low-energy (<5 keV) helium and oxygen ion flux changes indicate a strong field-aligned enhancement, which could be attributed to their outflows from the topside ionosphere. After the dipolarization onset, >5-keV electron and >30-keV ion fluxes exhibit a perpendicular anisotropy, which is most pronounced for 50–200-keV electrons. Both distributions become more isotropic with decreasing energy. A noticeable field-aligned anisotropy is seen for <5-keV ion fluxes. These statistical pitch angle-dependent electron and ion properties during dipolarizations may be explained by combining various processes, including adiabatic acceleration and/or transport, wave-particle interactions, and ion outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electron Mixing and Isotropization in the Exhaust of Asymmetric Magnetic Reconnection With a Guide Field

In this work, we investigate an exhaust crossing of asymmetric guide field reconnection observed by Magnetospheric Multiscale (MMS) mission at Earth's dayside magnetopause. One MMS spacecraft (MMS 4) observes multicomponent electron distributions, including two counterstreaming electron beams, while the other three MMS spacecraft, with a separation of ~30 km, record nearly isotropic electron distributions. As counterstreaming electrons are unstable for the electron two-stream instability, our observations suggest that the electrostatic waves generated by the fast-growing electron two-stream instability can contribute to the rapid isotropization of electron distributions in the reconnection exhaust, indicating that wave-particle interactions play an important role in electron dynamics.

58 GEOSCIENCES↗
DOE OSTI2020

The Effects of Field Line Curvature (FLC) Scattering on Ring Current Dynamics and Isotropic Boundary

Abstract In the ring current dynamics, various loss mechanisms contribute to the ring current decay, including losses to the upper atmosphere through particle precipitation. This study implements the field line curvature (FLC) scattering mechanism in a kinetic ring current model and investigates its role in precipitating ions into the ionosphere during the 17 March 2013 storm. Simulation results indicate that (1) the FLC scattering process exerts on energetic ions on the nightside where the magnetospheric configuration is more stretching. It is more effective on heavy ions (e.g., O + ). These ion losses thereafter lead to a faster recovery of the ring current. (2) The FLC‐associated ion precipitation mainly occurs in the outer region (L > 5 for protons and L > 4.5 for oxygen ions) on the nightside. The O + precipitation takes places in a wider region than protons although its intensity is much lower. Comparisons with POES observations suggest that more proton precipitation is needed in the inner region. This is probably caused by the less stretched configuration in the simulation that prevents more precipitation. It may also imply that other loss process is required in the model such as wave‐particle interactions. (3) The storm time precipitating proton flux of tens of keV due to the FLC scattering sometimes becomes comparable to that of electrons on the nightside, although electrons usually dominate the ionospheric energy deposition from the midnight eastward toward the dayside. (4) The FLC scattering process seems to be capable of explaining the formation of isotropic boundary in the ionosphere during the investigated event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Lightning Rings and Gravity Waves: Insights Into the Giant Eruption Plume From Tonga's Hunga Volcano on 15 January 2022

Abstract On 15 January 2022, Hunga Volcano in Tonga produced the most violent eruption in the modern satellite era, sending a water‐rich plume at least 58 km high. Using a combination of satellite‐ and ground‐based sensors, we investigate the astonishing rate of volcanic lightning (>2,600 flashes min −1 ) and what it reveals about the dynamics of the submarine eruption. In map view, lightning locations form radially expanding rings. We show that the initial lightning ring is co‐located with an internal gravity wave traveling >80 m s −1 in the stratospheric umbrella cloud. Buoyant oscillations of the plume's overshooting top generated the gravity waves, which enhanced turbulent particle interactions and triggered high‐current electrical discharges at unusually high altitudes. Our analysis attributes the intense lightning activity to an exceptional mass eruption rate (>5 × 10 9 kg s −1 ), rapidly expanding umbrella cloud, and entrainment of abundant seawater vaporized from magma‐water interaction at the submarine vent.

58 GEOSCIENCES↗
DOE OSTI2023

Outer Radiation Belt Flux and Phase Space Density Response to Sheath Regions: Van Allen Probes and GPS Observations

Turbulent and compressed sheath regions preceding interplanetary coronal mass ejections strongly impact electron dynamics in the outer radiation belt. Changes in electron flux can occur on timescales of tens of minutes, which are unlikely to be captured by a two-satellite mission. The recently released Global Positioning System (GPS) data set generally has shorter revisit times (at L ~ 4–8) owing to the large number of satellites in the constellation equipped with energetic particle detectors. Investigating electron fluxes at energies from 140 keV to 4 MeV and sheaths observed in 2012–2018, we show that the flux response to sheaths on a timescale of 6 hr, previously reported from Van Allen Probes (RBSP) data, is reproduced by GPS measurements. Furthermore, GPS data enables derivation of the response on a timescale of 30 min, which further confirms that the energy and L-shell dependent changes in electron flux are associated with the impact of the sheath. Sheath-driven loss is underestimated over longer timescales as the electrons recover during the ejecta. We additionally show the response of electron phase space density (PSD), which is a key quantity in identifying non-adiabatic loss from the system and electron energization through wave-particle interactions. The PSD response is calculated from both RBSP and GPS data for the 6 hr timescale, as well as from GPS data for the 30 min timescale. The response is divided based on the geoeffectiveness of the sheaths revealing that electrons are effectively accelerated only during geoeffective sheaths, while loss commonly occurs during all sheaths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Global Distribution of EMIC Waves and Its Association to Subauroral Proton Precipitation During the 27 May 2017 Storm: Modeling and Multipoint Observations

Recent simulation studies using the RAM-SCB model showed that proton precipitation contributes significantly to the total energy flux deposited into the subauroral ionosphere thereby affecting the magnetosphere-ionosphere coupling. Here, in this study, we use the BATS-R-US + RAM-SCB model to understand the evolution of ElectroMagnetic Ion Cyclotron (EMIC) waves in the inner magnetosphere, their correspondence to the proton precipitation into the subauroral ionosphere, and to assess the performance of the model in reproducing the EMIC wave-particle interactions. During the 27 May 2017 storm, Arase and RBSP-A satellites observed typical signatures of EMIC waves in the inner magnetosphere. Within this interval, Defense Meteorological Satellite Program (DMSP) and National Oceanic and Atmospheric Administration (NOAA)/MetOp satellites observed significant proton precipitation in the dusk-midnight sector. Simulation results show that H- and He-band EMIC waves are excited within regions of strong temperature anisotropy near the plasmapause. The simulated growth rates of EMIC waves show a similar trend to that of the EMIC wave power observed by the Arase and RBSP-A satellites, suggesting that the model can reproduce the EMIC wave activity qualitatively. The simulated H-band waves in the dusk sector are stronger than He-band waves possibly due to the presence of excess protons in the boundary conditions obtained from the BATS-R-US code. The precipitating proton fluxes reproduced by the simulation with EMIC waves are found to agree reasonably well with the DMSP and NOAA/MetOp satellite observations. It is suggested that EMIC wave scattering of ring current ions can account for proton precipitation observed by the DMSP and MetOp satellites during the 27 May 2017 storm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Giant Undulations Driven by Pitch‐Angle Scattering of Time Domain Structures Modulated by Plasmapause Surface Wave

Abstract Plasmapause surface waves (PSWs) near the plasmapause boundary are regarded to be the magnetospheric source of ionospheric auroral giant undulations (GUs) located at the equatorward boundary of diffuse aurora. However, the observational evidence of wave‐particle interaction connecting PSWs and GUs is absent. In this letter, we demonstrate GUs are driven by pitch‐angle scattering of time domain structures modulated by the PSWs, based on the conjugated ionospheric and magnetospheric observations. Specifically, ionospheric GUs are lighted by the pitch‐angle scattering of <1 keV thermal electron and ions and energetic ions with energy up to dozens of keV near the plasmapause. Further, the total fluxes during one PSW period and energy of scattered electron and ions determine the size and luminosity of GUs. Our research provides observational evidence that PSWs cause periodic electron precipitation via modulating the time domain structures rather than the previously predicted chorus or electron cyclotron harmonic waves.

Zhou, Yi‐Jia [Weihai Institute for Interdisciplina↗
DOE OSTI2018

Switchable geometric frustration in an artificial-spin-ice–superconductor heterosystem

Geometric frustration emerges when local interaction energies in an ordered lattice structure cannot be simultaneously minimized, resulting in a large number of degenerate states. The numerous degenerate configurations may lead to practical applications in microelectronics(1), such as data storage, memory and logic(2). However, it is difficult to achieve very high degeneracy, especially in a two-dimensional system(3,4). Here, we showcase in situ controllable geometric frustration with high degeneracy in a two-dimensional flux-quantum system. We create this in a superconducting thin film placed underneath a reconfigurable artificial-spin-ice structure(5). The tunable magnetic charges in the artificial-spin-ice strongly interact with the flux quanta in the superconductor, enabling switching between frustrated and crystallized flux quanta states. The different states have measurable effects on the superconducting critical current profile, which can be reconfigured by precise selection of the spin-ice magnetic state through the application of an external magnetic field. We demonstrate the applicability of these effects by realizing a reprogrammable flux quanta diode. The tailoring of the energy landscape of interacting 'particles' using artificial-spin-ices provides a new paradigm for the design of geometric frustration, which could illuminate a path to control new functionalities in other material systems, such as magnetic skyrmions(6), electrons and holes in two-dimensional materials(7,8), and topological insulators(9), as well as colloids in soft materials(10-13).

Wang, Yong-Lei (ORCID:0000000303917757)↗
DOE OSTI2018

Membrane-induced hydroelastic migration of a particle surfing its own wave

While coupling between fluid flow and soft elastic surfaces is common in biology and engineering, an analytical description is challenging as it often involves non-linear dynamics. Here we show using theory and experiments that a small particle moving along an elastic membrane through a viscous fluid is repelled from the membrane due to hydroelastic forces. The flow field produces an elastic disturbance in the membrane leading to particle–wave coupling. We derive an analytic expression for the particle trajectory and find that the normal migration velocity of the particle is quadratic in its speed and depends on a combination of the tension and bending resistance of the membrane. Experimentally, we measure the normal displacement of spheres sedimenting under gravity along a suspended elastic membrane and find quantitative agreement with the theoretical predictions with no fitting parameters. We experimentally demonstrate that the effect is strong enough for separation and sorting of particles on the basis of both their size and density. We discuss the significance of our results for particles interacting with biological membranes, and propose the use of our model for membrane elasticity measurements.

catalysis (homogeneous), solar (photovoltaic), bio↗