Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tensor factorizations”

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 73 records · Page 4

Mechanical form factors and densities of nonrelativistic fermions

The hadron physics community has been actively debating the interpretation of so-called mechanical properties of hadrons. Nonrelativistic quantum-mechanical systems like the hydrogen atom have been appealed to in these debates as analogies. Since such appeals are likely to continue, it is important to have Galilei-covariant expressions for matrix elements of the energy-momentum tensor. In this work, I obtain Galilei-covariant breakdowns of such matrix elements into mechanical form factors, with a special focus on spin-half states. I additionally study the spatial densities associated with these form factors, using the pilot wave interpretation to guide their breakdown into contributions from internal structure and from quantum-mechanical effects such as wave packet dispersion. For completeness, I also obtain nonrelativistic Breit frame densities.

form factors↗

Data-Driven Insights into the Structural Essence of Plasticity in High-Entropy Alloys

The heterogeneous mechanical response of a crystalline alloy with multiple principal elements was investigated using molecular dynamics simulations. The local configuration of the alloy in its quiescent state was characterized by the variables derived from the gyration tensor and the atomic electronegativity. A multivariate analysis identified the geometric and chemical factors that influenced the atomic packing variations. Further, upon straining, the non-affine displacement exhibited spatial heterogeneity. A statistical correlation was established between the local yield events and the specific features of the local configuration. Our findings, validated by the performance metrics analysis, provided a structural criterion for the instability mechanisms in high-entropy alloys (HEAs) and enhanced the understanding of their plasticity.

36 MATERIALS SCIENCE↗

Near-Efficient and Non-Asymptotic Multiway Inference

We establish non-asymptotic efficiency guarantees for tensor decomposition–based inference in count data models. Under a Poisson framework, we consider two related goals: (i) parametric inference , the estimation of the full distributional parameter tensor, and (ii) multiway analysis , the recovery of its canonical polyadic (CP) decomposition factors. Our main result shows that in the rank-one setting, a rank-constrained maximum-likelihood estimator achieves multiway analysis with variance matching the Cramér–Rao Lower Bound (CRLB) up to absolute constants and logarithmic factors. This provides a general framework for studying “near-efficient” multiway estimators in finite-sample settings. For higher ranks, we illustrate that our multiway estimator may not attain the CRLB; nevertheless, CP-based parametric inference remains nearly minimax optimal, with error bounds that improve on prior work by offering more favorable dependence on the CP rank. Numerical experiments corroborate near-efficiency in the rank-one case and highlight the efficiency gap in higher-rank scenarios.

97 MATHEMATICS AND COMPUTING↗

Improved Tensor Current Limit from 8 B 𝛽 Decay Including New Recoil-Order Calculations

A precision measurement of the 𝛽 + decay of 8 B was performed using the Beta-decay Paul Trap to determine the 𝛽−𝜈 angular correlation coefficient 𝑎 𝛽⁢𝜈 . The experimental results were combined with new ab initio symmetry-adapted no-core shell-model calculations to yield the second-most precise measurement from Gamow-Teller decays, 𝑎 𝛽⁢𝜈 = −0.3345 ± 0.001⁢9 stat ± 0.002⁢1 syst . This value agrees with the standard model value of −1/3 and improves uncertainties in 8 B by nearly a factor of 2. By combining results from 8 B and 8 Li , a tight limit on tensor current coupling to right-handed neutrinos was obtained. A recent global evaluation of all other precision 𝛽 decay studies suggested a nonzero value for right-handed neutrino coupling in contradiction with the standard model at just above 3⁢𝜎. Finally, the present results are of comparable sensitivity and do not support this finding.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hybrid simulations of collisionless ion tearing

Ion tearing in a thin collisionless neutral sheet is investigated using hybrid (particle ion, massless fluid electron) simulations in which dissipation effects are included through a model of the full electron pressure tensor in which self-consistently generated anisotropies are reduced via an isotropization time factor. The linear growth rates and ion dynamics are in qualitative agreement with linear theory and one-component simulations. The reconnection rate increases with the isotropization time.

Hesse, Michael↗

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗

A finite element tensor approach to plate buckling and postbuckling.

A practical finite element method is presented for geometrically nonlinear (von Karman) plate problems. Symmetric strain energy tensors provide an efficient formulation and solution, and allow the effects of initial imperfections to be considered. An unbalanced force iteration is employed, with the advantage that the time consuming Gauss decomposition needs to be performed only once. The convergence is controlled by an automatically computed relaxation factor, and this makes the method applicable in advanced stages of non-linearity. Several different triangular elements are used in explicitly deriving the required tensors.

Vos, R. G.↗

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗

Determination of Stress Coefficient Terms in Cracked Solids for Monoclinic Materials with Plane Symmetry at x3 = 0

Determination of all the coefficients in the crack tip field expansion for monoclinic materials under two-dimensional deformation is presented in this report. For monoclinic materials with a plane of material symmetry at x(sub 3) = 0, the in-plane deformation is decoupled from the anti-plane deformation. In the case of in-plane deformation, utilizing conservation laws of elasticity and Betti's reciprocal theorem, together with selected auxiliary fields, T-stress and third-order stress coefficients near the crack tip are evaluated first from path-independent line integrals. To determine the T-stress terms using the J-integral and Betti's reciprocal work theorem, auxiliary fields under a concentrated force and moment acting at the crack tip are used respectively. Through the use of Stroh formalism in anisotropic elasticity, analytical expressions for all the coefficients including the stress intensity factors are derived in a compact form that has surprisingly simple structure in terms of the Barnett-Lothe tensors, L. The solution forms for degenerated materials, orthotropic, and isotropic materials are presented.

Yuan, F. G.↗

Mass radius and D-term of atomic nuclei in relativistic mean field theory

Based on relativistic mean field theory for atomic nuclei, we compute the mass radius and other radii associated with the energy momentum tensor for dozens of spin-0 nuclei across the nuclear chart. We also compute the D-term of these nuclei, the forward limit of the gravitational form factor 𝐷⁡(𝑡=0)=𝐷. The dependence on the neutron number 𝑁 is systematically studied for calcium (Ca), nickel (Ni), zirconium (Zr), tin (Sn), and lead (Pb) isotopes. Remarkably, |𝐷| does not monotonically increase with 𝑁. Instead, it exhibits local maxima and minima when 𝑁 equals a magic number and even a submagic number. This results in characteristic kinks in the mass, scalar, tensor, and shear radii of these isotopes. Our work for the first time elucidates the strong sensitivity of the various mechanical properties of nuclei to the nuclear shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GPU-Accelerated Analytic Simulation of Sparse Ionization Signal Formation in Pixelated Projection Detector

This paper presents a GPU-accelerated simulation package, TRED, for next-generation neutrino detectors with pixelated charge readout, leveraging community-driven software ecosystems to ensure adaptability and extensibility. We introduce two generic contributions: (i) an effective-charge representation based on Gaussian quadrature rules, in which the linear- interpolation factors for the field response inside each voxel are absorbed into the effective charge, and (ii) a sparse, block- binned tensor representation that enables efficient FFT-based computation of induced signals on readout electrodes for sparsely activated detector volumes. The former captures structure inside a voxel without dense sampling, while the latter achieves low memory usage and scalable runtime, as demonstrated in bench- mark studies. The underlying data representation is applicable to large-scale detectors and to other computational problems involving sparse activity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ring Current-Electromagnetic Ion Cyclotron Waves Coupling

The effect of Electromagnetic Ion Cyclotron (EMIC) waves, generated by ion temperature anisotropy in Earth s ring current (RC), is the best known example of wave- particle interaction in the magnetosphere. Also, there is much controversy over the importance of EMIC waves on RC depletion. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt (RB) by EMIC wave scattering during a magnetic storm. That is why the calculation of EMIC waves must be a very critical part of the space weather studies. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch-angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) the EMIC waves saturation due to their modulation instability and LHW generation are included as an additional factor that contributes to this process; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗

Exploring Quantum State Preparation Using Tensor Networks and Sparse Wavefunction Simulations

The variational quantum eigenvalue solver is a powerful hybrid quantum-classical approach that has been suggested as a candidate method to run on near-term quantum hardware for computing ground state electronic energies of molecular systems. However, even for small molecules, the number of variational parameters and qubits required to minimize the electronic energy is beyond the reach of current quantum computers except for small basis sets. We explore a new paradigm for state preparation where we test how much of the optimization can be approximately prepared with classical computers to reduce the number of optimization steps performed using a quantum device. By adapting a recent algorithm for the factorized form of the UCC ansatz, we can study molecular electronic structure problems with up to 64 qubits. In addition, we also test a related approach of using tensor networks to optimize quantum circuits in order to benchmark various lattice models. We present results using these approaches and discuss strategies for incorporating these ideas into variational algorithms involving near-term quantum computers. Our results help demonstrate the strength of the UCC ansatz and address pressing questions about optimal initial parameterizations and circuit construction.

quantum computing↗

Variational principle for a prototype Rastall theory of gravitation

A prototype of Rastall's theory of gravity, in which the divergence of the energy-momentum tensor is proportional to the gradient of the scalar curvature, is shown to be derivable from a variational principle. Both the proportionality factor and the unrenormalized gravitational constant are found to be covariantly constant, but not necessarily constant. The prototype theory is, therefore, a gravitational theory with variable gravitational constant.

Smalley, L. L.↗

Complex deuteron NMR signals

To determine the spin polarization of deuterons, nuclear magnetic resonance (NMR) is used. This is necessary for polarized targets, such as for the upcoming A zz and b 1 experiment at Jefferson Lab. NMR measures the impedance of a solenoid around a deuterated sample. Although the impedance is a complex value, conventionally only the real part of the impedance has been used for this purpose. However, often the tune is not precisely real, meaning the signal has at least some small imaginary portion. This conventionally has been dealt with by an offset parameter, such as Dulya’s false asymmetry method. For vector polarization, this suffices, as the tuning is factored into the overall error of the results, and for a small phase angle doesn’t make much of a difference. However, for tensor polarization, the exact lineshape of the signal is quite significant, and treating the impedance as complex during analysis removes the need for an offset parameter. As a result, it also provides more accurate results, as the conventional false asymmetry method over- or underestimates polarization, depending on the sign of the phase angle.

McClellan, Michael [University of New Hampshire, D↗

A new determination of the halo luminosity density of the galaxy

The distribution of tangential velocities for arbitrary halo velocity dispersions and rotation rates are explicitly computed, assuming a particular form of the halo velocity distribution function near the sun. It is found that for a specific model of the halo velocity dispersion tensor that Schmidt's (1975) estimate of the mass carried by halo stars of 0.25-0.70 solar masses should be revised by a factor of between 1.0 and 2.0, depending on the rotation rate of the halo. Although obtaining a mass density in all stars from the stellar density in that mass range is extremely uncertain, it is found, in agreement with earlier work, that even under rather extreme assumptions about the remnants the observed stars and expected remnants fail to provide the theoretically expected mass. Comparison of the luminosity to the theoretical mass yields a mass to light ratio consistent with that of rich clusters of galaxies.

Richstone, D. O.↗

Ring Current Dynamic in the Presence of EMIC Waves

The effect of EMIC waves, generated by a positive ion temperature anisotropy on Earth s RC ions is one of the best known examples of wave-particle interaction in the magnetosphere and the most controversial mechanism of RC losses. Under certain conditions, relativistic electrons with energy 21 MeV can be removed from the outer RB by EMIC wave scattering during a magnetic storm much faster than by any other loss mechanisms. That is why the calculation of EMIC waves is a very critical part of the NASA LWS program. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch- angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) LHW are included as an additional factor that contributes to saturation process of EMIC waves; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗

Impossibility of obtaining time-independent, three-dimensional, spherically symmetric densities of confined systems of relativistically moving constituents

The quantum-mechanical definition of probability, the uncertainty principle, and Poincaré invariance provide strong basic restrictions on the ability to define spatial densities associated with form factors describing the properties of confined systems of relativistically moving constituents. Despite this, many papers ignore one or more of these restrictions. Here I show how to obtain time-independent , two-dimensional densities that are consistent with the stated restrictions. This is done using the light-front, infinite momentum frame formalism. Two-dimensional density interpretations of the axial-vector form factor and all three gravitational form factors is obtained. The resulting mass radius is smaller than the charge radius. Additionally, an expression of a two-dimensional mass density related to the trace of the energy momentum tensor is obtained. I also show that all known methods for finding three-dimensional densities—using the Breit frame, Abel transformations, Wigner distributions, and spherically symmetric wave packets with vanishing spatial extent—violate the basic restrictions in different ways. Furthermore, the use of the latter leads to densities that vanish almost everywhere in space as time increases from an initial value.

form factors↗