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 181 records · Page 10

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↗

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↗

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↗

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↗

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↗

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↗

Metal-insulator transition in transition metal dichalcogenide heterobilayer moiré superlattices

Moiré superlattices formed in two-dimensional semiconductor heterobilayers provide a new realization of Hubbard model physics in which the number of electrons per effective atom can be tuned at will. Here we report on an exact diagonalization study of the electronic properties of half-filled narrow moiré bands in which correlation strengths are varied by changing twist angles or interaction strengths. We construct a phase diagram for the bilayer, identifying where the metal-insulator phase transition occurs, estimating the sizes of the charge gaps in the insulating phase, and commenting on the nature of the transition and the importance of subdominant interaction parameters.

2-dimensional systems↗

π phase shift across stripes in a charge density wave system

Many strongly correlated materials are characterized by deeply intertwined charge and spin order. Besides their high superconducting transition temperatures, one of the central features of these complex patterns in cuprates is a phase shift which occurs across lines of decreased hole density. That is, when doped away from their antiferromagnetic phase, the additional charge is not distributed uniformly, but rather in “stripes.” The sublattices preferentially occupied by up and down spin are reversed across these stripes, a phenomenon referred to as a “π phase shift.” Many of the spin-charge patterns, including the π phase shift, are reproduced by density matrix renormalization group and quantum Monte Carlo calculations of simplified tight binding (repulsive Hubbard) models. In this paper we demonstrate that this sublattice reversal is generic by considering the corresponding phenomenon in the attractive Hubbard Hamiltonian, where a charge density wave phase forms at half filling. Here, we introduce charge stripes via an appropriate local chemical potential; measurements of charge correlation across the resulting lines of lowered density reveal a clear π phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on the group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multilayer fully connected neural network, which is trained from data sets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multiscale modeling of correlated electron systems.

36 MATERIALS SCIENCE↗

Signature of Many-Body Localization of Phonons in Strongly Disordered Superlattices

Many-body localization (MBL) has attracted significant research attention due to its immunity to thermalization, role in logarithmic entanglement entropy growth, and opportunities to reach exotic quantum orders. However, experimental realization of MBL in solid-state systems has remained challenging. Here we report evidence of a possible phonon MBL phase in disordered GaAs/AlAs superlattices. Through grazing-incidence inelastic X-ray scattering, we observe a strong deviation of the phonon population from equilibrium in disordered samples doped with ErAs nanodots at low temperature, signaling a departure from thermalization. This behavior occurs within finite phonon energy and wavevector windows, suggesting a localization-thermalization crossover. We support our observation by proposing a theoretical model for the effective phonon Hamiltonian in disordered superlattices, and showing that it can be mapped exactly to a disordered 1D Bose-Hubbard model with a known MBL phase. The large on-site repulsion qualitatively explains the experimental findings. Furthermore, our work provides mode-resolved experimental evidence of phonon localization, potentially extending the scope of MBL to disordered solid-state systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum scars of bosons with correlated hopping

Abstract Recent experiments on Rydberg atom arrays have found evidence of anomalously slow thermalization and persistent density oscillations, which have been interpreted as a many-body analog of the phenomenon of quantum scars. Periodic dynamics and atypical scarred eigenstates originate from a “hard” kinetic constraint: the neighboring Rydberg atoms cannot be simultaneously excited. Here we propose a realization of quantum many-body scars in a 1D bosonic lattice model with a “soft” constraint in the form of density-assisted hopping. We discuss the relation of this model to the standard Bose-Hubbard model and possible experimental realizations using ultracold atoms. We find that this model exhibits similar phenomenology to the Rydberg atom chain, including weakly entangled eigenstates at high energy densities and the presence of a large number of exact zero energy states, with distinct algebraic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗

Particle-hole asymmetric phases in doped twisted bilayer graphene

Twisted bilayer graphene (TBG) has emerged as a paradigmatic platform for exploring the interplay between strong interactions in a multi-band system with nearly flat bands, while offering unprecedented control over the filling fraction of electron/hole carriers. Despite much theoretical work, developing a comprehensive ab initio model for this system has proven challenging due to the inherent trade-off between accurately describing the band structure and incorporating the interactions within the Hamiltonian. Here, we constructed an extended multi-orbital Hubbard model, and performed Hartree-Fock (HF) calculations to explore its phase diagram across commensurate fillings. We found several nearly-degenerate insulating states at charge neutrality, all of which exhibit orbital orders. We also observed that most of the Hartree–Fock ground states exhibit a generalized valley Hund’s-like rule, resulting in valley polarization. Importantly, we show that the incorporation of the intra-valley and inter-valley exchange interactions is crucial to properly stabilize the ordered symmetry-broken states. In agreement with experiments, we find significant particle-hole asymmetry, which underscores the importance of using particle-hole asymmetric models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hamming distance and the onset of quantum criticality

Simulating models for quantum correlated matter unveils the inherent limitations of deterministic classical computations. In particular, in the case of quantum Monte Carlo methods, this is manifested by the emergence of negative weight configurations in the sampling, that is, the sign problem (SP). There have been several recent calculations which exploit the SP to locate underlying critical behavior. Here, utilizing a metric that quantifies phase-space ergodicity in such sampling, the Hamming distance, we suggest a significant advance on these ideas to extract the location of quantum critical points in various fermionic models, in spite of the presence of a severe SP. Combined with other methods, exact diagonalization in our case, it elucidates both the nature of the different phases as well as their location, as we demonstrate explicitly for the honeycomb and triangular Hubbard models, in both their U(1) and SU(2) forms. Furthermore, our approach exemplifies a possible path allowing the exploration of the phase diagram of a variety of fermionic quantum models hitherto considered to be impractical via quantum Monte Carlo simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Two-particle correlation functions in cluster perturbation theory: Hubbard spin susceptibilities

Cluster perturbation theory (CPT) is a computationally economic method commonly used to estimate the momentum- and energy-resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. Finally, we benchmark the method with the one-dimensional Fermi-Hubbard model by comparing with known results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnon heat transport in a two-dimensional Mott insulator

Whether or not anomalies in the thermal conductivity in insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity $\kappa$ and its relationship with the specific heat $c_v$ in the half-filled 2D single-band Hubbard model, using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. At low temperatures where the charge degrees of freedom are gapped-out and $c_v$ exhibits a clear magnon peak, we observe that thermal conductivity $\kappa$ also tends to form a peak at similar temperatures. Reducing temperature further produces a sharp upturn in $\kappa$, associated with an increasing mean-free path. Here we identify this as the high-temperature side of the anomalous peak in insulating cuprates, where the mean-free path eventually is cut-off by other scattering effects, including phonons, disorder, and physical size. Different scattering effects in our model are identified and analyzed in the thermal diffusivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Unconventional Superconductivity at the Crossover between Strong and Weak Electronic Interactions

High-temperature superconductivity emerges in many different quantum materials, often in regions of the phase diagram where the electronic kinetic energy is comparable to the electron-electron repulsion. Describing such intermediate-coupling regimes has proven challenging as standard perturbative approaches are inapplicable. In this study, we employ quantum Monte Carlo methods to solve a multiband Hubbard model that does not suffer from the sign problem and in which only repulsive interband interactions are present. In contrast to previous sign-problem-free studies, we treat magnetic, superconducting, and charge degrees of freedom on an equal footing. We find an antiferromagnetic dome accompanied by a metal-to-insulator crossover line in the intermediate-coupling regime, with a smaller superconducting dome appearing in the metallic region. Across the antiferromagnetic dome, the magnetic fluctuations change from overdamped in the metallic region to propagating in the insulating region. Our findings shed new light on the intertwining between superconductivity, magnetism, and charge correlations in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗