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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Analytical Failure Prediction Method Developed for Woven and Braided Composites

Historically, advances in aerospace engine performance and durability have been linked to improvements in materials. Recent developments in ceramic matrix composites (CMCs) have led to increased interest in CMCs to achieve revolutionary gains in engine performance. The use of CMCs promises many advantages for advanced turbomachinery engine development and may be especially beneficial for aerospace engines. The most beneficial aspects of CMC material may be its ability to maintain its strength to over 2500 F, its internal material damping, and its relatively low density. Ceramic matrix composites reinforced with two-dimensional woven and braided fabric preforms are being considered for NASA s next-generation reusable rocket turbomachinery applications (for example, see the preceding figure). However, the architecture of a textile composite is complex, and therefore, the parameters controlling its strength properties are numerous. This necessitates the development of engineering approaches that combine analytical methods with limited testing to provide effective, validated design analyses for the textile composite structures development.

Min, James B.↗
DOE OSTI2019

25 Years of Organic Electro-Optics and Future Prospects

Multi-scale modeling methods are employed to assist in the design of improved organic electro-optic (OEO) materials, including simulating the effects of the interaction of such materials with the electrodes of silicon organic hybrid (SOH) and plasmonic organic hybrid (POH) nanoscopic waveguide devices. Theoretical methods have promoted the evolution of OEO materials from low number density chromophore-polymer guest-host composites and chromophores covalently attached to traditional polymer matrices to custom designed high chromophore number density macromolecular materials with controlled intermolecular interactions and reduced matrix dimensionality. New OEO materials exhibit electrooptic activity of 500-600 pm/V in electrically-poled μm thick films and 90-350 pm/V in nm thick films in SOH and POH devices. The improvement in EO activity together with concentration of RF and optical fields in SOH and POH devices has resulted in a substantial improvement in device performance, including reducing voltage-length performance to 40 V-μm, energy efficiency on the order of a femtojoule/bit, device footprints on the order of 1 μm 2 and bandwidths of > 170 GHz for POH devices. Furthermore significant advances with respect to device insertion loss have also been achieved.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Quantification of Porosity in Ceramic Matrix Composites Using Thermography

Thermography is one of the Non-Destructive Evaluation methods used to qualitatively estimate the amount of porosity in as-manufactured ceramic matrix composite (CMC) parts. Through-transmission (TT) thermography allows for comparison of thermal diffusivity caused by differences in material density and porosity at different locations in the part. In this study, volume fraction of pores and its distribution in CMC samples containing different amounts of porosity was quantified using optical microscopy and X-ray computed tomography scans. TT thermography was done on the same specimens and images were analyzed to determine the distribution of thermal diffusivity. Comparison of porosity measurements and thermal diffusivity images shows an inverse correlation. A simple two-phase analytical model describing this relationship is derived for composites and its solution is shown to compare favorably with experimental data. In the absence of porosity data, the model can be used to estimate porosity directly from thermal diffusivity data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗
DOE OSTI2021

Single-Step Direct Laser Writing of Multimetal Oxygen Evolution Catalysts from Liquid Precursors

We investigate a laser direct-write method to synthesize and deposit metastable, mixed transition metal oxides and evaluate their performance as oxygen evolution reaction catalysts. This laser processing method enabled the rapid synthesis of diverse heterogeneous alloy and oxide catalysts directly from cost-effective solution precursors, including catalysts with a high density of nanocrystalline metal alloy inclusions within an amorphous oxide matrix. The nanoscale heterogenous structures of the synthesized catalysts were consistent with reactive force-field Monte Carlo calculations. By evaluating the impact of varying transition metal oxide composition ratios, we created a stable Fe 0.63 Co 0.19 Ni 0.18 O x /C catalyst with a Tafel slope of 38.23 mV dec -1 and overpotential of 247 mV, a performance similar to that of IrO 2 . Synthesized Fe 0.63 Co 0.19 Ni 0.18 O x /C and Fe 0.14 Co 0.46 Ni 0.40 O x /C catalysts were experimentally compared in terms of catalytic performance and structural characteristics to determine that higher iron content and a less crystalline structure in the secondary matrix decreases the charge transfer resistance and thus is beneficial for electrocatalytic activity. Furthermore, this conclusion is supported by density-functional theory calculations showing distorted active sites in ternary metal catalysts are key for lowering overpotentials for the oxygen evolution reaction.

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

Chiral spin liquid and quantum phase transition in the triangular-lattice Hofstadter-Hubbard model

Recent advances in moiré engineering motivate the study of lattice models of strongly correlated electrons subjected to substantial orbital magnetic flux. We analyze the triangular-lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette and a density of one electron per site, where a chiral spin liquid phase may exist between weak-coupling integer quantum Hall and strong-coupling 120° antiferromagnetic phases. Here, we use matrix product state methods and analytical arguments to investigate this model compactified to cylinders of finite circumference. We uncover a glide particle-hole symmetry operation which, we argue, is spontaneously broken at the quantum Hall to spin liquid transition on odd-circumference cylinders. We numerically verify the spontaneous symmetry breaking and further demonstrate that this transition is associated with algebraic long-range correlations of various spin-singlet, charge-neutral operators. For even-circumference cylinders, the transition becomes a crossover associated with a large correlation length that grows substantially with circumference. Our findings suggest that in the two-dimensional limit, the transition to a chiral spin liquid phase is continuous and features critical fluctuations of the current.

Divic, Stefan [University of Pennsylvania, Philade↗
DOE OSTI2025

Entanglement Witness for Indistinguishable Electrons Using Solid-State Spectroscopy

Characterizing entanglement in quantum materials is crucial for advancing next-generation quantum technologies. Despite recent strides in witnessing entanglement in magnetic materials with distinguishable spin modes, quantifying entanglement in systems formed by indistinguishable electrons remains a formidable challenge. To solve this problem, we introduce a method to extract various four-fermion correlations by analyzing the nonlinearity in resonant inelastic x-ray scattering spectra. These correlations constitute the primary components of the cumulant two-particle reduced density matrix. We further derive bounds for its eigenvalues and demonstrate the linear scaling with fermionic entanglement depth, providing a reliable witness for entanglement. Using the material-relevant strongly correlated models as examples, we show how this entanglement witness can efficiently quantify multipartite entanglement across different phase regions, highlighting its advantage over quantum Fisher information. Published by the American Physical Society 2025

Liu, Tongtong (ORCID:0000000295324061)↗
DOE OSTI2015

Canonical-ensemble state-averaged complete active space self-consistent field (SA-CASSCF) strategy for problems with more diabatic than adiabatic states: Charge-bond resonance in monomethine cyanines

This paper reviews basic results from a theory of the a priori classical probabilities (weights) in state-averaged complete active space self-consistent field (SA-CASSCF) models. It addresses how the classical probabilities limit the invariance of the self-consistency condition to transformations of the complete active space configuration interaction (CAS-CI) problem. Such transformations are of interest for choosing representations of the SA-CASSCF solution that are diabatic with respect to some interaction. I achieve the known result that a SA-CASSCF can be self-consistently transformed only within degenerate subspaces of the CAS-CI ensemble density matrix. For uniformly distributed (“microcanonical”) SA-CASSCF ensembles, self-consistency is invariant to any unitary CAS-CI transformation that acts locally on the ensemble support. Most SA-CASSCF applications in current literature are microcanonical. A problem with microcanonical SA-CASSCF models for problems with “more diabatic than adiabatic” states is described. The problem is that not all diabatic energies and couplings are self-consistently resolvable. A canonical-ensemble SA-CASSCF strategy is proposed to solve the problem. For canonical-ensemble SA-CASSCF, the equilibrated ensemble is a Boltzmann density matrix parametrized by its own CAS-CI Hamiltonian and a Lagrange multiplier acting as an inverse “temperature,” unrelated to the physical temperature. Like the convergence criterion for microcanonical-ensemble SA-CASSCF, the equilibration condition for canonical-ensemble SA-CASSCF is invariant to transformations that act locally on the ensemble CAS-CI density matrix. The advantage of a canonical-ensemble description is that more adiabatic states can be included in the support of the ensemble without running into convergence problems. The constraint on the dimensionality of the problem is relieved by the introduction of an energy constraint. The method is illustrated with a complete active space valence-bond (CASVB) analysis of the charge/bond resonance electronic structure of a monomethine cyanine: Michler’s hydrol blue. The diabatic CASVB representation is shown to vary weakly for “temperatures” corresponding to visible photon energies. Canonical-ensemble SA-CASSCF enables the resolution of energies and couplings for all covalent and ionic CASVB structures contributing to the SA-CASSCF ensemble. The CASVB solution describes resonance of charge- and bond-localized electronic structures interacting via bridge resonance superexchange. The resonance couplings can be separated into channels associated with either covalent charge delocalization or chemical bonding interactions, with the latter significantly stronger than the former.

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

Electron dynamics in complex environments with real-time time dependent density functional theory in a QM-MM framework

This article presents a time dependent density functional theory (TDDFT) implementation to propagate the Kohn-Sham equations in real time, including the effects of a molecular environment through a Quantum-Mechanics Molecular-Mechanics (QM-MM) hamiltonian. The code delivers an all-electron description employing Gaussian basis functions, and incorporates the Amber force-field in the QM-MM treatment. The most expensive parts of the computation, comprising the commutators between the hamiltonian and the density matrix—required to propagate the electron dynamics—, and the evaluation of the exchange-correlation energy, were migrated to the CUDA platform to run on graphics processing units, which remarkably accelerates the performance of the code. The method was validated by reproducing linear-response TDDFT results for the absorption spectra of several molecular species. Two different schemes were tested to propagate the quantum dynamics: (i) a leap-frog Verlet algorithm, and (ii) the Magnus expansion to first-order. These two approaches were confronted, to find that the Magnus scheme is more efficient by a factor of six in small molecules. Interestingly, the presence of iron was found to seriously limitate the length of the integration time step, due to the high frequencies associated with the core-electrons. This highlights the importance of pseudopotentials to alleviate the cost of the propagation of the inner states when heavy nuclei are present. Finally, the methodology was applied to investigate the shifts induced by the chemical environment on the most intense UV absorption bands of two model systems of general relevance: the formamide molecule in water solution, and the carboxy-heme group in Flavohemoglobin. In both cases, shifts of several nanometers are observed, consistently with the available experimental data.

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

Sci—Thur AM: YIS - 06: An EPID-based 3D patient dose verification method for SBRT-VMAT delivery

Purpose: Stereotactic body radiation therapy (SBRT) delivered via volumetric modulated arc therapy (VMAT) can strongly benefit from an in vivo patient dose verification due to the large doses per fraction. Electronic portal imaging devices (EPIDs) can be utilized as a patient dose dosimeter. In this work we present a physics-based model which utilizes on-treatment EPID images to reconstruct the dose delivered to an anthropomorphic phantom during SBRT-VMAT delivery. Methods: An SBRT linac beam was modeled using Monte Carlo methods and verified with measured data. Our dose reconstruction model back-projects EPID measured focal fluence upstream of the patient and adds a predicted extra-focal fluence component. This fluence is forward projected onto the patient's density matrix and convolved with dose kernels to calculate dose. The model was validated for two prostate, three lung, and two spine SBRT-VMAT treatments. Results were compared to the treatment planning system's calculation. Results: 2%/2 mm chi comparison calculations gave pass rates for the whole volume, infield, and high dose region respectively, and no lower than: 98%, 95%, 99% for the prostate plans, 99%, 92%, 85% for the lung plans, and 91%, 85%, 81% for the spine plans. A 3%/3mm calculation gave pass rates no lower than 99%, 94%, and 90% for all dose regions for the prostate, lung, and spine respectively. Conclusions: We have developed a physics-based model which calculates delivered dose to phantom (or patient) for SBRT-VMAT delivery using on treatment EPID images. The accuracy of the results has allowed us to test this model clinically.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Complex Langevin simulation of a random matrix model at nonzero chemical potential

In this study we test the complex Langevin algorithm for numerical simulations of a random matrix model of QCD with a first order phase transition to a phase of finite baryon density. We observe that a naive implementation of the algorithm leads to phase quenched results, which were also derived analytically in this article. We test several fixes for the convergence issues of the algorithm, in particular the method of gauge cooling, the shifted representation, the deformation technique and reweighted complex Langevin, but only the latter method reproduces the correct analytical results in the region where the quark mass is inside the domain of the eigenvalues. In order to shed more light on the issues of the methods we also apply them to a similar random matrix model with a milder sign problem and no phase transition, and in that case gauge cooling solves the convergence problems as was shown before in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Waste Characterization Using Gamma Ray Spectrometry with Automated Efficiency Optimization - 13404

Gamma ray spectrometry using High Purity Germanium (HPGe) detectors is commonly employed in assaying radioactive waste streams from a variety of sources: nuclear power plants, Department of Energy (DOE) laboratories, medical facilities, decontamination and decommissioning activities etc. The radioactive material is typically packaged in boxes or drums (for e.g. B-25 boxes or 208 liter drums) and assayed to identify and quantify radionuclides. Depending on the origin of the waste stream, the radionuclides could be special nuclear materials (SNM), fission products, or activation products. Efficiency calibration of the measurement geometry is a critical step in the achieving accurate quantification of radionuclide content. Due to the large size of the waste items, it is impractical and expensive to manufacture gamma ray standard sources for performing a measurement based calibration. For well over a decade, mathematical efficiency methods such as those in Canberra's In Situ Object Counting System (ISOCS) have been successfully employed in the efficiency calibration of gamma based waste assay systems. In the traditional ISOCS based calibrations, the user provides input data such as the dimensions of the waste item, the average density and fill height of the matrix, and matrix composition. As in measurement based calibrations, the user typically defines a homogeneous matrix with a uniform distribution of radioactivity. Actual waste containers can be quite nonuniform, however. Such simplifying assumptions in the efficiency calibration could lead to a large Total Measurement Uncertainty (TMU), thus limiting the amount of waste that can be disposed of as intermediate or low activity level waste. To improve the accuracy of radionuclide quantification, and reduce the TMU, Canberra has developed the capability to optimize the efficiency calibration using the ISOCS method. The optimization is based on benchmarking the efficiency shape and magnitude to the data available in the analyzed gamma ray spectra. Data from measurements of a given item in multiple counting geometries are among the powerful benchmarks that could be used in the optimization. Also, while assaying a waste stream with fission products and activation products emitting gamma lines of multiple energies, optimizing the efficiency on the basis of line activity consistency is very effective. In the present paper, the ISOCS- based optimization methodology is applied to measurement scenarios involving multiple counting geometries, and multi-gamma-line radionuclides. Results will be presented along with accuracy and precision estimates for each measurement. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Mechanical & morphological properties of attapulgite/NR composites: Effect of mixing time variation

The development of composite material based on attapulgite clay (ATP) as a filler and natural rubber (NR) matrices were prepared by combination of melt mixing and latex compounding methods. Sonication technique was chosen in this work to disperse the attapulgite suspension. 6 phr of attapulgite loading was fabricated using different time of mixing ranging from 30 minutes until 2 hours and sonication time was kept constant at 15 minutes. Then, co-coagulating HA latex with attapulgite clay suspension through latex compounding method produced the masterbatch. The masterbatch was compounded with natural rubber by melt mixing method. The mechanical and morphological characteristics were investigated in this work. From mechanical testing, M1 showed the highest value of tensile and tear strength. By comparing with M30 and M2, M1 shows high 300% tensile modulus and lower crosslink density. However, when the time of mixing was prolonged to 2 hours, the results for tensile strength, elongation at break and tear strength were decreased. This is due to flocculation of attapulgite particles. Sonication techniques also proved that the tensile strength and elongation at break of these three samples were higher compared to gum NR (NR) and attapulgite compounded with NR using a conventional method (in-situ 6). From field emission scanning electron microscope (FESEM) results, it revealed that M1 had good dispersion in the NR system. It is proved that the higher tensile strength was due to good dispersion of attapulgite clay in the NR matrix. It was also supported from crosslink density, which is lower than NR and in-situ 6 results. It showed that the penetration of toluene solvent into rubber compound was restricted. The optimum time, M1 give the best results, which can be compared to control the sample.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A universal density matrix functional from molecular orbital-based machine learning: Transferability across organic molecules

We address the degree to which machine learning (ML) can be used to accurately and transferably predict post-Hartree-Fock correlation energies. Refined strategies for feature design and selection are presented, and the molecular-orbital-based machine learning (MOB-ML) method is applied to several test systems. Strikingly, for the second-order Møller-Plessett perturbation theory, coupled cluster with singles and doubles (CCSD), and CCSD with perturbative triples levels of theory, it is shown that the thermally accessible (350 K) potential energy surface for a single water molecule can be described to within 1 mhartree using a model that is trained from only a single reference calculation at a randomized geometry. To explore the breadth of chemical diversity that can be described, MOB-ML is also applied to a new dataset of thermalized (350 K) geometries of 7211 organic models with up to seven heavy atoms. In comparison with the previously reported Δ-ML method, MOB-ML is shown to reach chemical accuracy with threefold fewer training geometries. Finally, a transferability test in which models trained for seven-heavy-atom systems are used to predict energies for thirteen-heavy-atom systems reveals that MOB-ML reaches chemical accuracy with 36-fold fewer training calculations than Δ-ML (140 vs 5000 training calculations).

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

An adaptive Hessian approximated stochastic gradient MCMC method

Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. Here, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation (SA) to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the samples and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix via SA, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and demonstrates better learning results compared to standard SG-MCMC updating rules. The approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. Numerical experiments are performed on several problems, including sampling from 2D correlated distribution, synthetic regression problems, and learning the numerical solutions of heterogeneous elliptic PDE. The numerical results demonstrate great improvement in both the convergence rate and accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Segmentation of RDX and TNT in X‐Ray Computed Tomography Reconstructions of Melt‐Cast Explosives

ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.

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

Precision study of the massive Schwinger model near quantum criticality

We perform a numerical analysis of the massive Schwinger model in the presence of a background electric field. Using the Density Matrix Renormalization Group approach, we efficiently compute the spectrum of the Schwinger model on a staggered lattice with up to 3000 qubits. As a result, we achieve a precise computation of the critical mass of the massive Schwinger model to five digits using four different “criticality criteria,” observing perfect agreement among them Additionally, we discuss the effect of a four-fermion operator deformation of the Schwinger model and compute the critical mass for various values of the deformation parameter.

Critical phenomena↗
NASA NTRS1994

Beginning the 21st century with advanced Automatic Parts Identification (API)

Under the direction of the NASA George C. Marshall Space Flight Center, Huntsville, Alabama, the development and commercialization of an advanced Automated Parts Indentification (API) system is being undertaken by Rockwell International Corporation. The new API system is based on a variable sized, machine-readable, two-dimensioanl matrix symbol that can be applied directly onto most metallic and nonmetallic materials using safe, permanent marking methods. Its checkerboard-like structure is the most space efficient of all symbologies. This high data-density symbology can be applied to products of different material sizes and geometries using application-dependent, computer-driven marking devices. The high fidelity markings produced by these devices can then be captured using a specially designed camera linked to any IBM-compatible computer. Applications of compressed symbology technology will reduce costs and improve quality, productivity, and processes in a wide variety of federal and commercial applications.

Schramm, Fred↗