Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum many body models”

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

Dynamical Complexity of Non-Gaussian Many-Body Systems with Dissipation

We characterize the dynamical state of many-body bosonic and fermionic many-body models with intersite Gaussian couplings, on-site non-Gaussian interactions, and local dissipation comprising incoherent particle loss, particle gain, and dephasing. We first establish that, for fermionic systems, if the dephasing noise is larger than the non-Gaussian interactions, irrespective of the Gaussian coupling strength, the system state is a convex combination of Gaussian states at all times. Furthermore, for bosonic systems, we show that if the particle loss and particle gain rates are larger than the Gaussian intersite couplings, the system remains in a separable state at all times. Building on this characterization, we establish that at noise rates above a threshold, there exists a classical algorithm that can efficiently sample from the system state of both the fermionic and bosonic models. Finally, we show that, unlike fermionic systems, bosonic systems can evolve into states that are not convex Gaussian even when the dissipation is much higher than the on-site non-Gaussianity. Similarly, unlike bosonic systems, fermionic systems can generate entanglement even with noise rates much larger than the intersite couplings.

Computational complexity↗
DOE OSTI2024

A precise electron EDM constraint on CP-odd heavy-quark Yukawas

CP-odd Higgs couplings to bottom and charm quarks arise in many extensions of the standard model and are of potential interest for electroweak baryogenesis. These couplings induce a contribution to the electron EDM. The experimental limit on the latter then leads to a strong bound on the CP-odd Higgs couplings. We point out that this bound receives large QCD corrections, even though it arises from a leptonic observable. We calculate the contribution of CP-odd Higgs couplings to the bottom and charm quarks in renormalisation-group improved perturbation theory at next-to-leading order in the strong interaction, thereby reducing the uncertainty to a few percent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Strongly correlated Chern insulators in magic-angle twisted bilayer graphene

Interactions between electrons and the topology of their energy bands can create unusual quantum phases of matter. Most topological electronic phases appear in systems with weak electron–electron interactions. The instances in which topological phases emerge only as a result of strong interactions are rare and mostly limited to those realized in intense magnetic fields. The discovery of flat electronic bands with topological character in magic-angle twisted bilayer graphene (MATBG) has created a unique opportunity to search for strongly correlated topological phases. Here we introduce a local spectroscopic technique using a scanning tunnelling microscope to detect a sequence of topological insulators in MATBG with Chern numbers C = ±1, ±2 and ±3, which form near filling factors of ±3, ±2 and ±1 electrons per moiré unit cell, respectively, and are stabilized by modest magnetic fields. One of the phases detected here (C = +1) was previously observed when the sublattice symmetry of MATBG was intentionally broken by a hexagonal boron nitride substrate, with interactions having a secondary role. We also demonstrate that strong electron–electron interactions alone can produce not only the previously observed phase, but also other unexpected Chern insulating phases in MATBG. The full sequence of phases that we observe can be understood by postulating that strong correlations favour breaking time-reversal symmetry to form Chern insulators that are stabilized by weak magnetic fields. Our findings illustrate that many-body correlations can create topological phases in moiré systems beyond those anticipated from weakly interacting models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Vibrational exciton nanoimaging of phases and domains in porphyrin nanocrystals

Much of the electronic transport, photophysical, or biological functions of molecular materials emerge from intermolecular interactions and associated nanoscale structure and morphology. However, competing phases, defects, and disorder give rise to confinement and many-body localization of the associated wavefunction, disturbing the performance of the material. Here, we employ vibrational excitons as a sensitive local probe of intermolecular coupling in hyperspectral infrared scattering scanning near-field optical microscopy (IR s-SNOM) with complementary small-angle X-ray scattering to map multiscale structure from molecular coupling to long-range order. In the model organic electronic material octaethyl porphyrin ruthenium(II) carbonyl (RuOEP), we observe the evolution of competing ordered and disordered phases, in nucleation, growth, and ripening of porphyrin nanocrystals. From measurement of vibrational exciton delocalization, we identify coexistence of ordered and disordered phases in RuOEP that extend down to the molecular scale. Even when reaching a high degree of macroscopic crystallinity, identify significant local disorder with correlation lengths of only a few nanometers. This minimally invasive approach of vibrational exciton nanospectroscopy and -imaging is generally applicable to provide the molecular-level insight into photoresponse and energy transport in organic photovoltaics, electronics, or proteins.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spacecraft Geometry Effects on Kinetic Impactor Missions

Abstract The DART (Double Asteroid Redirection Test) mission will impact a spacecraft on the secondary (Dimorphos) of the binary asteroid system Didymos in 2022 September, with the goal of altering the orbital period of Dimorphos about Didymos sufficiently to be observed from ground-based observations. Numerical impact modeling is a crucial component in understanding the outcome of the DART experiment, and while many have investigated the effects of target properties, such as material strength and porosity (which remain unknown), an often overlooked factor is the importance of accurately representing the spacecraft itself in such models. Most impact modeling to date has considered simple impactor geometries such as a solid uniform sphere, but in reality the spacecraft is a complex shape full of different components, open spaces, and thin walled structures. At a minimum, a simple solid representation underestimates the surface area of the impact: for a small body such as Dimorphos (approximately 160 m in diameter), the difference between a spacecraft spanning 20 m (including solar arrays) impacting and a sub-1 m idealized shape may be important. In this paper, we compare models impacting high-fidelity models of the spacecraft based on the CAD geometry with various simplified impactors, in order to assess the potential importance of this effect. We find that the difference between the simplest impactor geometries (such as a uniform sphere) and the real spacecraft is measurable, and has an interesting dependence on the material properties of the asteroid itself.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗
DOE OSTI2021

Potential energy surfaces for high-energy N + O 2 collisions

Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O 2 ( 3 Σ$^{–}_{g}$) with N( 4 S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N 4 and N 2 O 2 ). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.

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

Interaction-enhanced nesting in spin-fermion and Fermi-Hubbard models

The spin-fermion (SF) model postulates that the dominant coupling between low-energy fermions in near critical metals is mediated by collective spin fluctuations (paramagnons) peaked at the Néel wave vector, Q N , connecting hot spots on opposite sides of the Fermi surface. It has been argued that strong correlations at hot spots lead to a Fermi surface deformation (FSD) featuring flat regions and increased nesting. This conjecture was confirmed in the perturbative self-consistent calculations when the paramagnon propagator dependence on momentum deviation from Q N is given by χ − 1 ∝ | Δ q | . Using diagrammatic Monte Carlo (diagMC) technique we show that such a dependence holds only at temperatures orders of magnitude smaller than any other energy scale in the problem, indicating that a different mechanism may be at play. Instead, we find that a χ − 1 ∝ | Δ q | 2 dependence yields a robust finite- T scenario for achieving FSD. To link phenomenological and microscopic descriptions, we applied the connected determinant diagMC method to the ( t − t ′ ) Hubbard model and found that at large U / t > 5.5 before the formation of electron and hole pockets (i) the FSD defined as a maximum of the spectral function is not very pronounced; instead, it is the lines of zeros of the renormalized dispersion relation that deforms toward nesting, and (ii) the static spin susceptibility is well described by χ − 1 ∝ | Δ q | 2 . Flat FS regions yield a nontrivial scenario for realizing a non-Fermi liquid state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Beyond-mean-field approaches for nuclear neutrinoless double beta decay in the standard mechanism

Nuclear weak decays provide important probes to fundamental symmetries in nature. A precise description of these processes in atomic nuclei requires comprehensive knowledge on both the strong and weak interactions in the nuclear medium and on the dynamics of quantum many-body systems. In particular, an observation of the hypothetical double beta decay without emission of neutrinos (0νββ) would unambiguously demonstrate the Majorana nature of neutrinos and the existence of the lepton-number-violation process. It would also provide unique information on the ordering and absolute scale of neutrino masses. The next-generation tonne-scale experiments with sensitivity up to 10 28 years after a few years of running will probably provide a definite answer to these fundamental questions based on our current knowledge on the nuclear matrix element (NME), the precise determination of which is a challenge to nuclear theory. Beyond-mean-field approaches have been frequently adapted for the study of nuclear structure and decay throughout the nuclear chart for several decades. In this review, we summarize the status of beyond-mean-field calculations of the NMEs of 0νββ decay assuming the standard mechanism of an exchange of light Majorana neutrinos. In conclusion, the challenges and prospects in the extension and application of beyond-mean-field approaches for 0νββ decay are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

MICCO: An Enhanced Multi-GPU Scheduling Framework for Many-Body Correlation Functions

Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.

Wang, Qihan↗
DOE OSTI2024

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. Lastly, the methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Entanglement Structure of Non-Gaussian States and How to Measure It

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. Here, we present a protocol that constrains quantum states using experimentally measured correlation functions. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest nontrivial expansion quantitatively predicts early time thermalization dynamics, including signaling the onset of quantum chaos indicated by the entanglement Hamiltonian.

Fermi gases↗
DOE OSTI2021

Eigenstate Entanglement: Crossover from the Ground State to Volume Laws

For the typical quantum many-body systems that obey the eigenstate thermalization hypothesis (ETH), we argue that the entanglement entropy of (almost) all energy eigenstates is described by a single crossover function. The ETH implies that the crossover functions can be deduced from subsystem entropies of thermal ensembles and have universal properties. These functions capture the full crossover from the ground-state entanglement regime at low energies and small subsystem size (area or log-area law) to the extensive volume-law regime at high energies or large subsystem size. For critical one-dimensional systems, a universal scaling function follows from conformal field theory and can be adapted for nonlinear dispersions. We use it to also deduce the crossover scaling function for Fermi liquids in d > 1 dimensions. Here, the analytical results are complemented by numerics for large noninteracting systems of fermions in d ≤ 3 dimensions and have also been confirmed for bosonic systems and nonintegrable spin chains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Quenching of the quantum p -strength in light exotic nuclei

One-nucleon spectroscopic overlaps, their strength or spectroscopic factors (SFs), and nucleon removal cross sections for light nuclei in the mass range A≤12 were evaluated using the fully correlated Quantum Monte Carlo (QMC) wave functions (WFs). Harder (AV18+UX) and softer (NV2+3) bare interactions were used providing consistent results. The SFs were also taken from simple Shell Model (SM) calculations. This structure information was incorporated in the standard three-body Faddeev/Alt-Grassberger-Sandhas reaction formalism to evaluate the (p,pN) cross sections. The results further our understanding of the quenching of the quantum p-shell strength obtained from structure and reactions. We have found the ratios of the total QMC sums of SFs to the SM ones to be uniform and ∼ 3/4. The corresponding ratios of the sums Below Particle Threshold (BPT) of SFs, as well as of total cross sections, deviate strikingly from the uniform trend in some special cases. We find these ratios to be close to unit for 9 Li → 8 Li + n and 9 C → 8 B + p. In contrast, they are strongly reduced for 11 C → 10 C + n and 11 B → 10 Be + p mirror transitions, resulting from the quenching of the strength for low-lying 2 + final state transitions. The QMC theoretical cross section BPT for 11 C(p,pn) is about two times smaller than the experimental data.

(p,pN) reactions↗
DOE OSTI2025

Stabilizing Non-Abelian Topological Order Against Heralded Noise via Local Lindbladian Dynamics

An important open question for the current generation of highly controllable quantum devices is understanding which phases can be realized as stable steady states under local quantum dynamics. In this work, we show how robust steady-state phases with both Abelian and non-Abelian mixed-state topological order can be stabilized, in two spatial dimensions, against generic “heralded” noise using active dynamics that incorporate measurement and feedback, modeled as a fully local Lindblad master equation. These topologically ordered steady states are two-way connected to pure topologically ordered ground states using local quantum channels, and preserve quantum information for a time that is exponentially large in the system size. Specifically, we present explicit constructions of families of local Lindbladians for both Abelian (ℤ 2 ) and non-Abelian (𝐷 4 ) topological order whose steady states host mixed-state topological order when the noise is below a threshold strength. As the noise strength is increased, these models exhibit first-order transitions to intermediate mixed-state phases where they encode robust classical memories, followed by (first-order) transitions to a trivial steady state at high noise rates. When the noise is imperfectly heralded, steady-state order disappears but our active dynamics significantly enhances the lifetime of the encoded logical information. To carry out the numerical simulations for the non-Abelian 𝐷 4 case, we introduce a generalized stabilizer tableau formalism that permits efficient simulation of the non-Abelian Lindbladian dynamics.

Monte Carlo methods↗
DOE OSTI2025

One- and two-body current contributions to lepton-nucleus scattering

Modeling lepton-nucleus scattering with the accuracy required to extract neutrino-oscillation parameters from long- and short-baseline experiments necessitates retaining most quantum-mechanical effects. One such effect is the interference between one- and two-body current operators in the transition currents, which has been known to enhance the cross sections, especially in transverse kinematics. Here, we incorporate such interference in the spectral function formalism, which combines relativistic currents and kinematics with an accurate description of the initial target state. Our analysis of lepton-scattering off 12 C demonstrates that interference effects appreciably enhance the transverse electromagnetic response functions and the flux-folded neutrino-nucleus cross section, in both cases, improving the agreement with experimental data. We discuss the impact on the neutrino-oscillation program and the determination of nucleon axial form factors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗