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 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Stable Solid Electrolyte Interphase Layer Formed by Electrochemical Pretreatment of Gel Polymer Coating on Li Metal Anode for Lithium–Oxygen Batteries

Lithium (Li)-oxygen (O2) batteries (LOBs) exhibit the highest theoretical specific energy density among candidates of the next-generation energy storage systems, but the instability of Li metal anode (LMA), air electrode, and electrolyte largely limit the practical energy density of these batteries. Herein, we report an effective method to protect LMA against side reactions between LMA and crossover contaminants such as highly reactive oxygen moieties. A solid electrolyte interphase (SEI) layer rich in inorganic components was formed on a LMA coated with polyethylene oxide thin-film through an in-situ electrochemical pre-charging step under oxygen atmosphere. This uniformly distributed stiff SEI layer interacted with flexible polymer matrix and formed a submicron-sized gel-like polymer layer. This polymer supported SEI layer leads to much longer cycle life (130 vs. 65 cycles) as compared to those of pristine cell under the same testing conditions (1.0 mAh cm-2 at 0.2 mA cm-2 in the voltage range of 2 to 5 V without any catalysts). It is also very effective during low voltage (2 to 4.5 V) cycling with a redox mediator (0.1 and 0.15M (2,2,6,6-tetramethylpiperidin-1-yl) oxidanyl). Therefore, this approach can be used to stabilize LMA/electrolytes interphase and improve the cycle life of rechargeable LOBs.

lithium-oxygen batteries, lithium metal anodes, ar↗

DOE OSTI2022

Microstructural and physicochemical origins of electroless copper deposition on graphite enhanced by acid pretreatment

We report Acid treatment is the most widely used surface modification method for enhancing the electroless metal deposition (EMD) on carbon reinforcement materials (CRMs) for metal matrix nanocomposites. However, specific microstructural and physicochemical origins of the enhanced EMD on carbon surfaces by acid treatments have been rarely studied. Here, we investigated the effects of the nitric acid treatment on graphite, a prototypical combination of acid treatment and CRM, on the fidelity of Cu EMD and their structural and chemical origins. Complementary materials characterizations and density functional theory calculations revealed the acid-induced formation of broken C–C/C=C graphitic bonds and resulting surface micropores on graphite; this enabled a uniform dispersion of catalytic Sn/Pd nanoparticles during pre-EMD sensitization/activation processes via spontaneous binding of Sn and Pd ions and, consequently, a much more uniform Cu layer EMD compared to the untreated graphite. We proposed a general mechanism illustrating how the acid-induced microstructural and chemical modifications of carbon surface affected the spatial uniformity of catalytic metal reduction during EMD and, finally, the quality of deposited metal layer. The results clearly reveal the origins of the enhanced EMD on carbon materials by acid treatments, providing guidelines for optimizing EMD on general CRMs for high-performance metal matrix nanocomposites.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Global quark spin correlations in relativistic heavy ion collisions

The observation of the vector meson’s global spin alignment by the STAR Collaboration reveals that strong spin correlations may exist for quarks and antiquarks in relativistic heavy-ion collisions in the normal direction of the reaction plane. We propose a systematic method to describe such correlations in the quark matter. We classify them as local and long range quark spin correlations in the system. We show in particular that the effective quark spin correlations contain the genuine spin correlations originated directly from the dynamical process as well as those induced by averaging over other degrees of freedom. We also show that such correlations can be studied by measuring the vector meson’s spin density matrix and hyperon-hyperon and hyperon-anti-hyperon spin correlations. We present the relationships between these measurable quantities and spin correlations of quarks and antiquarks. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1989

Kaon-nucleus scattering

Two kinds of number density distributions of the nucleus, harmonic well and Woods-Saxon models, are used with the t-matrix that is taken from the scattering experiments to find a simple optical potential. The parameterized two body inputs, which are kaon-nucleon total cross sections, elastic slope parameters, and the ratio of the real to imaginary part of the forward elastic scattering amplitude, are shown. The eikonal approximation was chosen as the solution method to estimate the total and absorptive cross sections for the kaon-nucleus scattering.

Hong, Byungsik↗

DOE OSTI2021

Solving puzzles in deformed JT gravity: phase transitions and non-perturbative effects

Recent work has shown that certain deformations of the scalar potential in Jackiw-Teitelboim gravity can be written as double-scaled matrix models. However, some of the deformations exhibit an apparent breakdown of unitarity in the form of a negative spectral density at disc order. We show here that the source of the problem is the presence of a multi-valued solution of the leading order matrix model string equation. While for a class of deformations we fix the problem by identifying a first order phase transition, for others we show that the theory is both perturbatively and non-perturbatively inconsistent. Aspects of the phase structure of the deformations are mapped out, using methods known to supply a non-perturbative definition of undeformed JT gravity. Some features are in qualitative agreement with a semi-classical analysis of the phase structure of two-dimensional black holes in these deformed theories.

2D gravity↗

DOE OSTI2021

DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various density functional theory packages

We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2007

Sequestration of Single-Walled Carbon Nanotubes in a Polymer

Sequestration of single-walled carbon nanotubes (SWCNs) in a suitably chosen polymer is under investigation as a means of promoting the dissolution of the nanotubes into epoxies. The purpose of this investigation is to make it possible to utilize SWCNs as the reinforcing fibers in strong, lightweight epoxy-matrix/carbon-fiber composite materials. SWCNs are especially attractive for use as reinforcing fibers because of their stiffness and strength-to-weight ratio: Their Young s modulus has been calculated to be 1.2 TPa, their strength has been calculated to be as much as 100 times that of steel, and their mass density is only one-sixth that of steel. Bare SWCNs cannot be incorporated directly into composite materials of the types envisioned because they are not soluble in epoxies. Heretofore, SWCNS have been rendered soluble by chemically attaching various molecular chains to them, but such chemical attachments compromise their structural integrity. In the method now under investigation, carbon nanotubes are sequestered in molecules of poly(m-phenylenevinylene-co-2,5-dioctyloxy-p-phenylenevinylene) [PmPV]. The strength of the carbon nanotubes is preserved because they are not chemically bonded to the PmPV. This method exploits the tendency of PmPV molecules to wrap themselves around carbon nanotubes: the wrapping occurs partly because there exists a favorable interface between the conjugated face of a nanotube and the conjugated backbone of the polymer and partly because of the helical molecular structure of PmPV. The constituents attached to the polymer backbones (the side chains) render the PmPV-wrapped carbon nanotubes PmPV soluble in organic materials that, in turn, could be used to suspend the carbon nanotubes in epoxy precursors. At present, this method is being optimized: The side chains on the currently available form of PmPV are very nonpolar and unable to react with the epoxy resins and/or hardeners; as a consequence, SWCN/PmPV composites have been observed to precipitate out of epoxies while the epoxies were being cured. If the side chains of the PmPV molecules were functionalized to make them capable of reacting with the epoxy matrices, it might be possible to make progress toward practical applications. By bonding the side chains of the PmPV molecules to an epoxy matrix, one would form an PmPV conduit between the epoxy matrix and the carbon nanotubes sequestered in the PmPV. This conduit would transfer stresses from the epoxy matrix to the nanotubes. This proposed load-transfer mode is similar to that of the current practice in which silane groups are chemically bonded to both the epoxy matrices and the fibers (often glass fibers) in epoxymatrix/fiber composites.

Bley, Richard A.↗

DOE OSTI2022

Ab initio study on fcc Pr with correlation matrix renormalization theory

Here, we studied fcc praseodymium (Pr) with the ab initio correlation matrix renormalization theory (CMRT) explicitly calculating the Coulomb interactions among basis orbitals and without using any adjustable parameters to work for strongly correlated electron systems. We calculated its total energy in a paramagnetic ground state and studied the role of the correlated 4f electrons in the system. Good agreement was obtained between CMRT and experiments in the pressure volume dependence of the fcc phase. We also compared the CMRT results against other theoretical methods including local density approximation + dynamical mean-field theory and showed consistent results among them. Moreover, we found the normalized local charge fluctuation of the 4f electrons as the function of volume exhibits a clear slope change at the volume collapse region, indicating a switch in their correlation nature there.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Simple and Robust Carrier-Based PWM Technique for Single-Stage Three-phase Rectifier Indirect Matrix Converter

Indirect matrix converters (IMC) can offer high efficiency, high power density, and bidirectional power flow capability with a single-stage three-to single-phase AC-AC power conversion integrating a conventional three-phase active rectifier with a primary side full-bridge inverter of a dual active bridge converter. However, the IMCs generally require space vector modulation techniques, leading to complex sectors and dwell times calculations. A practical implementation of the modulation methods with a commonly used digital signal processor becomes even more challenging when considering power factor correction, harmonics reduction, and deadtime between sector transitions. This paper proposes a simple and robust carrier-based pulse width modulation (PWM) technique determined through a voltage-source inverter (VSI) to current-source inverter (CSI) to IMC conversion. This translation enables the adoption of all conventional VSI control methods with active and reactive power control in the IMCs through simple carrier and reference waveform comparison. The VSI vectors are first one-on-one matched with CSI vectors and translated to IMC modulation signals without calculating sectors and dwell times. Nine different carrier waveforms are proposed, and the converter performance is analyzed. Grid current THD and efficiency are estimated for each modulation method. The symmetric A and B methods had the lowest THD of the grid current, 2.80%, and the efficiency remained at approximately 98.3%, making it the most suitable carrier-based method proposed in this paper. The space vector PWM (SVPWM) modulation method was implemented to the DSP and generated the expected twelve gate signals.

Benson, Mikayla↗

DOE OSTI2024

Interfacial Engineering Using Covalent Organic Frameworks in Polymer Composites for High‐Temperature Electrostatic Energy Storage

Abstract The use of inorganic nanofillers has been an effective method to improve high‐temperature capacitive performance of dielectric polymers, though there are unmet challenges such as undesirable organic–inorganic compatibility, and low efficiencies and energy densities. Herein, a surface functionalization strategy using covalent organic frameworks (COFs) is employed to address such challenges in realizing high‐performing polymer composites. Specifically, core–shell structured nanoparticles, where ZrO 2 nanoparticles act as the core and a COF material forms the shell, are constructed and composited with the polyetherimide (PEI) matrix. The design leverages the high electron affinity ( E A ) of the outer COF shell to create energy traps, thereby capturing free charges and limiting electrical conduction. Concurrently, the low E A and wide bandgap of the ZrO 2 core introduce energy barriers to impede charge injection and migration. This orchestrated “energy level cascade” results in a marked reduction of leakage current and energy loss. The resulting polymer composite showcases an impressive discharged energy density of 6.21 J cm −3 at an efficiency above 90%, with a maximum discharged energy density reaching 7.43 J cm −3 at 150 °C. These performance metrics position the PEI/ZrO 2 @COF polymer composite to surpass or be on par with state‐of‐the‐art high‐temperature PEI composites and other advanced polymer dielectrics.

Xie, Zongliang↗

DOE OSTI2021

Effect of Ligand Polarity on the Internal Dipoles and Ferroelectric Distortion in BaTiO 3 Nanocubes

Here, surface adsorbates and surrounding matrix species have been demonstrated to affect the properties of nanoscale ferroelectrics and nanoscale ferroelectric composites; potentially counteracting performance losses that can occur in small particle sizes. In this work, the effects of nonpolar oleic acid (OA) and polar tetrafluoroborate (BF 4 - ) ligand capping on the surface of various sizes of BaTiO 3 nanocubes have been investigated with combined neutron diffraction and neutron pair distribution function (PDF), density functional theory (DFT), and ab initio molecular dynamics (AIMD) methods. The low real space PDF region provides an unobstructed view of rhombohedral (split short and long) Ti-O distances in BaTiO 3 nanocubes, mimicking the well-established order-disorder local structure found in bulk BaTiO 3 . Interestingly, the intermediate-range order in nanocubes is found to be orthorhombic, rather than tetragonal. It is concluded that polar ligands adsorbed at BaTiO 3 surfaces stabilize the correlation length scale of local rhombohedral distortions in ferroelectric nanoparticles relative to nonpolar ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Data assimilation in operator algebras

We develop an algebraic framework for sequential data assimilation of partially observed dynamical systems. In this framework, Bayesian data assimilation is embedded in a nonabelian operator algebra, which provides a representation of observables by multiplication operators and probability densities by density operators (quantum states). In the algebraic approach, the forecast step of data assimilation is represented by a quantum operation induced by the Koopman operator of the dynamical system. Moreover, the analysis step is described by a quantum effect, which generalizes the Bayesian observational update rule. Projecting this formulation to finite-dimensional matrix algebras leads to computational schemes that are i) automatically positivity-preserving and ii) amenable to consistent data-driven approximation using kernel methods for machine learning. Moreover, these methods are natural candidates for implementation on quantum computers. Applications to the Lorenz 96 multiscale system and the El Niño Southern Oscillation in a climate model show promising results in terms of forecast skill and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Coupled Gel/Seeded Growth of Cs2AgBiBr6 Single Crystals for Radiation Detection

Abstract— The all-inorganic lead free Cs2AgBiBr6 perovskite represents a highly promising environment-friendly radiation detector material. The first attempt to grow Cs2AgBiBr6 single crystals using the unique gel growth method is reported in this study. This technique enables controlled reagent diffusion in a gel medium at room temperature, fostering the growth of visible single crystals. It offers advantages over conventional methods, leading to fewer equilibrium and non-equilibrium defects and simplified execution. Interestingly, the as-grown crystals exhibited variations in size and colors within the gel matrix, indicating the influence of Bi solutes in Bi-poor and Bi-rich Cs2AgBiBr6 crystals. Once the crystals were confirmed by powder xray diffraction, a series of comprehensive material characterizations was performed. The resistivity of Cs2AgBiBr6 single crystals was measured at 1.12x107 Ω•cm. Trap density analysis yielded a value of 2.82 x1010 cm-3. The charge carrier mobility was determined to be 32.9 cm2V-1S-1, suggesting the suitability for radiation detector applications. Additional growth optimization is going on and further characterization will be reported during the presentation. Our work helps accelerate the development of lead-free perovskite materials toward radiation detector applications.

James, Ralph B.↗

NASA NTRS1993

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗

DOE OSTI2022

Additive manufacturing of complex objects using refractory matrix materials

A method for the manufacture of a three-dimensional object using a refractory matrix material is provided. The method includes the additive manufacture of a green body from a powder-based refractory matrix material followed by densification via chemical vapor infiltration (CVI). The refractory matrix material can be a refractory ceramic (e.g., silicon carbide, zirconium carbide, or graphite) or a refractory metal (e.g., molybdenum or tungsten). In one embodiment, the matrix material is deposited according to a binder-jet printing process to produce a green body having a complex geometry. The CVI process increases its density, provides a hermetic seal, and yields an object with mechanical integrity. The residual binder content dissociates and is removed from the green body prior to the start of the CVI process as temperatures increase in the CVI reactor. The CVI process selective deposits a fully dense coating on all internal and external surfaces of the finished object.

Terrani, Kurt A.↗

DOE OSTI2024

Additive manufacturing of complex objects using refractory matrix materials

Terrani, Kurt A.↗