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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermodynamics and magnetism in the two-dimensional to three-dimensional crossover of the Hubbard model

The realization of antiferromagnetic (AF) correlations in ultracold fermionic atoms on an optical lattice is a significant achievement. Experiments have been carried out in one, two, and three dimensions, and have also studied anisotropic configurations with stronger tunneling in some lattice directions. Such anisotropy is relevant to the physics of cuprate superconductors and other strongly correlated materials. Moreover, this anisotropy might be harnessed to enhance AF order. Here we numerically investigate, using the determinant quantum Monte Carlo method, a simple realization of anisotropy in the three-dimensional (3D) Hubbard model in which the tunneling between planes, $t_⊥$, is unequal to the intraplane tunneling $\textit{t}$. This model interpolates between the three-dimensional isotropic ($t_⊥ = t$) and two-dimensional (2D;$t_⊥$ = 0) systems. In this work, we show that at fixed interaction strength to tunneling ratio ($\textit{U/t}$), anisotropy can enhance the magnetic structure factor relative to both 2D and 3D results. However, this enhancement occurs at interaction strengths below those for which the Néel temperature $T_{\text{Néel}}$ is largest, in such a way that the structure factor cannot be made to exceed its value in isotropic 3D systems at the optimal $\textit{U/t}$. We characterize the 2D-3D crossover in terms of the magnetic structure factor, real space spin correlations, number of doubly occupied sites, and thermodynamic observables. An interesting implication of our results stems from the entropy's dependence on anisotropy. As the system evolves from 3D to 2D, the entropy at a fixed temperature increases. Correspondingly, at fixed entropy, the temperature will decrease going from 3D to 2D. This suggests a cooling protocol in which the dimensionality is adiabatically changed from 3D to 2D.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Spin-triplet pairing induced by near-neighbor attraction in the extended Hubbard model for cuprate chain

Abstract In quantum materials, the electronic interaction and the electron-phonon coupling are, in general, two essential ingredients, the combined impact of which may drive exotic phases. Recently, an anomalously strong electron-electron attraction, likely mediated by phonons, has been proposed in one-dimensional copper-oxide chain Ba 2− x Sr x CuO 3+ δ . Yet, it is unclear how this strong near-neighbor attraction V influences the superconductivity pairing in the system. Here we perform accurate many-body calculations to study the extended Hubbard model with on-site Coulomb repulsion U > 0 and near-neighbor attraction V < 0 that could well describe the cuprate chain and likely other similar transition-metal materials with both strong correlations and lattice effects. We find a rich quantum phase diagram containing an intriguing Tomonaga-Luttinger liquid phase — besides the spin density wave and various phase separation phases — that can host dominant spin-triplet pairing correlations and divergent superconductive susceptibility. Upon doping, the spin-triplet superconducting regime can be further broadened, offering a feasible mechanism to realize p -wave superconductivity in realistic cuprate chains.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tendencies of enhanced electronic nematicity in the Hubbard model and a comparison with Raman scattering on high-temperature superconductors

The pseudogap regime of the cuprate high-temperature superconductors is characterized by a variety of competing orders, the nature of which are still widely debated. Recent experiments have provided evidence for electron nematic order, in which the electron fluid breaks rotational symmetry while preserving translational invariance. Raman spectroscopy, with its ability to symmetry resolve low energy excitations, is a unique tool that can be used to assess nematic fluctuations and nematic ordering tendencies. Here we compare results from determinant quantum Monte Carlo simulations of the Hubbard model to experimental results from Raman spectroscopy in La 2-x Sr x CuO 4 , which show a prominent increase in the B 1g response around 10% hole doping as the temperature decreases, indicative of a rise in nematic fluctuations at low energy. Finally, our results support a picture of nematic fluctuations with B 1g symmetry occurring in underdoped cuprates, which may arise from melted stripes at elevated temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Measuring the Loschmidt Amplitude for Finite-Energy Properties of the Fermi-Hubbard Model on an Ion-Trap Quantum Computer

Calculating the equilibrium properties of condensed-matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude ⟨ ψ | e − i H ^ t | ψ ⟩ from initial states | ψ ⟩ and a time evolution under the Hamiltonian H ^ up to short times t . In this work, we study the operation of this algorithm on a present-day quantum computer. Specifically, we measure the Loschmidt amplitude for the Fermi-Hubbard model on a 16 -site ladder geometry (32 orbitals) on the Quantinuum H2-1 trapped-ion device. We assess the effect of noise on the Loschmidt amplitude and implement algorithm-specific error-mitigation techniques. By using a thus-motivated error model, we numerically analyze the influence of noise on the full operation of the quantum-classical algorithm by measuring expectation values of local observables at finite energies. Finally, we estimate the resources needed for scaling up the algorithm. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Ground-state properties of the Hubbard model in one and two dimensions from the Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Off-shell effective energy theory: A unified treatment of the Hubbard model from $\textit{d} = 1$ to $\textit{d = ∞}$

Here we propose an exact formalism, off-shell effective energy theory (OET), which provides a thermodynamic description of a generic quantum Hamiltonian. The OET is based on a partitioning of the Hamiltonian and a corresponding density matrix ansatz constructed from an off-shell extension of the equilibrium density matrix; and there are dual realizations based on a given partitioning. To approximate OET, we introduce the central point expansion (CPE), which is an expansion of the density matrix ansatz, and we renormalize the CPE using a standard expansion of the ground-state energy. In this work we showcase the OET for the one-band Hubbard model in $\textit{d}$ = 1 , 2, and $\infty$, using a partitioning between kinetic and potential energy, yielding two realizations denoted as $\mathcal{K}$ and $\mathcal{X}$. OET shows favorable agreement with exact or state-of-the-art results over all parameter space, and has a negligible computational cost. Physically, $\mathcal{K}$ describes the Fermi liquid, while $\mathcal{X}$ gives an analogous description of both the Luttinger liquid and the Mott insulator. Our approach should find broad applicability in lattice model Hamiltonians, in addition to real materials systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Experimental evaluation of the stabilization of the COT orbitals by 4f orbitals in COT 2 Ce using a Hubbard model

Significant orbital mixing is rare in lanthanide complexes because of the limited radial extent of the 4f orbitals, which results in a generally small stabilization due to 4f orbital interactions. Nevertheless, even a small amount of additional stabilization could enhance lanthanide separations. One lanthanide complex in which orbital mixing has been extensively studied both experimentally and computationally is cerocene, COT 2 Ce, where COT is cyclooctatetraene. This compound has a singlet ground state with a low-lying, triplet excited state. Previous fluorescence studies on trimethylsilyl-substituted cerocenes indicate the triplet state is 0.4 eV higher in energy than the singlet state. In addition, computational studies predict that the triplet is 0.3 to 1 eV higher in energy than the singlet. The synthesis of highly pure COT 2 Ce by Walter and Andersen allowed its physical properties to be accurately measured. Using these measurements, we evaluate the stabilization of the 4f orbitals using two, independent approaches. A Hubbard model is used to evaluate the stabilization of the ground state due to orbital mixing. This stabilization, which is also the singlet–triplet gap, is –0.29 eV using this model. Furthermore, this gap was also from the temperature independent paramagnetism of COT 2 Ce, which yielded a value of –0.32 eV.

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

Interaction-enhanced nesting in spin-fermion and Fermi-Hubbard models

The spin-fermion (SF) model postulates that the dominant coupling between low-energy fermions in near critical metals is mediated by collective spin fluctuations (paramagnons) peaked at the Néel wave vector, Q N , connecting hot spots on opposite sides of the Fermi surface. It has been argued that strong correlations at hot spots lead to a Fermi surface deformation (FSD) featuring flat regions and increased nesting. This conjecture was confirmed in the perturbative self-consistent calculations when the paramagnon propagator dependence on momentum deviation from Q N is given by χ − 1 ∝ | Δ q | . Using diagrammatic Monte Carlo (diagMC) technique we show that such a dependence holds only at temperatures orders of magnitude smaller than any other energy scale in the problem, indicating that a different mechanism may be at play. Instead, we find that a χ − 1 ∝ | Δ q | 2 dependence yields a robust finite- T scenario for achieving FSD. To link phenomenological and microscopic descriptions, we applied the connected determinant diagMC method to the ( t − t ′ ) Hubbard model and found that at large U / t > 5.5 before the formation of electron and hole pockets (i) the FSD defined as a maximum of the spectral function is not very pronounced; instead, it is the lines of zeros of the renormalized dispersion relation that deforms toward nesting, and (ii) the static spin susceptibility is well described by χ − 1 ∝ | Δ q | 2 . Flat FS regions yield a nontrivial scenario for realizing a non-Fermi liquid state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Constructing Hubbard models for the hydrogen chain using sliced-basis density matrix renormalization group

Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and here we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dipolar Bose-Hubbard model

We study a simple model of interacting bosons on a d-dimensional cubic lattice whose dynamics conserves both total boson number and total boson dipole moment. This model provides a simple framework in which several remarkable consequences of dipole conservation can be explored. As a function of chemical potential and hopping strength, the model can be tuned between gapped Mott insulating phases and various types of gapless condensates. The condensed phase realized at large hopping strengths, which we dub a Bose-Einstein insulator, is particularly interesting: despite having a Bose condensate, it is insulating, and despite being an insulator, it is compressible.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Attractive Su-Schrieffer-Heeger-Hubbard model on a square lattice away from half-filling

The Su-Schrieffer-Heeger (SSH) model, with bond phonons modulating electron tunneling, is a paradigmatic electron-phonon model that hosts an antiferromagnetic order to bond order transition at half-filling. In the presence of a repulsive Hubbard interaction, the antiferromagnetic phase is enhanced, but the phase transition remains first order. Here, we explore the physics of the SSH model with an attractive Hubbard interaction, which hosts an interesting interplay among charge order, s-wave pairing, and bond order. Using the numerically exact determinant quantum Monte Carlo method, we show that both charge order, present at weak electron-phonon coupling, and bond order, at large coupling, give way to s-wave pairing when the system is doped. Lastly, we demonstrate that the SSH electron-phonon interaction competes with the attractive Hubbard interaction and reduces the s-wave pairing correlation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Sublattice modulated superconductivity in the kagome Hubbard model

Here we identify a superconducting order featuring spatial pair modulations on the kagome lattice subject to onsite Hubbard U and nearest neighbor V interactions. Within our functional renormalization group analysis, this state appears with a concomitant d-wave superconducting (SC) instability at zero lattice momentum, where it distinguishes itself through intra-unit cell modulations of the pairing function thus breaking the discrete space group symmetry. The relative weight of the sublattice modulated superconductor (SMS) and d-wave SC is influenced by the absolute interaction strength and coupling ratio V/U. Parametrically adjacent to this domain at weak coupling, we find an intra-unit cell modulated vestigial charge density wave and an s-wave SC instability. Our study provides a microscopic setting and thorough description of this novel SMS arising within a translation symmetry broken background.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Enhanced Pair-Density-Wave Vertices in a Bilayer Hubbard Model at Half Filling

Motivated by the pair-density-wave (PDW) state found in the one-dimensional Kondo-Heisenberg chain, we report on a determinant quantum Monte Carlo study of pair fields for a two-dimensional half-filled Hubbard layer coupled to an itinerant, noninteracting layer with one electron per site. In a specific range of interlayer hopping, the pairing vertex associated with PDW order becomes more attractive than that for uniform d-wave pairing, although both remain subdominant to the leading antiferromagnetic correlations at half filling. In conclusion, our result sheds light on where one potentially may find a PDW state in such a model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗