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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Gravitational form factors of charmonia

We investigate the gravitational form factors of charmonium. Our method is based on a Hamiltonian formalism on the light front known as basis light-front quantization. The charmonium mass spectrum and light-front wave functions were obtained from diagonalizing an effective Hamiltonian that incorporates confinement from holographic QCD and one-gluon exchange interaction from light-front QCD. We proposed a quantum many-body approach to construct the hadronic matrix elements of the energy momentum tensor T + + and T + − , which are used to extract the gravitational form factors A ( Q 2 ) and D ( Q 2 ) . The obtained form factors satisfy the known constraints, e.g., the von Laue condition. From these quantities, we also extract the energy, pressure and light-front energy distributions of the system. We find that hadrons are multilayer systems. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2015

Symmetry breaking indication for supergravity inflation in light of the Planck 2015

Supergravity (SUGRA) theories with exact global U(1) symmetry or shift symmetry in Kähler potential provide natural frameworks for inflation. However, quadratic inflation is disfavoured by the new results on primordial tensor fluctuations from the Planck Collaboration. To be consistent with the new Planck data, we point out that the explicit symmetry breaking is needed, and study these two SUGRA inflation in detail. For SUGRA inflation with global U(1) symmetry, the symmetry breaking term leads to a trigonometric modulation on inflaton potential. Coefficient of the U(1) symmetry breaking term is of order 10{sup −2}, which is sufficient large to improve the inflationary predictions while its higher order corrections are negligible. Such models predict sizeable tensor fluctuations and highly agree with the Planck results. In particular, the model with a linear U(1) symmetry breaking term predicts the tensor-to-scalar ratio around r∼0.01 and running spectral index α{sub s∼} −0.004, which comfortably fit with the Planck observations. For SUGRA inflation with breaking shift symmetry, the inflaton potential is modulated by an exponential factor. The modulated linear and quadratic models are consistent with the Planck observations. In both types of models the tensor-to-scalar ratio can be of order 10{sup −2}, which will be tested by the near future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Primordial perturbations from dilaton-induced gauge fields

We study the primordial scalar and tensor perturbations in inflation scenario involving a spectator dilaton field. In our setup, the rolling spectator dilaton causes a tachyonic instability of gauge fields, leading to a copious production of gauge fields in the superhorizon regime, which generates additional scalar and tensor perturbations through gravitational interactions. Our prime concern is the possibility to enhance the tensor-to-scalar ratio r relative to the standard result, while satisfying the observational constraints. To this end, we allow the dilaton field to be stabilized before the end of inflation, but after the CMB scales exit the horizon. We show that for the inflaton slow roll parameter ε ∼> 10{sup −3}, the tensor-to-scalar ratio in our setup can be enhanced only by a factor of O(1) compared to the standard result. On the other hand, for smaller ε corresponding to a lower inflation energy scale, a much bigger enhancement can be achieved, so that our setup can give rise to an observably large r∼> 10{sup −2} even when ε|| 10{sup −3}. The tensor perturbation sourced by the spectator dilaton can have a strong scale dependence, and is generically red-tilted. We also discuss a specific model to realize our scenario, and identify the parameter region giving an observably large r for relatively low inflation energy scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

The principal component analysis method used with polynomial Chaos expansion to propagate uncertainties through critical transport problems

A critical neutron transport problem with random material properties is introduced. The total cross section and the average neutron multiplicity are assumed to be uncertain, characterized by the mean and variance with a log-normal distribution. The average neutron multiplicity and the total cross section are assumed to be uncorrected and the material properties for differing materials are also assumed to be uncorrected. The principal component analysis method is used to decompose the covariance matrix into eigenvalues and eigenvectors and then 'realizations' of the material properties can be computed. A simple Monte Carlo brute force sampling of the decomposed covariance matrix is employed to obtain a benchmark result for each test problem. In order to save computational time and to characterize the moments and probability density function of the multiplication factor the polynomial chaos expansion method is employed along with the stochastic collocation method. A Gauss-Hermite quadrature set is convolved into a multidimensional tensor product quadrature set and is successfully used to compute the polynomial chaos expansion coefficients of the multiplication factor. Finally, for a particular critical fuel pin assembly the appropriate number of random variables and polynomial expansion order are investigated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Experimental evaluation of electrical conductivity imaging of anisotropic brain tissues using a combination of diffusion tensor imaging and magnetic resonance electrical impedance tomography

Anisotropy of biological tissues is a low-frequency phenomenon that is associated with the function and structure of cell membranes. Imaging of anisotropic conductivity has potential for the analysis of interactions between electromagnetic fields and biological systems, such as the prediction of current pathways in electrical stimulation therapy. To improve application to the clinical environment, precise approaches are required to understand the exact responses inside the human body subjected to the stimulated currents. In this study, we experimentally evaluate the anisotropic conductivity tensor distribution of canine brain tissues, using a recently developed diffusion tensor-magnetic resonance electrical impedance tomography method. At low frequency, electrical conductivity of the biological tissues can be expressed as a product of the mobility and concentration of ions in the extracellular space. From diffusion tensor images of the brain, we can obtain directional information on diffusive movements of water molecules, which correspond to the mobility of ions. The position dependent scale factor, which provides information on ion concentration, was successfully calculated from the magnetic flux density, to obtain the equivalent conductivity tensor. By combining the information from both techniques, we can finally reconstruct the anisotropic conductivity tensor images of brain tissues. The reconstructed conductivity images better demonstrate the enhanced signal intensity in strongly anisotropic brain regions, compared with those resulting from previous methods using a global scale factor.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Emergent classical spacetime from microstates of an incipient black hole

Black holes have an enormous underlying space of microstates, but universal macroscopic physics characterized by mass, charge and angular momentum as well as a causally disconnected interior. This leads two related puzzles: (1) How does the effective factorization of interior and exterior degrees of freedom emerge in gravity?, and (2) How does the underlying degeneracy of states wind up having a geometric realization in the horizon area and in properties of the singularity? We explore these puzzles in the context of an incipient black hole in the AdS/CFT correspondence, the microstates of which are dual to half-BPS states of the $\mathcal{N}$ = 4 super-Yang-Mills theory. First, we construct a code subspace for this black hole and show how to organize it as a tensor product of a universal macroscopic piece (describing the exterior), and a factor corresponding to the microscopic degrees of freedom (describing the interior). We then study the classical phase space and symplectic form for low-energy excitations around the black hole. On the AdS side, we find that the symplectic form has a new physical degree of freedom at the stretched horizon of the black hole, reminiscent of soft hair, which is absent in the microstates. We explicitly show how such a soft mode emerges from the microscopic phase space in the dual CFT via a canonical transformation and how it encodes partial information about the microscopic degrees of freedom of the black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Trace anomaly form factors from lattice QCD

The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to Q 2 ∼ 4.3 GeV 2 and the nucleon up to Q 2 ∼ 1 GeV 2 . The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from ∼ 250 to ∼ 540 MeV , including the unitary point ∼ 340 MeV . We calculate the radius of the glue trace anomaly for the pion and the nucleon from the z expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion’s form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Inflation alternative via the gravitational field of a singularity

Abstract We explore the scenario that the observable Universe emerged from the vicinity of a negative mass ring singularity, and all content of the Universe travels at the same group velocity close to the speed of light on a geodesic trajectory along the axis of rotation of the singularity. In appropriate coordinate parametrization and evaluated on the trajectory, we find that the metric tensor in the vicinity of the trajectory exhibits a conformal scale factor a ( η ) with contraction and subsequent expansion properties that solve the horizon problem. We then introduce a static flow of gravitating radiation along the trajectory (perturbatively with respect to the mass scale of the singularity) to model a homogeneous radiation dominated Universe. Solving the Einstein field equations with a physically motivated ansatz of metric perturbation then reveals that the effective conformal scale factor indeed grows asymptotically with the same power law as expected in a conventional radiation dominated Universe.

Astronomy & Astrophysics↗
DOE OSTI2018

Plasma equilibrium with fast ion orbit width, pressure anisotropy, and toroidal flow effects

Here, we formulate the problem of tokamak plasma equilibrium including the toroidal flow and fast ion (or energetic particle, EP) pressure anisotropy and the finite drift orbit width (FOW) effects. The problem is formulated via the standard Grad-Shafranov equation (GShE) amended by the solvability condition which imposes physical constraints on allowed spacial dependencies of the anisotropic pressure. The GShE problem employs the pressure coupling scheme and includes the dominant diagonal terms and non-diagonal corrections to the standard pressure tensor. The anisotropic tensor elements are obtained via the distribution function represented in the factorized form via the constants of motion. Considered effects on the plasma equilibrium are estimated analytically, if possible, to understand their importance for GShE tokamak plasma problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonperturbative running of the tensor operator for N f = 3 QCD from the chirally rotated Schrödinger functional

We study the renormalization group (RG) running of the nonsinglet tensor operator, for N f = 3 QCD with Wilson fermions in a mixed action setup, with standard Schrödinger functional (SF) boundary conditions for sea quarks and chirally rotated Schrödinger functional ( χ SF ) boundary conditions for valence quarks. Based on a recursive finite-size scaling technique we compute nonperturbatively the tensor step-scaling function for an energy range between a hadronic scale and an electroweak scale, above which perturbation theory may be safely applied. Our result is expressed as the RG-running factor T RGI / [ T ( μ had ) ] R , where the numerator is the scale independent (renormalization group invariant—RGI) tensor operator and the denominator is its renormalized counterpart at a hadronic scale μ had = 233 ( 8 ) MeV in a given scheme. We determine the step-scaling function in four distinct renormalization schemes. We also compute the renormalization parameters of these schemes at μ had which, combined with the RG-running factor, gives the scheme-independent quantity Z T RGI ( g 0 2 ) in four schemes and for a range of bare gauge couplings in which large volume hadronic matrix element simulations are performed by the CLS consortium in N f = 2 + 1 QCD. All four results are compatible and also agree with a recent determination based on a unitary setup for Wilson quarks with Schrödinger functional boundary conditions [arXiv:2309.04314]. This provides a strong universality test. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2011

On the magnetostatics of chains of magnetic nanoparticles.

A novel approach is presented for the computation of the magnetostatic energy of straight and bent chains of identical, uniformly magnetized particles of arbitrary shape. The formalism relies on the concept of the magnetometric tensor field, and allows for closed form expressions for the magnetostatic energy, demagnetization factor, Young's modulus, and bending modulus of chains in terms of the shape amplitude of the particles. Analytical solutions are presented for straight chains of spheres, cubes, and cylinders, and for bent chains of spheres. Numerical results include chains of octahedra, tetrahedra, cuboctahedra, and bi-cones. The axial demagnetization factor for the bi-cone shape is derived in analytical form. An approximate energy expression, using the full shape-dependent interaction formalism for short separation distances, and the standard dipolar interaction expression for larger distances, is introduced.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Validating quantum-classical programming models with tensor network simulations

The exploration of hybrid quantum-classical algorithms and programming models on noisy near-term quantum hardware has begun. As hybrid programs scale towards classical intractability, validation and benchmarking are critical to understanding the utility of the hybrid computational model. In this paper, we demonstrate a newly developed quantum circuit simulator based on tensor network theory that enables intermediate-scale verification and validation of hybrid quantum-classical computing frameworks and programming models. We present our tensor-network quantum virtual machine (TNQVM) simulator which stores a multi-qubit wavefunction in a compressed (factorized) form as a matrix product state, thus enabling single-node simulations of larger qubit registers, as compared to brute-force statevector simulators. Our simulator is designed to be extensible in both the tensor network form and the classical hardware used to run the simulation (multicore, GPU, distributed). The extensibility of the TNQVM simulator with respect to the simulation hardware type is achieved via a pluggable interface for different numerical backends (e.g., ITensor and ExaTENSOR numerical libraries). We demonstrate the utility of our TNQVM quantum circuit simulator through the verification of randomized quantum circuits and the variational quantum eigensolver algorithm, both expressed within the eXtreme-scale ACCelerator (XACC) programming model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Implementation of a generalized coherent elastic scattering formulation for thermal neutron scattering analysis - 14487

The scattering interactions of thermal neutrons in matter are classified into inelastic and elastic components. In crystalline materials, coherent elastic scattering represents a dominant form of interaction. Traditionally, the treatment of this interaction in standard computer codes (e.g., the LEAPR module of the NJOY code system) was performed within the cubic approximation and restricted the analyst to a few materials that were hardwired into the code. In general, these restrictions can hinder the production of ENDF style thermal neutron scattering libraries for new materials outside the current portfolio of treated materials. In this work, a general formulation based on the exact theory of coherent elastic scattering was developed. This formulation is capable of processing the crystalline structure of any material. The cubic approximation is relaxed by accounting for the partial forces from asymmetric crystal structures in the calculation. The atom sites approximation is relaxed by utilizing a site-wise Debye-Waller factor. Specifically, the Debye-Waller factor is no longer regarded as a uniform constant across the crystal but treated as a matrix tensor that uses the integral of the partial phonon density of states. Therefore, the calculation of the coherent elastic scattering cross section is executed by applying the matrix-form of the Debye-Waller factor on each reciprocal space point and summing them up. Furthermore, the developed algorithm supports compounds with any number of elements and of any ratio to each other. This general formulation was successfully implemented and contributed to the production of ENDF style thermal neutron scattering libraries for silicon dioxide (α-SiO{sub 2}) and silicon carbide (3C-SiC). The needed input for the analysis was generated using density functional theory (DFT) atomistic simulations and dynamical matrix calculations for both materials. (authors)

42 ENGINEERING↗
DOE OSTI2012

Light-Front Holographic QCD

The relation between the hadronic short-distance constituent quark and gluon particle limit and the long-range confining domain is yet one of the most challenging aspects of particle physics due to the strong coupling nature of Quantum Chromodynamics, the fundamental theory of the strong interactions. The central question is how one can compute hadronic properties from first principles; i.e., directly from the QCD Lagrangian. The most successful theoretical approach thus far has been to quantize QCD on discrete lattices in Euclidean space-time. Lattice numerical results follow from computation of frame-dependent moments of distributions in Euclidean space and dynamical observables in Minkowski spacetime, such as the time-like hadronic form factors, are not amenable to Euclidean lattice computations. The Dyson-Schwinger methods have led to many important insights, such as the infrared fixed point behavior of the strong coupling constant, but in practice, the analyses are limited to ladder approximation in Landau gauge. Baryon spectroscopy and the excitation dynamics of nucleon resonances encoded in the nucleon transition form factors can provide fundamental insight into the strong-coupling dynamics of QCD. New theoretical tools are thus of primary interest for the interpretation of the results expected at the new mass scale and kinematic regions accessible to the JLab 12 GeV Upgrade Project. The AdS/CFT correspondence between gravity or string theory on a higher-dimensional anti-de Sitter (AdS) space and conformal field theories in physical space-time has led to a semiclassical approximation for strongly-coupled QCD, which provides physical insights into its nonperturbative dynamics. The correspondence is holographic in the sense that it determines a duality between theories in different number of space-time dimensions. This geometric approach leads in fact to a simple analytical and phenomenologically compelling nonperturbative approximation to the full light-front QCD Hamiltonian 'Light-Front Holography'. Light-Front Holography is in fact one of the most remarkable features of the AdS/CFT correspondence. The Hamiltonian equation of motion in the light-front (LF) is frame independent and has a structure similar to eigenmode equations in AdS space. This makes a direct connection of QCD with AdS/CFT methods possible. Remarkably, the AdS equations correspond to the kinetic energy terms of the partons inside a hadron, whereas the interaction terms build confinement and correspond to the truncation of AdS space in an effective dual gravity approximation. One can also study the gauge/gravity duality starting from the bound-state structure of hadrons in QCD quantized in the light-front. The LF Lorentz-invariant Hamiltonian equation for the relativistic bound-state system is P{sub {mu}}P{sup {mu}}|{psi}(P)> = (P{sup +}P{sup -} - P{sub {perpendicular}}{sup 2})|{psi}(P)> = M{sup 2}|{psi}(P)>, P{sup {+-}} = P{sup 0} {+-} P{sup 3}, where the LF time evolution operator P{sup -} is determined canonically from the QCD Lagrangian. To a first semiclassical approximation, where quantum loops and quark masses are not included, this leads to a LF Hamiltonian equation which describes the bound-state dynamics of light hadrons in terms of an invariant impact variable {zeta} which measures the separation of the partons within the hadron at equal light-front time {tau} = x{sup 0} + x{sup 3}. This allows us to identify the holographic variable z in AdS space with an impact variable {zeta}. The resulting Lorentz-invariant Schroedinger equation for general spin incorporates color confinement and is systematically improvable. Light-front holographic methods were originally introduced by matching the electromagnetic current matrix elements in AdS space with the corresponding expression using LF theory in physical space time. It was also shown that one obtains identical holographic mapping using the matrix elements of the energy-momentum tensor by perturbing the AdS metric around its static solution. A gravity dual to QCD is not known, but the mechanisms of confinement can be incorporated in the gauge/gravity correspondence by modifying the AdS geometry in the large infrared (IR) domain z {approx} 1 = {Lambda}{sub QCD}, which also sets the scale of the strong interactions. In this simplified approach we consider the propagation of hadronic modes in a fixed effective gravitational background asymptotic to AdS space, which encodes salient properties of the QCD dual theory, such as the ultraviolet (UV) conformal limit at the AdS boundary, as well as modifications of the background geometry in the large z IR region to describe confinement. The modified theory generates the point-like hard behavior expected from QCD, instead of the soft behavior characteristic of extended objects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

An infinite swampland of U(1) charge spectra in 6D supergravity theories

We analyze the anomaly constraints on 6D supergravity theories with a single abelian U(1) gauge factor. For theories with charges restricted to q = ±1, ±2 and no tensor multiplets, anomaly-free models match those models that can be realized from F-theory compactifications almost perfectly. For theories with tensor multiplets or with larger charges, the F-theory constraints are less well understood. We show, however, that there is an infinite class of distinct massless charge spectra in the “swampland” of theories that satisfy all known quantum consistency conditions but do not admit a realization through F-theory or any other known approach to string compactification. We also compare the spectra of charged matter in abelian theories with those that can be realized from breaking nonabelian SU(2) and higher rank gauge symmetries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Magnetizability and rotational g tensors for density fitted local second-order Møller-Plesset perturbation theory using gauge-including atomic orbitals

In this paper, we present theory and implementation of an efficient program for calculating magnetizabilities and rotational g tensors of closed-shell molecules at the level of local second-order Møller-Plesset perturbation theory (MP2) using London orbitals. Density fitting is employed to factorize the electron repulsion integrals with ordinary Gaussians as fitting functions. The presented program for the calculation of magnetizabilities and rotational g tensors is based on a previous implementation of NMR shielding tensors reported by S. Loibl and M. Schütz [J. Chem. Phys. 137, 084107 (2012)]. Extensive test calculations show (i) that the errors introduced by density fitting are negligible, and (ii) that the errors of the local approximation are still rather small, although larger than for nuclear magnetic resonance (NMR) shielding tensors. Electron correlation effects for magnetizabilities are tiny for most of the molecules considered here. MP2 appears to overestimate the correlation contribution of magnetizabilities such that it does not constitute an improvement over Hartree-Fock (when comparing to higher-order methods like CCSD(T)). For rotational g tensors the situation is different and MP2 provides a significant improvement in accuracy over Hartree-Fock. The computational performance of the new program was tested for two extended systems, the larger comprising about 2200 basis functions. It turns out that a magnetizability (or rotational g tensor) calculation takes about 1.5 times longer than a corresponding NMR shielding tensor calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗