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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Synthesis and Characterization of Polymeric Pore Former Microspheres using Modified Polymer Resin

Porosity in polymer parts is often desired to improve mechanical properties. Individual pores must be small enough that bulk mechanical properties are retained but significant enough to reduce the density of the material. Filler materials can be incorporated into polymer matrices and subsequently removed to create voids. The materials used to create the voids are known as pore-formers. The present research aims to create spherical pores in a silicone matrix by inserting and removing core materials or fillers post-cure. A method was developed to incorporate pore formers into an uncured polymer resin and remove the materials from the matrix. Polyethylene Glycol (PEG) and Polylactic Acid (PLA) microspheres are created and incorporated into silicone matrices as pore formers. Microspheres with controllable size distributions are formed via ultrasonic spray pyrolysis and emulsion formation. Microspheres are removed from the matrix by solvent extraction and calcination. PLA microspheres were modified by hydroxyl addition to the surface to study their solubility in water. PEG-PLA block copolymers are also synthesized and investigated, comparing their physical and chemical properties to the monomers to identify whether the block copolymer microspheres will be more easily removable than their monomeric microspheres. Moreover, other properties including thermal stability, rheology, and mechanical response are studied. The present dissertation provides a comprehensive understanding of the synthesis and characterization of PEG and PLA-based pore formers, toward the future development of feedstock for additive manufacturing porous structures.

Damptey, Ransford↗

DOE OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗

DOE OSTI2021

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Identifying Bounds of Inorganic Content in Solventless Processing of Hybrid Solid Electrolytes

Solid-state lithium batteries require safe, robust electrolytes to enable higher energy densities and improved safety over conventional cells. Hybrid polymer–ceramic electrolytes are a promising solution, combining the processability of polymers with the high ionic conductivity and mechanical strength of inorganic fillers. In this work, we demonstrate a solventless, UV-curing method to produce hybrid solid electrolytes using a poly(ethylene glycol) dimethyl ether (PEGDME)-based photocurable matrix incorporating Li 1.5 Al 0.5 Ge 1.5 (PO 4 ) 3 (LAGP) or Li 7 La 3 Zr 2 O 12 (LLZO) ceramic electrolyte. Inorganic filler loadings up to ∼55 wt.% could be successfully incorporated via this process which was the highest inorganic content at which the slurry remains processable and cured into a uniform film. The resulting UV-cured composite electrolytes remain flexible and exhibit room-temperature ionic conductivities on the order of 10 −4 S·cm −1 , along with notably improved lithium-ion transference numbers compared to conventional polymer electrolytes. Similar performance and processing limits were observed for both LAGP and LLZO, indicating that ceramic filler chemistry does not significantly affect the UV-curing process or the electrolyte's ion transport properties in this regime. Eliminating solvents from fabrication not only simplifies processing and mitigates environmental concerns but also enables higher solid contents that enhance mechanical strength and help suppress lithium dendrite formation. In conclusion, this scalable approach thus paves the way for manufacturing robust composite solid electrolytes for next-generation solid-state batteries (SSBs).

Batteries↗

DOE OSTI2024

Separation of Lipoproteins for Quantitative Analysis of 14 C-Labeled Lipid-Soluble Compounds by Accelerator Mass Spectrometry

To date, 14 C tracer studies using accelerator mass spectrometry (AMS) have not yet resolved lipid-soluble analytes into individual lipoprotein density subclasses. The objective of this work was to develop a reliable method for lipoprotein separation and quantitative recovery for biokinetic modeling purposes. The novel method developed provides the means for use of small volumes (10–200 µL) of frozen plasma as a starting material for continuous isopycnic lipoprotein separation within a carbon- and pH-stable analyte matrix, which, following post-separation fraction clean up, created samples suitable for highly accurate 14 C/ 12 C isotope ratio determinations by AMS. Manual aspiration achieved 99.2 ± 0.41% recovery of [5- 14 CH 3 ]-(2R, 4'R, 8'R)-α-tocopherol contained within 25 µL plasma recovered in triacylglycerol rich lipoproteins (TRL = Chylomicrons + VLDL), LDL, HDL, and infranatant (INF) from each of 10 different sampling times for one male and one female subject, n = 20 total samples. Small sample volumes of previously frozen plasma and high analyte recoveries make this an attractive method for AMS studies using newer, smaller footprint AMS equipment to develop genuine tracer analyses of lipophilic nutrients or compounds in all human age ranges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗

DOE OSTI2023

Constructing a Simulation Surrogate with Partially Observed Output

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. In conclusion, the proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

42 ENGINEERING↗

DOE OSTI2020

Impact of raw material surface oxide removal on dual band infrared optical properties of As 2 Se 3 chalcogenide glass

The manufacturing of low loss chalcogenide glasses (ChGs) for optoelectronic applications is ultimately defined by the concentration of impurities present in starting materials or imparted via processing. We describe a rapid method for purifying metallic starting materials in As 2 Se 3 glass where oxide reduction is correlated to optical and physical properties. Specifically, As-O reduction enhances the glass’ dual-band optical transparency proportional to the extent (13-fold reduction) of oxide reduction, and is accompanied by a change in density and hardness associated with changes in matrix bonding. A significant modification of the glass’ index and LWIR Abbe number is reported highlighting the significant impact purification has on material dispersion control required in optical designs.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

DOE OSTI2023

Witnessing Nonequilibrium Entanglement Dynamics in a Strongly Correlated Fermionic Chain

Many-body entanglement in condensed matter systems can be diagnosed from equilibrium response functions through the use of entanglement witnesses and operator-specific quantum bounds. Here, we investigate the applicability of this approach for detecting entangled states in quantum systems driven out of equilibrium. We use a multipartite entanglement witness, the quantum Fisher information, to study the dynamics of a paradigmatic fermion chain undergoing a time-dependent change of the Coulomb interaction. Our results show that the quantum Fisher information is able to witness distinct signatures of multipartite entanglement both near and far from equilibrium that are robust against decoherence. Here, we discuss implications of these findings for probing entanglement in light-driven quantum materials with time-resolved optical and x-ray scattering methods.

1-dimensional spin chains↗

DOE OSTI2022

Scalable, Infiltration-Free Ceramic Matrix Composite (SIF-CMC) manufacturing for Molten Salt Receiver: Investigation of particle-loaded preceramic resin

Due to their unique combination of high temperature, corrosion resistance, and toughness, Ceramic Matrix Composites (CMC) are the optimal materials to be used in the next generation of Concentrated Solar Power (CSP), which will operate at working temperatures up to 800 °C. However, today, CMC is too costly to be used in CSP applications due to its complex manufacturing process based on repeated infiltration-pyrolysis cycles. To solve this issue, PARC is developing a scalable, infiltration-free method (SIF) to manufacture CMC, which can reduce the production cost, increase batch consistency, and leverage all of the intrinsic advantages of current CMCs. The key to the SIF method is using a flowable preceramic resin with high loading of functional components. The formulation impregnates the fiber during curing and leads to high-density, low-shrinkage parts upon pyrolysis. In this paper, we demonstrate a proof-of-concept of the SIF process and report the relationships between key process parameters in the manufacturing process: functional material loading in the resin vs. resin’s flowability and morphology changes upon pyrolysis.

ceramic matrix composite, concentrated solar power↗

DOE OSTI2022

Design of a multi-functional gel polymer electrolyte with a 3D compact stacked polymer micro-sphere matrix for high-performance lithium metal batteries

We report lithium metal batteries (LMBs) are considered as promising high energy density batteries. However, they still suffer from poor cyclability due to the instability of solid electrolyte interphases (SEIs) and Li dendrite issues. Herein, a novel 3D densely packed multifunctional crosslinked gel polymer electrolyte (PHGPE) has been synthesized by in situ copolymerization of pentaerythritol tetraacrylate and 2-hydroxyethyl acrylate in a liquid electrolyte via a thermal initiation method. Due to intermolecular hydrogen bonding in the frameworks, the PHGPE exhibits a 3D compact-stacked and crosslinked micro-sphere structure, which benefits the construction of homogeneous Li + flux, enabling uniform Li plating/stripping and formation of robust SEIs onto the lithium electrode. Benefiting from the above advantages, symmetric Li cells in the PHGPE deliver excellent cyclabilities of 6000 h at a current density of 1 mA cm -2 (1 mA h cm -2 ) and 7700 h at 0.5 mA cm -2 (2 mA h cm -2 ). Moreover, the PHGPE has such a wide electrochemical window that it is compatible with both high-voltage LiNi 0.8 Co 0.1 Mn 0.1 O 2 and sulfur cathodes to realize long cycle life. Additionally, a LiFePO 4 based pouch size LMB also presents promising performance. Therefore, this study of the PHGPE has opened a new avenue for high-performance LMBs with excellent versatility and practicality.

25 ENERGY STORAGE↗

DOE OSTI2023

Microporous structure in silicone polymers with polylactic acid and polyethylene glycol microspheres: Synthesis and characterization

Porosity in polymer parts could improve mechanical properties while reducing weight and density. Individual pores must be small enough that bulk mechanical properties are retained but significant enough to reduce the density of the material. Filler materials can be incorporated into polymer matrices and removed to create voids. Here, this paper focuses on creating spherical pores in silicone matrices by inserting and removing core materials or fillers during post-cure. A method was developed to incorporate pore formers into an uncured polymer resin and remove the pore former materials from the matrix. Polyethylene glycol (PEG) and polylactic acid (PLA) microspheres were created and incorporated into silicone matrices as pore formers. Microspheres with controllable size distributions were formed via an emulsion formation. Microspheres are removed from the matrix by solvent extraction and calcination. Thermogravimetric Analysis (TGA), Dynamic Mechanical Analysis (DMA), and optical imaging characterize and determine the successful creation of micropores in the polymer matrices.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗

DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗