Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “many-body methods”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗
DOE OSTI2025

Impact of hole polaron formation on excitonic transitions in MgO from first principles

Here, we present a first-principles investigation of the excitonic properties of magnesia (MgO), an ionic insulator known to host hole polarons. We combine a density functional theory-based approach for structural relaxation in the presence of the hole and many-body perturbation theory to describe the excitonic properties. We determine that the hole polaron introduces new in-gap occupied states 0.6–0.8 eV above the valence band maximum that lead to two low-energy peaks in the optical spectrum. The predicted redshift of the lowest-energy transition due to polaron formation of 0.8 eV agrees well with the experimental Stokes shift of 0.8–0.9 eV. Analysis of the exciton wave function indicates that the electron-hole pair consists of a localized hole and delocalized electron, but that the wave function retains its Wannier-Mott character even in the presence of the hole polaron. Our study demonstrates that combining these previously established methods allows for a relatively computationally inexpensive approach to studying the exciton polaron in materials where only one charge carrier forms a polaron.

electronic structure↗
DOE OSTI2026

Correlation function metrology for warm dense matter: Recent developments and practical guidelines

X-ray Thomson scattering (XRTS) has emerged as a valuable diagnostic for matter under extreme conditions, as it captures the intricate many-body physics of the probed sample. Recent advances, such as the model-free temperature diagnostic of Dornheim et al. [Nat. Commun. 13 , 7911 (2022)], have demonstrated how much information can be extracted directly within the imaginary-time formalism. However, since the imaginary-time formalism is a concept often difficult to grasp, we provide here a systematic overview of its theoretical foundations and explicitly demonstrate its practical applications to temperature inference, including relevant subtleties. Furthermore, we present recent developments that enable the determination of the absolute normalization, Rayleigh weight, and density from XRTS measurements without reliance on uncontrolled model assumptions. Finally, we outline a unified workflow that guides the extraction of these key observables, offering a practical framework for applying the method to interpret experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Tunable Few-Layer van der Waals Crystals and Heterostructures as Emerging Energy and Quantum Materials (Final Technical Report)

2D and layered (van der Waals) semiconductors offer extraordinary opportunities for manipulating optically excited charge carriers, many-body excitations, and non-charge based quantum numbers. To date, research has focused on a limited group of materials, mostly transition metal dichalcogenides in the monolayer limit. Other van der Waals semiconductors, and especially few-layer to multilayer crystals and their heterostructures, carry large potential for the discovery of phenomena of interest for future energy and information technologies. But they remain largely unexplored, often due to a lack of access to high-quality materials and approaches for measuring their properties at the relevant scales. The goal of this project was to develop an EPSCoR-State/National Laboratory Partnership that addresses the challenges of preparing high-quality van der Waals semiconductors and of probing their structure, composition, and especially their optoelectronic and photonic properties, near the atomic scale using electron microscopy techniques. A central component of the project was the development of advanced methods for electron microscopy and electron-excited spectroscopy, taking advantage of unique samples as well as leading capabilities and expertise at the partner institutions. Efforts to advance leading-edge techniques was supported by ancillary developments, such as precision sample preparation for electron microscopy/spectroscopy, coordinated chemical imaging, and analytical electron microscopy. Experiments in materials synthesis and technique development were closely linked to theory and computation. The results obtained under this project yielded multifaceted benefits to the involved partners and their institutions, DOE-BES, the wider scientific community, and society at large, particularly in the State of Nebraska through dividends from knowledge and human capital generated under the project.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is #BQP-complete in the worst case. In this work, we introduce QFAMES (quantum filtering and analysis of multiplicities in eigenvalue spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

d -Electron Heavy Fermion Behavior in a Near-Room-Temperature Polar Metallic Ferrimagnet: A Case of Mn 5 SiC

Polycrystalline Mn 5 SiC was synthesized by using a high-temperature solid-state method. Mn 5 SiC adopts a polar space group (Cmc2 1 ) with six crystallographic Mn sites confirmed by X-ray and neutron diffraction, transmission electron microscopy, and second harmonic generation experiments. The complex crystal structure features edge-sharing trigonal prisms and icosahedra, as well as face/edge-sharing pentagonal prisms. Magnetic measurements indicate ferrimagnetic ordering with a transition temperature of 284 K. The ferrimagnetic structure (magnetic space group Cm’c’2 1 ) was further identified by powder neutron diffraction, where collinear Mn spins align along the crystallographic c-axis. The refined magnetic moment for each crystallographic Mn site at 4 K is 1.8(2), −2.42(9), −1.72(8), 0.51(6), 0.50(4), and 1.7(2) μ B . Density functional theory calculations confirm both the metallic behavior and the ferrimagnetic structure observed experimentally and further provide insight into the observed Mn moment dependence across crystallographic sites. The resistivity and specific heat measurements and density functional theory calculations reveal a substantially large Kadowaki–Woods ratio of 5 × 10 –5 μΩ·cm/(mJ/mol) 2 and a many-body renormalization factor of 5.5, indicating the unusual heavy Fermion behavior in such an itinerant magnetic metal.

Messegee, Zachary T. [George Mason University, Fai↗
DOE OSTI2026

Toward the “platinum standard” of quantum chemistry on quantum computers: Perturbative quadruple corrections in unitary coupled cluster theory

We propose a non-iterative, post-hoc correction to the unitary coupled cluster theory with the single, double, and triple excitations (UCCSDT) Ansatz, which considers the leading-order effects of neglected quadruple excitations. We present two ways to derive this correction, henceforth referred to as [Q-6], which leads to an improvement in the correlation energy shown to be truncated to sixth-order in many-body perturbation theory. Furthermore, a comparison between the UCC-based [Q-6] correction proposed in this work and analogous, “platinum standard” quadruple corrections proposed in conventional coupled cluster theory recognizes that [Q-6] is distinct from prior corrections since it is constructed entirely from internally connected components. Although trotterized (t) and full operator variants of UCCSDT exhibit errors in scans of small molecule potential energy surfaces that routinely exceed 1.6 mH, we find that t/UCCSDT[Q-6] is, nevertheless, able to achieve chemical accuracy as measured by the mean unsigned error.

Correlation energy↗
DOE OSTI2025

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

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

Coupled Lindblad Pseudomode Theory for Simulating Open Quantum Systems

Coupled Lindblad pseudomode theory is a promising approach for simulating non-Markovian quantum dynamics on both classical and quantum platforms, with dynamics that can be realized as a quantum channel. We provide theoretical evidence that the number of coupled pseudomodes only needs to scale as polylog⁡(𝑇/𝜖) in the simulation time 𝑇 and precision 𝜖. Inspired by the realization problem in control theory, we also develop a robust numerical algorithm for constructing the coupled modes that avoid the nonconvex optimization required by existing approaches. We demonstrate the effectiveness of our method by computing population dynamics and absorption spectra for the spin-boson model. Furthermore, this Letter provides a significant theoretical and computational improvement to the coupled Lindblad framework, which impacts a broad range of applications from classical simulations of quantum impurity problems to quantum simulations on near-term quantum platforms.

Anderson impurity model↗
DOE OSTI2026

Quantum state preparation with resolution refinement

We introduce a method called resolution refinement that allows one to bootstrap eigenstate preparation on a quantum computer. We first prepare an eigenstate of a low-resolution Hamiltonian using any method of choice. The eigenstate is then lifted to higher resolution and adiabatically evolved to produce the corresponding eigenstate of a higher-fidelity Hamiltonian. We give examples of resolution refinement applied to both single-particle basis states as well as a spatial lattice grid. For basis refinement, we compute few-body ground states of the Busch model for interacting particles in a harmonic trap in one dimension. For lattice refinement, we compute Hartree-Fock nuclear states for a central Woods-Saxon potential in three dimensions, and we compute bound states and continuum states in a multi-species Hubbard model of fermions in one dimension. In all cases, the method is efficient and requires an adiabatic evolution time that scales with the inverse of the energy gap times the square root of the system size. We show that this very favorable scaling arises from the fact that resolution refinement does not make large changes to the structure or energies of the low-energy eigenstates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Spontaneous magnon decays from nonrelativistic time-reversal symmetry breaking in altermagnets

Quasiparticles are central to condensed matter physics, but their stability can be undermined by quantum many-body interactions. Magnons, i.e., quasiparticles in quantum magnets, are particularly intriguing because their properties are governed by both real and spin space. While crystal symmetries may be low, spin interactions often remain approximately isotropic, limiting spontaneous magnon decay. Textbook wisdom holds that collinear Heisenberg magnets follow a dichotomy: ferromagnets host stable magnons, while antiferromagnetic magnons may decay depending on dispersion curvature. Up to now, relativistic spin-orbit coupling and noncollinear order that connect spin space to real space were shown to introduce more complex magnon instability mechanisms. Here, we show that even in nonrelativistic isotropic collinear systems, this conventional dichotomy is disrupted in altermagnets. Altermagnets, a newly identified class of collinear magnets, exhibit compensated spin order with nonrelativistic time-reversal symmetry breaking and even-parity band splitting. Using kinematic analysis, nonlinear spin-wave theory, and quantum simulations, we reveal that even weak band splitting opens a decay phase space, driving quasiparticle breakdown. Additionally, 𝑑-wave altermagnets form a rare “island of stability” at the Brillouin-zone center. Furthermore, our findings establish a quasiparticle stability trichotomy in collinear Heisenberg magnets and position altermagnets as a promising platform for unconventional spin dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Symmetry dilemmas in quantum computing for chemistry: A comprehensive analysis

Symmetry adaptation, universality, and gate efficiency are central but often competing requirements in quantum algorithms for electronic structure and many-body physics. For example, fully symmetry-adapted universal operator pools typically generate long and deep quantum circuits; gate-efficient universal operator pools generally break symmetries; and gate-efficient, fully symmetry-adapted operator pools may not be universal. In this work, we analyze such symmetry dilemmas both theoretically and numerically. On the theory side, we prove that the popular, gate-efficient operator pool consisting of singlet spin-adapted singles and perfect-pairing doubles is not universal when spatial symmetry is enforced. To demonstrate the strengths and weaknesses of the three types of pools, we perform numerical simulations using an adaptive algorithm paired with operator pools that are (i) fully symmetry-adapted and universal, (ii) fully symmetry-adapted and non-universal, and (iii) breaking a single symmetry and universal. Our numerical simulations encompass three physically relevant scenarios in which the target state is (i) the global ground state, (ii) the ground state crossed by a state differing in multiple symmetry properties, and (iii) the ground state crossed by a state differing in a single symmetry property. Our results show when symmetry-breaking but universal pools can be used safely, when enforcing at least one distinguishing symmetry suffices, and when a particular symmetry must be rigorously preserved to avoid variational collapse. Together, the formal and numerical analyses provide a practical guide for designing and benchmarking symmetry-adapted operator pools that balance universality, resource requirements, and robust state targeting in quantum simulations for chemistry.

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

Variations in the charge radii of indium isotopes between 𝑁=52 and 82

Measurements of the 5⁢𝑝 2 𝑃 3/2 → 9⁢𝑠 2 𝑆 1/2 and 5⁢𝑝 2 𝑃 1/2 → 8⁢𝑠 2 𝑆 1/2 transitions in the indium atom, combined with new atomic physics calculations, were used to extract the changes in mean-square nuclear charge radii, 𝛿⁡⟨𝑟 2 ⟩, of the indium (𝑍 = 49) isotopes 101–111 In , 113–123 In , and 125–131 In . With a proton hole in the closed nuclear shell of 𝑍 = 50, indium provides a detailed study of the effect of unpaired nucleons adjacent to the proton-shell closure, allowing investigation into the charge radii for isotopes between the two major neutron-shell closures at 𝑁 = 50 and 𝑁 = 82. A study of the variations in charge radii between neighboring isotopes with neutron number (the ‘odd-even staggering') is presented and provides further insight and challenges for the theoretical description of the size of proton-hole nuclei. Two nuclear theories, density functional theory and the valence-space in-medium similarity renormalization group method, were employed to interpret the data. The new information obtained in this work provides valuable insights into the successes and shortcomings of the theoretical approaches employed.

ab initio calculations↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2025

Taming the virtual space for incremental full configuration interaction

Incremental full configuration interaction (iFCI) closely approximates the FCI limit with polynomial cost through a many-body expansion of the correlation energy, providing highly accurate total energies within a given basis set. To extend iFCI beyond previous basis set limitations, this work introduces a novel natural orbital (NO) screening approach, incremental NO full configuration interaction (iNO-FCI). By consideration of the importance of virtual orbital selection in the convergence of iFCI, iNO-FCI maximizes the consistency between orbitals selected for each correlated body. iNO-FCI employs a principle of cancellation of errors and ensures that the same set of virtual NOs is used for interdependent terms. Here, this strategy significantly reduces computational cost without compromising precision. Computational savings of up to 95% are demonstrated, allowing access to larger basis sets that were previously computationally prohibitive. iNO-FCI is herein introduced and benchmarked for several difficult test cases involving double-bond dissociation, biradical systems, conjugated π systems, and the spin gap of a Cu-based transition metal complex.

Correlation energy↗
DOE OSTI2025

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
DOE OSTI2025

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗
DOE OSTI2026

Efficient Measurement-Driven Eigenenergy Estimation with Classical Shadows

Quantum algorithms exploiting real-time evolution under a target Hamiltonian have demonstrated remarkable efficiency in extracting key spectral information. However, the broader potential of these methods, particularly beyond ground-state calculations, is underexplored. In this work, we introduce the framework of multiobservable dynamic mode decomposition (MODMD), which combines the observable dynamic mode decomposition (DMD), a measurement-driven eigensolver tailored for near-term implementation, with classical shadow tomography. MODMD leverages random scrambling in the classical shadow technique to construct, with exponentially reduced resource requirements, a signal subspace that encodes rich spectral information. Notably, we replace typical Hadamard-test circuits with a protocol designed to predict low-rank observables, thereby broadening the use of classical shadow tomography for predicting many low-rank observables. We establish theoretical guarantees on the spectral approximation from MODMD, taking into account distinct sources of error. In the ideal case, we prove that the spectral error scales as exp (−Δ⁢𝐸⁢𝑡 max ), where Δ⁢𝐸 is the Hamiltonian spectral gap and 𝑡 max is the maximal simulation time. This analysis provides a rigorous justification of the rapid convergence observed across simulations. To demonstrate the utility of our framework, we consider its application to fundamental tasks, such as determining the low-lying, i.e., ground or excited, energies of representative many-body systems. Our work paves the path for efficient designs of measurement-driven algorithms on near-term and early fault-tolerant quantum devices.

quantum algorithms & computation↗