Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2025

Formation of Composite SiC/SiC Joints by Embedded Wire Chemical Vapor Deposition

The joining of ceramic monoliths or composites to date has primarily been limited to the formation of brittle monolithic joints using heterogeneous (dissimilar) materials, similar to brazing in metals. The development of a damage-tolerant joint layer by SiC fiber reinforcements is demonstrated here. Tube workpieces made of SiC fiber-SiC matrix composite are joined using a nonwoven SiC fiber mat densified by embedded wire chemical vapor deposition (EWCVD), creating a fiber-reinforced weld-like joint by homogeneous joining. EWCVD uses a localized heating method to target deposition and growth to the joint region specifically, while minimizing thermal damage to the surrounding composite tube material. X-ray computed tomography (XCT) is used to nondestructively characterize as-made joints for relative density, adhesion, and composition. In situ XCT analysis during mechanical testing revealed crack deflections in the bonding layer, which indicates a toughening mechanism typical of ceramic matrix composite phase. Gas permeation testing of these proof-of-concept composite joints identified relatively high leak rates in comparison to fully coated SiC/SiC composite tube workpieces. In conclusion, the novelty of the composite joining method and current technology challenges, including gas permeability, are discussed in comparison with traditional ceramic joints and materials.

SiC↗
DOE OSTI2025

Towards exact finite temperature electronic structure in solids and molecules (Final Technical Report)

This report describes the University of Iowa portion of a project that is now continuing at Michigan State University. We are developing novel methods, algorithms, and software to enable simulations of molecules and materials at high temperature. This is a key challenge in chemistry and materials science. By refining an approach called Density Matrix Quantum Monte Carlo (DMQMC), we developed faster and more accurate ways to conduct these simulations. These advances will help us understand how temperature affects the behavior of electrons, chemical bonds, and phase transitions in solids and molecules. These breakthroughs are especially important for applications where light and heat drive chemical reactions, superconductivity, and materials used in energy and sensing. In addition, this project involved the development of the open-source HANDE-QMC software package, supporting the broader community in benchmarking and developing finite-temperature electronic structure methods.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Controlled expansion for pairing in a polarized band with strong repulsion

Can strong repulsive interactions be shown to give rise to pairing in a controlled way? We find that for a single flavor polarized band, there is a small expansion parameter in the low density limit, once the Bloch wavefunction form factor is taken into account. A perturbative expansion is possible, even if the interaction is much stronger than the Fermi energy ε F . As a matter of principle, our work shows analytically how strong pairing can emerge from strong repulsion. Here, we illustrate our method with two examples: a 2D Dirac model and a 1D tight binding model with two orbitals. In the latter case, using density matrix renormalization group, we show that the analytical theory indeed guided us to discover the parameter regime where p-wave pairing with order-1 strength is dominant.

Ferromagnetic superconductors↗
DOE OSTI2020

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗
DOE OSTI2021

Application of HRGS for forensic characterization of uranium oxides, pure uranium metals and uranium alloys

A nondestructive iterative method for uranium-bearing material characterization with HRGS developed earlier in (J. appl. Rad. Isotopes, 166 (2020) 109433) is applied to determine matrix densities and uranium isotope masses of a variety of uranium materials, namely uranium ore, UO 2 and U 3 O 8 powders, fuel elements in the form of UO2 microspheres, uranium metal and uranium alloys. It is shown that U 3 O 8 powders with uranium mass fraction of about 84% can be distinguished from the powders of UO 2 with uranium mass fraction of about 87%; uranium products in the form of liquid or loose powder with matrix density of 0.5-2.0 g/cm 3 can be distinguished from uranium products in the form of compacted fuel elements with matrix density of 6.0-10.0 g/cm 3 and from pure metal uranium and uranium alloys with matrix density of 14.0- 19.0 g/cm 3 . In fuel microspheres based on UO2 the uranium mass fraction 88.02% measured by HRGS is consistent, within the measurement uncertainties, with the results of isotope dilution mass spectrometry 87.76±0.64 % and also is confirmed by X-ray diffraction technique. The uranium mass fraction of the uranium ore estimated as 0.08% by HRGS is consistent, within the measurement uncertainties, with the value 0.09±0.01% determined with WDXRF. Densities of two different uranium metal samples, estimated as 18.42 g/cm 3 and 19.33 g/cm 3 by HRGS are consistent with values 18.24±0.55 g/cm 3 and 18.86±0.59 g/cm 3 , respectively, obtained by the gas pycnometry technique.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Accelerated Coupled Monte Carlo-Thermal Hydraulic Calculations using a Hybrid GTF-Diffusion-based Prediction Block: First Results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work [1, 2] investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections’ variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Reconstructing thermal quantum quench dynamics from pure states

Simulating the nonequilibrium dynamics of thermal states is a fundamental problem across scales from high-energy to condensed-matter physics. Quantum computers may provide a way to solve this problem efficiently. Preparing a thermal state on a quantum computer is challenging, but there exist methods to circumvent this by computing a weighted sum of time-dependent matrix elements in a convenient basis. Further, while the number of basis states can be large, in this paper we show that it can be reduced by simulating only the largest density matrix elements by weight, capturing the density matrix to a specified precision. Leveraging Hamiltonian symmetries enables further reductions. This approach paves the way to more accurate thermal-state dynamics simulations on near-term quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Polarization Dependent Excitation and High Harmonic Generation from Intense Mid-IR Laser Pulses in ZnO

The generation of high order harmonics from femtosecond mid-IR laser pulses in ZnO has shown great potential to reveal new insight into the ultrafast electron dynamics on a few femtosecond timescale. In this work we report on the experimental investigation of photoluminescence and high-order harmonic generation (HHG) in a ZnO single crystal and polycrystalline thin film irradiated with intense femtosecond mid-IR laser pulses. The ellipticity dependence of the HHG process is experimentally studied up to the 17th harmonic order for various driving laser wavelengths in the spectral range 3–4 µm. Interband Zener tunneling is found to exhibit a significant excitation efficiency drop for circularly polarized strong-field pump pulses. For higher harmonics with energies larger than the bandgap, the measured ellipticity dependence can be quantitatively described by numerical simulations based on the density matrix equations. The ellipticity dependence of the below and above ZnO band gap harmonics as a function of the laser wavelength provides an efficient method for distinguishing the dominant HHG mechanism for different harmonic orders.

tunneling excitation↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

Synthesis of core/shell nanocrystals with ordered intermetallic single-atom alloy layers for nitrate electroreduction to ammonia

Structurally ordered intermetallic nanocrystals (NCs) and single-atom catalysts (SACs) are two emerging catalytic motifs for sustainable chemical production and energy conversion. However, both have synthetic limitations which can lead to the aggregation of NCs or metal atoms. Single-atom alloys (SAAs), which contain isolated metal atoms in a host metal, can overcome the aggregation concern because of the thermodynamic stabilization of single atoms on host metal surfaces. We report a direct solution-phase synthesis of Cu/CuAu core/shell NCs with tunable SAA layers. This synthesis can be extended to other Cu/CuM (M = Pt, Pd) systems, in which M atoms are isolated in the copper host. Using this method, the density of SAAs on a copper surface can be controlled, resulting in both low and high densities of single atoms. Alloying gold into the copper matrix introduced ligand effects that optimized the chemisorption of *NO 3 and *N. As a result, the densely packed Cu/CuAu material demonstrated a high selectivity toward NH 3 from the electrocatalytic nitrate reduction reaction with an 85.5% Faradaic efficiency while maintaining a high yield rate of 8.47 mol h -1 g -1 . This work advances the design of atomically precise catalytic sites by creating core/shell NCs with SAA atomic layers, opening an avenue for broad catalytic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Determining the oxidation behavior of matrix graphite

This work presents the oxidation behavior of matrix graphite in air. Matrix graphite, graphite powder/flakes bonded by a small amount of non-graphitic carbon, surrounds coated fuel particles in order to form cylindrical fuel compacts (in prismatic core designs) or spheres (in pebble-bed reactor designs). This work focuses on oxidation tests conducted on two matrix graphite materials, one provided by Kairos Power and the other A3 matrix graphite. Some of the tests followed American Society for Testing and Materials (ASTM) oxidation testing standards using a vertical furnace system and others were performed in a thermogravimetric analyzer (TGA). It was determined that, at temperatures of 450 °C–700 °C, the oxidation rate of the Kairos matrix graphite follows the Arrhenius equation. In comparison with A3 matrix graphite, the Kairos matrix graphite shows better oxidation resistance at high temperatures (≥550 °C), but also a higher oxidation rate at low temperatures. Both the A3 matrix graphite and the Kairos matrix graphite materials may experience preferential oxidation of the partially graphitized binder. An oxygen penetration gradient was also observed when using the three characterization methods (i.e., optical microscope, x-ray tomography [XCT], and density profile by the lathe) enlisted in this research. In conclusion, the oxygen penetration depth increases with decreasing isothermal oxidation temperature, while the center of the oxidized samples (10% weight loss) remains almost untouched even at 500 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

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

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Entanglement and Bell inequalities with boosted t t ¯

The Large Hadron Collider provides a unique opportunity to study quantum entanglement and violation of Bell inequalities at the highest energy available today. In this paper, we will investigate these quantum correlations with top quark pair production, which represents a system of two-qubits. The semileptonic top pair channel has a factor of 6 increase in statistics and easier reconstruction with respect to the dileptonic channel. Although measuring the spin polarization of the hadronic top quark is known to be challenging, our study indicates that it is feasible to reconstruct the spin density matrix of the two-qubit system using an optimal hadronic polarimeter. This is achieved with the aid of jet substructure techniques and NN-inspired reconstruction methods, which improve the mapping between subjets and quarks. We find that entanglement can already be observed at more than the 5 σ level with existing data, and violation of Bell inequalities may be probed above the 4 σ level at the HL-LHC with 3 ab − 1 of data. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Well-designed Crosslinked Polymer Electrolyte Enables High Ionic Conductivity and Enhanced Salt Solvation

A new facile single-step method to fabricate crosslinked polymer electrolyte membranes consisting of branched poly(ethyleneimine), (PEI) and poly(ethylene oxide), (PEO) is demonstrated. The membranes exhibit excellent ionic conductivity (1.2 × 10 –3 S cm –1 at 80 °C) with minimal addition of plasticizer (20 wt%). The amine functional group in the PEI-PEO crosslinked matrix provides Lewis basic and hydrogen bonding characteristics that facilitate the dissolution of lithium salt and enables a higher cation transport number than a PEO crosslinked matrix. The glass transition temperature, degree of crystallinity, and room temperature storage modulus increases with decreasing crosslink density and increasing ratio of free amines. The resultant ionic conductivity and mechanical strength can be flexibly tailored by varying the molar ratio of free amine moieties in the crosslinked PEI-PEO matrix. This study provides an improved synthesis method, in-depth characterization, and fundamental insights on the effect of free amine moieties on the transport properties of a highly conductive gel polymer electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗