Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix”

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

Density-matrix approach for sequential dissociative double ionization of molecules

Intense femtosecond infrared (IR) laser pulses have been used in recent years to study the breakup dynamics of molecules. However, observables such as kinetic energy release and branching ratios of molecular fragments are often difficult to predict by theory, making information of the molecular dynamics difficult to retrieve. In this work, we develop a simple model for sequential double ionization of molecules based on a density-matrix approach. The model describes tunneling ionization of the neutral and the ion as well as laser couplings between the ionic states simultaneously. Population of different doubly charged states can be obtained at a low computational cost. Furthermore, we applied our model to N 2 and obtained a good agreement on the kinetic energy release spectrum with a previous experiment. This theoretical development could open up the opportunities to use intense short IR laser pulses with coincidence measurement to probe molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Large classes of quantum scarred Hamiltonians from matrix product states

Motivated by the existence of exact many-body quantum scars in the Affleck-Kennedy-Lieb-Tasaki (AKLT) chain, in this work we explore the connection between matrix product state (MPS) wave functions and many-body quantum scarred Hamiltonians. We provide a method to systematically search for and construct parent Hamiltonians with towers of exact eigenstates composed of quasiparticles on top of an MPS wave function. These exact eigenstates have low entanglement in spite of being in the middle of the spectrum, thus violating the strong eigenstate thermalization hypothesis. Using our approach, we recover the AKLT chain starting from the MPS of its ground state, and we derive the most general nearest-neighbor Hamiltonian that shares the AKLT quasiparticle tower of exact eigenstates. We further apply this formalism to other simple MPS wave functions, and derive families of Hamiltonians that exhibit AKLT-like quantum scars. As a consequence, we also construct a scar-preserving deformation that connects the AKLT chain to the integrable spin-1 pure biquadratic model. Finally, we also derive other families of Hamiltonians that exhibit types of exact quantum scars, including a $\textit{U}$(1)-invariant perturbed Potts model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model

We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Role of Coulomb interaction in the phase formation of fcc Ce: Correlation matrix renormalization theory

The effect of electronic Coulomb interaction on the phase formation of fcc Ce lattice is investigated by full ab initio calculations without adjustable Coulomb U and J parameters using the Gutzwiller wavefunction-based correlation matrix renormalization theory (CMRT). Its total energy and pressure as a function of volume agree reasonably well with existing DFT+Gutzwiller calculations and experiments, indicating correct capture of electronic correlation and screening effects within the CMRT formalism. Here, a stable phase is found in line with the experimental α-Ce phase, and a lurking phase is identified supposedly linked with the experimental γ-Ce phase. A criterion based on the local 4f electron charge fluctuation is introduced to confirm the distinct electronic correlation natures of both phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Root-N Krylov-space correction vectors for spectral functions with the density matrix renormalization group

In this work, we propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov-space decomposition to compute the correction vectors. Our approach entails the calculation of the root-N (N=2 is the standard square root) of the Hamiltonian propagator using Krylov-space decomposition and repeating this procedure N times to obtain the actual correction vector. We show that our method greatly alleviates the burden of keeping a large bond dimension at large target frequencies, a problem found with conventional correction-vector DMRG, whereas achieving better computational performance at large N. We apply our method to spin and charge spectral functions of t-J and Hubbard models in the challenging two-leg ladder geometry and provide evidence that the root-N approach reaches a much improved spectral resolution compared to the conventional correction vector.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ab initio study on fcc Pr with correlation matrix renormalization theory

Here, we studied fcc praseodymium (Pr) with the ab initio correlation matrix renormalization theory (CMRT) explicitly calculating the Coulomb interactions among basis orbitals and without using any adjustable parameters to work for strongly correlated electron systems. We calculated its total energy in a paramagnetic ground state and studied the role of the correlated 4f electrons in the system. Good agreement was obtained between CMRT and experiments in the pressure volume dependence of the fcc phase. We also compared the CMRT results against other theoretical methods including local density approximation + dynamical mean-field theory and showed consistent results among them. Moreover, we found the normalized local charge fluctuation of the 4f electrons as the function of volume exhibits a clear slope change at the volume collapse region, indicating a switch in their correlation nature there.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Energy transfer in random-matrix ensembles of Floquet Hamiltonians

Here, we explore the statistical properties of energy transfer in ensembles of doubly driven random-matrix Floquet Hamiltonians based on universal symmetry arguments. The energy-pumping efficiency distribution P⁡($\overline{E}$) is associated with the Hamiltonian parameter ensemble and the eigenvalue statistics of the Floquet operator. For specific Hamiltonian ensembles, P⁡($\overline{E}$) undergoes a transition which cannot be associated with a symmetry breaking of the instantaneous Hamiltonian. The Floquet eigenvalue spacing distribution indicates the considered ensembles constitute generic nonintegrable Hamiltonian families. As a step towards Hamiltonian engineering, we develop a machine-learning classifier to understand the relative parameter importance in resulting high-conversion efficiency. We propose random Floquet Hamiltonians as a general framework to investigate frequency conversion effects in a class of generic dynamical processes beyond adiabatic pumps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ab initio description of bcc iron with correlation matrix renormalization theory

We applied the ab initio spin-polarized correlation matrix renormalization theory to the ferromagnetic state of the bulk bcc iron. We showed that it was capable of reproducing the equilibrium physical properties and the pressure-volume curve in good comparison with experiments. We then focused on the analysis of its local electronic correlations. By exploiting different local fluctuation-related physical quantities as measures of electronic correlation within target orbits, we elucidated the different roles of t 2g and e g states in both spin channels and presented compelling evidence to showcase this distinction in their electronic correlation.

3-dimensional systems↗
DOE OSTI2020

Spectroscopic factors, overlaps, and isospin symmetry from an R -matrix point of view

Background: Spectroscopic factors, overlaps, and isospin symmetry are often used in conjunction with single-particle wave functions for the phenomenological analysis of nuclear structure and reactions. Many differing prescriptions for connecting these quantities to physically relevant asymptotic normalization constants or widths are available in the literature, but their relationship and degree of validity are not always clear. Purpose: This paper derives relationships among the above quantities of interest using well-defined methodology and starting assumptions. Method: R -matrix theory is used as the primary tool to interoperate among the quantities of interest to this work. Particular attention is paid to effects arising from beyond the nuclear surface, where isospin symmetry is strongly violated. Results: Relationships between the quantities of interest are derived. Example applications of these methods to mirror levels in nucleon + 12 C, nucleon + 16 O, and nucleon + 26 Al are presented. A new approach to multilevel mirror symmetry is derived and applied to the first three 2 + states of 18 O and 18 Ne. Conclusions: In this work, the relationship between the quantities of interest is clarified and certain procedures are recommended. It is found that the asymptotic normalization constant of the second 2 + state in 18 Ne deduced from the mirror state in 18 O is significantly larger than found in previous work. This finding has the effect of increasing the 17 F($\textit{p,γ}$) 18 Ne reaction rate in novae.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. In this study, we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measurement of spin density matrix elements in Λ ( 1520 ) photoproduction at 8.2–8.8 GeV

We report on the measurement of spin density matrix elements of the $\Lambda(1520)$ in the photoproduction reaction $\gamma p\rightarrow \Lambda(1520)K^+$, via its subsequent decay to $K^{-}p$. The measurement was performed as part of the GlueX experimental program in Hall D at Jefferson Lab using a linearly polarized photon beam with $E_\gamma =$ 8.2-8.8 GeV. These are the first such measurements in this photon energy range. Results are presented in bins of momentum transfer squared, $-(t-t_\text{0})$. We compare the results with a Reggeon exchange model and determine that natural exchange amplitudes are dominant in $\Lambda(1520)$ photoproduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Density matrix renormalization group description of the island of inversion isotopes F 28 – 33

Recent experiments have confirmed that the neutron-rich isotopes 28,29 F belong to the so-called island of inversion (IOI), a region of the nuclear chart around Z = 10 and N = 20 where nuclear structure deviates from the standard shell model predictions due to deformation and continuum effects. However, while the general principles leading to the IOI are relatively well understood, the details of the low-lying structure of the exotic fluorine isotopes 28–33 F are basically unknown. In this study, we perform large-scale shell model calculations including continuum states to investigate the properties of the neutron-rich isotopes 25–33 F, from a core of 24 O and using an effective two-body interaction with a small number of adjustable parameters in the central and tensor channels. We develop two models adjusted on experimentally confirmed states in 25,26 O and 25–27 F based on different assumptions concerning the positions of the neutron 0d 3/2 and 1p 3/2 shells, and solve the many-body problem using the density matrix renormalization group (DMRG) method for open quantum systems in an sd–fp model space. We obtain the first detailed spectroscopy of 25–33F in the continuum and show how the interplay between continuum effects and deformation explains the recent data on 28,29 F. Several deformed one- and two-neutron halo states are predicted in 29,31 F, and we provide some information about the possible structure of the heaviest fluorine isotopes. We also suggest several experimental studies of interest to constraint models and test the present predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Statistical analysis for the neutrinoless double-β-decay matrix element of 48 Ca

Neutrinoless double-β-decay (0⁢vββ) nuclear matrix elements (NME) are the object of many theoretical calculation methods, and are very important for analysis and guidance of a large number of experimental efforts. However, there are large discrepancies between the NME values provided by different methods. Here, in this paper, we propose a statistical analysis of the 48 Ca 0v⁢ββ NME using the interacting shell model, emphasizing the range of the NME probable values and their correlations with observables that can be obtained from the existing nuclear data. Based on this statistical analysis with three independent effective Hamiltonians, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (0.45–0.95) at 90% confidence level, and a mean value of 0.68.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Predicting the neutrinoless double-β-decay matrix element of 136 Xe using a statistical approach

Calculation of the nuclear matrix elements (NMEs) for double-β decay is of paramount importance for guiding experiments and for analyzing and interpreting the experimental data, especially for the search of the neutrinoless double β decay mode (0v⁢ββ). However, there are currently still large differences between the NME values calculated by different methods, hence a quantification of their uncertainties is very much required. Here, in this paper, we propose a statistical analysis of 0⁢v⁢ββ NME for the 136 Xe isotope, based on the interacting shell model, but using three independent effective Hamiltonians, emphasizing the range of the NMEs' most probable values and its correlations with observables that can be obtained from the existing nuclear data. Consequently, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (1.55–2.65) at 90% confidence level, with a mean value of 1.99 and a standard deviation of 0.37.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Measurement of spin-density matrix elements in ρ ( 770 ) production with a linearly polarized photon beam at E γ = 8.2 – 8.8 GeV

The GlueX experiment at Jefferson Lab studies photoproduction of mesons using linearly polarized 8.5GeV photons impinging on a hydrogen target which is contained within a detector with near-complete coverage for charged and neutral particles. We present measurements of spin-density matrix elements for the photoproduction of the vector meson rho(770). The statistical precision achieved exceeds that of previous experiments for polarized photoproduction in this energy range by orders of magnitude. Here we confirm a high degree of s-channel helicity conservation at small squared four-momentum transfer t and are able to extract the t-dependence of natural and unnatural-parity exchange contributions to the production process in detail. We confirm the dominance of natural-parity exchange over the full t range. We also find that helicity amplitudes in which the helicity of the incident photon and the photoproduced rho(770) differ by two units are negligible for -t<0.5GeV 2 /c 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Probabilistic neural networks for improved analyses with phenomenological R -matrix

Here we present a method for measurement analyses based on probabilistic deep neural networks that provide several advantages over conventional analyses with phenomenological models. These include predicting physical quantities directly from data, the rapid generation of statistically robust uncertainties, and the ability to bypass some parameters that may induce ambiguities and complications in data analysis. As deep learning methods make predictions through “black boxes,” the uncertainty quantification is typically challenging. We use a probabilistic framework that provides thorough uncertainty quantification and is straightforward to follow in practice. With the network architecture based on the Transformer, we demonstrate the current method for predicting nuclear resonance parameters from scattering data using the phenomenological R-matrix model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗