Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A photoemission moments model using density functional and transfer matrix methods applied to coating layers on surfaces: Theory

Recent experimental measurements of a bulk material covered with a small number of graphene layers reported by Yamaguchi et al. [NPJ 2D Mater. Appl. 1, 12 (2017)] (on bialkali) and Liu et al.[Appl. Phys. Lett. 110, 041607 (2017)] (on copper) and the needs of emission models in beam optics codes have lead to substantial changes in a Moments model of photoemission. The changes account for (i) a barrier profile and density of states factor based on density functional theory (DFT) evaluations, (ii) a Drude-Lorentz model of the optical constants and laser penetration depth, and (iii) a transmission probability evaluated by an Airy Transfer Matrix Approach. Importantly, the DFT results lead to a surface barrier profile of a shape similar to both resonant barriers and reflectionless wells: the associated quantum mechanical transmission probabilities are shown to be comparable to those recently required to enable the Moments (and Three Step) model to match experimental data but for reasons very different than the assumption by conventional wisdom that a barrier is responsible. The substantial modifications of the Moments model components, motivated by computational materials methods, are developed. The results prepare the Moments model for use in treating heterostructures and discrete energy level systems (e.g., quantum dots) proposed for decoupling the opposing metrics of performance that undermine the performance of advanced light sources like the x-ray Free Electron Laser. The consequences of the modified components on quan-tum yield, emittance, and emission models needed by beam optics codes are discussed. Published by AIP Publishing. https://doi.org/10.1063/1.5008600

36 MATERIALS SCIENCE↗
DOE OSTI2015

Communication: Towards ab initio self-energy embedding theory in quantum chemistry

The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green’s function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.

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

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Increasing the power density when using inert matrix fuels to reduce production of transuranics

Reducing the production of transuranics is a goal of most advanced nuclear fuel cycles. One way to do this is to recycle the transuranics into the same reactors that are currently producing them using an inert matrix fuel. In previous work we have modeled such a reactor where 72%, of the core is comprised of standard enriched uranium fuel pins, with the remaining 28% fuel made from Yttria stabilized zirconium, in which transuranics are loaded. A key feature of this core is that all of the transuranics produced by the uranium fuel assemblies are later burned in inert matrix fuel assemblies. It has been shown that this system can achieve reductions in transuranic waste of more than 86%. The disadvantage of such a system is that the core power rating must be significantly lower than a standard pressurized water reactor. One reason for the lower power is that high burnup of the uranium fuel precludes a critical level of reactivity at the end of the campaign. Increasing the uranium enrichment and changing the pin pitch are two ways to increase burnup while maintaining criticality. In this paper we use MCNPX and a linear reactivity model to quantify the effect of these two parameters on the end of campaign reactivity. Importantly, we show that in the region of our proposed reactor, enrichment increases core reactivity by 0.02 per percent uranium 235 and pin pitch increases reactivity by 0.02 per mm. Reactivity is lost at a rate of 0.005 per MWd/kgIHM uranium burnup. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

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↗
DOE OSTI2014

Average local ionization energy generalized to correlated wavefunctions

The average local ionization energy function introduced by Politzer and co-workers [Can. J. Chem. 68, 1440 (1990)] as a descriptor of chemical reactivity has a limited utility because it is defined only for one-determinantal self-consistent-field methods such as the Hartree–Fock theory and the Kohn–Sham density-functional scheme. We reinterpret the negative of the average local ionization energy as the average total energy of an electron at a given point and, by rewriting this quantity in terms of reduced density matrices, arrive at its natural generalization to correlated wavefunctions. The generalized average local electron energy turns out to be the diagonal part of the coordinate representation of the generalized Fock operator divided by the electron density; it reduces to the original definition in terms of canonical orbitals and their eigenvalues for one-determinantal wavefunctions. The discussion is illustrated with calculations on selected atoms and molecules at various levels of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

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↗
DOE OSTI2014

Transition properties from the Hermitian formulation of the coupled cluster polarization propagator

Theory of one-electron transition density matrices has been formulated within the time-independent coupled cluster method for the polarization propagator [R. Moszynski, P. S. Żuchowski, and B. Jeziorski, Coll. Czech. Chem. Commun. 70, 1109 (2005)]. Working expressions have been obtained and implemented with the coupled cluster method limited to single, double, and linear triple excitations (CC3). Selected dipole and quadrupole transition probabilities of the alkali earth atoms, computed with the new transition density matrices are compared to the experimental data. Good agreement between theory and experiment is found. The results obtained with the new approach are of the same quality as the results obtained with the linear response coupled cluster theory. The one-electron density matrices for the ground state in the CC3 approximation have also been implemented. The dipole moments for a few representative diatomic molecules have been computed with several variants of the new approach, and the results are discussed to choose the approximation with the best balance between the accuracy and computational efficiency.

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

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↗
DOE OSTI2017

Cosmic Web of Galaxies in the COSMOS Field: Public Catalog and Different Quenching for Centrals and Satellites

We use a mass complete (log(M/M{sub ⊙}) ⩾ 9.6) sample of galaxies with accurate photometric redshifts in the COSMOS field to construct the density field and the cosmic web to z = 1.2. The comic web extraction relies on the density field Hessian matrix and breaks the density field into clusters, filaments, and the field. We provide the density field and cosmic web measures to the community. We show that at z ≲ 0.8, the median star formation rate (SFR) in the cosmic web gradually declines from the field to clusters and this decline is especially sharp for satellites (∼1 dex versus ∼0.5 dex for centrals). However, at z ≳ 0.8, the trend flattens out for the overall galaxy population and satellites. For star-forming (SF) galaxies only, the median SFR is constant at z ≳ 0.5 but declines by ∼0.3–0.4 dex from the field to clusters for satellites and centrals at z ≲ 0.5. We argue that for satellites, the main role of the cosmic web environment is to control their SF fraction, whereas for centrals, it is mainly to control their overall SFR at z ≲ 0.5 and to set their fraction at z ≳ 0.5. We suggest that most satellites experience a rapid quenching mechanism as they fall from the field into clusters through filaments, whereas centrals mostly undergo a slow environmental quenching at z ≲ 0.5 and a fast mechanism at higher redshifts. Our preliminary results highlight the importance of the large-scale cosmic web on galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

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↗
NASA NTRS1991

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.↗
DOE OSTI2021

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↗
DOE OSTI2022

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↗
NASA NTRS1971

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.↗
DOE OSTI2020

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↗
DOE OSTI2025

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 (↗