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 361 records · Page 20

Bound Constrained Partial DifferentialEquation Inverse Problem Solution by theSemi-Smooth Newton Method

We present the mathematical derivation, software implementation details, and computational results for a semi-smooth Newton method applied to two inverse problems governed by partial differential equations with bound constraints. The two problems share mathematical structural similarities to density-based topology optimization problems. The semi-smooth Newton method provides a mesh independent solution computation for the two test problems. A key step is that the complementarity part of the necessary optimality conditions are reformulated with the use of a complementarity functionφsuch that the complementarity conditions are satisfied if and only if a zero of a nonsmooth function has been obtained. The modular finite element package MFEM is utilized for the software implementation. In addition we constructed a matrix-free Operator to enable the use of efficient Krylov subspace IterativeSolver of MFEM for the solution of our two target problems.

97 MATHEMATICS AND COMPUTING↗

Microstructures and mechanical properties of a modified 9Cr ferritic-martensitic steel in the as-built condition after additive manufacturing

A newly developed nano-structured high-Mn 9Cr ferritic-martensitic (FM) steel designed for additive manufacturing (Additive-manufactured Nanostructured Alloy, ANA) was fabricated via the direct energy deposition method. The as-printed ANA2 showed tensile properties (strength and elongation) and Charpy impact toughness (upper-shelf energy and ductile-brittle transition temperature) comparable to oxide-dispersion-strengthened (ODS) Eurofer and PM2000. Compared with wrought P91, the as-built ANA2 showed significantly higher yield strength but inferior Charpy impact toughness and ductility. Microstructural characterization revealed uniformly distributed sub-micron pores, an extremely high dislocation density, and a high density of ultrafine nano-structured precipitates in the matrix. The microstructures indicated the dislocations as the dominant strengthening element and the ultrafine precipitates as the primary contribution to estimated irradiation sink density. This work demonstrates the ability to exploit the unique consolidation properties of additive manufacturing to fabricate steels which marry the benefits of both ODS alloys and conventionally processed FM steels.

36 MATERIALS SCIENCE↗

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗

Extracting the electronic structure of light elements in bulk materials through a Compton scattering method in the readily accessible hard x-ray regime

Here, our Compton profile measurements of Ti and TiH 2 using readily available hard X-ray radiation at 27.5 keV, detected by both a Hitachi Vortex silicon-drift detector and a high-resolution superconducting transition-edge sensor array, are found to be in excellent accord with state-of-the-art density functional theory based calculations. The spherically averaged difference between the Compton profiles of TiH 2 and Ti is well described by an inverted parabola, supporting an itinerant behavior of the electron gas screening the protons in the Ti matrix. Our experimental approach, validated by two different detectors, extends the applicability of Compton scattering technique to the readily accessible hard x-ray regime (below 30 keV). Our study suggests possibilities for experiments at low-flux bending magnet synchrotron beamlines and paves the way for the development of tabletop Compton experiments with x-ray tubes.

36 MATERIALS SCIENCE↗

Microstructural and nanomechanical studies of PVD Cr coatings on SiC for LWR fuel cladding applications

The microstructure and nanomechanical properties of pure Cr coatings deposited on SiC with various physical vapor deposition (PVD) methods have been investigated for the mitigation of hydrothermal corrosion of SiC-SiC f (SiC fiber-reinforced SiC matrix) composite fuel cladding in light water reactors (LWR). Cr coatings (4–7 μm thick) were deposited on unbiased SiC substrates using six variants of magnetron sputtering processes: (i) standard (S-DCMS), (ii) pulsed (P-DCMS), (iii) ion-assisted (I-DCMS), and (iv) pulsed ion-assisted direct current magnetron sputtering (PI-DCMS), (v) high-power impulse magnetron sputtering (HiPIMS), and (vi) bipolar HiPIMS (B-HiPIMS). Microstructural characterization and nanoindentation testing were used to evaluate effects of deposition technique on coating microstructure and nanomechanical properties. The application of positive reverse pulses to the sputter target induced coarsening of the columnar grain structure, while ion bombardment of substrates during deposition promoted densification of the coating microstructure. Both effects are exemplified in B-HiPIMS deposition, resulting in a high-density microstructure with compressive residual stress. Material pile-up around nano-indents and plastic work during nanoindentation were analyzed to understand deformation behavior of the coatings. Finally, the results suggest the B-HiPIMS process to be among the promising methods for the deposition of Cr coatings on SiC-SiC f for LWR cladding application.

36 MATERIALS SCIENCE↗

A highly stable Cu(OH) 2 -Poly(vinyl alcohol) nanocomposite membrane for dramatically enhanced direct borohydride fuel cell performance

Nano-additive aggregation, limited performance, and unclear modification mechanisms are the main obstacles in developing nanocomposite anion exchange membranes (AEMs). In this work, for the first time, an effective and highly stable Poly(vinyl alcohol) (PVA)-based AEMs with dispersive Cu(OH)$_2$ nanoclusters (Cu-AEMs) are prepared by a simple and eco-friendly three-step method: ‘CuCl$_2$ doping-casting-KOH immersing’. The doped Cu$^{2+}$ ions chemically combine with OH- ions to form anionic conductive Cu(OH)2 nanoclusters intermediated by attaching resins. The PVA skeletons wrap around Cu(OH)$_2$ while the hydroxyl groups expose to bulk water, forming Cu(OH)$_2$-PVA complex, which avoids nano-additive aggregation, increases anionic channels, and strengthens additive-matrix connection. A direct borohydride fuel cell using Cu-AEM with 0.56 wt % CuCl$_2$ possess the highest power density of 403.3 mW cm$^{-2}$ at 60 °C and a life span of over 200 h. The high-performance and durability come from the unique structure of Cu(OH)$_2$ nanoclusters-PVA complex. ‘Vehicle’ theory is considered to be the dominant mechanism for enhancing such nanocomposite Cu-AEMs. This work demonstrates a new concept for preparing the stable AEMs toward high-performance fuel cells. The synthetic chemistry involved can be broadly extended for fabricating versatile AEMs.

25 ENERGY STORAGE↗

Alpha Shape Analysis (ASA) Framework for Post- Clustering Property Determination in Atom Probe Tomographic Data

While application of clustering algorithms to atom probe tomography data have enabled quantification of solute clusters in terms of number density, size, and subcomposition there exist other properties (e.g., volume, surface area, and composition) that are better determined by defining an interface between the cluster and the surrounding matrix. The limitation in composition results from an ion selection step where the expected matrix ion types are omitted from the cluster search algorithm to enhance the contrast between the matrix and cluster and to reduce the complexity of the search. Previously, composition determination within solute clusters has utilized a secondary envelopment and erosion step on top of conventional methods such as maximum separation. In this work, we present a novel stochastic method that combines the particle identification fidelity of a conventional clustering algorithm with the analytical flexibility of mesh-based approaches through the generation of alpha shapes for each identified cluster. Furthermore, the corresponding mesh accounts for concave components of the clusters and determines the volume and surface area of the clusters; additionally, the mesh boundary is utilized to update the total composition according to the internal ions.

42 ENGINEERING↗

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗

Surface Functionalized Barium Titanate Nanoparticles: A Combined Experimental and Computational Study

Barium titanate (BTO) nanoparticles show great potential for use in electrostatic capacitors with high energy density. This includes both polymer composite and sintered capacitors. However, questions about the nanoparticles' size distribution, amount of agglomeration, and surface ligand effect on performance properties remain. Reducing particle agglomeration is a crucial step to understanding the properties of nanoscale particles, as agglomeration has significant effects on the composite dielectric constant. BTO surface functionalization using phosphonic acids is known reduce BTO nanoparticle agglomeration. We explore solution synthesized 10 nm BTO particles with tert-butylphosphonic acid ligands. Recent methods to quantifying agglomeration using an epoxy matrix before imaging shows that tert-butylphosphonic acid ligands reduce BTO agglomeration by 33%. Thermometric, spectroscopic, and computational methods provide confirmation of ligand binding and provide evidence of multiple ligand binding modes on the BTO particle surface.

36 MATERIALS SCIENCE↗

Effect of off-diagonal elements in the Wannier Hamiltonian on DFT + DMFT for low-symmetry materials: Study of Li 2 MnO 3

Here, we study the effect of the off-diagonal elements of the Wannier Hamiltonian on the electronic structure of the low-symmetry material Li 2 MnO 3 ( C2/m ), using dynamical mean field theory calculations with a continuous-time quantum Monte Carlo impurity solver. The presence of significant off-diagonal elements leads to a pronounced suppression of the energy gap. The off-diagonal elements are largest when the Wannier projection is used based on the global coordinate, and they remain substantial even with the projection using the local coordinate close to the direction of Mn-O bonds. We show that the energy gap is enhanced by the diagonalization of the Mn d block in the full p-d Hamiltonian with the application of a unitary rotation matrix. Additionally, the inclusion of small double counting energy is crucial for achieving the experimental gap by reducing p-d hybridization. Furthermore, we establish the efficiency of a low-energy (d-only basis) model for studying the electronic structure of Li 2 MnO 3 , as the Wannier basis represents a hybridized state of Mn d and O p orbitals. These findings suggest an appropriate approach for investigating low-symmetry materials using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. We also find that the antiferromagnetic ground state $\Gamma$ 2u is stable with U ≤ 2 eV within density functional theory+$U$ calculations, which is much smaller than the widely used U = 5 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structural response reconstruction using a system-equivalent singular vector basis

Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

42 ENGINEERING↗

Thermal effects on the electronic properties of sodium electride under high pressures

Sodium is a simple metal at ambient conditions, while it transits to an electride phase at pressures above ~160 GPa along the room temperature isotherm. Here we explore the thermal effects on the electronic properties of the $hP4$ phase of sodium along the $ρ$ = 5.872 g/cm 3 isochore. We quantitatively classify this phase as an insulator based on the criterion of nearsightedness of the one-particle density matrix. Ab initio calculations suggest that the band gap of this insulator decreases with increasing temperature along the isochore, primarily because of ionic distortions; culminating in an insulator-to-metal transition upon melting at $T_m \approx 2100$ $\hbox{K}$. This transition is accompanied by residual electronic localization (in d orbitals) in the form of dynamic electron bubbles and a change in hybridization form $p–d$ to $s–p$ upon melting. This transition is explored by tracking the electronic and electro-optical properties along the isochore under consideration.

36 MATERIALS SCIENCE↗

Conformal collider physics meets LHC data

The remarkably high energies of the Large Hadron Collider (LHC) have allowed for the first measurements of the shapes and scalings of multipoint correlators of energy flow operators, ⟨ Ψ | E ( n → 1 ) E ( n → 2 ) ⋯ E ( n → k ) | Ψ ⟩ , providing new insights into the Lorentzian dynamics of quantum chromodynamics (QCD). In this letter, we use recent advances in effective field theory to derive a rigorous factorization theorem for the light-ray density matrix, ρ = | Ψ ⟩ ⟨ Ψ | , inside high transverse momentum jets at the LHC. Using the light-ray operator product expansion, the scaling behavior of multipoint correlators can be computed from the expectation value of the twist-2 spin- J light-ray operators, O [ J ] , in this state, Tr [ ρ O [ J ] ] . We compute the light-ray density matrix at next-to-leading order, and combine this with results for the next-to-leading logarithmic scaling behavior of the correlators up to six-points, comparing with CMS open data. This theoretical accuracy allows us to resolve the quantum scaling dimensions of QCD light-ray operators inside jets at the LHC. Our factorization theorem for the light-ray density matrix at the LHC completes the link between recent developments in the study of energy correlators and LHC phenomenology, opening the door to a wide variety of precision jet substructure studies. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Soft X-ray tomography reveals variations in B. subtilis biofilm structure upon tasA deletion

Bacterial biofilms are complex cell communities within a self-produced extracellular matrix, crucial in various fields but challenging to analyze in 3D. We developed a "biofilm-in-capillary" growth method compatible with full-rotation soft X-ray tomography, enabling high-resolution 3D imaging of bacterial cells and their matrix during biofilm formation. This approach offers 50 nm isotropic spatial resolution, rapid imaging, and quantitative native analysis of biofilm structure. Using Bacillus subtilis biofilms, we detected coherent alignment and chaining of wild-type cells towards the oxygen-rich capillary tip. In contrast, the ΔtasA genetic knock-out showed a loss of cellular orientation and changes in the extracellular matrix. Adding TasA protein to the ΔtasA strain restored matrix density and led to cell assembly compaction, but without the chaining observed in wild-type biofilms. This scalable and transferable approach opens new avenues for examining biofilm structure and function across various species, including mixed biofilms, and response to genetic and environmental factors.

59 BASIC BIOLOGICAL SCIENCES↗

Investigation of Ga 2 O 3 as a new transparent conductive oxide for photovoltaics applications

This small innovation project intended to leverage recent activity and advances in ultrawide-bandgap oxide power electronics toward the investigation of Ga 2 O 3 as a transparent conductive oxide (TCO) for use in photovoltaics (PV). Ga 2 O 3 ’s theoretical advantage over incumbent TCOs is its very large bandgap of 4.8 eV, ensuring optical transparency of photons with λ ≥ 260 nm, effectively the full terrestrial solar spectrum. At the time of proposal writing, literature on Ga 2 O 3 as a TCO/optical material was relatively sparse — what little existed was focused on for UV sensors and/or emitters, with none related to PV. As such, this project was intended to help determine the potential of Ga 2 O 3 as a PV-oriented TCO by investigating its deposition using tools common to the PV and electronics industries — atomic layer deposition (ALD) and RF sputtering (RS), both of which are already used in PV manufacturing — and its resultant optical and electronic properties. Final project goals were to test its application to Si and III-V solar cells. Controllable film thickness and excellent uniformity was established for both methods deposition methods, with ALD providing higher precision for thinner films and RS more effective for growing thicker films. The as-deposited films were found to be amorphous in nature. Spectroscopic ellipsometry (SE) confirmed bandgaps of at least 4.8 eV and non-parasitic absorption across both AM0 and AM1.5G spectra; a refractive index approaching the expected value of 1.8 was also observed, with some degree of tunability based on process parameters. Using this initial optical data, a transfer matrix model, which interfaces with in-house EQE and LIV models, was developed to simulate Ga 2 O 3 optical effects on various solar cells, included potential in antireflection coatings (ARCs). Despite promising optical properties, the resultant resistivity / conductivity metrics did not meet expectations, regardless of deposition process. Although thick RS-deposited films, using both undoped and 1 at% Ge (n-type) doped sintered targets — demonstrated high net carrier concentrations (via C-V), and low specific resistance Ohmic contacts were demonstrated, transmission line measurements (TLM) showed very high resistivity. Further analysis indicated high vertical film conductivity, but very low lateral conductivity. Deeper characterization revealed a relatively high density of nano/microcrystalline inclusions that were ostensibly the source of the vertical conductivity, with the amorphous matrix serving as an effective insulator (likely due to high concentrations of electronic trap states). Despite extensive work to increase film polycrystallinity — demonstrated via high temperature annealing — sufficient lateral conductivity for TCO use was not achieved. The final phase of the project shifted focus toward investigation of ALD-deposited Ga 2 O 3 ’s potential as a passivant and/or passivating contact for both Si and III-V solar cells, with expectation of performance similar to Al 2 O 3 . However, initial rounds of testing using our baseline ALD process yielded minimal (but not quite zero) passivation of both GaInP and Si surfaces. Dielectric passivation is well-known to be highly process sensitive; given the unoptimized nature of the process used, this work was deemed inconclusive. The final outlook with respect to feasibility of Ga 2 O 3 as a PV-oriented TCO is, ultimately, still uncertain. The work performed in this project confirmed the optical properties and deposition methods, but the achieved electrical properties do not meet technological needs; literature reports in recent years are somewhat inconsistent and potentially untrustworthy, but generally appear to be in line with our results. At the very least it is clear that, should Ga 2 O 3 still be under consideration for PV TCO and/or selective contact applications, a significant amount of optimization and study is still needed.

14 SOLAR ENERGY↗

Entanglement and many-body effects in collective neutrino oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings, such as supernovae, where the neutrino density is large. In this regime, neutrino-neutrino interactions are important and simulations in the mean-field approximation show evidence for collective oscillations occurring at timescales much shorter than those associated with vacuum oscillations. In this work, we study the out-of-equilibrium dynamics of a corresponding spin model using matrix product states and show how collective bipolar oscillations can be triggered by many-body correlations if appropriate initial conditions are present. We find entanglement entropies scaling at most logarithmically in the system size suggesting that classical tensor network methods could be efficient in describing collective neutrino dynamics more generally. These observation provide a clear path forward, not only to increase the accuracy of current simulations, but also to elucidate the mechanism behind collective flavor oscillations without resorting to the mean-field approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗