Fermionic Reduced Density Low-Rank Matrix Completion, Noise Filtering, and Measurement Reduction in Quantum Simulations
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “density matrix”
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.
A multisite decomposition of the tensor network path integrals
Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Ion Exchange Technology Development in Support of the Urine Processor Assembly
The urine processor assembly (UPA) on the International Space Station (ISS) recovers water from urine via a vacuum distillation process. The distillation occurs in a rotating distillation assembly (DA) where the urine is heated and subjected to sub-ambient pressure. As water is removed, the original organics, salts, and minerals in the urine become more concentrated and result in urine brine. Eventually, water removal will concentrate the urine brine to super saturation of individual constituents, and precipitation occurs. Under typical UPA DA operating conditions, calcium sulfate or gypsum is the first chemical to precipitate in substantial quantity. During preflight testing with ground urine, the UPA achieved 85% water recovery without precipitation. However, on ISS, it is possible that crewmember urine can be significantly more concentrated relative to urine from ground donors. As a result, gypsum precipitated in the DA when operating at water recovery rates at or near 85%, causing the failure and subsequent re14 NASA Tech Briefs, September 2013 placement of the DA. Later investigations have demonstrated that an excess of calcium and sulfate will cause precipitation at water recovery rates greater than 70%. The source of the excess calcium is likely physiological in nature, via crewmembers' bone loss, while the excess sulfate is primarily due to the sulfuric acid component of the urine pretreatment. To prevent gypsum precipitation in the UPA, the Precipitation Prevention Project (PPP) team has focused on removing the calcium ion from pretreated urine, using ion exchange resins as calcium removal agents. The selectivity and effectiveness of ion exchange resins are determined by such factors as the mobility of the liquid phase through the polymer matrix, the density of functional groups, type of functional groups bound to the matrix, and the chemical characteristics of the liquid phase (pH, oxidation potential, and ionic strength). Previous experience with ion exchange resins has demonstrated that the most effective implementation for an ion exchange resin is a cartridge, or column, in which the resin is contained. Based on the results of equilibrium and sub-scale dynamic column testing, a possible solution for mitigating the calcium precipitation issue on the ISS has been identified. From an original pool of 13 ion exchange resins, two candidates have been identified that demonstrate substantial calcium removal on the sub-scale. The dramatic reduction in resin performance from published calcium uptake demonstrates the need for thorough evaluation of resins at the low pH and strong oxidizing environment present in the UPA. Chemical variations in the influent (calcium concentrations and pretreatment dosing) appear to have a noticeable impact on the calcium capacity of the resin. Low calcium concentrations and high pretreatment dosing will likely result in a decrease in calcium capacity. Conversely, low pre trea t - ment dosing will likely result in an increase in calcium capacity. In contrast, investigations at a variety of flow rates, length-to-diameter ratios, resin volumes, and flow regimes (continuous versus pulsed) show that changes in physical parameters do not have substantial impacts on resin performance in the very low specific velocity ranges of interest. This result is particularly useful because most commercial applications at higher specific velocities do show a relatively strong relationship between flow and capacity. The lack of a strong relationship will allow more flexibility in the implementation of an ion exchange bed for flight. Verification of subscale tests with flight-scale resin beds is recommended prior to implementation in the on-orbit UPA.
Mitchell, Julie↗
On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH
We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Accelerating the density-functional tight-binding method using graphical processing units
Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Low-density, high-strength intermetallic matrix composites by XD (trademark) synthesis
A feasibility study was conducted to evaluate the potential of particulate composites based on low-density, L1(sub 2) trialuminide matrices for high-temperature applications. The compounds evaluated included Al22Fe3Ti8 (as a multiphase matrix), Al67Ti25Cr8, and Al66Ti25Mn9. The reinforcement consisted of TiB2 particulates. The TiB2 composites were processed by ingot and powder metallurgy techniques. Microstructural characterization and mechanical testing were performed in the hot-pressed and hot-isostatic-pressed condition. The casting were sectioned and isothermally forged into pancakes. All the materials were tested in compression as a function of temperature, and at high temperatures as a function of strain rate. The test results are discussed.
Kumar, K. S.↗
Multiscale characterization and representation of variability in ceramic matrix composites
Low density, high strength, and high creep and oxidation resistance properties of ceramic matrix composites (CMCs) make them an ideal choice for use in extreme environments in space and military applications. This paper presents a detailed characterization study of structural and manufacturing flaws in Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs at different length-scales. Energy-dispersive spectroscopy (EDS) is used for the chemical characterization of the material’s elemental constituents. High-resolution multiscale graphs obtained from scanning electron microscope (SEM) and confocal laser scanning microscope (LSM) are used to characterize the distribution and morphology of defects at different length scales. This is followed by the classification and quantification of the common manufacturing defects. An image processing algorithm based on the image segmentation process is developed to quantify the variability of various scale-dependent architectural parameters. Finally, a three-dimensional stochastic representative volume element (SRVE) generation algorithm is developed to provide precise representations of material textures at multiple length scales. The developed algorithm accurately accounts for material features and flaws based on a range of multiscale structural and defects characterization results.
36 MATERIALS SCIENCE↗
237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes
A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Granular two-phase insulation systems
Easily prepared system, consisting of matrix of hollow zirconia microspheres containing dispersed tungsten powder, produces minimum-cost, prototype test specimen. Combination represents basic concept of highly reflective dispersed phase in low density insulative matrix and is stable at 2200 K. Other combinations of materials are suggested.
Carpenter, H. W.↗
Oxide dispersoid coherency of a ferritic-martensitic 12Cr oxide-dispersion-strengthened alloy under self-ion irradiation
Ferritic-martensitic oxide-dispersion-strengthened (ODS) alloy has shown excellent mechanical property and high radiation tolerance. However, the stability of dispersoids during displacive irradiation in the individual ferritic and martensitic phases is unclear. In this work, the correlation among dispersoid coherency, size, density, and matrix phase are studied in dual-phase 12Cr ODS after 100 peak displacements-per-atom (dpa) irradiation at 475°C, using 3.5 MeV Fe 2+ self-ions. The size and density changes of coherent and incoherent dispersoids were analyzed as a function of irradiation depth. The average dispersoid size decreased after irradiation in both phases, and the large incoherent dispersoids in the tempered martensite phase underwent a more dramatic change than those in the ferrite phase. The dispersoid density significantly increased in the ferrite phase within the irradiated region, mostly resulting in an increase in coherent dispersoid density. On the other hand, only a small change in density was observed in the tempered martensite phase for both coherent and incoherent dispersoids. This study shows that dispersoid evolution is dramatically different in ferrite and tempered martensite phases.
36 MATERIALS SCIENCE↗
Elucidating the impact of fiber source on polypropylene/hemp composite performance for the automotive industry
Given their high strength-to-weight ratio, there is an ever-increasing volume of plastics being used in the automotive industry as plastics aid in the charge to lightweight vehicles for improved fuel efficiency. However, these plastics are often landfilled at their end-of-life, which has given rise to the demand for sustainable materials and waste management alternatives compared to purely synthetic systems. Natural fiber composites have been explored as a viable material option to reduce the environmental impact of plastic use in automobiles while simultaneously ensuring the part performance is not sacrificed. Herein, we explored the use of hemp/polypropylene (PP) composites in which US-sourced hemp is compared to internationally sourced and industrially available hemp. There appears to be a minimal impact on the composite properties regardless of fiber sourcing, and the addition of a natural filler to the PP matrix results in up to a 367% increase in Young's modulus, 126% increase in heat deflection temperature, and comparable water uptake performance. It should be noted that the natural filler addition does increase density by up to 13% due to the higher density of natural fibers compared to the low-density PP matrix. A modified rule of mixtures calculation revealed that the composite materials produced in this study demonstrated good agreement with analytical modeling. Finally, a screening analysis was performed exploring the transportation of hemp fibers, and the results build a strong case for regionalized manufacturing of automotive parts.
Hubbard, Amber M. [Oak Ridge National Laboratory (↗
High-performance strategies for the recent MRSF-TDDFT in GAMESS
Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy
We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems
Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.
Calculus of variations↗
Ethynylated aromatics as high temperature matrix resins
Difunctional and trifunctional arylacetylenes were used as monomers to form thermoset matrix resin composites. Composites can be hot pressed at 180 C to react 80 percent of the acetylene groups. Crosslinking is completed by postcuring at 350 C. The postcured resins are thermally stable to nominally 460 C in air. As a result of their high crosslink density, the matrix exhibits brittle failure when unaxial composites are tested in tension. Failure of both uniaxial tensile and flexural specimens occurs in shear at the fiber matrix interface. Tensile fracture stresses for 0 deg composites fabricated with 60 v/o Celion 6K graphite fiber were 827 MPa. The strain to failure was 0.5 percent. Composites fabricated with 8 harness satin Celion cloth (Fiberite 1133) and tested in tension also failed in shear at tensile stresses of 413 MPa.
Hurwitz, Frances I.↗
Ethynylated aromatics as high temperature matrix resins
Difunctional and trifunctional arylacetylenes were used as monomers to form thermoset matrix resin composites. Composites can be hot-pressed at 180 C to react 80 percent of the acetylene groups. Crosslinking is completed by postcuring at 350 C. The postcured resins are thermally stable to nominally 460 C in air. As a result of their high crosslink density, the matrix exhibits brittle failure when uniaxial composites are tested in tension. Failure of both uniaixial tensile and flexural specimens occurs in shear at the fiber-matrix interface. Tensile fracture stresses for 0-deg composites fabricated with 60 v/o Celion 6K graphite fiber were 827 MPa. The strain to failure was 0.5 percent. Composites fabricated with 8 harness satin Celion cloth (Fiberite 1133) and tested in tension also failed in shear at tensile stresses of 413 MPa.
Hurwitz, F. I.↗
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.↗
Design and Analysis of Battery Thermal Management Systems
Thermal management of battery cell packs is a critically needed technology. The purpose of this work is to design new and improved Battery Thermal Management Systems (BTMS) for use in electric airplanes. The BTMS should be 3D printable and hold twelve to sixteen 18650 batteries. To minimize aircraft’s weight, an actively air-cooled battery pack was chosen to avoid the excess weight of water-cooled and phase change material (PCM) cooled packs. As a result, significant changes were made to remove all metal in the pack and replace it with lower density polymer matrix composites (PMCs) which can be additively manufactured. Two air cooled battery pack configurations (traditional propeller fans & bladeless fans) were designed, modeled, and compared. These packs were first modeled in SolidWorks 2021 3D CAD, then imported into COMSOL MultiPhysics to be studied using the “Heat Transfer in Solids and Fluids” module. The design with the bladeless fans eliminated the need to use high conductivity heavy metal to remove unwanted thermal energy. These bladeless fans were designed entirely out of PMCs. This thermal pack design weighs 0.04 kg less than that of the traditional propeller fans design and has increased the battery pack energy density by 8.25 Wh/kg.
Thermal Management System for Battery Packs↗