SEARCH · Engineering Papers

Results for “Matrix 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Genomic selection in algae with biphasic lifecycles: A Saccharina latissima (sugar kelp) case study

Introduction Sugar kelp ( Saccharina latissima ) has a biphasic life cycle, allowing selection on both thediploid sporophytes (SPs) and haploid gametophytes (GPs). Methods We trained a genomic selection (GS) model from farm-tested SP phenotypic data and used a mixed-ploidy additive relationship matrix to predict GP breeding values. Topranked GPs were used to make crosses for further farm evaluation. The relationship matrix included 866 individuals: a) founder SPs sampled from the wild; b) progeny GPs from founders; c) Farm-tested SPs crossed from b); and d) progeny GPs from farm-tested SPs. The complete pedigree-based relationship matrix was estimated for all individuals. A subset of founder SPs ( n = 58) and GPs ( n = 276) were genotyped with Diversity Array Technology and whole genome sequencing, respectively. We evaluated GS prediction accuracy via cross validation for SPs tested on farm in 2019 and 2020 using a basic GBLUP model. We also estimated the general combining ability (GCA) and specific combining ability (SCA) variances of parental GPs. A total of 11 yield-related and morphology traits were evaluated. Results The cross validation accuracies for dry weight per meter ( r ranged from 0.16 to 0.35) and wet weight per meter ( r ranged 0.19 to 0.35) were comparable to GS accuracy for yield traits in terrestrial crops. For morphology traits, cross validation accuracy exceeded 0.18 in all scenarios except for blade thickness in the second year. Accuracy in a third validation year (2021) was 0.31 for dry weight per meter over a confirmation set of 87 individuals. Discussion Our findings indicate that progress can be made in sugar kelp breeding by using genomic selection.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Status of SAS4A/SASSYS-1 Software Development and Application (FY2024)

SAS4A/SASSYS-1 is a simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. With its origin as SAS1A in the late 1960s, the SAS series of codes has been under continuous use and development for over fifty years. It has been identified as a critical element of safety analysis capabilities for the U.S. Department of Energy and is utilized within industry to perform the transient safety analyses required to support the licensing of Liquid Metal-cooled Fast Reactors (LMFRs). This report summarizes the code development and update activities carried out during FY2024. In FY2024, programmatic activities focused on key improvements to software useability, such as enhanced user interfaces for reactivity feedback modeling, improvements in stability/useability of the Code Manual, and improvements to the acceptance testing infrastructure, including automation of acceptance testing and generation of the Acceptance Testing Report. To support end user applications, an open training was held, a semi-public forum was maintained, and a practical benchmarking and validation matrix was developed which allowed limitations of existing testing capabilities to be assessed. The existing fuel models were also enhanced with improved modeling capabilities and testing for the oxide and annular fuel models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

In this work, we present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system’s governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh–Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Comparative study of He 6 β -decay based on different similarity-renormalization-group evolved chiral interactions

We report on a study of the Gamow-Teller matrix element contributing to ${}^6$He $\beta$-decay with similarity renormalization group (SRG) versions of momentum- and configuration-space two-nucleon interactions. These interactions are derived from two different formulations of chiral effective field theory ($\chi$EFT) -- without and with the explicit inclusion of $\Delta$-isobars. We consider evolution parameters $\Lambda_{\rm SRG}$ in the range between 1.2 and 2.0 fm$^{-1}$ and, for the $\Delta$-less case, also the unevolved (bare) interaction. The axial current contains one- and two-body terms, consistently derived at tree level (no loops) in the two distinct $\chi$EFT formulations we have adopted here. The ${}^6$He and ${}^6$Li ground-state wave functions are obtained from hyperspherical-harmonics (HH) solutions of the nuclear many-body problem. In $A\,$=$\,6$ systems, the HH method is limited at present to treat only two-body interactions and non-SRG evolved currents. Our results exhibit a significant dependence on $\Lambda_{\text{SRG}}$ of the contributions associated with two-body currents, suggesting that a consistent SRG-evolution of these is needed in order to obtain reliable estimates. We also show that the contributions from one-pion-exchange currents depend strongly on the model (chiral) interactions and on the momentum- or configuration-space cutoffs used to regularize them. These results might prove helpful in clarifying the origin of the sign difference recently found in No-Core-Shell-Model and Quantum Monte Carlo calculations of the ${}^6$He Gamow-Teller matrix element.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Tensor renormalization group for fermions

Abstract We review the basic ideas of the tensor renormalization group method and show how they can be applied for lattice field theory models involving relativistic fermions and Grassmann variables in arbitrary dimensions. We discuss recent progress for entanglement filtering, loop optimization, bond-weighting techniques and matrix product decompositions for Grassmann tensor networks. The new methods are tested with two-dimensional Wilson–Majorana fermions and multi-flavor Gross–Neveu models. We show that the methods can also be applied to the fermionic Hubbard model in 1+1 and 2+1 dimensions.

Physics↗

DOE OSTI2021

An SYK-inspired model with density–density interactions: Spectral & wave function statistics, Green’s function and phase diagram

The Sachdev–Ye–Kitaev (SYK) model is a rare example of a strongly-interacting system that is analytically tractable. Tractability arises because the model is largely structureless by design and therefore artificial: while the interaction is restricted to two-body terms, interaction matrix elements are “randomized” and therefore the corresponding interaction operator does not commute with the local density. Unlike conventional density–density-type interactions, the SYK-interaction is, in this sense, not integrable. We here investigate a variant of the (complex) SYK model, which restores this integrability. It features a randomized single-body term and a density–density-type interaction. We present numerical investigations suggesting that the model exhibits two integrable phases separated by several intermediate phases including a chaotic one. The chaotic phase carries several characteristic SYK-signatures including in the spectral statistics and the frequency scaling of the Green’s function and therefore should be adiabatically connected to the non-Fermi liquid phase of the original SYK model. Thus, our model Hamiltonian provides a bridge from the SYK-model towards microscopic realism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Interpretable and flexible non-intrusive reduced-order models using reproducing kernel Hilbert spaces

This paper develops an interpretable, non-intrusive reduced-order modeling technique using regularized kernel interpolation. Existing non-intrusive approaches approximate the dynamics of a reduced-order model (ROM) by solving a data-driven least-squares regression problem for low-dimensional matrix operators. Our approach instead leverages regularized kernel interpolation, which yields an optimal approximation of the ROM dynamics from a user-defined reproducing kernel Hilbert space. We show that our kernel-based approach can produce interpretable ROMs whose structure mirrors full-order model structure by embedding judiciously chosen feature maps into the kernel. The approach is flexible and allows a combination of informed structure through feature maps and closure terms via more general nonlinear terms in the kernel. We also derive a computable a posteriori error bound that combines standard error estimates for intrusive projection-based ROMs and kernel interpolants. In conclusion, the approach is demonstrated in several numerical experiments that include comparisons to operator inference using both proper orthogonal decomposition and quadratic manifold dimension reduction.

Data-driven model reduction↗

DOE OSTI2021

An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars

The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Phenomenology of the Prethermal Many-Body Localized Regime

The dynamical phase diagram of interacting disordered systems has seen substantial revision over the past few years. Theory must now account for a large prethermal many-body localized regime in which thermalization is extremely slow, but not completely arrested. We derive a quantitative description of these dynamics in short-ranged one-dimensional systems using a model of successive many-body resonances. Here, the model explains the decay timescale of mean autocorrelators, the functional form of the decay—a stretched exponential—and relates the value of the stretch exponent to the broad distribution of resonance timescales. The Jacobi method of matrix diagonalization provides numerical access to this distribution, as well as a conceptual framework for our analysis. The resonance model correctly predicts the stretch exponents for several models in the literature. Successive resonances may also underlie slow thermalization in strongly disordered systems in higher dimensions, or with long-range interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Decentralized Low-Rank State Estimation for Power Distribution Systems

This article considers the low-observability state estimation problem in power distribution networks and develops a decentralized state estimation algorithm leveraging the matrix completion methodology. Matrix completion has been shown to be an effective technique in state estimation that exploits the low dimensionality of the power system measurements to recover missing information. This technique can utilize an approximate (linear) load flow model, or it can be used with no physical models in a network where no information about the topology or line admittance is available. The direct application of matrix completion algorithms requires solving a semi-definite programming (SDP) problem, which becomes computationally challenging for large networks. We therefore develop a decentralized algorithm that capitalizes on the popular proximal alternating direction method of multipliers (proximal ADMM). The method allows us to distribute the computation among different areas of the network, leading to a scalable algorithm. By doing all computations at individual control areas and only communicating with neighboring areas, the algorithm eliminates the need for data to be sent to a central processing unit and thus increases efficiency and contributes to the goal of autonomous control of distribution networks. We illustrate the advantages of the proposed algorithm numerically using standard IEEE test cases.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2024

Altered post-fracture systemic bone loss in a mouse model of osteocyte dysfunction

Femur fracture leads to loss of bone at uninjured skeletal sites, which may increase risk of subsequent fracture. Osteocytes, the most abundant bone cells, can directly resorb bone matrix and regulate osteoclast and osteoblast activity, but their role in systemic bone loss after fracture remains poorly understood. In this study we used a transgenic (TG+) mouse model that overexpresses human B-cell lymphoma 2 (BCL-2) in osteoblasts and osteocytes. This causes enhanced osteoblast proliferation, followed by disruption in lacunar-canalicular connectivity and massive osteocyte death by 10 wk of age. We hypothesized that reduced viable osteocyte density would decrease the magnitude of systemic bone loss after femur fracture, reduce perilacunar remodeling, and alter callus formation. Bone remodeling was assessed using serum biomarkers of bone formation and resorption at 5 d post-fracture. We used micro-computed tomography, high resolution x-ray microscopy, mechanical testing, and Raman spectroscopy to quantify the magnitude of systemic bone loss, as well as changes in osteocyte lacunar volume, bone strength, and bone composition 2 wk post-fracture. Fracture was associated with a reduction in circulating markers of bone resorption in non-transgenic (TG-) animals. TG+ mice exhibited high bone mass in the limbs, greater cortical elastic modulus and reduced post-yield displacement. After fracture, TG+ mice lost less trabecular bone than TG- mice, but conversely TG+ mice exhibited trends toward a lower yield point and reduced femoral cortical thickness after fracture, though these were not statistically significant. Lacunar density was greater in TG+ mice, but fracture did not alter lacunar volume in TG+ or TG- mice. These findings suggest that osteocytes potentially play a significant role in the post-traumatic systemic response to fracture, though the effects differ between trabecular and cortical bone.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Coexistence of non-Abelian chiral spin liquid and magnetic order in a spin-1 antiferromagnet

We study the ground-state properties of a spin-1 Heisenberg model on a square lattice with the first- and second-nearest-neighbor antiferromagnetic couplings $J_1$ and $J_2$ and a three-spin scalar chirality term $J_\chi$. Using the density matrix renormalization group calculation, we map out a global phase diagram including various magnetic order phases and an emergent quantum spin liquid phase. The nature of the spin liquid is identified as a bosonic non-Abelian Moore-Read state from the fingerprint of the entanglement spectra and identification of a full set of topological sectors. We further unveil a stripe magnetic order coexisting with this spin liquid. Furthermore, our results not only establish a rare example of non-Abelian spin liquids in simple spin systems but also demonstrate the coexistence of fractionalized excitations and magnetic order beyond mean-field descriptions.

2-dimensional systems↗

DOE OSTI2020

Source identification by non-negative matrix factorization combined with semi-supervised clustering

Machine-learning methods and apparatus are provided to solve blind source separation problems with an unknown number of sources and having a signal propagation model with features such as wave-like propagation, medium-dependent velocity, attenuation, diffusion, and/or advection, between sources and sensors. In exemplary embodiments, multiple trials of non-negative matrix factorization are performed for a fixed number of sources, with selection criteria applied to determine successful trials. A semi-supervised clustering procedure is applied to trial results, and the clustering results are evaluated for robustness using measures for reconstruction quality and cluster separation. The number of sources is determined by comparing these measures for different trial numbers of sources. Source locations and parameters of the signal propagation model can also be determined. Disclosed methods are applicable to a wide range of spatial problems including chemical dispersal, pressure transients, and electromagnetic signals, and also to non-spatial problems such as cancer mutation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Source identification by non-negative matrix factorization combined with semi-supervised clustering

Alexandrov, Boian S.↗

DOE OSTI2026

Simplex‐based model for nanoparticle grain identification in four‐dimensional scanning transmission electron microscopy data

Grain identification in polycrystalline nanoparticles, for example, determining which crystal phases are present at each spatial location, is fundamental to materials characterisation. This is particularly challenging when grains overlap extensively, as commonly occurs in four-dimensional scanning transmission electron microscopy (4D-STEM) datasets. We propose a simplex-based model (SBM) in which each simplex vertex represents the diffraction pattern (DP) of a pure grain, and the simplex edges and interior represent overlapping grains. Our SBM grain identification algorithm operates on the Bragg disk (BD) data matrix distilled from the 4D-STEM data to identify the grain membership at each scan position, together with a BD feature matrix whose columns represent the DPs for each constituent grain, which is important for identifying the crystal structure of each grain. We solve the model using a two-stage algorithm. In Stage 1, we adapt a linear mixing algorithm to estimate an initial BD feature matrix whose columns represent DPs of potentially overlapping grains. Our Stage 2 algorithm incorporates sparsity considerations to transform the initial BD feature matrix so that its columns represent DPs of pure grains. Using simulated datasets with various grain configurations, we demonstrate that SBM recovers both the BD feature matrix and membership maps more accurately than existing methods, even when a grain lacks any pure region and completely overlaps with other grains.

4D-STEM segmentation↗