Search NASA⌕ Search

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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Charge correlations suppress unconventional pairing in the Holstein model

In a recent work by Schrodi et al. [Phys. Rev. B. 104, L140506 (2021)], the authors find an unconventional superconducting state with a sign-changing order parameter using the Migdal-Eliashberg theory, including the first vertex correction. This unconventional solution arises despite using an isotropic bare electron-phonon coupling in the Hamiltonian. We examine this claim using hybrid quantum Monte Carlo for a single-band Holstein model with a cuprate-like noninteracting band structure and identical parameters to Schrodi et al.. Our Monte Carlo results for these parameters suggest that unconventional pairing correlations do not exceed their noninteracting values at any carrier concentration we have checked. Instead, strong charge-density-wave correlations persist at the lowest accessible temperatures for dilute and nearly half-filled bands. As a result, we present arguments for how vertex-corrected Migdal-Eliashberg calculation schemes can lead to uncontrolled results in the presence of Fermi surface nesting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multistep Coulomb excitation of Ni 64 : Shape coexistence and nature of low-spin excitations

Here, the structure of 64 Ni, the heaviest stable Ni isotope, has been investigated via high-statistics, multistep safe Coulomb excitation to search for shape coexistence, a phenomenon recently observed in neutron-rich 66 Ni and 70 Ni as well as in doubly magic, N = 40, 68 Ni. The study was motivated by recent, state-of-the-art Monte Carlo shell-model calculations (MCSM), where a Hamiltonian with effective interactions incorporating the monopole tensor force predicts the existence of shape coexistence, also in the lower-mass 62,64 Ni isotopes. A set of transition and static E2 matrix elements for both yrast and near-yrast structures was extracted from the differential Coulomb excitation cross sections. From comparisons between the new results and MCSM as well as other shell-model calculations, a clearer picture of the structure of 64 Ni emerges. Specifically, the low-spin states are shown to be dominated by proton and neutron excitations mainly within the fp shell, with minimal contribution from the g 9/2 shape-driving neutron orbital. The agreement between experimental data and MCSM results indicates a small oblate deformation for the $0$$^{+}_{2}$ level and a spherical shape for the $0$$^{+}_{3}$ state. In addition, the small upper limit determined for the B(E2) probability of a transition associated with the decay of the recently observed 3463-keV, $0$$^{+}_{4}$ state agrees with its proposed assignment to a prolate shape, herewith providing first evidence for triple shape coexistence in a stable Ni isotope.

59 ≤ A ≤ 89↗
DOE OSTI2021

From trivial to topological paramagnets: The case of Z 2 and Z 2 3 symmetries in two dimensions

Using quantum Monte Carlo simulations, we map out the phase diagram of Hamiltonians interpolating between trivial and nontrivial bosonic symmetry-protected topological phases, protected by Z 2 and Z 2 3 symmetries, in two dimensions. In all cases, we find that the trivial and the topological phases are separated by an intermediate phase in which the protecting symmetry is spontaneously broken. Depending on the model, we identify a variety of magnetic orders on the triangular lattice, including ferromagnetism, 3 × 3 order, and stripe orders (both commensurate and incommensurate). Critical properties are determined through a finite-size scaling analysis. Possible scenarios regarding the nature of the phase transitions are discussed.

2-dimensional systems↗
DOE OSTI2020

Charge density waves on a half-filled decorated honeycomb lattice

Tight binding models like the Hubbard Hamiltonian are most often explored in the context of uniform intersite hopping $\textit{t}$. The electron-electron interactions, if sufficiently large compared to this translationally invariant $\textit{t}$, can give rise to ordered magnetic phases and Mott insulator transitions, especially at commensurate filling. The more complex situation of nonuniform $\textit{t}$ has been studied within a number of situations, perhaps most prominently in multiband geometries where there is a natural distinction of hopping between orbitals of different degree of overlap. In this paper we explore related questions arising from the interplay of multiple kinetic energy scales and electron-phonon interactions. Specifically, we use determinant quantum Monte Carlo (DQMC) to solve the half-filled Holstein Hamiltonian on a “decorated honeycomb lattice,” consisting of hexagons with internal hopping $\textit{t}$ coupled together by $\textit{t'}$. This modulation of the hopping introduces a gap in the Dirac spectrum and affects the nature of the topological phases. Here, we determine the range of $\textit{t/t'}$ values which support a charge density wave phase about the Dirac point of uniform hopping $\textit{t = t'}$, as well as the critical transition temperature $T_c$. The QMC simulations are compared with the results of mean field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗
DOE OSTI2022

Phase diagram of the Su-Schrieffer-Heeger-Hubbard model on a square lattice

The Hubbard and Su-Schrieffer-Heeger (SSH) Hamiltonians are iconic models for understanding the qualitative effects of electron-electron and electron-phonon interactions, respectively. In the two-dimensional square lattice Hubbard model at half filling, the on-site Coulomb repulsion U between up and down electrons induces antiferromagnetic (AFM) order and a Mott insulating phase. On the other hand, for the SSH model, there is an AFM phase when the electron-phonon coupling λ is less than a critical value λ c and a bond order wave when λ > λ c . In this Letter, we perform numerical studies on the square lattice optical Su-Schrieffer-Heeger-Hubbard Hamiltonian, which combines both interactions. We use the determinant quantum Monte Carlo method which does not suffer from the fermionic sign problem at half filling. We map out the phase diagram and find that it exhibits a direct first-order transition between an antiferromagnetic phase and a bond-ordered wave as λ increases. Here, the AFM phase is characterized by two different regions. At smaller λ the behavior is similar to that of the pure Hubbard model; the other region, while maintaining long-range AFM order, exhibits larger kinetic energies and double occupancy, i.e., larger quantum fluctuations, similar to the AFM phase found in the pure SSH model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Magnetic, charge, and bond order in the two-dimensional Su-Schrieffer-Heeger-Holstein model

Most nonperturbative numerical studies of electron-phonon interactions focus on model Hamiltonians where the electrons interact with a phonon branch via a single type of microscopic mechanism. Two commonly explored couplings in this context are the Holstein and Su-Schrieffer-Heeger (SSH) interactions, which describe phonons modulating the on-site energy and intersite electron hopping, respectively. Many materials, however, have multiple phonon branches that can each interact with electronic degrees of freedom in different ways. We present here a determinant quantum Monte Carlo study of the half-filled two-dimensional (bond) SSH-Holstein Hamiltonian, where electrons couple to different phonon branches via either the Holstein or SSH mechanism. As a result, we map the model's phase diagram and determine the nature of the transitions between charge-density wave, bond-order wave, and antiferromagnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Evidence for deconfined U(1) gauge theory at the transition between toric code and double semion

Building on quantum Monte Carlo simulations, we study the phase diagram of a one-parameter Hamiltonian interpolating between trivial and topological Ising paramagnets in two dimensions, which are dual to the toric code and the double semion. We discover an intermediate phase with stripe order which spontaneously breaks the protecting Ising symmetry. Remarkably, we find evidence that this intervening phase is gapless due to the incommensurability of the stripe pattern and that it is dual to a $U(1)$ gauge theory exhibiting Cantor deconfinement.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Inhomogeneous mean-field approach to collective excitations near the superfluid–Mott glass transition

Highlights: • Unconventional collective excitations across the superfluid–Mott glass transition. • Higgs mode localized in superfluid and insulating phases. • Broad uncritical amplitude response function at quantum phase transition. • Goldstone mode delocalizes in the superfluid phase. We develop an inhomogeneous quantum mean-field approach to the behavior of collective excitations across the superfluid–Mott glass quantum phase transition in two dimensions, complementing recent quantum Monte Carlo simulations (Puschmann et al. 2020). In quadratic (Gaussian) approximation, the Goldstone (phase) and Higgs (amplitude) modes completely decouple. Each is described by a disordered Bogoliubov Hamiltonian which can be solved by an inhomogeneous multi-mode Bogoliubov transformation. We find that the Higgs mode is spatially localized in both phases. The corresponding scalar spectral function shows a broad peak that is noncritical in the sense that its peak frequency does not soften but remains nonzero across the quantum phase transition. In contrast, the lowest-energy Goldstone mode delocalizes in the superfluid phase, leading to a zero-frequency spectral peak. We compare these findings to the results of the quantum Monte Carlo simulations. We also relate them to general results on the localization of bosonic excitations, and we discuss the limits and generality of our approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗
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

Contributing to Sunny open source code

The Sunny package is an existing open source software to support the simulation of magnetic spin systems. It is available online under the MIT license here: https://github.com/SunnySuite/Sunny.jl The package aims to include the following features: - Analysis for symmetry-allowed interactions - Monte Carlo algorithms for simulating equilibrium configurations. - Schemes for integrating spin dynamics. - Inference of model Hamiltonians - Support for generalized SU(N) spins - Tools for interpreting experimental data, such as diffraction data, magnetization curves, and heat capacity data

Barros, Kipton↗
DOE OSTI2024

Evaluating a quantum-classical quantum Monte Carlo algorithm with Matchgate shadows

Solving the electronic structure problem of molecules and solids to high accuracy is a major challenge in quantum chemistry and condensed matter physics. The rapid emergence and development of quantum computers offer a promising route to systematically tackle this problem. Recent work by [Huggins et al ., Nature (London) 603 , 416 (2022)] proposed a hybrid quantum-classical quantum Monte Carlo (QC-QMC) algorithm using Clifford shadows to determine the ground state of a Fermionic Hamiltonian. This approach displayed inherent noise resilience and the potential for improved accuracy compared to its purely classical counterpart. Nevertheless, the use of Clifford shadows introduces an exponentially scaling postprocessing cost. In this work, we investigate an improved QC-QMC scheme utilizing the recently developed Matchgate shadows technique [Commun. Math. Phys. 404 , 629 (2023)], which removes the aforementioned exponential bottleneck. We observe from experiments on quantum hardware that the use of Matchgate shadows in QC-QMC is inherently noise robust. We show that this noise resilience has a more subtle origin than in the case of Clifford shadows. Nevertheless, we find that classical postprocessing, while asymptotically efficient, requires hours of runtime on thousands of classical CPUs for even the smallest chemical systems, presenting a major challenge to the scalability of the algorithm.

Monte Carlo methods↗
DOE OSTI2026

Constraining Hamiltonians from chiral effective field theory with neutron-star data

Multi-messenger observations of neutron stars (NSs) and their mergers have placed strong constraints on the dense-matter equation of state (EOS). The EOS, in turn, depends on microscopic nuclear interactions that are described by nuclear Hamiltonians. These Hamiltonians are commonly derived within chiral effective field theory (EFT). Ideally, multi-messenger observations of NSs could be used to directly inform our understanding of EFT interactions, but such a direct inference necessitates millions of model evaluations. This is computationally prohibitive because each evaluation requires us to calculate the EOS from a Hamiltonian by solving the quantum many-body problem with methods such as auxiliary-field diffusion Monte Carlo (AFDMC), which provides very accurate and precise solutions but at a significant computational cost. Additionally, we need to solve the stellar structure equations for each EOS which further slows down each model evaluation by a few seconds. In this work, we combine emulators for AFDMC calculations of neutron matter, built using parametric matrix models, and for the stellar structure equations, built using multilayer perceptron neural networks, with the PyCBC data-analysis framework to enable a direct inference of coupling constants in an EFT Hamiltonian using multi-messenger observations of NSs. We find that astrophysical data can provide informative constraints on two-nucleon couplings despite the high densities probed in NS interiors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2024

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗