Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Theory”

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 1,117 records · Page 62

Data from a four-day long microcosm experiment addressing the destabilization of artificial mineral-associated organic matter by model root exudates embedded in a soil matrix from the Rocky Mountain Biological Laboratory (Gothic, CO, USA), 2019

This dataset provides data collected during a four-day long laboratory soil microcosm experiment testing the efficacy of root exudate-driven mineral-associated organic matter destabilization. This dataset contains four data files in comma-separate values (*.csv). The files provide the metadata and the experimental results on microbial respiration, MAOM-derived respiration, and sequential mineral-extractions. This data was used to produce the figures in Bölscher et al., 2026. The results of the experiment can be found in the open access article Bölscher et al., 2026 (https://doi.org/10.1016/j.soilbio.2026.110276). Abstract: Mineral-associated organic matter (MAOM) is often considered stable, but root exudates can destabilize MAOM via various pathways. Theory and model system studies suggest that direct MAOM destabilization by strong ligands, like oxalic acid, or reducing agents, like catechol, is more effective than indirect, microbial-mediated MAOM destabilization, stimulated by less reactive compounds like glucose. Here, we demonstrate that the presence of a soil matrix alters the efficacy of exudate-driven MAOM destabilization pathways. Glucose and catechol destabilized significantly greater amounts of MAOM from ferrihydrite and aluminum hydroxide (Al (OH)3) embedded in a soil matrix than oxalic acid. Our findings indicate that indirect, microbial-mediated MAOM destabilization may play a larger role than direct MAOM destabilization in soil environments.

Destabilization↗

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalous dimensions from gauge couplings in SMEFT with right-handed neutrinos

Standard Model Neutrino Effective Field Theory (SMNEFT) is an effective theory with Standard Model (SM) gauge-invariant operators constructed only from SM and right-handed neutrino fields. For the full set of dimension-six SMNEFT operators, we present the gauge coupling terms of the one-loop anomalous dimension matrix for renormalization group evolution (RGE) of the Wilson coefficients between a new physics scale and the electroweak scale. We find that the SMNEFT operators can be divided into five subsets which are closed under RGE. Our results apply for both Dirac and Majorana neutrinos. We also discuss the operator mixing pattern numerically and comment on some interesting phenomenological implications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

Subleading corrections to the S 3 free energy of necklace quiver theories dual to massive IIA

We investigate the S 3 free energy of N = 3 Chern-Simons-matter quiver gauge theories with gauge group U(N) r (r ≥ 2) where the sum of Chern-Simons levels does not vanish, beyond the leading order in the large-N limit. We take two different approaches to explore the sub-leading structures of the free energy. First we evaluate the matrix integral for the partition function in the ’t Hooft limit using a saddle point approximation. Second we use an ideal Fermi-gas model to compute the same partition function, but in the limit of fixed Chern-Simons levels. The resulting expressions for the free energy F = – log Z are then compared in the overlapping parameter regime. The Fermi-gas approach also hints at a universal 1/6 log N correction to the free energy. Since the quiver gauge theories we consider are dual to massive Type IIA theory, we expect the sub-leading correction of the planar free energy in the large ’t Hooft parameter limit to match higher-derivative corrections to the tree-level holographic dual free energy, which have not yet been fully investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three-dimensional imaging of pion using lattice QCD: generalized parton distributions

In this work, we report a lattice calculation of x-dependent valence pion generalized parton distributions (GPDs) at zero skewness with multiple values of the momentum transfer −t. The calculations are based on an N f = 2 + 1 gauge ensemble of highly improved staggered quarks with Wilson-Clover valence fermion. The lattice spacing is 0.04 fm, and the pion valence mass is tuned to be 300 MeV. We determine the Lorentz-invariant amplitudes of the quasi-GPD matrix elements for both symmetric and asymmetric momenta transfers with similar values and show the equivalence of both frames. Then, focusing on the asymmetric frame, we utilize a hybrid scheme to renormalize the quasi-GPD matrix elements obtained from the lattice calculations. After the Fourier transforms, the quasi-GPDs are then matched to the light-cone GPDs within the framework of large momentum effective theory with improved matching, including the next-to-next-to-leading order perturbative corrections, and leading renormalon and renormalization group resummations. We also present the 3-dimensional image of the pion in impact-parameter space through the Fourier transform of the momentum transfer −t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Single-soft emissions for amplitudes with two colored particles at three loops

We compute the three-loop correction to the universal single-soft emission current for the case of scattering amplitudes with two additional color-charged partons. We present results valid for QCD and $\mathcal{N}$ = 4 super-symmetric Yang-Mills theory. To achieve our results we develop a new integrand expansion technique for scattering amplitudes in the presence of soft emissions. Furthermore, we obtain contributions from single final-state parton matrix elements to the Higgs boson and Drell-Yan production cross section at next-to-next-to-next-to-next-to leading order (N 4 LO) in perturbative QCD in the threshold limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A microstructure-based modeling approach to predict the mechanical properties of Zr alloy with hydride precipitates

In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.

Kulkarni, Shank S.↗

Correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field methods

The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Directional Laplacian Centrality for Cyber Situational Awareness

Cyber operations is drowning in diverse, high-volume, multi-source data. To get a full picture of current operations and identify malicious events and actors, analysts must see through data generated by a mix of human activity and benign automated processes. Although many monitoring and alert systems exist, they typically use signature-based detection methods. We introduce a general method rooted in spectral graph theory to discover patterns and anomalies without a priori knowledge of signatures. We derive and propose a new graph-theoretic centrality measure based on the derivative of the graph Laplacian matrix in the direction of a vertex. To build intuition about our measure, we show how it identifies the most central vertices in standard network datasets and compare to other graph centrality measures. Finally, we focus our attention on studying its effectiveness in identifying important IP addresses in network flow data. Using both real and synthetic network flow data, we conduct several experiments to test our measure’s sensitivity to two types of injected attack profiles and show that vertices participating in injected attack profiles exhibit noticeable changes in our centrality measures, even when the injected anomalies are relatively small, and in the presence of simulated network dynamics.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

nuclear-score-maximization v1.0

This software library presents efficient and multithreaded implementations of matrix low rank approximation via column selection in C++17 code. The algorithms are described in Fornace, Mark, and Michael Lindsey. "Column and row subset selection using nuclear scores: algorithms and theory for Nystro m approximation, CUR decomposition, and graph Laplacian reduction." arXiv preprint arXiv:2407.01698 (2024). The presented methods are by-and-large ver novel, have provable approximation guarantees, multiple use-cases, and exhibit higher quality approximations on a variety of studied examples.

Fornace, Mark↗

Viscosity of rock-ice mixtures and applications to the evolution of icy satellites

Theory and experiments are used to establish lower and upper bounds on the ratio of actual viscosity to pure ice viscosity for a suspension of rock particles in a water ice matrix. A rheological model for rock-ice mixtures is described, establishing bounds for the range of possible viscosity enhancements provided by a suspension of silicate spheres in an ice matrix. A parametrized thermal convection model is described and used to determine a criterion for criticality, defined as the heat flow and/or silicate volume fraction for which the satellite temperature profile intercepts the melting curve of water ice. The consequences of achieving this critical state are examined, and it is shown that under certain circumstances a 'runaway' differentiation can occur in which the silicates settle to form a core and extensive melting of water ice takes place, the latent heat being supplied by the gravitational energy of differentiation. A possible application of these results to Ganymede and Callisto is described.

Friedson, A. J.↗

Radiative transfer theory for active remote sensing of a layer of nonspherical particles

The radiative transfer theory is applied to calculate the scattering by a layer of randomly positioned and oriented nonspherical particles. The scattering amplitude functions of each individual particle are calculated with Waterman's T matrix method, which utilizes vector spherical wave functions for expansion of incident, scattered, and surface fields. The orientation of the particles is described by a probability density function of the Eulerian angles of rotation. A rotation matrix is used to relate the T matrix of the principal frame to that of the natural frame of the particle. The extinction matrix and phase matrix of the radiative transfer equations are expressed in terms of the T matrix elements. The extinction matrix for nonspherical particles is generally nondiagonal. There are only two attenuation rates in a specified direction of propagation. The radiative transfer equations are solved by an iterative method to first order in albedo. Numerical results are illustrated as functions of incidence angle and frequency with applications to active remote sensing.

Tsang, L.↗

Extended aeroelastic analysis for helicopter rotors with prescribed hub motion and blade appended penduluum vibration absorbers

The mathematical development for the expanded capabilities of the G400 rotor aeroelastic analysis was examined. The G400PA expanded analysis simulates the dynamics of all conventional rotors, blade pendulum vibration absorbers, and the higher harmonic excitations resulting from prescribed vibratory hub motions and higher harmonic blade pitch control. The methodology for modeling the unsteady stalled airloads of two dimensional airfoils is discussed. Formulations for calculating the rotor impedance matrix appropriate to the higher harmonic blade excitations are outlined. This impedance matrix, and the associated vibratory hub loads, are the rotor dynamic characteristic elements for use in the simplified coupled rotor/fuselage vibration analysis (SIMVIB). Updates to the development of the original G400 theory, program documentation, user instructions and information are presented.

Bielawa, R. L.↗

The Ohio State 1991 geopotential and sea surface topography harmonic coefficient models

The computation is described of a geopotential model to deg 360, a sea surface topography model to deg 10/15, and adjusted Geosat orbits for the first year of the exact repeat mission (ERM). This study started from the GEM-T2 potential coefficient model and it's error covariance matrix and Geosat orbits (for 22 ERMs) computed by Haines et al. using the GEM-T2 model. The first step followed the general procedures which use a radial orbit error theory originally developed by English. The Geosat data was processed to find corrections to the a priori geopotential model, corrections to a radial orbit error model for 76 Geosat arcs, and coefficients of a harmonic representation of the sea surface topography. The second stage of the analysis took place by doing a combination of the GEM-T2 coefficients with 30 deg gravity data derived from surface gravity data and anomalies obtained from altimeter data. The analysis has shown how a high degree spherical harmonic model can be determined combining the best aspects of two different analysis techniques. The error analysis was described that has led to the accuracy estimates for all the coefficients to deg 360. Significant work is needed to improve the modeling effort.

Rapp, Richard H.↗

Application of vector-valued rational approximations to the matrix eigenvalue problem and connections with Krylov subspace methods

Let F(z) be a vectored-valued function F: C approaches C sup N, which is analytic at z=0 and meromorphic in a neighborhood of z=0, and let its Maclaurin series be given. We use vector-valued rational approximation procedures for F(z) that are based on its Maclaurin series in conjunction with power iterations to develop bona fide generalizations of the power method for an arbitrary N X N matrix that may be diagonalizable or not. These generalizations can be used to obtain simultaneously several of the largest distinct eigenvalues and the corresponding invariant subspaces, and present a detailed convergence theory for them. In addition, it is shown that the generalized power methods of this work are equivalent to some Krylov subspace methods, among them the methods of Arnoldi and Lanczos. Thus, the theory provides a set of completely new results and constructions for these Krylov subspace methods. This theory suggests at the same time a new mode of usage for these Krylov subspace methods that were observed to possess computational advantages over their common mode of usage.

Sidi, Avram↗

A User''s Guide to the Zwikker-Kosten Transmission Line Code (ZKTL)

This user's guide documents updates to the Zwikker-Kosten Transmission Line Code (ZKTL). This code was developed for analyzing new liner concepts developed to provide increased sound absorption. Contiguous arrays of multi-degree-of-freedom (MDOF) liner elements serve as the model for these liner configurations, and Zwikker and Kosten's theory of sound propagation in channels is used to predict the surface impedance. Transmission matrices for the various liner elements incorporate both analytical and semi-empirical methods. This allows standard matrix techniques to be employed in the code to systematically calculate the composite impedance due to the individual liner elements. The ZKTL code consists of four independent subroutines: 1. Single channel impedance calculation - linear version (SCIC) 2. Single channel impedance calculation - nonlinear version (SCICNL) 3. Multi-channel, multi-segment, multi-layer impedance calculation - linear version (MCMSML) 4. Multi-channel, multi-segment, multi-layer impedance calculation - nonlinear version (MCMSMLNL) Detailed examples, comments, and explanations for each liner impedance computation module are included. Also contained in the guide are depictions of the interactive execution, input files and output files.

Kelly, J. J.↗

Multiscale Failure Analysis of a 3D Woven Unit Cell Containing Defects

Multiscale failure simulations have been performed for a Three-dimensional woven composite unit cell considering five, or more, length scales spanning the woven composite mesoscale to the sub-microscale voids. The multiscale recursive micromechanics approach, which enables recursive integration of general micromechanics theories over an arbitrary number of length scales, has been employed within the NASA Multiscale Analysis Tool. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories, and considers failure in the constituent materials using a simple damage model. Baseline results, containing distributed voids, are compared to uniaxial experimental data for an AS4 carbon fiber/ RTM6 epoxy matrix 3D orthogonal woven composite with good agreement in terms of global stiffness and global failure stress. The simulations demonstrate that the 3D woven composite exhibits damage tolerance through sustaining increasing axial load far beyond the first initiation of damage. The multiscale model is used to examine the nonlinear response of the material to other loading conditions. Case studies, motivated by X-ray computed tomography data, are presented on the effects of manufacturing induced voids and cracks.

3D woven↗