Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extended 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 55 records · Page 3

Extending the Gutzwiller approximation to intersite interactions

In this work, we develop an extension of the Gutzwiller approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping, and Coulomb-assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multiband Hubbard models in the presence of only local interactions. This is accomplished by combining the 1 / z expansion—where z is the coordination number, and only the leading-order terms contribute in the limit of infinite dimensions—with a P R † P R - I expansion, where P R is the Gutzwiller projector on the site R . Furthermore, the method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single-band Hubbard model. Similarly to the usual Brinkman-Rice mechanism, we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wave function is enhanced with respect to the uncorrelated Slater determinant wave function.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A random-sampling method as an efficient alternative to variational Monte Carlo for solving Gutzwiller wavefunctions

Abstract We present a random-sampling (RS) method for evaluating expectation values of physical quantities using the variational approach. We demonstrate that the RS method is computationally more efficient than the variational Monte Carlo method using the Gutzwiller wavefunctions applied on single-band Hubbard models as an example. Non-local constraints can also been easily implemented in the current scheme that capture the essential physics in the limit of strong on-site repulsion. In addition, we extend the RS method to study the antiferromagnetic states with multiple variational parameters for 1D and 2D Hubbard models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Statistical mechanics of light elements at high pressure. IV - A model free energy for the metallic phase

A large quantity of data on the thermodynamic properties of hydrogen-helium metallic liquids have been obtained in extended computer calculations in which a Monte Carlo code essentially identical to that described by Hubbard (1972) was used. A model free energy for metallic hydrogen with a relatively small mass fraction of helium is discussed, taking into account the definition of variables, a procedure for choosing the free energy, values for the fitting parameters, and the evaluation of the entropy constants. Possibilities concerning a use of the obtained data in studies of the interiors of the outer planets are briefly considered.

Dewitt, H. E.↗

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↗

Transformer-based operator learning framework for self-energy in strongly correlated systems

We introduce Σ-Attention, a transformer-based operator-learning framework for approximating the self-energy operator of strongly correlated electronic systems. By creating a batched dataset that combines results from three complementary approaches, i.e., many-body perturbation theory, strong-coupling expansion, and exact diagonalization, each effective in specific parameter regimes, Σ-Attention is applied to learn an accurate approximation for the self-energy operator that is valid across a wide range of parameter regimes. This hybrid strategy leverages the strengths of existing methods while relying on the transformer's ability to generalize beyond individual limitations. More importantly, the scalability of the transformer architecture allows the learned self-energy to be extended to systems with larger sizes, leading to much improved computational scaling. Using the one-dimensional Hubbard model, we demonstrate that Σ-Attention can accurately predict the Matsubara Green's function of large systems with a wide range of coupling strength. Our framework offers a promising and scalable pathway for studying strongly correlated systems with many possible generalizations.

Zhu, Yuanran↗

Localization in Energy Materials (Final Project Report)

The last year of this award we have continued our research on using quantum machine learning to identify phase transitions. By combining quantum machine learning with quantum computing, we extended a hybrid classical-quantum algorithm to capture the metal-insulator quantum phase transition of the Hubbard model. We have also continued our studies of non-equilibrium dynamics of interacting disordered systems. In particular, we studied the non-equilibrium transient dynamics of a system described by the Anderson-Hubbard model following an interaction and disorder quench, and the effect of disorder on the out-of-time-order correlator on the Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Auxiliary field Quantum Monte Carlo for dilute neutrons on the lattice

We employ constrained path Auxiliary Field Quantum Monte Carlo (AFQMC) in the pursuit of studying physical nuclear systems using a lattice formalism. Since AFQMC has been widely used in the study of condensed-matter systems such as the Hubbard model, we benchmark our method against published results for both one- and two-dimensional Hubbard model calculations. We then turn our attention to cold atomic and nuclear systems. We use an onsite contact interaction that can be tuned in order to reproduce the known scattering length and effective range of a given interaction. Developing this machinery allows us to extend our calculations to study nuclear systems within a lattice formalism. In conclusion, we perform initial calculations for a range of nuclear systems from two- to few-body neutron systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Downfolding Complex Materials Problems Onto Model Hamiltonians for Quantum Computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios Markos Alvertis↗

Downfolding complex materials problems onto model Hamiltonians for quantum computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios M. Alvertis↗

A path integral ground state Monte Carlo algorithm for entanglement of lattice bosons

A ground state path integral quantum Monte Carlo algorithm is introduced that allows for the study of entanglement in lattice bosons at zero temperature. The Rényi entanglement entropy between spatial subregions is explored across the phase diagram of the one dimensional Bose-Hubbard model for systems consisting of up to L=256 L = 256 sites at unit-filling without any restrictions on site occupancy, far beyond the reach of exact diagonalization. The favorable scaling of the algorithm is demonstrated through a further measurement of the Rényi entanglement entropy at the two dimensional superfluid-insulator critical point for large system sizes, confirming the existence of the expected entanglement boundary law in the ground state. The Rényi estimator is extended to measure the symmetry resolved entanglement that is operationally accessible as a resource for experimentally relevant lattice gases with fixed total particle number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Emergent Superconductivity and Competing Charge Orders in Hole-Doped Square-Lattice t – J Model

The square-lattice Hubbard and closely related t-J models are considered as basic paradigms for understanding strong correlation effects and unconventional superconductivity (SC). Recent large-scale density matrix renormalization group simulations on the extended t-J model have identified d-wave SC on the electron-doped side (with the next-nearest-neighbor hopping t 2 > 0) but a dominant charge density wave (CDW) order on the hole-doped side (t 2 < 0), which is inconsistent with the SC of hole-doped cuprate compounds. We re-examine the ground-state phase diagram of the extended t-J model by employing the state-of-the-art density matrix renormalization group calculations with much enhanced bond dimensions, allowing more accurate determination of the ground state. On six-leg cylinders, while different CDW phases are identified on the hole-doped side for the doping range δ = 1/16 – 1/8, a SC phase emerges at a lower doping regime, with algebraically decaying pairing correlations and d-wave symmetry. On the wider eight-leg systems, the d-wave SC also emerges on the hole-doped side at the optimal 1=8 doping, demonstrating the winning of SC over CDW by increasing the system width. Furthermore, our results not only suggest a new path to SC in general t-J model through weakening the competing charge orders, but also provide a unified understanding on the SC of both hole- and electron-doped cuprate superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Cosmic rays and the emission line regions of active galactic nuclei

The effects that the synchrotron emitting relativistic electrons could have on the emission line regions which characterize active nuclei are discussed. Detailed models of both the inner, dense, broad line region and the outer, lower density, narrow line region are presented, together with the first models of the optically emitting gas often found within extended radio lobes. If the relativistic gas which produces the synchrotron radio emission is mixed with the emission line region gas then significant changes in the emission line spectrum will result. The effects of the synchrotron emitting electrons on filaments in the Crab Nebula are discussed in an appendix, along with a comparison between the experimental calculations, which employ the mean escape probability formalism, and recent Hubbard and Puetter models.

Ferland, G. J.↗

Cosmic rays and the emission-line regions of active galactic nuclei

The effects that the synchrotron emitting relativistic electrons could have on the emission line regions which characterize active nuclei are discussed. Detailed models of both the inner, dense, broad line region and the outer, lower density, narrow line region are presented, together with the first models of the optically emitting gas often found within extended radio lobes. If the relativistic gas which produces the synchrotron radio emission is mixed with the emission line region gas then significant changes in the emission line spectrum will result. The effects of the synchrotron emitting electrons on filaments in the Crab Nebula are discussed in an appendix, along with a comparison between the experimental calculations, which employ the mean escape probability formalism, and recent Hubbard and Puetter models.

Ferland, G. J.↗

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exact solution of two simple non-equilibrium electron-phonon and electron-electron coupled systems: The atomic limit of the Holstein-Hubbard model and the generalized Hatsugai-Komoto model

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales—for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections, but strongly correlated electrons require a more complex description, and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. Then, we examine the photoemission spectra of this model in and out of equilibrium. Next, we use similar methods to exactly solve an extended version of the Hatsugai-Komoto model and examine its behavior in and out of equilibrium. These calculations lead us to propose using the first moment of the photoemission spectra to signal non-equilibrium changes in electron-electron and electron-phonon couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗