SEARCH · Engineering Papers

Results for “Density matrix method”

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.

124 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Electrolytic gold plating, stripping, and ion transport dynamics through a solid-state iodide perovskite

The pronounced electrochemical reactivity between halide perovskites and metal electrodes can introduce mobile extrinsic metal ions which can cause device instability or enable novel functionalities. Here we systematically investigate the kinetics of gold cation (Au + ) migration in indium tin oxide (ITO)/methylammonium lead triiodide (MAPbI 3 )/Au model devices under long-term potentiostatic biasing. Scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) analyses reveal that Au + ions, electrochemically generated at the Au anode, traverse the perovskite layer with diffusion coefficients on the order of 10 −11 to 10 −10 cm 2 s −1 and are subsequently reduced at the cathode as Au 0 clusters, resembling metal plating behavior in electrolytic cells and solid-state batteries during charging. Furthermore, reversing the applied bias strips the plated Au 0 , revealing reversibility suitable for bipolar resistive switching devices and providing direct evidence of the electrochemical and ionic nature of Au transport within the perovskite matrix. Quantitatively determining diffusion coefficients and ion concentrations provides foundational inputs for future drift-diffusion modelling opportunities and allows us to relate our findings to implications on long term operation of devices like photovoltaic modules. These results clearly demonstrate the solid-state electrochemical nature of perovskite devices, highlight methods to be more quantitative about ion transport properties, provide and emphasize the importance of disentangling electro-, photo-, photoelectrochemical processes for understanding device performance and unlocking new functionalities.

14 SOLAR ENERGY↗

DOE OSTI2024

In-Plane Liftout and Push-to-Pull for In Situ Mechanical Testing of Irradiated Inconel X-750

A streamlined sample preparation method for nanomechanical testing is needed to improve the quality of specimens, reduce the cost, and increase the versatility of specimen fabrication. This work outlines an in-plane liftout focused ion beam (FIB) fabrication procedure to prepare electron-transparent specimens for in situ transmission electron microscopy (TEM) nanomechanical testing. Ion etching and electron backscatter diffraction (EBSD) techniques were used to lift out a [110] oriented grain from a neutron-irradiated bulk X-750 alloy. Careful control of voltages and currents ensured precision. Top surface thinning sweeps prevented resurfacing and redeposition while dog-bone geometries were shaped with a 1:4 gauge width-to-milling pattern diameter ratio. Nanotensile testing in the TEM with a picoindenter allowed for the estimation of an ultimate tensile strength of 2.41 GPa, and inspection revealed a high density of bubbles in the X-750 matrix. The proposed fabrication procedure is significant for preparing samples from radioactive materials, studying complex structures that are orientation-dependent, and analyzing desired planar areas.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Quantum-classical embedding via ghost Gutzwiller approximation for enhanced simulations of correlated electron systems

Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. Here, this work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.

Chen, I-Chi [Ames Laboratory (AMES), Ames, IA (Uni↗

DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Extracting off-diagonal order from diagonal basis measurements

Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multiscale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-dimensional extended Hubbard model with an off-diagonal order parameter, sandwiched between diagonal charge and spin density wave phases, using only diagonal descriptors. We study the model directly in the thermodynamic limit using the recently developed variational uniform matrix product states algorithm, and draw our samples from degenerate ground states related by global spin rotations, emulating the projective measurements that are accessible in experiments. Our results will have important implications for the study of exotic phases using optical lattice experiments. Published by the American Physical Society 2024

1-dimensional systems↗

DOE OSTI2024

Fiber Sorbents – A Versatile Platform for Sorption-Based Gas Separations

Increasing demand for high-purity fine chemicals and a drive for process intensification of large-scale separations have driven significant work on the development of highly engineered porous materials with promise for sorption-based separations. While sorptive separations in porous materials offer energy-efficient alternatives to longstanding thermal-based methods, the particulate nature of many of these sorbents has sometimes limited their large-scale deployment in high-throughput applications such as gas separations, for which the necessary high feed flow rates and gas velocities accrue prohibitive operational costs. These processability limitations have been historically addressed through powder shaping methods aimed at the fabrication of structured sorbent contactors based on pellets, beads or monoliths, commonly obtained as extrudates. These structures overcome limitations such as elevated pressure drops commonly recorded across powder adsorption beds but often accrue thermal limitations arising from elevated particle density and aggregation, which ultimately cap their maximum separation performance. Furthermore, the harsh mechanical strain to which powder particles are subjected during contactor fabrication, in the form of extrusion/compression forces, can result in partial pore occlusion and framework degradation, further limiting their performance. Here, we present the development of porous fiber sorbents as an alternative sorbent contactor design capable of addressing sorbent processability limitations while enabling an array of performance-maximizing heat integration capabilities. This new sorbent form factor leverages pre-existing know-how from hollow fiber spinning to produce fiber-shaped sorbent contactors through the phase inversion of known polymers in a process known as dry-jet/wet quenching. The process of phase inversion allows microporous sorbent particles to be latched onto a macroporous polymer matrix under mild processing conditions, thus making it compatible with soft porous materials prone to amorphization under traditional pelletization conditions. Sorbent fibers can be created with different geometries through control of the spinning apparatus and process, offering the possibility to produce monolithic and hollow fibers alike, the latter of which can be integrated with thermalization fluid flows. In this Account, we summarize our progress in the field of fiber sorbents from both design and application standpoints. We further guide the reader through the evolution of this field from the early inceptive work on zeolite hollow fibers to recent developments on MOF fibers. We highlight the versatile nature of fiber sorbents, both from the composition, fabrication and structure points of view, and further demonstrate how fiber sorbents offer alternative paths in tackling new and challenging chemical separation challenges like direct air capture (DAC), with a final perspective on the future of the field.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗

DOE OSTI2026

Testing of High S Matrix Glasses to Expand DFHLW Glass Compositional Ranges (Rev.1)

Gaps in glass composition-property data for direct-feed high-level waste (DFHLW) have recently been identified. One such gap is the region of high sulfur solubility since previous, pretreated, high-level wastes contained very little sulfur. Filling this data gap will significantly broaden the range of process flowsheet options including minimal washing and will allow for optimized waste loading in DFHLW glasses. This report summarizes the data collected during the characterization of the DFHLW High S Glass Matrix (HS24). A glass matrix of 50 glass compositions was developed to evenly cover the DFHLW composition region for high sulfur glass. Matrix glasses were designed to expand the composition region outside the current component concentration and property limits so as to reduce uncertainties at the limits. The 50 matrix glasses were fabricated and tested for properties important to the success of DFHLW vitrification including: compositions, canister centerline cooling (CCC) crystallinity and isothermal crystallinity, density, viscosity, electrical conductivity (EC), product consistency test (PCT) response, toxicity, and sulfate solubility. Melter materials corrosion testing is reported elsewhere. These glasses were intentionally designed to have high SO 3 solubilities (0.7 to 2.2 SO 3 wt%) in compositional regions that had not been previously explored. Forty-eight glasses showed the measured SO 3 content retained >80% of the target SO 3 and the densities of all the glasses ranged from 2.49 g·cm -3 to 2.74 g·cm -3 . While the model predicted nepheline formation in 5 glasses, one of the tested 50 CCC glasses formed nepheline, and 35 glasses formed Cr-containing phases such as spinels and eskolaite. Only five glasses were amorphous after CCC treatment where 44 glasses with detectable crystals contained =10 wt% crystals and only one glass had > 10 wt% crystals. None of the glasses exceeded the allowable T 2% for spinel crystal formation at 950 ºC (i.e., no glasses had >2 wt% spinel at 950 ºC) during isothermal crystal fraction tests where 10 glasses showed no crystalline phases at or below 950 ºC. All the glasses (except one which failed being slightly lower than the target) satisfied the SO 3 constraint while 98 glasses did not meet the viscosity constraints and 4 failed the EC constraints. Six quenched (Q) and six CCC glasses failed the Defense Waste Processing Facility (DWPF) Environmental Assessment (EA) glass PCT threshold and 3 Q and 4 CCC failed the PCT design constraint. One glass exceeded the WTP delisting limits for Cr via EPA Method 1311 (i.e., Toxicity Characteristic Leaching Procedure, TCLP). It should be emphasized that some of these glasses were specifically designed to approach or even exceed certain property constraints, as filling data gaps in these regions will provide the greatest benefit for future model development by improving accuracy and reducing uncertainties. These insights will ultimately support the development of more robust glass formulation strategies, enabling higher waste loading, reducing operational risks, and expanding the processing envelope.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Thermal nature of confining strings

We investigate the quantum statistical properties of the confining string connecting a static fermion-antifermion pair in the massive Schwinger model. By analyzing the reduced density matrix of the subsystem located in between the fermion and antifermion, we demonstrate that as the interfermion separation approaches the string-breaking distance, the overlap between the microscopic density matrix and an effective thermal density matrix exhibits a pronounced, narrow peak, approaching unity at the onset of string breaking. This behavior reveals that the confining flux tube evolves toward a genuinely thermal state as the separation between the charges grows, even in the absence of an external heat bath. In other words, one cannot tell whether a reduced state of the subsystem arises from a surrounding heat bath or from entanglement with the rest of the system. The entanglement spectrum near the critical string-breaking distance exhibits a rapid transition from the dominance of a single state describing the confining electric string towards a strongly entangled state containing virtual fermion-antifermion pairs. Our findings establish a quantitative link between confinement, entanglement, and emergent thermality, and suggest that string breaking corresponds to a microscopic thermalization transition within the flux tube.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Second harmonic generation in topological insulators under quantizing magnetic fields

We theoretically investigate the second harmonic generation (SHG) of topological insulator surface states in a perpendicular magnetic field. Our theory is based on the microscopic expression of the second-order magneto-optical conductivity developed from the density matrix formalism, taking into account hexagonal warping effects on the surface states' band structure. Using numerically exact Landau-level energies and wave functions including hexagonal warping, we calculate the spectrum of SHG conductivities under normal incidence for different values of magnetic field and chemical potential. The imaginary parts of the SHG conductivities show prominent resonant peaks corresponding to one-photon and two-photon inter-Landau-level transitions. Treating the hexagonal warping term perturbatively, these transitions are clarified analytically within a perturbation theory from which approximate selection rules for the allowable optical transitions for SHG are determined. Furthermore, our results show extremely high SHG susceptibility that is easily tunable by magnetic field and doping level for topological surface states in the far-infrared regime, exceeding that of many conventional nonlinear materials. This work highlights the key role of hexagonal warping effects in generating second-order optical responses and provides insights into the nonlinear magneto-optical properties of the topological insulators.

Landau levels↗

DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗

DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗