Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor product”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Distributed-Memory Sparse Deep Neural Network Inference Using Global Arrays

Partitioned Global Address Space (PGAS) models exhibit tremendous promise in developing efficient and productive distributed-memory parallel applications. They have been used extensively in scientific computations due to conveniently offering a ``shared-memory''-like model and convenient interfaces that separate communication with synchronization. Traditionally, PGAS communication models have been applied to dense/contiguously distributed data, but most modern applications depict varied levels of sparsity. Existing PGAS models require certain adaptations to support distributed sparse computations, since associated computations often require matrix arithmetic, in addition to data movement. The Global Arrays toolkit from Pacific Northwest National Laboratory (PNNL) is one of the earliest PGAS models to combine one-sided data communication and distributed matrix operations and is still used in the popular NWChem quantum chemistry suite. Recently, we have expanded the Global Arrays toolkit to support common sparse operations, like sparse matrix-dense matrix multiplies (SpMM), sparse matrix-sparse matrix multiplication (SpGEMM) and Sampled Dense-Dense Matrix Multiplication (SDDMM). As it turns out, these operations are the bedrock of sparse Deep Learning (DL); sparse deep neural networks and Graph Neural Networks (GNNs) have gained increasing attention recently in achieving speedups on training and inference with reduced memory footprints. Unlike scientific applications in High Performance Computing (HPC), modern (distributed-memory capable) DL toolkits often rely on non-standardized and closed-source vendor software optimizations, creating challenges in software-hardware co-design at scale. Our goal is to support a variety of distributed-memory sparse matrix operations and helper functions in the newly created Sparse Global Arrays (SGA), such that it is possible to build portable and productive Machine Learning scenarios for algorithm/software and hardware codesign purposes. Contemporary data-parallel schemes for training/inference are undergoing a major overhaul since model replication limits scalability and causes resource inefficiencies. As such, we have adopted tensor parallelism in decomposing the model and inputs, to mitigate memory issues. Current implementation is built on top of MPI and uses CPUs to maximize the portability across the platforms.

Distributed computing, machine learning↗
DOE OSTI2021

Entanglement and many-body effects in collective neutrino oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings, such as supernovae, where the neutrino density is large. In this regime, neutrino-neutrino interactions are important and simulations in the mean-field approximation show evidence for collective oscillations occurring at timescales much shorter than those associated with vacuum oscillations. In this work, we study the out-of-equilibrium dynamics of a corresponding spin model using matrix product states and show how collective bipolar oscillations can be triggered by many-body correlations if appropriate initial conditions are present. We find entanglement entropies scaling at most logarithmically in the system size suggesting that classical tensor network methods could be efficient in describing collective neutrino dynamics more generally. These observation provide a clear path forward, not only to increase the accuracy of current simulations, but also to elucidate the mechanism behind collective flavor oscillations without resorting to the mean-field approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Overview of SoLID

SoLID spectrometer was proposed to fully exploit the potential of JLab 12 GeV energy upgrade. It is a large acceptance detector which can handle very high luminosity. An overview of the rich physics program will be given, which includes number of planned measurements: a multi-dimensional mapping of semi-inclusive DIS asymmetries for tomography of the nucleon in momentum space in the high-x region and precision determination of the tensor charge; a measurement of parity-violating DIS to provide a precision test of the Standard Model, reaching a sensitivity to new physics at 10-20 TeV level; a precision measurement of J/psi photo- and electro-production cross sections in the threshold region to probe the strong color fields in the nucleon and to study the origin of the proton mass. The current status and the plan of the project will be discussed.

Chen, Jian-Ping↗
DOE OSTI2025

A quantum eigenvalue solver based on tensor networks

Electronic ground states are of central importance in chemical simulations, but have remained beyond the reach of efficient classical algorithms except in cases of weak electron correlation or one-dimensional spatial geometry. We introduce a hybrid quantum-classical eigenvalue solver that constructs a wavefunction ansatz from a linear combination of matrix product states in rotated orbital bases, enabling the characterization of strongly correlated ground states with arbitrary spatial geometry. The energy is converged via a gradient-free generalized sweep algorithm based on quantum subspace diagonalization, with a potentially exponential speedup in the off-diagonal matrix element contractions upon translation into compact quantum circuits of linear depth in the number of qubits. Chemical accuracy is attained in numerical experiments for both a stretched water molecule and an octahedral arrangement of hydrogen atoms, achieving substantially better correlation energies compared to a unitary coupled-cluster benchmark, with orders of magnitude reductions in quantum resource estimates and a surprisingly high tolerance to shot noise. This proof-of-concept study suggests a promising new avenue for scaling up simulations of strongly correlated chemical systems on near-term quantum hardware.

chemistry↗
DOE OSTI2023

Quasi-Fermi liquid behavior in a one-dimensional system of interacting spinless fermions

We present numerical evidence for a paradigm in one-dimensional interacting fermion systems, whose phenomenology has traits of both Luttinger liquids and Fermi liquids. This state, dubbed a quasi-Fermi liquid, possesses a discontinuity in its fermion occupation number at the Fermi momentum. The excitation spectrum presents particlelike quasiparticles and absence of holelike quasiparticles, giving rise instead to edge singularities. Such a state is realized in a one-dimensional spinless fermion lattice Hamiltonian by fine-tuning the interactions to a regime where they become irrelevant in the renormalization group sense. We show, using uniform infinite matrix products states and finite-entanglement scaling analysis, that the system ground state is characterized by a Luttinger parameter K = 1 and a discontinuous jump in the fermion occupation number. We support the characterization with calculations of the spectral function that show a particle-hole asymmetry reflected in the existence of well-defined Landau quasiparticles above the Fermi level and edge singularities without the associated quasiparticles below. Furthermore, these results indicate that the quasi-Fermi liquid paradigm can be realized beyond the low-energy perturbative realm.

1-dimensional systems↗
DOE OSTI2022

A universal method to compare parts from STEP files

Abstract Model Based Definition (MBD) captures the complete specification of a part in digital form and leverages (at least) the universal “Standard for the Exchange of Product” (STEP) file format. MBD has revolutionized manufacturing due to time and cost savings associated with containing all engineering data within a single digital source. This work presents a novel method to transform digital definitions in any given STEP file into a tensor-like structure that is unique for each part and can be used to regenerate the original STEP file completely. Resulting STEP tensors are amenable to part comparison based on various part specifications in a general and straightforward manner. Here, part similarity is evaluated among sets of parts according to specific geometry, material composition, and design intent. Importantly, specification similarity can be quantified using only the tensors’ structure. As such, this approach is not limited to families of geometric shapes, part types, or fabrication methods; nor does it require any prior knowledge about the parts being compared.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Conformal conserved currents in embedding space

We study conformal conserved currents in arbitrary irreducible representations of the Lorentz group using the embedding space formalism. With the help of the operator product expansion, we first show that conservation conditions can be fully investigated by considering only two- and three-point correlation functions. We then find an explicitly conformally-covariant differential operator in embedding space that implements conservation based on the standard position space operator product expansion differential operator ∂μ, although the latter does not uplift to embedding space covariantly. The differential operator in embedding space that imposes conservation is the same differential operator $D_{ijA}$ used in the operator product expansion in embedding space. We provide several examples including conserved currents in irreducible representations that are not symmetric and traceless. With an eye on four-point conformal bootstrap equations for four conserved vector currents $\langle JJJJ\rangle$ and four energy-momentum tensors $\langle TTTT\rangle$, we mostly focus on conservation conditions for $\langle JJ\mathcal{O}\rangle$ and $\langle TT\mathcal{O}\rangle$. Finally, we reproduce and extend the consequences of conformal Ward identities at coincident points by determining three-point coefficients in terms of charges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Separable shape tensors for aerodynamic design

Airfoil shape design is a classical problem in engineering and manufacturing. In this work, we combine principled physics-based considerations for the shape design problem with modern computational techniques using a data-driven approach. Modern and traditional analyses of two-dimensional (2D) and three-dimensional (3D) aerodynamic shapes reveal a flow-based sensitivity to specific deformations that can be represented generally by affine transformations (rotation, scaling, shearing, and translation). We present a novel representation of shapes that decouples affine-style deformations over a submanifold and a product submanifold principally of the Grassmannian. As an analytic generative model, the separable representation, informed by a database of physically relevant airfoils, offers: (i) a rich set of novel 2D airfoil deformations not previously captured in the data, (ii) an improved low-dimensional parameter domain for inferential statistics informing design/manufacturing, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal 2D shapes.

42 ENGINEERING↗
DOE OSTI2023

Gravitational wave probes of massive gauge bosons at the cosmological collider

We extend the reach of the "cosmological collider" for massive gauge boson production during inflation from the CMB scales to the interferometer scales. Considering a Chern-Simons coupling between the gauge bosons and the pseudoscalar inflaton, one of the transverse gauge modes is efficiently produced and its inverse decay leaves an imprint in the primordial scalar and tensor perturbations. We study the correlation functions of these perturbations and derive the updated constraints on the parameter space from CMB observables. We then extrapolate the tensor power spectrum to smaller scales consistently taking into account the impact of the gauge field on inflationary dynamics. Our results show that the presence of massive gauge fields during inflation can be detected from characteristic gravitational wave signals encompassing the whole range of current and planned interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Drell-Yan angular lepton distributions at small x from TMD factorization.

The Drell-Yan process is studied in the framework of TMD factorization in the Sudakov region s » Q 2 » \( {q}_{\perp}^2 \) corresponding to recent LHC experiments with Q 2 of order of mass of Z-boson and transverse momentum of DY pair ~ few tens GeV. The DY hadronic tensors are expressed in terms of quark and quark-gluon TMDs with \( \frac{1}{Q^2} \) and \( \frac{1}{N_c^2} \) accuracy. It is demonstrated that in the leading order in N c the higher-twist quark-quark-gluon TMDs reduce to leading-twist TMDs due to QCD equation of motion. The resulting hadronic tensors depend on two leading-twist TMDs: f 1 responsible for total DY cross section, and Boer-Mulders function \( {h}_1^{\perp } \) . The corresponding qualitative and semi-quantitative predictions seem to agree with LHC data on five angular coefficients A 0 – A 4 of DY pair production. The remaining three coefficients A 5 – A 7 are determined by quark-quark-gluon TMDs multiplied by extra \( \frac{1}{N_c} \) so they appear to be relatively small in accordance with LHC results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Triviality of quantum trajectories close to a directed percolation transition

Here we study quantum circuits consisting of unitary gates, projective measurements, and control operations that steer the system toward a pure absorbing state. Two types of phase transition occur as the rate of these control operations is increased: a measurement-induced entanglement transition, and a directed percolation transition into the absorbing state (taken here to be a product state). In this work, we show analytically that these transitions are generically distinct, with the quantum trajectories becoming disentangled before the absorbing state transition is reached, and we analyze their critical properties. We introduce a simple class of models where the measurements in each quantum trajectory define an effective tensor network (ETN)—a subgraph of the initial spacetime graph where nontrivial time evolution takes place. By analyzing the entanglement properties of the ETN, we show that the entanglement and absorbing-state transitions coincide only in the limit of the infinite local Hilbert-space dimension. Focusing on a Clifford model which allows numerical simulations for large system sizes, we verify our predictions and study the finite-size crossover between the two transitions at large local Hilbert space dimension. We give evidence that the entanglement transition is governed by the same fixed point as in hybrid circuits without feedback.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Models of $J/\varPsi$ photo-production reactions on the nucleon

The J/Ψ photo-production reactions on the nucleon can provide information on the roles of gluons in determining the J/Ψ-nucleon (J/Ψ-N) interactions and the structure of the nucleon. The information on the J/Ψ-N interactions is needed to test lattice QCD (LQCD) calculations and to understand the nucleon resonances such as N*(P c ) recently reported by the LHCb Collaboration. In addition, it is also needed to investigate the production of nuclei with hidden charms and to extract the gluon distributions in nuclei. The main purpose of this article is to review six models of the γ+p→J/Ψ+p reaction which have been and can be applied to analyze the data from Thomas Jefferson National Accelerator Facility (JLab). Here, the formulae for each model are given and used to obtain the results to show the extent to which the available data can be described. The models presented include the Pomeron-exchange model of Donnachie and Landshoff (Pom-DL) and its extensions to include J/Ψ-N potentials extracted from LQCD (Pom-pot) and to also use the constituent quark model (CQM) to account for the quark substructure of J/Ψ (Pom-CQM). The other three models are developed from applying the perturbative QCD approach to calculate the two-gluon exchange using the generalized parton distribution (GPD) of the nucleon (GPD-based), two- and three-gluon exchanges using the parton distribution of the nucleon (2g+3g), and the exchanges of scalar (0 ++ ) and tensor (2 ++ ) glueballs within the holographic formulation (holog). The results of investigating the excitation of the nucleon resonances N*(P c ) in the γ+p→J/Ψ+p reactions are also given. We demonstrate that the differences between these six models can be unambiguously distinguished and the N* can be better studied by using the forthcoming JLab data at large |t| and at energies very near the J/Ψ production threshold. Possible improvements of the considered models are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Flexible Power Method for Solving Infinite Dimensional Tensor Eigenvalue Problems

We propose a flexible power method for computing the leftmost, i.e., algebraically smallest, eigenvalue of an infinite dimensional tensor eigenvalue problem, $H x = \lambda x$, where the infinite dimensional symmetric matrix $H$ exhibits a translational invariant structure. We assume the smallest eigenvalue of $H$ is simple and apply a power iteration of $e^{-H}$ with the eigenvector represented in a compact way as a translational invariant infinite Tensor Ring (iTR). Hence, the infinite dimensional eigenvector can be represented by a finite number of iTR cores of finite rank. In order to implement this power iteration, we use a small parameter $t$ so that the infinite matrix-vector operation $e^{-Ht}x$ can efficiently be approximated by the Lie product formula, also known as Suzuki--Trotter splitting, and we employ a low rank approximation through a truncated singular value decomposition on the iTR cores in order to keep the cost of subsequent power iterations bounded. We also use an efficient way for computing the iTR Rayleigh quotient and introduce a finite size iTR residual which is used to monitor the convergence of the Rayleigh quotient and to modify the timestep $t$. In this paper, we discuss 2 different implementations of the flexible power algorithm and illustrate the automatic timestep adaption approach for several numerical examples.

Beeumen, Roel Van↗
DOE OSTI2023

Conceptual design of DIII-D experiments to diagnose the lifetime of spin polarized fuel

Abstract In magnetic fusion experiments, the cross sections for the D-T and D- 3 He fusion reactions are increased by as much as 50% if the fuel remains spin polarized parallel to the magnetic field. The goal of this study is to assess the feasibility of lifetime measurements of spin polarization, in magnetic fusion relevant conditions, on the DIII-D tokamak using relative changes in charged fusion product (CFP) loss measurements that depend upon the differential fusion cross section d σ /dΩ. Relative measurements that capture changes in the escaping CFP pitch, poloidal, and energy distributions are studied in two realistic TRANSP calculated plasma scenarios: (a) vector-polarized 3 He and D pellets are injected into a hot hydrogen plasma to produce thermonuclear reactions, and (b) a tensor-polarized deuterium pellet is injected into an L-mode hydrogen background plasma that includes unpolarized 3 He neutral beam injection. Ideal CFP signals in both scenarios show substantial pitch sensitivity to polarization for 14.7 MeV proton detection at a poloidal angle of − 56 ∘ (on the outer wall in the ion ∇ B direction), and strong sensitivity to polarization for 3.6 MeV alpha flux detection by an array of poloidal detectors. Energy-resolved measurements of 14.7 MeV protons are also sensitive to the degree of polarization for the − 56 ∘ port in the beam-plasma scenario. A realistic assessment of CFP signals in the thermonuclear scenario show count rates in the range of 2 × 10 4 cps for pitch-resolved proton detection and 2 × 10 5 cps for alpha flux measurements. Reduced chi-squared χ r 2 calculations show polarization lifetime measurements are feasible for either proton or alpha detection of both enhanced or suppressed polarization for the thermonuclear scenario. Measurements of gamma rays produced in the weak D + 3 He → γ + 5 Li branch complement the CFP measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Search for charged-lepton flavor violation in top quark production and decay in pp collisions at $ \sqrt{s} $ = 13 TeV

Results are presented from a search for charged-lepton flavor violating (CLFV) interactions in top quark production and decay in pp collisions at a center-of-mass energy of 13 TeV. The events are required to contain one oppositely charged electron-muon pair in the final state, along with at least one jet identified as originating from a bottom quark. The data correspond to an integrated luminosity of 138 fb$^{-1}$, collected by the CMS experiment at the LHC. This analysis includes both the production (q → eμt) and decay (t → eμq) modes of the top quark through CLFV interactions, with q referring to a u or c quark. These interactions are parametrized using an effective field theory approach. With no significant excess over the standard model expectation, the results are interpreted in terms of vector-, scalar-, and tensor-like CLFV four-fermion effective interactions. Finally, observed exclusion limits are set at 95% confidence levels on the respective branching fractions of a top quark to an eμ pair and an up (charm) quark of 0.13 × 10$^{-6}$ (1.31 × 10$^{-6}$), 0.07 × 10$^{-6}$ (0.89 × 10$^{-6}$), and 0.25 × 10$^{-6}$ (2.59 × 10$^{-6}$) for vector, scalar, and tensor CLFV interactions, respectively.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Semi-inclusive deep-inelastic scattering on a polarized spin-1 target. I. Cross section and spin observables

We develop the theoretical framework for semi-inclusive deep-inelastic scattering on a polarized spin-1 target and apply it to scattering on the polarized deuteron with spectator nucleon tagging. In Part I (this article) we present the general form of the semi-inclusive cross section and polarization observables for the spin-1 target. A relativistically covariant formulation in terms of four-vectors and invariant polarization parameters is employed. The target polarization is described by a spin density matrix with vector and tensor polarization. The spin and azimuthal angle dependence of the semi-inclusive cross section is derived and parametrized in terms of invariant structure functions. To validate the result, the structure functions are expressed as photon-target helicity amplitudes with known symmetry properties. The expressions presented here are kinematic (no assumptions about particle production dynamics) and valid in all regions of the deep-inelastic final state (current and target fragmentation regions). In Part II (following article), we consider deep-inelastic scattering on the polarized deuteron with spectator nucleon tagging as a special case of target fragmentation. The semi-inclusive structure functions are computed by separating nuclear and hadronic structure, and the polarization observables are explored as functions of the tagged nucleon momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗