Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tensor”

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

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

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

Open-Shell Tensor Hypercontraction

The extension of least-squares tensor hypercontracted second- and third-order Møller–Plesset perturbation theory (LS-THC-MP2 and LS-THC-MP3) to open-shell systems is an important development due to the scaling reduction afforded by THC and the ubiquity of molecular ions, radicals, and other open-shell reactive species. The complexity of wavefunction-based quantum-chemical methods such as Møller–Plesset and coupled cluster theory is reflected in the steep scaling of the computational costs with the molecular size. The least-squares tensor hypercontraction (LS-THC) method is an efficient, single-step factorization for the two-electron integral tensor but can also be used to factorize the double excitation amplitudes, leading to significant scaling reduction. Herein we extend this promising method to open-shell variants of LS-THC-MP2 and -MP3 by using diagrammatic techniques and explicit spin summation. The accuracy of the resulting methods for open-shell species is benchmarked on standard test systems such as regular alkanes as well as realistic systems involving bond breaking, radical stabilization, and other effects. We find that open-shell LS-THC-MPn methods exhibit errors highly comparable to those produced by closed-shell LS-THC-MPn and are highly insensitive to particular chemical interactions, geometries, or even moderate spin contamination.

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

Holographic codes from hyperinvariant tensor networks

Holographic quantum-error correcting codes are models of bulk/boundary dualities such as the anti-de Sitter/conformal field theory (AdS/CFT) correspondence, where a higher-dimensional bulk geometry is associated with the code’s logical degrees of freedom. Previous discrete holographic codes based on tensor networks have reproduced the general code properties expected from continuum AdS/CFT, such as complementary recovery. However, the boundary states of such tensor networks typically do not exhibit the expected correlation functions of CFT boundary states. In this work, we show that a new class of exact holographic codes, extending the previously proposed hyperinvariant tensor networks into quantum codes, produce the correct boundary correlation functions. This approach yields a dictionary between logical states in the bulk and the critical renormalization group flow of boundary states. Furthermore, these codes exhibit a state-dependent breakdown of complementary recovery as expected from AdS/CFT under small quantum gravity corrections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Determination of the Raman polarizability tensor in the optically anisotropic crystal potassium dihydrogen phosphate and its deuterated analog

Abstract The Raman tensor of the dominant A 1 modes of the nonlinear optical crystalline material potassium dihydrogen phosphate and its 70% deuterated analog have been ascertained. Challenges in determining the A 1 mode tensor element values based on previous reports have been resolved using a specially designed experimental setup that makes use of spherical crystal samples. This novel experimental design enabled the determination of measurement artifacts, including polarization rotation of the pump and/or scattered light propagating through the sample and the contribution of additional overlapping phonon modes, which have hindered previous efforts. Results confirmed that the polarization tensor is diagonal, and matrix elements were determined with high accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning of 27Al NMR electric field gradient tensors for crystalline structures from DFT

NMR crystallography has emerged as a promising technique for the determination and refinement of atomic coordinates in crystal structures. The crystal structure of compounds containing quadrupolar nuclei, such as 27Al, can be improved by directly comparing solid-state NMR measurements to DFT computations of the electric field gradient (EFG) tensor. The non-negligible computational cost of these first-principles calculations limits the applicability of this method to all but the most well-defined structures. We developed a fast, low-cost machine learning model to predict EFG parameters based on local structural motifs and elemental parameters. We computed 8081 EFG tensors from 1681 27Al crystalline solids using DFT and benchmarked them against 105 experimentally measured 27Al sites. Surprisingly, simple local geometric features dominate the predictive performance of the resulting random-forest model, yielding an R2 value of 0.98 and an RMSE of 0.61 MHz for CQ, the quadrupolar coupling constant. This model accuracy should enable pre-refining future structural assignments before finally validating with first-principles calculations. Such a catalogue of 27Al NMR tensors can serve as a tool for researchers assigning complex NMR spectra influenced by the nuclear electric quadrupole interaction.

Sun, He↗
DOE OSTI2014

Eckart frame vibration-rotation Hamiltonians: Contravariant metric tensor

Eckart frame is a unique embedding in the theory of molecular vibrations and rotations. It is defined by the condition that the Coriolis coupling of the reference structure of the molecule is zero for every choice of the shape coordinates. It is far from trivial to set up Eckart kinetic energy operators (KEOs), when the shape of the molecule is described by curvilinear coordinates. In order to obtain the KEO, one needs to set up the corresponding contravariant metric tensor. Here, I derive explicitly the Eckart frame rotational measuring vectors. Their inner products with themselves give the rotational elements, and their inner products with the vibrational measuring vectors (which, in the absence of constraints, are the mass-weighted gradients of the shape coordinates) give the Coriolis elements of the contravariant metric tensor. The vibrational elements are given as the inner products of the vibrational measuring vectors with themselves, and these elements do not depend on the choice of the body-frame. The present approach has the advantage that it does not depend on any particular choice of the shape coordinates, but it can be used in conjunction with all shape coordinates. Furthermore, it does not involve evaluation of covariant metric tensors, chain rules of derivation, or numerical differentiation, and it can be easily modified if there are constraints on the shape of the molecule. Both the planar and non-planar reference structures are accounted for. The present method is particular suitable for numerical work. Its computational implementation is outlined in an example, where I discuss how to evaluate vibration-rotation energies and eigenfunctions of a general N-atomic molecule, the shape of which is described by a set of local polyspherical coordinates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Development of a vector-tensor system to measure the absolute magnetic flux density and its gradient in magnetically shielded rooms

Several experiments in fundamental physics demand an environment of very low, homogeneous, and stable magnetic fields. For the magnetic characterization of such environments, we present a portable SQUID system that measures the absolute magnetic flux density vector and the gradient tensor. This vector-tensor system contains 13 integrated low-critical temperature (LTc) superconducting quantum interference devices (SQUIDs) inside a small cylindrical liquid helium Dewar with a height of 31 cm and 37 cm in diameter. The achievable resolution depends on the flux density of the field under investigation and its temporal drift. Inside a seven-layer mu-metal shield, an accuracy better than ±23 pT for the components of the static magnetic field vector and ±2 pT/cm for each of the nine components of the gradient tensor is reached by using the shifting method.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

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↗
DOE OSTI2012

POWER ANISOTROPY IN THE MAGNETIC FIELD POWER SPECTRAL TENSOR OF SOLAR WIND TURBULENCE

We observe the anisotropy of the power spectral tensor of magnetic field fluctuations in the fast solar wind for the first time. In heliocentric RTN coordinates, the power in each element of the tensor has a unique dependence on the angle between the magnetic field and velocity of the solar wind ({theta}{sub B}) and the angle of the vector in the plane perpendicular to the velocity ({phi}{sub B}). We derive the geometrical effect of the high speed flow of the solar wind past the spacecraft on the power spectrum in the frame of the plasma P(k) to arrive at the observed power spectrum P(f, {theta}{sub B}, {phi}{sub B}) based on a scalar field description of turbulence theory. This allows us to predict the variation in the {phi}{sub B} direction and compare it to the data. We then transform the observations from RTN coordinates to magnetic-field-aligned coordinates. The observed reduced power spectral tensor matches the theoretical predictions we derive in both RTN and field-aligned coordinates, which means that the local magnetic field we calculate with wavelet envelope functions is an accurate representation of the physical axis of symmetry for the turbulence and implies that on average the turbulence is axisymmetric. We also show that we can separate the dominant toroidal component of the turbulence from the smaller but significant poloidal component and that these have different power anisotropy. We also conclude that the magnetic helicity is anisotropic and mostly two dimensional, arising from wavevectors largely confined to the plane perpendicular to B .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Canonical single field slow-roll inflation with a non-monotonic tensor-to-scalar ratio

We take a pragmatic, model independent approach to single field slow-roll canonical inflation by imposing conditions, not on the potential, but on the slow-roll parameter ε(φ) and its derivatives ε'(φ) and ε''(φ), thereby extracting general conditions on the tensor-to-scalar ratio r and the running n {sub sk} at φ {sub H} where the perturbations are produced, some 50–60 e -folds before the end of inflation. We find quite generally that for models where ε(φ) develops a maximum, a relatively large r is most likely accompanied by a positive running while a negligible tensor-to-scalar ratio implies negative running. The definitive answer, however, is given in terms of the slow-roll parameter ξ{sub 2}(φ). To accommodate a large tensor-to-scalar ratio that meets the limiting values allowed by the Planck data, we study a non-monotonic ε(φ) decreasing during most part of inflation. Since at φ {sub H} the slow-roll parameter ε(φ) is increasing, we thus require that ε(φ) develops a maximum for φ > φ {sub H} after which ε(φ) decrease to small values where most e -folds are produced. The end of inflation might occur trough a hybrid mechanism and a small field excursion Δφ {sub e} ≡ |φ {sub H} −φ {sub e} | is obtained with a sufficiently thin profile for ε(φ) which, however, should not conflict with the second slow-roll parameter η(φ). As a consequence of this analysis we find bounds for Δφ {sub e} , r {sub H} and for the scalar spectral index n {sub sH} . Finally we provide examples where these considerations are explicitly realised.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Tensor Minkowski Functionals: first application to the CMB

Tensor Minkowski Functionals (TMFs) are tensor generalizations of the usual Minkowski Functionals which are scalar quantities. We introduce them here for use in cosmological analysis, in particular to analyze the Cosmic Microwave Background (CMB) radiation. They encapsulate information about the shapes of structures and the orientation of distributions of structures. We focus on one of the TMFs, namely W {sub 2}{sup 1,1}, which is the (1,1) rank tensor generalization of the genus. The ratio of the eigenvalues of the average of W {sub 2}{sup 1,1} over all structures, α, encodes the net orientation of the structures; and the average of the ratios of the eigenvalues of W {sub 2}{sup 1,1} for each structure, β, encodes the net intrinsic anisotropy of the structures. We have developed a code that computes W {sub 2}{sup 1,1}, and from it α and β, for a set of structures on the 2-dimensional Euclidean plane. We use it to compute α and β as functions of chosen threshold levels for simulated Gaussian and isotropic CMB temperature and E mode fields. We obtain the value of α to be one for both temperature and E mode, which means that we recover the statistical isotropy of density fluctuations that we input in the simulations. We find that the standard ΛCDM model predicts a charateristic shape of β for temperature and E mode as a function of the threshold, and the average over thresholds is β∼ 0.62 for temperature and β∼ 0.63 for E mode. Accurate measurements of α and β can be used to test the standard model of cosmology and to search for deviations from it. For this purpose we compute α and β for temperature and E mode data of various data sets from PLANCK mission. We compare the values measured from observed data with those obtained from simulations to which instrument beam and noise characteristics of the 44GHz frequency channel have been added (which are provided as part of the PLANCK data release). We find very good agreement of β and α between all PLANCK temperature data sets with ΛCDM expectations. E mode data show good agreement for β but α for all data sets deviate from ΛCDM predictions higher than 3−σ. It is most likely that the deviations are probing the anisotropy of the noise field and beam characteristics of the detector rather than the true E mode signal since for 44GHz the signal-to-noise ratio is well below one. This will be further investigated after the full PLANCK data becomes publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Potential for a Tensor Asymmetry Azz Measurement in the x > 1 Region at Jefferson Lab

The tensor asymmetry A(zz) in the quasi-elastic region through the tensor polarized D(e, e')X channel is sensitive to the nucleon-nucleon potential. Previous measurements of A(zz) have been used to extract b(1) in the DIS region and T(20) in the elastic region. In the quasielastic region, A(zz) can be used to compare light cone calculations with variation nucleon- nucleon calculations, and is an important quantity to determine for understanding tensor effects, such as the dominance of pn correlations in nuclei. In the quasi-elastic region, A(zz) was first calculated in 1988 by Frankfurt and Strikman using the Hamada-Johnstone and Reid soft-core wave functions [1]. Recent calculations by M. Sargsian revisit A(zz) in the x > 1 range using virtual-nucleon and light-cone methods, which differ by up to a factor of two [2]. Discussed in these proceedings, a study has been completed that determines the feasibility of measuring A(zz) in the quasi-elastic x > 1 region at Jefferson Lab's Hall C.

Long, Elena↗
DOE OSTI2019

Towards quantum machine learning with tensor networks

Machine learning is a promising application of quantum computing, but challenges remain for implementation today because near-term devices have a limited number of physical qubits and high error rates. Motivated by the usefulness of tensor networks for machine learning in the classical context, we propose quantum computing approaches to both discriminative and generative learning, with circuits based on tree and matrix product state tensor networks, that could already have benefits with such near-term devices. The result is a unified framework in which classical and quantum computing can benefit from the same theoretical and algorithmic developments, and the same model can be trained classically then transferred to the quantum setting for additional optimization. Tensor network circuits can also provide qubit-efficient schemes in which, depending on the architecture, the number of physical qubits required scales only logarithmically with, or independently of the input or output data sizes. Here, we demonstrate our proposals with numerical experiments, training a discriminative model to perform handwriting recognition using a hybrid quantum-classical optimization procedure that could be carried out on quantum hardware today, and testing the noise resilience of the trained model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coarse-grain cluster analysis of tensors with application to climate biome identification

A tensor provides a concise way to codify the interdependence of complex data. Treating a tensor as a d-way array, each entry records the interaction between the different indices. Clustering provides a way to parse the complexity of the data into more readily understandable information. Clustering methods are heavily dependent on the algorithm of choice, as well as the chosen hyperparameters of the algorithm. However, their sensitivity to data scales is largely unknown. In this work, we apply the discrete wavelet transform to analyze the effects of coarse-graining on clustering tensor data. We are particularly interested in understanding how scale affects clustering of the Earth's climate system. The discrete wavelet transform allows classification of the Earth's climate across a multitude of spatial-temporal scales. The discrete wavelet transform is used to produce an ensemble of classification estimates, as opposed to a single classification. Each element of the ensemble is a clustering at a different spatial-temporal scale. Information theoretic approaches are used to identify important scale lengths in clustering the L15 Climate Datset. We also discover a sub-collection of the ensemble that spans the majority of the variance observed, allowing for efficient consensus clustering techniques that can be used to identify climate biomes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Seismic moment tensor classification using elliptical distribution functions on the hypersphere

Discrimination of underground explosions from naturally occurring earthquakes and other anthropogenic sources is one of the fundamental challenges of nuclear explosion monitoring. In an operational setting, the number of events that can be thoroughly investigated by analysts is limited by available resources. The capability to rapidly screen out events that can be robustly identified as not being explosions is, therefore, of great potential benefit. Nevertheless, possible mis-classification of explosions as earthquakes currently limits the use of screening methods for verification of test-ban treaties. Moment tensors provide a physics-based classification tool for the characterization of different seismic sources and have enabled the advent of new techniques for discriminating between earthquakes and explosions. Following normalization and projection of their six-degree vectors onto the hypersphere, existing screening approaches use spherically symmetric metrics to determine whether any new moment tensor may have been an explosion. Here, we show that populations of moment tensors for both earthquakes and explosions are anisotropically distributed on the hypersphere. Distributions possessing elliptical symmetry, such as the scaled von Mises–Fisher distribution, therefore provide a better description of these populations than the existing spherically symmetric models. We describe a method that uses these elliptical distributions in combination with a Bayesian classifier to achieve successful classification rates of 99 per cent for explosions and 98 per cent for earthquakes using existing catalogues of events from the western United States. The 1983 May 5 Crowdie underground nuclear test and 2018 July 20 DAG-1 deep-borehole chemical explosion are the only two explosions out of 140 that are incorrectly classified. Application of the method to the 2006–2017 nuclear tests in the Democratic People’s Republic of Korea yields 100 per cent identification rates and we provide a simple routine MTid for general usage. The approach provides a means to rapidly assess the likelihood of an event being an explosion and can be built into monitoring workflows that rely on simultaneously assessing multiple different discrimination metrics.

58 GEOSCIENCES↗
DOE OSTI2018

Analysis of seismic moment tensors and in situ stress during Enhanced Geothermal System development at The Geysers geothermal field, California

Seismicity in the vicinity of an Enhanced Geothermal System (EGS) demonstration project in the Northwest Geysers at the Prati 32 (P-32) injection well and the Prati State 31 (PS-31) production well is investigated to determine earthquake focal mechanisms, moment magnitudes and in situ stress of events that occur before and during reservoir stimulation starting 2011 October 6. We present a catalogue of 167 waveform-based seismic moment tensors ranging in moment magnitude from 0.6 to 3.9. The moment tensor catalogue is subsequently used to invert for the stress tensor and to investigate possible temporal stress changes resulting from fluid injection. We find an approximate 15-degree counter-clockwise rotation of the least compressive stress σ3 during injection. More remarkable is a change in orientation of the maximum compressive stress σ1 from subhorizontal to vertical when injection operations temporarily cease. The orientation of σ1 returns to subhorizontal when injection operations resume. It is found that there was a systematic reduction in the stress shape factor, R, as the injected water volume increases, indicating an evolution towards a more transtensional stress state.

58 GEOSCIENCES↗
DOE OSTI2026

Elastic and resonance structures of the nucleon from the hadronic tensor in lattice QCD: Implications for neutrino-nucleon scattering and hadron physics

We compute the Euclidean hadronic tensor from charge density operators and extract elastic and resonance structures by employing exponential fits to the four-point function correlator, as well as a Bayesian reconstruction inverse algorithm to obtain the corresponding spectral density for qualitative comparison. We present the determination of the nucleon’s Sachs electric form factor using the hadronic tensor formalism and verify that it is consistent with that from the conventional three-point function calculation. Beyond the elastic peak, we observe a structure located approximately 0.5–0.7 GeV above the nucleon mass in the Bayesian reconstruction. This structure is interpreted as a mixture of the Roper resonance [𝑁⁡(1440)], and states with both positive and negative parities in this mass region, as well as multihadron states. Assuming the observed structure is dominated by 𝐽 𝑃 = 1/2 ± states, we extract the transition electric form factor 𝐺$^*_𝐸$⁡(𝑄 2 ) and the corresponding longitudinal helicity amplitude 𝑆 1/2 ⁡(𝑄 2 ), and compare them with those determined from the CLAS experimental data of nucleon-to-Roper transition. Although fitting to the four-point correlation function or using the inverse algorithm does not resolve individual resonances, it nevertheless enables the determination of total inclusive lepton–nucleon scattering cross sections in appropriate energy bins. This lattice QCD calculation presents the first major step toward studying the inclusive 𝑁 → 𝑋 contributions within the hadronic tensor formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Symmetry analysis of tensors in the honeycomb lattice of edge-sharing octahedra

In this work, we obtain the most general forms of rank-2 and rank-3 tensors allowed by the crystal symmetries of the honeycomb lattice of edge-sharing octahedra for crystals belonging to different crystallographic point groups, including the monoclinic point group 2/$\textit{m}$ and the trigonal (or rhombohedral) point group $\bar{3}$. Our results are relevant for two-dimensional materials, such as α-RuCl 3 , CrI 3 , and the honeycomb iridates. We focus on the magnetic-field-dependent thermal conductivity tensor $κ_{ij}$ (H), which describes a system’s longitudinal and thermal Hall responses, for the cases when the magnetic field is applied along high-symmetry directions, perpendicular to the plane and in the plane. We highlight some unexpected results, such as the equality of fully longitudinal components to partially transverse components in rank-3 tensors for

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗