Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix theories”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl 2 and OsOCl 2

The van der Waals oxide dichlorides MOX 2 (M=V, Ta, Nb, Ru, and Os; X=halogen element), with different electronic densities, are attracting considerable attention. Ferroelectricity, spin-singlet formation, and orbital-selective Peierls phases were reported in this family with d 1 or d 2 electronic configurations, all believed to be caused by the strongly anisotropic electronic orbital degree of freedom. Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl 2 and OsOCl 2 with d 4 electronic configurations. Different from a previous study using VOI 2 with d 1 configuration, these systems with 4d4 or 5d4 do not exhibit a ferroelectric instability along the a axis. Due to the fully occupied dxy orbital in RuOCl 2 and OsOCl 2 , the Peierls instability distortion disappears along the b axis, leading to an undistorted Immm phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the a axis. For this reason, the materials of our focus can be regarded as “effective one-dimensional” systems even when they apparently have a dominant two-dimensional lattice geometry. The large crystal-field splitting energy (between dxz/yz and dxy orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the J = 0 singlet state in MOCl 2 (M=Ru or Os) with electronic density n = 4, resulting in a spin-1 system. Moreover, we find staggered antiferromagnetic order with π wave vector along the M-O chain direction (a axis) while the magnetic coupling along the b axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the t 2g orbitals for the Os atoms to construct a multiorbital Hubbard model for the M-O chains. Staggered AFM with ↑|-|↓-↑|-|↓ spin structure dominates in our density matrix renormalization group calculations, in agreement with density functional theory calculations. Our results for RuOCl 2 and OsOCl 2 provide guidance to experimentalists and theorists working on this interesting family of oxide dichlorides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Few-nucleon matrix elements in pionless effective field theory in a finite volume

Here, pionless effective field theory in a finite volume (FVEFT $\notπ$ ) is investigated as a framework for the analysis of multinucleon spectra and matrix elements calculated in lattice QCD (LQCD). By combining FVEFT π with the stochastic variational method, the spectra of nuclei with atomic number A ∈ {2, 3} are matched to existing finite-volume LQCD calculations at heavier-than-physical quark masses corresponding to a pion mass m $\notπ$ = 806 MeV , thereby enabling infinite-volume binding energies to be determined using infinite-volume variational calculations. Based on the variational wave functions that are constructed in this approach, the finite-volume matrix elements of various local operators are computed in FVEFT $\notπ$ and matched to LQCD calculations of the corresponding QCD operators in the same volume, thereby determining the relevant one- and two-body effective field theory counterterms and enabling an extrapolation of the LQCD matrix elements to infinite volume. As examples, the scalar, tensor, and axial matrix elements are considered, in addition to the magnetic moments and the isovector longitudinal momentum fraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dynamical coupled-channel models for hadron dynamics

Dynamical coupled-channel (DCC) approaches parametrize the interactions and dynamics of two and more hadrons and their response to different electroweak probes. The inclusion of unitarity, three-body channels, and other properties from scattering theory allows for a reliable extraction of resonance spectra and their properties from data. Here, we review the formalism and application of the ANL-Osaka, the Juelich-Bonn-Washington, and other DCC approaches in the context of light baryon resonances from meson, (virtual) photon, and neutrino-induced reactions, as well as production reactions, strange baryons, light mesons, heavy meson systems, exotics, and baryon-baryon interactions. Finally, we also provide a connection of the formalism to study finite-volume spectra obtained in Lattice QCD, and review applications involving modern statistical and machine learning tools.

Amplitude analysis↗
DOE OSTI2026

Thermal nature of confining strings

We investigate the quantum statistical properties of the confining string connecting a static fermion-antifermion pair in the massive Schwinger model. By analyzing the reduced density matrix of the subsystem located in between the fermion and antifermion, we demonstrate that as the interfermion separation approaches the string-breaking distance, the overlap between the microscopic density matrix and an effective thermal density matrix exhibits a pronounced, narrow peak, approaching unity at the onset of string breaking. This behavior reveals that the confining flux tube evolves toward a genuinely thermal state as the separation between the charges grows, even in the absence of an external heat bath. In other words, one cannot tell whether a reduced state of the subsystem arises from a surrounding heat bath or from entanglement with the rest of the system. The entanglement spectrum near the critical string-breaking distance exhibits a rapid transition from the dominance of a single state describing the confining electric string towards a strongly entangled state containing virtual fermion-antifermion pairs. Our findings establish a quantitative link between confinement, entanglement, and emergent thermality, and suggest that string breaking corresponds to a microscopic thermalization transition within the flux tube.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Heisenberg spin chain as a worldsheet coordinate for lightcone quantized string

Although the energy spectrum of the Heisenberg spin chain on a circle defined by H = 1 4 Σ k = 1 M σ k x σ k + 1 x + σ k y σ k + 1 y + Δ σ k z σ k + 1 z is well known for any fixed M, the boundary conditions vary according to whether Mϵ 4N + r, where r = -1,0,1,2, and also according to the parity of the number of overturned spins in the state, In string theory all these cases must be allowed because interactions involve a string with M spins breaking into strings with M 1 < M and M - M 1 spins (or vice versa). We organize the energy spectrum and degeneracies of H in the case Δ = 0 where the system is equivalent to a system of free fermions. In spite of the multiplicity of special cases, in the limit M → ∞ the spectrum is that of a free compactified worldsheet field. Such a field can be interpreted as a compact transverse string coordinate x(σ) ≡ x(σ) + R 0 . We construct the bosonization formulas explicitly in all separate cases, and for each sector give the Virasoro conformal generators in both fermionic and bosonic formulations. Furthermore from calculations in the literature for selected classes of excited states, there is strong evidence that the only change for Δ ≠ 0 is a change in the compactification radius R 0 → R Δ . As Δ→ -1 this radius goes to infinity, giving a concrete example of noncompact space emerging from a discrete dynamical system. Finally we apply our work to construct the three string vertex implied by a string whose bosonic coordinates emerge from this mechanism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

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

Diabatic Valence-Hole States in the C 2 Molecule: “Putting Humpty Dumpty Together Again”

Despite the long history of spectroscopic studies of the C$_2$ molecule, fundamental questions about its chemical bonding are still being hotly debated. The complex electronic structure of C$_2$ is a consequence of its dense manifold of near-degenerate, low-lying electronic states. A global multi-state diabatic model is proposed here to disentangle the numerous configuration interactions within four symmetry manifolds of C$_2$ ($^{1}\Pi_g$, $^{3}\Pi_g$, $^{1}\Sigma_u^+$, and $^{3}\Sigma_u^+$). The key concept of our model is the existence of two "valence-hole" configurations, $2\sigma_g^22\sigma_u^11\pi_{u}^33\sigma_g^2$ for $^{1,3}\Pi_g$ states and $2\sigma_g^22\sigma_u^11\pi_{u}^43\sigma_g^1$ for $^{1,3}\Sigma_u^+$ states that derive from $3\sigma_g\leftarrow2\sigma_u$ electron promotion. The lowest-energy state from each of the four C$_2$ symmetry species is dominated by this type of valence-hole configuration at its equilibrium internuclear separation. Finally, as a result of their large binding energy (nominal bond order of 3) and correlation with the 2s$^2$2p$^2$+2s2p$^3$ separated-atom configurations, the presence of these valence-hole configurations has a profound impact on the $global$ electronic structure and unimolecular dynamics of C$_2$.

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

Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding

The optical spectra of neutral oxygen vacancies (F 0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. Here, to estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet–singlet decay pathway.

embedding↗
DOE OSTI2025

Strongly Correlated States of Transition Metal Spin Defects: The Case of an Iron Impurity in Aluminum Nitride

We investigate the electronic properties of an exemplar transition metal impurity in an insulator, with the goal of accurately describing strongly correlated defect states. Here, we consider iron in aluminum nitride, a material of interest for hybrid quantum technologies, and we carry out calculations with quantum embedding methods, density matrix embedding theory (DMET) and quantum defect embedding theory (QDET), and with spin-flip time-dependent density functional theory (TDDFT). We show that both DMET and QDET accurately describe the ground state and low-lying excited states of the defect and that TDDFT yields photoluminescence spectra in agreement with experiments. In addition, we provide a detailed discussion of the convergence of our results as a function of the active space used in the embedding methods, thus defining a protocol to obtain converged data directly comparable with experiments.

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

New Apparatus for Neutron Capture Measurements on Extra Small Radioactive Samples: The DICER Instrument at LANSCE

We report the neutron capture cross section, the probability per unit area a neutron is absorbed by a nucleus followed by the emission of gamma ray(s), is an important quantity in several fields such as astrophysics, nuclear criticality safety, radiochemical diagnostics, nuclear medicine, nuclear forensics and nuclear security. Accurate knowledge of how neutrons interact with matter is a crucial component to understand physical processes in depth. A great example is the understanding of how our universe was created, a quest that requires accurate modeling in which neutron capture cross sections are a key ingredient

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Hardness of Observing Strong-to-Weak Symmetry Breaking

Spontaneous symmetry breaking (SSB) is the cornerstone of our understanding of quantum phases of matter. Recent works have generalized this concept to the domain of mixed states in open quantum systems, where symmetries can be realized in two distinct ways dubbed strong and weak. Novel intrinsically mixed phases of quantum matter can then be defined by the spontaneous breaking of strong symmetry down to weak symmetry. However, proposed order parameters for strong-to-weak SSB (based on mixed-state fidelities or purities) seem to require exponentially many copies of the state, raising the question: Is it possible to efficiently detect strong-to-weak SSB, in general? In this work, we answer this question negatively in the paradigmatic cases of ℤ2 and U(1) symmetries. We construct ensembles of pseudorandom mixed states that do not break the strong symmetry yet are computationally indistinguishable from states that do. This rules out the existence of efficient state-agnostic protocols to detect strong-to-weak SSB.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Black Hole Airy Tail

In Jackiw-Teitelboim (JT) gravity, which is dual to a random matrix ensemble, the annealed entropy differs from the quenched entropy at low temperatures and goes negative. However, computing the quenched entropy in JT gravity requires a replica limit that is poorly understood. To circumvent this, we define an intermediate quantity called the semiquenched entropy, which has the positivity properties of the quenched entropy, while requiring a much simpler replica trick. We compute this in JT gravity in different regimes using (i) a bulk calculation involving wormholes corresponding to the Airy limit of the dual matrix integral and (ii) a boundary calculation involving one-eigenvalue instanton saddles proposed by Hernández-Cuenca, demonstrating consistency between these two calculations in their common regime of validity. We also clarify why similar one-eigenvalue instanton saddles cannot be used to compute the quenched entropy due to a breakdown of the saddle-point approximation for the one-eigenvalue instanton in the replica limit. Our results show how to use the gravitational path integral to prove that black holes in JT gravity have isolated ground states and to study their properties.

Gauge-gravity dualities↗
DOE OSTI2025

Deviations from the Porter-Thomas Distribution due to Nonstatistical 𝛾 Decay below the 150 Nd Neutron Separation Threshold

We introduce a new method for the study of fluctuations of partial transition widths based on nuclear resonance fluorescence experiments with quasimonochromatic linearly polarized photon beams below particle separation thresholds. It is based on the average branching of decays of 𝐽=1 states of an even-even nucleus to the 2$^{+}_{1}$ state in comparison to the ground state. Between 5 and 7 MeV, a constant average branching ratio for 𝛾 decays from 1 − states of 0.490(16) is observed for the nuclide 150 Nd. Assuming 𝜒 2 -distributed partial transition widths, this average branching ratio is related to a degree of freedom of 𝜈 = 1.93⁢(12), rejecting the validity of the Porter-Thomas distribution, requiring 𝜈 = 1. The observed deviation can be explained by nonstatistical effects in the 𝛾-decay behavior with contributions in the range of 9.4(10)% up to 94(10)%.

150 ≤ A ≤ 189↗
DOE OSTI2023

Multipartite entanglement in the one-dimensional spin- 1 2 Heisenberg antiferromagnet

Multipartite entanglement refers to the simultaneous entanglement between multiple subsystems of a many-body quantum system. While multipartite entanglement can be difficult to quantify analytically, it is known that it can be witnessed through the quantum Fisher information (QFI), a quantity that can also be related to dynamical Kubo response functions. In this work, we first show that the finite temperature QFI can generally be expressed in terms of a static structure factor of the system, plus a correction that vanishes as T→0. This implies that the static structure factor witnesses multipartite entanglement near quantum critical points at temperatures below a characteristic energy scale of the system. Therefore, in systems with a known static structure factor, we can deduce finite temperature scaling of multipartite entanglement and low temperature entanglement depth without knowledge of the full dynamical response function of the system. This is particularly useful to study 1D quantum critical systems in which sub-power-law divergences can dominate entanglement growth, where the conventional scaling theory of the QFI breaks down. Furthermore, the 1D spin- 1 2 antiferromagnetic Heisenberg model is an important example of such a system, and we show that multipartite entanglement in the Heisenberg chain diverges nontrivially as ~ln(1/T) 3/2 . We verify these predictions with calculations of the QFI using conformal field theory and matrix product state simulations. Finally, we discuss the implications of our results for experiments to probe entanglement in quantum materials, comparing to neutron scattering data in KCuF 3 , a material well described by the Heisenberg chain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Slow thermalization and subdiffusion in $U$(1) conserving Floquet random circuits

Random quantum circuits are paradigmatic models of minimally structured and analytically tractable chaotic dynamics. We study a family of Floquet unitary circuits with Haar random U(1) charge conserving dynamics; the minimal such model has nearest-neighbor gates acting on spin-1/2 qubits, and a single layer of even/odd gates repeated periodically in time. We find that this minimal model is not robustly thermalizing at numerically accessible system sizes, and displays slow subdiffusive dynamics for long times. We map out the thermalization dynamics in a broader parameter space of charge conserving circuits, and understand the origin of the slow dynamics in terms of proximate localized and integrable regimes in parameter space. In contrast, we find that small extensions to the minimal model are sufficient to achieve robust thermalization; these include (i) increasing the interaction range to three-site gates, (ii) increasing the local Hilbert space dimension by appending an additional unconstrained qubit to the conserved charge on each site, or (iii) using a larger Floquet period comprised of two independent layers of gates. Furthermore, our results should inform future studies of charge conserving circuits, which are relevant for a wide range of topical theoretical questions.

1-dimensional spin chains↗