SEARCH · Search NASA

Results for “ANALOG 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Real-time simulation of asymmetry generation in fermion-bubble collisions

Motivated by the out-of-equilibrium dynamics during an early-Universe first-order phase transition, we perform real-time simulations of fermion-bubble scattering in 1 + 1 dimensions. This nonequilibrium process can generate a charge-conjugation C asymmetry outside the bubble wall, induced by the complex fermion mass profile. The resulting C asymmetry is the 1 + 1 -dimensional analog of the C P asymmetry in 3 + 1 dimensions, a key ingredient in baryon asymmetry generation at the electroweak scale. Using tensor network methods, we track the real-time evolution of the C asymmetry in the charge density as the fermion interacts with the bubble wall, a regime inaccessible to analytic calculations. We further introduce two observables to quantify the asymmetry in the asymptotic region where reflected particles are well separated from the scattering point: one based on the net charge outside the bubble wall, and the other on the spatial displacement between the reflected particle and antiparticle wave packets. Our study represents a first step toward nonperturbative, real-time computations of C P asymmetry in 3 + 1 dimensions for electroweak baryogenesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Emulating 2D Materials with Magnons

Spin waves (magnons) in two-dimensional (2D) materials have received increasing interest due to their unique states and potential for tunability. However, many interesting features of these systems, including Dirac points and topological states, occur at high frequencies, where experimental probes are limited. Here, we study a crystal formed by patterning a hexagonal array of holes in a perpendicularly magnetized thin film. Through simulation, we find that the magnonic band structure imitates that of graphene, but additionally has some kagomelike character and includes a few flat bands. Surprisingly, its nature can be understood using a nine-band tight-binding Hamiltonian. This clear analogy to 2D materials enables band-gap engineering in 2D, topological magnons along 1D phase boundaries, and spectrally isolated modes at 0D point defects. Interestingly, the 1D phase boundaries allow access to the valley degree of freedom through a magnonic analog of the quantum valley Hall insulator. These approaches can be extended to other magnonic systems, but are potentially more general due to the simplicity of the model, which resembles existing results from electron, phonon, photon, and cold-atom systems. This finding brings the physics of spin waves in 2D materials to more experimentally accessible scales, augments it, and outlines a few principles for controlling magnonic states.

Ferrimagnets↗

DOE OSTI2026

Wasserstein Normalized Autoencoder for Anomaly Detection in ProtoDUNE Vertical-Drift Detector

ProtoDUNE Vertical Drift needs a selective triggering algorithm. The detector sits on Earth's surface, so cosmic activity dominates its data. Our goal in this paper is to trigger on neutrino events more robustly than the current deployed Analog-to-Digital Converter Simple Window (ADCSW) model and, eventually, search for signals of Beyond Standard Model (BSM) physics at DUNE as our ultimate North Star objective. As a step towards this goal, we evaluate a Wasserstein Normalized Autoencoder (WNAE) on simulated collection-plane only windows of shape $1\times10\times10$ where Neutrinos act as our BSM-proxy and Cosmic-ray Muons serve as our learned background. The network parameters are fitted using only cosmic-ray muon events as background in order to maintain an unsupervised pipeline. Training uses finite-step Langevin $x^-$ samples, positive-sample reconstruction energy, and an empirical sliced $2$-Wasserstein objective to learn a normalized Boltzmann energy model. We then calibrate on a nominal $5\,\mathrm{Hz}$ operating threshold calculated from cosmic validation data. Both WNAE and ADCSW accept 311 of 194,083 held-out cosmic background events at this $5\,\mathrm{Hz}$ threshold. We found that WNAE accepts 9,677 of 34,634 neutrino-proxy events $(27.9\pm0.24)\%$, compared with 10,076 $(29.1\pm0.24)\%$ for ADCSW, an observed WNAE-minus-ADCSW difference of $-1.15\%$. At another nominal $2\,\mathrm{Hz}$ target threshold, the corresponding efficiencies are $(20.5\pm0.22)\%$ and $(22.6\pm0.22)\%$, respectively. Of the WNAE-selected neutrino proxies at $5\,\mathrm{Hz}$, $(20.8\pm0.4)\%$ of the classified neutrino-proxy events are unique to WNAE, where the uncertainty is an absolute binomial standard error of $0.4\%$.

Zheng, Jake [U. Chicago (main)] (ORCID:00090002189↗

DOE OSTI2026

Wasserstein Normalized Autoencoder for Anomaly Detection in ProtoDUNE Vertical-Drift Detector

Zheng, Jake [Chicago U.] (ORCID:0009000218901379)↗

DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Sensitivity of nEXO to 136 Xe charged-current interactions: Background-free searches for solar neutrinos and fermionic dark matter

We study the sensitivity of nEXO to solar neutrino charged-current interactions, 𝜈 𝑒 + 136 Xe → 136 Cs*+𝑒 − , as well as analogous interactions predicted by models of fermionic dark matter. Due to the recently observed low-lying isomeric states of 136 Cs, these interactions will create a time-delayed coincident signal observable in the scintillation channel. Here we develop a detailed Monte Carlo simulation of scintillation emission, propagation, and detection in the nEXO detector to model these signals under different assumptions about the timing resolution of the photosensor readout. We show this correlated signal can be used to achieve background discrimination on the order of 10 −9 , enabling nEXO to make background-free measurements of solar neutrinos above the reaction threshold of 0.668 MeV. We project that nEXO could measure the flux of neutrinos from the carbon-nitrogen-oxygen cycle with a statistical uncertainty of 25%, thus contributing a novel and competitive measurement toward addressing the solar metallicity problem. Additionally, nEXO could measure the mean energy of the 7 Be neutrinos with a precision of 𝜎 ≤ 1.5 keV and could determine the survival probability of 7 Be and pep solar 𝜈𝑒 with precision comparable to the state of the art. These quantities are sensitive to the Sun’s core temperature and to nonstandard neutrino interactions, respectively. Here, the strong background suppression would allow nEXO to search for charged-current interactions of fermionic dark matter in the mass range 𝑚 𝜒 = 0.668–7 MeV with a sensitivity up to three orders of magnitude better than current limits.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Wholesale Electricity Analysis via Simulation & Learning Experiments (WEASLE): Platform Development and Pilot Competition

This document reports the development of the Wholesale Electricity Analysis via Simulation and Learning Experiments (WEASLE) platform and the pilot competition that was conducted to test the platform. Due to the increasing reliance on variable renewable energy resources for bulk power, the pilot competition, called the Energy Storage Participation Algorithm Competition (ESPA-Comp), was used to test the effect of various market designs on storage utilization and market efficiency. Basic details of the platform are provided, including an overview of the market clearing engine, the battery dispatch and degradation models, electric grid topology and resource mix, and software architecture. Two market designs were tested: a two-settlement market analogous to typical ISO design today, and a multi-settlement market that allows additional forward-trading periods during the real-time market. Results from the pilot competition show that the storage bidding problem is nontrivial and is well suited for future challenges. We find that: 1) all four teams utilized different approaches to the bidding problem, 2) different methodological approaches led to substantially different offer behaviors, 3) resource profits are clustered by team and methodological approach, 4) simulated offers reduced market surplus by about 0.5%, 5) substantially different prices between two-settlement and multi-settlement markets albeit minimal difference in overall market surplus.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Toward a quantitative description of solvation structure: a framework for differential solution scattering measurements

Appreciating that the role of the solute–solvent and other outer-sphere interactions is essential for understanding chemistry and chemical dynamics in solution, experimental approaches are needed to address the structural consequences of these interactions, complementing condensed-matter simulations and coarse-grained theories. High-energy X-ray scattering (HEXS) combined with pair distribution function analysis presents the opportunity to probe these structures directly and to develop quantitative, atomistic models of molecular systems in situ in the solution phase. However, at concentrations relevant to solution-phase chemistry, the total scattering signal is dominated by the bulk solvent, prompting researchers to adopt a differential approach to eliminate this unwanted background. Though similar approaches are well established in quantitative structural studies of macromolecules in solution by small- and wide-angle X-ray scattering (SAXS/WAXS), analogous studies in the HEXS regime—where sub-ångström spatial resolution is achieved—remain underdeveloped, in part due to the lack of a rigorous theoretical description of the experiment. To address this, herein we develop a framework for differential solution scattering experiments conducted at high energies, which includes concepts of the solvent-excluded volume introduced to describe SAXS/WAXS data, as well as concepts from the time-resolved X-ray scattering community. Our theory is supported by numerical simulations and experiment and paves the way for establishing quantitative methods to determine the atomic structures of small molecules in solution with resolution approaching that of crystallography.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Safety First: Stability and Dissipation of Line-tied Force-free Flux Tubes in Magnetized Coronae

Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes ( m = 1) and efficient magnetic energy dissipation develop for plasma safety factors q ≲ 1, where q is the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes ( m > 1) can distort equilibrium flux tubes for q > 1 but induce significantly less dissipation. In our analysis, the characteristic pitch μ ˜ 0 of flux-tube field lines determines the growth rate ( ∝ μ ˜ 0 3 ) and minimum wavelength of the kink instability ( ∝ μ ˜ 0 − 1 ). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given μ ˜ 0 . By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗

DOE OSTI2024

A nonlinear journey from structural phase transitions to quantum annealing

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the “$\phi$ 4 model”) with weak inter-chain coupling, we explore connections between the driven variants of the two systems. Here, we argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical $\phi$ 4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical $\phi$ 4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗

DOE OSTI2025

Transport of Delocalized Excitons through DNA-Based Molecular Photonic Wires

Molecular photonic wires conduct electronic energy via their rapid transport properties. In photosynthesis, nature achieves efficient transport across large distances using delocalized excitons, generated by strong excitonic coupling between chromophores. How, or even whether, delocalization facilitates long-distance energy transport in synthetic systems has been challenging to experimentally test and optimize. Thus, far, studies have been limited to strongly coupled, heterogeneous chromophore aggregates or weakly coupled chromophore monomers. Here, in this work, we employed DNA nanostructures to engineer molecular photonic wires constructed from a series of excitonically coupled indocarbocyanine chromophores─achieving the intermediate and strong coupling regimes. Using time-resolved fluorescence spectroscopy and complementary simulations, we demonstrated that an intermediate intermolecular electronic coupling (∼ k B T ) enables up to 40% faster exciton transport as compared to strongly coupled chromophores. The delocalized excitons generated in the intermediate coupling regime exhibited properties conducive to rapid diffusivity, similar to their monomeric counterparts. Thus, intermediate excitonic coupling, analogous to natural systems, achieves long-distance exciton transport with the high chromophore density required for energy capture.

DNA origami↗

DOE OSTI2026

Final report: Insights from ARM observations into aerosol processing and transport by extratropical cyclones and aerosol effects on cyclone clouds

This project advanced understanding of aerosol–cloud interactions in extratropical cyclones using observations from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) observatory and ACE-ENA field campaigns, along with satellite datasets and numerical modeling. We have quantified the susceptibility of frontal clouds to aerosols, showing that this susceptibility is similar to that in non-frontal clouds despite substantial (factor 2) differences in droplet concentrations between these cloud types, and also significant differences in cloud albedo. A key finding was that Aitken-mode (here, 60-100nm diameter) aerosols play a disproportionately important role in droplet activation in frontal clouds, and concentrations of these particles are frequently strongly biased in climate models. Activation of aerosols depends strongly on updraft speeds, and via collaboration with Argonne National Laboratory we contributed to quantifying these better with the ARM radar wind profiler at the ENA site. Our studies in the North Atlantic have analogs in the Southern Ocean, and to draw these out we examined data from the MARCUS and SOCRATES field campaigns. We identified new particle formation within extratropical cyclones as a potentially important source of cloud condensation nuclei in these pristine marine environments. By combining ARM observations with perturbed parameter ensemble simulations, we demonstrated that surface observations can effectively constrain aerosol–cloud adjustments and reduce uncertainty in simulated cloud liquid water path responses by approximately 15%, yielding stronger constraints on historical aerosol cooling effects. The work further established precipitation processes as a dominant control on aerosol–cloud interactions and Earth system predictability.

Gordon, Hamish [Department of Chemical Engineering↗

DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Quasiaperiodic grain boundary phases of Σ⁢5 tilt grain boundaries in refractory metals

We report ground-state structures and phase transitions in Σ⁢5⁢[001] tilt grain boundaries (GBs) in body-centered-cubic (bcc) refractory metals Nb, Ta, Mo, and W. Σ⁢5 tilt GBs have been extensively investigated over the past several decades, with their ground-state structure—composed of kite-shaped structural units—previously thought to be well understood. By performing a rigorous GB structure search that optimizes the number of atoms in the boundary core, we predict different quasiaperiodic “split kite” phases analogous to those previously found in GBs in face-centered-cubic metals. Furthermore, our results suggest that complex aperiodic phases of GBs appear to be a general phenomenon, as validated through density functional theory calculations. Moreover, the atoms in the split kite phase demonstrate distinct collective diffusion dynamics. Phase-contrast image simulations of split kites show better agreement with experimental observations, offering an alternative explanation for previous microscopy results and motivating future atomically resolved imaging of the GB structure.

Body-centered cubic↗

DOE OSTI2024

Topological polarons in halide perovskites

Halide perovskites emerged as a revolutionary family of high-quality semiconductors for solar energy harvesting and energy-efficient lighting. There is mounting evidence that the exceptional optoelectronic properties of these materials could stem from unconventional electron–phonon couplings, and it has been suggested that the formation of polarons and self-trapped excitons could be key to understanding such properties. By performing first-principles simulations across the length scales, here we show that halide perovskites harbor a uniquely rich variety of polaronic species, including small polarons, large polarons, and charge density waves, and we explain a variety of experimental observations. We find that these emergent quasiparticles support topologically nontrivial phonon fields with quantized topological charge, making them nonmagnetic analog of the helical Bloch points found in magnetic skyrmion lattices.

Science & Technology - Other Topics↗

DOE OSTI2025

Quantifying the impact of precision errors on quantum approximate optimization algorithms

The quantum approximate optimization algorithm (QAOA) is a hybrid quantum-classical algorithm that seeks to achieve approximate solutions to optimization problems by iteratively alternating between intervals of controlled quantum evolution. Here, we examine the effect of analog precision errors on QAOA performance from the perspective of both algorithmic training and performance guarantees. Leveraging cumulant expansions, we recast the faulty QAOA as a control problem in which precision errors are expressed as multiplicative control noise and derive bounds on the performance of QAOA. We show using both analytical techniques and numerical simulations that fixed precision implementations of QAOA circuits are subject to an exponential degradation in performance dependent upon the number of optimal QAOA layers and magnitude of the precision error. Despite this significant reduction, we show that it is possible to mitigate precision errors in QAOA via digitization of the variational parameters at the cost of increasing circuit depth.

quantum algorithms↗