Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 271 records · Page 15

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.↗

Quasinormal modes for non-minimally coupled scalar fields in regular black hole spacetimes: Grey-body factors, area spectrum and shadow radius

Black hole perturbation theory is a useful approach to study interactions between black holes and fundamental fields. A particular class of black hole solutions arising out of modification of Einstein’s general theory of relativity are regular black holes (RBHs) which can be constructed using a nonlinear electrodynamic Lagrangian. Because of their importance, we are interested in studying the behavior of three kinds of such RBHs under perturbations generated by an external field. Indeed, we investigate the quasinormal modes (QNMs) of a massive scalar field propagating near the RBHs which is non-minimally coupled to the Ricci scalar-tensor of background geometry. We will attempt to find the low-lying quasinormal frequencies of the perturbations by using WKB approximation. We shall also study the relationship between the QNMs and some characteristic properties of black holes such as grey-body factors, area quantization, and shadow radius.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A continuous analogue of the tensor-train decomposition

We develop new approximation algorithms and data structures for representing and computing with multivariate functions using the functional tensor-train (FT), a continuous extension of the tensor-train (TT) decomposition. The FT represents functions using a tensor-train ansatz by replacing the three-dimensional TT cores with univariate matrix-valued functions. The main contribution of this paper is a framework to compute the FT that employs adaptive approximations of univariate fibers, and that is not connected to any tensorized discretization. The algorithm can be coupled with any univariate linear or nonlinear approximation procedure. We demonstrate that this approach can generate multivariate function approximations that are several orders of magnitude more accurate, for the same cost, than those based on the conventional approach of compressing the coefficient tensor of a tensor-product basis. Our approach is in the spirit of other continuous computation packages such as Chebfun, and yields an algorithm which requires the computation of “continuous” matrix factorizations such as the LU and QR decompositions of vector-valued functions. To support these developments, we describe continuous versions of an approximate maximum-volume cross approximation algorithm and of a rounding algorithm that re-approximates an FT by one of lower ranks. We demonstrate that our technique improves accuracy and robustness, compared to TT and quantics-TT approaches with fixed parameterizations, of high-dimensional integration, differentiation, and approximation of functions with local features such as discontinuities and other nonlinearities.

42 ENGINEERING↗

Gravitational form factors of light mesons

Here, we calculate the gravitational form factors of the pion, $\sigma$ meson, and $\rho$ meson in the Nambu-Jona-Lasinio model of quantum chromodynamics. The canonical energy-momentum tensor (EMT) is used in their derivation, allowing the possibility of an antisymmetric contribution when the hadron has intrinsic spin. We show that the asymmetric graviton vertex arising from the canonical EMT satisfies a simpler Ward-Takahashi identity (WTI) than the symmetric graviton vertex of the Belinfante EMT. The necessity of fully dressing the graviton vertex through the relevant Bethe-Salpeter equation is demonstrated for observing both the WTI and a low-energy pion theorem. Lastly, we calculate static moments of the meson EMT decompositions, obtaining predictions for the meson mass radii. We find light cone mass radii of 0.27 fm for the pion, 0.32 fm for the $\sigma$, and 0.32 fm for the $\rho$. For the pion and rho, these are smaller than the light cone charge radii, respectively 0.51 fm and 0.45 fm, while we have a $\sigma$ charge radius of zero. Our light cone pion mass radius agrees with a phenomenological extraction from KEKB data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Polycrystalline Mg{sub 2}Si thin films: A theoretical investigation of their electronic transport properties

The electronic structures and thermoelectric properties of a polycrystalline Mg{sub 2}Si thin film have been investigated by first-principle density-functional theory (DFT) and Boltzmann transport theory calculations within the constant-relaxation time approximation. The polycrystalline thin film has been simulated by assembling three types of slabs each having the orientation (001), (110) or (111) with a thickness of about 18 Å. The effect of applying the relaxation procedure to the thin film induces disorder in the structure that has been ascertained by calculating radial distribution functions. For the calculations of the thermoelectric properties, the energy gap has been fixed at the experimental value of 0.74 eV. The thermoelectric properties, namely the Seebeck coefficient, the electrical conductivity and the power factor, have been determined at three temperatures of 350 K, 600 K and 900 K with respect to both the energy levels and the p-type and n-type doping levels. The best Seebeck coefficient is obtained at 350 K: the S{sub yy} component of the tensor amounts to about ±1000 μV K{sup −1}, depending on the type of charge carriers. However, the electrical conductivity is much too small which results in low values of the figure of merit ZT. Structure–property relationship correlations based on directional radial distribution functions allow us to tentatively draw some explanations regarding the anisotropy of the electrical conductivity. Finally, the low ZT values obtained for the polycrystalline Mg{sub 2}Si thin film are paralleled with those recently reported in the literature for bulk chalcogenide glasses. - Graphical abstract: Structure of the polycrystalline thin film of Mg{sub 2}Si. - Author-Highlights: • Polycrystalline Mg{sub 2}Si film has been modelled by DFT approach. • Thermoelectric properties have been evaluated by semi-classical Boltzmann theory. • The structure was found to be slightly disordered after relaxation. • The highest value of Seebeck coefficient reaches −1000 μV/K at 350 K. • The ZT value is estimated between 0.0175 at 350 K and 0.2 at 900 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

On the prior dependence of constraints on the tensor-to-scalar ratio

We investigate the prior dependence of constraints on cosmic tensor perturbations. Commonly imposed is the strong prior of the single-field inflationary consistency equation, relating the tensor spectral index n{sub T} to the tensor-to-scalar ratio r. Dropping it leads to significantly different constraints on n{sub T}, with both positive and negative values allowed with comparable likelihood, and substantially increases the upper limit on r on scales k = 0.01 Mpc{sup -1} to 0.05 Mpc{sup -1}, by a factor of ten or more. Even if the consistency equation is adopted, a uniform prior on r on one scale does not correspond to a uniform one on another; constraints therefore depend on the pivot scale chosen. We assess the size of this effect and determine the optimal scale for constraining the tensor amplitude, both with and without the consistency relation.

79 ASTRONOMY AND ASTROPHYSICS↗

Twisted bilayer graphene. V. Exact analytic many-body excitations in Coulomb Hamiltonians: Charge gap, Goldstone modes, and absence of Cooper pairing

We find exact analytic expressions for the energies and wave functions of the charged and neutral excitations above the exact ground states (at rational filling per unit cell) of projected Coulomb Hamiltonians in twisted bilayer graphene. Our exact expressions are valid for any form of the Coulomb interaction and any form of A A and A B / B A tunneling. The single charge excitation energy is a convolution of the Coulomb potential with a quantum geometric tensor of the TBG bands. The neutral excitations are (high-symmetry group) magnons, and their dispersion is analytically calculated in terms of the form factors of the active bands in TBG. The two-charge excitation energy and wave functions are also obtained, and a sufficient condition on the graphene eigenstates for obtaining a Cooper pair from Coulomb interactions is obtained. For the actual TBG bands at the first magic angle, we can analytically show that the Cooper pair binding energy is zero in all such projected Coulomb models, implying that either phonons and/or nonzero kinetic energy are needed for superconductivity. Since Vafek and Kang [Phys. Rev. Lett. 125, 257602 (2020)] showed that the kinetic energy bounds on the superexchange energy are less 10 -3 in Coulomb units, the phonon mechanism becomes then very likely. If nonetheless the superconductivity is due to kinetic terms which render the bands nonflat, one prediction of our theory is that the highest T c would not occur at the highest DOS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Neutrino Distributions in Phase Space for the Rotating Core-collapse Supernova Simulated with a Boltzmann-neutrino-radiation-hydrodynamics Code

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, in this work we simulate the collapse of a rotating core of the progenitor with a zero-age- main-sequence mass of 11.2 M ⊙ and a shellular rotation of $1\,\mathrm{rad}\,{{\rm{s}}}^{-1}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches ~20% in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum Monte Carlo methods for nuclear physics

Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Monte Carlo methods for nuclear physics

Quantum Monte Carlo methods have proved valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments, and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. The nuclear interactions and currents are reviewed along with a description of the continuum quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. A variety of results are presented, including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. Low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars are also described. Furthermore, a coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of the 2021 Milford, Utah earthquake swarm: Enhanced earthquake catalog and migration patterns

An earthquake swarm occurred in Spring 2021 in south-central Utah near the town of Milford. The University of Utah Seismograph Stations located 125 earthquakes between March 19 and May 10 with magnitudes ranging from 0.5 to 3.2. We implement a matched-filter technique in order to identify additional earthquakes that went undetected during the routine network location. The 125 network-located earthquakes are used as templates and are cross-correlated with continuous data for the dates Feb 17—June 10. This time period corresponds to approximately 1 month before the earthquake swarm began through 1 month after it ended. For the matched-filter analysis, we rely heavily on station FOR1, which is located within 5 km of most template events. Four other stations within 20–30 km of most template events provide a supplement to the closest station. The matched-filter implementation results in the detection of over 600 earthquakes in addition to the original 125 catalog events. This is one of the largest swarms ever recorded in Utah, and no previous large swarms have been recorded in this location. We use HypoDD to obtain relative double-difference locations of the catalog events. Both routine locations and HypoDD relocations of the catalog events suggest a fault dipping west beneath the mountain range, opposite of typical Basin and Range normal faults that dip beneath alluvial valleys. Moment tensors for the largest five events show normal faulting consistent with the west-dipping fault seen in the seismicity. Hydrothermal features in the area, including a geothermal power plant, suggest that fluids may be a contributing factor to the earthquake swarm triggering. We examine the role of fluids by exploring bounds on diffusion parameters and investigating spatial migration characteristics of the swarm seismicity. We conclude that this swarm is the result of heterogenous stress conditions in a prefractured region.

58 GEOSCIENCES↗

Amplitude/operator basis in chiral perturbation theory

We establish a systematic construction of the on-shell amplitude/operator basis for Chiral Perturbation Theory (ChPT) in D = 4 spacetime dimensions and with an arbitrary number of flavors N f . For kinematic factors, we employ spinor-helicity variables to construct the soft blocks, which are local amplitudes satisfying the Adler’s zero condition, as well as to take into account the reduction in the kinematic basis due to the Gram determinant, which arises at O(p 10 ) when the number of multiplicity N in an amplitude becomes large: N > D. For flavor factors, we include group-theoretic relations at small N f , N f ≤ N, which decreases the flavor basis. The result is obtained by adapting the Young tensor method of constructing the operator basis for generic effective field theories to the case of non-linearly realized symmetries. Working in the massless quark limit, we present purely mesonic operators for both even- and odd-parity at O(p 6 ) and O(p 8 ) for N = 6 and arbitrary N f , and establish a direct correspondence between the amplitude basis and the operator basis. Furthermore, the redundancy due to the Gram determinant is studied at O(p 10 ) for N = 6, 8 and 10.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tensor Force Manifestations in Ab Initio Study of the {sup 2}H(d,{gamma}){sup 4}He, {sup 2}H(d,p){sup 3}H, and {sup 2}H(d,n){sup 3}He Reactions

The {sup 2}H(d,p){sup 3}H, {sup 2}H(d,n){sup 3}He, and {sup 2}H(d,{gamma}){sup 4}He reactions are studied at low energies in a multichannel ab initio model that takes into account the distortions of the nuclei. The internal wave functions of these nuclei are given by the stochastic variational method with the AV8{sup '} realistic interaction and a phenomenological three-body force included to reproduce the two-body thresholds. The obtained astrophysical S factors are all in very good agreement with the experiment. The most important channels for both transfer and radiative capture are identified by comparing to calculations with an effective central force. They are all found to dominate thanks to the tensor force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleon Isovector Axial Form Factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1-flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron $Nπ$ states that make significant contributions. Our final results are $g_A$ = 1.292(53) stat (24) sys for the axial charge; $g_S$ = 1.085(50) stat (103) sys and $g_T$ = 0.991(21) stat (10) sys for the scalar and tensor charges; $\langle{r^2_A}\rangle$ = 0.439(56) sta t(34) sys fm 2 for the mean squared axial charge radius, $g^*_P$ = 9.03(47) stat (42) sys for the induced pseudoscalar charge; and $g_{πNN}$ = 14.14(81) stat (85) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor $G_A(Q^2)$ over the range 0 ≤ $Q^2$ ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of $G_A(Q^2)$ from the $ν$-deuterium scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleon isovector axial form factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1 -flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron N π states that make significant contributions. Our final results are g A = 1.292 ( 53 ) stat ( 24 ) sys for the axial charge; g S = 1.085 ( 50 ) stat ( 103 ) sys and g T = 0.991 ( 21 ) stat ( 10 ) sys for the scalar and tensor charges; ⟨ r A 2 ⟩ = 0.439 ( 56 ) stat ( 34 ) sys fm 2 for the mean squared axial charge radius, g P * = 9.03 ( 47 ) stat ( 42 ) sys for the induced pseudoscalar charge; and g π N N = 14.14 ( 81 ) stat ( 85 ) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor G A ( Q 2 ) over the range 0 ≤ Q 2 ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of G A ( Q 2 ) from the ν -deuterium scattering data. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Results from the Atacama B-mode Search (ABS) experiment

We report the Atacama B-mode Search is an experiment designed to measure the cosmic microwave background polarization at large angular scales (0ℓ>4). It observes at 145 GHz from a site at 5,190 m elevation in northern Chile. The noise equivalent polarization temperature, or NEQ, is 41 μK√s. One of the unique features of ABS is its use of a rapidly rotating ambient-temperature half-wave plate (HWP) {as the first optical element}. {The HWP spins} at 2.55 Hz to modulate the incident polarized signal at frequencies above where instrument white noise dominates over atmospheric fluctuations and other sources of low-frequency noise. We report here on the analysis of data from a 2,400 deg 2 region of sky. We perform a blind analysis to reduce potential bias. After unblinding, we find agreement with the Planck TE and EE measurements on the same region of sky, {with a derived calibration factor of 00.89 ± 0.1}. We marginally detect polarized dust emission {(at 3.2 σ for EE and 2.2 σ for BB)} and give an upper limit on the tensor-to-scalar ratio of r<2.3 (95% confidence level) with the equivalent of 100 on-sky days of observation. Finally, we also present a new measurement of the polarization of Tau A and introduce new methods for calibration and data analysis associated with HWP-based observations.

79 ASTRONOMY AND ASTROPHYSICS↗