Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY MATRIX”

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 577 records · Page 32

Higher Hall conductivity from a single wave function: Obstructions to symmetry-preserving gapped edge of (2+1)-dimensional topological order

A (2+1)D topologically ordered phase with U(1) symmetry may or may not have a symmetric gapped edge state, even if both thermal and electric Hall conductivity are vanishing. It has recently been discovered that there are “higher” versions of Hall conductivity valid for fermionic fractional quantum Hall (FQH) states that obstruct symmetry-preserving gapped edge states beyond thermal and electric Hall conductivity. In this paper, we show that one can extract higher Hall conductivity from a single wave function of an FQH state, by evaluating the expectation value of the “partial rotation” unitary, which is a combination of partial spatial rotation and a U(1) phase rotation. This result is verified numerically with the fermionic Laughlin state with 𝜈=1/3 and 1/5, as well as the non-Abelian Moore-Read state. Together with topological entanglement entropy, we prove that the expectation values of the partial rotation completely determine if a bosonic/fermionic Abelian topological order with U(1) symmetry has a symmetry-preserving gappable edge state or not. We also show that thermal and electric Hall conductivity of Abelian topological order can be extracted by partial rotations. Even in non-Abelian FQH states, partial rotation provides the Lieb-Schultz-Mattis type theorem constraining the low-energy spectrum of the bulk-boundary system. The generalization of higher Hall conductivity to the case with Lie group symmetry is also presented.

2-dimensional systems↗

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗

A retrospective review of von Neumann’s analysis of hidden variables in quantum mechanics

This article reviews the history of J. von Neumann’s analysis of hidden variables in quantum mechanics and the subsequent analysis by others. In his book The Mathematical Foundations of Quantum Mechanics , published in 1932, von Neumann performed an analysis of the consequences of introducing hidden parameters (hidden variables) into quantum mechanics. He arrived at two principal conclusions: first, hidden variables cannot be incorporated into the existing theory of quantum mechanics without major modifications, and second, if they did exist, the theory would have already failed in situations where it has been successfully applied. This analysis has been taken as an “incorrect proof” against the existence of hidden variables, possibly due to a mistranslation of the German word prufen . von Neumann’s so-called proof isn’t even wrong as such a proof does not exist, but it is an examination of the limitations imposed by internal consistency of the Hilbert space formulation of the theory. One of the earliest attempts to eliminate uncertainty, by D. Bohm, requires a major modification of quantum mechanics (observables are not represented by Hermitian operators), which supports von Neumann’s first principal conclusion. However, testing the Bohm theory requires constructing a physically impossible initial state. As such, the theory has no experimental consequences, so W. Pauli referred to it as an “uncashable check”. As there are no observable consequences, the Bohm theory is possibly a counterexample to von Neumann’s second conclusion that hidden variables in particular would have already led to a failure of the theory.

density matrix↗

Gamma ray Imaging Based Methods to Improve the Accuracy of Uranium Holdup Quantification

Deposits in commercial uranium processing facilities of interest to nuclear safeguards are predominantly low enriched uranium (LEU). The shapes and sizes of holdup deposits can vary significantly, thus making it difficult to model the deposits accurately for calibration purposes. This presents a challenge to currently employed methods such as the Generalized Geometry Holdup (GGH) that rely on quantifying 235U mass by simplifying deposit shapes as a point, a line, or an area. In this work, gamma ray imaging using high energy resolution Germanium Gamma Imagers (GeGI™) is employed to determine the distribution of uranium inside the source containment. Two types of imaging methods are employed; coded-aperture imaging and Compton imaging. Gamma-ray emissions seen in different regions of the image are quantified using an inverse gamma-ray transport solver being developed at ORNL. Using measured data, the inverse method will be used to solve for unknown source parameters such as source matrix thickness, density, and attenuation due to container wall and shielding. The intrinsic efficiency of the detector is determined based on a library of built-in response functions and key dimensions of the detector (e.g., thickness and radius of the detector crystal). Measurements are performed in both imaging modalities using uranium sources of well-known masses and enrichments configured inside mocked-up holdup fixtures. Results from the two imaging modalities are intercompared. The paper presents preliminary results of the imaging measurements and the progress made in the computational methods to quantify uranium mass.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Two-stage combustor

A two-stage combustor having as constituent parts: a partial oxidation reactor, which catalytically converts a hydrocarbon fuel and a first supply of oxidant into a gaseous partial oxidation product; and a deep oxidation reactor having a premixer plenum fluidly connected to a porous heat spreader, which converts the gaseous partial oxidation product to deep oxidation products. In one embodiment, the premixer plenum provides an empty space wherein combustion occurs in flame mode. In a second embodiment, the premixer plenum contains a high pore density foam matrix, absent catalyst, which facilitates holding a flameless combustion downstream within the porous heat spreader. In both embodiments heat produced during combustion is transmitted from the heat spreader to an associated heat acceptor, such as a heater head of a Stirling engine.

Roychoudhury, Subir↗

Localized Phonon Densities of States at Grain Boundaries in Silicon

Since it is now possible to record vibrational spectra at nanometer scales in the electron microscope, it is of interest to explore whether extended defects in crystals such as dislocations or grain boundaries will result in measurable changes of the phonon densities of states (dos) that are reflected in the spectra. Phonon densities of states were calculated for a set of high angle grain boundaries in silicon. The boundaries are modeled by supercells with up to 160 atoms, and the vibrational densities of states were calculated by taking the Fourier transform of the velocity–velocity autocorrelation function from molecular dynamics simulations with larger supercells doubled in all three directions. In selected cases, the results were checked on the original supercells by comparison with the densities of states obtained by diagonalizing the dynamical matrix calculated using density functional theory. Near the core of the grain boundary, the height of the optic phonon peak in the dos at 60 meV was suppressed relative to features due to acoustic phonons that are largely unchanged relative to their bulk values. This can be attributed to the variation in the strength of bonds in grain boundary core regions where there is a range of bond lengths.

Materials Science↗

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↗

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)↗

Machine Learning–Based Condition Monitoring of a Circulating Water System of a Canadian Nuclear Plant

With the need to maintain long-term reliable energy using nuclear power plants, there is an underlying demand to ensure that the maintenance of plant components and systems is also done in an efficient and cost-effective manner. One way to achieve this is by moving from time-based maintenance to condition-based maintenance. The research presented in this paper focuses on applying statistical and machine-learning-based methods to capture anomalies within data for fault detection to further develop into condition monitoring. This paper focuses on system data for a circulating water system (CWS) of a pressurized heavy-water reactor for detecting anomalies. The different methodologies used for detecting and capturing anomalies in the CWS data are matrix profile, density-based spatial clustering of applications with noise (DBSCAN), and support vector machines (SVMs). Matrix profile and DBSCAN are used to distinguish between normal data and anomalous data. This paper presents a hybrid method using DBSCAN and SVM when a portion of the data is used for DBSCAN to generate clusters. This portion of data is then used to train the SVM along with the clusters generated by DBSCAN as output. SVM is then tested on unseen data as a predictive tool, which can work in real time to categorize data points as either normal or anomalous. This paper presents results that show the high accuracies of DBSCAN and SVM in capturing anomalies within the data for a CWS for fault detection. Thus, the maintenance plan would be focused on component condition rather than a time-based schedule by switching to an automated system to identify and predict faults within a CWS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High Critical Current Density in Internally-oxidized Nb 3 Sn Superconductors and its Origin

In this paper, we report achievement of non-Cu-matrix critical current density $(J_c)$ in internally-oxidized Nb 3 Sn superconductors that surpasses the best state-of-the-art Nb 3 Sn and the $J_c$ specification of the Future Circular Collider. We then explore the mechanism of the improvement, which was believed to be caused by refined grain size. However, our experiment showed that even with similar grain size, internally-oxidized wires had significantly higher $J_c$ than non-oxidized wires. This, combined with transmission electron microscopy and atom probe tomography studies which showed a high density of intragranular nano-precipitates with suitable sizes as pinning centers in internally-oxidized wires, demonstrates a significant contribution from the nano-precipitates.

43 PARTICLE ACCELERATORS↗

High susceptibility to adiabatic shear banding and high dynamic strength in tungsten heavy alloys with a high-entropy alloy matrix

The susceptibility to adiabatic shear banding (ASB) of tungsten heavy alloys (WHAs) is a key factor to determine their penetration performance for kinetic energy penetrators (KEPs). This work reports a novel WHA with strong susceptibility to ASB that uses a precipitation-hardening high-entropy alloy (HEA) as the binding phase of W particles. The results find that nanoscale L12 precipitates and high dislocation density in the matrix contribute a high dynamic yield strength of ~2300 MPa. Also, a narrow ASB with a width of ~12 μm forms in the novel WHA, which is only 6% of the width of ASB in traditional WHAs like W-NiFe, exhibiting a great susceptibility to ASB. The redissolution of precipitates into the matrix because of temperature rise during shearing could sharply soften the shearing location, which is suggested to be a very important factor of destabilizing mechanism to promote the development of ASB.

36 MATERIALS SCIENCE↗

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↗

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of Cooling Rates on Self-Reinforced Polyethylene Composites via Multiscale Experimental Characterizations

This study investigates cooling rate effects on self-reinforced polyethylene composites (PE-SRCs), focusing on macroscopic thermal and mechanical properties, and micro- to sub-micron morphologies. Unlike conventional fiber-reinforced composites (e.g., carbon and glass fiber based), cooling rate effects on SRCs remain less studied. PE-SRCs were fabricated from ultrahigh molecular weight polyethylene (UHMWPE) fabric and high-density polyethylene (HDPE) matrix with cooling rates from 0.76°C/min to 425°C/min. Tensile properties, including modulus (~9.5 GPa), strength (~360 MPa), and ultimate elongation (~6%), showed no significant cooling rate dependence. However, thermal analysis revealed melting temperature of HDPE matrix decreased from ~130°C to ~126.5°C with total crystallinity reducing from 84% to 75%. Wide-angle X-ray scattering confirmed decreasing crystallinity but at lower values (67%–59%), as thermal analysis might overestimate crystallinity. Small-angle X-ray scattering showed long period decreased from 24.3 to 17.4 nm, attributed to thinner lamellae and closer packing, suggesting limited crystal growth and faster nucleation. While cooling rates alter melting point and crystalline structures of HDPE matrix, these changes do not translate to substantial impact on macroscale mechanical properties of PE-SRCs. This insensitivity might be related to low interfacial shear strength at PE/PE interfaces, enabling process optimization while maintaining mechanical performance.

Cooling rate↗

Hemp fiber reinforced dual dynamic network vitrimer biocomposites with direct incorporation of amino silane

Natural fiber composites are inexpensive and renewable alternatives to traditional fiber reinforced polymers (FRPs). However, existing natural fiber composites primarily rely on petrochemical based thermosetting matrices, which are difficult to recycle due to their stable crosslinked network structures. To address this challenge, it is desirable to develop natural fiber composites with inherently recyclable biobased matrices. In this study, we developed a dual dynamic network vitrimer matrix from hempseed oil and limonene derivatives and demonstrated its application in hemp fiber reinforced composites. To improve the weak interface between hemp fiber and the polymer matrix, a common challenge in the realm of natural fiber composites, we directly incorporated amino silane into the vitrimer matrix. The amino silane participated in the polymer network, increasing the crosslink density and toughening the matrix, as evidenced by a significantly improved impact strength from 3.5 kJ/m 2 to 10.3 kJ/m 2 . Moreover, the incorporation of amino silane resulted in a lower water absorption by 11% in 7-day soaking for the composites, demonstrating improved fiber/matrix interfacial interaction. Furthermore, our biobased vitrimer matrix exhibits both imine and hydroxy-ester dynamic bonds, which enable the recycling of the biocomposite through a mild and cost-effective aminolysis process (100°C, 3h, ambient pressure). The decomposed polymer matrix was successfully reused as a polyol for polyurethane adhesives, and the surface morphology of recovered hemp fibers was analyzed and compared with that of the original fibers. Furthermore, these findings will help broaden the use of vitrimers for practical natural fiber composite applications, raise awareness of the problem of recycling natural fiber composite waste, and shed light on the interfacial challenge of natural fiber reinforced vitrimer composites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of processing temperature, pressure, and fiber volume fraction on mechanical and morphological behaviors of fully-recyclable uni-directional thermoplastic polymer-fiber-reinforced polymers

This work explores a type of composite called thermoplastic polymer-fiber-reinforced polymers (PFRPs), often referred to as self-reinforced composites (SRCs). A representative PFRP was exemplified using unidirectional (UD) ultra-high-molecular-weight polyethylene (UHMWPE) fibers embedded in a high-density polyethylene (HDPE) matrix. The effects of compression molding temperature and pressure on the mechanical and morphological behaviors of the filament-wound PFRPs with various fiber volume fractions (V f ) were experimentally investigated. The results elucidate the evolution of morphologies and tensile properties of the PFRPs due to thermal melting, fiber misalignment from pressure, and (V f )-induced structural variance, which has not been comprehensively reported yet. The highest specific tensile strength and modulus of the PFRP laminae reach 600 MPa/(g/cm 3 ) and 31 GPa/(g/cm 3 ), respectively. These properties are comparable to glass-/aramid-fiber-reinforced polymers (GFRPs, GFRTPs, AFRPs, and AFRTPs), with PFRPs exhibiting better ductility (specific strain at peak load ≈ 4%/(g/cm 3 )) than other common polymer composites. The motivation for this work was the high recyclability of PFRPs, which can be recycled by melting both the fibers and the matrix, and then reshaped them for re-manufacturing composites to maximize the efficiency in material reuse. This process simplifies the implementation of closed-loop recycling, re-manufacturing, and reuse to support sustainability in composites. This work aims to contribute to advancing thermoplastic PFRPs for their potential applications in various industries.

36 MATERIALS SCIENCE↗