Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HUBBARD MODEL”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fluctuations analysis of spin susceptibility: Néel ordering revisited in dynamical mean field theory

Here in this paper, we revisit the antiferromagnetic (AF) phase diagram of the single-band three-dimensional half-filled Hubbard model on a simple cubic lattice studied within the dynamical mean field theory (DMFT). Although this problem has been investigated extensively in the literature, a comprehensive understanding of the impact of the different one- and, in particular, two-particle local correlation functions of DMFT on the AF transition temperature is still missing. We have, hence, performed a fluctuation analysis of T N with respect to different local bosonic fluctuations (charge, spin, particle-particle) contained in the two-particle vertex of DMFT. Our results indicate that, beyond weak coupling, the screening of the DMFT vertex by local fluctuations leads to an enhancement of T N with respect to a random phase approximation (RPA) like calculation where this vertex is replaced by the bare interaction. The overall suppression of T N in DMFT with respect to RPA is then solely due to the incoherence introduced by the DMFT self-energy in the one-particle Green's functions. This illustrates the Janus-faced role of the local moment formation in the DMFT solution of the Hubbard model, which leads to completely opposite effects in the one- and two-particle correlation functions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ferromagnetism and doublon localization in a Wannier-Hubbard chain

We derive a “Wannier-Hubbard” model consisting of an array of overlapping atomic orbitals interacting via a local Coulomb interaction. Transforming to an orthogonal Wannier basis set, the resulting Hamiltonian displays long-range hopping and interactions, with new terms such as correlated hopping and ferromagnetic direct exchange, among others. We numerically study the one-dimensional version of the model at half-filling using the density matrix renormalization group method, unveiling a rich phase diagram as a function of the interaction U and the overlap s with metallic and ferromagnetic phases separated by a ferrimagnetic region. Our results indicate a path toward understanding emergent phases under pressure and beyond standard model Hamiltonians.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermodynamics of correlated electrons in a magnetic field

The Hofstadter–Hubbard model captures the physics of strongly correlated electrons in an applied magnetic field, which is relevant to many recent experiments on Moiré materials. Few large-scale, numerically exact simulations exists for this model. In this work, we simulate the Hubbard–Hofstadter model using the determinant quantum Monte Carlo (DQMC) algorithm. We report the field and Hubbard interaction strength dependence of charge compressibility, fermion sign, local moment, magnetic structure factor, and specific heat. The gross structure of magnetic Bloch bands and band gaps determined by the non-interacting Hofstadter spectrum is preserved in the presence of U. Incompressible regions of the phase diagram have improved fermion sign. At half filling and intermediate and larger couplings, a strong orbital magnetic field delocalizes electrons and reduces the effect of Hubbard U on thermodynamic properties of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mitigating the Sign Problem through Basis Rotations

Quantum Monte Carlo simulations of quantum many-body systems are plagued by the Fermion sign problem. The computational complexity of simulating Fermions scales exponentially in the projection time β and system size. The sign problem is basis dependent and an improved basis, for fixed errors, leads to exponentially quicker simulations. Herein, we show how to use sign-free quantum Monte Carlo simulations to optimize over the choice of basis on large two-dimensional systems. We numerically illustrate these techniques decreasing the “badness” of the sign problem by optimizing over single-particle basis rotations on one- and two-dimensional Hubbard systems. We find a generic rotation which improves the average sign of the Hubbard model for a wide range of U and densities for L×4 systems. In one example improvement, the average sign (and hence simulation cost at fixed accuracy) for the 16×4 Hubbard model at U/t=4 and n=0.75 increases by exp[8.64(6)β]. For typical projection times of β ≳ 100, this accelerates such simulation by many orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A random-sampling method as an efficient alternative to variational Monte Carlo for solving Gutzwiller wavefunctions

Abstract We present a random-sampling (RS) method for evaluating expectation values of physical quantities using the variational approach. We demonstrate that the RS method is computationally more efficient than the variational Monte Carlo method using the Gutzwiller wavefunctions applied on single-band Hubbard models as an example. Non-local constraints can also been easily implemented in the current scheme that capture the essential physics in the limit of strong on-site repulsion. In addition, we extend the RS method to study the antiferromagnetic states with multiple variational parameters for 1D and 2D Hubbard models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Optimizing the critical temperature and superfluid density of a metal-superconductor bilayer

A promising path to realizing higher superconducting transition temperatures 𝑇 c is the strategic engineering of artificial heterostructures. For example, quantum materials could, in principle, be coupled with other materials to produce a more robust superconducting state. Here, in this work, we add numerical support to the hypothesis that a strongly interacting superconductor weakened by phase fluctuations can boost its 𝑇 c by hybridizing the system with a metal. Using determinant quantum Monte Carlo, we simulate a two-dimensional bilayer composed of an attractive Hubbard model and a metallic layer in two regimes of the interaction strength −|𝑈|. In the strongly interacting regime, we find that increasing the interlayer hybridization 𝑡 ⊥ results in a nonmonotonic enhancement of 𝑇 c , with an optimal value comparable to the maximum 𝑇 c observed in the single-layer attractive Hubbard model, confirming trends inferred from other approaches. In the intermediate coupling regime, when −|𝑈| is close to the value associated with the maximum 𝑇 c of the single-layer model, increasing 𝑡 ⊥ tends to decrease 𝑇 c , implying that the correlated layer was already optimally tuned. Importantly, we demonstrate that the mechanism behind these trends is related to enhancement in the superfluid stiffness, as was initially proposed by Kivelson [Phys. B: Condens. Matter 318, 61 (2002)].

Zhang, Yutan [Univ. of California, Davis, CA (Unit↗
DOE OSTI2025

Beyond-Hubbard Pairing in a Cuprate Ladder

The Hubbard model is believed to capture the essential physics of cuprate superconductors. However, recent theoretical studies suggest that it fails to reproduce a robust and homogeneous superconducting ground state. Here, using resonant inelastic x-ray scattering and density matrix renormalization group calculations, we show that magnetic excitations in the prototypical cuprate ladder Sr 14 Cu 24 O 41 are inconsistent with those of a simple Hubbard model. The magnetic response of hole carriers, contributing to an emergent branch of spin-flip excitations, is strongly suppressed. This effect is the consequence of strong d -wavelike pairing, enhanced by nearly an order of magnitude through a large nearest-neighbor attractive interaction and persisting up to at least 260 K. The close connection between the physics of cuprate ladders and that of the two-dimensional compounds suggests that such an enhanced hole pairing may be a universal feature of superconducting cuprates. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

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

Disentangling quantum matter with measurements

Measurements destroy entanglement. Building on ideas used to study ‘quantum disentangled liquids,’ we explore the use of this effect to characterize states of matter. In this work, we focus on systems with multiple components, such as charge and spin in a Hubbard model or local moments and conduction electrons in a Kondo lattice model. In such systems, measurements of (a subset of) one of the components can leave behind a quantum state of the other that is easy to understand, for example in terms of scaling of entanglement entropy of subregions. We bound the outcome of this protocol, for any choice of measurement, in terms of more standard information-theoretic quantities. We apply this quantum disentangling protocol to several problems of physical interest, including gapless topological phases, heavy fermions, and scar states in the Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Imaging local discharge cascades for correlated electrons in WS 2 /WSe 2 moiré superlattices

Transition metal dichalcogenide moiré heterostructures provide an ideal platform to explore the extended Hubbard model, where long-range Coulomb interactions play a critical role in determining strongly correlated electron states. This has led to experimental observations of Mott insulator states at half filling as well as a variety of extended Wigner crystal states at different fractional fillings. However, a microscopic understanding of these emerging quantum phases is still lacking. Here we describe a scanning tunnelling microscopy (STM) technique for the local sensing and manipulation of correlated electrons in a gated WS 2 /WSe 2 moiré superlattice, which enables the experimental extraction of fundamental extended Hubbard model parameters. We demonstrate that the charge state of the local moiré sites can be imaged by their influence on the STM tunnelling current. In addition to imaging, we are also able to manipulate the charge state of correlated electrons. When we ramp the bias on the STM tip, there is a local discharge cascade of correlated electrons in the moiré superlattice, which allows us to estimate the nearest-neighbour Coulomb interaction. Two-dimensional mapping of the moiré electron charge states also enables us to determine the on-site energy fluctuations at different moiré sites. Furthermore, our technique should be broadly applicable to many semiconductor moiré systems, offering a powerful tool for the microscopic characterization and control of strongly correlated states in moiré superlattices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Doping Dependence of 2-Spinon Excitations in the Doped 1D Cuprate Ba 2 ⁢CuO 3+δ

Recent photoemission experiments on the quasi-one-dimensional Ba-based cuprates suggest that doped holes experience an attractive potential not captured using the simple Hubbard model. This observation has garnered significant attention due to its potential relevance to Cooper pair formation in high-T c cuprate superconductors. To scrutinize this assertion, we examined signatures of such an attractive potential in doped 1D cuprates Ba 2 ⁢CuO 3+δ by measuring the dispersion of the 2-spinon excitations using Cu L 3 -edge resonant inelastic x-ray scattering (RIXS). Upon doping, the 2-spinon excitations appear to weaken, with a shift of the minimal position corresponding to the nesting vector of the Fermi points, q F . Notably, we find that the energy scale of the 2-spinons near the Brillouin zone boundary is substantially softened compared to that predicted by the Hubbard model in one dimension. Finally, such a discrepancy implies missing ingredients, which lends support for the presence of an additional attractive potential between holes.

1-dimensional spin chains↗
DOE OSTI2021

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Location and thermal evolution of the pseudogap due to spin fluctuations

We study pseudogap behavior in a metal near a spin density wave (SDW) instability due to thermal magnetic fluctuations. We consider the t–t' Hubbard model on a square lattice at a finite doping, at intermediate coupling strength, and analyze the thermal evolution of the electron spectral function between a SDW ordered state at low temperatures and a normal Fermi liquid at high temperatures. We argue that for proper description of the pseudogap one needs to sum up infinite series of diagrams for both the fermionic self-energy and the SDW order parameter in the SDW state or the magnetic correlation length in the paramagnetic state. We use the eikonal approach to sum up an infinite series of diagrammatic contributions from thermal fluctuations. Earlier studies found that in the SDW state the spectral function A k (ω) of a hot fermion at a finite T is exponentially small below the energy scale Δ(T), which scales with SDW order and vanishes at the ordering temperature T N , and has a hump at a larger frequency Δ PG , comparable to the zero-temperature SDW gap Δ(T = 0). We argue that the hump, which we associate with the pseudogap, survives in some T range above T N . We show that this range is split by regions of strong and weak pseudogap behavior. In the first region, Δ PG is weakly temperature dependent, despite the fact that it comes from thermal fluctuations. Such a behavior has been seen in numerical studies of the Hubbard model. We show that to obtain it one needs to go beyond the one-loop approximation and sum up the infinite series of diagrams. In the second regime, Δ PG decreases with increasing T and eventually vanishes. Here, we further argue that a magnetic pseudogap at a finite T emerges only if the ground state is magnetically ordered. We present the phase diagram and apply the results to high-T c cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
DOE OSTI2026

Improving Spectral Resolution from Real-time Evolution for Correlated Systems

Abstract The quality of numerically simulated spectra using real-time evolution methods for strongly correlated systems is affected by both the length of simulation time and the system size, limiting resolution in both frequency and momentum. In this work, we propose a computationally cheap, linear autoregressive machine learning-based framework to extend short-time and short-distance results over a wider range. We use the proposed method to extend the lesser Green’s function for both the Hubbard model and the much more computationally challenging Hubbard-extended Holstein model. This technique significantly improves both the frequency and momentum resolution of the single-particle removal spectrum $${\mathcal{A}}(k,\omega )$$ A ( k , ω ) , allowing the observation of otherwise obscured spectral features due to electron-phonon coupling.

Tang, Ta↗
DOE OSTI2023

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗