Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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

Coupled-Cluster Calculations of Neutrinoless Double- β Decay in Ca 48

We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-$\beta$ decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-$\beta$ decay of $^{48}$Ca with a quenching factor deduced from two-body currents in recent ab initio calculation of the Ikeda sum rule in $^{48}$Ca [Gysbers et al., Nat. Phys. 15, 428 (2019)].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

pyDecay: A CRAM-Based Isotope Decay Solver

This module (pyDecay) implements a Chebyshev Rational Approximation Method (CRAM) for solving isotope decay equations, based on the work of M. Pusa. It provides a numerically stable and efficient method for evaluating the matrix exponential involved in nuclear decay calculations. This implementation of CRAM relies on the incomplete partial factorization (IPF) algorithm published by Pusa [3], with corrections noted by Romano et al.

Skutnik, SteveEugene [Oak Ridge National Laborator↗
DOE OSTI2022

Spectral form factor of a quantum spin glass

It is widely expected that systems which fully thermalize are chaotic in the sense of exhibiting random-matrix statistics of their energy level spacings, whereas integrable systems exhibit Poissonian statistics. In this paper, we investigate a third class: spin glasses. These systems are partially chaotic but do not achieve full thermalization due to large free energy barriers. We examine the level spacing statistics of a canonical infinite-range quantum spin glass, the quantum p-spherical model, using an analytic path integral approach. We find statistics consistent with a direct sum of independent random matrices, and show that the number of such matrices is equal to the number of distinct metastable configurations — the exponential of the spin glass “complexity” as obtained from the quantum Thouless-Anderson-Palmer equations. We also consider the statistical properties of the complexity itself and identify a set of contributions to the path integral which suggest a Poissonian distribution for the number of metastable configurations. Our results show that level spacing statistics can probe the ergodicity-breaking in quantum spin glasses and provide a way to generalize the notion of spin glass complexity beyond models with a semi-classical limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The S-matrix of 2D type 0B string theory. Part II. D-instanton effects

We study the effect of D-instantons on closed string scattering amplitudes in the two-dimensional type 0B string theory from the worldsheet perspective. We find that the contribution from a pair of D-instanton and anti-D-instanton to the closed string reflection amplitude, with a suitable contour prescription for the integration over the D-instanton moduli space, agrees with the corresponding leading non-perturbative corrections in the proposed dual matrix quantum mechanics. This analysis is further extended to thermal observables defined at finite temperature. The infrared divergence in charged processes is understood through the measure factor for charged D-instantons, and can be treated with spacetime dimensional regularization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Comparison of δ 13 C analyses of individual foraminifer ( Orbulina universa ) shells by secondary ion mass spectrometry and gas source mass spectrometry

Rationale: The use of secondary ion mass spectrometry (SIMS) to perform micrometer-scale in situ carbon isotope (δ 13 C) analyses of shells of marine microfossils called planktic foraminifers holds promise to explore calcification and ecological processes. The potential of this technique, however, cannot be realized without comparison to traditional whole-shell δ 13 C values measured by gas source mass spectrometry (GSMS). Methods: Paired SIMS and GSMS δ 13 C values measured from final chamber fragments of the same shell of the planktic foraminifer Orbulina universa are compared. The SIMS–GSMS δ 13 C differences (Δ 13 C SIMS-GSMS ) were determined via paired analysis of hydrogen peroxide-cleaned fragments of modern cultured specimens and of fossil specimens from deep-sea sediments that were either untreated, sonicated, and cleaned with hydrogen peroxide or vacuum roasted. After treatment, fragments were analyzed by a CAMECA IMS 1280 SIMS instrument and either a ThermoScientific MAT-253 or a Fisons Optima isotope ratio mass spectrometer (GSMS). Results: Paired analyses of cleaned fragments of cultured specimens (n = 7) yield no SIMS–GSMS δ 13 C difference. However, paired analyses of untreated (n = 18) and cleaned (n = 12) fragments of fossil shells yield average Δ 13 C SIMS-GSMS values of 0.8‰ and 0.6‰ (±0.2‰, 2 SE), respectively, while vacuum roasting of fossil shell fragments (n = 11) removes the SIMS–GSMS δ 13 C difference. Conclusions: The noted Δ 13 C SIMS-GSMS values are most likely due to matrix effects causing sample–standard mismatch for SIMS analyses but may also be a combination of other factors such as SIMS measurement of chemically bound water. The volume of material analyzed via SIMS is ~10 5 times smaller than that analyzed by GSMS; hence, the extent to which these Δ 13 C SIMS-GSMS values represent differences in analyte or instrument factors remains unclear.

58 GEOSCIENCES↗
DOE OSTI2022

Randomized Cholesky Preconditioning for Graph Partitioning Applications

Graph partitioning has emerged as an area of interest due to its use in various applications in computational research. One way to partition a graph is to solve for the eigenvectors of the corresponding graph Laplacian matrix. This project focuses on the eigensolver LOBPCG and the evaluation of a new preconditioner: Randomized Cholesky Factorization (rchol). This proconditioner was tested for its speed and accuracy against other well-known preconditioners for the method. After experiments were run on several known test matrices, rchol appears to be a better preconditioner for structured matrices. This research was sponsored by National Nuclear Security Administration Minority Serving Institutions Internship Program (NNSA-MSIIP) and completed at host facility Sandia National Laboratories. As such, after discussion of the research project itself, this report contains a brief reflection on experience gained as a result of participating in the NNSA-MSIIP.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Binning is Sinning (Supernova Version): The Impact of Self-calibration in Cosmological Analyses with Type Ia Supernovae

Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia supernovae (SNe Ia) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic uncertainty vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of ~ 1.5 × with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach indeed yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is solely degenerate with the cosmological model in redshift space or whether it can be described by additional correlations between supernova properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the data set, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically classified samples, we must first address the issue that binning is required in currently used photometric analysis methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗
DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

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

Low-energy physics of isotropic spin-1 chains in the critical and Haldane phases

Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

In this work, we study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu L-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped t-t'-J model for a set of t' values. We find evidence for the breakdown of spin-charge separation as |t'| increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.

1D spin chains↗
DOE OSTI2023

Brief Announcement: Communication Optimal Sparse LU Factorization for Planar Matrices

We introduce a new parallel algorithm for solving sparse LU factorization of planar matrices, which commonly arise in the finite element method for 2D PDEs. Existing scalable methods, such as the multifrontal approach with subtree-to-subcube mapping by Gupta et al. [1] and right-looking with 3D mapping by Sao et al. [2] fail to achieve optimal communication costs for these matrices. Our new algorithm combines 3D mapping and subtree-to-subcube mapping to minimize communication costs while allowing trade-offs between extra memory and reduced communication. We demonstrate that our proposed algorithm attains the communication lower bound up to a factor of O(log log n) in the memory-optimal case and up to a factor of O(log P) in the memory-independent case for an n-dimensional planar sparse matrix on P processors.

Sao, Piyush↗
DOE OSTI2023

Advanced Ab Initio Methods for Nuclear Structure (Final Report)

Over the past decade, there has been enormous progress in the description of nuclear structure from first principles, using interactions from Chiral Effective Field Theory that are rooted in Quantum Chromodynamics, the fundamental theory of strong interactions, and many-body methods that solve the Schrödinger equation with systematically improvable approximations. Amongst those are the family of In-Medium Similarity Renormalization Group (IMSRG) framework developed by the PI and his co-workers. While these methods scale polynomially in the size N of the single-particle basis, computational efforts still grows dramatically as we increase the degrees of freedom for the nucleons by relaxing symmetries or introducing continuum couplings (see below), or as we push to improved truncations to provide precise inputs for experimental efforts, in particular in fundamental symmetry searches. In order to address this growing computational cost, one focus area of this award was the exploration of compression and factorization methods. The key to success or failure is the presence of low-rank structures within the matrix elements of NN and 3N interactions or the IMSRG evolution operator, and a means to reformulate the method that will let us exploit them.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

FLAG review 2024

We review lattice results related to pion, kaon, 𝐷-meson, 𝐵-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor 𝑓+⁡(0) arising in the semileptonic 𝐾 → 𝜋 transition at zero momentum transfer, as well as the decay-constant ratio 𝑓 𝐾 ⁡/𝑓 𝜋 and its consequences for the Cabibbo–Kobayashi–Maskawa (CKM) matrix elements 𝑉 𝑢⁢𝑠 and 𝑉 𝑢⁢𝑑 . We review the determination of the 𝐵 𝐾 parameter of neutral kaon mixing as well as the additional four 𝐵 parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for 𝑚 𝑐 and 𝑚 𝑏 as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant 𝛼 𝑠 . We review the determinations of nucleon charges from the matrix elements of both isovector and flavor-diagonal axial, scalar and tensor local quark bilinears, and momentum fraction, helicity moment and the transversity moment from one-link quark bilinears. We also review determinations of scale-setting quantities. Finally, in this review we have added a new section on the general definition of the low-energy limit of the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Control of Excitonic Energy Transfer in RGB Quantum Dot:Polymer Composites for Tunable White Emission

Tint-controlled white light is crucial for both illumination systems and display applications. Here, in this study, we demonstrate solution-processed quantum-dot light-emitting diodes (QD-LEDs) featuring a red–green–blue (RGB) QD–poly(methyl methacrylate) (PMMA) composite emissive layer (EML) for tunable white electroluminescence (EL). In this composite EML, PMMA functions as a dispersion matrix that modulates the interdot spacing (d), thereby controlling Förster resonance energy transfer (FRET) between QDs. By adjustment of the PMMA content, the balance of R, G, and B emissions is controlled, enabling systematic and continuous tuning of the EL color from greenish to reddish white at a fixed RGB ratio and constant driving bias. Time-resolved photoluminescence measurements confirm that the variation in the exciton lifetime with the PMMA fraction is the primary factor for tuning the EL color. Notably, nearly pure white EL with CIE coordinates close to (0.33, 0.33) is achieved using a diluted PMMA matrix without significant degradation of the electrical properties. Our results demonstrate d as an independent design parameter for decoupling color tuning from RGB composition and electrical operation, providing a versatile design framework for high-quality white- or tint-controlled QD-LEDs toward advanced solid-state lighting and display technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Method of Moments Wide Band Adaptive Rational Interpolation Method for High-Quality Factor Resonant Cavities

A new adaptive rational interpolation method is proposed to obtain the wideband frequency response of a resonant cavity simulated with the method of moments (MoM). This interpolation method uses both the Loewner matrix to construct a rational expression for the solution vector of MoM’s matrix system and an error estimator generated by the solution vectors and their derivatives. This error estimator is implemented in the adaptive procedure to gain a minimum set of frequencies and solution vectors required in the interpolation. The resulting set of frequencies and solution vectors is applied to interpolate other system variables, such as shielding effectiveness and input impedance. Here numerical results of a slotted cylindrical cavity supporting high-quality factor resonances are presented, showing that the new rational interpolation method is accurate and efficient in interpolating the complicated resonant response of the solution vector functions.

42 ENGINEERING↗
DOE OSTI2022

Duality and form factors in the thermally deformed two-dimensional tricritical Ising model

The thermal deformation of the critical point action of the 2D tricritical Ising model gives rise to an exact scattering theory with seven massive excitations based on the exceptional E_7 E 7 Lie algebra. The high and low temperature phases of this model are related by duality. This duality guarantees that the leading and sub-leading magnetisation operators, \sigma(x) σ ( x ) and \sigma'(x) σ ′ ( x ) , in either phase are accompanied by associated disorder operators, \mu(x) μ ( x ) and \mu'(x) μ ′ ( x ) . Working specifically in the high temperature phase, we write down the sets of bootstrap equations for these four operators. For \sigma(x) σ ( x ) and \sigma'(x) σ ′ ( x ) , the equations are identical in form and are parameterised by the values of the one-particle form factors of the two lightest \mathbb{Z}_2 ℤ 2 odd particles. Similarly, the equations for \mu(x) μ ( x ) and \mu'(x) μ ′ ( x ) have identical form and are parameterised by two elementary form factors. Using the clustering property, we show that these four sets of solutions are eventually not independent; instead, the parameters of the solutions for \sigma(x)/\sigma'(x) σ ( x ) / σ ′ ( x ) are fixed in terms of those for \mu(x)/\mu'(x) μ ( x ) / μ ′ ( x ) . We use the truncated conformal space approach to confirm numerically the derived expressions of the matrix elements as well as the validity of the \Delta Δ -sum rule as applied to the off-critical correlators. We employ the derived form factors of the order and disorder operators to compute the exact dynamical structure factors of the theory, a set of quantities with a rich spectroscopy which may be directly tested in future inelastic neutron or Raman scattering experiments.

Cortes Cubero, Axel↗
DOE OSTI2023

Irradiation Effects on Stability of δ-UZr2 phase in U-50 wt% Zr Alloy

U-50wt%Zr is a candidate metallic nuclear fuel with potential application in light water reactors due to its excellent thermal properties and high radiation tolerance. The Zr-rich UZr fuels possess greater swelling resistance and fission gas release characteristics compared with U-rich UZr fuels. In this current study, the δ-phase U-50wt%Zr is proton irradiated at room temperature to 1 displacement per atom (dpa) to provide insights on phase stability under irradiation conditions. High resolution characterization of Transmission Electron Microscopy (TEM) and Atom Probe Tomography (APT) characterization techniques are used to elucidate microstructural changes due to irradiation. TEM and APT results show highly oriented bcc β-Zr-rich platelet precipitates nucleating adjacent to α-U phases inside the UZr 2 matrix. Formation of this platelet morphology is characteristic of Widmänstatten structure which can be attributed to a variety of factors such as differences in thermal expansion coefficient between the phases, grain size, alloy composition, and cooling rate. The phases present are distinctly different than those observed through in situ annealing, but irradiation accelerates diffusion and phase separation kinetics. These microstructural changes in U-50wt%Zr are different from those achieved by pure thermodynamic or high temperature heavy ion irradiation experiments. Finally, our work, together with previous ones, highlight the necessity to study the U-Zr phase diagram under non-equilibrium thermodynamics conditions to support this material's deployment as a viable nuclear fuel form.

36 MATERIALS SCIENCE↗