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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

The structure of the planets Jupiter and Saturn

Planetary models for Jupiter and Saturn are computed using a fourth-order theory and an equation of state for the molecular hydrogen and helium planetary envelopes taken from the Monte Carlo calculations of Slattery and Hubbard (1976). Models for Jupiter are found that have a small amount of heavy elements either mixed with hydrogen and helium throughout the interior of the planet or concentrated in a small dense core. Saturn is modeled with a solar-composition hydrogen and helium envelope and a small dense core. It is concluded that the molecular equation of state linked with suitable interior equations of state can produce Jovian models which satisfy the observational data. The planetary models show that the enrichment of heavy elements (relative to solar composition) is approximately 3 times for Jupiter and 10 times for Saturn.

Slattery, W. L.↗
DOE OSTI2024

Magnetic phase diagram of a two-orbital model for bilayer nickelates with varying doping

Motivated by the recently discovered high-T c bilayer nickelate superconductor La 3⁢ Ni 2 ⁢O 7 , we comprehensively research a bilayer 2×2×2 cluster for different electronic densities n by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying n. Based on the spin structure factor S(q), we have obtained a rich magnetic phase diagram in the plane defined by n and U/W, at fixed Hund coupling, where U is the Hubbard strength and W the bandwidth. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. At half-filling, n=2 (two electrons per Ni site, corresponding to N=16 electrons), the canonical superexchange interaction leads to a robust G-AFM state (π,π,π) with antiferromagnetic couplings both in-plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the “half-empty” and “half-full” mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At n=1.5 (or N=12), density corresponding to La 3 ⁢Ni 2 ⁢O 7 , we obtained the “Stripe 2” ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the z axis) in the 2×2×2 cluster. In addition, we obtained a much stronger AFM coupling along the z axis than the magnetic coupling in the xy plane. The random-phase approximation calculations with varying n give very similar results as Lanczos, even though both techniques are based on quite different procedures. Additionally, a state with q/π=(0.6,0.6,1) close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to n=1.25, possibly responsible for the E-phase SDW recently observed in experiments. In conclusion, our predictions can be tested by chemically doping La 3 ⁢Ni 2 ⁢O 7 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Bootstrap embedding for interacting electrons in phonon coherent-state mean field

Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Islam, Shariful [North Carolina State University, ↗
NASA NTRSDate not supplied

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗
DOE OSTI2022

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗
DOE OSTI2022

Quasiparticle excitations in a one-dimensional interacting topological insulator: Application for dopant-based quantum simulation

Here, we study the effects of electron-electron interactions on the charge excitation spectrum of the spinful Su-Schrieffer-Heeger (SSH) model, a prototype of a one-dimensional bulk obstructed topological insulator. In view of recent progress in the fabrication of dopant-based quantum simulators we focus on experimentally detectable signatures of interacting topology in finite lattices. To this end we use Lanczos-based exact diagonalization to calculate the single-particle spectral function in real space which generalizes the local density of states to interacting systems. Its spatial and spectral resolution allows for the direct investigation and identification of edge states. By studying the noninteracting limit, we demonstrate that the topological in-gap states on the boundary are robust against both finite-size effects as well as random bond and onsite disorder which suggests the feasibility of simulating the SSH model in engineered dopant arrays in silicon. While edge excitations become zero-energy spinlike for any finite interaction strength, our analysis of the spectral function shows that the single-particle charge excitations are gapped out on the boundary. Despite the loss of topological protection, we find that these edge excitations are quasiparticlelike as long as they remain within the bulk gap. Above a critical interaction strength of U c ≈5t these quasiparticles on the boundary lose their coherence which is explained by the merging of edge and bulk states. This is in contrast to the many-body edge excitations which survive the limit of strong coupling, as established in the literature. Our findings show that for moderate repulsive interactions the nontrivial phase of the interacting SSH model can be detected through remnant signatures of topological single-particle states using single-particle local measurement techniques such as scanning tunneling spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Oxygen magnetic polarization, nodes in spin density, and zigzag spin order in oxides

Recent studies on Ba 2 CoO 4 (BCO) and SrRuO 3 (SRO) have unveiled a variety of intriguing phenomena, such as magnetic polarization on oxygens, unexpected nodes in the spin density profile along bonds, and unusual zigzag spin patterns in triangular lattices. Here, using simple model calculations supplemented by density functional theory we explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the transition metal (TM) to oxygen bonds. Our simple model also allows us to explain the net polarization on oxygen as originated from the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Overall, our results are not limited to BCO and SRO, but they are generic and qualitatively predict the net polarization expected on any ligands, according to the spin order of the closest TM atoms and the number of intermediate ligand atoms. Finally, we propose that a robust easy-axis anisotropy would suppress the competing 120° antiferromagnetic order to stabilize the zigzag pattern order as ground state in a triangular lattice. Our generic predictions should be applicable to any other compound with characteristics similar to those of BCO and SRO.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La 4 ⁢Ni 3 ⁢O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4⁢ Ni 3 ⁢O 10 , analogous to the dimers of La 3 ⁢Ni 2 ⁢O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3⁢z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 ⁢Ni 2 ⁢O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3⁢ Ni 2 ⁢O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 ⁢Ni 3 ⁢O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Complexity Phase Diagram for Interacting and Long-Range Bosonic Hamiltonians

Here, we classify phases of a bosonic lattice model based on the computational complexity of classically simulating the system. We show that the system transitions from being classically simulable to classically hard to simulate as it evolves in time, extending previous results to include on-site number-conserving interactions and long-range hopping. Specifically, we construct a complexity phase diagram with easy and hard “phases” and derive analytic bounds on the location of the phase boundary with respect to the evolution time and the degree of locality. We find that the location of the phase transition is intimately related to upper bounds on the spread of quantum correlations and protocols to transfer quantum information. Remarkably, although the location of the transition point is unchanged by on-site interactions, the nature of the transition point does change. Specifically, we find that there are two kinds of transitions, sharp and coarse, broadly corresponding to interacting and noninteracting bosons, respectively. Our Letter motivates future studies of complexity in many-body systems and its interplay with the associated physical phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Adsorption-Induced Surface Magnetism

We report the emergence of adsorption-induced magnetism from heterohelicene molecules on a nonmagnetic Cu(100) surface. Spin-polarized low-energy electron microscopy measurements reveal spin-dependent electron reflectivity for enantiopure 7,12,17-trioxa[11]helicene (TO[11]H) monolayers, indicating the formation of a spin-polarized state localized in the topmost copper layer. Control experiments on clean Cu(100) and TO[11]H on highly oriented pyrolytic graphite show no such effect, excluding artifacts and chirality-induced spin selectivity as origins. Spin-polarized density functional theory calculations with hybrid functionals attribute the magnetism to strong chemisorption, which induces hybridization between the molecular HOMO and copper s- and d-states, driving asymmetric spin-polarized charge redistribution at the interface. An extended Newns-Anderson-Grimley model incorporating on-site Coulomb repulsion in Cu d-orbitals reproduces the emergence of interfacial spin polarization above a threshold interaction strength, highlighting the key roles of hybridization parameters and Coulomb correlation. These findings reveal a mechanism for inducing magnetism at molecule-metal interfaces without inherently magnetic components, offering avenues for engineering spin-polarized states in organic-inorganic hybrid systems.

Hubbard model↗
DOE OSTI2025

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Design Rules for Open-Shell Molecular Wires: Insights from Correlated Many-Body Transport

Open-shell π-conjugated systems have emerged as promising molecular wires that can sustain unusually high low-bias conductance over tens of nanometers. However, predictive design rules remain limited because near-degeneracy and spin polarization give rise to correlated electronic states that are not reliably captured by standard density functional theory and Landauer transport descriptions. Here, we establish physically transparent design rules for long-range, weakly length-dependent charge transport based on the interplay between bond-length alternation and electron–electron interactions. Using a fully correlated transport framework that combines density matrix renormalization group calculations with nonequilibrium Green’s function embedding, we show that weak dimerization maintains spatially extended edge states, while moderate on-site interaction stabilizes open-shell character without excessively separating transport-relevant resonances. This balance yields zero-bias transmission that is largely insensitive to molecular length, providing a route to connect chemically tunable structure to correlated transport beyond mean-field descriptions.

Charge transport↗
DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

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

Ground-state phase diagram of the t-t ′ -J model

Significance Theoretical efforts to understand the superconducting cuprates have focused on a few simplified but challenging models: notably the Hubbard and the t - t ′ - J models. Using density matrix renormalization group (DMRG) simulations on width-8 cylinders, we study the ground-state phase diagram of the t - t ′- J model. The next nearest-neighbor hopping t ′ allows one to distinguish between hole ( t ′ < 0 ) and electron doping ( t ′ > 0 ). We find for the magnetic, charge, and single-particle properties that this model roughly captures the differences between the electron- and hole-doped cuprates. However, the superconducting properties differ. For hole doping ( t ′ < 0 ), there is an absence of superconductivity, while for electron doping ( t ′ > 0), we find d -wave superconductivity coexisting with antiferromagnetism and a secondary triplet p -wave pairing component.

Jiang, Shengtao↗