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 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 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 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 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 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↗
NASA NTRS1989

Numerical simulations - Some results for the 2- and 3-D Hubbard models and a 2-D electron phonon model

Numerical simulations on the half-filled three-dimensional Hubbard model clearly show the onset of Neel order. Simulations of the two-dimensional electron-phonon Holstein model show the competition between the formation of a Peierls-CDW state and a superconducting state. However, the behavior of the partly filled two-dimensional Hubbard model is more difficult to determine. At half-filling, the antiferromagnetic correlations grow as T is reduced. Doping away from half-filling suppresses these correlations, and it is found that there is a weak attractive pairing interaction in the d-wave channel. However, the strength of the pair field susceptibility is weak at the temperatures and lattice sizes that have been simulated, and the nature of the low-temperature state of the nearly half-filled Hubbard model remains open.

Scalapino, D. J.↗
DOE OSTI2025

Unit-density SU(3) Fermi-Hubbard model with spin-flavor imbalance

The advent of ultracold alkaline-earth atoms in optical lattices has established a platform for investigating correlated quantum matter with SU⁡(𝑁) symmetry, offering highly tunable model parameters that allow experiments to access phenomena that are unavailable in conventional materials. Understanding the ground-state physics of SU⁡(𝑁) Fermi-Hubbard models away from the Heisenberg limit and from the spin-flavor balanced setting is important, as examining the flavor imbalance reveals new physics in Fermi-Hubbard models and shows how SU⁡(𝑁) phases react to practical experimental imperfections in optical lattices. Here, in this study, mean-field phase diagrams are presented for the unit-density SU(3) Fermi-Hubbard model at two sets of flavor densities, ($\frac{1}{3}$ − 𝛿, $\frac{1}{3}$ + 𝛿, $\frac{1}{3}$) and ($\frac{1}{4}$ − 𝛿, $\frac{1}{4}$ +𝛿, $\frac{1}{2}$), with the flavor imbalance introduced as 𝛿. Novel phases are identified at moderate interaction strengths for both densities, and their robustness is investigated in the presence of flavor imbalance. Furthermore, we provide microscopic explanations of the phases found and their stability. Analysis of thermal ensembles of random mean-field solutions indicate that, at temperatures accessible in state-of-the-art cold atom experiments, some spin orders are hard for conventional scattering or local observable measurements to detect, but can be more accessible with quantum gas microscopy in optical lattice experiments. This work also shows that nesting and Mottness, intertwined in the usual SU(2) Hubbard model in stark contrast to generic materials, can be tuned in the SU(3) model and play distinct roles. The resulting phase diagrams not only deepen our understanding of SU⁡(𝑁) Fermi-Hubbard models but also inform future experimental search for new phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Moiré Kanamori-Hubbard model for transition metal dichalcogenide homobilayers

Ab initio and continuum model studies predicted that the Γ valley transition metal dichalcogenide (TMD) homobilayers could simulate the conventional multiorbital Hubbard model on the moiré honeycomb lattice. Here, in this work, we perform the Wannierization starting from the continuum model and show that a more general moiré Kanamori-Hubbard model emerges, beyond the extensively studied standard multiorbital Hubbard model, which can be used to investigate the many-body physics in the Γ valley TMD homobilayers. Using the unrestricted Hartree-Fock and Lanczos techniques, we study these half-filled multiorbital moiré bands. By constructing the phase diagrams we predict the presence of an antiferromagnetic state and in addition we found unexpected and dominant states, such as a S = 1 ferromagnetic insulator and a charge density wave state. Our theoretical predictions made using this model can be tested in future experiments on the Γ valley TMD homobilayers.

2-dimensional systems↗
DOE OSTI2023

Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model

The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU⁡(2⁢N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Finite-temperature kinetic ferromagnetism in the square-lattice Hubbard model

While the exact phase diagram of the Fermi-Hubbard model remains poorly understood despite decades of progress, nearly 60 years ago Nagaoka proved that a single dopant in an otherwise halffilled Hubbard system can bring about ferromagnetism through kinetic means. The phenomenon was recently observed with ultracold atoms in triangular optical lattices. Here, we explore the kinetic ferromagnetism within the square lattice Hubbard model and its strong-coupling counterpart, the t − J model, at finite temperatures in the thermodynamic limit via numerical linked-cluster expansions. We find evidence of ferromagnetic Nagaoka polarons at dopings up to ∼ 30% away from half filling for a variety of interaction strengths and at temperatures as low as 0.2 of the hopping energy. We map out the boundaries of this phase through analyzing various correlation functions.

Newby, Robin C. [San José State University] (ORCID↗
DOE OSTI2024

Exploring thermal equilibria of the Fermi-Hubbard model with variational quantum algorithms

Here, this study investigates the thermal properties of the repulsive Fermi-Hubbard model with chemical potential using variational quantum algorithms, crucial in comprehending particle behaviour within lattices at heightened temperatures in condensed matter systems. Conventional computational methods encounter challenges, especially in managing chemical potential, prompting exploration into Hamiltonian approaches. Despite the promise of quantum algorithms, their efficacy is hampered by coherence limitations when simulating extended imaginary time evolution sequences. To overcome these constraints, this research focuses on optimizing variational quantum algorithms to probe the thermal properties of the Fermi-Hubbard model. Physics-inspired circuit designs are tailored to alleviate coherence constraints, facilitating a more comprehensive exploration of materials at elevated temperatures. Our study demonstrates the potential of variational algorithms in simulating the thermal properties of the Fermi-Hubbard model while acknowledging limitations stemming from error sources in quantum devices and encountering barren plateaus.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Universal thermodynamics of an SU ( N ) Fermi-Hubbard model

The SU(2) symmetric Fermi-Hubbard model (FHM) plays an essential role in strongly correlated fermionic many-body systems. In the one particle per site and strongly interacting limit U/t $\gg$ 1, it is effectively described by the Heisenberg Hamiltonian. In this limit, enlarging the spin and extending the typical SU(2) symmetry to SU(N) has been predicted to give exotic phases of matter in the ground state, with a complicated dependence on N. This raises the question of what—if any—are the finite-temperature signatures of these phases, especially in the currently experimentally relevant regime near or above the superexchange energy. We explore this question for thermodynamic observables by numerically calculating the thermodynamics of the SU(N) FHM in the two-dimensional square lattice near densities of one particle per site, using determinant quantum Monte Carlo and numerical linked cluster expansion. Interestingly, we find that for temperatures above the superexchange energy, where the correlation length is short, the energy, number of on-site pairs, and kinetic energy are universal functions of N. Although the physics in the regime studied is well beyond what can be captured by low-order high-temperature series, we show that an analytic description of the scaling is possible in terms of only one- and two-site calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
DOE OSTI2023

Enhanced superconductivity by near-neighbor attraction in the doped extended Hubbard model

A recent experiment has unveiled an anomalously strong electron-electron attraction in the one-dimensional copper-oxide chain Ba 2–x Sr x CuO 3+δ . While the effect of the near-neighbor electron attraction V in the one-dimensional extended Hubbard chain has been examined recently, its effect in the Hubbard model beyond the one-dimensional chain remains unclear. Here, in this work, we report a density matrix renormalization group study of the extended Hubbard model on long four-leg cylinders on the square lattice. We find that the near-neighbor electron attraction V can notably enhance the long-distance superconducting correlations while simultaneously suppressing the charge-density-wave (CDW) correlations. Specifically, for a modestly strong electron attraction, the superconducting correlations become dominant over the CDW correlations with a Luttinger exponent K sc ~1 and strong divergent superconducting susceptibility. Our results provide a promising way to realize long-range superconductivity in the doped Hubbard model. The relevance of our numerical results to cuprate materials is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantitative assessment of the universal thermopower in the Hubbard model

As primarily an electronic observable, the room-temperature thermopower S in cuprates provides possibilities for a quantitative assessment of the Hubbard model. Using determinant quantum Monte Carlo, we demonstrate agreement between Hubbard model calculations and experimentally measured room-temperature S across multiple cuprate families, both qualitatively in terms of the doping dependence and quantitatively in terms of magnitude. We observe an upturn in S with decreasing temperatures, which possesses a slope comparable to that observed experimentally in cuprates. From our calculations, the doping at which S changes sign occurs in close proximity to a vanishing temperature dependence of the chemical potential at fixed density. Our results emphasize the importance of interaction effects in the systematic assessment of the thermopower S in cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗