Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level density”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development of cathode-electrolyte-interphase for safer lithium batteries

Accompanied by the adoption of aggressive cathodes to continuously improve batteries energy density, enhancing their safety is becoming increasingly urgent for the electric vehicle development. In-situ controllable formation of robust cathode-electrolyte interphase (CEI) with high inorganic content seems to be the most promising strategy to address the thermal runaway concerns. Moreover, the in-situ formation strategy via suitable electrolyte replacement or electrolyte additives is extremely simple yet effective, especially for industrial manufacture of batteries. Here, this paper briefly reviews recent advanced CEIs formed by conventional carbonate-based electrolytes, fluorinated electrolytes, concentrated electrolytes, and solid state electrolytes. The focus on the thermal stability of the cathodes after CEI modification, and at the same time conduct safety tests at the material, cell, and module levels for comprehensive evaluation are encouraged. The review will provide inspiration for future developments in battery safety and push forward the practical applications of newly developed high-energy density batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Dark Energy Survey Year 1 Results: Cosmological constraints from cluster abundances and weak lensing

In this work we perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Electronic structure, magnetic properties, and pairing tendencies of the copper-based honeycomb lattice Na 2 Cu 2 TeO 6

Spin-1/2 chains with alternating antiferromagnetic (AFM) and ferromagnetic (FM) couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density-matrix renormalization-group (DMRG) methods, we have systematically studied the dimerized chain system Na 2 Cu 2 TeO 6 with a d 9 electronic configuration. Near the Fermi level, in the nonmagnetic phase the dominant states are mainly contributed by the Cu 3d x 2 -y 2 orbitals highly hybridized with the O 2p orbitals, leading to an “effective” single-orbital low-energy model. By calculating the relevant hoping amplitudes, we explain the size and sign of the exchange interactions in Na 2 Cu 2 TeO 6 . In addition, a single-orbital Hubbard model is constructed for this dimerized chain system where the quantum fluctuations are taken into account. Both AFM and FM couplings (leading to an ↑|-|↓-|↓-↑| state) along the chain were found in our DMRG and Lanczos calculations, in agreement with density functional theory and neutron-scattering results. The hole pairing binding energy ΔE is predicted to be negative at Hubbard U~ 11eV, suggesting incipient pairing tendencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Effects of SRPPF ARS Effluent Variations on Magnesium Oxysulfate Solidification

A magnesium oxysulfate (MOS) grout formulation has been identified to solidify the liquid effluent from the Savannah River Plutonium Processing Facility (SRPPF) Aqueous Recovery System (ARS). The formulation uses reactive, light burnt MgO (28 wt%), anhydrous MgSO 4 (10 wt%), and dead burnt MgO (62 wt%), resulting in a grout with good mixability, acceptable density, acceptable setting times, and a leachate pH around 9.4, within the assumed Waste Isolation Pilot Plant (WIPP) brine pH range. A series of tests were then designed to observe how this formula reacted to variations in the liquid effluent. Liquid effluent temperature, the caustic pH level, concentration of NaNO 3 , and the presence of neutralization products, trace metals, and/or sodium sulfate (all potential residuals found in the effluent from upstream processing), were evaluated to see how each impacts the grout mixing time, peak temperature, setting time, presence of bleed water, leachate pH, and density. The standard mix has an average mixing time of 18 minutes, peak temperature of 100.5 °C, and set time of 60 minutes. Increasing the temperature of the liquid effluent increased the reaction rate of the mix and reduced mixing and setting times. Between 30 and 40°C the change was relatively small, raising no more than 6 °C compared to the standard mix maximum temperature. Increasing the pH of the liquid effluent had no significant impact on the mix. The addition of neutralization products and trace metals had an overall impact of decreasing the reactivity of the mix. The addition of more ions in the liquid effluent, such as an increase in the concentration of NaNO 3 and the addition of Na 2 SO 4 , generally decreased the reactivity of the mix. None of the variations to the liquid effluent hindered solidification as indicated by the lack of bleed water found on all the samples. The leachate pH values for all mixes tested did not significantly vary from the expected pH of 9.4 and none of the density values fell below 1.8 g/cm 3 . Overall, changes to the liquid effluent were found to have a minimal impact on the formulation suggesting the grout's ability to properly form despite increased temperatures and the presence of residuals typically found in the liquid effluent from the ARS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Evaluation of High Level Waste Sludge Processing Behavior

The U.S. Department of Energy’s (DOE) Hanford Site has 177 underground storage tanks that contain wastes from past nuclear fuel reprocessing and waste-management operations. Over 20% of this waste is in the form of an insoluble sludge that will require slurry modification before its transfer to the Waste Treatment and Immobilization Plant (WTP). Specific WTP acceptance criteria for waste feed delivery describe the physical and chemical characteristics of the waste that must be met before the waste is transferred to the WTP. One challenging requirement relates to the undissolved solids (UDS) composition in a waste feed because the waste contains solid particles that settle, and their concentration and relative proportion can change during the transfer of the waste in individual batches. A key uncertainty is the ability to transfer and mix wastes with large variations in UDS concentrations and resulting settling rates. To address this uncertainty, a number of small scale mixing and settling tests have been conducted to determine the mobilization performance of variable chemistry simulants. Comparison of the size and density of the particulate for each simulant to that of southeast area Hanford sludge was made using metrics for particle mobilization, suspension, settling, and pipeline transfer where dependance on particle size and density may be different, including: 1. Settling velocity, 2. Critical shear stress for erosion, 3. Just-suspended impeller speed, and 4. Pipeline critical transport velocity. Existing high-level waste sludge data has shown the effect that increasing Al concentration has on resulting settled solids. This differential settling of particles in the sludge has the possibility of resulting in solids segregation during feed preparation and uneven particle distribution during pipeline transportation or mixer jet pump operations. Understanding the predictive capabilities of HLW solids settling and transport as well as potential remedies for addressing disparate sludge behaviors can help provide technical guidance during HLW flowsheet planning.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗
DOE OSTI2023

Functional group effect of chemically modified microcrystalline methyl cellulose on thermoplastic polyurethane composites

Microcrystalline cellulose (MCC) is a promising bio-based filler for lightweight yet mechanically high-performance eco-polymer composites because of its low density and high biocompatibility. However, intermolecular hydrogen bonding among MCC is stronger than the polymer–filler interactions, which deteriorate the mechanical properties of the composites. Herein, we investigated the effects of functional groups on the mechanical properties of composites by scrutinizing chemically modified microcrystalline methyl cellulose (m-MMC) with three different substitution levels of hydroxyl group to hydroxypropyl (HP) group: no-, low-, and high-level substitution (HP-0, HP-low, and HP-high). The degree of HP substitution of m-MMC was quantitatively measured by CP/MAS NMR analysis. The relatively bulky HP groups interrupted the filler–filler intermolecular interactions and reduced the crystallinity and density of m-MMC, as evident from X-ray diffractometer and pycnometer data, respectively. For scalable production of the composites, the m-MMC were compounded with thermoplastic polyurethane (TPU) by a twin-screw extruder at concentrations between 0.5 and 10 wt%. Despite its low filler concentration, the toughness of m-MMC/TPU composites was remarkably enhanced, up to 28% (229.2 to 294.4 MJ/m 3 ) at 0.5 wt% loading of HP-low, owing to the enhanced polymer–filler interactions. In conclusion, the fundamental understanding on structure–property relationships will provide insights for designing of mechanically robust yet eco-friendly polymer composites.

composites↗
DOE OSTI2023

Spin-polarization anisotropy controlled by bending in tungsten diselenide nanoribbons and tunable excitonic states

A WSe 2 monolayer shows many interesting properties due to its spin-orbit coupling induced spin splitting in bands around the Fermi level and the spin-valley configuration. The orientation of the spin polarization in the relevant bands is crucial for the nature of exciton states and the optical valley selectivity. In this work, we studied the WSe 2 nanoribbons under different mechanical bending curvatures and electron/hole doping with density functional theory and their optical absorption and excitonic states with many-body perturbation GW and BSE (Bethe-Salpeter equation) methods. We found that the WSe 2 nanoribbons can exhibit an enhanced SOC effect and a spatially varying spin polarization in bands around the Fermi level under bending. The spin-polarization can show an anisotropy (or asymmetry) in those nearly degenerate bands, leading to a controllable magnetism via bending and electron/hole doping to the nanoribbons, suggesting a potential application in compact and controllable magnetic nanodevices and spintronics. The optical absorption spectrum of the nanoribbon presents a large tunability with bending within the near infrared region of about 0.4 to 1.5 eV, showing an enhanced absorption at a large bending condition. Finally, the exciton states generally show mixed or various spin configuration in the electron and hole pairs that are controlled by bending, potentially useful for applications in spin-based quantum information processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic structure and magnetic tendencies of trilayer La 4 Ni 3 O 10 under pressure: Structural transition, molecular orbitals, and layer differentiation

Motivated by the recent observation of superconductivity in the pressurized trilayer Ruddlesden-Popper (RP) nickelate La 4 Ni 3 O 10 , we explore its structural, electronic, and magnetic properties as a function of hydrostatic pressure from first-principles calculations. We find that an orthorhombic (monoclinic)-to-tetragonal transition under pressure takes place concomitantly with the onset of superconductivity. The electronic structure of La 4 Ni 3 O 10 can be understood using a molecular trimer basis wherein n molecular subbands arise as the d z 2 orbitals hybridize strongly along the c axis within the trilayer. The magnetic tendencies indicate that the ground state at ambient pressure is formed by nonmagnetic inner planes and stripe-ordered outer planes that are antiferromagnetically coupled along the c axis, resulting in an unusual ↑, 0, ↓ stacking that is consistent with the spin density wave model previously suggested by neutron diffraction. Such a state is destabilized at the pressure where superconductivity arises. Despite the presence of d z 2 states at the Fermi level, the d x 2 –y 2 orbitals also play a key role in the electronic structure of La4Ni3O10. Finally, this active role of the d x 2 –y 2 states in the low-energy physics of the trilayer RP nickelate, together with the distinct electronic behavior of the inner and outer planes, resembles the physics of multilayer cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electron Bernstein wave conversion of high-field side injected X-modes in QUEST

This paper presents a detailed design of the Q-shu University experimental steady state spherical tokamak's (QUEST's) high-field side (HFS) injection system for electron Bernstein wave (EBW) excitation and the results of an experimental comparison of the HFS eXtraordinary X-mode and low-field side (LFS) ordinary O-mode injection of 8.2 GHz radio frequency (RF) power. Waveguides, as an alternative to mirror polarizers for transmitting RF X-mode power from LFS to HFS for EBW conversion, were used instead of the installation of an RF mirror. Testing of LFS-to-HFS RF power transmission at 8.2 GHz, using an RG-50-type vacuum waveguide in a bench-scale device filled with SF6 gas at 0.03 Mpa, revealed that an RF power of 10.8 kW could traverse the fundamental electron cyclotron resonance layer for 60 s without breakdown. The short-length, open-ended waveguide antenna used in the HFS injection-induced wave diffraction reduced the efficiency of power delivery to the upper hybrid resonance (UHR) by approximately 7% at an electron temperature of 50 eV. The HFS injection was able to produce brighter camera images than the standard LFS injection. The location of the UHR, as estimated by measuring the density with an interferometer, agreed with its location as measured by plasma radiation low-field, side-edge positions shown by fast camera imaging. This indicates that the plasma was produced by mode-converted EBW. The HFS injection had an absorption efficiency of 96%, compared to 40% for LFS. A greater fluctuation of floating potential adjustable to the lower hybrid wave (LHW) was observed in the HFS case by installing a Langmuir probe, confirming that EBW conversion efficiency was higher in the HFS case. Moreover, after setting the poloidal field to B PF = 7.6 mT, plasma current (I P ) in the HFS peaked at 1.3 kA, as opposed to 0.3 kA for LFS, despite LFS injection having a total power of 55 kW, compared to 40 kW for HFS. However, as the impurity level was comparatively high, it is believed that this I P is dominated by pressure-drive, which makes it difficult to analyze EBWCD. Finally, the line-integrated density in the HFS injection peaked at 1.6 × 10 18 m -2 , compared to 8 × 10 17 m -2 in the LFS one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Coherent Phonon-Induced Modulation of Charge Transfer in 2D Hybrid Perovskites

Electron–phonon interactions play an essential role in charge transport and transfer processes in semiconductors. For most structures, tailoring electron–phonon interactions for specific functionality remains elusive. Here, it is shown that, in hybrid perovskites, coherent phonon modes can be used to manipulate charge transfer. In the 2D double perovskite, (AE2T) 2 AgBiI 8 (AE2T: 5,5“-diylbis(amino-ethyl)-(2,2”-(2)thiophene)), the valence band maximum derived from the [Ag 0.5 Bi 0.5 I4] 2– framework lies in close proximity to the AE2T-derived HOMO level, thereby forming a type-II heterostructure. During transient absorption spectroscopy, pulsed excitation creates sustained coherent phonon modes, which periodically modulate the associated electronic levels. Thus, the energy offset at the organic–inorganic interface also oscillates periodically, providing a unique opportunity for modulation of interfacial charge transfer. Density-functional theory corroborates the mechanism and identifies specific phonon modes as likely drivers of the coherent charge transfer. These observations are a striking example of how electron–phonon interactions can be used to manipulate fundamentally important charge and energy transfer processes in hybrid perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Topotactic Reduction‐Driven Crystal Field Excitations in Brownmillerite Manganite Thin Films

Topotactic reduction of perovskite oxides offers a powerful approach for discovering novel phenomena, such as superconducting infinite-layer nickelates and polar metallicity, and is commonly accompanied by the emergence of multiple valence states and/or complex crystal fields of transition metals. However, understanding the complex interplay between crystal chemistry, electronic structure, and physical properties at the spin- and orbital-resolved levels in these reduced systems remains elusive. Here, in this study, x-ray absorption spectroscopy, resonant inelastic x-ray scattering (RIXS), and density functional theory calculations are used to uncover topotactic metal-insulator transition and orbital-specific crystal field excitations in brownmillerite La 0.67 Ca 0.33 MnO 2.5 thin films. The Mn valence states are found to be Mn 2+ /Mn 3+ , along with their corresponding populations at octahedral and tetrahedral sites, which effectively weaken the Mn-O hybridization compared to the parent perovskite phase. As a result, La 0.67 Ca 0.33 MnO 2.5 films exhibit an antiferromagnetic insulating ground state. Moreover, by combining the RIXS measurements on selected single-valence manganites, specifically MnO, LaMnO 3 , and CaMnO 3 , with orbital- and spin-resolved density-of-states calculations, the study identifies the dd excitations of octahedrally and tetrahedrally coordinated Mn 2+ /Mn 3+ ions, directly linking the microscopic electronic structure to the macroscopic magnetic/electrical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enabling Deformable and Stretchable Batteries

Abstract Deformable energy storage devices are needed to power next‐generation wearable electronics that interface intimately with human skin. Currently, deformable energy storage devices demonstrate poor performance compared to their rigid lithium‐ion counterparts, forcing wearable manufacturers to design their devices around bulky battery compartments. However, technological advances to create deformable batteries at the component and device level have yielded continuous improvement in stretchable batteries over the last five years. In this Essay, the major strategies at the component and device level that have been successfully employed to create stretchable batteries are reviewed. The outstanding challenges facing deformable energy storage are also discussed, namely, energy density, packaging, delamination, device integration, and manufacturing. This Essay will give researchers who are interested in contributing to the development of deformable batteries a cursory understanding of the most successful strategies to date, and provide insights into the most important directions to pursue in the future.

Mackanic, David G.↗
DOE OSTI2021

Temperature and radiation effects on brittle versus ductile fracture behavior in miscible phase boundaries: insight from atomistic simulations

Temperature- and irradiation-assisted failure mechanisms in miscible phase boundaries are clarified via atomistic calculations. We first establish the temperature-dependent brittle-to-ductile transition in U–Zr miscible phase boundaries. Our results confirm that these boundaries are mostly brittle at low temperatures and ductile at elevated temperatures. We then investigate the changes induced by irradiation on the fracture mechanisms in such phase boundaries. The irradiation-induced defect accumulation follows a multi-stage process that starts with the accumulation of isolated small dislocation loops before transitioning to the saturation and growth of larger dislocation loops and end up with a reorganization into forest dislocations. The accumulation of loops is the primary feature to participate in the delineation between brittle and ductile interfacial fracture in irradiated phase boundaries. At low damage levels, radiation defect interactions with the crack tip are limited and U–Zr miscible boundaries fail through the emission of dislocations ahead of the crack tip followed by brittle cleavage in agreement with the classical Griffith’s criterion for crack stability. At higher damage levels, the failure mode transitions from brittle crack growth to ductile void growth. In this case, the microcrack tip is blunted by the high density of pre-existing, radiation-induced defects in the vicinity of the crack. This interaction leads to the development and growth of a cavity at the interface as opposed to interfacial cleavage.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A spin-wave magnetometer with a positive feedback

In this work, we demonstrate experimentally the operation of a spin-wave magnetometer integrated into a circuit with a positive feedback. The circuit consists of the passive magnetic and active electric parts. The magnetic part includes a sensing element, which is a magnetic cross junction made of Y 3 Fe 2 (FeO 4 ) 3 . The electric part includes a non-linear amplifier and a phase shifter. The electric and magnetic parts are connected via micrometer size antennae. Spin waves are excited by two of these antennae while the output inductive voltage produced by the interfering spin waves is detected by the third antenna. Spin waves propagating in the orthogonal arms of the cross can accumulate significantly different phase shifts, depending on the direction and the strength of the external magnetic field. The output inductive voltage reaches its maximum in the case of constructive spin wave interference. The positive feedback provides further signal amplification. It appears possible to enhance the response function, compared to the passive circuits without a feedback, by a factor of ×100 without an increase in the noise level. The experimental data show a prominent response to the external magnetic field variation, exceeding 5 x 10 3 V/T. The intrinsic noise spectral density of the device can be as low as 10 -16 V 2 /Hz. The estimated sensitivity of the prototype device is 2 x 10 -12 T/√Hz at room temperature. We argue that spin-wave magnetometers can potentially be as sensitive as SQUIDs while operating at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2021

Excitonic density wave and spin-valley superfluid in bilayer transition metal dichalcogenide

Artificial moiré superlattices in 2 d van der Waals heterostructures are a new venue for realizing and controlling correlated electronic phenomena. Recently, twisted bilayer WSe 2 emerged as a new robust moiré system hosting a correlated insulator at moiré half-filling over a range of twist angle. In this work, we present a theory of this insulating state as an excitonic density wave due to intervalley electron–hole pairing. We show that exciton condensation is strongly enhanced by a van Hove singularity near the Fermi level. Our theory explains the remarkable sensitivity of the insulating gap to the vertical electric field. In contrast, the gap is weakly reduced by a perpendicular magnetic field, with quadratic dependence at low field. The different responses to electric and magnetic field can be understood in terms of pair-breaking versus non-pair-breaking effects in a BCS analog of the system. We further predict superfluid spin transport in this electrical insulator, which can be detected by optical spin injection and spatial-temporal imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quantum-classical embedding via ghost Gutzwiller approximation for enhanced simulations of correlated electron systems

Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. Here, this work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.

Chen, I-Chi [Ames Laboratory (AMES), Ames, IA (Uni↗