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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluation of PbTe and SnTe as Ohmic Contact Layers in CdTe Solar Cell Devices

For solar cells based on CdTe, the choice of a suitable back contact material is limited by CdTe's deep work function. Here, we explored p-type PbTe and SnTe as ohmic contacts to CdTe. These contact layers were grown on single crystal CdTe substrates by molecular beam epitaxy, and the valence band offset between film and substrate was measured using X-ray photoemission spectroscopy. Polycrystalline device structures were also grown by sublimation to assess performance improvements. Doping was achieved in PbTe by thallium incorporation. Only the highest Tl doped PbTe resulted in a desirable band alignment with the CdTe, forming an electron reflector and no hole barrier. Time-resolved photoluminescence measurements also revealed significant photocarrier lifetime improvements for only the highest doped PbTe. Consequentially, devices incorporating the highest doped PbTe layers showed increased power conversion efficiency, primarily from increased fill factor. Doping of the PbTe was measured via Hall effect with variable magnetic field, which was required due to the formation of an n-type parasitic interface layer. To properly interpret the variable field Hall measurements, we derived an expression for the magnetic field-dependent conductivity tensor of an L-valley semiconductor.

anisotropic band valley↗
DOE OSTI2018

High-resolution sub-Doppler infrared spectroscopy of atmospherically relevant Criegee precursor CH 2 I radicals: CH 2 stretch vibrations and “charge-sloshing” dynamics

The combination of a pulsed supersonic slit-discharge source and single-mode difference frequency direct absorption infrared spectroscopy permit first high resolution infrared study of the iodomethyl (CH 2 I) radical, with the CH 2 I radical species generated in a slit jet Ne/He discharge and cooled to 16 K in the supersonic expansion. Dual laser beam detection and collisional collimation in the slit expansion yield sub-Doppler linewidths (60 MHz), an absolute frequency calibration of 13 MHz, and absorbance sensitivities within a factor of two of the shot-noise limit. Fully rovibrationally resolved direct absorption spectra of the CH 2 symmetric stretch mode (ν 2 ) are obtained and fitted to a Watson asymmetric top Hamiltonian with electron spin-rotation coupling, providing precision rotational constants and spin-rotation tensor elements for the vibrationally excited state. Assessment of the asymmetric top rotational constants confirms a vibrationally averaged planar geometry in both the ground- and first-excited vibrational levels. Sub-Doppler resolution permits additional nuclear spin hyperfine structures to be observed, with splittings in excellent agreement with microwave measurements on the ground state. Spectroscopic data on CH 2 I facilitate systematic comparison with previous studies of halogen-substituted methyl radicals, with the periodic trends strongly correlated with the electronegativity of the halogen atom. Surprisingly, we do not observe any asymmetric CH 2 stretch transitions, despite S/N ≈ 25:1 on strongest lines in the corresponding symmetric CH 2 stretch manifold. This dramatic reversal of the more typical 3:1 antisymmetric/symmetric CH 2 stretch intensity ratio signals a vibrational transition moment poorly described by simple “bond-dipole” models. Instead, the data suggest that this anomalous intensity ratio arises from “charge sloshing” dynamics in the highly polar carbon-iodine bond, as supported by ab initio electron differential density plots and indeed consistent with observations in other halomethyl radicals and protonated cluster ions.

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

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

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

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Embedding DBI inflation in scalar-tensor theory

The Dirac-Born-Infeld (DBI) action has been widely studied as an interesting example of a model of k-inflation in which the sound speed of the cosmological perturbations differs from unity. In this article we consider a scalar-tensor theory in which the matter component is a field with a DBI action. Transforming to the Einstein frame, we explore the effect of the resulting coupling on the background dynamics of the fields and the first-order perturbations. We find that the coupling forces the scalar field into the minimum of its effective potential. While the additional scalar field contributes significantly to the energy density during inflation, the dynamics are determined by the DBI field, which has the interesting effect of increasing the number of efolds of inflation and decreasing the boost factor of the DBI field. Focusing on this case, we show, with the benefit of numerical examples, that the power spectrum of the primordial perturbations is determined by the behaviour of the perturbations of the modified DBI field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

M5-brane and D-brane scattering amplitudes

We present tree-level $\mathcal{n}$-particle on-shell scattering amplitudes of various brane theories with 16 conserved supercharges. These include the world-volume theory of a probe D3-brane or D5-brane in 10D Minkowski spacetime as well as a probe M5-brane in 11D Minkowski spacetime, which describes self interactions of an abelian tensor supermultiplet with 6D (2, 0) supersymmetry. Twistor-string-like formulas are proposed for tree-level scattering amplitudes of all multiplicities for each of these theories. The R symmetry of the D3-brane theory is shown to be SU(4) × U(1), and the U(1) factor implies that its amplitudes are helicity conserving. Each of 6D theories (D5-brane and M5-brane) reduces to the D3-brane theory by dimensional reduction. Conclusively, as special cases of the general M5-brane amplitudes, we present compact formulas for examples involving only the self-dual B field with $\mathcal{n}$ = 4, 6, 8.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Tensor-structured coupled cluster theory

We derive and implement a new way of solving coupled cluster equations with lower computational scaling. Our method is based on the decomposition of both amplitudes and two electron integrals, using a combination of tensor hypercontraction and canonical polyadic decomposition. While the original theory scales as O ( N 6 ) with respect to the number of basis functions, we demonstrate numerically that we achieve sub-millihartree difference from the original theory with O ( N 4 ) scaling. This is accomplished by solving directly for the factors that decompose the cluster operator. Here, the proposed scheme is quite general and can be easily extended to other many-body methods.

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

Faster Johnson–Lindenstrauss transforms via Kronecker products

The Kronecker product is an important matrix operation with a wide range of applications in signal processing, graph theory, quantum computing and deep learning. In this work, we introduce a generalization of the fast Johnson–Lindenstrauss projection for embedding vectors with Kronecker product structure, the Kronecker fast Johnson–Lindenstrauss transform (KFJLT). The KFJLT reduces the embedding cost by an exponential factor of the standard fast Johnson–Lindenstrauss transform’s cost when applied to vectors with Kronecker structure, by avoiding explicitly forming the full Kronecker products. Here, we prove that this computational gain comes with only a small price in embedding power: consider a finite set of $p$ points in a tensor product of $d$ constituent Euclidean spaces $\bigotimes _{k=d}^{1}{\mathbb{R}}^{n_k}$, and let $N = \prod _{k=1}^{d}n_k$. With high probability, a random KFJLT matrix of dimension $m \times N$ embeds the set of points up to multiplicative distortion $(1\pm \varepsilon )$ provided $m \gtrsim \varepsilon ^{-2} \, \log ^{2d - 1} (p) \, \log N$. We conclude by describing a direct application of the KFJLT to the efficient solution of large-scale Kronecker-structured least squares problems for fitting the CP tensor decomposition.

Kronecker structure↗
DOE OSTI2015

Similarity and differences between the radion and Higgs boson production and decay processes involving off-shell fermions

The radion is a scalar particle that occurs in brane world models and interacts with the trace of the energy–momentum tensor of the Standard Model (SM). The radion–SM fermion interaction Lagrangian differs from the Higgs boson–fermion interaction Lagrangian for off-shell fermions. It is shown that all additional, as compared to the Higgs boson, contributions to the amplitudes of radion production and decay processes involving off-shell fermions are canceled out for both massless and massive fermions. Thus, additional terms in the interaction Lagrangian do not change properties of these processes for the radion and the Higgs boson, except for the general normalization factors. This similarity is a consequence of gauge invariance for the processes with production of gauge bosons. When an additional scalar particle is produced, there are no apparent reasons for the above cancellation, as confirmed, for example, by the process with production of two scalar particles, which features an additional contribution of the radion in comparison with the Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Determining Gyration Tensor of Orienting Macromolecules through Their Scattering Signature

A method is presented for quantitatively evaluating the shape and size of deformed particles in dispersion from their two-dimensional anisotropic spectra by small-angle scattering. By means of real spherical harmonic expansion, we derive analytical expressions of the gyration tensor R in terms of experimentally measured anisotropic scattering functions, yielding a tensorial extension of the Guinier law. We demonstrate the usefulness of this approach by a model study of an affinely deformed Gaussian chain. We further show that radius of gyration R g is the source term of intraparticle structure factor at the mean-field limit, and from this perspective, we address the connection between R and conformation asphericity. As a result, the developed method not only facilitates quantitative scattering studies of deforming materials, but also provides insightful information regarding their deformation behavior at the molecular level based on the symmetric properties of real spherical harmonics.

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

Displacement Fields and Self-Energies of Circular and Polygonal Dislocation Loops in Homogeneous and Layered Anisotropic Solids

There are large classes of materials problems that involve the solutions of stress, displacement, and strain energy of dislocation loops in elastically anisotropic solids, including increasingly detailed investigations of the generation and evolution of irradiation induced defect clusters ranging in sizes from the micro- to meso-scopic length scales. Based on a two-dimensional Fourier transform and Stroh formalism that are ideal for homogeneous and layered anisotropic solids, we have developed robust and computationally efficient methods to calculate the displacement fields for circular and polygonal dislocation loops. Using the homogeneous nature of the Green tensor of order -1, we have shown that the displacement and stress fields of dislocation loops can be obtained by numerical quadrature of a line integral. In addition, it is shown that the sextuple integrals associated with the strain energy of loops can be represented by the product of a pre-factor containing elastic anisotropy effects and a universal term that is singular and equal to that for elastic isotropic case. Furthermore, we have found that the self-energy pre-factor of prismatic loops is identical to the effective modulus of normal contact, and the pre-factor of shear loops differs from the effective indentation modulus in shear by only a few percent. These results provide a convenient method for examining dislocation reaction energetic and efficient procedures for numerical computation of local displacements and stresses of dislocation loops, both of which play integral roles in quantitative defect analyses within combined experimental–theoretical investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Dissipative axial inflation

We analyze in detail the background cosmological evolution of a scalar field coupled to a massless abelian gauge field through an axial term φ/ f {sub γ} F ∼ F , such as in the case of an axion. Gauge fields in this case are known to experience tachyonic growth and therefore can backreact on the background as an effective dissipation into radiation energy density ρ{sub R}, which can lead to inflation without the need of a flat potential. We analyze the system, for momenta k smaller than the cutoff f {sub γ}, including the backreaction numerically. We consider the evolution from a given static initial condition and explicitly show that, if f {sub γ} is smaller than the field excursion φ{sub 0} by about a factor of at least O (20), there is a friction effect which turns on before the field can fall down and which can then lead to a very long stage of inflation with a generic potential. In addition we find superimposed oscillations, which would get imprinted on any kind of perturbations, scalars and tensors. Such oscillations have a period of 4–5 efolds and an amplitude which is typically less than a few percent and decreases linearly with f {sub γ}. We also stress that the curvature perturbation on uniform density slices should be sensitive to slow-roll parameters related to ρ{sub R} rather than φ-dot {sup 2}/2 and we discuss the existence of friction terms acting on the perturbations, although we postpone a calculation of the power spectrum and of non-gaussianity to future work and we simply define and compute suitable slow roll parameters. Finally we stress that this scenario may be realized in the axion case, if the coupling 1/ f {sub γ} to U(1) (photons) is much larger than the coupling 1/ f {sub G} to non-abelian gauge fields (gluons), since the latter sets the range of the potential and therefore the maximal allowed φ{sub 0∼} f {sub G}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Photonic band structures of one-dimensional photonic crystals doped with plasma

The photonic band structures (PBSs) of oblique incidence propagation in one-dimensional plasma-doped photonic crystals (PCs) are investigated carefully. When the lattice constant of plasma-doped PCs is less than the incident wavelength, the PC becomes anisotropic. Therefore, the dielectric constant of PC is converted into a complex tensor dielectric constant. This determines the PBSs of PCs. In the present paper, one-dimensional PCs are taken as an example to study both normal and absorption PBSs. Using both the effective medium approximation and the transfer matrix method, we can derive the dispersion relation for PCs. The dependence of the plasma filling factor on the effective dielectric constant and PBSs is calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fermionic systems for quantum information people

The operator algebra of fermionic modes is isomorphic to that of qubits, the difference between them is twofold: the embedding of subalgebras corresponding to mode subsets and multiqubit subsystems on the one hand, and the parity superselection in the fermionic case on the other. We discuss these two fundamental differences extensively, and illustrate these through the Jordan–Wigner representation in a coherent, self-contained, pedagogical way, from the point of view of quantum information theory. Our perspective leads us to develop useful new tools for the treatment of fermionic systems, such as the fermionic (quasi-)tensor product, fermionic canonical embedding, fermionic partial trace, fermionic products of maps and fermionic embeddings of maps. We formulate these by direct, easily applicable formulas, without mode permutations, for arbitrary partitionings of the modes. It is also shown that fermionic reduced states can be calculated by the fermionic partial trace, containing the proper phase factors. We also consider variants of the notions of fermionic mode correlation and entanglement, which can be endowed with the usual, local operation based motivation, if the parity superselection rule is imposed. We also elucidate some other fundamental points, related to joint map extensions, which make the parity superselection inevitable in the description of fermionic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2023

BitGNN: Unlocking the Performance Potential of Binary Graph Neural Networks on GPUs

Graph Neural Networks (GNNs) have shown compelling results in many graph-based learning tasks. They are, however, time-consuming. Recent work has shown a promising direction in improving GNN speed and shrinking the size — network binarization, which binarizes network values and operations. Prior work, however, mainly focused on algorithm designs, leaving it open on how to fully materialize the performance potential. This work fills the gap by proposing techniques to best map binary GNNs and their computations to fit the nature of bit manipulations, optimizations and algorithms to maximize BSpMM kernel efficiency, and solutions to other factors influencing the end-to-end time on GPUs. Results on real-world graphs show that the proposed techniques outperform state of-the-art binary GNN implementations by 21-67× with little accuracy loss.

Chen, Jou-An↗