Search NASA⌕ Search

SEARCH · Search NASA

Results for “many-body interactions”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum Electrodynamics Coupled-Cluster at Scale: High-Performance Implementation for Complex Systems

Coupled-cluster theory (CC) is a highly accurate and versatile method for simulating complex interactions within quantum systems. The extension of CC theory to model mixed electron-photon processes with quantum electrodynamics (QED) has improved our capability to predict cavity-modified chemistry, a field where photons are used as cost-effective and eco-friendly alternatives to catalyze/inhibit chemical reactions. However, calculations with CC methods, even without incorporating QED effects, are often prohibitively expensive. Simulations of larger systems require scalable infrastructures that exist for traditional CC methods but not for QED-CC methods. As such, we present a GPU-enabled, high-performance, open-source implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) within the ExaChem quantum chemistry software package. ExaChem relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure: a parallel heterogeneous tensor library designed to achieve scalable performance on modern heterogeneous supercomputing platforms. Furthermore, we discuss theoretical foundations, algorithmic details, and numerical benchmarks to showcase the larger systems that ExaChem can simulate and how the integration of photonic degrees-of-freedom alters their ground-state properties.

Basis sets↗
DOE OSTI2026

Ferromagnetic Insulator to Metal Transition in Noncentrosymmetric Graphene Nanoribbons

Engineering sublattice imbalance within the unit cell of bottom-up synthesized graphene nanoribbons (GNRs) represents a versatile tool for realizing custom-tailored quantum nanomaterials. The interaction between low-energy zero-modes (ZMs) not only contributes to frontier bands but can form the basis for magnetically ordered phases. Here, we present the bottom-up synthesis of a noncentrosymmetric GNR that places all ZMs on the majority sublattice sites. Scanning tunneling microscopy and spectroscopy reveal that strong electron–electron correlations, leading to the Stoner magnetic instability, drive the system into a ferromagnetically ordered insulating ground state featuring a sizable band gap of E g ∼1.2 eV. At higher temperatures, a chemical transformation induces an insulator-to-metal transition that quenches the ferromagnetic order. Tight-binding (TB), density functional theory, and GW calculations corroborate our experimental observations. Furthermore, this work showcases how control over molecular symmetry, sublattice polarization, and ZM hybridization in bottom-up synthesized nanographenes can open a path to the exploration of many-body physics in rationally designed quantum materials.

Band structure↗
DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Physics-tailored machine learning reveals unexpected physics in dusty plasmas

Dusty plasma is a mixture of ions, electrons, and macroscopic charged particles that is commonly found in space and planetary environments. The particles interact through Coulomb forces mediated by the surrounding plasma, and as a result, the effective forces between particles can be nonconservative and nonreciprocal. Machine learning (ML) models are a promising route to learn these complex forces, yet their structure should match the underlying physical constraints to provide useful insight. Here, we demonstrate and experimentally validate an ML approach that incorporates physical intuition to infer force laws in a laboratory dusty plasma. Trained on 3D particle trajectories, the model accounts for inherent symmetries, nonidentical particles, and learns the effective nonreciprocal forces between particles with exquisite accuracy (R 2 > 0.99). We validate the model by inferring particle masses in two independent yet consistent ways. The model’s accuracy enables precise measurements of particle charge and screening length, identifying large deviations from common theoretical assumptions. Our ability to identify unknown physics from experimental data demonstrates how ML-powered approaches can guide new routes of scientific discovery in many-body systems. Furthermore, we anticipate our ML approach to be a starting point for inferring laws from dynamics in a wide range of many-body systems, from colloids to living organisms.

Science & Technology - Other Topics↗
DOE OSTI2025

Alternating and Gaussian Fermionic Isometric Tensor Network States

Isometric tensor networks in two dimensions enable efficient and accurate study of quantum many-body states, yet the effect of the isometric restriction on the represented quantum states is not fully understood. We address this question in two main contributions. First, we introduce an improved variant of isometric tensor network states (isoTNS) in two dimensions, where the isometric arrows on the columns of the network alternate between pointing upward and downward; hence the name alternating isometric tensor network states. Second, we introduce a numerical tool—the isometric Gaussian fermionic TNS (isoGfTNS)—that incorporates isometric constraints into the framework of Gaussian fermionic tensor network states. We demonstrate in numerous ways that alternating isoTNSs represent many-body ground states of two-dimensional quantum systems significantly better than the original isoTNSs. First, we show that the entanglement in an isoTNS is mediated along the isometric arrows and that alternating isoTNSs mediate entanglement more efficiently than conventional isoTNSs. Second, alternating isoTNSs correspond to a deeper, and thus more representative, sequential-circuit construction of depth 𝒪⁢(𝐿𝑥 ⋅𝐿𝑦) compared to the original isoTNSs of depth 𝒪⁢(𝐿𝑥 +𝐿𝑦). Third, using the Gaussian framework and gradient-based energy minimization, we provide numerical evidence of better bond-dimension scaling and variational energy of alternating isoGfTNSs for ground states of various free-fermionic models, including the Fermi surface, the band insulator, and the 𝑝𝑥 +𝑖⁢𝑝𝑦 mean-field superconductor. Finally, benchmarking on the transverse-field Ising model, we demonstrate that an alternating isoTNS provides substantially improved performance and stability relative to the original isoTNS for the ground-state search algorithm in interacting systems.

Wu, Yantao [Chinese Academy of Sciences, Beijing (↗
DOE OSTI2025

Probing the Kitaev honeycomb model on a neutral-atom quantum computer

Quantum simulations of many-body systems are among the most promising applications of quantum computers. In particular, models based on strongly correlated fermions are central to our understanding of quantum chemistry and materials problems, and can lead to exotic, topological phases of matter. However, owing to the non-local nature of fermions, such models are challenging to simulate with qubit devices. Here we realize a digital quantum simulation architecture for two-dimensional fermionic systems based on reconfigurable atom arrays. We utilize a fermion-to-qubit mapping based on Kitaev’s model on a honeycomb lattice, in which fermionic statistics are encoded using long-range entangled states. We prepare these states efficiently using measurement and feedforward, realize subsequent fermionic evolution through Floquet engineering with tunable entangling gates interspersed with atom rearrangement, and improve results with built-in error detection. Leveraging this fermion description of the Kitaev spin model, we efficiently prepare topological states across its complex phase diagram and verify the non-Abelian spin-liquid phase by evaluating an odd Chern number. We further explore this two-dimensional fermion system by realizing tunable dynamics and directly probing fermion exchange statistics. Finally, we simulate strong interactions and study the dynamics of the Fermi–Hubbard model on a square lattice. These results pave the way for digital quantum simulations of complex fermionic systems for materials science, chemistry and high-energy physics.

atomic and molecular physics↗
DOE OSTI2025

Second-order renormalized Hamiltonian of Yukawa theory

Using the renormalization group procedure for effective particles we calculate the effective Hamiltonians in the theory of a fermion field coupled to a scalar field via the Yukawa interaction. The theory is renormalized by the addition of counterterms. Necessary counterterms are determined by computing matrix elements of the effective Hamiltonian. All calculations are performed up to the second order in the expansion in powers of the coupling constant. Renormalized effective Hamiltonians are well-defined symmetric forms acting in the Fock space as opposed to the renormalized bare Hamiltonian, which is not well defined without regularization. We introduce computational techniques that should streamline higher-order calculations and may be of independent interest.

Ab initio calculations↗
DOE OSTI2025

New Class of Three-Nucleon Forces and Their Implications

Here, we identify a new class of three-nucleon forces that arises in the low-energy effective theory of nuclear interactions including pions. We estimate their contribution to the energy of neutron and nuclear matter and find that it can be as important as the leading-order three-nucleon forces previously considered in the literature. The magnitude of this force is set by the strength of the coupling of pions to two nucleons and is presently not well constrained by experiments. The implications for nuclei, nuclear matter, and the equation of state of neutron matter are briefly discussed.

chiral perturbation theory↗
DOE OSTI2025

Energy transfer between localized emitters in photonic cavities from first principles

Radiative and nonradiative resonant couplings between defects are ubiquitous phenomena in photonic devices used in classical and quantum information technology applications. In this work, we present a first-principles approach to enable quantitative predictions of the energy transfer between defects in photonic cavities, beyond the dipole-dipole approximation and including the many-body nature of the electronic states. As an example, we discuss the energy transfer from a dipolelike emitter to an 𝐹 center in MgO in a spherical cavity. We show that the cavity can be used to controllably enhance or suppress specific spin-flip and spin-conserving transitions. Specifically, we predict that an ∼10–100 enhancement in the resonant energy transfer rate can be gained in the case of the 𝐹 center in MgO at ∼10 nm distances from a dipolar source, using rather moderate cavity with quality factor 𝑄 ∼ 400. We also show that a similar suppression in the transfer rate can be achieved by off-tuning the cavity resonance relative to the emitter transition energy. The framework presented here is general and readily applicable to a wide range of devices where localized emitters are embedded in microspheres, core-shell nanoparticles, and dielectric Mie resonators. Hence, our approach paves the way to predict how to control energy transfer in quantum memories and in ultrahigh-density optical memories, and in a variety of quantum information platforms.

First-principles calculations↗
DOE OSTI2025

Neutron stars and the dense matter equation of state

The past two decades have witnessed tremendous progress in understanding the properties of neutron stars, their maximum mass and radii, and the properties of the dense matter in their cores, made possible by electromagnetic observations of neutron stars and the detection of gravitational waves from their mergers. These observations have provided novel constraints on neutron-star structure that are intimately related to the properties of dense neutron-rich matter described by the nuclear equation of state. Nevertheless, constraining the equation of state over the wide range of densities probed by astrophysical observations is still challenging, as the physics involved is broad and the system spans many orders of magnitude in density. Here, theoretical approaches to calculate and model the neutron-star equation of state in various regimes of densities are reviewed, and the related consequent properties of neutron stars are discussed. How the equation of state at low densities can be calculated from nuclear interactions that are constrained and benchmarked by nuclear experiments is described. Neutron-star observations, with a particular emphasis on information provided by gravitational-wave signals and electromagnetic observations, are reviewed. Finally, future challenges and opportunities in the field are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Finite-range pairing in nuclear density functional theory

Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. Here, we show that a folding radius of about 1 fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.

Nuclear density functional theory↗
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

Efficient use of quantum computers for collider physics

Most observables at particle colliders involve physics at a wide variety of distance scales. Due to asymptotic freedom of the strong interaction, the physics at short distances can be calculated reliably using perturbative techniques, while long distance physics is non-perturbative in nature. Factorization theorems separate the contributions from different scales, allowing to identify the pieces that can be determined perturbatively from those that require non-perturbative information, and if the non-perturbative pieces can be reliably determined, one can use experimental measurements to extract the short distance effects, sensitive to possible new physics. Without the ability to compute the non-perturbative ingredients from first principles one typically identifies observables for which the non-perturbative information is universal in the sense that it can be extracted from some experimental observables and then used to predict other observables. In this paper we argue that the future ability to use quantum computers to calculate non-perturbative matrix elements from first principles will allow to make predictions for observables with non-universal non-perturbative long-distance physics.

Algorithms and Theoretical Developments↗
DOE OSTI2025

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
DOE OSTI2025

Deconstruction of the Anisotropic Magnetic Interactions from Spin‐Entangled Optical Excitations in van der Waals Antiferromagnets

Magneto-optical excitations in antiferromagnetic d systems can originate from a multiplicity of light-spin and spin-spin interactions, as the light and spin degrees of freedom can be entangled. This is exemplified in van der Waals systems with attendant strong anisotropy between in-plane and out-of-plane directions, such as MnPS 3 and NiPS 3 films studied here. The rich interplay between the magnetic ordering and sub-bandgap optical transitions poses a challenge to resolve the mechanisms driving spin-entangled optical transitions, as well as the single-particle bandgap itself. Here, a high-fidelity ab initio theory is applied to find a realistic estimation of the bandgap by elucidating the atom- and orbital-resolved contributions to the fundamental sub-bands. It is further demonstrated that the spin-entangled excitations, observable as photoluminescence and absorption resonances, originate from an on-site spin-flip transition confined to a magnetic atom (Mn or Ni). The evolution of the spin-flip transition in a magnetic field is used to deduce the effective exchange coupling and anisotropy constants.

2D materials↗
DOE OSTI2025

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗
DOE OSTI2025

Two-channel multi-impurity Kondo model: RKKY-induced criticality

We show that the Ruderman-Kittel-Kasua-Yoisida interaction between overscreened spins in two-channel Kondo impurity systems is a relevant perturbation when the number of impurities N is greater than 3 driving the system to a new quantum critical point with anomalous dimensions $\frac{1}{(𝑁+1)}$ for the spin operator and the Sommerfeld coefficient of the specific heat scales as 𝛾 ∼ 𝑇 −$\frac{3}{𝑁+1}$ . The critical point universal properties are relevant to many strong correlation problems, such as impurity placed in a Majorana metal and the multichannel Kondo lattice model of heavy fermion materials. In conclusion, we discuss relevance of our results for cluster DMFT studies of quantum criticality.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Solving reaction dynamics with quantum computing algorithms

The description of quantum many-body dynamics is extremely challenging on classical computers, as it can involve many degrees of freedom. However, the time evolution of quantum states is a natural application for quantum computers that are designed to efficiently perform unitary transformations. Here, in this paper, we study quantum algorithms for response functions, relevant for describing different reactions governed by linear response. We focus on nuclear-physics applications and consider a qubit-efficient mapping on the lattice, which can efficiently represent the large volumes required for realistic scattering simulations. For the case of a contact interaction, we develop an algorithm for time evolution based on the Trotter approximation that scales logarithmically with the lattice size and is combined with quantum phase estimation. We eventually focus on the nuclear two-body system and a typical response function relevant for electron scattering as an example. We also investigate ground-state preparation and examine the total circuit depth required for a realistic calculation and the hardware noise level required to interpret the signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗