Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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↗

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↗

A Cartesian parametrization for the numerical analysis of material instability

We examine four parametrizations of the unit sphere in the context of material stability analysis by means of the singularity of the acoustic tensor. We then propose a Cartesian parametrization for vectors that lie a cube of side length two and use these vectors in lieu of unit normals to test for the loss of the ellipticity condition. This parametrization is then used to construct a tensor akin to the acoustic tensor. It is shown that both of these tensors become singular at the same time and in the same planes in the presence of a material instability. Furthermore, the performance of the Cartesian parametrization is compared against the other parametrizations, with the results of these comparisons showing that in general, the Cartesian parametrization is more robust and more numerically efficient than the others.

36 MATERIALS SCIENCE↗

Results and forecasts on cosmic inflation from weak lensing

Here, we highlight the role of weak lensing measurements from current and upcoming stage-IV imaging surveys in the search for cosmic inflation, specifically in measuring the scalar spectral index n s . To do so, we combine the Dark Energy Survey 3 years of observation weak lensing and clustering data with Bicep/Keck, Planck, and Sloan Digital Sky Survey data in r⁢Λ⁢ CDM (cold dark matter) where r is the tensor-to-scalar ratio. While there is no significant improvement in constraining power, we obtain a 1⁢σ shift on n s . Additionally, we forecast a weak lensing and clustering data vector from the 10-year Legacy Survey of Space and Time by the Vera C. Rubin Observatory and show its combination with current data would improve their n s constraints by 25% in r⁢Λ⁢ CDM.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

Interactive color graphical postprocessing as a unifying influence in numerical analysis research

The general capabilities and the hardware environment of an interactive color graphical postprocessor for three-dimensional numerical analysis are described as well as the synergistic influence of this aspect of interactive graphics on the research environment. Areas of practical concern include finite element analysis, fracture propagation, surface and framed structures. The interactive color graphical postprocessor, POSTPRO3D, is capable of displaying the scalar, vector, and tensor results of finite difference, finite element, and boundary element analysis. In the case of fatigue cracking in steel crane-rail girders, it was shown that POSTPRO3D offers a high level of flexibility, interactivity and features which minimize response time.

Abel, J. F.↗

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MAGNETIC HELICITY AND ENERGY SPECTRA OF A SOLAR ACTIVE REGION

We compute for the first time the magnetic helicity and energy spectra of the solar active region NOAA 11158 during 2011 February 11-15 at 20° southern heliographic latitude using observational photospheric vector magnetograms. We adopt the isotropic representation of the Fourier-transformed two-point correlation tensor of the magnetic field. The sign of the magnetic helicity turns out to be predominantly positive at all wavenumbers. This sign is consistent with what is theoretically expected for the southern hemisphere. The magnetic helicity normalized to its theoretical maximum value, here referred to as relative helicity, is around 4% and strongest at intermediate wavenumbers of k ≈ 0.4 Mm{sup –1}, corresponding to a scale of 2π/k ≈ 16 Mm. The same sign and a similar value are also found for the relative current helicity evaluated in real space based on the vertical components of magnetic field and current density. The modulus of the magnetic helicity spectrum shows a k {sup –11/3} power law at large wavenumbers, which implies a k {sup –5/3} spectrum for the modulus of the current helicity. A k {sup –5/3} spectrum is also obtained for the magnetic energy. The energy spectra evaluated separately from the horizontal and vertical fields agree for wavenumbers below 3 Mm{sup –1}, corresponding to scales above 2 Mm. This gives some justification to our assumption of isotropy and places limits resulting from possible instrumental artifacts at small scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic fluid moments closure for a magnetized plasma with collisions

A novel method aimed at a kinetic moments closure for a magnetized plasma with collisions is proposed. The velocity distribution function for each species is expanded in 8 Gaussian Radial Basis Functions (GRBFs) which are essentially shifted Maxwellians at eight representative 3D-velocity points of drift. The vector of 8 fluid moments (for particle density, 3 particle fluxes, total energy density, and 3 energy fluxes) has an 8x8 analytic linear matrix relation to the vector of 8 GRBF density weighs in 3D-real space. The 8 fluid moments with sources for each species are advanced in time while the 8 GRBF weighs are determined from the 8x8 inverse matrix. The two closure moments (for the stress tensor and the energy weighted stress tensor) are linearly determined from the GRBF weighs. Most importantly the velocity moments of the nonlinear Coulomb Fokker-Planck collision operator [Rosenbluth et al, Phys. Rev 107, 1957] are evaluated from the GRBF weights. Generalization from 8 to 12, 16, 20 .., in an energy weighted moment hierarchy is straightforward. The electric field follows from a generalized vorticity (quasi-neutral charge conservation) equation. A strong drift ordering approximation can be applied to eliminate any spuriously unstable high frequency cyclotron motions. A novel weak drift ordering two-time step scheme avoids the vorticity equation by following ion cyclotron motion in time to get the electric field with ion gyroaveraging. Inclusion of low-beta magnetic perturbations is straightforward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic quantum effects in anti-de Sitter spacetime

The two-point functions of the vector potential and of the field tensor for the electromagnetic field in the background of anti-de Sitter (AdS) spacetime are evaluated in an arbitrary number of spatial dimensions. First, we consider the two-point functions in the boundary-free geometry and then generalize the results in the presence of a reflecting boundary parallel to the AdS horizon. By using the expressions for the two-point functions of the field tensor, we investigate the vacuum expectation values of the electric field squared and of the energy–momentum tensor. Simple asymptotic expressions are provided near the AdS boundary and horizon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermal Radiation Transport with Tensor Trains

We present a novel tensor network algorithm to solve the time-dependent, gray thermal radiation transport equation. The method invokes a tensor train (TT) decomposition for the specific intensity. The efficiency of this approach is dictated by the rank of the decomposition. When the solution is “low rank,” the memory footprint of the specific intensity solution vector may be significantly compressed. The algorithm, following a step-then-truncate approach of a traditional discrete ordinates method, operates directly on the compressed state vector, thereby enabling large speedups for low-rank solutions. To achieve these speedups, we rely on a recently developed rounding approach based on the Gram-SVD. We detail how familiar S N algorithms for (gray) thermal transport can be mapped to this TT framework and present several numerical examples testing both the optically thick and thin regimes. The TT framework finds low-rank structure and supplies up to ≃60× speedups and ≃1000× compressions for problems demanding large angle counts, thereby enabling previously intractable SN calculations and supplying a promising avenue to mitigate ray effects.

79 ASTRONOMY AND ASTROPHYSICS↗

Finite-strain large-deflection elastic-viscoplastic finite-element transient response analysis of structures

A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.

Rodal, J. J. A.↗

Modified Einstein and Navier–Stokes Equations

The appearance of inertial rest mass-energy is associated with the kinematic slowing-down of time and with the vortex state of the elementary massive space with zero integral of its kinetic and potential energies. An analog of the Einstein equation is found for moving densities of a non-empty metric space in the concept of the Einstein–Infeld material field. The vector consequences of this tensor equation for a metric medium of overlapping elementary carriers of continuous mass-energies allow us to modify the Navier–Stokes equation under inertial motion of the matter of the nonlocal field in the nonrelativistic limit. The nonlocality of massenergy generates kinematic accelerations of feedback to Newtonian acceleration, which impedes asymptotic divergence of energy fluxes. Stabilization of inertial media by dynamic Bernoulli pressure corresponds to nonlocal self-organization of Einstein–Infeld non-empty space and invalidates Newtonian localization of masses in empty space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cosmology for quadratic gravity in generalized Weyl geometry

A class of vector-tensor theories arises naturally in the framework of quadratic gravity in spacetimes with linear vector distortion. Requiring the absence of ghosts for the vector field imposes an interesting condition on the allowed connections with vector distortion: the resulting one-parameter family of connections generalises the usual Weyl geometry with polar torsion. The cosmology of this class of theories is studied, focusing on isotropic solutions wherein the vector field is dominated by the temporal component. De Sitter attractors are found and inhomogeneous perturbations around such backgrounds are analysed. In particular, further constraints on the models are imposed by excluding pathologies in the scalar, vector and tensor fluctuations. Various exact background solutions are presented, describing a constant and an evolving dark energy, a bounce and a self-tuning de Sitter phase. However, the latter two scenarios are not viable under a closer scrutiny.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

IIB supergravity and the E 6(6) covariant vector-tensor hierarchy

IIB supergravity is reformulated with a manifest local USp(8) invariance that makes the embedding of five-dimensional maximal supergravities transparent. In this formulation the ten-dimensional theory exhibits all the 27 one-form fields and 22 of the 27 two-form fields that are required by the vector-tensor hierarchy of the five-dimensional theory. The missing 5 two-form fields must transform in the same representation as a descendant of the ten-dimensional ‘dual graviton’. The invariant E 6(6) symmetric tensor that appears in the vector-tensor hierarchy is reproduced. Generalized vielbeine are derived from the supersymmetry transformations of the vector fields, as well as consistent expressions for the USp(8) covariant fermion fields. Implications are further discussed for the consistency of the truncation of IIB supergravity compactified on the five-sphere to maximal gauged supergravity in five space-time dimensions with an SO(6) gauge group.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Parity-violating CMB correlators with non-decaying statistical anisotropy

We examine the effect induced on cosmological correlators by the simultaneous breaking of parity and of statistical isotropy. As an example of this, we compute the scalar-scalar, scalar-tensor, tensor-tensor and scalar-scalar-scalar cosmological correlators in presence of the coupling L = f(φ) ( − 1/4 F{sup 2} + γ/4 F ∼F ) between the inflaton φ and a vector field with vacuum expectation value A. For a suitably chosen function f, the energy in the vector field ρ{sub A} does not decay during inflation. This results in nearly scale-invariant signatures of broken statistical isotropy and parity. Specifically, we find that the scalar-scalar correlator of primordial curvature perturbations includes a quadrupolar anisotropy, P{sub ζ}(k) = P(k)[1+g{sub *}( k-circumflex ⋅Â){sup 2}], and a (angle-averaged) scalar bispectrum that is a linear combination of the first 3 Legendre polynomials, B{sub ζ}(k{sub 1}, k{sub 2}, k{sub 3}) = ∑{sub L} c{sub L} P{sub L} ( k-circumflex {sub 1} ⋅ k-circumflex {sub 2}) P(k{sub 1}) P(k{sub 2}) + 2 perms , with c{sub 0}:c{sub 1}:c{sub 2}=2-3:1 (c{sub 1}≠0 is a consequence of parity violation, corresponding to the constant 0γ ≠ ). The latter is one of the main results of this paper, which provides for the first time a clear example of an inflationary model where a non-negligible c{sub 1} contribution to the bispectrum is generated. The scalar-tensor and tensor-tensor correlators induce characteristic signatures in the Cosmic Microwave Background temperature anisotropies (T) and polarization (E/B modes); namely, non-diagonal contributions to (a{sub ℓ1m1}a{sup *}{sub ℓ2m2}), with |ℓ{sub 1} − ℓ{sub 2}| = 1 in TT, TE, EE and BB, and |ℓ{sub 1} − ℓ{sub 2}| = 2 in TB and EB. The latest CMB bounds on the scalar observables (g{sub *}, c{sub 0}, c{sub 1} and c{sub 2}), translate into the upper limit ρ{sub A} / ρ{sub φ} ∼< 10{sup −9} at 0γ=. We find that the upper limit on the vector energy density becomes much more stringent as γ grows.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗

Adiabatic fast passage spin manipulation measurements in solid polarized targets

Adiabatic fast passage (AFP) is a rapid method for reversing nuclear polarization and manipulating spin populations in polarized solid targets, avoiding the long repolarization times associated with dynamic nuclear polarization (DNP). We report AFP measurements in a 5 T, 1 K polarized-target system for irradiated 15 NH 3 , irradiated 14 ND 3 , and butanol-based materials prepared either with TEMPO doping or by irradiation. We also present a joint manipulated-lineshape analysis for spin-1 targets and demonstrate that vector and tensor polarizations can be extracted from AFP-manipulated deuteron NMR spectra even when the populations are not described by a single Boltzmann spin temperature. In conclusion, we report a reproducible polarization- and direction-dependent AFP response in a large irradiated 15 NH 3 sample. These ammonia results are presented as empirical observations under the specific sample–coil conditions of the experiment, with possible circuit-mediated mechanisms such as radiation damping or superradiant behavior discussed but not assigned as a definitive cause.

Adiabatic fast passage↗