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

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2022

An efficient method to propagate model uncertainty when inverting seismic data for time domain seismic moment tensors

SUMMARY We present a computationally efficient method to approximately propagate uncertainty when linearly inverting seismic data for point source, time variable moment tensor components. The method is based on the assumption that the data residual, given by the difference between the observed seismic data and the data predicated by a linear inversion, contains the effects of both data and model uncertainty. Our method uses a distribution of data residuals, added directly to the data, in a pseudo-Monte Carlo scheme. Using the assumption that the data residual is a stochastic process, we use the well-known Karhunen–Loève (KL) theorem to construct a distribution of data residuals, where the required basis functions are constructed using Fourier series. The Fourier series are scaled by a product of a random variable and the real-valued spectral amplitudes of the original data residual’s spectrum. Thus, the Fourier series and spectral amplitudes are eigenfunction-eigenvalue pairs used in the KL-based construction of data residual distribution. Using tests with synthetic data, we show that our method compares closely with a Finite Difference Monte Carlo (FDMC) method that we presented previously. More importantly, the method presented here is computationally several orders of magnitude faster than our previous FDMC method, and requires no a priori assumptions of model and/or data uncertainty.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2023

Moment tensor event identification for collapses

SUMMARY We introduce a seismic identification method for collapse events using moment tensors (MTs). We start by computing full (six-element) MT solutions for 43 identified collapse events from around the world, and statistically characterizing the population on the MT hypersphere. We then test a large data set of over 1000 full MTs for the western U.S. against the distribution of collapses using a MT-based identification method similarly as used for testing explosions. Known collapses and explosions are readily identified, along with other anomalous events in the Geysers and central California coast. Misidentification rates are determined for various screening angles with optimal misidentification rates between earthquakes and collapses on the order of 3 per cent. The method is demonstrated to be very effective at identifying non-earthquake sources with a 97–98 per cent accuracy. It is likely to be transportable to other regions, and can be used for event identification anywhere full MT solutions are routinely calculated.

58 GEOSCIENCES↗
DOE OSTI2025

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗
DOE OSTI2022

Bias on tensor-to-scalar ratio inference with estimated covariance matrices

ABSTRACT We investigate simulation-based bandpower covariance matrices commonly used in cosmological parameter inferences such as the estimation of the tensor-to-scalar ratio r. We find that upper limits on r can be biased low by tens of per cent. The underestimation of the upper limit is most severe when the number of simulation realizations is similar to the number of observables. Convergence of the covariance-matrix estimation can require a number of simulations an order of magnitude larger than the number of observables, which could mean $\mathcal {O}(10\ 000)$ simulations. This is found to be caused by an additional scatter in the posterior probability of r due to Monte Carlo noise in the estimated bandpower covariance matrix, in particular, by spurious non-zero off-diagonal elements. We show that matrix conditioning can be a viable mitigation strategy in the case that legitimate covariance assumptions can be made.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mixed-precision iterative refinement using tensor cores on GPUs to accelerate solution of linear systems

Double-precision floating-point arithmetic (FP64) has been the de facto standard for engineering and scientific simulations for several decades. Problem complexity and the sheer volume of data coming from various instruments and sensors motivate researchers to mix and match various approaches to optimize compute resources, including different levels of floating-point precision. In recent years, machine learning has motivated hardware support for half-precision floating-point arithmetic. A primary challenge in high-performance computing is to leverage reduced-precision and mixed-precision hardware. We show how the FP16/FP32 Tensor Cores on NVIDIA GPUs can be exploited to accelerate the solution of linear systems of equations Ax = b without sacrificing numerical stability. The techniques we employ include multiprecision LU factorization, the preconditioned generalized minimal residual algorithm (GMRES), and scaling and auto-adaptive rounding to avoid overflow. We also show how to efficiently handle systems with multiple right-hand sides. On the NVIDIA Quadro GV100 (Volta) GPU, we achieve a 4×-5× performance increase and 5× better energy efficiency versus the standard FP64 implementation while maintaining an FP64 level of numerical stability.

GMRES↗
DOE OSTI2025

Emergent Nonthermal Fluid from Jets in the Massive Schwinger Model Using Tensor Networks

We analyze the correlation between the energy, momentum, and spatial entanglement produced by two luminal jets in the massive Schwinger model. Using tensor network methods, we show that for m/g > 1/π, in the vicinity of the strong- to weak-coupling transition, a nearly perfect and chargeless effective fluid behavior appears around the midrapidity region with a universal energy-pressure relationship. The evolution of energy and pressure is strongly correlated with the rise of the spatial entanglement entropy, indicating a key role of quantum dynamics. Some of these observations may be used to analyze high multiplicity jet fragmentation events, energy-energy and energy-charge correlators at current collider energies.

Ab initio calculations↗
DOE OSTI2022

Global dipole symmetry, compact Lifshitz theory, tensor gauge theory, and fractons

We study field theories with global dipole symmetries and gauge dipole symmetries. The famous Lifshitz theory is an example of a theory with a global dipole symmetry. We study in detail its 1+1d version with a compact field. When this global symmetry is promoted to a $U(1)$ dipole gauge symmetry, the corresponding gauge field is a tensor gauge field. This theory is known to lead to fractons. In order to resolve various subtleties in the precise meaning of these global or gauge symmetries, we place these 1+1d theories on a lattice and then take the continuum limit. Interestingly, the continuum limit is not unique. Different limits lead to different continuum theories, whose operators, defects, global symmetries, etc. are different. We also consider a lattice gauge theory with a $\mathbb Z_N$ dipole gauge group. Surprisingly, several physical observables, such as the ground state degeneracy and the mobility of defects depend sensitively on the number of sites in the lattice. Our analysis forces us to think carefully about global symmetries that do not act on the standard Hilbert space of the theory, but only on the Hilbert space in the presence of defects. We refer to them as time-like global symmetries and discuss them in detail. Furthermore, these time-like global symmetries allow us to phrase the mobility restrictions of defects (including those of fractons) as a consequence of a global symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

(2+1)-dimensional compact Lifshitz theory, tensor gauge theory, and fractons

The (2+1)-dimensional continuum Lifshitz theory of a free compact scalar field plays a prominent role in a variety of quantum systems in condensed matter physics and high energy physics. It is known that in compact space, it has an infinite ground-state degeneracy. In order to understand this theory better, we consider two candidate lattice regularizations of it using the modified Villain formalism. We show that these two lattice theories have significantly different global symmetries (including a dipole global symmetry), anomalies, ground-state degeneracies, and dualities. In particular, one of them is self-dual. Given these theories and their global symmetries, we can couple them to corresponding gauge theories. Here, these are two different U(1) tensor gauge theories. The resulting models have excitations with restricted mobility, i.e., fractons. Finally, we give an exact lattice realization of the fracton and lineon-elasticity dualities for the Lifshitz theory, and scalar and vector charge gauge theories.

2-dimensional systems↗
DOE OSTI2022

Improved limits on the tensor-to-scalar ratio using BICEP and P l a n c k data

In this study, we present constraints on the tensor-to-scalar ratio r using a combination of BICEP/Keck 2018 (BK18) and Planck PR4 data allowing us to fit for r consistently with the six parameters of the Λ CDM model. We discuss the sensitivity of constraints on r to uncertainties in the Λ CDM parameters as defined by the Planck data. In particular, we are able to derive a constraint on the reionization optical depth τ and thus propagate its uncertainty into the posterior distribution for r . While Planck sensitivity to r is slightly lower than the current ground-based measurements, the combination of Planck with BK18 and baryon-acoustic-oscillation data yields results consistent with r = 0 and tightens the constraint to r < 0.032 at 95% confidence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Phase diagram of generalized XY model using the tensor renormalization group

We use the higher-order tensor renormalization group method to study the two-dimensional generalized XY model that admits integer and half-integer vortices. This model is the deformation of the classical XY model and has a rich phase structure consisting of nematic, ferromagnetic, and disordered phases and three transition lines belonging to the Berezinskii-Kosterlitz-Thouless and Ising class. We explore the model for a wide range of temperatures, T , and the deformation parameter, Δ , and compute specific heat along with integer and half-integer magnetic susceptibility, finding both Berezinskii-Kosterlitz-Thouless-like and Ising-like transitions and the region where they meet. Published by the American Physical Society 2024

2-dimensional systems↗
DOE OSTI2024

Gravitational tensor-monopole moment of the hydrogen atom to order O ( α )

We calculate the gravitational tensor-monopole moment of the momentum-current density T i j in the ground state of the hydrogen atom to order O ( α ) in quantum electrodynamics (QED). The result is τ H / τ 0 − 1 = 4 α 3 π ( ln α 2 − 0.028 ) = − 3.06 × 10 − 2 where τ 0 = ℏ 2 / 4 m e is the leading-order moment. The physics of the next-to-leading-order correction is similar to that of the famous Lamb shift for energy levels. Published by the American Physical Society 2024

Ji, Xiangdong↗
DOE OSTI2025

Discrete generative diffusion models without stochastic differential equations: A tensor network approach

Diffusion models (DMs) are a class of generative machine learning methods that sample a target distribution by transforming samples of a trivial (often Gaussian) distribution using a learned stochastic differential equation. In standard DMs, this is done by learning a “score function” that reverses the effect of adding diffusive noise to the distribution of interest. Here we consider the generalisation of DMs to lattice systems with discrete degrees of freedom, and where noise is added via Markov chain jump dynamics. We show how to use tensor networks (TNs) to efficiently define and sample such “discrete diffusion models” (DDMs) without explicitly having to solve a stochastic differential equation. We show the following: (i) by parametrising the data and evolution operators as TNs, the denoising dynamics can be represented exactly; (ii) the auto-regressive nature of TNs allows to generate samples efficiently and without bias; (iii) for sampling Boltzmann-like distributions, TNs allow to construct an efficient learning scheme that integrates well with Monte Carlo. We illustrate this approach to study the equilibrium of two models with non-trivial thermodynamics, the d = 1 constrained Fredkin chain and the d = 2 Ising model. Published by the American Physical Society 2025

Causer, Luke (ORCID:0000000194243473)↗
DOE OSTI2022

Improved Limit on Tensor Currents in the Weak Interaction from Li 8 β Decay

Here, the electroweak interaction in the standard model is described by a pure vector-axial-vector structure, though any Lorentz-invariant component could contribute. In this Letter, we present the most precise measurement of tensor currents in the low-energy regime by examining the $β-\bar{ν}$ correlation of trapped 8 Li ions with the Beta-decay Paul Trap. We find $a_{β\barν}$ = - 0.3325 ± 0.0013 stat ± 0.0019 syst at 1σ for the case of coupling to right-handed neutrinos (C T = -C$^{'}_{T}$), which is consistent with the standard model prediction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Angular Correlations in the β Decay of B 8 : First Tensor-Current Limits from a Mirror-Nucleus Pair

We present the first measurement of the α–β–ν angular correlation in the Gamow-Teller β + decay of 8 B. This was accomplished using the Beta-decay Paul Trap, expanding on our previous work on the β – decay of 8 Li. The 8 B result is consistent with the V – A electroweak interaction of the standard model and, on its own, provides a limit on the exotic right-handed tensor current relative to the axial-vector current of |C T /C A | 2 < 0.013 at the 95.5% confidence level. This represents the first high-precision angular correlation measurements in mirror decays and was made possible through the use of an ion trap. Here, by combining this 8 B result with our previous 8 Li results, we demonstrate a new pathway for increased precision in searches for exotic currents.

6 ≤ A ≤ 19↗
DOE OSTI2024

High-Energy Collision of Quarks and Mesons in the Schwinger Model: From Tensor Networks to Circuit QED

With the aim of studying nonperturbative out-of-equilibrium dynamics of high-energy particle collisions on quantum simulators, we investigate the scattering dynamics of lattice quantum electrodynamics in 1 + 1 dimensions. Working in the bosonized formulation of the model and in the thermodynamic limit, we use uniform-matrix-product-state tensor networks to construct multiparticle wave-packet states, evolve them in time, and detect outgoing particles post collision. This facilitates the numerical simulation of scattering experiments in both confined and deconfined regimes of the model at different energies, giving rise to rich phenomenology, including inelastic production of quark and meson states, meson disintegration, and dynamical string formation and breaking. We obtain elastic and inelastic scattering cross sections, together with time-resolved momentum and position distributions of the outgoing particles. Furthermore, we propose an analog circuit-QED implementation of the scattering process that is native to the platform, requires minimal ingredients and approximations, and enables practical schemes for particle wave-packet preparation and evolution. Furthermore, this study highlights the role of classical and quantum simulation in enhancing our understanding of scattering processes in quantum field theories in real time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Improved Tensor Current Limit from 8 B 𝛽 Decay Including New Recoil-Order Calculations

A precision measurement of the 𝛽 + decay of 8 B was performed using the Beta-decay Paul Trap to determine the 𝛽−𝜈 angular correlation coefficient 𝑎 𝛽⁢𝜈 . The experimental results were combined with new ab initio symmetry-adapted no-core shell-model calculations to yield the second-most precise measurement from Gamow-Teller decays, 𝑎 𝛽⁢𝜈 = −0.3345 ± 0.001⁢9 stat ± 0.002⁢1 syst . This value agrees with the standard model value of −1/3 and improves uncertainties in 8 B by nearly a factor of 2. By combining results from 8 B and 8 Li , a tight limit on tensor current coupling to right-handed neutrinos was obtained. A recent global evaluation of all other precision 𝛽 decay studies suggested a nonzero value for right-handed neutrino coupling in contradiction with the standard model at just above 3⁢𝜎. Finally, the present results are of comparable sensitivity and do not support this finding.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Weak Ferromagnetism in Altermagnets from Alternating 𝑔-Tensor Anisotropy

Altermagnets are magnetic materials with antiferromagnetic spin ordering but exhibit ferromagnetic properties. Understanding the microscopic origin of the latter is a central problem. Ferromagnetlike properties such as the anomalous Hall effect are linked with weak ferromagnetism, whose microscopic origin in altermagnets remains unclear however. We show theoretically that the alternating 𝑔-tensor anisotropy in altermagnets can induce weak orbital ferromagnetism even when the Dzyaloshinskii-Moriya interaction is forbidden. We demonstrate this mechanism to explain weak ferromagnetism for both collinear and noncollinear spin altermagnets. Our findings provide new insights into the origin of weak ferromagnetism and suggest orbital-based ways for manipulating magnetic configurations in altermagnets.

Altermagnetism↗