SEARCH · Search NASA

Results for “Matrix exponential”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Kinematic modeling of the Galaxy. I - The Yale Bright Star Catalogue

A code has been developed to investigate the kinematic properties of stellar populations of the Galaxy by modeling the distributions of proper motion and line-of-sight velocity. The coordinate frames of reference are discussed in detail and the general matrix transformations valid beyond the solar neighborhood are given. As a first application, the kinematic properties of the magnitude-selected sample of 3993 stars in the Bright Star Catalogue with V less than 6.0 and (B-V) greater than 0 are selected. The observed kinematic distributions are compared with a conventional model that has an exponential disk with scale heights, components of velocity dispersion ellipsoid, and asymmetric drift taken from the literature and defined as functions of spectral type. The comparisons suggest that a significant fraction of the Bright Star Catalogue belongs to a relatively young population that is not properly represented in the model.

Ratnatunga, Kavan U.↗

DOE OSTI2024

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Characterizing random-singlet state in two-dimensional frustrated quantum magnets and implications for the double perovskite Sr 2 CuTe 1 – x W x O 6

Motivated by the experimental observation of a nonmagnetic phase in compounds with frustration and disorder, we study the ground state of a spin-1/2 square-lattice Heisenberg model with randomly distributed nearest-neighbor J 1 and next-nearest-neighbor J 2 couplings. By using the density matrix renormalization group (DMRG) calculation on a cylinder system with a circumference of up to ten lattice sites, we identify a disordered phase between the Néel and stripe magnetic phase with growing J 2 /J 1 in the presence of strong bond randomness. The vanished spin-freezing parameter indicates the absence of spin-glass order. The large-scale DMRG results unveil the size-scaling behaviors of the spin-freezing parameter, the power-law decay of the average spin correlation, and the exponential decay of the typical spin correlation, which all agree with the corresponding behavior in the one-dimensional random-singlet (RS) state and characterize the RS nature of this disordered phase. The DMRG simulation also provides insights and opportunities for characterizing a class of nonmagnetic states in two-dimensional frustrated magnets with disorder. Here, we also compare with existing experiments and suggest more measurements for understanding the spin-liquid-like behaviors in the double perovskite Sr 2 CuTe 1–x W x O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Quantum error correction from complexity in Brownian SYK

We study the robustness of quantum error correction in a one-parameter ensemble of codes generated by the Brownian SYK model, where the parameter quantifies the encoding complexity. The robustness of error correction by a quantum code is upper bounded by the “mutual purity” of a certain entangled state between the code subspace and environment in the isometric extension of the error channel, where the mutual purity of a density matrix ρAB is the difference $\mathcal{F}$ p ($A : B$) ≡ $\mathrm{T}$r $p^{2}_{AB}$ - $\mathrm{T}$r $p^{2}_{A}$ $\mathrm{T}$r $p^{2}_{B}$. We show that when the encoding complexity is small, the mutual purity is O(1) for the erasure of a small number of qubits (i.e., the encoding is fragile). However, this quantity decays exponentially, becoming O(1/N) for O(log N) encoding complexity. Further, at polynomial encoding complexity, the mutual purity saturates to a plateau of O(e -N ). We also find a hierarchy of complexity scales associated to a tower of subleading contributions to the mutual purity that quantitatively, but not qualitatively, adjust our error correction bound as encoding complexity increases. In the AdS/CFT context, our results suggest that any portion of the entanglement wedge of a general boundary subregion A with sufficiently high encoding complexity is robustly protected against low-rank errors acting on A with no prior access to the encoding map. From the bulk point of view, we expect such bulk degrees of freedom to be causally inaccessible from the region A despite being encoded in it.

1/N expansion↗

DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

DOE OSTI2026

Deuterated Xylene-Based Neutron Spectrometry

Deuterated scintillators, such as deuterated xylene (C8D10) have been theorized to offer better pulse shape discrimination (PSD) and neutron spectrometry performance than traditional protiated scintillators due to the anisotropy of neutron scattering on deuterium. However, the light output responses of said deuterated scintillators to monoenergetic has not been characterized in its entirety at both high and low energies. Several time-of-flight experiments were conducted at the Idaho national Laboratory to characterize a pair of EJ-301D deuterated xylene detectors. The low energy detector response to quasi-monoenergetic neutrons was measured using a Cf-252 source, and the high energy detector response was measured using a D-T source and angles of 20, 45, and 60 deg, corresponding to theoretical neutron energies of 13.2, 10.4, and 8.2 MeV. During post-processing, the neutron and gamma bands were clearly able to be separated, even at low light outputs, solidifying the superior PSD capability of deuterated xylene. An exponential of fit of LO=0.63E-2.22*(1-exp(-0.25E)), which has good agreement with previously measured data, was found as the light output response from deuteron recoils. Additionally, measurements of AmBe, D-T, and D-D sources were unfolded using the light output response curve, which resulted in favorable agreement with theoretical values. Further testing of the response matrix by unfolding more sources (such as PuBe) needs to be done, but the unfolding accuracy for currently available sources suggests that deuterated xylene can be used to measure and unfold both fission, monoenergetic, and continuous spectra.

scintillator↗

NASA NTRS1993

A transfer matrix approach to vibration localization in mistuned blade assemblies

A study of mode localization in mistuned bladed disks is performed using transfer matrices. The transfer matrix approach yields the free response of a general, mono-coupled, perfectly cyclic assembly in closed form. A mistuned structure is represented by random transfer matrices, and the expansion of these matrices in terms of the small mistuning parameter leads to the definition of a measure of sensitivity to mistuning. An approximation of the localization factor, the spatially averaged rate of exponential attenuation per blade-disk sector, is obtained through perturbation techniques in the limits of high and low sensitivity. The methodology is applied to a common model of a bladed disk and the results verified by Monte Carlo simulations. The easily calculated sensitivity measure may prove to be a valuable design tool due to its system-independent quantification of mistuning effects such as mode localization.

Ottarson, Gisli↗

NASA NTRS1983

Dissociative recombination of O2(+), NO(+) and N2(+)

A new L(2) approach for the calculation of the threshold molecular capture width needed for the determination of DR cross sections was developed. The widths are calculated with Fermi's golden rule by substituting Rydberg orbitals for the free electron continuum coulomb orbital. It is shown that the calculated width converges exponentially as the effective principal quantum number of the Rydberg orbital increases. The threshold capture width is then easily obtained. Since atmospheric recombination involves very low energy electrons, the threshold capture widths are essential to the calculation of DR cross sections for the atmospheric species studied here. The approach described makes use of bound state computer codes already in use. A program that collects width matrix elements over CI wavefunctions for the initial and final states is described.

Guberman, S. L.↗

DOE OSTI2021

Gleipnir: toward practical error analysis for Quantum programs

Practical error analysis is essential for the design, optimization, and evaluation of Noisy Intermediate-Scale Quantum(NISQ) computing. However, bounding errors in quantum programs is a grand challenge, because the effects of quantum errors depend on exponentially large quantum states. In this work, we present Gleipnir, a novel methodology toward practically computing verified error bounds in quantum programs. Gleipnir introduces the (ρ,δ)-diamond norm, an error metric constrained by a quantum predicate consisting of the approximate state ρ and its distance δ to the ideal state ρ. This predicate (ρ,δ) can be computed adaptively using tensor networks based on the Matrix Product States. Gleipnir features a lightweight logic for reasoning about error bounds in noisy quantum programs, based on the (ρ,δ)-diamond norm metric. Furthermore, our experimental results show that Gleipnir is able to efficiently generate tight error bounds for real-world quantum programs with 10 to 100 qubits, and can be used to evaluate the error mitigation performance of quantum compiler transformations.

Tao, Runzhou↗

DOE OSTI2022

Exponential Water Uptake in Ionomer Membranes Results from Polymer Plasticization

Water content is the most influential factor in the performance of ion-exchange membranes (IEMs), controlling their mechanical rigidity, ionic conductivity, and degradation rates. Membranes absorb water exponentially at high water activity (a w ), making that region of the sorption isotherm the most influential on membrane properties. Plasticization of the polymer by water has been proposed to cause this exponential uptake at high a w . However, an integrated microscopic picture of the structure, thermodynamic, and mechanical properties of IEMs as a function of a w has remained elusive. In this work, we use large-scale molecular simulations validated with experimental measurements to compute the sorption isotherms, Young’s modulus, polymer dynamics, and structure of IEMs. The simulations unveil that the exponential increase in water uptake (WU) coincides in all cases with the glass to rubber transition of the membrane, as measured through its Young’s modulus and segmental polymer dynamics. Functionalization of the polymer with alkyl groups further contributes to its plasticization, increasing the WU at a given a w and ion-exchange capacity (IEC). The alkyl chains act synergistically with water to plasticize the polymer matrix and allow water penetration in the membrane. The simulations reveal that the width of the water channels depends on the ratio λ of water to ions in the polymer but is independent of its IEC. We conclude that differences in the polymer matrix–not the water channels–are responsible for the distinct uptake response of IEMs to the thermodynamic driving force of water activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Multiphase complete exchange: A theoretical analysis

Complete Exchange requires each of N processors to send a unique message to each of the remaining N-1 processors. For a circuit switched hypercube with N = 2(sub d) processors, the Direct and Standard algorithms for Complete Exchange are optimal for very large and very small message sizes, respectively. For intermediate sizes, a hybrid Multiphase algorithm is better. This carries out Direct exchanges on a set of subcubes whose dimensions are a partition of the integer d. The best such algorithm for a given message size m could hitherto only be found by enumerating all partitions of d. The Multiphase algorithm is analyzed assuming a high performance communication network. It is proved that only algorithms corresponding to equipartitions of d (partitions in which the maximum and minimum elements differ by at most 1) can possibly be optimal. The run times of these algorithms plotted against m form a hull of optimality. It is proved that, although there is an exponential number of partitions, (1) the number of faces on this hull is Theta(square root of d), (2) the hull can be found in theta(square root of d) time, and (3) once it has been found, the optimal algorithm for any given m can be found in Theta(log d) time. These results provide a very fast technique for minimizing communication overhead in many important applications, such as matrix transpose, Fast Fourier transform, and ADI.

Bokhari, Shahid H.↗

DOE OSTI2022

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ0, the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ0 relates to a herd immunity threshold (HIT), which is given by 1−1/ℛ0. When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21 January 2020 to 21 June 2020. Our ℛ0 estimates characterize the infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we found that no state had achieved herd immunity as of 20 September 2021.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Laser Disdrometer Quantities (LDQUANTS) and Video Disdrometer Quantities (VDISQUANTS) Value-Added Products Report

Disdrometers are useful instruments that measure the distribution of raindrop sizes (drop size distribution: DSD) and the associated rainfall rates/accumulation as those raindrops fall to the ground. These DSDs are a quantity of interest to members of the modeler and observational communities. However, use of disdrometer data for model evaluation, radar monitoring, or other activities requires careful quality control and processing for key DSD properties of interest (e.g., the number concentration of drops) to ensure appropriate physical (scattering, fall speed) assumptions. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Laser Disdrometer Quantities and Video Disdrometer Quantities Value-Added Product (LDQUANTS/VDISQUANTS VAP) uses standard methods from Tokay et al. (e.g., 2013, 2014) to filter drops with unrealistic fall speeds. Further, it estimates several microphysical/geophysical quantities of parameterized DSDs (gamma or exponential assumption type fitting methods) as in previous disdrometer studies and ARM long-term efforts (e.g., Testud et al. 2001, Giangrande et al. 2014, Thompson et al. 2015). Disdrometers are also beneficial for cross-checks with other instrumentation, including rain gauges and radars. To support research interests and related radar monitoring activities, this product calculates radar equivalent quantities, including dual polarization radar quantities (e.g., Z, Differential Reflectivity ZDR, etc.), using T-Matrix scattering and additional wavelength, temperature, and drop shape assumptions (e.g., Thurai et al. 2007). All of these efforts allow disdrometer data sets to become more useful and easily handled by modeling and observational studies, or routine instrumentation checks.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Multi-energy reconstructions, central electron temperature measurements, and early detection of the birth and growth of runaway electrons using a versatile soft x-ray pinhole camera at MST

A multi-energy soft x-ray pinhole camera has been designed, built, and deployed at the Madison Symmetric Torus to aid the study of particle and thermal transport, as well as MHD stability physics. This novel imaging diagnostic technique employs a pixelated x-ray detector in which the lower energy threshold for photon detection can be adjusted independently on each pixel. The detector of choice is a PILATUS3 100 K with a 450 μm thick silicon sensor and nearly 100 000 pixels sensitive to photon energies between 1.6 and 30 keV. An ensemble of cubic spline smoothing functions has been applied to the line-integrated data for each time-frame and energy-range, obtaining a reduced standard-deviation when compared to that dominated by photon-noise. The multi-energy local emissivity profiles are obtained from a 1D matrix-based Abel-inversion procedure. Central values of Te can be obtained by modeling the slope of the continuum radiation from ratios of the inverted radial emissivity profiles over multiple energy ranges with no a priori assumptions of plasma profiles, magnetic field reconstruction constraints, high-density limitations, or need of shot-to-shot reproducibility. In tokamak plasmas, a novel application has recently been tested for early detection, 1D imaging, and study of the birth, exponential growth, and saturation of runaway electrons at energies comparable to 100 × Te,0; thus, early results are also presented.

Delgado-Aparicio, L. F. (ORCID:0000000167394380)↗

DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Celestial amplitudes from UV to IR

Celestial amplitudes represent 4D scattering of particles in boost, rather than the usual energy-momentum, eigenstates and hence are sensitive to both UV and IR physics. We show that known UV and IR properties of quantum gravity translate into powerful constraints on the analytic structure of celestial amplitudes. For example the soft UV behavior of quantum gravity is shown to imply that the exact four-particle scattering amplitude is meromorphic in the complex boost weight plane with poles confined to even integers on the negative real axis. Would-be poles on the positive real axis from UV asymptotics are shown to be erased by a flat space analog of the AdS resolution of the bulk point singularity. The residues of the poles on the negative axis are identified with operator coefficients in the IR effective action. Far along the real positive axis, the scattering is argued to grow exponentially according to the black hole area law. Exclusive amplitudes are shown to simply factorize into conformally hard and conformally soft factors. The soft factor contains all IR divergences and is given by a celestial current algebra correlator of Goldstone bosons from spontaneously broken asymptotic symmetries. The hard factor describes the scattering of hard particles together with the boost-eigenstate clouds of soft photons or gravitons required by asymptotic symmetries. These provide an IR safe S-matrix for the scattering of hard particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1997

Algorithms for Performance, Dependability, and Performability Evaluation using Stochastic Activity Networks

Modeling tools and technologies are important for aerospace development. At the University of Illinois, we have worked on advancing the state of the art in modeling by Markov reward models in two important areas: reducing the memory necessary to numerically solve systems represented as stochastic activity networks and other stochastic Petri net extensions while still obtaining solutions in a reasonable amount of time, and finding numerically stable and memory-efficient methods to solve for the reward accumulated during a finite mission time. A long standing problem when modeling with high level formalisms such as stochastic activity networks is the so-called state space explosion, where the number of states increases exponentially with size of the high level model. Thus, the corresponding Markov model becomes prohibitively large and solution is constrained by the the size of primary memory. To reduce the memory necessary to numerically solve complex systems, we propose new methods that can tolerate such large state spaces that do not require any special structure in the model (as many other techniques do). First, we develop methods that generate row and columns of the state transition-rate-matrix on-the-fly, eliminating the need to explicitly store the matrix at all. Next, we introduce a new iterative solution method, called modified adaptive Gauss-Seidel, that exhibits locality in its use of data from the state transition-rate-matrix, permitting us to cache portions of the matrix and hence reduce the solution time. Finally, we develop a new memory and computationally efficient technique for Gauss-Seidel based solvers that avoids the need for generating rows of A in order to solve Ax = b. This is a significant performance improvement for on-the-fly methods as well as other recent solution techniques based on Kronecker operators. Taken together, these new results show that one can solve very large models without any special structure.

Deavours, Daniel D.↗

DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗