Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Avoiding ergodicity problems in lattice discretizations of the Hubbard model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. Here, we show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

First-principles calculations↗
DOE OSTI2025

Twisting the Hubbard model into the momentum-mixing Hatsugai–Kohmoto model

The Hubbard model is a standard theoretical tool for studying materials with strong electron–electron interactions, such as cuprate superconductors. Unfortunately, interaction-driven phenomena, such as a transition into the strongly correlated Mott insulator phase, are difficult to treat with established theoretical techniques. However, the exactly solvable Hatsugai–Kohmoto model displays similar Mott physics. In this work we show how the Hatsugai–Kohmoto model can be deformed continuously into the Hubbard model. The trick is to systematically reintroduce all the momentum mixing that the original Hatsugai–Kohmoto model omits. This can be accomplished by grouping n momenta into a cell and hybridizing them, resulting in the momentum-mixing Hatsugai–Kohmoto model. We recover the Bethe ansatz ground-state energy of the one-dimensional Hubbard model to within 1% from only ten mixed momenta. Overall, the convergence scales as 1/n2 as opposed to the inverse linear behaviour of standard finite-cluster techniques. Our results for a square lattice reproduce all the known features from state-of-the-art simulations also with only a few mixed momenta. Consequently, we believe that the momentum-mixing Hatsugai–Kohmoto model offers an alternative tool for strongly correlated quantum matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Stripe order in the underdoped region of the two-dimensional Hubbard model

Competing inhomogeneous orders are a central feature of correlated electron materials, including the high-temperature superconductors. The two-dimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the latest numerical methods in exhaustive simulations, we uncover the ordering in the underdoped ground state. We find a stripe order that has a highly compressible wavelength on an energy scale of a few kelvin, with wavelength fluctuations coupled to pairing order. The favored filled stripe order is different from that seen in real materials. Our results demonstrate the power of modern numerical methods to solve microscopic models, even in challenging settings.

Zheng, Bo-Xiao↗
DOE OSTI2019

Rotationally invariant slave-boson and density matrix embedding theory: Unified framework and comparative study on the one-dimensional and two-dimensional Hubbard model [Rotationally invariant slave-boson and density matrix embedding theory: A unified framework and a comparative study on the 1D and 2D Hubbard Model]

We present detailed benchmark ground-state calculations of the one- and two-dimensional Hubbard model utilizing the cluster extensions of the rotationally invariant slave-boson mean-field theory and the density matrix embedding theory. Our analysis shows that the overall accuracy and the performance of these two methods are very similar. Furthermore, we propose a unified computational framework that allows us to implement both of these techniques on the same footing. In conclusion, this provides us with a different line of interpretation and paves the ways for developing systematically distinct generalizations of these complementary approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

How correlations change the magnetic structure factor of the kagome Hubbard model

The kagome Hubbard model (KHM) is a paradigmatic example of a frustrated two-dimensional model. While its strongly correlated regime, described by a Heisenberg model, is of topical interest due to its enigmatic prospective spin-liquid ground state, the weakly and moderately correlated regimes remain largely unexplored. Motivated by the rapidly growing number of metallic kagome materials (e.g., Mn 3 Sn, Fe 3 Sn 2 , FeSn, Co 3 Sn 2 S 2 , Gd 3 Ru 4 Al 12 , and AV 3 Sb 5 with A = K, Rb, Cs), here we study the respective regimes of the KHM by means of three complementary numerical methods: the dynamical mean-field theory, the dynamical vertex approximation, and determinant quantum Monte Carlo. In contrast to the archetypal square lattice, we find no tendencies toward magnetic ordering, as magnetic correlations remain short-range. Nevertheless, the magnetic correlations undergo a remarkable crossover as the system approaches the metal-to-insulator transition. The Mott transition itself does not affect the magnetic correlations. Our equal-time and dynamical structure factors can be used as a reference for inelastic neutron scattering experiments on the growing family of metallic kagome materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Kekulé valence bond order in an extended Hubbard model on the honeycomb lattice with possible applications to twisted bilayer graphene

Using large-scale quantum Monte Carlo simulations, we exactly solve a model of fermions hopping on the honeycomb lattice with cluster charge interactions, which has been proposed as an effective model with possible application to twisted bilayer graphene near half-filling. Here, we find an interaction driven semimetal to insulator transition to an insulating phase consisting of a valence bond solid with Kekulé pattern. Finite size scaling reveals that the phase transition of the semimetal to Kekulé valence bond solid phase is continuous and belongs to chiral XY universality class.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Disordered Supersolids in the Extended Bose-Hubbard Model

The extended Bose-Hubbard model captures the essential properties of a wide variety of physical systems including ultracold atoms and molecules in optical lattices, Josephson junction arrays, and certain narrow band superconductors. It exhibits a rich phase diagram including a supersolid phase where a lattice solid coexists with a superfluid. We use quantum Monte Carlo to study the supersolid part of the phase diagram of the extended Bose-Hubbard model on the simple cubic lattice. We add disorder to the extended Bose-Hubbard model and find that the maximum critical temperature for the supersolid phase tends to be suppressed by disorder. But we also find a narrow parameter window in which the supersolid critical temperature is enhanced by disorder. Our results show that supersolids survive a moderate amount of spatial disorder and thermal fluctuations in the simple cubic lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Publisher Correction: Disordered Supersolids in the Extended Bose-Hubbard Model

The extended Bose-Hubbard model captures the essential properties of a wide variety of physical systems including ultracold atoms and molecules in optical lattices, Josephson junction arrays, and certain narrow band superconductors. It exhibits a rich phase diagram including a supersolid phase where a lattice solid coexists with a superfluid. We use quantum Monte Carlo to study the supersolid part of the phase diagram of the extended Bose-Hubbard model on the simple cubic lattice. We add disorder to the extended Bose-Hubbard model and find that the maximum critical temperature for the supersolid phase tends to be suppressed by disorder. But we also find a narrow parameter window in which the supersolid critical temperature is enhanced by disorder. Our results show that supersolids survive a moderate amount of spatial disorder and thermal fluctuations in the simple cubic lattice.

Lin, Fei↗
DOE OSTI2026

Quantum utility in simulating the real-time dynamics of the Fermi–Hubbard model using superconducting quantum computers

The Fermi–Hubbard model is a fundamental model in condensed matter physics that describes strongly correlated electrons. On the other hand, quantum computers are emerging as powerful tools for exploring the complex dynamics of these quantum many-body systems. In this work, we demonstrate the quantum simulation of the one-dimensional Fermi–Hubbard model using IBM's superconducting quantum computers, employing over 100 qubits. We introduce a first-order Trotterization scheme and extend it to an optimized second-order Trotterization for the time evolution in the Fermi–Hubbard model, specifically tailored for the limited qubit connectivity of quantum architectures, such as IBM's platforms. Notably, both Trotterization approaches are scalable and maintain a constant circuit depth at each Trotter step, regardless of the qubit count, enabling us to precisely investigate the relaxation dynamics in the Fermi–Hubbard model by measuring the expectation value of the Néel observable (staggered magnetization) for time-evolved quantum states. Lastly, our successful measurement of expectation values in such large-scale quantum many-body systems, especially at longer time scales with larger entanglement, highlights the quantum utility of superconducting quantum platforms over conventional classical approximation methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effective hopping in holographic Bose and Fermi-Hubbard models

In this paper, we analyze a proposed gravity dual to a SU(N) Bose-Hubbard model, as well as construct a holographic dual of a SU(N) Fermi-Hubbard model from D-branes in string theory. In both cases, the SU(N) is dynamical, i.e. the hopping degrees of freedom are strongly coupled to SU(N) gauge bosons which themselves are strongly interacting. The vacuum expectation value (VEV) of the hopping term (i.e. the hopping energy) is analyzed in the gravity dual as a function of the bulk mass of the field dual to the hopping term, as well as of the coupling constants of the model. The bulk mass controls the anomalous dimension (i.e. the critical exponent) of the hopping term in the SU(N) Bose-Hubbard model. We compare the hopping energy to the corresponding result in a numerical simulation of the ungauged SU(N) Bose-Hubbard model. We find agreement when the hopping parameter is smaller than the other couplings. Our analysis shows that the kinetic energy increases as the bulk mass increases, due to increased contributions from the IR. The holographic Bose-Hubbard model is then compared with the string theory construction of a SU(N) Fermi-Hubbard model. The string theory construction makes it possible to describe fluctuations around a half-filled state in the supergravity limit, which map to O(1) occupation number fluctuations in the Fermi-Hubbard model at half filling. Finally, the VEV of the Bose-Hubbard model is shown to agree with the one of the fermionic Hubbard model with the help of a two-site version of the Jordan-Wigner transformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Effects of the next-nearest-neighbor hopping on the low-dimensional Hubbard model: ferromagnetism, antiferromagnetism, and superconductivity

The Hubbard model has attracted considerable interest due to its prototypical role in describing strongly interacting electronic systems, such as high-critical-temperature superconductors as well as many novel quantum materials. By introducing next-nearest-neighbor (NNN) hoppings to the Hubbard model, the phase diagram becomes richer, and fascinating phenomena arise in both, one-dimensional chains and square lattices, such as: antiferromagnetism, ferromagnetism, superconductivity (SC), as well as charge orders, among others. Moreover, NNN hoppings play a fundamental role in understanding effects of doping on magnetism and pairing orders in strongly interacting regimes. In this article, we review the recent progress in understanding the different competing phases of this model in one and two dimensions from a computational perspective. In conclusion, we comment on the pressing technical challenges, illustrate the controversial results concerning the emergence of the SC phase, and conclude with our perspectives on future explorations.

Hubbard model↗
DOE OSTI2022

Estimation of biquadratic and bicubic Heisenberg effective couplings from multiorbital Hubbard models

Abstract We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion U , from weak to strong coupling, and multiple ratios of the Hund coupling J H to U . The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large U , we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a S = 1 Heisenberg model where by symmetry both nearest-neighbor ( S i ⋅ S j ) and ( S i ⋅ S j ) 2 are allowed, with respective couplings J 1 and J 2 . For the case of three orbitals per site, the mapping is into a S = 3/2 Heisenberg model with ( S i ⋅ S j ), ( S i ⋅ S j ) 2 , and ( S i ⋅ S j ) 3 terms, and respective couplings J 1 , J 2 , and J 3 . The strength of these coupling constants in the Heisenberg models depend on the U , J H , and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J 2 / J 1 and J 3 / J 1 can be. We show that those ratios appear rather limited and, as a qualitative guidance, we conclude that J 2 / J 1 is less than 0.4 and J 3 / J 1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J 2 / J 1 and J 3 / J 1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Superconductivity in the doped Hubbard model and its interplay with next-nearest hopping t ′

The Hubbard model is widely believed to contain the essential ingredients of high-temperature superconductivity. However, proving definitively that the model supports superconductivity is challenging. Here, we report a large-scale density matrix renormalization group study of the lightly doped Hubbard model on four-leg cylinders at hole doping concentration δ = 12.5%. We reveal a delicate interplay between superconductivity and charge density wave and spin density wave orders tunable via next-nearest neighbor hopping t ′. For finite t ′, the ground state is consistent with a Luther-Emery liquid with power-law superconducting and charge density wave correlations associated with half-filled charge stripes. In contrast, for t ′ = 0, superconducting correlations fall off exponentially, whereas charge density and spin density modulations are dominant. Our results indicate that a route to robust long-range superconductivity involves destabilizing insulating charge stripes in the doped Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Finite-temperature tensor network study of the Hubbard model on an infinite square lattice

Here, the Hubbard model is a longstanding problem in the theory of strongly correlated electrons and a very active one in the experiments with ultracold fermionic atoms. Motivated by current and prospective quantum simulations, we apply a two-dimensional tensor network—an infinite projected entangled pair state—evolved in imaginary time by the neighborhood tensor update algorithm working directly in the thermodynamic limit. With U⁡(1)×U⁡(1) symmetry and the bond dimensions up to 29, we generate thermal states down to the temperature of 0.17 times the hopping rate. We obtain results for spin and charge correlators, unaffected by boundary effects. The spin correlators—measurable in prospective ultracold atoms experiments attempting to approach the thermodynamic limit—provide evidence of disruption of the antiferromagnetic background with mobile holes in a slightly doped Hubbard model. The charge correlators reveal the presence of hole-doublon pairs near half filling and signatures of hole-hole repulsion on doping. We also obtain specific heat in the slightly doped regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗