SEARCH · Search NASA

Results for “Hamiltonian Monte Carlo”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Quantum Monte Carlo studies of a trimer scaling function with microscopic two- and three-body interactions

Here, we present an energy scaling function to predict, in a specific range, the energy of bosonic trimers with large scattering lengths and finite range interactions, which is validated by quantum Monte Carlo calculations using microscopic Hamiltonians with two- and three-body potentials. The proposed scaling function depends on the scattering length, effective range, and a reference energy, which we chose as the trimer energy at unitarity. We obtained the scaling function as a limit cycle from the solution of the renormalized zero-range model with effective range corrections. We proposed a simple parametrization of the energy scaling function. Besides the intrinsic interest in theoretical and experimental investigations, this scaling function allows one to probe Efimov physics with only the trimer ground states, which may open opportunities to identify Efimov trimers whenever access to excited states is limited.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Unconventional skyrmions in an interfacial frustrated triangular lattice

Magnetic skyrmions are spin topological textures of potential interest in spintronics-related data storage and processing devices. Here, in this work, we show the emergence of unconventional skyrmions in a geometrically frustrated triangular lattice on an inversion-symmetry-breaking two-dimensional electron gas substrate. Starting with a classical double-exchange mechanism, this generic interface induces exotic skyrmionic and unique noncoplanar magnetic states not observed in the equivalent square lattice interface. We study the model by deriving an effective spin Hamiltonian. Large scale classical Monte Carlo simulations provide a quantitative evidence for the emergence of these exotic magnetic states. We found that these chiral magnetic states exhibit a substantial and nonzero topological Hall conductivity. As potential material candidates, we propose Cr/MoS 2 , Fe/MoS 2 , and Fe/WSe 2 interfaces because they have the requisite underlying triangular lattice structures and large spin-orbit coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Perfect short-range ordered alloy with line-compound-like properties in the ZnSnN 2 :ZnO system

We present a new solid-state material phase which is a disordered solid solution but offers many ordered line-compound features. The emergent physical phenomena are rooted in the perfect short-range order which conserves the local octet rule. We model the dual-sublattice-mixed semiconductor alloy (ZnSnN 2 ) 1-x (ZnO) 2x using first-principles calculations, Monte-Carlo simulations with a model Hamiltonian, and an extension of the regular solution model by incorporating short-range order. We demonstrate that this unique solid solution, occurring at a “magic” composition, can provide an electronically pristine character without disorder-induced charge localization and, therefore, a superior carrier transport similar to ordered phases. Interestingly, this phase shows singularities that are absent in the conventional solid-solution models, such as the regular solution and band-gap bowing model. Thermodynamically, this alloy phase has a sharply reduced enthalpy at its composition (like a line compound), but it still requires the entropy from long-range disorder to be stabilized at experimentally accessible temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

SmoQyDQMC.jl: A flexible implementation of determinant quantum Monte Carlo for Hubbard and electron-phonon interactions

We introduce the SmoQyDQMC.jl package, a Julia implementation of the determinant quantum Monte Carlo algorithm. SmoQyDQMC.jl supports generalized tight-binding Hamiltonians with on-site Hubbard and generalized electron-phonon ( e e -ph) interactions, including non-linear e e -ph coupling and anharmonic lattice potentials. Our implementation uses hybrid Monte Carlo methods with exact forces for sampling the phonon fields, enabling efficient simulation of low-energy phonon branches, including acoustic phonons. The SmoQyDQMC.jl package also uses a flexible scripting interface, allowing users to adapt it to different workflows and interface with other software packages in the Julia ecosystem. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDQMC.jl or installed using the Julia package manager. The online documentation, including examples, can be obtained from our document page at https://smoqysuite.github.io/SmoQyDQMC.jl/stable/.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Codebase release r0.3 for SmoQyDQMC.jl

Cohen-Stead, Benjamin (ORCID:0000000279156280)↗

DOE OSTI2021

Quantum Monte Carlo Simulations of the 2D Su-Schrieffer-Heeger Model

Over the past several years, a new generation of quantum simulations has greatly expanded our understanding of charge density wave phase transitions in Hamiltonians with coupling between local phonon modes and the on-site charge density. A quite different, and interesting, case is one in which the phonons live on the bonds, and hence modulate the electron hopping. This example, described by the Su-Schrieffer-Heeger (SSH) Hamiltonian, has so far only been studied with quantum Monte Carlo in one dimension. Here we present results for the 2D SSH model, show that a bond ordered wave (BOW) insulator is present in the ground state at half filling, and argue that a critical value of the electron-phonon coupling is required for its onset, in contradistinction with the 1D case where BOW exists for any nonzero coupling. We determine the precise nature of the bond ordering pattern, which has hitherto been controversial, and the critical transition temperature, which is associated with a spontaneous breaking of Ζ 4 symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Monte Carlo simulation of order-disorder transition in refractory high entropy alloys: A data-driven approach

High entropy alloys (HEAs) are a series of novel materials that demonstrate many exceptional mechanical properties. To understand the origin of these attractive properties, it is important to investigate the thermodynamics and elucidate the evolution of various chemical phases. In this work, we introduce a machine learning approach to construct the effective Hamiltonian and study the thermodynamics of HEAs through canonical Monte Carlo simulation. The main characteristic of our method is to use pairwise interactions between atoms as the machine learning features and systematically improve the representativeness of the dataset using active learning from the Monte Carlo samples. We find this method produces highly robust and accurate effective Hamiltonians that give less than 0.1 mRy test error for all the three refractory HEAs: MoNbTaW, MoNbTaVW, and MoNbTaTiW. Using replica exchange to speed up the MC simulation, we calculated the specific heats and short-range order parameters in a wide range of temperatures. For all the studied materials, we find there are two major order-disorder transitions occurring respectively at T 1 and T 2 , where T 1 is near room temperature but T 2 is much higher. Further, we demonstrate that the transition at T 1 is caused by W and Nb while the one at T 2 is caused by the other elements. By comparing with experiments, the results demonstrate that the order-disorder transitions have a profound impact on the strength and ductility of HEAs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Emergent Glassy Behavior in a Kagome Rydberg Atom Array

Here, we present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids—as well as the crossover towards a trivial paramagnetic phase—are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Variational calculation of closed-shell nuclei up to A=40 [Slides]

Compute ground state properties of closed-shell nuclei up to A = 40 with quantum Monte Carlo. Analyze the behavior of phenomenological Hamiltonians in medium-heavy nuclei. Study the high-momentum components of the nuclear wave function and potential in-medium modifications of the nucleon form factors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Magnetic ground state and perturbations of the distorted kagome Ising metal TmAgGe

Here, we present the magnetic orders and excitations of the distorted kagome intermetallic magnet TmAgGe. Using neutron single crystal diffraction we identify the propagation vectors k = ($\frac{1}{2}$ 0 0) and k = (0 0 0) and determine the magnetic structures of the zero-field and magnetic field-induced phases for H along the $\mathcal{a}$ and [–110] crystal directions. We determine the experimental magnetic field- temperature (H, T)-phase diagram and reproduce it by Monte Carlo simulations of an effective spin exchange Hamiltonian for one distorted kagome layer. Our model includes a strong axial single-ion anisotropy and significantly smaller exchange couplings, which span up to the third-nearest neighbors within the layer. Single crystal inelastic neutron scattering (INS) measurements reveal an almost flat, only weakly dispersive mode around 7 meV that we use alongside bulk magnetization data to deduce the crystal-electric field (CEF) scheme for the Tm 3+ ions. Random phase approximation (RPA) calculations based on the determined CEF wave functions of the two lowest quasidoublets enable an estimation of the interlayer coupling that is compatible with the experimental INS spectra. No evidence for low-energy spin waves associated to the magnetic order was found, which is consistent with the strongly Ising nature of the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

First-Principles Statistical Mechanics Study of Magnetic Fluctuations and Order–Disorder in the Spinel LiNi 0.5 Mn 1.5 O 4 Cathode

While significant magnetic interactions exist in lithium transition metal oxides, commonly used as Li-ion cathodes, the interplay between magnetic couplings, disorder, and redox processes remains poorly understood. In this work, we focus on the high-voltage spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) cathode as a model system on which to apply a computational framework that uses first principles-based statistical mechanics methods to predict the finite temperature magnetic properties of materials and provide insights into the complex interplay between magnetic and chemical degrees of freedom. Density functional theory calculations on multiple distinct Ni–Mn orderings within the LNMO system, including the ordered ground-state structure (space group P4332), reveal a preference for a ferrimagnetic arrangement of the Ni and Mn sublattices due to strong antiferromagnetic superexchange interactions between neighboring Mn 4+ and Ni 2+ ions and ferromagnetic Mn–Mn and Ni–Ni couplings, as revealed by magnetic cluster expansions. These results are consistent with qualitative predictions using the Goodenough-Kanamori-Anderson rules. Simulations of the finite temperature magnetic properties of LNMO are conducted using Metropolis Monte Carlo. We find that a “semiclassical” Monte Carlo sampling method based on the Heisenberg Hamiltonian accurately predicts experimental magnetic transition temperatures observed in magnetometry measurements. This study highlights the importance of a robust computational toolkit that accurately captures the complex chemomagnetic interactions and predicts finite temperature magnetic behavior to help analyze experimental magnetic and magnetic resonance spectroscopy data acquired ex situ and operando.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Deducing subnanometer cluster size and shape distributions of heterogeneous supported catalysts

Abstract Infrared (IR) spectra of adsorbate vibrational modes are sensitive to adsorbate/metal interactions, accurate, and easily obtainable in-situ or operando. While they are the gold standards for characterizing single-crystals and large nanoparticles, analogous spectra for highly dispersed heterogeneous catalysts consisting of single-atoms and ultra-small clusters are lacking. Here, we combine data-based approaches with physics-driven surrogate models to generate synthetic IR spectra from first-principles. We bypass the vast combinatorial space of clusters by determining viable, low-energy structures using machine-learned Hamiltonians, genetic algorithm optimization, and grand canonical Monte Carlo calculations. We obtain first-principles vibrations on this tractable ensemble and generate single-cluster primary spectra analogous to pure component gas-phase IR spectra. With such spectra as standards, we predict cluster size distributions from computational and experimental data, demonstrated in the case of CO adsorption on Pd/CeO 2 (111) catalysts, and quantify uncertainty using Bayesian Inference. We discuss extensions for characterizing complex materials towards closing the materials gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Floating Block Method for Quantum Monte Carlo Simulations

Quantum Monte Carlo simulations are powerful and versatile tools for the quantum many-body problem. In addition to the usual calculations of energies and eigenstate observables, quantum Monte Carlo simulations can in principle be used to build fast and accurate many-body emulators using eigenvector continuation or design time-dependent Hamiltonians for adiabatic quantum computing. Furthermore, these new applications require something that is missing from the published literature, an efficient quantum Monte Carlo scheme for computing the inner product of ground state eigenvectors corresponding to different Hamiltonians. In this work, we introduce an algorithm called the floating block method, which solves the problem by performing Euclidean time evolution with two different Hamiltonians and interleaving the corresponding time blocks. We use the floating block method and nuclear lattice simulations to build eigenvector continuation emulators for energies of 4 He, 8 Be, 12 C, and 16 O nuclei over a range of local and nonlocal interaction couplings. From the emulator data, we identify the quantum phase transition line from a Bose gas of alpha particles to a nuclear liquid.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

CASM — A software package for first-principles based study of multicomponent crystalline solids

CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Emulators for Scarce and Noisy Data: Application to Auxiliary-Field Diffusion Monte Carlo for Neutron Matter

Understanding the equation of state (EOS) of pure neutron matter is necessary for interpreting multimessenger observations of neutron stars. Reliable data analyses of these observations require well-quantified uncertainties for the EOS input, ideally propagating uncertainties from nuclear interactions directly to the EOS. This, however, requires calculations of the EOS for a prohibitively larger number of nuclear Hamiltonians, solving the nuclear many-body problem for each one. Quantum Monte Carlo methods, such as auxiliary-field diffusion Monte Carlo (AFDMC), provide precise and accurate results for the neutron matter EOS, but they are very computationally expensive, making them unsuitable for the fast evaluations necessary for uncertainty propagation. Here, we employ parametric matrix models to develop fast emulators for AFDMC calculations of neutron matter and use them to directly propagate uncertainties of coupling constants in the Hamiltonian to the EOS. As these uncertainties include estimates of the effective field theory truncation uncertainty, this approach provides robust uncertainty estimates for use in astrophysical data analyses. In conclusion, this Letter will enable novel applications such as using astrophysical observations to put constraints on coupling constants for nuclear interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Quantum Monte Carlo Methods for Astrophysical Applications

In recent years, new astrophysical observations have provided a wealth of exciting input for nuclear physics. For example, the observations of two-solar-mass neutron stars put strong constraints on possible phase transitions to exotic phases in strongly interacting matter at high densities. Furthermore, the recent observation of a neutron-star merger in both the electromagnetic spectrum and gravitational waves has provided compelling evidence that neutron-star mergers are an important site for the production of extremely neutron-rich nuclei within the r-process. In the coming years, an abundance of new observations is expected, which will continue to provide crucial constraints on the nuclear physics of these events. To reliably analyze such astrophysical observations and extract information on nuclear physics, it is very important that a consistent approach to nuclear systems is used. Such an approach consists of a precise and accurate method to solve the nuclear many-body problem in nuclei and nuclear matter, combined with modern nuclear Hamiltonians that allow to estimate the theoretical uncertainties. Quantum Monte Carlo methods are ideally suited for such an approach and have been successfully used to describe atomic nuclei and nuclear matter. In this contribution, I will present a detailed description of Quantum Monte Carlo methods focusing on the application of these methods to astrophysical problems. In particular, I will discuss how to use Quantum Monte Carlo methods to describe nuclear matter of relevance to the physics of neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Impurity lattice Monte Carlo for hypernuclei

Abstract We consider the problem of including $$\varLambda $$ Λ hyperons into the ab initio framework of nuclear lattice effective field theory. In order to avoid large sign oscillations in Monte Carlo simulations, we make use of the fact that the number of hyperons is typically small compared to the number of nucleons in the hypernuclei of interest. This allows us to use the impurity lattice Monte Carlo method, where the minority species of fermions in the full nuclear Hamiltonian is integrated out and treated as a worldline in Euclidean projection time. The majority fermions (nucleons) are treated as explicit degrees of freedom, with their mutual interactions described by auxiliary fields. This is the first application of the impurity lattice Monte Carlo method to systems where the majority particles are interacting. Here, we show how the impurity Monte Carlo method can be applied to compute the binding energies of the light hypernuclei. In this exploratory work we use spin-independent nucleon–nucleon and hyperon–nucleon interactions to test the computational power of the method. We find that the computational effort scales approximately linearly in the number of nucleons. The results are very promising for future studies of larger hypernuclear systems using chiral effective field theory and realistic hyperon–nucleon interactions, as well as applications to other quantum many-body systems.

Frame, Dillon↗

DOE OSTI2022

Investigation of the monopole magneto-chemical potential in spin ices using capacitive torque magnetometry

The single-ion anisotropy and magnetic interactions in spin-ice systems give rise to unusual non-collinear spin textures, such as Pauling states and magnetic monopoles. The effective spin correlation strength (J eff ) determines the relative energies of the different spin-ice states. With this work, we display the capability of capacitive torque magnetometry in characterizing the magneto-chemical potential associated with monopole formation. We build a magnetic phase diagram of Ho 2 Ti 2 O 7 , and show that the magneto-chemical potential depends on the spin sublattice (α or β), i.e., the Pauling state, involved in the transition. Monte Carlo simulations using the dipolar-spin-ice Hamiltonian support our findings of a sublattice-dependent magneto-chemical potential, but the model underestimates the J eff for the β-sublattice. Additional simulations, including next-nearest neighbor interactions (J 2 ), show that long-range exchange terms in the Hamiltonian are needed to describe the measurements. This demonstrates that torque magnetometry provides a sensitive test for J eff and the spin-spin interactions that contribute to it.

36 MATERIALS SCIENCE↗