Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

DIGS: deep inference of galaxy spectra with neural posterior estimation

Abstract With the advent of billion-galaxy surveys with complex data, the need of the hour is to efficiently model galaxy spectral energy distributions (SEDs) with robust uncertainty quantification. The combination of simulation-based inference (SBI) and amortized neural posterior estimation (NPE) has been successfully used to analyse simulated and real galaxy photometry both precisely and efficiently. In this work, we utilise this combination and build on existing literature to analyse simulated noisy galaxy spectra. Here, we demonstrate a proof-of-concept study of spectra that is (a) an efficient analysis of galaxy SEDs and inference of galaxy parameters with physically interpretable uncertainties; and (b) amortized calculations of posterior distributions of said galaxy parameters at the modest cost of a few galaxy fits with Markov chain Monte Carlo (MCMC) methods. We utilise the SED generator and inference framework Prospector to generate simulated spectra, and train a dataset of 2 × 10 6 spectra (corresponding to a five-parameter SED model) with NPE. We show that SBI—with its combination of fast and amortized posterior estimations—is capable of inferring accurate galaxy stellar masses and metallicities. Our uncertainty constraints are comparable to or moderately weaker than traditional inverse-modelling with Bayesian MCMC methods (e.g. 0.17 and 0.26 dex in stellar mass and metallicity for a given galaxy, respectively). We also find that our inference framework conducts rapid SED inference (0.9–1.2 × 10 5 galaxy spectra via SBI/NPE at the cost of 1 MCMC-based fit). With this work, we set the stage for further work that focuses of SED fitting of galaxy spectra with SBI, in the era of JWST galaxy survey programs and the wide-field Roman Space Telescope spectroscopic surveys.

spectroscopy↗
DOE OSTI2025

Bayesian event categorization matrix approach for explosion monitoring

Current efforts to correctly categorize natural events from suspected explosion sources with data that is collected by ground- or space-based sensors presents historical challenges that remain unaddressed by the Event Categorization Matrix (ECM) model. Smaller historical events (lower yield explosions) may have data available from fewer measurement techniques than are available today, and therefore, a historical event record can lack a complete set of discriminants. The covariance structures can also differ between such observations of event (source-type) categories. Both obstacles are problematic for the classic ECM model. Our work addresses this gap and presents a Bayesian update to the previous ECM model, termed the Bayesian Event Categorization Matrix model, which can be trained on partial observations and does not rely on a pooled covariance structure. We further augment the ECM model with Bayesian Decision Theory so that false negative or false positive rates of an event categorization can be reduced in an intuitive manner. To demonstrate improved categorization rates for the Bayesian Event Categorization Matrix model, we compare an array of Bayesian and classic models with multiple performance metrics using Monte Carlo experiments. We use both synthetic and real data. Our Bayesian models show consistent gains in overall accuracy and lower false negative rates relative to the classic ECM model. Here, we propose future avenues to improve Bayesian Event Categorization Matrix models’ decision making and predictive capability.

58 GEOSCIENCES↗
DOE OSTI2023

The GMRT archive atomic gas survey – II. Mass modelling and dark matter halo properties across late-type spirals

Studying the kinematics and mass modelling of galaxies from H i 21 cm data provides valuable insights into the properties of both the baryonic components and the dark matter halo in nearby galaxies. Despite many observational studies, mass modelling of galaxies remains challenging due to different limitations. For example, most of the previous studies involving mass modelling are based on rotation curves derived from 2D velocity fields from H i or H α spectroscopic observation which are often affected by beam smearing and projection effect. However, kinematic modelling done by fitting the ‘Tilted ring model’ to 3D data cube is not affected by these issues. In this study, we present and compare 3D kinematic modelling of a pilot sample of 11 galaxies from the GMRT archive atomic gas survey (GARCIA) using two different publicly available pipelines. We model the observed H i rotation curve using 3.6-μm infrared data and SDSS r-band data for stellar contribution, H i surface density profile for gas, and Navarro–Frenk–White profile for dark matter halo; and employ the Markov chain Monte Carlo optimization method for parameter estimation. Further, to validate our analysis, we revisit important scaling relations, e.g. the M gas –M star relation, M star –M halo relation, M gas –M halo relation and Baryonic Tully–Fisher relation. Here, the scaling relations from our analysis are broadly consistent with that reported in the literature. A larger sample of galaxies from GARCIA in the near future will allow studying these scaling relations in greater details.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measurement of the branching fraction of Λ + 𝑐 →𝑝⁢𝐾 0 𝑆⁢ 𝜋 0 at Belle

We report a precise measurement of the ratio of branching fractions ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾$^0_𝑆$⁢𝜋 0 )/ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ) using 980 fb −1 of 𝑒 + ⁢𝑒 − data from the Belle experiment. We obtain a value of ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾$^0_𝑆$⁢𝜋 0 )/ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ) =0.339 ±0.002 ±0.009, where the first and second uncertainties are statistical and systematic, respectively. This Belle result is consistent with the previous measurement from the CLEO experiment but has a fivefold improvement in precision. By combining our result with the world average ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ), we obtain the absolute branching fraction ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾$^0_𝑆$⁢𝜋 0 ) =(2.12 ± 0.01 ±0.05 ±0.10)%, where the uncertainties are statistical, systematic, and the uncertainty in the absolute branching fraction scale ℬ⁡(Λ$^+_𝑐$ → 𝑝⁢𝐾 − ⁢𝜋 + ), respectively. This measurement can shed light on hadronic decay mechanisms in charmed baryon decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Determination of proton PDF uncertainties with Markov chain Monte Carlo

We present an analysis of parton distribution functions (PDFs) of the proton using Markov chain Monte Carlo (MCMC) methods. The MCMC approach naturally implements Bayes’ theorem and, thus, provides a means to directly sample the underlying probability distribution—in this case, the probability distribution of the PDF parameters. This allows for a straightforward propagation of the resulting uncertainties into any PDF-dependent observable, preserving their simple probabilistic interpretation. In our analysis we include a broad set of deep inelastic scattering data from HERA, BCDMS and NMC experiments along with the Drell-Yan, 𝑊 and 𝑍 boson data from LHC and Tevatron experiments, which combined with theoretical calculations at next-to-next-to-leading order in QCD allow for realistic determination of PDFs. The main focus of this analysis is to explore alternative methods for PDF uncertainty estimation that are more firmly grounded in statistical principles. We show that the flexibility of the Bayes framework, allowing one, e.g., to account for non-Gaussianity or inconsistencies of datasets, is crucial to extract realistic uncertainties when such assumptions are not fulfilled. We also demonstrate that MCMC allows one to determine the Δ⁢𝜒 2 value corresponding to a given confidence level in the sample, which can, in turn, be used as a statistically well-founded tolerance criterion used in the Hessian method, thus addressing one of its main long-standing drawbacks.

Risse, Peter Clemens [Universität Münster (Germany↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
DOE OSTI2026

Baryon–antibaryon photoproduction cross sections off the proton

The GlueX experiment at Jefferson Lab has observed 𝑝⁢$\overline{p}$ and, for the first time, Λ$\overlineΛ$ and 𝑝$\overlineΛ$ photoproduction from a proton target at photon energies up to 11.6 GeV. The angular distributions are forward peaked for all produced pairs, consistent with Regge-like 𝑡-channel exchange. Asymmetric wide-angle antibaryon distributions show the presence of additional processes. In a phenomenological model, we find consistency with a double-𝑡-channel exchange process where antibaryons are created only at the middle vertex. The model matches all observed distributions with a small number of free parameters. In the hyperon channels, we observe a clear distinction between photoproduction of the Λ$\overlineΛ$ and 𝑝$\overlineΛ$ systems but general similarity to the 𝑝$\overline{p}$ system. We report both total cross sections and cross sections differential with respect to momentum transfer and the invariant masses of the created particle pairs. No narrow resonant structures were found in these reaction channels. The suppression of 𝑠$\overline{s}$ quark pairs relative to 𝑑$\overline{d}$ quark pairs is similar to what has been seen in other reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

External Magnetic Field Suppression of Carbon Diffusion in Iron

External magnetic fields reduce diffusion of carbon in BCC iron, but the physical mechanism is not understood. Using DFT calculations with magnetic moments sampled from a Heisenberg model, we calculate diffusivities of carbon in iron at high temperatures and with field. Our model reproduces the measured suppression of diffusivity from field. Here, we find that increasing magnetic disorder flattens the electron density of states compared with the ferromagnetic case, which distorts the octahedral cages around carbon, lowering the activation barrier to diffusion; an applied field reverses these trends.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stretched exponential magnetic relaxation dynamics in artificial square ice revealed through x-ray photon correlation spectroscopy

X-ray photon correlation scattering measurements are undertaken on a thermally active artificial spin ice based on the square lattice, referred to as artificial square ice, to probe the fluctuation timescales as a function of temperature as the system passes through the paramagnetic-antiferromagnetic phase transition, which belongs to the two-dimensional Ising universality class. In the paramagnetic regime, a single exponential timescale is seen, whereas at and below the critical temperature, a stretched exponential decorrelation is observed, with the stretching exponent decreasing from unity down to below one-half as the temperature reduces. This trend is confirmed by kinetic Monte Carlo simulations of a simplified point-dipolar square ice system, and is in agreement with past theoretical work on the kinetic Ising model where stretched exponential relaxation due to equilibrium domain wall dynamics below the critical temperature was found.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Indirect Measurement of the 23 Na(p, γ) 24 Mg Direct Capture Reaction Rate via ( 3 He,d) Spectroscopy

The cross section of the $^{23}\text{Na}(p,γ)^{24}\text{Mg}$ reaction is dominated by direct capture at low energies relevant for stellar burning. Such cross sections can be constrained using spectroscopic factors($C^2S$) or asymptotic normalization coefficients(ANCs) from transfer reactions. In this work, the $^{23}\text{Na}(^3\text{He},d)^{24}\text{Mg}$ reaction was measured at $E_{lab}=21$ MeV to extract spectroscopic factors for $^{24}\text{Mg}$ states with excitation energies in $E_x=7 - 12$ MeV using the Enge split-pole spectrograph at the Triangle Universities Nuclear Laboratory. A new non-resonant astrophysical S factor and the direct capture reaction rate for the $^{23}\text{Na}(p,γ)$ reaction are calculated and presented based on this measurement. The new rate at $T<0.04$ GK is 43$\%$ smaller than in previous studies. Rigorous treatments of uncertainties are presented using a Bayesian Markov Chain Monte Carlo (MCMC) method. Finally, sources of uncertainties for computing the direct capture cross section are also discussed in detail.

electromagnetic moments↗
DOE OSTI2022

Efficient quantum state tomography with mode-assisted training

Neural networks (NNs) representing quantum states are typically trained using Markov chain Monte Carlo based methods. However, unless specifically designed, such samplers only consist of local moves, making the slow-mixing problem prominent even for extremely simple quantum states. Here, we propose to use mode-assisted training that provides global information via the modes of the NN distribution. Further, applied to quantum state tomography using restricted Boltzmann machines, this method improves the quality of reconstructed quantum states by orders of magnitude. The method is applicable to other types of NNs and may efficiently tackle problems previously unmanageable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Classical analog of quantum models in synthetic dimensions

We introduce a classical analog of quantum matter in ultracold molecule-synthetic or Rydberg atom-synthetic dimensions, by extending the Potts model to include interactions J 1 between atoms adjacent in both real and synthetic space and studying its finite-temperature properties. For intermediate values of J 1 , the resulting phases and phase diagrams are similar to those of the clock and Villain models, in which three phases emerge. There exists a sheet phase analogous to that found in quantum synthetic dimension models between the high-temperature disordered phase and the low-temperature ferromagnetic phase. Furthermore, we also employ machine learning to uncover nontrivial features of the phase diagram using the learning by confusion approach, which is able to discern several successive phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum parton shower with kinematics

Parton showers, which can efficiently incorporate quantum interference effects, have been shown to be run efficiently on quantum computers. However, so far these quantum parton showers have not included the full kinematical information required to reconstruct an event, which in classical parton showers requires the use of a veto algorithm. Here, in this work, we show that adding one extra assumption about the discretization of the evolution variable allows one to construct a quantum veto algorithm, which reproduces the full quantum interference in the event, and allows one to include kinematical effects. We finally show that for certain initial states the quantum interference effects generated in this veto algorithm are classically tractable, such that an efficient classical algorithm can be devised.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Generalizing Similarity Laws for Radio-Frequency Discharge Plasmas across Nonlinear Transition Regimes

We generalize similarity theory based on the scaling and solution invariance of the Boltzmann equation, coupled with the Poisson equation, and demonstrate similarity laws for radio-frequency (rf) discharge plasmas across three nonlinear transitional regimes, namely, the alpha-gamma mode transition, the stochastic-Ohmic-heating mode transition, and the bounce-resonance-heating mode transition. Fundamental plasma parameters, e.g., the electron power absorption, under similar discharge conditions are examined via fully kinetic particle-in-cell simulations, and electron-kinetic invariance is exemplified in similar rf discharge plasmas. The results unambiguously confirm the applicability of similarity laws for rf plasmas in extended operating regimes, and strengthen the foundations and framework of similarity physics with universality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Origin of mean-field behavior in an elastic Ising model

Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. Here, we demonstrate how the long-range interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean-field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean-field theory yields similarly accurate results.

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

Exotic bonding interactions and coexistence of chemically distinct periodic lattice distortions in the charge density wave compound TaTe2

The superstructure of TaTe 2 is examined by total high-energy x-ray scattering and large-scale structure modeling over a broad temperature range. At room temperature, it features double zigzag chains of Ta atoms exhibiting metallic bonding. The chains are sandwiched between planes of weakly interacting Te atoms. Below 170 K, the chains appear broken to “butterflylike” clusters and single zigzag chains of covalently bonded Te atoms emerge while Ta-Te interactions remain largely ionic. This is a rare example of a charge density wave compound stabilized by bonding that is a mix of covalent, ionic, and metallic character giving rise to chemically distinct periodic lattice distortions. We argue that the formation of periodic lattice distortions in the Te sublattice in addition to those in the Ta sublattice favors charge delocalization, which may explain the unusually diminished resistivity and increased magnetic susceptibility exhibited by the low-temperature charge density wave phase of TaTe 2 . The effect may be common for 5d transition metal chalcogenides containing Te because of the extended electronic orbitals of the constituent atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗