SEARCH · Search NASA

Results for “Direct simulation Monte Carlo”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Deep probabilistic direction prediction in 3D with applications to directional dark matter detectors

Abstract We present the first method to probabilistically predict 3D direction in a deep neural network model. The probabilistic predictions are modeled as a heteroscedastic von Mises-Fisher distribution on the sphere S 2 , giving a simple way to quantify aleatoric uncertainty. This approach generalizes the cosine distance loss which is a special case of our loss function when the uncertainty is assumed to be uniform across samples. We develop approximations required to make the likelihood function and gradient calculations stable. The method is applied to the task of predicting the 3D directions of electrons, the most complex signal in a class of experimental particle physics detectors designed to demonstrate the particle nature of dark matter and study solar neutrinos. Using simulated Monte Carlo data, the initial direction of recoiling electrons is inferred from their tortuous trajectories, as captured by the 3D detectors. For 40 keV electrons in a 70% He 30% CO 2 gas mixture at STP, the new approach achieves a mean cosine distance of 0.104 (26 ∘ ) compared to 0.556 (64 ∘ ) achieved by a non-machine learning algorithm. We show that the model is well-calibrated and accuracy can be increased further by removing samples with high predicted uncertainty. This advancement in probabilistic 3D directional learning could increase the sensitivity of directional dark matter detectors.

Computer Science↗

DOE OSTI2024

Data Science Shows that Entropy Correlates with Accelerated Zeolite Crystallization in Monte Carlo Simulations

We have performed a data science study of Monte Carlo simulation trajectories to understand factors that can accelerate formation of zeolite nanoporous crystals, a process that can take days or even weeks. In previous work, Monte Carlo simulations predicted and experiments confirmed that using a secondary organic structure-directing agent (OSDA) accelerates crystallization of all-silica LTA zeolite, with experiments finding a three-fold speedup [PCCP 24, 142-148 (2022)]. However, it remains unclear what physical factors cause the speed-up. Here, we apply data science to analyze the simulation trajectories to discover what drives accelerated zeolite crystallization in Monte Carlo going from a one-OSDA synthesis (1OSDA) to a two-OSDA version (2OSDA). We encoded simulation snapshots using the Smooth Overlap of Atomic Positions approach, which represents all 2- and 3-body correlations within a given cutoff distance. Principal component analyses failed to discriminate datasets of structures from 1OSDA and 2OSDA simulations, while the Support Vector Machine (SVM) approach succeeded at classifying such structures with an area-under-curve (AUC) score of 0.99 (where AUC = 1 is a perfect classification) with all 3-body correlations, and as high as 0.94 with only 2-body correlations. SVM decision functions reveal relatively broad / narrow histograms for 1OSDA / 2OSDA datasets, suggesting that the two simulations differ strongly in information heterogeneity. Informed by these results, we performed pair (2-body) entropy calculations during crystallization, resulting in entropy differences that semi-quantitatively account for the speedup observed in the previous Monte Carlo simulations. We conclude that altering synthesis conditions in ways that substantially changes the entropy of labile silica networks may accelerate zeolite crystallization, and we discuss possible approaches for achieving such acceleration.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗

DOE OSTI2025

Discovery of Stacking Heterogeneity, Layer Buckling, and Residual Water in COF-999-NH 2 and Implications on CO 2 Capture

Covalent organic frameworks (COFs), with their modular architectures and tunable functionalities, provide a versatile platform to design sorbents for the direct capture of CO 2 from air. Here, for this work, we combined density functional theory, molecular dynamics, and grand canonical Monte Carlo simulations with experiment to understand structural factors for furthering COF-999-NH 2 ’s performance as the precursor to COF-999 for direct air CO 2 capture. Small energy differences among laterally shifted stackings suggest intrinsic stacking heterogeneity. The simulations show pronounced layer buckling coupled to extensive amine–nitrile hydrogen bonding and persistent pore water, which initiates undesired polymerization and undermines uptake. The predicted presence of water is confirmed by subsequent experiments. These insights point to a single, actionable design rule: exclude retained water by introducing hydrophobic pore environments to maximize the CO 2 capture efficiency.

adsorption↗

DOE OSTI2024

Solar reflection of dark matter with dark-photon mediators

Abstract We consider the scattering of low-mass halo dark-matter particles in the hot plasma of the Sun, focusing on dark matter that interact with ordinary matter through a dark-photon mediator. The resulting “solar-reflected” dark matter (SRDM) component contains high-velocity particles, which significantly extend the sensitivity of terrestrial direct-detection experiments to sub-MeV dark-matter masses. We use a detailed Monte Carlo simulation to model the propagation and scattering of dark-matter particles in the Sun, including thermal effects, with special emphasis on ultralight dark-photon mediators. We study the properties of the SRDM flux, obtain exclusion limits from various direct-detection experiments, and provide projections for future experiments, focusing especially on those with silicon and xenon targets. We find that proposed future experiments with xenon and silicon targets can probe the entire “freeze-in benchmark”, in which dark matter is coupled to an ultralight dark photon, including dark-matter masses as low as 𝒪(keV). Our simulations and SRDM fluxes are publicly available.

Astronomy & Astrophysics↗

DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Transforming E-Waste Into Strategic Resources: Techno-Economic Analysis of Gallium and By-products Recovery from LEDs via Bioleaching

The growing demand for gallium in optoelectronics and renewable energy applications raises concerns about supply security and production sustainability. This study evaluates the techno-economic feasibility of recovering gallium and by-products (copper and nickel) from waste GaN-based LEDs via bioleaching. A process flowsheet encompassing transportation, robotic disassembly, ball milling, bioleaching, solvent extraction/electrowinning, and refining was modeled. Based on mass balance analysis, more than 53 tons of LED waste are required annually to yield 1 kg of gallium alongside substantial copper and nickel co-products. Preliminary techno-economic analysis (TEA) shows an average total cost (ATC) of 6.84 USD/kg metal when costs are allocated by mass-weighted economic value (market price) fraction, corresponding to 6.75 USD/kg for copper, 15.91 USD/kg for nickel, and 470.95 USD/kg for gallium. For gallium, direct operational costs account for more than 70% of the cost. Monte Carlo simulations further quantify cost uncertainties under market price fluctuations. This work represents the first TEA of gallium recovery from GaN-based LEDs and highlights potential pathways for future cost reduction.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Modeling hyperbranched polymer formation via ATRP using dissipative particle dynamics

Hyperbranched polymers (HBPs) offer distinguishing, advantageous properties that arise from their distinctive complex topology. One of the effective approaches to the synthesis of hyperbranched structures involves the use of a branching initiator (inibramer) that is activated only after incorporation into a polymer chain. There remain, however, challenges in determining and characterizing the structures of the synthesized HBPs. Dissipative particle dynamics (DPD) was used to probe the effects of inibramer concentration, solvent concentration, and inibramer reactivity on the kinetics, molecular weight, and dispersity of HBPs. Additionally, DPD allows for direct observation of branched structures, which was not possible in previously reported Monte Carlo type simulations. It was found that higher inibramer concentrations led to faster monomer consumption while forming more dendritic structures with fewer defects. Additionally, high dispersities characteristic of HBP systems were found to originate from asymmetric propagation rates between inibramer-inibramer and inibramer-monomer reactions.

Biswas, Santidan↗

DOE OSTI2026

Global Λ hyperon polarization in low-energy heavy ion collisions: A scenario without vorticity

Since its discovery, global polarization of the Λ hyperon in heavy-ion collisions has been firmly established and is widely attributed to the large vorticity generated in the rotating quark-gluon plasma. In contrast, nearly fifty years after the first observation of unexpectedly large transverse Λ polarization in unpolarized hadron collisions, its underlying mechanism remains an open and long-standing puzzle, despite being observed across a broad range of collision systems. Although these two phenomena exhibit notable similarities, they are generally regarded as arising from distinct physical origins. In this work, we propose a direct connection between Λ global polarization in heavy-ion collisions and the long-standing transverse polarization observed in unpolarized collision systems. We demonstrate that the alignment between the Λ production plane and the reaction plane, driven by directed flow, can transfer transverse polarization into the measured global polarization signal. Realistic Monte Carlo simulations of Au+Au collisions at √𝑠 NN =3 GeV indicate that this mechanism can generate a sizable global polarization, accounting for approximately 23% ±6% of the magnitude reported by the STAR Collaboration. Our results establish, for the first time, a quantitative link between these two well-known phenomena and have important implications for the interpretation of Λ global polarization measurements in low-energy heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Review of top quark mass measurements in CMS

The top quark mass is one of the most intriguing parameters of the standard model (SM). Its value indicates a Yukawa coupling close to unity, and the resulting strong ties to Higgs physics make the top quark mass a crucial ingredient for understanding essential aspects of the electroweak sector of the SM. This review offers the first comprehensive overview of the top quark mass measurements performed by the CMS Collaboration using the data collected at centre-of-mass energies of 7, 8, and 13 TeV. Moreover, a detailed description of the top quark event reconstruction is provided and dedicated studies of the dominant uncertainties in the modelling of the signal processes are discussed. The interpretation of the experimental results on the top quark mass in terms of the SM Lagrangian parameter is challenging and is a focus of an ongoing discussion in the theory community. The CMS Collaboration has performed two main types of top quark mass measurements, addressing this challenge from different perspectives: highly precise ‘direct’ measurements, based on reconstructed top quark decay products and relying exclusively on Monte-Carlo simulations, as well as ‘indirect’ measurements, where the simulations are employed to determine parton-level cross sections that are compared to fixed-order perturbative calculations. Recent mass extractions using Lorentz-boosted top quarks open a new avenue of measurements based on top quark decay products contained in a single particle jet, with promising prospects for accurate theoretical interpretations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Determining hexavalent chromium transport properties in alkaline nuclear waste using nuclear magnetic resonance spectroscopy

This study focuses on the transport properties of hexavalent chromium, specifically the chromate anion, to improve predictive models and environmental remediation strategies for Cr(VI) migration. Using 53 Cr Nuclear Magnetic Resonance (NMR) spectroscopy, the research quantifies chromate in multicomponent electrolytes replicating nuclear waste conditions at the Hanford Site in Washington State. The consistency of the 53 Cr NMR signal integral with chromate concentration, despite varying matrix compositions, establishes it as a reliable concentration indicator. The transport properties of chromate in an alkaline solution were assessed using relaxation-based measurements via saturation recovery and Carr-Purcell-Meiboom-Gill experiments, determining spin-lattice and spin-spin relaxation times. These measurements, combined with the Bloembergen-Purcell-Pound equation, helped estimate the rotational correlation time and the 53 Cr self-diffusion coefficient using Stokes-Einstein-Debye and Stokes-Einstein equations. Direct measurements were obtained through pulsed field gradient stimulated echo 53 Cr NMR spectroscopy. Monte Carlo simulations further estimated uncertainty propagation. The results enhance comprehension of chromate transport and highlight prospects for identifying transport properties of NMR-active nuclei, traditionally considered unreachable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗

DOE OSTI2026

Design and simulation of a muon detector to characterize geological overburden

This study presents the design, construction, and simulation of a mobile muon detector tailored for geological overburden characterization. The detector employs plastic scintillator paddles with silicon photomultipliers (SiPMs) and a QuarkNet data acquisition system, offering a portable solution suitable for remote field deployment. The simulator’s modular aluminum frame allows for adjustable geometry and directional sensitivity, while its battery system supports over a week of autonomous operation. Preliminary experimental tests confirmed that its muon flux measurements were consistent with theoretical expectations. A comprehensive simulation framework using Geant4 and CORSIKA was developed to model detector response and overburden effects. Analytical and Monte Carlo methods were used to assess quadrant resolution and infer muon directionality. This work lays the foundation for future overburden mapping and supports the development of reconstruction algorithms for geological applications.

72 - PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗