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 325 records · Page 18

Hydrodynamics of nonintegrable systems from a relaxation-time approximation

We develop a general kinetic theory framework to describe the hydrodynamics of strongly interacting, nonequilibrium quantum systems in which integrability is weakly broken, leaving a few residual conserved quantities. This framework is based on a generalized relaxation-time approximation; it gives a simple, but surprisingly accurate, prescription for computing nonequilibrium transport even in strongly interacting systems. We validate the predictions of this approximation against matrix product operator calculations on chaotic quantum spin chains, finding surprisingly good agreement. Lastly, we show that despite its simplicity, our framework can capture phenomena distinctive to strongly interacting systems, such as widely separated charge and energy diffusion constants.

1-dimensional spin chains↗

Fission Product Evaluation (Quarterly Report)

The LLNL team has developed a theoretical framework to compute the fission fragment charge and mass distributions of odd-mass nuclei, which can be formed in neutron-induced reactions on even-even nuclei. For example, n+ 238 U leads to the compound nucleus 239 U which has 147 neutrons. Odd-mass and odd-odd nuclei are spin-polarized: the spin of the nucleus is not equal to 0 as in even-even nuclei (a consequence of pairing correlations) but is in fact set by a single particle. In practice, that single particle can occupy different quantum-mechanical configurations, each with different spin. This affects the deformation properties of the nucleus and, as a result, the distribution of fission fragments. These distributions require two inputs for the calculation: a potential energy surface (= the total energy of the nucleus as a function of the deformation of the nuclear shape, often parametrized with 2 degrees of freedom) and the collective inertia tensor associated with each point of the PES. In an odd-mass or odd-odd nucleus, each quantity depends on the configuration of the single particle. We computed of the order of 20 to 30 different such configurations to extract the PES for spins K=1/2, 3/2, 5/2, 7/2, and 9/2 of 239 U. We then computed the primary yields for each such PES, as shown in Fig. 1 below. These yields were weighed with the spin distribution associated with the reaction 238 U+n and the resulting distribution used as input to FREYA to compute the independent yields shown in Fig. 2. These yields are then used as inputs to the FIER code to extract the cumulative yields and compare with LLNL-LANL measurements at TUNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalable learning of potentials to predict time-dependent Hartree–Fock dynamics

We propose a framework to learn the time-dependent Hartree–Fock (TDHF) inter-electronic potential of a molecule from its electron density dynamics. Although the entire TDHF Hamiltonian, including the inter-electronic potential, can be computed from first principles, we use this problem as a testbed to develop strategies that can be applied to learn a priori unknown terms that arise in other methods/approaches to quantum dynamics, e.g., emerging problems such as learning exchange–correlation potentials for time-dependent density functional theory. We develop, train, and test three models of the TDHF inter-electronic potential, each parameterized by a four-index tensor of size up to 60 × 60 × 60 × 60. Two of the models preserve Hermitian symmetry, while one model preserves an eight-fold permutation symmetry that implies Hermitian symmetry. Across seven different molecular systems, we find that accounting for the deeper eight-fold symmetry leads to the best-performing model across three metrics: training efficiency, test set predictive power, and direct comparison of true and learned inter-electronic potentials. All three models, when trained on ensembles of field-free trajectories, generate accurate electron dynamics predictions even in a field-on regime that lies outside the training set. To enable our models to scale to large molecular systems, we derive expressions for Jacobian-vector products that enable iterative, matrix-free training.

97 MATHEMATICS AND COMPUTING↗

Taking Measure of the Milky Way

We intend to use SIM to make definitive measurements of fundamental structural and dynamical parameters of the Milky Way. The important niche in dynamical parameter space afforded by SIM can be exploited to resolve, with unprecedented precision, a number of classical problems of Galactic astronomy. In addition, we have developed new tests of the Galactic mass distribution specifically designed for data with the special properties of SIM products. Our proposed suite of experiments will utilize the SIM Astrometric Grid as well as complementary observations of star clusters and other strategically-selected, distant "test particles" for a definitive characterization of the major components (bulge, disk, halo, satellite system) of the Milky Way. Specifically, our goals will be: 1) The determination of two fundamental parameters that play a central role in virtually every problem in Galactic astronomy, namely (a) the solar distance to the center of the Milky Way, R(sub 0); (b) the solar angular velocity around the Galactic: center, omega(sub 0). 2) The measurement of fundamental dynamical properties of the Milky Way, among them (a) the pattern speed of the central bar (b) the rotation field and velocity-dispersion tensor in the disk (c) the kinematics (mean rotational velocity and velocity dispersion tensor) of the halo as a function of position 3. The definition of the mass distribution of the Galaxy, which is dominated by the presence of dark matter. We intend to measure (a) the relative contribution of the disk and halo to the gravitational potential (b) the local volume and surface mass density of the disk (c) the shape, mass and extent of the dark halo of the Milky Way out to 250 kpc.

Majewski, Steven R.↗

Performance Evaluations of Noisy Approximate Quantum Fourier Arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route towards a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted sum optimization problems in data processing and machine learning and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance, and identify the most effective approximation depths for quantum add and quantum multiply within the given context. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges - relevant to signed, unsigned, modular and non-modular versions - that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze performance trends in our results and speculate on possible future development within this computational paradigm.

97 MATHEMATICS AND COMPUTING↗

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Spectral and Statistical Properties of Turbulence Generated by a Vortex/Blade-Tip Interaction

The perpendicular interaction of a streamwise vortex with the tip of a lifting blade was studied in incompressible flow to provide information useful to the accurate prediction of helicopter rotor noise and the understanding of vortex dominated turbulent flows. The vortex passed 0.3 chord lengths to the suction side of the blade tip, providing a weak interaction. Single and two-point turbulence measurements were made using sub-miniature four sensor hot-wire probes 15 chord lengths downstream of the blade trailing edge; revealing the mean velocity and Reynolds stress tensor distributions of the turbulence, as well as its spanwise length scales as a function of frequency. The single point measurements show the flow downstream of the blade to be dominated by the interaction of the original tip vortex and the vortex shed by the blade. These vortices rotate about each other under their mutual induction, winding up the turbulent wakes of the blades. This interaction between the vortices appears to be the source of new turbulence in their cores and in the region between them. This turbulence appears to be responsible for some decay in the core of the original vortex, not seen when the blade is removed. The region between the vortices is not only a region of comparatively large stresses, but also one of intense turbulence production. Velocity autospectra measured near its center suggests the presence quasi-periodic large eddies with axes roughly parallel to a line joining the vortex cores. Detailed two-point measurements were made on a series of spanwise cuts through the flow so as to reveal the turbulence scales as they would be seen along the span of an intersecting airfoil. The measurements were made over a range of probe separations that enabled them to be analyzed not only in terms of coherence and phase spectra but also in terms of wave-number frequency (kappa-omega) spectra, computed by transforming the measured cross-spectra with respect to the spanwise separation of the probes. These data clearly show the influence of the coherent eddies in the spiral wake and the turbulent region between the cores. These eddies produce distinct peaks in the upwash velocity kappa-omega spectra, and strong anisotropy manifested both in the decay of the kappa-omega spectrum at larger wave-numbers and in differences between the kappa-omega spectra of different components. None of these features are represented in the von Karman spectrum for isotropic turbulence that is often used in broadband noise computations. Wave-number frequency spectra measured in the cores appear to show some evidence that the turbulence outside sets tip core waves, as has previously been hypothesized. These spectra also provide for the first time a truly objective method for distinguishing velocity fluctuations produced by core wandering from other motions.

Devenport, William J.↗

A two-loop four-point form factor at function level

Recently, the maximally-helicity-violating four-point form factor for the chiral stress-energy tensor in planar $\mathcal{N}$ = 4 super Yang-Mills was computed to three loops at the level of the symbol associated with multiple polylogarithms. It exhibits antipodal self-duality, or invariance under the combined action of a kinematic map and reversing the ordering of letters in the symbol. Here we lift the two-loop form factor from symbol level to function level. We provide an iterated representation of the function’s derivatives (coproducts). In order to do so, we find a three-parameter limit of the five-parameter phase space where the symbol’s letters are all rational. We also use function-level information about dihedral symmetries and the soft, collinear, and factorization limits, as well as limits governed by the form-factor operator product expansion (FFOPE). We provide plots of the remainder function on several kinematic slices, and show that the result is compatible with the FFOPE data. We further verify that antipodal self-duality is valid at two loops beyond the level of the symbol.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Novel Solid Form of Erlotinib: Synthesis by Heterogeneous Complexation and Characterization by NMR Crystallography

We describe the synthesis of a novel complex of the anticancer "active pharmaceutical ingredient erlotinib (ERL) via heterogeneous nucleation on polished zinc tiles. The resulting product, ERL 2 ·ZnCl 2 , is characterized by single-crystal X-ray diffraction, multinuclear solid-state NMR (ssNMR) spectroscopy, and density functional theory (DFT) calculations. Also characterized are the hydrochloride salt (ERL·HCl) and monohydrate free base (ERL·H 2 O) forms of erlotinib. 13 C ssNMR spectroscopy is useful for site-by-site assignment and rapid fingerprinting, while also providing preliminary structural interpretations, such as the number of molecules in the asymmetric unit. 35 Cl ssNMR can readily differentiate between the chloride ions in ERL·HCl and the covalently bonded chlorine in ERL 2 ·ZnCl 2 . 15 N ssNMR proves to be critical here because of the large isotropic chemical shift differences between ERL·H 2 O, ERL·HCl, and ERL 2 ·ZnCl 2 . The 15 N chemical shift tensors are linked directly to differences in structure and bonding with the aid of DFT calculations. Together, these results demonstrate the utility of multinuclear NMR crystallography for the characterization of solid forms of APIs, especially when other analytical techniques face significant challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

WUS256: An Adjoint Waveform Tomography Model of the Crust and Upper Mantle of the Western United States for Improved Waveform Simulations

Abstract We report a new model (WUS256) of radially anisotropic seismic wavespeeds of the crust and upper mantle of the western United States (WUS) obtained from adjoint waveform tomography for the purpose of improving synthetic waveform fits to observed data. WUS256 is based on inversion of over 94,000 waveforms from 72 earthquakes recorded by nearly 3,400 stations. We started with the SPiRaL global model (Simmons et al., 2021, https://doi.org/10.1093/gji/ggab277 ) and waveforms in the period band of 50–120 s. We followed a conservative multiscale inversion approach with eight stages and 256 total inversion iterations which enabled monotonic misfit reduction to 20‐s minimum‐period waves. WUS256 relied on time‐frequency (TF) phase misfits and a trust region limited memory Broyden–Fletcher–Goldfarb–Shanno (L‐BFGS) optimization. Hessian‐vector products were used to qualitatively assess model resolution. Results indicate that WUS256 has good coverage of the continental regions to depths of about 150 km and is able to resolve features on lateral scales of about 200 km. We quantify waveform fits by the reduction in TF and normalized amplitude difference misfits between WUS256 and the SPiRaL starting model. WUS256 significantly improves waveform fits with misfit reduction 64% for both inversion and validation data sets compared to the SPiRaL starting model and shows even better fits compared to other models. Waveform fits illustrate that WUS256 reproduces body‐waves, fundamental mode surface waves as well as late arriving dispersed and/or scattered short period surface waves. The improvement in waveform fit indicates that WUS256 can be used to reproduce path effects on regional complete waveforms and moment tensor inversions.

58 GEOSCIENCES↗

Far-from-equilibrium slow modes and momentum anisotropy in expanding plasma

The momentum distribution of particle production in heavy-ion collisions encodes information about thermalization processes in the early-stage quark-gluon plasma. We use kinetic theory to study the far-from-equilibrium evolution of an expanding plasma with an anisotropic momentumspace distribution. We identify slow and fast degrees of freedom in the far-from-equilibrium plasma from the evolution of moments of this distribution. At late times, the slow modes correspond to hydrodynamic degrees of freedom and are naturally gapped from the fast modes by the inverse of the relaxation time, $τ$$^{–1}_{R}$. At early times, however, there are an infinite number of slow modes with a gap inversely proportional to time, τ –1 . From the evolution of the slow modes we generalize the paradigm of the far-from-equilibrium attractor to vector and tensor components of the energy momentum tensor, and even to higher moments of the distribution function that are not part of the hydrodynamic evolution. We predict that initial-state momentum anisotropy decays slowly in the far-from-equilibrium phase and may persist until the relaxation time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for charged-lepton flavour violation in top quark interactions with an up-type quark, a muon, and a $τ$ lepton in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for charged-lepton flavour violation (CLFV) in top quark (t) production and decay is presented. The search uses proton-proton collision data corresponding to 138 fb$^{-1}$ collected with the CMS experiment at $\sqrt{s}$ = 13 TeV. The signal consists of the production of a single top quark via a CLFV interaction or top quark pair production followed by a CLFV decay. The analysis selects events containing a hadronically decaying $τ$ lepton and a muon of opposite electric charge, as well as at least three jets, one of which is identified as originating from the fragmentation of a bottom quark. Machine learning classification techniques are used to distinguish signal from standard model background events. The results of this search are consistent with the standard model expectations. The upper limits at 95% confidence level on the branching fraction $\mathcal{B}$ for CLFV top quark decays to a muon, a $τ$ lepton, and an up or a charm quark are set at $\mathcal{B}$(t $\to$ $μτ$u) $\lt$ (0.04, 0.08, and 0.12) $\times$ 10$^{-6}$, and $\mathcal{B}$(t $\to$ $μτ$c) $\lt$ (0.81, 1.71, and 2.05) $\times$ 10$^{-6}$ for scalar, vector, and tensor-like operators, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Far-from-equilibrium slow modes and momentum anisotropy in an expanding plasma

The momentum distribution of particle production in heavy-ion collisions encodes information about thermalization processes in the early-stage quark-gluon plasma. We use kinetic theory to study the far-from-equilibrium evolution of an expanding plasma with an anisotropic momentum-space distribution. We identify slow and fast degrees of freedom in the far-from-equilibrium plasma from the evolution of moments of this distribution. At late times, the slow modes correspond to hydrodynamic degrees of freedom and are naturally gapped from the fast modes by the inverse of the relaxation time, τ R − 1 . At early times, however, there are an infinite number of slow modes with a gap inversely proportional to time, τ − 1 . From the evolution of the slow modes we generalize the paradigm of the far-from-equilibrium attractor to vector and tensor components of the energy-momentum tensor, and even to higher moments of the distribution function that are not part of the hydrodynamic evolution. We predict that initial-state momentum anisotropy decays slowly in the far-from-equilibrium phase and may persist until the relaxation time. Published by the American Physical Society 2024

Brewer, Jasmine (ORCID:0000000230840663)↗

Visualizing the SEI formation between lithium metal and solid-state electrolyte

The solid electrolyte interphase (SEI) is regarded as the most important factor affecting the durability of lithium-metal anode in all-solid-state batteries (ASSBs). Despite its significance, the nucleation and growth mechanism of SEI is not yet well understood. Here, we elucidate the thermodynamics and kinetics governing SEI formation at the Li|β-Li 3 PS 4 interface at the atomic scale via thermodynamic phase equilibrium analysis and machine-learning-potential-assisted molecular dynamics (MD) simulations. An accurate moment tensor potential using the machine-learning method is developed for a reactive model of Li|β-Li 3 PS 4 . This potential enabled us to perform large-scale MD simulations with the model size expanded to the experimental dimensions (~40 nm) while maintaining the same level of accuracy as density functional theory calculations. The results reveal a four-stage evolution process at the Li|β-Li 3 PS 4 interface, namely (i) fast ion diffusion, (ii) nucleation, (iii) Li 2 S growth, and (iv) stabilization. Notably, we demonstrate that the SEI can be categorized into crystalline and amorphous regions. The simulated SEI thickness, structure, and composition closely match experimental findings, validating the accuracy of the MD simulations. We further disclose the significant impact of ion diffusion kinetic limitations on the phase formation and crystallization of interfacial products. Furthermore, we shed light on the detailed potential energy (PE) distribution of lithium along the direction perpendicular to the interface. Finally, this information is crucial for better understanding interfacial ion mobility.

25 ENERGY STORAGE↗

Performance evaluations of signed and unsigned noisy approximate quantum Fourier arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route toward a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted-sum optimization problems in data processing and machine learning, and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting-qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance and identify the most effective approximation depths for unsigned quantum addition and quantum multiplication within the given context. We then perform a similar analysis of signed addition and compare to the unsigned results. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges—relevant to signed, unsigned, modular and non-modular versions—that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze the performance trends in our results and speculate on possible future developments within this computational paradigm.

Computational models↗

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Extracting off-diagonal order from diagonal basis measurements

Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multiscale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-dimensional extended Hubbard model with an off-diagonal order parameter, sandwiched between diagonal charge and spin density wave phases, using only diagonal descriptors. We study the model directly in the thermodynamic limit using the recently developed variational uniform matrix product states algorithm, and draw our samples from degenerate ground states related by global spin rotations, emulating the projective measurements that are accessible in experiments. Our results will have important implications for the study of exotic phases using optical lattice experiments. Published by the American Physical Society 2024

1-dimensional systems↗

Structural Studies of Alloyed and Nanoparticle Transition Metal Dichalcogenides by Selenium-77 Solid-State Nuclear Magnetic Resonance Spectroscopy

Layered transition metal dichalcogenides (TMDCs) such as MoS 2 , MoSe 2 , and WSe 2 are under intense investigation because they are atomically thin semiconductors with photophysical properties that can be tuned by changing their composition or morphology. Mechanochemical processing has been proposed as a method to obtain alloyed TMDCs in the series Mo 1–x W x S y Se 2–y (x = 0–1; y = 0–2). However, elucidating the chemical transformations occurring at the atomic scale following mechanochemical processing can be challenging because the products are often amorphous or microcrystalline. To address this challenge, we probe TMDC mixing and alloying by using a combination of powder X-ray diffraction, Raman spectroscopy, diffuse reflectance spectroscopy, 77 Se solid-state nuclear magnetic resonance (SSNMR) spectroscopy, and planewave density functional theory (DFT) calculations. The nature of the milling material and reaction atmosphere are shown to be essential factors in limiting the formation of undesired oxide byproducts. We demonstrate acquisition of 77 Se SSNMR spectra using different combinations of Carr-Purcell Meiboom-Gill acquisition (CPMG) pulse sequences, magic angle spinning (MAS), and MAS dynamic nuclear polarization. Further, the combination of SSNMR with the other characterization methods clearly demonstrates that high energy impact ball milling induces molecular level alloying of Mo, W and chalcogen atoms in the family Mo 1–x W x S y Se 2–y . Gauge including projector augmented wave DFT calculations yield accurate 77 Se chemical shift (CS) tensor components. 77 Se SSNMR spectroscopy was also applied to study the structure of WSe 2 nanocrystals intercalated with ethylenediamine. The intercalated WSe 2 nanocrystals exhibit a more positive isotropic 77 Se CS as compared to bulk WSe 2 , however, the 77 Se CS anisotropy is the same, confirming the WSe 2 layers have a similar structure as in their bulk counterparts.

36 MATERIALS SCIENCE↗