Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Imaging elliptically polarized infrared near-fields on nanoparticles by strong-field dissociation of functional surface groups

Abstract We investigate the strong-field ion emission from the surface of isolated silica nanoparticles aerosolized from an alcoholic solution, and demonstrate the applicability of the recently reported near-field imaging at 720 nm [Rupp et al., Nat. Comm., 10(1):4655, 2019] to longer wavelength (2 $$\mu $$ μ m) and polarizations with arbitrary ellipticity. Based on the experimental observations, we discuss the validity of a previously introduced semi-classical model, which is based on near-field driven charge generation by a Monte-Carlo approach and classical propagation. We furthermore clarify the role of the solvent in the surface composition of the nanoparticles in the interaction region. We find that upon injection of the nanoparticles into the vacuum, the alcoholic solvent evaporates on millisecond time scales, and that the generated ions originate predominantly from covalent bonds with the silica surface rather than from physisorbed solvent molecules. These findings have important implications for the development of future theoretical models of the strong-field ion emission from silica nanoparticles, and the application of near-field imaging and reaction dynamics of functional groups on isolated nanoparticles. Graphical abstract

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hydroxide films on mica form charge-stabilized microphases that circumvent nucleation barriers

Crystal nucleation is facilitated by transient, nanoscale fluctuations that are extraordinarily difficult to observe. Here, we use high-speed atomic force microscopy to directly observe the growth of an aluminum hydroxide film from an aqueous solution and characterize the dynamically fluctuating nanostructures that precede its formation. Nanoscale cluster distributions and fluctuation dynamics show many similarities to the predictions of classical nucleation theory, but the cluster energy landscape deviates from classical expectations. Kinetic Monte Carlo simulations show that these deviations can arise from electrostatic interactions between the clusters and the underlying substrate, which drive microphase separation to create a nanostructured surface phase. This phase can evolve seamlessly from a low-coverage state of fluctuating clusters into a high-coverage nanostructured network, allowing the film to grow without having to overcome classical nucleation barriers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Measurement of the 1-jettiness event shape observable in deep-inelastic electron-proton scattering at HERA

The H1 Collaboration reports the first measurement of the 1-jettiness event shape observable $τ^b_1$ in neutral-current deep-inelastic electron-proton scattering (DIS). The observable $τ^b_1$ is equivalent to a thrust observable defined in the Breit frame. The data sample was collected at the HERA ep collider in the years 2003–2007 with center-of-mass energy of $\sqrt{s}$ = 319 GeV, corresponding to an integrated luminosity of 351.1 pb -1 . Triple differential cross sections are provided as a function of $τ^b_1$, event virtuality $Q^2$, and inelasticity y, in the kinematic region $Q^2$ > 150 GeV 2 . Single differential cross section are provided as a function of $τ^b_1$ in a limited kinematic range. Double differential cross sections are measured, in contrast, integrated over $τ^b_1$ and represent the inclusive neutral-current DIS cross section measured as a function of $Q^2$ and y. The data are compared to a variety of predictions and include classical and modern Monte Carlo event generators, predictions in fixed-order perturbative QCD where calculations up to $\mathcal{O}$($α^3_s$) are available for $τ^b_1$ or inclusive DIS, and resummed predictions at next-to-leading logarithmic accuracy matched to fixed order predictions at $\mathcal{O}$($α^2_s$). These comparisons reveal sensitivity of the 1-jettiness observable to QCD parton shower and resummation effects, as well as the modeling of hadronization and fragmentation. Within their range of validity, the fixed-order predictions provide a good description of the data. Monte Carlo event generators are predictive over the full measured range and hence their underlying models and parameters can be constrained by comparing to the presented data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo

Stochastic reaction network models are often used to explain and predict the dynamics of gene regulation in single cells.These models usually involve several parameters, such as the kinetic rates of chemical reactions, that are not directly measurable and must be inferred from experimental data. Bayesian inference provides a rigorous probabilistic frame-work for identifying these parameters by finding a posterior parameter distribution that captures their uncertainty.Traditional computational methods for solving inference problems such as Markov Chain Monte Carlo methods based on classical Metropolis-Hastings algorithm involve numerous serial evaluations of the likelihood function, which in turn requires expensive forward solutions of the chemical master equation (CME). We propose an alternate approach based on a multifidelity extension of the Sequential Tempered Markov Chain Monte Carlo (ST-MCMC) sampler. This algorithm is built upon Sequential Monte Carlo and solves the Bayesian inference problem by decomposing it into a sequence of efficiently solved subproblems that gradually increase both model fidelity and the influence of the observed data. We reformulate the finite state projection (FSP) algorithm, a well-known method for solving the CME, to produce a hierarchy of surrogate master equations to be used in this multifidelity scheme. To determine the appropriate fidelity, we introduce a novel information-theoretic criteria that seeks to extract the most information about the ultimate Bayesian posterior from each model in the hierarchy without inducing significant bias. This novel sampling scheme is tested with high performance computing resources using biologically relevant problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The utility of Bayesian data reconciliation for separations

Data reconciliation methods for separation processes typically rely on classical statistical approaches to generate estimates of true mass flow rates from measurements. Knowledge regarding the uncertainty of these estimates has value in decision making, but is often not acquired. Bayesian approaches intrinsically quantify uncertainty; however, literature for Bayesian data reconciliation of separation processes is scarce. This publication outlines two Bayesian data reconciliation models and provides details for how the models were implemented for the BayesMassBal (V 1.0.0) software package written in R. To demonstrate the advantages of this approach for data reconciliation, the models were first applied to simulated data and then compared to a classical model through a Monte Carlo experiment. In this example, the Bayesian models were found to provide more accurate estimates of the simulated data, while also providing quantitative information on the estimate uncertainty. To demonstrate the use of the technique in a practical problem, the models were also applied to real data collected from a pilot-scale rare earth solvent extraction process. Here, this publication provides a small window into how Bayesian methods can be used for data reconciliation, but findings suggest Bayesian data reconciliation models for separation processes have distinct advantages over classical alternatives.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Composite Qdrift-product formulas for quantum and classical simulations in real and imaginary time

Recent study has shown that it can be advantageous to implement a composite channel that partitions the Hamiltonian H for a given simulation problem into subsets A and B such that H = A + B , where the terms in A are simulated with a Trotter-Suzuki channel and the B terms are randomly sampled via the Qdrift algorithm. Here we extend Qdrift and composite product formulas to imaginary time, formulating candidate classical algorithms for quantum Monte Carlo calculations. We upper bound the induced Schatten- 1 → 1 norm on both imaginary-time Qdrift and composite channels. Another recent result demonstrated that simulations of lattice Hamiltonians containing geometrically local interactions can be improved using a Lieb-Robinson argument to decompose H into subsets that contain only terms supported on that subset of the lattice. Here, we provide a quantum algorithm by unifying this result with the composite approach into “local composite channels” and we upper bound the diamond distance. We provide exact numerical simulations of algorithmic cost by counting the number of gates of the form e − i H j t and e − H j β to meet a certain error tolerance ε . In doing so, we optimize the partitioning into sets A and B using gradient boosted tree models from machine learning. These numerical studies are important given that product formulas have been historically known to outperform analytic upper bounds. We show constant factor advantages for a variety of interesting Hamiltonians, the maximum of which is a ≈ 20 -fold speedup that occurs in the simulation of Jellium. Published by the American Physical Society 2024

Pocrnic, Matthew (ORCID:0000000203089376)↗
DOE OSTI2022

Essential role of magnetic frustration in the phase diagrams of doped cobaltites

Doped perovskite cobaltites (e.g., La 1–x Sr x CoO 3 ) have been extensively studied for their spin-state physics, electronic inhomogeneity, and insulator-metal transitions. Ferromagnetically interacting spin-state polarons emerge at low x in the phase diagram of these compounds, eventually yielding long-range ferromagnetism. The onset of long-range ferromagnetism (x≈0.18) is substantially delayed relative to polaron percolation (x≈0.05), however, generating a troubling inconsistency. Here, in this study, Monte Carlo simulations of a disordered classical spin model are used to establish that previously ignored magnetic frustration is responsible for this effect, enabling faithful reproduction of the magnetic phase diagram.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Efficient Probabilistic Visualization of Local Divergence of 2D Vector Fields with Independent Gaussian Uncertainty

This work focuses on visualizing uncertainty of local divergence of two-dimensional vector fields. Divergence is one of the fundamental attributes of fluid flows, as it can help domain scientists analyze potential positions of sources (positive divergence) and sinks (negative divergence) in the flow. However, uncertainty inherent in vector field data can lead to erroneous divergence computations, adversely impacting downstream analysis. While Monte Carlo (MC) sampling is a classical approach for estimating divergence uncertainty, it suffers from slow convergence and poor scalability with increasing data size and sample counts. Thus, we present a two-fold contribution that tackles the challenges of slow convergence and limited scalability of the MC approach. (1) We derive a closed-form approach for highly efficient and accurate uncertainty visualization of local divergence, assuming independently Gaussian-distributed vector uncertainties. (2) We further integrate our approach into Viskores, a platform-portable parallel library, to accelerate uncertainty visualization. In our results, we demonstrate significantly enhanced efficiency and accuracy of our serial analytical (speed-up up to 1946×) and parallel Viskores (speed-up up to 19698×) algorithms over the classical serial MC approach. We also demonstrate qualitative improvements of our probabilistic divergence visualizations over traditional mean-field visualization, which disregards uncertainty. We validate the accuracy and efficiency of our methods on wind forecast and ocean simulation datasets.

Ouermi, Timbwaoga [University of Utah]↗
DOE OSTI2026

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2022

Equation of State Optimization and Uncertainty Quantization: Implementation in the LANL EOS Production Code OpenSesame

We detail the approaches of particle swarm optimization and Bayesian inference through Markov chain Monte Carlo for equation of state development. This work includes formulation of the modeling for the equation of state, numeric optimization of the parametric models via particle swarm optimization, and generation of probability distributions of equations of state from Markov chain Monte Carlo.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hydrogen Burning of 29 Si and Its Impact on Presolar Stardust Grains from Classical Novae

Abstract Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small N ( 12 C)/ N ( 13 C) and N ( 14 N)/ N ( 15 N) number abundance ratios compared to solar values, and abundance excesses in 30 Si relative to 29 Si, as previously predicted by models of classical novae. We report on a direct measurement of the 29 Si(p, γ ) 30 P reaction, which strongly impacts simulated δ 29 Si values from classical novae. Our new experimental 29 Si(p, γ ) 30 P thermonuclear reaction rate differs from previous results by up to 50% in the classical nova temperature range ( T = 100–400 MK), while the rate uncertainty is reduced by up to a factor of 3. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate δ 29 Si values with much reduced uncertainties. Our results establish δ 29 Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that δ 30 Si values from nova simulations are currently not a useful diagnostic tool unless the large uncertainty of the 30 P(p, γ ) 31 S reaction rate can be significantly reduced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗
DOE OSTI2021

Local structure and distortions of mixed methane-carbon dioxide hydrates

A vast source of methane is found in gas hydrate deposits, which form naturally dispersed throughout ocean sediments and arctic permafrost. Methane may be obtained from hydrates by exchange with hydrocarbon byproduct carbon dioxide. It is imperative for the development of safe methane extraction and carbon dioxide sequestration to understand how methane and carbon dioxide co-occupy the same hydrate structure. Pair distribution functions (PDFs) provide atomic-scale structural insight into intermolecular interactions in methane and carbon dioxide hydrates. We present experimental neutron PDFs of methane, carbon dioxide and mixed methane-carbon dioxide hydrates at 10 K analyzed with complementing classical molecular dynamics simulations and Reverse Monte Carlo fitting. Mixed hydrate, which forms during the exchange process, is more locally disordered than methane or carbon dioxide hydrates. The behavior of mixed gas species cannot be interpolated from properties of pure compounds, and PDF measurements provide important understanding of how the guest composition impacts overall order in the hydrate structure.

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

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗
DOE OSTI2023

X-ray induced electron and ion fragmentation dynamics in IBr

Characterization of the inner-shell decay processes in molecules containing heavy elements is key to understanding x-ray damage of molecules and materials and for medical applications with Auger-electron-emitting radionuclides. The 1s hole states of heavy atoms can be produced by absorption of tunable x rays and the resulting vacancy decays characterized by recording emitted photons, electrons, and ions. The 1s hole states in heavy elements have large x-ray fluorescence yields that transfer the hole to intermediate electron shells that then decay by sequential Auger-electron transitions that increase the ion’s charge state until the final state is reached. In molecules, the charge is spread across the atomic sites, resulting in dissociation to energetic atomic ions. We have used x-ray/ion coincidence spectroscopy to measure charge states and energies of Iq+ and Brq′+ atomic ions following 1s ionization at the I and Br K-edges of IBr. We present the charge states and kinetic energies of the two correlated fragment ions associated with core-excited states produced during the various steps of the cascades. Here, to understand the dynamics leading to the ion data, we develop a computational model that combines Monte-Carlo/Molecular-Dynamics (MC/MD) simulations with a classical over-the-barrier model to track inner-shell cascades and redistribution of electrons in valence orbitals and nuclear motion of fragments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗