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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1988

Quantum mechanical algebraic variational methods for inelastic and reactive molecular collisions

The quantum mechanical problem of reactive or nonreactive scattering of atoms and molecules is formulated in terms of square-integrable basis sets with variational expressions for the reactance matrix. Several formulations involving expansions of the wave function (the Schwinger variational principle) or amplitude density (a generalization of the Newton variational principle), single-channel or multichannel distortion potentials, and primitive or contracted basis functions are presented and tested. The test results, for inelastic and reactive atom-diatom collisions, suggest that the methods may be useful for a variety of collision calculations and may allow the accurate quantal treatment of systems for which other available methods would be prohibitively expensive.

Schwenke, David W.↗
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↗
NASA NTRS2000

Optical Control of Intersubband Absorption in a Multiple Quantum Well-Embedded Semiconductor Microcravity

Optical intersubband response of a multiple quantum well (MQW)-embedded microcavity driven by a coherent pump field is studied theoretically. The n-type doped MQW structure with three subbands in the conduction band is sandwiched between a semi-infinite medium and a distributed Bragg reflector (DBR). A strong pump field couples the two upper subbands and a weak field probes the two lower subbands. To describe the optical response of the MQW-embedded microcavity, we adopt a semi-classical nonlocal response theory. Taking into account the pump-probe interaction, we derive the probe-induced current density associated with intersubband transitions from the single-particle density-matrix formalism. By incorporating the current density into the Maxwell equation, we solve the probe local field exactly by means of Green's function technique and the transfer-matrix method. We obtain an exact expression for the probe absorption coefficient of the microcavity. For a GaAs/Al(sub x)Ga(sub 1-x)As MQW structure sandwiched between a GaAs/AlAs DBR and vacuum, we performed numerical calculations of the probe absorption spectra for different parameters such as pump intensity, pump detuning, and cavity length. We find that the probe spectrum is strongly dependent on these parameters. In particular, we find that the combination of the cavity effect and the Autler-Townes effect results in a triplet in the optical spectrum of the MQW system. The optical absorption peak value and its location can be feasibly controlled by varying the pump intensity and detuning.

Liu, Ansheng↗
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↗
NASA NTRS2021

A Sulfur Dioxide Covariance-Based Retrieval Algorithm (COBRA): Application to TROPOMI Reveals New Emission Sources

Sensitive and accurate detection of sulfur dioxide (SO 2 ) from space is important for monitoring and estimating global sulfur emissions. Inspired by detection methods applied in the thermal infrared, we present here a new scheme to retrieve SO 2 columns from satellite observations of ultraviolet back-scattered radiances. The retrieval is based on a measurement error covariance matrix to fully represent the SO 2 -free radiance variability, so that the SO 2 slant column density is the only retrieved parameter of the algorithm. We demonstrate this approach, named COBRA, on measurements from the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Sentinel-5 Precursor (S-5P) satellite. We show that the method reduces significantly both the noise and biases present in the current TROPOMI operational DOAS SO 2 retrievals. The performance of this technique is also benchmarked against that of the Principal Component Algorithm (PCA) approach. We find that the quality of the data is similar and even slightly better with the proposed COBRA approach. The ability of the algorithm to retrieve SO 2 accurately is also further supported by comparison with ground-based observations. We illustrate the great sensitivity of the method with a high-resolution global SO 2 map, considering two and a half years of TROPOMI data. In addition to the known sources, we detect many new SO 2 emission hotspots worldwide. For the largest sources, we use the COBRA data to estimate SO 2 emission rates. Results are comparable to other recently published TROPOMI-based SO 2 emissions estimates, but the associated uncertainties are significantly lower than with the operational data. Next, for a limited number of weak sources, we demonstrate the potential of our data for quantifying SO 2 emissions with a detection limit of about 8 kt yr -1 , a factor of 4 better than the emissions derived from the Ozone Monitoring Instrument (OMI). We anticipate that the systematic use of our TROPOMI COBRA SO 2 column data set at a global scale will allow identifying and quantifying missing sources, and help improving SO 2 emission inventories.

SO2↗
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↗
NASA NTRS1974

Reaction formulation for radiation and scattering from plates, corner reflectors and dielectric-coated cylinders

The reaction concept is employed to formulate an integral equation for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders. The surface-current density on the conducting surface is expanded with subsectional bases. The dielectric layer is modeled with polarization currents radiating in free space. Maxwell's equation and the boundary conditions are employed to express the polarization-current distribution in terms of the surface-current density on the conducting surface. By enforcing reaction tests with an array of electric test sources, the moment method is employed to reduce the integral equation to a matrix equation. Inversion of the matrix equation yields the current distribution, and the scattered field is then obtained by integrating the current distribution. The theory, computer program and numerical results are presented for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders.

Wang, N. N.↗
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↗
NASA NTRS2003

Calculation of Free-Free Opacities

Free-free absorption is an important contribution to the opacity for radiation transport through hot materials Temperatures can be as high as several keV, such that it becomes a computational challenge to solve the Schrodinger equation efficiently for rapidly oscillating continuum functions for high angular momenta. Several groups\footnots, including ours, have studied the phase amplitude solution (PAS) of the Schrodinger equation, in which one solves equations for the wave function amplitude and phase, which are: smooth functions of the electron energy. It is also important to have an accurate Schroudinger benchmark for the development of the PAS method. We present results for dipole matrix elements, Gaunt factors, and cross sections for the absorption of radiation at various energies for Cs XIX at temperature=100 eV and density=0.187 g/cc for our newly developed PAS and Schrodinger benchmark.

Bhatia, A. K.↗