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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗
DOE OSTI2025

High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo

We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius–Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. Here, we compare with previous simulations and experimental measurements.

08 HYDROGEN↗
DOE OSTI2024

Disentangling the physics of the attractive Hubbard model as a fully interacting model of fermions via the accessible and symmetry-resolved entanglement entropies

The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states like the entanglement entropy can provide a complementary and clarifying lens on electronic behavior. One particularly useful measure that can be used to distinguish between quantum and classical sources of entanglement is the accessible entanglement, the entanglement available as a quantum resource for systems subject to conservation laws, such as fixed particle number, due to superselection rules. In this work, we introduce an algorithm and demonstrate how to compute accessible and symmetry-resolved entanglements for interacting fermion systems. This is accomplished by combining an incremental version of the swap algorithm with a recursive auxiliary field quantum Monte Carlo algorithm recently developed by the authors. We apply these tools to study the pairing and charge density waves exhibited in the paradigmatic attractive Hubbard model via entanglement. We find that the particle and spin symmetry-resolved entanglements and their related full probability distribution functions show very clear—and unique—signatures of the underlying electronic behavior even when those features are less pronounced in conventional correlation functions. Altogether, this work provides a systematic means of characterizing the entanglement within quantum systems that can grant a deeper understanding of the complicated electronic behavior that underlies quantum phase transitions and crossovers in many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2023

Application of the variational autoencoder to detect the critical points of the anisotropic Ising model

We generalize the previous study on the application of variational autoencoders to the two-dimensional Ising model to a system with anisotropy. Due to the self-duality property of the system, the critical points can be located exactly for the entire range of anisotropic coupling. This presents an excellent test bed for the validity of using a variational autoencoder to characterize an anisotropic classical model. Furthermore, we reproduce the phase diagram for a wide range of anisotropic couplings and temperatures via a variational autoencoder without the explicit construction of an order parameter. Considering that the partition function of ($d$ + 1)-dimensional anisotropic models can be mapped to that of the $d$-dimensional quantum spin models, the present study provides numerical evidence that a variational autoencoder can be applied to analyze quantum systems via the quantum Monte Carlo method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Preparation for quantum simulation of the ( 1 + 1 ) -dimensional O(3) nonlinear σ model using cold atoms

The 1+1D O(3) non-linear σ-model is a model system for future quantum lattice simulations of other asymptotically-free theories, such as non-Abelian gauge theories. We find that utilizing dimensional reduction can make efficient use of two-dimensional layouts presently available on cold atom quantum simulators. A new definition of the renormalized coupling is introduced, which is applicable to systems with open boundary conditions and can be measured using analog quantum simulators. Monte Carlo and tensor network calculations are performed to determine the quantum resources required to reproduce perturbative short-distance observables. In particular, we show that a rectangular array of 48 Rydberg atoms with existing quantum hardware capabilities should be able to adiabatically prepare low-energy states of the perturbatively-matched theory. Furthermore, these states can then be used to simulate non-perturbative observables in the continuum limit that lie beyond the reach of classical computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Toward Discord : Code for Simulating Continuous Spin Systems

A new computational tool to simulate classical spin systems with frustrated crystal structures is presented. Complementary single- and cluster-spin flip algorithms are implemented to calculate the diffuse scattering patterns, spin-pair correlations, and thermodynamic quantities. Test cases of geometrically frustrated kagome, pyrochlore, and cubic systems are detailed. Two recent scientific cases are also shown here. This new method, together with recent developments of the rmc-discord package (https://github.com/zjmorgan/rmc-discord), represent integrated and strategic step in a complete forward and reverse Monte Carlo framework discord.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stochastic renormalization group and gradient flow

A non-perturbative and continuous definition of RG transformations as stochastic processes is proposed, inspired by the observation that the functional RG equations for effective Boltzmann factors may be interpreted as Fokker-Planck equations. The result implies a new approach to Monte Carlo RG that is amenable to lattice simulation. Long-distance correlations of the effective theory are shown to approach gradient-flowed correlations, which are simpler to measure. The Markov property of the stochastic RG transformation implies an RG scaling formula which allows for the measurement of anomalous dimensions when transcribed into gradient flow expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Photometry of Outer Solar System Objects from the Dark Energy Survey. I. Photometric Methods, Light-curve Distributions, and Trans-Neptunian Binaries

We report the methods of and initial scientific inferences from the extraction of precision photometric information for the >800 trans-Neptunian objects (TNOs) discovered in the images of the Dark Energy Survey (DES). Scene-modeling photometry is used to obtain shot-noise-limited flux measures for each exposure of each TNO, with background sources subtracted. Comparison of double-source fits to the pixel data with single-source fits are used to identify and characterize two binary TNO systems. A Markov Chain Monte Carlo method samples the joint likelihood of the intrinsic colors of each source as well as the amplitude of its flux variation, given the time series of multiband flux measurements and their uncertainties. A catalog of these colors and light-curve amplitudes A is included with this publication. We show how to assign a likelihood to the distribution q(A) of light-curve amplitudes in any subpopulation. Using this method, we find decisive evidence (i.e., evidence ratio <0.01) that cold classical (CC) TNOs with absolute magnitude 6 < H$_{r}$ < 8.2 are more variable than the hot classical (HC) population of the same H$_{r}$, reinforcing theories that the former form in situ and the latter arise from a different physical population. Resonant and scattering TNOs in this H$_{r}$ range have variability consistent with either the HCs or CCs. DES TNOs with H$_{r}$ < 6 are seen to be decisively less variable than higher-H$_{r}$ members of any dynamical group, as expected. More surprising is that detached TNOs are decisively less variable than scattering TNOs, which requires them to have distinct source regions or some subsequent differential processing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Eliashberg theory of phonon-mediated superconductivity — When it is valid and how it breaks down

We analyze the validity of Eliashberg theory of phonon-mediated superconductivity in 2D systems in light of recent extensive Monte-Carlo studies of the Holstein model. Conventional wisdom says that Eliashberg theory is applicable as long as vertex corrections remain small. For small ratio of the phonon energy Ω 0 and the Fermi energy E F , this condition is supposed to hold even when the dimensionless electron–phonon coupling λ is larger than one, i.e., in the strong coupling regime. A comparison between various quantities computed in the Migdal approximation and those computed by Quantum Monte Carlo prove that this belief is wrong, and we identify analytically some of the ways in which this breakdown occurs for various ‘‘normal state’’ properties at λ = λ cr , where λ cr = O(1). Here, the breakdown occurs at temperatures high enough that neither superconducting nor charge–density wave correlations extend over any significant range of distances, so it cannot be associated with the onset of an instability toward any of the relevant ordered ground-states — rather it is associated with the local physics of classical bipolaron formation. Still, we show that certain properties, including the superconducting Tc and the superconducting gap structure below T c , can be accurately inferred from the strong-coupling limit of Eliashberg theory at λ ≤ λ cr .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

MLPF: efficient machine-learned particle-flow reconstruction using graph neural networks

In general-purpose particle detectors, the particle-flow algorithm may be used to reconstruct a comprehensive particle-level view of the event by combining information from the calorimeters and the trackers, significantly improving the detector resolution for jets and the missing transverse momentum. In view of the planned high-luminosity upgrade of the CERN Large Hadron Collider (LHC), it is necessary to revisit existing reconstruction algorithms and ensure that both the physics and computational performance are sufficient in an environment with many simultaneous proton–proton interactions (pileup). Machine learning may offer a prospect for computationally efficient event reconstruction that is well-suited to heterogeneous computing platforms, while significantly improving the reconstruction quality over rule-based algorithms for granular detectors. We introduce MLPF, a novel, end-to-end trainable, machine-learned particle-flow algorithm based on parallelizable, computationally efficient, and scalable graph neural network optimized using a multi-task objective on simulated events. We report the physics and computational performance of the MLPF algorithm on a Monte Carlo dataset of top quark–antiquark pairs produced in proton–proton collisions in conditions similar to those expected for the high-luminosity LHC. The MLPF algorithm improves the physics response with respect to a rule-based benchmark algorithm and demonstrates computationally scalable particle-flow reconstruction in a high-pileup environment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2022

Phase diagram of a distorted kagome antiferromagnet and application to Y-kapellasite

Abstract We investigate the magnetism of a previously unexplored distorted spin-1/2 kagome model consisting of three symmetry-inequivalent nearest-neighbor antiferromagnetic Heisenberg couplings J ⬡ , J , and $$J^{\prime}$$ J ′ , and uncover a rich ground state phase diagram even at the classical level. Using analytical arguments and numerical techniques we identify a collinear $$\overrightarrow{Q}=0$$ Q → = 0 magnetic phase, two unusual non-collinear coplanar $$\overrightarrow{Q}=(1/3,1/3)$$ Q → = ( 1 / 3 , 1 / 3 ) phases and a classical spin liquid phase with a degenerate manifold of non-coplanar ground states, resembling the jammed spin liquid phase found in the context of a bond-disordered kagome antiferromagnet. We further show with density functional theory calculations that the recently synthesized Y-kapellasite Y 3 Cu 9 (OH) 19 Cl 8 is a realization of this model and predict its ground state to lie in the region of $$\overrightarrow{Q}=(1/3,1/3)$$ Q → = ( 1 / 3 , 1 / 3 ) order, which remains stable even after the inclusion of quantum fluctuation effects within variational Monte Carlo and pseudofermion functional renormalization group. The presented model opens a new direction in the study of kagome antiferromagnets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deconstructing magnetization noise: Degeneracies, phases, and mobile fractionalized excitations in tetris artificial spin ice

Direct detection of spontaneous spin fluctuations, or “magnetization noise,” is emerging as a powerful means of revealing and studying magnetic excitations in both natural and artificial frustrated magnets. Depending on the lattice and nature of the frustration, these excitations can often be described as fractionalized quasiparticles possessing an effective magnetic charge. Here, by combining ultrasensitive optical detection of thermodynamic magnetization noise with Monte Carlo simulations, we reveal emergent regimes of magnetic excitations in artificial “tetris ice.” A marked increase of the intrinsic noise at certain applied magnetic fields heralds the spontaneous proliferation of fractionalized excitations, which can diffuse independently, without cost in energy, along specific quasi-1D spin chains in the tetris ice lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Monolayer CrCl 3 as an Ideal Test Bed for the Universality Classes of 2D Magnetism

The monolayer halides CrX 3 (X=Cl, Br, I) attract significant attention for realizing 2D magnets with genuine long-range order (LRO), challenging the Mermin-Wagner theorem. In this work, we show that monolayer CrCl 3 has the unique benefit of exhibiting tunable magnetic anisotropy upon applying a compressive strain. This opens the possibility to use CrCl 3 for producing and studying both ferromagnetic and antiferromagnetic 2D Ising-type LRO as well as the Berezinskii-Kosterlitz-Thouless (BKT) regime of 2D magnetism with quasi-LRO. Using state-of-the-art density functional theory, we explain how realistic compressive strain could be used to tune the monolayer's magnetic properties so that it could exhibit any of these phases. Building on large-scale quantum Monte Carlo simulations, we compute the phase diagram of strained CrCl 3 , as well as the magnon spectrum with spin-wave theory. Our results highlight the eminent suitability of monolayer CrCl 3 to achieve very high BKT transition temperatures, around 50 K, due to their singular dependence on the weak easy-plane anisotropy of the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Clock model interpolation and symmetry breaking in O(2) models

The $q$-state clock model is a classical spin model that corresponds to the Ising model when $q=2$ and to the $XY$ model when $q\to\infty$. The integer-$q$ clock model has been studied extensively and has been shown to have a single phase transition when $q=2$,$3$,$4$ and two phase transitions when $q>4$.We define an extended $q$-state clock model that reduces to the ordinary $q$-state clock model when $q$ is an integer and otherwise is a continuous interpolation of the clock model to noninteger $q$. We investigate this class of clock models in 2D using Monte Carlo (MC) and tensor renormalization group (TRG) methods, and we find that the model with noninteger $q$ has a crossover and a second-order phase transition. We also define an extended-$O(2)$ model (with a parameter $\gamma$) that reduces to the $XY$ model when $\gamma=0$ and to the extended $q$-state clock model when $\gamma\to\infty$, and we begin to outline the phase diagram of this model. These models with noninteger $q$ serve as a testbed to study symmetry breaking in situations corresponding to quantum simulators where experimental parameters can be tuned continuously.

Hostetler, Leon↗
DOE OSTI2022

Structural Properties of Nanometer‐Sized Gold Nanoparticles on a Silicon Substrate

Properties of gold nanoparticles (Au NPs) supported on Si substrates, which have been the focus of recent attention, depend sensitively on the size, structure, and coordination of Au NPs as well as on the processes of interdiffusion and chemical bonding at the Au–Si interface. There is a great need, therefore, for developing structural techniques capable of imaging such nanosystems in atomic detail. Here, it is discussed how the atomic structure of nanometer‐sized Au NPs deposited on a silicon surface can be determined by using high‐energy X‐ray diffraction measurements in combination with an analysis based on atomic‐pair‐distribution functions and advanced reverse Monte Carlo methods. This approach delivers simple and incisive pictures of the structure of Au clusters on Si and reveals direct structural evidence for the emergence of a Au‐silicide interface following room temperature deposition of Au NPs on Si. This study will contribute to designing Au/Si nanostructures with desirable properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗