Search NASA⌕ Search

SEARCH · Search NASA

Results for “Many-body calculations”

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

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗
DOE OSTI2026

Chiral interactions and superfluidity in the calcium isotopic chain

We perform ab initio calculations of three-point mass differences in the odd- and even-mass 39−49Ca isotopes to probe nuclear superfluidity via empirical neutron pairing gaps. We also quantify the sensitivity of those gaps to the parameters of the interaction at mean-field level. Recent studies employing accurate chiral nuclear interactions have found these gaps to be too small. We show that experimental values can be reproduced at mean-field level by substantially increasing the attraction of the singlet 𝑆 -wave two-nucleon contact interaction, but doing so induces an unphysical bound state of the dineutron. The sensitivity of these predictions to the full calibration of the nuclear interaction is then studied by performing Bayesian posterior sampling in a delta-full chiral effective field theory at third chiral order. We find that pairing gaps remain largely unaffected, leaving the explanation of nuclear superfluidity as a future task for improved many-body modeling and refined interactions at higher chiral orders.

Hagen, Gaute [ORNL] (ORCID:0000000160191687)↗
DOE OSTI2025

Global optimization of harmonic oscillator basis in covariant density functional theory

The present investigation focuses on the improvement of the accuracy of the description of binding energies within moderately sized fermionic basis. Using the solutions corresponding to infinite fermionic basis it was shown that in the case of meson exchange (ME) covariant energy density functionals (CEDFs) the global accuracy of the description of binding energies in the finite $N_F$ = 16 - 20 bases can be drastically (by a factor ranging from ~3 up to ~9 dependent on the functional and $N_F$) improved by a global optimization of oscillator frequency of the basis. This is a consequence of the unique feature of the ME functionals in which with increasing fermionic basis size fermionic and mesonic energies approach the exact (infinite basis) solution from above and below, respectively. As a consequence, an optimal oscillator frequency $\hbar\omega_0$ of the basis can be defined which provides an accurate reproduction of exact total binding energies by the ones calculated in truncated basis. This leads to a very high accuracy of the calculations in moderately sized $N_F=20$ basis when mass dependent oscillator frequency is used: global rms differences $\delta B_{rms}$ between the binding energies calculated in infinite and truncated bases are only 0.025 MeV and 0.031 MeV for the NL5(Z) and DD-MEZ functionals, respectively. Optimized values of the oscillator frequency $\hbar\omega_0$ are provided for three major classes of CEDFs, i.e. for density dependent meson exchange functionals, nonlinear meson exchange ones and point coupling functionals.

Binding energy & masses↗
DOE OSTI2026

Large-scale calculations of 𝛽-decay rates and implications for 𝑟-process nucleosynthesis

Nuclear 𝛽 decay is a key element of the astrophysical rapid neutron capture process (𝑟 process). In this work, we present state-of-the-art global 𝛽-decay calculations based on the quantified relativistic nuclear energy density functional theory and the deformed proton-neutron quasiparticle random-phase approximation. Our analysis considers contributions from allowed and first-forbidden transitions. We used two point-coupling functionals with carefully calibrated time-odd terms and isoscalar pairing strength. The new calculations display consistent results for both employed functionals, especially near the neutron drip line, suggesting slower 𝛽 decays past the 𝑁=126 neutron shell closure than in commonly used 𝛽-decay models. The new rates, along with the existing rates based on the recent nonrelativistic global calculations, are found to slow down the synthesis of heavy elements in the 𝑟 process and significantly reduce the contribution of neutron-induced fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ionization of Free and Confined Atoms and Ions (Final Technical Report)

Over the almost 22-year duration of this theoretical research project, a great deal of progress has been made in the understanding of the ionization of free and confined atomic and molecular systems along with a recent particular emphasis upon studies of the ultrafast time-delay phenomenon in the photoemission process. Our methods have emphasized the role of many-body interactions and relativistic effects in the various processes we have studied. The advances that we have made over the course of this project have been in the areas nondipole effects in the photoionization process, inner-shell photoemission from negative atomic ions, excited state photoionization of neutral atoms and positive atomic ions, Cooper minima in quadrupole photoionization, relativistic interaction in the photoionization of atoms and negative atomic ions, photoionization of the fullerene molecule C 60 along with the photoionization of atomic and molecular systems confined in the C 60 molecule, confinement resonances in endohedral atom photoemission, interchannel coupling effects in photoionization, Wigner time delay in free and confined atoms and ions, angular effects in Wigner time delay in atomic photoionization, photoionization of superheavy atoms and the influence of core relaxation on the photoionization of atomic systems, among others. In the course of this project, we have had extensive collaboration with experimentalists where we have explained the physics of numerous measurements and pointed the way to new laboratory investigations with our calculations. We have also collaborated with a number of theorists as well. Thus, owing to these many collaborations with scientists not funded by DOE, we have leveraged the DOE support rather considerably. To illustrate the depth and breadth of the theoretical investigations, listed below are the 138 refereed and invited publications that cite DOE support. These publications detail the many contributions made over the period of this project.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Puzzling 𝐵⁡(𝐸⁢2;0 + → 2 + ) strength in the proton dripline nucleus 36 Ca

Recent measurements of the 𝐸⁢2 transition rate from the ground state to the first 2 + excited state of the proton dripline nucleus 36 Ca show an unusual pattern when compared to its isotopic neighbor 38 Ca: despite having a higher 𝐸 𝑥 ⁡(2$^{+}_{1}$) excitation energy, the 𝐵⁡(𝐸⁢2;0$^{+}_{1}$ → 2$^{+}_{1}$) rate in 36 Ca is larger. The question that naturally arises is to what extent this observation can be attributed to the unbound character of the 2$^{+}_{1}$ state. To understand the influence of the continuum space on the low-energy properties of 36 Ca, we carried out Gamow shell model calculations that can account for the continuum coupling effects associated with the occupation of unbound 𝑓⁡𝑝 shells. We found that in the threshold 2 + state, 36 Ca is spatially diffused, which impacts the observed 𝐵⁡(𝐸⁢2) trend.

20 ≤ A ≤ 38↗
DOE OSTI2025

Satellites, core hole excitations, and spin-resolved electronic structure in the spectroscopy of half-metallic CrO 2

Photoelectron satellites—the structures appearing on the low kinetic or high binding-energy side of the “main” or “elastic” photopeak—betray the complex many-body interactions set in motion by the sudden creation of the core hole. In this work, we demonstrate, using the technologically important ferromagnetic half-metal CrO 2 , how such satellites can manifest themselves in other core-level spectroscopies of the material and how they can reveal important details pertinent to its electronic structure. Specifically, we identify a fluorescence satellite in the Cr 𝐿 3 resonant x-ray-emission spectra that radiates at a constant emission energy across the Cr 𝐿 3 x-ray edge with energy ≈1.3 eV above the ordinary valence fluorescence. Here, we provide evidence that this feature arises from the valence recombination of the Cr 2⁢𝑝 core hole “dressed” by the same shakeup charge-transfer process present in both the Cr x-ray photoelectron and the Cr x-ray absorption spectra with its energy uniquely measuring the exchange splitting of the Cr 3⁢𝑑 level. Further analysis of the x-ray emission data reveals three additional features that radiate at constant loss energy that are attributed to combinations of Cr 3⁢𝑑⁢(𝑡 2⁢𝑔 ) → Cr 3⁢𝑑⁢(𝑡 2⁢𝑔 ), charge-transfer O 2⁢𝑝→Cr 3⁢𝑑, and crystal-field Cr⁢ 3⁢𝑑⁡(𝑡 2⁢𝑔 )→Cr⁢ 3⁢𝑑⁡(𝑒 𝑔 ) excitations. These assignments and their energies are supported by density-functional theory calculations, the accuracy of which we demonstrate by hard x-ray valence-photoemission measurements. Atomic multiplet calculations, which include crystal-field effects, help interpret x-ray photoelectron and x-ray absorption spectra of the covalently mixed Cr ion. Resonant Cr K-𝐿 2,3 ⁢𝐿 2,3 Auger-electron emission spectra support a ligand-to-metal nature of the charge-transfer process while highlighting the charge sensitivity differences between photon-in/electron-out and photon-in/photon-out spectroscopies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Nuclear Responses with Neural-Network Quantum States

We introduce a variational Monte Carlo framework that combines neural-network quantum states with the Lorentz integral transform technique to compute the dynamical properties of self-bound quantum many-body systems in continuous Hilbert spaces. While broadly applicable to various quantum systems, including atoms and molecules, in this initial application we focus on the photoabsorption cross section of light nuclei, where benchmarks against numerically exact techniques are available. Our accurate theoretical predictions are complemented by robust uncertainty quantification, enabling meaningful comparisons with experiments. Here, we demonstrate that a relatively simple nuclear Hamiltonian—based on a leading-order pionless EFT expansion and known to accurately reproduce ground-state energies of nuclei with 𝐴 ≤ 40—also provides a reliable description of the photoabsorption cross section.

Ab initio calculations↗
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

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗
DOE OSTI2026

Tunable Electronic Energy Level Alignment and Exciton Diversity in Organic–Inorganic van der Waals Heterostructures

van der Waals stacking of two-dimensional (2D) materials offers a powerful platform for engineering material interfaces with tailored electronic and optical properties. While most van der Waals multilayers have featured inorganic monolayers, incorporating molecular monolayers introduces additional degrees of tunability and functionality. Here, in this study, we investigate hybrid bilayers composed of atomically thin perylene-based molecular crystals interfaced with monolayer transition metal dichalcogenides (TMDs), specifically MoS 2 and WS 2 . Using the ab initio many-body perturbation theory within the GW approximation and the Bethe-Salpeter equation approach, we predict emergent properties beyond those of the isolated constituent systems. Notably, we find substantial renormalization of monolayer molecular crystal band gap due to TMD-induced polarization. Furthermore, by varying the TMD monolayer, we demonstrate tuning of the energy level alignment of the bilayer and subsequent control over a diversity of lowest-energy excitons, which include strongly bound hybrid excitons and long-lived charge-transfer excitons. These findings establish organic-inorganic van der Waals heterostructures as a promising class of materials for tunable optoelectronic devices and quantum excitonic phenomena, expanding the design space for low-dimensional systems.

GW-BSE calculations↗
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

Charge Radii Measurements of Exotic Tin Isotopes in the Proximity of 𝑁 =50 and 𝑁 =82

We report nuclear charge radii for the isotopes 104–134 Sn , measured using two different collinear laser spectroscopy techniques at ISOLDE-CERN. These measurements clarify the archlike trend in charge radii along the isotopic chain and reveal an odd-even staggering that is more pronounced near the 𝑁 =50 and 𝑁 =82 shell closures. The observed local trends are well described by both nuclear density functional theory and valence space in-medium similarity renormalization group calculations. Both theories predict appreciable contributions from beyond-mean-field correlations to the charge radii of the neutron-deficient tin isotopes. The models, however, fall short of reproducing the magnitude of the known 𝐵⁡(𝐸⁢2) transition probabilities, highlighting the remaining challenges in achieving a unified description of both ground-state properties and collective phenomena.

laser spectroscopy↗
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↗
DOE OSTI2025

Applications of the modified Hulthén-Kohn method for bound and scattering states

We adapt the Hulthén–Kohn method suggested by Efros [Phys. Rev. C 99, 034620 (2019)] for calculating various observables in the continuum and discrete spectrum using two-body interactions in single- and coupled channel systems. We explore the convergence of phase shifts and wave functions as well as the location of S-matrix poles which enables obtaining both resonance and bound state parameters. We find that employing a harmonic oscillator basis, together with an interaction smoothing scheme introduced by Gyarmati et al. [Nucl. Phys. A 326, 119 (1979)], and adopting approximate bound-state solutions for the short-range components of basis wave functions lead to good convergence even with restricted oscillator quanta accessible for modern no-core shell model codes. The adapted Efros method will facilitate ab initio many-body nuclear structure applications.

Nuclear reactions↗