Search NASA⌕ Search

SEARCH · Search NASA

Results for “tensor times same vector”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Quantum corrections for spinning particles in de Sitter

We compute the one-loop quantum corrections to the gravitational potentials of a spinning point particle in a de Sitter background, due to the vacuum polarisation induced by conformal fields in an effective field theory approach. We consider arbitrary conformal field theories, assuming only that the theory contains a large number N of fields in order to separate their contribution from the one induced by virtual gravitons. The corrections are described in a gauge-invariant way, classifying the induced metric perturbations around the de Sitter background according to their behaviour under transformations on equal-time hypersurfaces. There are six gauge-invariant modes: two scalar Bardeen potentials, one transverse vector and one transverse traceless tensor, of which one scalar and the vector couple to the spinning particle. The quantum corrections consist of three different parts: a generalisation of the flat-space correction, which is only significant at distances of the order of the Planck length; a constant correction depending on the undetermined parameters of the renormalised effective action; and a term which grows logarithmically with the distance from the particle. This last term is the most interesting, and when resummed gives a modified power law, enhancing the gravitational force at large distances. As a check on the accuracy of our calculation, we recover the linearised Kerr-de Sitter metric in the classical limit and the flat-space quantum correction in the limit of vanishing Hubble constant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Legendre submanifolds in contact manifolds as attractors and geometric nonequilibrium thermodynamics

It has been proposed that equilibrium thermodynamics is described on Legendre submanifolds in contact geometry. It is shown in this paper that Legendre submanifolds embedded in a contact manifold can be expressed as attractors in phase space for a certain class of contact Hamiltonian vector fields. By giving a physical interpretation that points outside the Legendre submanifold can represent nonequilibrium states of thermodynamic variables, in addition to that points of a given Legendre submanifold can represent equilibrium states of the variables, this class of contact Hamiltonian vector fields is physically interpreted as a class of relaxation processes, in which thermodynamic variables achieve an equilibrium state from a nonequilibrium state through a time evolution, a typical nonequilibrium phenomenon. Geometric properties of such vector fields on contact manifolds are characterized after introducing a metric tensor field on a contact manifold. It is also shown that a contact manifold and a strictly convex function induce a lower dimensional dually flat space used in information geometry where a geometrization of equilibrium statistical mechanics is constructed. Legendre duality on contact manifolds is explicitly stated throughout.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Phenomenology in minimal theory of massive gravity

We investigate the minimal theory of massive gravity (MTMG) recently introduced. After reviewing the original construction based on its Hamiltonian in the vielbein formalism, we reformulate it in terms of its Lagrangian in both the vielbein and the metric formalisms. It then becomes obvious that, unlike previous attempts in the literature of Lorentz-violating massive gravity, not only the potential but also the kinetic structure of the action is modified from the de Rham-Gabadadze-Tolley (dRGT) massive gravity theory. We confirm that the number of physical degrees of freedom in MTMG is two at fully nonlinear level. This proves the absence of various possible pathologies such as superluminality, acausality and strong coupling. Afterwards, we discuss the phenomenology of MTMG in the presence of a dust fluid. We find that on a flat homogeneous and isotropic background we have two branches. One of them (self-accelerating branch) naturally leads to acceleration without the genuine cosmological constant or dark energy. For this branch both the scalar and the vector modes behave exactly as in general relativity (GR). The phenomenology of this branch differs from GR in the tensor modes sector, as the tensor modes acquire a non-zero mass. Hence, MTMG serves as a stable nonlinear completion of the self-accelerating cosmological solution found originally in dRGT theory. The other branch (normal branch) has a dynamics which depends on the time-dependent fiducial metric. For the normal branch, the scalar mode sector, even though as in GR only one scalar mode is present (due to the dust fluid), differs from the one in GR, and, in general, structure formation will follow a different phenomenology. The tensor modes will be massive, whereas the vector modes, for both branches, will have the same phenomenology as in GR.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Viscous corrections to electromagnetic emissivities in QCD

We provide a general framework for the derivation of the hydrodynamical corrections to the QCD electromagnetic emissivities in a viscous fluid. Here, assuming that the emission times are short in comparison to the fluid evolution time, we show that the leading corrections in the fluid gradients are controlled by the bulk and shear tensors times pertinent response functions involving the energy-momentum tensor. In a hadronic fluid phase, we detail these contributions using spectral functions. Using the vector dominance approximation, we show that the bulk viscosity correction to the photon rate is sizable, while the shear viscosity is negligible for about all frequencies. In the partonic phase near the transition temperature, we provide an assessment of the viscous corrections to the photon and dilepton emissions, using a nonperturbative quark-gluon plasma with soft thermal gluonic corrections in the form of operators of leading mass dimension. Again, the thermal bulk viscosity corrections are found to be larger than the thermal shear viscosity corrections at all energies for both the photon and dilepton in the partonic phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗
DOE OSTI2019

Anomaly detection in scientific data using joint statistical moments

We propose an anomaly detection method for multi-variate scientific data based on analysis of high-order joint moments. Using kurtosis as a reliable measure of outliers, we suggest that principal kurtosis vectors, by analogy to principal component analysis (PCA) vectors, signify the principal directions along which outliers appear. The inception of an anomaly, then, manifests as a change in the principal values and vectors of kurtosis. Obtaining the principal kurtosis vectors requires decomposing a fourth order joint cumulant tensor for which we use a simple, computationally less expensive approach that involves performing a singular value decomposition (SVD) over the matricized tensor. We demonstrate the efficacy of this approach on synthetic data, and develop an algorithm to identify the occurrence of a spatial and/or temporal anomalous event in scientific phenomena. The algorithm decomposes the data into several spatial sub-domains and time steps to identify regions with such events. Feature moment metrics, based on the alignments of the principal kurtosis vectors, are computed at each sub-domain and time step for all features to quantify their relative importance towards the overall kurtosis in the data. Accordingly, spatial and temporal anomaly metrics for each sub-domain are proposed using the Hellinger distance of the feature moment metric distribution from a suitable nominal distribution. Finally, we apply the algorithm to two turbulent auto-ignition combustion cases and demonstrate that the anomaly metrics reliably capture the occurrence of auto-ignition in relevant spatial sub-domains at the right time steps.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Time-reversal symmetry, anomalies, and dualities in (2+1)$d$

We study continuum quantum field theories in 2+1 dimensions with time-reversal symmetry \cal T . The standard relation {\cal T}^2=(-1)^F is satisfied on all the “perturbative operators” i.e. polynomials in the fundamental fields and their derivatives. However, we find that it is often the case that acting on more complicated operators {\cal T}^2=(-1)^F {\cal M} with \cal M a non-trivial global symmetry. For example, acting on monopole operators, \cal M could be \pm1 ± 1 depending on the magnetic charge. We study in detail U(1) U ( 1 ) gauge theories with fermions of various charges. Such a modification of the time-reversal algebra happens when the number of odd charge fermions is 2 ~{\rm mod }~4 , e.g. in QED with two fermions. Our work also clarifies the dynamics of QED with fermions of higher charges. In particular, we argue that the long-distance behavior of QED with a single fermion of charge 2 2 is a free theory consisting of a Dirac fermion and a decoupled topological quantum field theory. The extension to an arbitrary even charge is straightforward. The generalization of these abelian theories to SO(N) S O ( N ) gauge theories with fermions in the vector or in two-index tensor representations leads to new results and new consistency conditions on previously suggested scenarios for the dynamics of these theories. Among these new results is a surprising non-abelian symmetry involving time-reversal.

Cordova, Clay↗
DOE OSTI2023

Phonon Hall Viscosity of Ionic Crystals

When time-reversal symmetry is broken, the low-energy description of acoustic lattice dynamics allows for a dissipationless component of the viscosity tensor, the phonon Hall viscosity, which captures how phonon chirality grows with the wave vector. In this work, we show that, in ionic crystals, a phonon Hall viscosity contribution is produced by the Lorentz forces on moving ions. We calculate typical values of the Lorentz force contribution to the Hall viscosity using a simple square lattice toy model, and we compare it with literature estimates of the strengths of other Hall-viscosity mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2025

Observation of an Axial-Vector State in the Study of the Decay ψ ( 3686 ) → ϕ η η ′

Using ( 2712.4 ± 14.3 ) × 10 6 ψ ( 3686 ) events collected with the BESIII detector at BEPCII, a partial wave analysis of the decay ψ ( 3686 ) → ϕ η η ′ is performed with the covariant tensor approach. In addition to the established states h 1 ( 1900 ) and ϕ ( 2170 ) , an axial-vector state with a mass near 2.3 GeV / c 2 is observed for the first time. Its mass and width are measured to be 2316 ± 9 stat ± 3 0 syst MeV / c 2 and 89 ± 1 5 stat ± 2 6 syst MeV , respectively. The product branching fractions of B [ ψ ( 3686 ) → X ( 2300 ) η ′ ] B [ X ( 2300 ) → ϕ η ] and B [ ψ ( 3686 ) → X ( 2300 ) η ] B [ X ( 2300 ) → ϕ η ′ ] are determined to be ( 4.8 ± 1.3 stat ± 0.7 syst ) × 10 − 6 and ( 2.2 ± 0.7 stat ± 0.7 syst ) × 10 − 6 , respectively. The branching fraction B [ ψ ( 3686 ) → ϕ η η ′ ] is measured for the first time to be ( 3.14 ± 0.1 7 stat ± 0.2 4 syst ) × 10 − 5 . The first uncertainties are statistical and the second are systematic. Published by the American Physical Society 2025

Ablikim, M.↗
DOE OSTI2017

Spin-density fluctuations and the fluctuation-dissipation theorem in 3 d ferromagnetic metals

Spatial and time scales of spin-density fluctuations (SDFs) were analyzed in 3d ferromagnets using ab initio linear-response calculations of complete wave-vector and energy dependence of the dynamic spin susceptibility tensor. We demonstrate that SDFs are spread continuously over the entire Brillouin zone and while the majority of them reside within the 3d bandwidth, a significant amount comes from much higher energies. A validity of the adiabatic approximation in spin dynamics is discussed. The SDF spectrum is shown to have two main constituents: a minor low-energy spin-wave contribution and a much larger high-energy component from more localized excitations. Furthermore, using the fluctuation-dissipation theorem, the on-site spin correlator and the related effective fluctuating moment were properly evaluated and their universal dependence on the 3d band population is further discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Comment on “The Mössbauer rotor experiment and the general theory of relativity” [Ann. Physics 368 (2016) 258–266]

Highlights: •Our approach is considered from the viewpoint of arbitrary moving continuum. •Continuum characteristics are determined by a deformation tensor, a stress tensor, a strain velocity tensor, a spin tensor (Sing (1963)) and the first curvature vectors of the world lines of medium particles. •All these physical values are located on the hypersurfaces orthogonal to the world lines of medium particles and specify the “physical” space. •The time counting from the “physical” space permits to obtain the red shift of the frequency at the end of the tube that coincides with Corda result (Corda (2016)).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Transition matrices and orbitals from reduced density matrix theory

In this contribution, we report two different methodologies for characterizing the electronic structure reorganization occurring when a chromophore undergoes an electronic transition. For the first method, we start by setting the theoretical background necessary to the reinterpretation through simple tensor analysis of (i) the transition density matrix and (ii) the natural transition orbitals in the scope of reduced density matrix theory. This novel interpretation is made more clear thanks to a short compendium of the one-particle reduced density matrix theory in a Fock space. The formalism is further applied to two different classes of excited states calculation methods, both requiring a single-determinant reference, that express an excited state as a hole-particle mono-excited configurations expansion, to which particle-hole correlation is coupled (time-dependent Hartree-Fock/time-dependent density functional theory) or not (configuration interaction single/Tamm-Dancoff approximation). For the second methodology presented in this paper, we introduce a novel and complementary concept related to electronic transitions with the canonical transition density matrix and the canonical transition orbitals. Their expression actually reflects the electronic cloud polarisation in the orbital space with a decomposition based on the actual contribution of one-particle excitations from occupied canonical orbitals to virtual ones. This approach validates our novel interpretation of the transition density matrix elements in terms of the Euclidean norm of elementary transition vectors in a linear tensor space. A proper use of these new concepts leads to the conclusion that despite the different principles underlying their construction, they provide two equivalent excited states topological analyses. This connexion is evidenced through simple illustrations of (in)organic dyes electronic transitions analysis.

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

A general theory of linear cosmological perturbations: scalar-tensor and vector-tensor theories

We present a method for parametrizing linear cosmological perturbations of theories of gravity, around homogeneous and isotropic backgrounds. The method is sufficiently general and systematic that it can be applied to theories with any degrees of freedom (DoFs) and arbitrary gauge symmetries. In this paper, we focus on scalar-tensor and vector-tensor theories, invariant under linear coordinate transformations. In the case of scalar-tensor theories, we use our framework to recover the simple parametrizations of linearized Horndeski and ''Beyond Horndeski'' theories, and also find higher-derivative corrections. In the case of vector-tensor theories, we first construct the most general quadratic action for perturbations that leads to second-order equations of motion, which propagates two scalar DoFs. Then we specialize to the case in which the vector field is time-like (à la Einstein-Aether gravity), where the theory only propagates one scalar DoF. As a result, we identify the complete forms of the quadratic actions for perturbations, and the number of free parameters that need to be defined, to cosmologically characterize these two broad classes of theories.

79 ASTRONOMY AND ASTROPHYSICS↗