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
NASA NTRS2003

Ultraviolet and Visible Emission Mechanisms in Astrophysics

The project involved the study of ultraviolet (UV) and visible emission mechanisms in astrophysical and atmospheric environments. In many situations, the emission is a direct consequence of a charge transferring collision of an ion with a neutral with capture of an electron to an excited state of the product ion. The process is also important in establishing the ionization and thermal balance of an astrophysical plasma. As little of the necessary collision data are available, the main thrust of the project was the calculation of total and state-selective charge transfer cross sections and rate coefficients for a very large number of collision systems. The data was computed using modern explicit techniques including the molecular-orbital close-coupling (MOCC), classical trajectory Monte Carlo (CTMC), and continuum distorted wave (CDW) methods. Estimates were also made in some instances using the multichannel Landau-Zener (MCLZ) and classical over-the-barrier (COB) models. Much of the data which has been computed has been formatted for inclusion in a charge transfer database on the World Wide Web (cfadc.phy.ornl.gov/astro/ps/data/). A considerable amount of data has been generated during the lifetime of the grant. Some of it has not been analyzed, but it will be as soon as possible, the data placed on our website, and papers ultimately written.

Stancil, Phillip C.↗
DOE OSTI2023

A high-performance implementation of atomistic spin dynamics simulations on x86 CPUs

Atomistic spin dynamics simulations provide valuable information about the energy spectrum of magnetic materials in different phases, allowing one to identify instabilities and the nature of their excitations. However, the time cost of evaluating the dynamical correlation function S(q,t) increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. In this work, here, we propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel, which results in further performance improvements of 44\% - 71\% on several relatively large lattice sizes when compared to the implementation that uses the GEMM subroutine in OpenBLAS, which is the state-of-the-art open source library for Basic Linear Algebra Subroutine (BLAS).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core

In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Double differential distributions of e-emission in ionization of N 2 by 3,4 and 5 keV electron impact

In this paper, we report the measurement of the absolute double differential cross sections (DDCS) of secondary electrons emitted due to the ionization of N2 molecule in collisions with fast electrons having energies between 3 and 5 keV. The emitted electrons with energies from 1-500 eV have been measured for different forward and backward emission angles. The measured DDCS have been compared with the state-of-the-art first Born approximation with correct boundary condition (CB1) model calculations as well as with the classical trajectory Monte Carlo (CTMC) method. From the measured DDCS, the single differential cross sections (SDCS) as a function of the emission energies have been computed and eventually the total ionization cross sections (TCS) have been derived. The TCS values are also compared with a semi-empirical calculation, namely, the CSP-ic (complex scattering potential-ionization contribution) model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Dipolar focusing in laser-assisted positronium formation

Abstract We consider positronium formation in collisions of positrons with excited hydrogen atoms H(n) in an infrared laser field theoretically. This process is assisted by the dipolar focusing effect: a positron moving in a superposition of a laser field and the dipolar field can approach the atomic target even if its trajectory starts with a very large impact parameter, leading to a significant enhancement of the Ps formation cross section. The classical trajectory Monte Carlo method, which is justified for n ⩾ 3 , allows efficient calculation of this enhancement. A similar effect can occur in collisions of positrons with other atoms in excited states, which can lead to improvements in the efficiency of positronium formation.

Optics↗
DOE OSTI2025

Enhanced quantum radiation with flying-focus laser pulses

The emission of a photon by an electron in an intense laser field is one of the most fundamental processes in electrodynamics and underlies the many applications that utilize high-energy photon beams. This process is typically studied for electrons colliding head-on with a stationary-focus laser pulse. Here, we show that the energy lost by electrons in the quantum regime and the yield of emitted photons can be substantially increased by replacing a stationary-focus pulse with an equal-energy flying-focus pulse whose focus co-propagates with the electrons. Furthermore, these advantages of the flying focus result from the energy loss and the photon yield scaling more favorably with the interaction time than the laser intensity in the quantum regime, with the latter also holding in the classical regime. Monte Carlo simulations of electrons colliding with equal-energy stationary and flying-focus laser pulses demonstrate these advantages.

Compton scattering↗
DOE OSTI2020

Effect of the orientation of Rydberg atoms on their collisional ionization cross section

Collisional ionization between two Rydberg atoms in relative motion is examined. A classical trajectory Monte Carlo method is used to determine the cross sections associated with Penning ionization. The dependence of the ionization cross section on the magnitude and the direction of orbital angular momentum of the electrons and the direction of the Laplace-Runge-Lenz vector of the electrons is studied. For a given magnitude of angular momentum, there can exist a difference of a factor of up to ~2.5 in the ionization cross section between the orientations with the highest and the lowest ionization cross section. As a result, the case of exchange ionization is examined and its dependence on the magnitude of angular momentum is studied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Waveform-dependent air fluorescence from neutral and ionic nitrogen molecules

Laser-induced air fluorescence in the ultraviolet regime is primarily attributed to transitions between the C and B states in excited neutral N 2 molecules and between the B and X states in N$^+_2$ ions. However, the mechanism underlying the former remains contentious, as direct population to the C state by light fields is forbidden by electron spin constraints. In this work, we investigate the mechanism of air fluorescence from excited neutral N 2 molecules by carrier-envelope phase–stabilized sub–4 femtosecond pulses. Our results show that fluorescence from N$^+_2$ ions reaches a maximum with cosine-like pulses, while fluorescence from excited neutral N 2 molecules peaks with sine-like pulses. In addition, by scanning the chirp of the driving pulse, we find that ionic fluorescence is maximized with chirp-free pulses, whereas neutral fluorescence favors negatively chirped pulses. These observations, supported by classical trajectory Monte Carlo simulations, support the mechanism of intersystem crossing from excited spin-singlet states, which are populated via recollision-induced strong-field excitation.

Science & Technology - Other Topics↗
DOE OSTI2022

All-optical nanoscopic spatial control of molecular reaction yields on nanoparticles

Molecular adsorbate reactions on nanoparticles play a fundamental role in areas such as nano-photocatalysis, atmospheric, and astrochemistry. They can be induced, enhanced, and controlled by field localization and enhancement on the nanoparticle surface. In particular, the ability to perform highly controlled near-field-mediated reactions is key to deepening our understanding of surface photoactivity on nanosystems. Here, using reaction nanoscopy, we experimentally demonstrate all-optical nanoscopic control of surface reaction yields by tailoring the near fields on nanoparticles with waveform-controlled linear and bicircular two-color laser pulses, respectively. We observe site-selective proton emission from the dissociative ionization of adsorbate molecules on Si O 2 nanoparticles as a function of the polarization and relative phase of the two-color pulses. The angularly resolved close-to-uniform mapping between the surface reaction yields and the measured ion momentum enables the observation and spatial control of molecular reactions on the nanoparticle surface with nanoscopic resolution. The experimental results are modeled and reproduced qualitatively by classical trajectory Monte Carlo simulations. Our work paves the way toward reliable all-optical control of photocatalytic chemical reactions on nanoscale surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

The Oort cloud in transition

The evolution of theoretical and empirical models of the Oort cloud (OC) since it was first proposed by Oort in 1950 is traced, and the main features of current models are discussed, in a general review. Consideration is given to work on the classical OC (Monte Carlo simulations of OC evolution, population and mass estimates, and OC perturbation by passing stars and giant molecular clouds), models of a massive inner OC (simulations of planetesimal-swarm evolution in the Uranus-Neptune zone and IRAS observations of circumstellar dust shells), evidence for random and/or periodic comet showers, and the possible role of the Galactic missing mass. The current OC model comprises an almost spherical outer (10,000-100,000-AU) cloud of mass 7-8 earth mass and population (1.4-2.3) x 10 to the 12th, and a flat disklike inner (40-10,000 AU) cloud of mass 100-200 earth mass and population (1-10) x 10 to the 13th.

Weissman, P. R.↗
NASA NTRS2002

Charge Transfer Calculations and Database for Astrophysics

A variety of theoretical approaches, having different but overlapping energy ranges of applicability, are applied to investigate charge transfer processes for collisions of atomic ions with atoms and molecules. The methods include quantal molecular-orbital close-coupling, classical trajectory Monte Carlo, and continuum distorted wave methods. Recent collision systems studied include S(+4) + H, S(+4) + He, N(+7) + He, H2O, CO, and CO2, O(+q)(q = 1 - 8) + H, H2, and S(+q)(q = 1 - 16) + H2. The database effort is concentrating on astrophysically important reactions of atomic ions X(+q)(X=H-Zn, q=1-4, and selected higher charges) with H, He, various metal atoms, H2, and other selected molecular targets. Existing data, much of it produced by us, has been compiled and critically evaluated. Data for many reactions missing in the literature are estimated using the multichannel Landau-Zener approximation. Fits to cross sections and rate coefficients using standard functions are provided as well as tabulations of the raw data. The database is available on the World Wide Web at cfadc.phy.ornl.gov/astro/ps/data.

Wang, J. G.↗
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↗
NASA NTRS2019

Characterizing Lunar Polar Volatiles at the Working Scale: Going from Exploration Goals to Mission Requirements

The economic evaluation of natural resources depends on the accuracy of resource distribution estimates. On Earth such estimates are necessary in making decisions about opening new mines or in planning future investment for operating mines or industrial deposits. A frequently discussed lunar resource is water ice, however, we are only at the first stages of understanding its potential as a resource. In particular, we currently do not have a sufficient understanding of the distribution of water or its form at the scales it would be extracted and processed, that is, the “working scale”. Here the “working scale” is defined to be the scales at which sufficient material can be processed to meet some basic demand (for example, 100s of square meters), and the anticipated heterogeneity in the water distribution across those scales (scales <5 - 10s of meters). Several mission concepts have been developed to better understand lunar water, motivated by both scientific and exploration goals. This paper provides an analysis of the number and distribution of observations needed to provide the necessary next steps in lunar water ISRU. We use a combination of Monte Carlo studies and classic geostatistical approaches to go from the exploration goal of “understand the distribution of water” to quantification of specific mission sampling requirements.

Colaprete, A.↗
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↗