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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Exact results for nonlinear Drude weights in the spin-1/2 XXZ chain

Nonlinear Drude weight (NLDW) is a generalization of the linear Drude weight, which characterizes the nonlinear transport in quantum many-body systems. We investigate these weights for the spin-1/2 XXZ chain in the critical regime. The effects of the Dzyaloshinskii-Moriya interaction and an external magnetic field are also studied. Solving the Bethe equations numerically, we obtain these weights for very large system sizes and identify parameter regimes where the weights diverge in the thermodynamic limit. These divergences appear in all the orders studied in this Letter and can be regarded as a generic feature of the NLDWs. We study the origin of these divergences and reveal that they result from nonanalytic finite-size corrections to the ground-state energy. As a result, we compute closed-form expressions for several weights in the thermodynamic limit and find excellent agreement with the numerical results.

1-dimensional systems↗
DOE OSTI2020

Mixed-State Entanglement from Local Randomized Measurements

In this work, we propose a method for detecting bipartite entanglement in a many-body mixed state based on estimating moments of the partially transposed density matrix. The estimates are obtained by performing local random measurements on the state, followed by postprocessing using the classical shadows framework. Our method can be applied to any quantum system with single-qubit control. We provide a detailed analysis of the required number of experimental runs, and demonstrate the protocol using existing experimental data [Brydges et al., Science 364, 260 (2019)].

1-dimensional spin chains↗
DOE OSTI2021

Distinguishing localization from chaos: Challenges in finite-size systems

Highlights: • Provides an overview of the current understanding of the many-body localization phase transition. • Discusses the subtleties of finite-size scaling near this transition • Assesses the implications of these subtleties for numerical studies of spin chains. • Explores scaling of diagnostics in models with known localization transitions. • Presents suggestions for future numerical work. We re-examine attempts to study the many-body localization transition using measures that are physically natural on the ergodic/quantum chaotic regime of the phase diagram. Using simple scaling arguments and an analysis of various models for which rigorous results are available, we find that these measures can be particularly adversely affected by the strong finite-size effects observed in nearly all numerical studies of many-body localization. This severely impacts their utility in probing the transition and the localized phase. In light of this analysis, we discuss a recent study (Šuntajs et al., 2020) of the behaviour of the Thouless energy and level repulsion in disordered spin chains, and its implications for the question of whether MBL is a true phase of matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

Straintronics with single-layer MoS 2 : A quantum Monte Carlo study

Using state-of-the-art quantum Monte Carlo (QMC) methods, we study straintronic properties of a single MoS 2 monolayer. 2D MoS 2 is a quintessential straintronic material for which many experiments have been performed. First, we determine the equilibrium atomic structure which is not known experimentally and is strictly needed to correctly determine the straintronic properties. That enables us to precisely analyze the quasiparticle band gaps for any applied biaxial strain, which we describe by a bivariate paraboloid function of lattice constant and internal structural parameter. Using the fixed-node QMC calculations fitted by analytical formulas, we localize the following excited state crossings between the direct, K → K , and indirect Γ → X and K → K /2 excitations. Based on this highly accurate many-body treatment, we predict a gauge factor of 136 meV/% for the K → K transition and a fairly narrow window of ≈ 2.8 % from compressive to mildly tensile strains, accounting for only ≈ 0.3 eV band gap change maintaining the direct character of the gap. Consequently, we suggest that, compared to other 2D materials, such as phosphorene, there is only a limited straintronic tuneability in this material often studied for its straintronic properties. QMC results are compared to results of standard DFT modeling, which reveal insights into the corresponding inaccuracies and therefore open a window for educated use of rapid DFT approaches. Published by the American Physical Society 2024

Huang, Y.↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2020

Transient Studies of Nonequilibrium Transport in Two-Dimensional Semiconductors

The overarching objective of research under DOE support has been to investigate emergent carrier phenomena in different two-dimensional (2D) systems, under conditions where they are driven far from equilibrium. The description of nonequilibrium transport has long represented one of the most enduring challenges in condensed-matter physics, in large part due to the need to invoke the role of complex many-body phenomena (electron-phonon and electron-electron coupling) in this regime. In this project, we have explored aspects of this problem in a variety of different two-dimensional systems. These include the modulation-doped two-dimensional electrons gas, or 2DEG, in GaAs/AlGaAs heterostructures; the electron and hole systems in graphene and its related hetero-structures, and; the strongly-confined electron layers in transition-metal dichalcogenide (MoS 2 and WS 2 ) and trichalcogenide (TiS 3 ) materials. Significant understanding has been achieved in a number of different areas, as we elaborate upon in subsequent sections of this proposal: Nanoscale spectroscopy of defect-related barriers in MoS2 Negative differential conductance due to inter-valley scattering in WS 2 under high fields Incipient localization due to scattering from frozen-phonon disorder in quantum point contacts High-field transient transport in graphene/hexagonal BN heterostructures Universal voltage scaling of quantum corrections in graphene under nonequilibrium A defining feature of much of this research has been an approach to problem solving in which novel experimental techniques (such as transient electrical pulsing and scanning-probe microscopy) are combined with the development (in collaboration with colleagues from condensed-matter theory) of sophisticated models, a strategy that we propose to continue in our future research. In recent years, we have developed an increasing capacity for handling atomically-thin materials, and their van der Waals systems, a capacity that is reflected in the balance of topics listed above.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

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

Neural Network Water Model Based on the MB-Pol Many-Body Potential

The MB-pol many-body potential accurately predicts many properties of water, including cluster, liquid phase, and vapor–liquid equilibrium properties, but its high computational cost can make applying it in large-scale simulations quite challenging. Here, in order to address this limitation, we developed a “deep potential” neural network (DPMD) model based on the MB-pol potential for water. We find that a DPMD model trained on mostly liquid configurations yields a good description of the bulk liquid phase but severely underpredicts vapor–liquid coexistence densities. By contrast, adding cluster configurations to the neural network training set leads to a good agreement for the vapor coexistence densities. Liquid phase densities under supercooled conditions are also represented well, even though they were not included in the training set. These results confirm that neural network models can combine accuracy and transferability if sufficient attention is given to the construction of a representative training set for the target system.

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

Krylov complexity in mixed phase space

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

chaos & nonlinear dynamics↗
DOE OSTI2025

Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation

Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025

Wang, Ruoxi↗
DOE OSTI2022

Second-order magnetic responses in quantum magnets: Magnetization under ac magnetic fields

Here we investigate second-order magnetic responses of quantum magnets against ac magnetic fields. We focus on the case where the z component of the spin is conserved in the unperturbed Hamiltonian and the driving field is applied in the xy plane. We find that linearly polarized driving fields induce a second-harmonic response, while circularly polarized fields generate only a zero-frequency response, leading to a magnetization with a direction determined by the helicity. Employing an unbiased numerical method, we demonstrate the nonlinear magnetic effect driven by the circularly polarized field in the XXZ model and show that the magnitude of the magnetization can be predicted by the dynamical spin structure factor in the linear response regime.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generalized Hydrodynamics: A Perspective

Conventional hydrodynamics describes systems with few long-lived excitations. In one dimension, however, many experimentally relevant systems feature a large number of long-lived excitations and conserved quantities even at high temperature, because they are proximate to integrable limits. Such models cannot be treated using conventional hydrodynamics. The framework of generalized hydrodynamics (GHD) was recently developed to treat the dynamics of one-dimensional models: It combines ideas from integrability, hydrodynamics, and kinetic theory to come up with a quantitative theory of transport. GHD has successfully settled several long-standing questions about one-dimensional transport; it has also been leveraged to study dynamical questions beyond the transport of conserved quantities and to systems that are not integrable. In this article, we introduce the main ideas and predictions of GHD, survey some of the most recent theoretical extensions and experimental tests of the GHD framework, and discuss some open questions in transport that the GHD perspective might elucidate. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Directional Superradiance in a Driven Ultracold Atomic Gas in Free Space

Ultracold atomic systems are among the most promising platforms that have the potential to shed light on the complex behavior of many-body quantum systems. One prominent example is the case of a dense ensemble illuminated by a strong coherent drive while interacting via dipole-dipole interactions. Despite being subjected to intense investigations, this system retains many open questions. A recent experiment carried out in a pencil-shaped geometry [Ferioli Nat. Phys. 19, 1345 (2023)] has reported measurements that have seemed consistent with the emergence of strong collective effects in the form of a “superradiant” phase transition in free space, when looking at the light-emission properties in the forward direction. Motivated by the experimental observations, we carry out a systematic theoretical analysis of the steady-state properties of the system as a function of the driving strength and atom number N . We observe signatures of collective effects in the weak-driving regime, which disappear with increasing drive strength as the system evolves into a single-particle-like mixed state comprised of randomly aligned dipoles. Although the steady state features some similarities to the reported superradiant-to-normal nonequilibrium transition, also known as cooperative resonance fluorescence, we observe significant qualitative and quantitative differences, including a different scaling of the critical drive parameter (from N to N ). We validate the applicability of a mean-field treatment to capture the steady-state dynamics under currently accessible conditions. Furthermore, we develop a simple theoretical model that explains the scaling properties by accounting for interaction-induced inhomogeneous effects and spontaneous emission, which are intrinsic features of interacting disordered arrays in free space. Published by the American Physical Society 2024

Agarwal, Sanaa (ORCID:0000000336863804)↗
DOE OSTI2025

Equation-of-motion internally contracted multireference unitary coupled-cluster theory

The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Unitary coupled-channel three-body amplitude with pions and kaons

Three-body dynamics above threshold is required for the reliable extraction of many amplitudes and resonances from experiment and lattice QCD. The S-matrix principle of unitarity can be used to construct dynamical coupled-channel approaches in which three particles scatter off each other, rearranging two-body subsystems by particle exchange. This paper reports the development of a three-body coupled-channel, amplitude including pions and kaons. The unequal-mass amplitude contains two-body S- and P-wave subsystems (“isobars”) of all isospins, 𝐼 = 0,1/2,1,3/2,2 , and it also allows for transitions within a given isobar. The 𝑓 0 ⁡(500)⁢(𝜎),𝑓 0 ⁡(980),𝜌⁡(700),𝐾$^{*}_{0}$⁡(700)⁢(𝜅), and 𝐾*⁡(892) resonances are included, apart from repulsive isobars. Different methods to evaluate the amplitude for physical momenta are discussed. Production amplitudes for 𝑎 1 quantum numbers are shown as a proof of principle for the numerical implementation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Complex paths around the sign problem

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of physics, from QCD and nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half filling, and dynamical and nonequilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called sign problem. Here, a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it is reviewed. The mathematical ideas underpinning this approach, as well as the algorithms developed thus far, are described together with nontrivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗