Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Theory”

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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions for cometary and planetary applications, as predicted by mixed quantum-classical theory

We present new calculations of collision cross sections for state-to-state transitions between the rotational states in an H 2 O + H 2 O system, which are used to generate a new database of collisional rate coefficients for cometary and planetary applications. Calculations were carried out using a mixed quantum-classical theory approach that is implemented in the code MQCT. The large basis set of rotational states used in these calculations permits us to predict thermally averaged cross sections for 441 transitions in para- and ortho-H 2 O in a broad range of temperatures. It is found that all state-to-state transitions in the H 2 O + H 2 O system split into two well-defined groups, one with higher cross-section values and lower energy transfer, which corresponds to the dipole-dipole driven processes. The other group has smaller cross sections and higher energy transfer, driven by higher-order interaction terms. We present a detailed analysis of the theoretical error bars, and we symmetrized the state-to-state transition matrixes to ensure that excitation and quenching processes for each transition satisfy the principle of microscopic reversibility. We also compare our results with other data available from the literature for H 2 O + H 2 O collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Modeling of fission gas diffusion and release for Gd 2 O 3 doped UO 2

Uranium dioxide (UO 2 ) is the primary nuclear fuel in light water reactors, and its excess neutronic reactivity can be controlled by adding burnable absorbers, such as Gd 2 O 3 . This burnable absorber has a large neutron absorption cross-section, lowering the high reactivity of the reactor's initial fuel load. However, there needs to be more understanding of how added Gd 2 O 3 influences the properties of UO 2 under irradiation. To understand the behavior of defects and fission gas in the UO 2 /Gd 2 O 3 system under irradiation, we use cluster dynamics modeling supported by density functional theory calculations. First, we calculate the formation energies of Gd point and cluster defects, and evaluate the temperature-dependent defect concentrations using the defect formation energies and entropies. We show that Gd is soluble in UO 2 , introducing a negative charge in the system. Using this information, we adapted the cluster dynamics code Centipede to model the influence of Gd on U self-diffusion and Xe diffusion in UO 2 with 10 wt% Gd 2 O 3 . Also, we analyzed the Xe diffusion as a function of Gd 2 O 3 concentration, showing that the Xe diffusivity is decreased, which means that the athermal diffusivity due to electronic stopping persists at higher temperatures. In conclusion, the decrease in Xe diffusion means that more Xe stays in the matrix, decreasing the Xe release, and lowering its influence of fission gas release on the thermomechanical properties of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Entropy production in model colloidal suspensions under shear via the fluctuation theorem

Dissipative systems often exhibit novel and unexpected properties. This is, for instance, the case of simple liquids, which, when subjected to shear and after reaching a steady state, can exhibit a negative entropy production over finite length scales and timescales. This result, among others, is captured by nonequilibrium relations known as fluctuation theorems. Using nonequilibrium molecular dynamics simulations, we examine how, by fine-tuning the properties of the components of a complex fluid, we can steer the nonequilibrium response of the fluid. More specifically, we show how we control the nonequilibrium probability distribution for the shear stress and, in turn, how often states with a negative entropy production can occur. To achieve this, we start by characterizing how the size for the liquid matrix impacts the probability of observing negative entropy states, as well as the timescale over which these can be observed. We then measure how the addition of larger particles to this liquid matrix, i.e., simulating a model colloidal suspension, results in an increase in the occurrence of such states. As a result, this suggests how modifications in the composition of the mixture and in the properties of its components lead to an increase in the probability of observing states of negative entropy production and, thus, for the system to run in reverse.

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

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)—a selected configuration interaction technique for vibrational structure theory—has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm−1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Chemistry↗
DOE OSTI2025

Decoherence dynamics in molecular qubits: Exponential, Gaussian and beyond

In this article, we examine how the structure of system–bath interactions can determine commonly encountered temporal decoherence patterns, such as Gaussian and exponential decay, in molecular and other qubits coupled to a thermal bosonic bath. The analysis, based on a pure dephasing picture that admits analytical treatment, shows that decoherence, in general, is neither purely Gaussian nor exponential but rather the exponential of oscillatory functions, with periods determined by the bath’s frequencies. For initially unentangled qubit-bath states, Gaussian decay is always present at early times. It becomes increasingly dominant with increasing temperature, qubit–bath interaction strength, and bath correlation time. Initial system–bath entanglement that arises due to displacement in the position of the bath states preserves the Gaussian decay. By contrast, strict exponential decay arises only in very specific models that we isolate. However, it becomes dominant for times longer than the bath correlation time or for early times when there is initial entanglement due to momentum displacement of the bath states. For molecular electronic decoherence, the long-time exponential regime plays a limited role as it emerges after most coherence is lost. Thus, the Gaussian decay provides a more suitable (albeit imperfect) model of such decoherence. Furthermore, we discuss the connection between electronic decoherence dynamics and electronic spectroscopic line shape theory, where Gaussian spectral peaks correspond to Gaussian coherence decay and Lorentzian peaks correspond to exponential coherence decay. We find that Gaussian spectral peaks, usually associated with inhomogeneous broadening, can emerge from the entangling unitary system–bath dynamics even when there is no inhomogeneity in the initial conditions.

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

Polarized TMD fragmentation functions for 𝐽/𝜓 production

We calculate the matching, at leading order, of the transverse momentum-dependent fragmentation functions (TMDFFs) for light quarks and gluons fragmenting to a 𝐽/𝜓 onto polarized nonrelativistic QCD (NRQCD) TMDFFs. The NRQCD TMDFFs have an operator product expansion in terms of nonperturbative NRQCD production matrix elements. Using the results we obtain, we make predictions for the light quark fragmentation contribution to the production of polarized 𝐽/𝜓 in semi-inclusive deep inelastic scattering (SIDIS) both for unpolarized and longitudinally polarized beams. These results are an important contribution to polarized 𝐽/𝜓 production in SIDIS and thus are needed for comparison with experiments at the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Opportunistic CP violation

In the electroweak sector of the Standard Model, CP violation arises through a very particular interplay between the three quark generations, as described by the Cabibbo-Kobayashi-Maskawa (CKM) mechanism and the single Jarlskog invariant J 4 . Once generalized to the Standard Model Effective Field Theory (SMEFT), this peculiar pattern gets modified by higher-dimensional operators, whose associated Wilson coefficients are usually split into CP-even and odd parts. However, CP violation at dimension four, i.e., at the lowest order in the EFT expansion, blurs this distinction: any Wilson coefficient can interfere with J 4 and mediate CP violation. In this paper, we study such interferences at first order in the SMEFT expansion, $\mathcal{O}(1/Λ^2)$, and we capture their associated parameter space via a set of 1551 linear CP-odd flavor invariants. This construction describes both new, genuinely CP-violating quantities as well as the interference between J 4 and CP-conserving ones. We call this latter possibility opportunistic CP violation. Relying on an appropriate extension of the matrix rank to Taylor expansions, which we dub Taylor rank, we define a procedure to organize the invariants in terms of their magnitude, so as to retain only the relevant ones at a given precision. We explore how this characterization changes when different assumptions are made on the flavor structure of the SMEFT coefficients. Interestingly, some of the CP-odd invariants turn out to be less suppressed than J 4 , even when they capture opportunistic CPV, demonstrating that CP-violation in the SM, at dimension 4, is accidentally small.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Ising meson spectroscopy on a noisy digital quantum simulator

Abstract Quantum simulation has the potential to be an indispensable technique for the investigation of non-perturbative phenomena in strongly-interacting quantum field theories (QFTs). In the modern quantum era, with Noisy Intermediate Scale Quantum (NISQ) simulators widely available and larger-scale quantum machines on the horizon, it is natural to ask: what non-perturbative QFT problems can be solved with the existing quantum hardware? We show that existing noisy quantum machines can be used to analyze the energy spectrum of several strongly-interacting 1+1D QFTs, which exhibit non-perturbative effects like ‘quark confinement’ and ‘false vacuum decay’. We perform quench experiments on IBM’s quantum simulators to compute the energy spectrum of 1+1D quantum Ising model with a longitudinal field. Our results demonstrate that digital quantum simulation in the NISQ era has the potential to be a viable alternative to numerical techniques such as density matrix renormalization group or the truncated conformal space methods for analyzing QFTs.

Science & Technology - Other Topics↗
NASA NTRS1987

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗
DOE OSTI2026

Optimal Transfer Operators in Algebraic Two-Level Methods for Nonsymmetric and Indefinite Problems

Consider an algebraic two-level method applied to the 𝑛-dimensional linear system 𝐴⁢𝒙 = 𝒃 using fine-space preconditioner (i.e., “relaxation” or “smoother”) 𝑀, with 𝑀 ≈ 𝐴, restriction and interpolation 𝑅 and 𝑃, and algebraic coarse-space operator 𝐴 𝑐 : = 𝑅 ∗ ⁢𝐴⁢𝑃. Then, what are the best possible transfer operators 𝑅 and 𝑃 of a given dimension 𝑛 𝑐 < 𝑛? Brannick et al. [12] showed that when 𝐴 and 𝑀 are Hermitian positive definite (HPD), the optimal interpolation is such that its range contains the 𝑛 𝑐 smallest generalized eigenvectors of the matrix pencil (𝐴, 𝑀). Recently, in Ali et al. [5] we generalized this framework to the non-HPD setting, by considering both right (interpolation) and left (restriction) generalized eigenvectors of (𝐴, 𝑀) and defining corresponding nonsymmetric transfer operators {𝑅#, 𝑃#}. Tight convergence bounds for {𝑅#, 𝑃#} are derived in spectral radius, as well as a proof of pseudo-optimality. Note, {𝑅#, 𝑃#} are typically complex valued, which is not practical for real-valued problems. Here, in this work, we build on [5], first characterizing all inner products in which the coarse-space correction defined by {𝑅#, 𝑃#} is orthogonal. We then develop tight two-level convergence bounds in these norms, and prove that the underlying transfer operators {𝑅#, 𝑃#} are genuinely optimal. As a special case, our theory both recovers and extends the HPD results from [12]. Finally, we show how to construct optimal, real-valued transfer operators in the case of that 𝐴 and 𝑀 are real valued, but are not HPD. Numerical examples arising from a discretized advection-reaction equation, wave-equation, and Stokes equations are used to verify and illustrate the theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Analog Systems for Edge Optimization

Over the past decade, analog computing has the subject of substantial research interest providing a path toward improved computational efficiency in the post-Dennard era. Analog matrix vector multiplication (MVM) accelerators provide a popular approach given the ubiquity of MVM operations in numerous applications. However, historically analog computing systems can struggle with applications requiring high precision due to the inherent susceptibility of these systems to analog non-idealities. Therefore, prior work on analog systems has focused either on applications known to be tolerant of limited precision (e.g., neural network inference), or using expensive techniques to emulate high-precision using many analog MVM operations. In this work, we propose an alternative approach. Motivated by recent advances in inexact nonlinear solvers and optimizers, we explore the potential of co-designing optimization algorithms which can take full advantage of the fundamentally inexact analog MVM operations. To enable these co-designed algorithms we also develop a general mathematical theory of the precision and energy efficiency of analog operations, and a new system architecture for tightly-coupled analog and digital computation. Finally, we examine the applicability of analog computing to a wider class of symmetric positive definite systems and find potential in using analog operations as a sparse approximate inverse preconditioner. With these core innovations, this project provides a path toward effectively implementing optimization algorithms on power-constrained autonomous and semi-autonomous systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Shedding Light on Primary Donors in Photosynthetic Reaction Centers

Chlorophylls (Chl)s exist in a variety of flavors and are ubiquitous in both the energy and electron transfer processes of photosynthesis. The functions they perform often occur on the ultrafast (fs–ns) time scale and until recently, these have been difficult to measure in real time. Further, the complexity of the binding pockets and the resulting protein-matrix effects that alter the respective electronic properties have rendered theoretical modeling of these states difficult. Recent advances in experimental methodology, computational modeling, and emergence of new reaction center (RC) structures have renewed interest in these processes and allowed researchers to elucidate previously ambiguous functions of Chls and related pheophytins. This is complemented by a wealth of experimental data obtained from decades of prior research. Studying the electronic properties of Chl molecules has advanced our understanding of both the nature of the primary charge separation and subsequent electron transfer processes of RCs. In this review, we examine the structures of primary electron donors in Type I and Type II RCs in relation to the vast body of spectroscopic research that has been performed on them to date. Further, we present density functional theory calculations on each oxidized primary donor to study both their electronic properties and our ability to model experimental spectroscopic data. This allows us to directly compare the electronic properties of hetero- and homodimeric RCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Prediction of thermal cycling induced cracking in polymer matrix composites

This report summarizes the work done in the period February 1993 through July 1993 on the 'Prediction of Thermal Cycling Induced Cracking In Polymer Matrix Composites' program. An oral presentation of this work was given to Langley personnel in September of 1993. This document was prepared for archival purposes. Progress studies have been performed on the effects of spatial variations in material strength. Qualitative agreement was found with observed patterns of crack distribution. These results were presented to NASA Langley personnel in November 1992. The analytical methodology developed by Prof. McManus in the summer of 1992 (under an ASEE fellowship) has been generalized. A method for predicting matrix cracking due to decreasing temperatures and/or thermal cycling in all plies of an arbitrary laminate has been implemented as a computer code. The code also predicts changes in properties due to the cracking. Experimental progressive cracking studies on a variety of laminates were carried out at Langley Research Center. Results were correlated to predictions using the new methods. Results were initially mixed. This motivated an exploration of the configuration of cracks within laminates. A crack configuration study was carried out by cutting and/or sanding specimens in order to examine the distribution of cracks within the specimens. These investigations were supplemented by dye-penetrant enhanced X-ray photographs. The behavior of thin plies was found to be different from the behavior of thicker plies (or ply groups) on which existing theories are based. Significant edge effects were also noted, which caused the traditional metric of microcracking (count of cracks on a polished edge) to be very inaccurate in some cases. With edge and configuration taken into account, rough agreement with predictions was achieved. All results to date were reviewed with NASA Langley personnel in September 1993.

Mcmanus, Hugh L.↗
NASA NTRS2016

Scattering of Gaussian Beams by Disordered Particulate Media

A frequently observed characteristic of electromagnetic scattering by a disordered particulate medium is the absence of pronounced speckles in angular patterns of the scattered light. It is known that such diffuse speckle-free scattering patterns can be caused by averaging over randomly changing particle positions and/or over a finite spectral range. To get further insight into the possible physical causes of the absence of speckles, we use the numerically exact superposition T-matrix solver of the Maxwell equations and analyze the scattering of plane-wave and Gaussian beams by representative multi-sphere groups. We show that phase and amplitude variations across an incident Gaussian beam do not serve to extinguish the pronounced speckle pattern typical of plane-wave illumination of a fixed multi-particle group. Averaging over random particle positions and/or over a finite spectral range is still required to generate the classical diffuse speckle-free regime.

radiative transfer theory↗
DOE OSTI2020

The black hole weak gravity conjecture with multiple charges

We study the effect of higher-derivative corrections on asymptotically flat, fourdimensional, dyonic black holes in low-energy models of gravity coupled to N U(1) gauge fields. For large extremal black holes, the leading O (1/Q 2 ) correction to the extremality bound is calculated from the most general low-energy effective action containing operators with up to four derivatives. Motivated by the multi-charge generalization of the Weak Gravity Conjecture, we analyze the necessary kinematic conditions for an asymptotically large extremal black hole to decay into a multi-particle state of extremal black holes. In the large black hole regime, we show that the convex hull condition degenerates to the requirement that a certain quartic form constructed from the Wilson coefficients of the fourderivative effective operators, is everywhere positive. Using on-shell unitarity methods, we show that higher-derivative operators are renormalized at one-loop only if they generate local, on-shell matrix elements that are invariant tensors of the electromagnetic duality group U(N). The one-loop logarithmic running of the four-derivative Wilson coefficients is calculated and shown to imply the positivity of the extremality form at some finite value of Q 2 . This result generalizes an argument recently given by Charles [1], and shows that under the given assumptions the multi-charge Weak Gravity Conjecture is not a Swampland criterion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Doping dependence of the electron-phonon coupling in two families of bilayer superconducting cuprates

While electron-phonon coupling (EPC) is crucial for Cooper pairing in conventional superconductors, its role in high-Tc superconducting cuprates is debated. Here, using resonant inelastic x-ray scattering at the oxygen K edge, we study the EPC in Bi 2 Sr 2 CaCu 2 O 8+δ (Bi2212) and Nd 1+x Ba 2-x Cu 3 O 7-δ (NBCO) at different doping levels ranging from heavily underdoped (p=0.07) to overdoped (p=0.21). We analyze the data with a localized Lang-Firsov model that allows for the coherent excitations of two phonon modes. While electronic band dispersion effects are non-negligible, we are able to perform a study of the relative values of EPC matrix elements in these cuprate families. In the case of NBCO, the choice of the excitation energy allows us to disentangle modes related to the CuO chains and the CuO 2 planes. Combining the results from the two families, we find the EPC strength decreases with doping at q∥=(-0.25, 0) r.l.u., but has a nonmonotonic trend as a function of doping at smaller momenta. This behavior is attributed to the screening effect of charge carriers. We also find that the phonon intensity is enhanced in the vicinity of the charge-density-wave excitations while the extracted EPC strength appears to be less sensitive to their proximity. By performing a comparative study of two cuprate families, we are able to identify general trends in the EPC for the cuprates and provide experimental input to theories invoking a synergistic role for this interaction in d -wave pairing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1998

Effects of Fiber Coating Composition on Mechanical Behavior of Silicon Carbide Fiber-Reinforced Celsian Composites

Celsian matrix composites reinforced with Hi-Nicalon fibers, precoated with a dual layer of BN/SiC by chemical vapor deposition in two separate batches, were fabricated. Mechanical properties of the composites were measured in three-point flexure. Despite supposedly identical processing, the composite panels fabricated with fibers coated in two batches exhibited substantially different mechanical behavior. The first matrix cracking stresses (sigma(sub mc)) of the composites reinforced with fibers coated in batch 1 and batch 2 were 436 and 122 MPa, respectively. This large difference in sigma(sub mc) was attributed to differences in fiber sliding stresses(tau(sub friction)), 121.2+/-48.7 and 10.4+/-3.1 MPa, respectively, for the two composites as determined by the fiber push-in method. Such a large difference in values of tau(sub friction) for the two composites was found to be due to the difference in the compositions of the interface coatings. Scanning Auger microprobe analysis revealed the presence of carbon layers between the fiber and BN, and also between the BN and SiC coatings in the composite showing lower tau(sub friction). This resulted in lower sigma(sub mc) in agreement with the ACK theory. The ultimate strengths of the two composites, 904 and 759 MPa, depended mainly on the fiber volume fraction and were not significantly effected by tau(sub friction) values, as expected. The poor reproducibility of the fiber coating composition between the two batches was judged to be the primary source of the large differences in performance of the two composites.

Bansal, Narottam P.↗

Refine your search to explore more results.