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

Emergence of p+ip superconductivity in two-dimensional doped Dirac systems

Searching for the $p + ip$ superconducting (SC) state has become a fascinating subject in condensed matter physics recently, as a dream application awaiting in topological quantum computation. Unfortunately, so far there is no universal principle for realizing $p + ip$ in generic solid-state systems. Here we report a theoretical discovery of a $p + ip$ SC ground state (coexisting with ferromagnetic order) in the honeycomb lattice Hubbard model in the extremely strong-coupling limit (e.g., infinite U) at low doping $(δ < 0.2)$, by using both the state-of-art Grassmann tensor product state approach and a continuum quantum field theory approach. Our discovery suggests a mechanism for the $p + ip$ SC state in generic strongly correlated systems based on spin-charge separation and the charge current–current coupling scenario, which opens a door towards experimental realization. The $p + ip$ SC state has an instability towards a potential non-Fermi liquid with a large but finite U . Nevertheless, by applying an in-plane Zeeman field, such a $p + ip$ SC state can be stabilized with finite U in a very wide range of doping. Relevant realistic materials are also proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nonequilibrium optical response of a one-dimensional Mott insulator

We define, compute, and analyze the nonequilibrium differential optical conductivity of the one-dimensional extended Hubbard model at half-filling after applying a pump pulse, using the time-dependent density matrix renormalization group method. The melting of the Mott insulator is accompanied by a suppression of the local magnetic moment and ensuing photogeneration of doublon-holon pairs. The differential optical conductivity reveals (i) mid-gap states related to parity-forbidden optical states and (ii) strong renormalization and hybridization of the excitonic resonance and the absorption band, yielding a Fano resonance. We offer evidence and interpret such a resonance as a signature of nonequilibrium optical excitations resembling excitonic strings, (bi)excitons, and unbound doublon-holon pairs, depending on the magnitude of the intersite Coulomb repulsion. We discuss our results in the context of pump and probe spectroscopy experiments on organic Mott insulators.

1-dimensional systems↗
DOE OSTI2021

Magnetic properties of alternating Hubbard ladders

We explore the Hubbard Hamiltonian on ladders where the number of sites per rung alternates between two and three. These geometries are bipartite with nonequal or equal number of sites on the two sublattices. Thus they share a key feature of the Hubbard model in a class of lattices which Lieb has shown analytically to exhibit long-range ferrimagnetic order while being amenable to powerful numeric approaches developed for quasi-one-dimensional geometries. The density matrix renormalization group (DMRG) method is used to obtain the groundstate properties, e.g., excitation gaps, charge and spin densities as well as their correlation functions at half filling. We show the existence of long-range ferrimagnetic order in the one-dimensional ladder geometries. Our work provides detailed quantitative results which complement the general theorem of Lieb for generalized bipartite lattices. It also addresses the issue of how the alternation between quasi-long-range order and spin liquid behavior for uniform ladders with odd and even numbers of legs might be affected by a regular alternation pattern.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗
DOE OSTI2021

Electronic structure and two-band superconductivity in unconventional high- T c cuprates Ba 2 CuO 3+δ

We report that the recently discovered cuprate superconductor Ba 2 CuO 3+δ exhibits a high T c ≃73 K at δ ≃ 0.2. The polycrystal grown under high pressure has a structure similar to La 2 CuO 4 but with dramatically different lattice parameters due to the CuO 6 octahedron compression. The crystal field in the compressed Ba 2 CuO 4 leads to an inverted Cu 3d e g complex with the d x 2 -y 2 orbital sitting below the d 3z 2 - r 2 and an electronic structure highly unusual compared to the conventional cuprates. We construct a two-orbital Hubbard model for the Cu d 9 state at hole doping x = 2δ and study the orbital-dependent strong correlation and superconductivity. For the undoped case at x = 0 , we found that strong correlation drives an orbital-polarized Mott-insulating state with the spin-1/2 moment of the localized d 3z 2 - r 2 orbital. In contrast to the single-band cuprates where superconductivity is suppressed in the overdoped regime, hole doping the two-orbital Mott insulator leads to orbital-dependent correlations and the robust spin and orbital exchange interactions produce a high-T c antiphase d-wave superconductor even in the heavily doped regime at x = 0.4 . We conjecture that Ba 2 CuO 3+δ realizes mixtures of such heavily hole-doped superconducting Ba 2 CuO 4 and disordered Ba 2 CuO 3 chains in a single-layer or predominately separated bilayer structure. Our findings suggest that unconventional cuprates with liberated orbitals as doped two-band Mott insulators can be a direction for realizing high-T c superconductivity with enhanced transition temperature T c .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Circuit complexity near critical points

Here, we consider the Bose–Hubbard model in two and three spatial dimensions and numerically compute the quantum circuit complexity of the ground state in the Mott insulator and superfluid phases using a mean field approximation with additional quadratic fluctuations. After mapping to a qubit system, the result is given by the complexity associated with a Bogoliubov transformation applied to the reference state taken to be the mean field ground state. In particular, the complexity has peaks at the O(2) critical points where the system can be described by a relativistic quantum field theory. Given that we use a Gaussian approximation, near criticality the numerical results agree with a free field theory calculation. To go beyond the Gaussian approximation we use general scaling arguments that imply that, as we approach the critical point t → t c , there is a non-analytic behavior in the complexity c 2 (t) of the form |c 2 (t) – c 2 (t c )| ~ |t – t c | νd , up to possible logarithmic corrections. Here d is the number of spatial dimensions and ν is the usual critical exponent for the correlation length ξ ~ |t – t c | –ν . As a check, for d = 2 this agrees with the numerical computation if we use the Gaussian critical exponent $v=\frac{1}{2}$. Finally, using AdS/CFT methods, we study higher dimensional examples and confirm this scaling argument with non-Gaussian exponent ν for strongly interacting theories that have a gravity dual.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Metal-insulator transition and magnetism of SU(3) fermions in the square lattice

We study the SU(3) symmetric Fermi-Hubbard model (FHM) in the square lattice at 1/3-filling using numerically exact determinant quantum Monte Carlo and numerical linked-cluster expansion techniques. We present the different regimes of the model in the T–U plane, which are characterized by local and short-range correlations, and capture signatures of the metal-insulator transition and magnetic crossovers. These signatures are detected as the temperature scales characterizing the rise of the compressibility, and an interaction-dependent change in the sign of the diagonal spin-spin correlation function. The analysis of the compressibility estimates the location of the metal-insulator quantum critical point at U c /t ~ 6, and provides a temperature scale for observing Mott physics at finite T. Furthermore, from the analysis of the spin-spin correlation function we observe that for U/t ≳ 6 and T ~ J = 4⁢t 2 /U there is a development of a short-range two-sublattice (2SL) antiferromagnetic structure, as well as an emerging three-sublattice (3SL) antiferromagnetic structure as the temperature is lowered below T/J ≲ 0.57. This crossover from 2SL to 3SL magnetic ordering agrees with Heisenberg limit predictions, and has observable effects on the density of on-site pairs. Finally, we describe how the features of the regimes in the T–U plane can be explored with alkaline-earth-like atoms in optical lattices with currently achieved experimental techniques and temperatures. Furthermore, the results discussed in this paper provide a starting point for the exploration of the SU(3) FHM upon doping.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Spin-orbital liquids and insulator-metal transitions on the pyrochlore lattice

The two orbital Hubbard model, with the electrons additionally coupled to a complex magnetic background, arises in the pyrochlore molybdates. The background involves local moments Hund's coupled to the electrons, driving double exchange ferromagnetism, and antiferromagnetic (AF) tendency arising from competing superexchange. The key scales include the Hubbard repulsion and the superexchange, both of which can be tuned in these materials. They control the phase transition from a ferromagnetic metal to a spin glass metal and then to a spin glass (Mott) insulator. We provide a comprehensive description of the ground state of this model using an unrestricted Hartree-Fock scheme implemented via a simulated annealing procedure and establish the metal-insulator transition line for varying Hubbard interaction and superexchange. The electrons see an effective disorder, due to orbital frustration, already in the ferromagnetic phase. The disorder is further enhanced by antiferromagnetic coupling and the resulting magnetic disorder. As a result, increasing AF coupling shifts the metal-insulator transition to lower Hubbard interaction and gives it an additional "Anderson" character. In conclusion, we provide detailed results on the magnetic and orbital correlations, the density of states, and the optical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Superconductivity, charge density waves, and bipolarons in the Holstein model

The electron-phonon ( e-ph ) interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for e-ph coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless e-ph coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate e-ph coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Doping-driven resistive collapse of the Mott insulator in a minimal model for VO 2

Here we study the resistive collapse of the Mott insulator state in the dimer Hubbard model. This minimal model has been used to describe the physics of VO 2 , and should be relevant to other strongly correlated materials. It incorporates the physics of correlated dimers and the explicit competition between on-site Coulomb repulsion and magnetic exchange interactions. Our results unveil that between the Mott insulator at half filling and the Fermi liquid metal at high doping there is an intermediate bad metallic phase with exotic features such as a pseudogap, orbital selectivity, and a first-order metal-metal transition. The model is solved within dynamical mean field theory by means of quantum Monte Carlo, which provides the numerically exact solution of the model in the limit of large lattice dimensionality. This model can be considered as a minimal one that captures exotic phenomena associated to the physics of a doped Mott insulator, shading light on their basic physical mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tensor renormalization group for fermions

Abstract We review the basic ideas of the tensor renormalization group method and show how they can be applied for lattice field theory models involving relativistic fermions and Grassmann variables in arbitrary dimensions. We discuss recent progress for entanglement filtering, loop optimization, bond-weighting techniques and matrix product decompositions for Grassmann tensor networks. The new methods are tested with two-dimensional Wilson–Majorana fermions and multi-flavor Gross–Neveu models. We show that the methods can also be applied to the fermionic Hubbard model in 1+1 and 2+1 dimensions.

Physics↗
DOE OSTI2024

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantum Criticality of Antiferromagnetism and Superconductivity with Relativity

We study a quantum phase transition from a massless to massive Dirac fermion phase in a new two-dimensional bipartite lattice model of electrons that is amenable to sign-free quantum Monte Carlo simulations. Importantly, interactions in our model are not only invariant under SU(2) symmetries of spin and charge like the Hubbard model, but they also preserve an Ising-like electron spin-charge flip symmetry. From unbiased fermion bag Monte Carlo simulations with up to 2304 sites, we show that the massive fermion phase spontaneously breaks this Ising symmetry, picking either antiferromagnetism or superconductivity, and that the transition at which both orders are simultaneously quantum critical belongs to a new “chiral spin-charge symmetric” universality class. We explain our observations using effective potential and renormalization group calculations within the framework of a continuum field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nuclear two point correlation functions on a quantum computer

The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local, and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this paper we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hund bands in spectra of multiorbital systems

Spectroscopy experiments are routinely used to characterize the behavior of strongly correlated systems. An in-depth understanding of the different spectral features is thus essential. Here, we show that the spectrum of the multiorbital Hubbard model exhibits unique Hund bands that occur at energies given only by the Hund coupling J H , as distinct from the Hubbard satellites following the interaction U. We focus on experimentally relevant single-particle and optical spectra that we calculate for a model related to iron chalcogenide ladders. The calculations are performed via the density-matrix renormalization group and Lanczos methods. Finally, the generality of the implications is verified by considering a generic multiorbital model within dynamical mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A beginner's guide to non-abelian iPEPS for correlated fermions

Infinite projected entangled pair states (iPEPS) have emerged as a powerful tool for studying interacting two-dimensional fermionic systems. In this review, we discuss the iPEPS construction and some basic properties of this tensor network (TN) ansatz. Special focus is put on (i) a gentle introduction of the diagrammatic TN representations forming the basis for deriving the complex numerical algorithm, and (ii) the technical advance of fully exploiting non-abelian symmetries for fermionic iPEPS treatments of multi-band lattice models. The exploitation of non-abelian symmetries substantially increases the performance of the algorithm, enabling the treatment of fermionic systems up to a bond dimension D=24 D = 24 on a square lattice. A variety of complex two-dimensional (2D) models thus become numerically accessible. Here, we present first promising results for two types of multi-band Hubbard models, one with 2 2 bands of spinful fermions of \mathrm{SU}(2)_\mathrm{spin} \otimes \mathrm{SU}(2)_\mathrm{orb} S U ( 2 ) s p i n ⊗ S U ( 2 ) o r b symmetry, the other with 3 3 flavors of spinless fermions of \mathrm{SU}(3)_\mathrm{flavor} S U ( 3 ) f l a v o r symmetry.

Bruognolo, Benedikt↗
DOE OSTI2022

Prediction of orbital-selective Mott phases and block magnetic states in the quasi-one-dimensional iron chain Ce 2 O 2 FeSe 2 under hole and electron doping

The recent detailed study of quasi-one-dimensional iron-based ladders, with the 3d iron electronic density n=6, has unveiled surprises, such as orbital-selective phases. However, similar studies for n=6 iron chains are still rare. Here a three-orbital electronic Hubbard model was constructed to study the magnetic and electronic properties of the quasi-one-dimensional n=6 iron chain Ce 2 O 2 FeSe 2 , with focus on the effect of doping. Specifically, introducing the Hubbard U and Hund J H couplings and studying the model via the density matrix renormalization group, we report the ground-state phase diagram varying the electronic density away from n=6. For the realistic Hund coupling J H /U=1/4, several electronic phases were obtained, including a metal, orbital-selective Mott, and Mott insulating phases. Doping away from the parent phase, the competition of many tendencies leads to a variety of magnetic states, such as ferromagnetism, as well as several antiferromagnetic and magnetic “block” phases. In the hole-doping region, two different interesting orbital-selective Mott phases were found: OSMP1 (with one localized orbital and two itinerant orbitals) and OSMP2 (with two localized orbitals and one itinerant orbital). Moreover, charge disproportionation phenomena were found in special doping regions. We argue that our predictions can be tested by simple modifications in the original chemical formula of Ce 2 O 2 FeSe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1972

Role of phonons in the metal-insulator phase transition.

Review, for the transition series oxides, of the Mattis and Lander model, which is one of electrons interacting with lattice vibrations (electron and phonon interaction). The model displays superconducting, insulating, and metallic phases. Its basic properties evolve from a finite crystallographic distortion associated with a dominant phonon mode and the splitting of the Brillouin zone into two subzones, a property of simple cubic and body centered cubic lattices. The order of the metal-insulator phase transition is examined. The basic model has a second-order phase transition and the effects of additional mechanisms on the model are calculated. The way in which these mechanisms affect the magnetically ordered transition series oxides as described by the Hubbard model is discussed.

Langer, W. D.↗