Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Derived Matrices”

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

Refinements for Bragg coherent X-ray diffraction imaging: electron backscatter diffraction alignment and strain field computation

Bragg coherent X-ray diffraction imaging (BCDI) allows the 3D measurement of lattice strain along the scattering vector for specific microcrystals. If at least three linearly independent reflections are measured, the 3D variation of the full lattice strain tensor within the microcrystal can be recovered. However, this requires knowledge of the crystal orientation, which is typically attained via estimates based on crystal geometry or synchrotron microbeam Laue diffraction measurements. Presented here is an alternative method to determine the crystal orientation for BCDI measurements using electron backscatter diffraction (EBSD) to align Fe–Ni and Co–Fe alloy microcrystals on three different substrates. The orientation matrix is calculated from EBSD Euler angles and compared with the orientation determined using microbeam Laue diffraction. The average angular mismatch between the orientation matrices is less than ∼6°, which is reasonable for the search for Bragg reflections. The use of an orientation matrix derived from EBSD is demonstrated to align and measure five reflections for a single Fe–Ni microcrystal via multi-reflection BCDI. Using this data set, a refined strain field computation based on the gradient of the complex exponential of the phase is developed. This approach is shown to increase accuracy, especially in the presence of dislocations. The results demonstrate the feasibility of using EBSD to pre-align BCDI samples and the application of more efficient approaches to determine the full lattice strain tensor with greater accuracy.

42 ENGINEERING↗
DOE OSTI2023

Examining the justification for the introduction of a fermion localization function

Becke and Edgecombe suggested in 1990 a theoretical tool to describe electron localization in atoms and molecules, an idea which was borrowed by a large number of nuclear theorists since 2011 to describe nucleon localization in nuclear systems. I argue here that these arguments are highly questionable and cannot be used in interacting systems where effects beyond the naive mean field or the simple Hartree-Fock framework are important and the inclusion of correlations induced by particle interactions is necessary in order to introduce such a localization function. Furthermore, I also describe several aspects of the exchange and irreducible two-body density matrices, which depend on the character and strength of the two-particle interaction and which can be useful in justifying the derivation of an appropriate energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Assessing Gas-Phase Ion Reactivity of 50 Elements with NO and the Direct Application for 239 Pu in Complex Matrices Using ICP-MS/MS

Understanding the reactivity of metal cations with various reaction gases in ICP-MS/MS is important to determine the best gas to use for a given analyte/interference pair. In this study, nitric oxide (NO) was investigated as the reaction gas following previous experimental designs. The reactions with 50 elements were investigated to examine periodic trends in reactivity, validate theoretical modeling of reaction enthalpies as a method to screen reactant gases, and provide a baseline for potential in-line gas separation methods. ICP-MS/MS studies involving actinides are typically limited to Th, U and Pu, with analyses of Np and Am rarely reported in the literature. To date, only two previous methods have investigated the use of NO in ICP-MS/MS analyses. To showcase the utility of NO, a method was developed to measure 239 Pu in the presence of environmental matrix constituent and other actinides, like what could be expected from post-detonation debris, with no chemical separation prior to analysis. 239 Pu + was reacted to form 239 Pu 16 O + , eliminating interferences derived from the sample matrix by measuring 239 Pu + intensity at m/z = 255 ( 239 Pu 16 O + ). To validate NO for 238 U 1 H + interference removal in environmental matrices, standard reference materials were diluted to 1 mg per gram of solution and spiked to 0.05 pg/g of 239 Pu and 1 µg/g 238 U (Pu/U = 5 x 10 -8 ). Measured 239 Pu concentrations were within 6% of the spiked value. In conclusion, these results demonstrate that reliable 239 Pu measurements can be made at levels relevant to nuclear forensics without the need for extensive chemical matrix separation prior to analysis.

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

Geometric conditions for saturating the data processing inequality

Abstract The data processing inequality (DPI) is a scalar inequality satisfied by distinguishability measures on density matrices. For some distinguishability measures, saturation of the scalar DPI implies an operator equation relating the arguments of the measure. These results are typically derived using functional analytic techniques. In a complementary approach, we use geometric techniques to derive a formula that gives an operator equation from DPI saturation for any distinguishability measure; moreover, for a broad class of distinguishability measures, the derived operator equation is sufficient to imply saturation as well. Our operator equation coincides with known results for the sandwiched Rényi relative entropies, and gives new results for α - z Rényi relative entropies and a family of of quantum f -divergences, which we compute explicitly.

Physics↗
DOE OSTI2026

BICEP/Keck. XX. Component-separated Maps of the Polarized Cosmic Microwave Background and Thermal Dust Emission Using Planck and BICEP/Keck Observations through the 2018 Observing Season

We present component-separated polarization maps of the cosmic microwave background (CMB) and Galactic thermal dust emission, derived using data from the BICEP/Keck experiments through the 2018 observing season and Planck. By employing a maximum-likelihood method that utilizes observing matrices, we produce unbiased maps of the CMB and dust signals. We outline the computational challenges and demonstrate an efficient implementation of the component map estimator. We show methods to compute and characterize power spectra of these maps, opening up an alternative way to infer the tensor-to-scalar ratio from our data. We compare the results of this map-based separation method with the baseline BICEP/Keck analysis. Our analysis demonstrates consistency between the two methods, finding an 84% correlation between the pipelines.

cosmic inflation↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On the Derivation of Quasi-Newton Formulas for Optimization in Function Spaces

Newton’s method is usually preferred when solving optimization problems due to its superior convergence properties compared to gradient-based or derivative-free optimization algorithms. However, deriving and computing second-order derivatives needed by Newton’s method often is not trivial and, in some cases, not possible. In such cases quasi-Newton algorithms are a great alternative. In this paper, we provide a new derivation of well-known quasi-Newton formulas in an infinite-dimensional Hilbert space setting. Furthermore, it is known that quasi-Newton update formulas are solutions to certain variational problems over the space of symmetric matrices. In this paper, we formulate similar variational problems over the space of bounded symmetric operators in Hilbert spaces. By changing the constraints of the variational problem we obtain updates (for the Hessian and Hessian inverse) not only for the Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton method but also for Davidon–Fletcher–Powell (DFP), Symmetric Rank One (SR1), and Powell-Symmetric-Broyden (PSB). In addition, for an inverse problem governed by a partial differential equation (PDE), we derive DFP and BFGS “structured” secant formulas that explicitly use the derivative of the regularization and only approximates the second derivative of the misfit term. We show numerical results that demonstrate the desired mesh-independence property and superior performance of the resulting quasi-Newton methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The geometry of the modular bootstrap

Abstract We explore the geometry behind the modular bootstrap and its image in the space of Taylor coefficients of the torus partition function. In the first part, we identify the geometry as an intersection of planes with the convex hull of moment curves onR + ⊗ℤ, with boundaries characterized by the total positivity of generalized Hankel matrices. We phrase the Hankel constraints as a semi-definite program, which has several advantages, such as the validity of bounds irrespective of spin truncation. We derive bounds on the gap, twist-gap, and the space of Taylor coefficients themselves. We find that if the gap is above$$ {\Delta }_{\textrm{gap}}^{\ast } $$ ∆ gap ∗ , where$$ \frac{c-1}{12}<{\Delta}_{\textrm{gap}}^{\ast }<\frac{c}{12} $$ c − 1 12 < Δ gap ∗ < c 12 , all coefficients become bounded on both sides and kinks develop in the space. In the second part, we propose an analytic method of imposing the integrality condition for the degeneracy number in the spinless bootstrap, which leads to a non-convex geometry. We find that even at very low derivative order this condition rules out regions otherwise allowed by bootstraps at high derivative order.

Physics↗
DOE OSTI2022

Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing

Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Classifying topology in photonic crystal slabs with radiative environments

Abstract In the recent years, photonic Chern materials have attracted substantial interest as they feature topological edge states that are robust against disorder, promising to realize defect-agnostic integrated photonic crystal slab devices. However, the out-of-plane radiative losses in those photonic Chern slabs has been previously neglected, yielding limited accuracy for predictions of these systems’ topological protection. Here, we develop a general framework for measuring the topological protection in photonic systems, such as in photonic crystal slabs, while accounting for in-plane and out-of-plane radiative losses. Our approach relies on the spectral localizer that combines the position and Hamiltonian matrices of the system to draw a real-picture of the system’s topology. This operator-based approach to topology allows us to use an effective Hamiltonian directly derived from the full-wave Maxwell equations after discretization via finite-elements method (FEM), resulting in the full account of all the system’s physical processes. As the spectral FEM-localizer is constructed solely from FEM discretization of the system’s master equation, the proposed framework is applicable to any physical system and is compatible with commonly used FEM software. Moving forward, we anticipate the generality of the method to aid in the topological classification of a broad range of complex physical systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

N -representability violations in truncated equation-of-motion coupled-cluster methods

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. Furthermore, these observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

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

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Generating accurate density matrices on the tangent space of a Grassmann manifold

Interpolating a density matrix from a set of known density matrices is not a trivial task. This is because a linear combination of density matrices does not necessarily correspond to another density matrix. In this Communication, density matrices are examined as objects of a Grassmann manifold. Although this manifold is not a vector space, its tangent space is a vector space. As a result, one can map the density matrices on this manifold to their corresponding vectors in the tangent space and then perform interpolations on that tangent space. The resulting interpolated vector can be mapped back to the Grassmann manifold, which can then be utilized (1) as an optimal initial guess for a self-consistent field (SCF) calculation or (2) to derive energy directly without time-consuming SCF iterations. Such a promising approach is denoted as Grassmann interpolation (G-Int). The hydrogen molecule has been used to illustrate that the described interpolated method in this work preserves the essential attributes of a density matrix. For phosphorus mononitride and ferrocene, it was demonstrated numerically that reference points for the definition of the corresponding tangent spaces can be chosen arbitrarily. In addition, the interpolated density matrices provide a superior and essentially converged initial guess for an SCF calculation to make the SCF procedure itself unnecessary. Finally, this accurate, efficient, robust, and systematically improved G-Int strategy has been used for the first time to generate highly accurate potential energy surfaces with fine details for the difficult case, ferrocene.

Chemistry↗
DOE OSTI2026

Carbamoylation as an Effective Tool in the Analysis of the Soman Nerve Agent Marker Pinacolyl Alcohol in Soil Matrices by EI-GC-MS and LC-HRMS

Pinacolyl alcohol (PA) is a Schedule 2 chemical commonly featured in most proficiency tests (PTs) administered by the Organisation for the Prohibition of Chemical Weapons (OPCW) due to its direct as well as forensic link to the nerve agent Soman. Therefore, its detection by Chemical Weapons Convention (CWC) inspection teams during on-site investigations is a strong indicator of the past or latent presence of Soman in the environment. Its small molecular weight (102), early elution time, and poor ionization profile make PA a challenging analyte to detect particularly at low concentrations (∼1–10 μg/g). In this work, 1,1′-carbonyldimidazole (CDI) has been used to effectively modify PA for the first time in two different soil matrices (Virginia type A soil and silt sediment) at two separate concentrations (1 and 10 μg/g) for its subsequent detection by EI-GC-MS and LC-HRMS methods. For the EI-GC-MS analysis, the PA carbamate derivative (PIC) exhibits improved chromatography relative to PA such as improved peak shape, increased molecular weight (196 for PIC and 102 for PA) and increased retention time (16.9 min for PIC and ∼4.1 min for PA). In addition, the derivatization also improves the detection of PA by LC-HRMS as the PIC product possesses protonation sites (i.e., imidazole ring) relative to none exhibited by PA. More importantly, the carbamoylation proceeds under mild conditions (55 °C, no base) and rapidly (3 h), characteristics that make it an appealing protocol for the analysis of PA during OPCW PTs or real case scenarios particularly in instances where it is present at low concentrations. It is anticipated that the protocol can be applied to the forensic analysis of this important Soman marker in various environmental matrices.

Chemistry↗
DOE OSTI2021

Compact representations of structured BFGS matrices

For general large-scale optimization problems compact representations exist in which recursive quasi-Newton update formulas are represented as compact matrix factorizations. For problems in which the objective function contains additional structure, recent structured quasi-Newton methods exploit available second-derivative information and approximate unavailable second derivatives. Here, this article develops the compact representations of two structured Broyden-Fletcher-Goldfarb-Shanno update formulas. The compact representations enable efficient limited memory and initialization strategies. Two limited memory line search algorithms are described for which extensive numerical results demonstrate the efficacy of the algorithms, including comparisons to IPOPT on large machine learning problems, and to L-BFGS on a real world large scale ptychographic imaging application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dispersive sum rules in AdS 2

Dispersion relations for S-matrices and CFT correlators translate UV consistency into bounds on IR observables. In this note, we construct dispersive sum rules for 1D CFTs. We use them to prove bounds on higher-derivative couplings in weakly-coupled non-gravitational EFTs in AdS2. At the leading order in the bulk-point limit, the bounds agree with the flat-space result. We compute the leading universal effect of finite AdS radius on the bounds. Along the way, we give an explicit formula for anomalous dimensions in general higher-derivative contact Witten diagrams in AdS 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗