Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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 703 records · Page 39

A structured framework for predicting sustainable aviation fuel properties using liquid-phase FTIR and machine learning

Sustainable aviation fuels have the potential to improve efficiency, reduce emissions, and enhance energy security. To help identify viable sustainable aviation fuels and accelerate research, machine learning models have been developed to predict relevant physicochemical properties. However, many models have limited applicability, leverage data from complex analytical techniques with confined spectral ranges, or use feature decomposition methods that offer limited interpretability. Using liquid-phase Fourier Transform Infrared (FTIR) spectra, this study presents a structured method for creating accurate and interpretable property prediction models for neat molecules, aviation fuels, and blends. Liquid FTIR spectra can be collected quickly and consistently, offering high reliability, sensitivity, and component specificity using less than 2 ml of sample. The method first decomposes FTIR spectra into fundamental building blocks using non-negative matrix factorization (NMF) to enable scientific analysis of FTIR spectra attributes and fuel properties. The NMF features are then used to create five ensemble models for predicting final boiling point, flash point, freezing point, density at 15°C, and kinematic viscosity at -20°C. All models were trained using experimental property data from neat molecules, aviation fuels, and blends. The models accurately predict key properties across a broad range of neat molecules and representative fuels and blends, while enabling interpretation of relationships between compositional elements, such as functional groups or chemical classes, and their resulting properties. This demonstrates strong potential to support sustainable aviation fuel research and development. The models and data are available on an interactive web tool.

Fourier transform infrared spectroscopy↗

Efficient optimal reconstruction of linear fields and band-powers from cosmological data

We present an efficient implementation of Wiener filtering of real-space linear field and optimal quadratic estimator of its power spectrum Band-powers. We first recast the field reconstruction into an optimization problem, which we solve using quasi-Newton optimization. We then recast the power spectrum estimation into the field marginalization problem, from which we obtain an expression that depends on the field reconstruction solution and a determinant term. We develop a novel simulation based method for the latter. We extend the simulations formalism to provide the covariance matrix for the power spectrum. We develop a flexible framework that can be used on a variety of cosmological fields and present results for a variety of test cases, using simulated examples of projected density fields, projected shear maps from galaxy lensing, and observed Cosmic Microwave Background (CMB) temperature anisotropies, with a wide range of map incompleteness and variable noise. For smaller cases where direct numerical inversion is possible, we show that our solution matches that created by direct Wiener Filtering at a fraction of the overall computation cost. Even more significant reduction of computational is achieved by this implementation of optimal quadratic estimator due to the fast evaluation of the Hessian matrix. This technique allows for accurate map and power spectrum reconstruction with complex masks and nontrivial noise properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Entanglement of exact excited states of Affleck-Kennedy-Lieb-Tasaki models: Exact results, many-body scars, and violation of the strong eigenstate thermalization hypothesis

We obtain multiple exact results on the entanglement of the exact excited states of nonintegrable models we introduced in Phys. Rev. B 98, 235155 (2018). We first discuss a general formalism to analytically compute the entanglement spectra of exact excited states using matrix product states and matrix product operators and illustrate the method by reproducing a general result on single-mode excitations. We then apply this technique to analytically obtain the entanglement spectra of the infinite tower of states of the spin-S AKLT models in the zero and finite energy density limits. We show that in the zero energy density limit, the entanglement spectra of the tower of states are multiple shifted copies of the ground-state entanglement spectrum in the thermodynamic limit. We show that such a resemblance is destroyed at any nonzero energy density. Furthermore, the entanglement entropy S of the states of the tower that are in the bulk of the spectrum is subthermal S ∝ log L, as opposed to a volume law S ∝ L, thus indicating a violation of the strong eigenstate thermalization hypothesis (ETH). These states are examples of what are now called many-body scars. Lastly, we analytically study the finite-size effects and symmetry-protected degeneracies in the entanglement spectra of the excited states, extending the existing theory.

1-dimensional spin chains↗

A Study of a Surrogate Nuclear Thermal Propulsion Fuel Element

We report a surrogate fuel, Zirconium Carbide/Titanium Nitride has been prepared using a simple method to obtain a more uniform distribution of titanium nitride in the zirconium carbide matrix. Powders were sintered using Spark Plasma Sintering at 1750oC and 1850oC using a modified sintering profile. This resulted in a fairly uniform distribution of titanium nitride in the zirconium carbide matrix with a density of 95% of theoretical at 1850oC. Scanning electron microscopy, indentation hardness and x-ray diffraction results are given.

36 MATERIALS SCIENCE↗

On differentiating the probability of error in multipopular feature selection

A method of linear feature selection for n dimensional observation vectors which belong to one of m populations is presented. Each population has a known apriori probability and is described by a known multivariate normal density function. Specifically we consider the problem of finding a k x n matrix B of rank k (k n) for which the transformed probability of misclassification is minimized. Providing that the transformed a posterior probabilities are distinct theoretical results are obtained which, for the case k = l, give rise to a numerically tractable formula for the derivative of the probability of misclassification. It is shown that for the two population problem this condition is also necessary. The dependence of the minimum probability of error on the a priori probabilities is investigated. The minimum probability of error satisfies a uniform Lipschitz condition with respect to the a priori probabilities.

Peters, B. C.↗

Matrix cracking in laminated composites under monotonic and cyclic loadings

An analytical model based on the internal state variable (ISV) concept and the strain energy method is proposed for characterizing the monotonic and cyclic response of laminated composites containing matrix cracks. A modified constitution is formulated for angle-ply laminates under general in-plane mechanical loading and constant temperature change. A monotonic matrix cracking criterion is developed for predicting the crack density in cross-ply laminates as a function of the applied laminate axial stress. An initial formulation for a cyclic matrix cracking criterion for cross-ply laminates is also discussed. For the monotonic loading case, a number of experimental data and well-known models are compared with the present study for validating the practical applicability of the ISV approach.

Allen, David H.↗

Single Atomic Iron Site Catalysts via Benign Aqueous Synthesis for Durability Improvement in Proton Exchange Membrane Fuel Cells

Atomically-dispersed iron-nitrogen-carbon (Fe–N–C) catalysts have arisen as promising candidates for replacing the costly precious metal catalysts in fuel cells but still face some grand challenges, such as insufficient site density and durability. Herein, we report a self-assembly method in an aqueous solution to develop an atomically-dispersed iron catalyst with high oxygen reduction reaction (ORR) activity and stability in acidic electrolytes. As determined by high-resolution transmission electron microscopy (HR-TEM), X-ray absorption spectroscopy (XAS), and high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM), this benign aqueous synthesis strategy facilitates the formation of homogeneous atomic nitrogen-coordinated iron sites embedded in a popcorn-like porous graphitic carbon matrix. These catalyst properties contribute to the improved ORR kinetic current density and mass transport. By controlling synthesis chemistry, the correlation between structure and property is systematically investigated. In this work, the iron content is the most critical material property and can regulate site density and graphitic carbon structures in the catalyst, impacting catalytic activity and stability. The enhanced performance and durability were examined in both acidic aqueous electrolytes and membrane electrode assemblies.

25 ENERGY STORAGE↗

Benchmark tests and spin adaptation for the particle-particle random phase approximation

The particle-particle random phase approximation (pp-RPA) provides an approximation to the correlation energy in density functional theory via the adiabatic connection [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013)]. It has virtually no delocalization error nor static correlation error for single-bond systems. However, with its formal O(N{sup 6}) scaling, the pp-RPA is computationally expensive. In this paper, we implement a spin-separated and spin-adapted pp-RPA algorithm, which reduces the computational cost by a substantial factor. We then perform benchmark tests on the G2/97 enthalpies of formation database, DBH24 reaction barrier database, and four test sets for non-bonded interactions (HB6/04, CT7/04, DI6/04, and WI9/04). For the G2/97 database, the pp-RPA gives a significantly smaller mean absolute error (8.3 kcal/mol) than the direct particle-hole RPA (ph-RPA) (22.7 kcal/mol). Furthermore, the error in the pp-RPA is nearly constant with the number of atoms in a molecule, while the error in the ph-RPA increases. For chemical reactions involving typical organic closed-shell molecules, pp- and ph-RPA both give accurate reaction energies. Similarly, both RPAs perform well for reaction barriers and nonbonded interactions. These results suggest that the pp-RPA gives reliable energies in chemical applications. The adiabatic connection formalism based on pairing matrix fluctuation is therefore expected to lead to widely applicable and accurate density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How many invariant polynomials are needed to decide local unitary equivalence of qubit states?

Given L-qubit states with the fixed spectra of reduced one-qubit density matrices, we find a formula for the minimal number of invariant polynomials needed for solving local unitary (LU) equivalence problem, that is, problem of deciding if two states can be connected by local unitary operations. Interestingly, this number is not the same for every collection of the spectra. Some spectra require less polynomials to solve LU equivalence problem than others. The result is obtained using geometric methods, i.e., by calculating the dimensions of reduced spaces, stemming from the symplectic reduction procedure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Relaxation dynamics in quantum dissipative systems: The microscopic effect of intramolecular vibrational energy redistribution

We investigate the effect of inter-mode coupling on the vibrational relaxation dynamics of molecules in weak dissipative environments. The simulations are performed within the reduced density matrix formalism in the Markovian regime, assuming a Lindblad form for the system-bath interaction. The prototypical two-dimensional model system representing two CO molecules approaching a Cu(100) surface is adapted from an ab initio potential, while the diatom-diatom vibrational coupling strength is systematically varied. In the weak system-bath coupling limit and at low temperatures, only first order non-adiabatic uni-modal coupling terms contribute to surface-mediated vibrational relaxation. Since dissipative dynamics is non-unitary, the choice of representation will affect the evolution of the reduced density matrix. Two alternative representations for computing the relaxation rates and the associated operators are thus compared: the fully coupled spectral basis, and a factorizable ansatz. The former is well-established and serves as a benchmark for the solution of Liouville-von Neumann equation. In the latter, a contracted grid basis of potential-optimized discrete variable representation is tailored to incorporate most of the inter-mode coupling, while the Lindblad operators are represented as tensor products of one-dimensional operators, for consistency. This procedure results in a marked reduction of the grid size and in a much more advantageous scaling of the computational cost with respect to the increase of the dimensionality of the system. The factorizable method is found to provide an accurate description of the dissipative quantum dynamics of the model system, specifically of the time evolution of the state populations and of the probability density distribution of the molecular wave packet. The influence of intra-molecular vibrational energy redistribution appears to be properly taken into account by the new model on the whole range of coupling strengths. It demontrates that most of the mode mixing during relaxation is due to the potential part of the Hamiltonian and not to the coupling among relaxation operators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Silicon in U-10Mo Alloy

This document details a method for evaluating the effect of silicon impurity content on U-10Mo alloys. Silicon concentration in U-10Mo alloys has been shown to impact the following: volume fraction of precipitate phases, effective density of the final alloy, and 235-U enrichment in the gamma-UMo matrix. This report presents a model for calculating these quantities as a function of Silicon concentration, which along with fuel foil characterization data, will serve as a reference for quality control of the U-10Mo final alloy Si content. Additionally, detailed characterization using scanning electron microscope imaging, transmission electron microscope diffraction, and atom probe tomography showed that Silicon impurities present in U-10Mo alloys form a Si-rich precipitate phase.

36 MATERIALS SCIENCE↗

An Accurate Quartic Force Field, Fundamental Frequencies, and Binding Energy for the High Energy Density Material T(d)N4

The CCSD(T) method has been used to compute a highly accurate quartic force field and fundamental frequencies for all N-14 and N-15 isotopomers of the high energy density material T(sub d)N(sub 4). The computed fundamental frequencies show beyond doubt that the bands observed in a matrix isolation experiment by Radziszewski and coworkers are not due to different isotopomers of T(sub d)N(sub 4). The most sophisticated thermochemical calculations to date yield a N(sub 4) -> 2N(sub 2) heat of reaction of 182.22 +/- 0.5 kcal/mol at 0 K (180.64 +/- 0.5 at 298 K). It is hoped that the data reported herein will aid in the ultimate detection of T(sub d)N(sub 4).

Lee, Timothy J.↗

Synthesis, structure, growth and characterization of a novel organic NLO single crystal: Morpholin-4-ium p-aminobenzoate

Highlights: ► A new organic NLO crystal morpholin-4-ium p-aminobenzoate has been grown for the first time. ► The structure is reported for the first time in the literature. ► Thermal, optical and SHG studies suggest its suitability for various NLO applications. -- Abstract: The title compound, morpholin-4-ium p-aminobenzoate (MPABA)(C{sub 4}H{sub 10}NO{sup +},C{sub 7}H{sub 6}NO{sub 2}{sup −}), has been synthesized for the first time by the addition of morpholine with 4-aminobenzoic acid in equi-molar ratio and good quality single crystals have been grown by solution growth technique using methanol as a solvent. The molecular structure of the compound was solved and refined by Direct Methods using SHELXS97 and full-matrix least-squares technique using SHELXL97, respectively. MPABA crystallizes in a monoclinic system with unit cell parameters, a = 5.948(5) Å, b = 18.033(4) Å, c = 10.577(5) Å, β = 90.40(1)° and non-centrosymmetric space group Cc. The experimentally measured density and chemical compositions were found to be in good agreement with the theoretical values. The phases and functional groups of MPABA have been identified and confirmed through powder X-ray diffraction and Fourier transform infrared (FTIR) studies, respectively. The thermal stability and decomposition details were studied through TG/DTA thermograms. The UV–visible transmission spectra were recorded for the grown crystal and its NLO characteristic was explored by powder second harmonic generation (SHG) studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exchange-correlation energy from pairing matrix fluctuation and the particle-particle random phase approximation

Despite their unmatched success for many applications, commonly used local, semi-local, and hybrid density functionals still face challenges when it comes to describing long-range interactions, static correlation, and electron delocalization. Density functionals of both the occupied and virtual orbitals are able to address these problems. The particle-hole (ph-) Random Phase Approximation (RPA), a functional of occupied and virtual orbitals, has recently known a revival within the density functional theory community. Following up on an idea introduced in our recent communication [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013)], we formulate more general adiabatic connections for the correlation energy in terms of pairing matrix fluctuations described by the particle-particle (pp-) propagator. With numerical examples of the pp-RPA, the lowest-order approximation to the pp-propagator, we illustrate the potential of density functional approximations based on pairing matrix fluctuations. The pp-RPA is size-extensive, self-interaction free, fully anti-symmetric, describes the strong static correlation limit in H{sub 2}, and eliminates delocalization errors in H{sub 2}{sup +} and other single-bond systems. It gives surprisingly good non-bonded interaction energies – competitive with the ph-RPA – with the correct R{sup −6} asymptotic decay as a function of the separation R, which we argue is mainly attributable to its correct second-order energy term. While the pp-RPA tends to underestimate absolute correlation energies, it gives good relative energies: much better atomization energies than the ph-RPA, as it has no tendency to underbind, and reaction energies of similar quality. The adiabatic connection in terms of pairing matrix fluctuation paves the way for promising new density functional approximations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parallel-in-Time Simulation of Lindblad's Equation

Constructing fast quantum logic gates is critical to building a scalable quantum computer. We consider a qudit, a quantum version of a bit that can take an arbitrary number of states, coupled with a cavity. In this project, we wish to force the qudit to reach the 0-state, for any possible initial state. The coupled system changes in time according to Lindblad’s equation, an ordinary differential equation on the density matrix of the quantum system. Lindblad’s equation contains some parameters that we can control, so-called control functions. We seek control functions which force the qudit to the 0-state within 2 microseconds, which is much faster than what is currently done in practice. The search method is gradient descent, a numerical optimization method that uses gradient information to iteratively improve the control parameters. My contribution to this project is an attempt to speed up the computation of the gradient. It currently takes about 40 seconds to compute the gradient which involves solving a set of ODEs sequentially. Current supercomputers have thousands of cores, but sequential computations can only make use of 1 core at a time. We wish to divide up the work better, so that we can use many more cores at once. To this end, we have implemented the Multigrid Reduction in Time (MGRIT) algorithm. We perform a systematic parameter search on how to best apply this algorithm. Results indicate a 25 percent speed up for solving Lindblad’s equation and determining how close the final state is the 0-state.

97 MATHEMATICS AND COMPUTING↗

Structure of Active Sites of Fe-N-C Nano-Catalysts for Alkaline Exchange Membrane Fuel Cells

Platinum group metal-free (PGM-free) catalysts based on transition metal-nitrogen-carbon nanomaterials have been studied by a combination of ex situ and in situ synchrotron X-ray spectroscopy techniques; high-resolution Transmission Electron Microscope (TEM); Mößbauer spectroscopy combined with electrochemical methods and Density Functional Theory (DFT) modeling/theoretical approaches. The main objective of this study was to correlate the HO2− generation with the chemical nature and surface availability of active sites in iron-nitrogen-carbon (Fe-N-C) catalysts derived by sacrificial support method (SSM). These nanomaterials present a carbonaceous matrix with nitrogen-doped sites and atomically dispersed and; in some cases; iron and nanoparticles embedded in the carbonaceous matrix. Fe-N-C oxygen reduction reaction electrocatalysts were synthesized by varying several synthetic parameters to obtain nanomaterials with different composition and morphology. Combining spectroscopy, microscopy and electrochemical reactivity allowed the building of structure-to-properties correlations which demonstrate the contributions of these moieties to the catalyst activity, and mechanistically assign the active sites to individual reaction steps. Associated with Fe-Nx motive and the presence of Fe metallic particles in the electrocatalysts showed the clear differences in the variation of composition; processing and treatment conditions of SSM. From the results of material characterization; catalytic activity and theoretical studies; Fe metallic particles (coated with carbon) are main contributors into the HO2− generation.

Kishi, Hirofumi↗

Up-and-coming IMCs

While the good oxidation and environmental resistance, high melting points, and comparatively low densities of such ordered intermetallics as Ti3Al, NiAl, FeAl, and NbAl3 render them good candidates for advanced aerospace structures, their poor toughness at low temperatures and low strength at elevated temperatures have prompted the development of fiber-reinforced intermetallic-matrix composites (IMCs) with more balanced characteristics. Fabrication methods for continuous-fiber IMCs under development include the P/M 'powder cloth' method, the foil/fiber method, and thermal spraying. The ultimate success of IMCs depends on fibers truly compatible with the matrix materials.

Bowman, Randy↗

Constructing LDPC Codes from Loop-Free Encoding Modules

A method of constructing certain low-density parity-check (LDPC) codes by use of relatively simple loop-free coding modules has been developed. The subclasses of LDPC codes to which the method applies includes accumulate-repeat-accumulate (ARA) codes, accumulate-repeat-check-accumulate codes, and the codes described in Accumulate-Repeat-Accumulate-Accumulate Codes (NPO-41305), NASA Tech Briefs, Vol. 31, No. 9 (September 2007), page 90. All of the affected codes can be characterized as serial/parallel (hybrid) concatenations of such relatively simple modules as accumulators, repetition codes, differentiators, and punctured single-parity check codes. These are error-correcting codes suitable for use in a variety of wireless data-communication systems that include noisy channels. These codes can also be characterized as hybrid turbolike codes that have projected graph or protograph representations (for example see figure); these characteristics make it possible to design high-speed iterative decoders that utilize belief-propagation algorithms. The present method comprises two related submethods for constructing LDPC codes from simple loop-free modules with circulant permutations. The first submethod is an iterative encoding method based on the erasure-decoding algorithm. The computations required by this method are well organized because they involve a parity-check matrix having a block-circulant structure. The second submethod involves the use of block-circulant generator matrices. The encoders of this method are very similar to those of recursive convolutional codes. Some encoders according to this second submethod have been implemented in a small field-programmable gate array that operates at a speed of 100 megasymbols per second. By use of density evolution (a computational- simulation technique for analyzing performances of LDPC codes), it has been shown through some examples that as the block size goes to infinity, low iterative decoding thresholds close to channel capacity limits can be achieved for the codes of the type in question having low maximum variable node degrees. The decoding thresholds in these examples are lower than those of the best-known unstructured irregular LDPC codes constrained to have the same maximum node degrees. Furthermore, the present method enables the construction of codes of any desired rate with thresholds that stay uniformly close to their respective channel capacity thresholds.

Divsalar, Dariush↗