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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

In-situ ion irradiation of fission products in a spent UO 2 fuel

This study investigates the behavior of fission gas bubbles and five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) in spent uranium dioxide (UO 2 ) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles and 5MPs in UO 2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. In conclusion, this study provides critical insights into the dynamic behavior of fission products under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.

Fission products↗
DOE OSTI2020

Axions in neutron star mergers

Supernovae and cooling neutron stars have long been used to constrain the properties of axions, such as their mass and interactions with nucleons and other Standard Model particles. We investigate the prospects of using neutron star mergers as a similar location where axions can be probed in the future. We examine the impact axions would have on mergers, considering both the possibility that they free-stream through the dense nuclear matter and the case where they are trapped. We calculate the mean free path of axions in merger conditions, and find that they would free-stream through the merger in all thermodynamic conditions. In contrast to previous calculations, we integrate over the entire phase space while using a relativistic treatment of the nucleons, assuming the matrix element is momentum-independent. In particular, we use a relativistic mean field theory to describe the nucleons, taking into account the precipitous decrease in the effective mass of the nucleons as density increases above nuclear saturation density. We find that within current constraints on the axion-neutron coupling, axions could cool nuclear matter on timescales relevant to neutron star mergers. Our results may be regarded as first steps aimed at understanding how axions affect merger simulations and potentially interface with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Analysis of the neutron matter equation of state and the symmetry energy up to fourth order of chiral effective field theory

We present predictions for the neutron matter equation of state, from leading to fourth order of chiral effective field theory, using recently developed, accurate chiral nucleon-nucleon potentials. For the many-body method, we employ the nonperturbative particle-particle ladder approximation, that is, we solve the G-matrix equation. Furthermore, we find the impact of subleading three-neutron forces to be mild and attractive. We also show order-by-order predictions for the symmetry energy, and discuss its density dependence in relation to empirical constraints. For the nuclear matter equation of state, in this work we adopt an empirical parametrization with good saturation properties. This is to highlight, specifically, the energy and pressure in neutron matter, particularly when comparing with empirical constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

High-temperature lean Cu alloys with Cr-to-Nb atomic ratio of 2

Two Cu-Cr-Nb alloys, denoted as alloy 1 (comprising Cu-0.89 at% Cr-0.42 at% Nb) and alloy 2 (comprising Cu-1.84 at% Cr-0.99 at% Nb), were produced through a series of manufacturing processes including vacuum induction melting, melt spinning, consolidation, brazing, and baking, with both alloys aimed at achieving a nominal Cr-to-Nb atomic ratio of 2. Microstructural characterization using transmission electron microscopy and X-ray diffraction identified the cubic C15 Laves-phase Cr 2 Nb as the dominant precipitate in both alloys, cross-validated by thermodynamic calculations and atomistic simulation-based density functional theory (DFT). Besides cubic C15 Cr 2 Nb, hexagonal C14-phase Cr 2 Nb and α-BiF3 cubic structured Cr 3 Nb were also observed in the alloys, including a coherent interface formed between the Cr 3 Nb precipitate and the Cu matrix. The hardness of the alloys increases, and the electrical conductivity decreases with increasing alloying addition content; two practical equations described the trends. Further DFT simulations revealed that the electrical conductivity (conductance) of the Cu/Cr 2 Nb interface is an order of magnitude higher than the intrinsic Cu high-angle grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Coherent gluon radiation: beyond leading-log accuracy

Abstract Results are presented for the medium-induced, soft coherent radiation spectrum for all 2 → 2 partonic channels in QCD, at leading-order inα s but beyond leading logarithmic accuracy. The general formula is valid in the full kinematic range of the underlying process, and reduces to previous results in special cases. The soft gluon radiation spectrum is expressed in terms of thecolor density matrixspecific to each channel, quantifying the entanglement between the color components of the 2 → 2 production amplitude. Beyond the leading logarithm, the spectrum depends explicitly on the off-diagonal elements of this matrix, owing to the soft gluon’s ability to probe the internal color structure of the parton pair.

Physics↗
DOE OSTI2024

Theoretical Prediction of Thermal Expansion Anisotropy for Y 2 Si 2 O 7 Environmental Barrier Coatings Using a Deep Neural Network Potential and Comparison to Experiment

Environmental barrier coatings (EBCs) are an enabling technology for silicon carbide (SiC)-based ceramic matrix composites (CMCs) in extreme environments such as gas turbine engines. However, the development of new coating systems is hindered by the large design space and difficulty in predicting the properties for these materials. Density Functional Theory (DFT) has successfully been used to model and predict some thermodynamic and thermo-mechanical properties of high-temperature ceramics for EBCs, although these calculations are challenging due to their high computational costs. In this work, we use machine learning to train a deep neural network potential (DNP) for Y 2 Si 2 O 7 , which is then applied to calculate the thermodynamic and thermo-mechanical properties at near-DFT accuracy much faster and using less computational resources than DFT. We use this DNP to predict the phonon-based thermodynamic properties of Y 2 Si 2 O 7 with good agreement to DFT and experiments. We also utilize the DNP to calculate the anisotropic, lattice direction-dependent coefficients of thermal expansion (CTEs) for Y 2 Si 2 O 7 . Molecular dynamics trajectories using the DNP correctly demonstrate the accurate prediction of the anisotropy of the CTE in good agreement with the diffraction experiments. In the future, this DNP could be applied to accelerate additional property calculations for Y 2 Si 2 O 7 compared to DFT or experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

0D-1D Hybrid Silicon Nanocomposite as Lithium-Ion Batteries Anodes

Lithium ion batteries (LIBs) are the enabling technology for many of the societal changes that are expected to happen in the following years. Among all the challenges for which LIBs are the key, vehicle electrification is one of the most crucial. Current battery materials cannot provide the required power densities for such applications and therefore, it makes necessary to develop new materials. Silicon is one of the proposed as next generation battery materials, but still there are challenges to overcome. Poor capacity retention is one of those drawbacks, and because it is tightly related with its high capacity, it is a problem rather difficult to address with common and scalable fabrication processes. Here we show that combining 0D and 1D silicon nanostructures, high capacity and stability can be achieved even using standard electrode fabrication processes. Capacities as high as 1200 mAh/g for more than 500 cycles at high current densities (2 A/g) were achieved with the produced hybrid 0D/1D electrodes. In this research, it was shown that while 0D nanostructures provide good strain relaxation capabilities, 1D nanomaterials contribute with enhanced cohesion and conductive matrix integrity.

0D↗
DOE OSTI2023

Structure–property relations in graphitic pebbles for nuclear applications

This work presents an analytical approach for holistically characterizing graphitic matrix pebbles for nuclear applications whereby the macrostructure, microstructure, and thermophysical properties of pebbles are determined. A systematic sectioning method was applied to several pebbles to describe the regional properties of the samples. Intact matrix-only spheres and sections of spheres fabricated by Kairos Power were characterized via optical imaging, x-ray computed tomography, x-ray diffraction, and ellipsometry to determine 2D and 3D macrostructure and anisotropy. The thermophysical properties of these materials were determined via measurements of density, specific heat, thermal expansion, and thermal diffusivity. The results of this study indicate that the pebble fabrication methods and their resultant effect on microstructure have a nontrivial effect on thermophysical properties, confirming the importance of robust characterization of these components. A discussion of the characterization approach and its applicability to nuclear fuel development activities is also included.

42 ENGINEERING↗
DOE OSTI2023

Progression of creep deformation from grain boundaries to grain interior in Al-Cu-Mn-Zr alloys

Creep mechanisms are studied in θ'-Al 2 Cu-strengthened Al-Cu-Mn-Zr alloys at 300 and 350°C for (i) ACMZ, a base alloy without further alloying elements and (ii) RR350, a commercial alloy with additions of Ni and Co forming distinct grain-boundary precipitates. At high stresses, creep is dominated by dislocations bypassing θ' precipitates within grains via the Orowan mechanism, as evidenced by (i) very high stress exponent (n~20-25) and (ii) α-Al and θ' lattice strains (measured via in-situ neutron diffraction) evolving during creep in a manner consistent with load transfer from the plastically-deforming α-Al matrix to elastically-deforming θ' precipitates. At intermediate stresses, both alloys exhibit a n~3 regime, where α-Al and θ' lattice strains scale near-linearly with applied stress while remaining largely unaffected by strain accumulation, indicating that Orowan looping or dislocation pile-up around θ' is now inactive within the grains. Rather, dislocation motion occurs solely in θ'-precipitate-free zones (θ'-PFZ) where high dislocation densities are observed via TEM after creep deformation. Plastic flow at θ'-PFZ and/or localized pipe diffusion are expected to enable grain-boundary sliding (GBS), which is proposed as the rate-limiting mechanism in the n~3 regime. Ni/Co-rich precipitates at RR350 grain-boundaries, with negligible θ'-PFZ around them, share load (as determined via neutron diffraction) with the α-Al matrix more effectively than θ-Al 2 Cu precipitates at ACMZ grain-boundaries, with wide surrounding θ'-PFZ. So, high creep resistance in the n~3 GBS regime of RR350 is enabled by coarsening-resistant grain-boundary precipitates, forming without concomitant development of weak θ'-PFZ, which effectively share load with the grains.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An adaptive Hessian approximated stochastic gradient MCMC method

Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. Here, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation (SA) to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the samples and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix via SA, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and demonstrates better learning results compared to standard SG-MCMC updating rules. The approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. Numerical experiments are performed on several problems, including sampling from 2D correlated distribution, synthetic regression problems, and learning the numerical solutions of heterogeneous elliptic PDE. The numerical results demonstrate great improvement in both the convergence rate and accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Compatibility of mesophase pitch and linear low-density polyethylene for low-cost carbon fiber

Mesophase pitch and linear low-density polyethylene (LLDPE) were blended by shear mixing and extruded producing precursor fibers which were then converted to carbon fibers (CFs). LLDPE additives agglomerate within the isotropic spheres rather than uniformly dispersing throughout the entire matrix. The blends possess different thermal properties and precursor tensile properties evident by two-phase tensile stress nature not present in precursor fibers produced from pure pitch. The resultant CFs were oxidized at 280 °C for 4 h and carbonized at 1000 °C for 15 min. This process created pitch CFs with a tensile strength of 1.23 GPa and a Young’s modulus of 99 GPa. All CFs produced from blends of pitch and LLDPE possessed a tensile strength of 0.63 GPa and less. These CFs showed a microstructure slightly different from the pure pitch CFs with fractured CF cross sections revealing voids present in all blended CFs due to burnt off LLDPE stemming from carbonization. Further study into compatibilizers and new conversion processes is required to successfully blend LLDPE with mesophase pitch for the purpose of producing a new low-cost CF.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A High-Throughput Method to Define Additive Manufacturing Process Parameters: Application to Haynes 282

This paper demonstrates how an analytical and experimental method can be used to rapidly define the 45 additive manufacturing settings for a new alloy where the process parameters were previously unknown. A nickel-based superalloy, Haynes 282, was chosen for the analysis. An experimental matrix of focused processing parameters was predicted with a dimensionless number and 100 samples were printed using the Selective Laser Melting technique. High-throughput measurements validated the predicted process conditions needed to achieve desired density and hardness. The whole process was completed in 16 hours. The new technique was confirmed with analytical processing maps adopted by the metal additive manufacturing community. With the predicted set of process parameters, microstructural characterizations and tensile testing were used to confirm the predictions. The resultant as-fabricated microstructures have refined length scales of both microsegregation and secondary phase distributions. As a result, the mechanical properties were comparable within the predicted processing window and exhibited high strength and high ductility.

Laser Material Processing↗
DOE OSTI2021

Development and characteristics of ultra high-performance lightweight cementitious composites (UHP-LCCs)

Highlights: • An ultra high-performance lightweight cementitious composite (UHP-LCC) was designed. • The UHP-LCCs had ultra high strength (>120 MPa) and very low density (3 ). • The UHP-LCCs had high structural efficiency, good functional properties and superior durability. • The strategies and mechanisms on achieving UHP-LCCs with excellent performance were proposed. High strength and light weight are two recent opposite development trends of concrete. This study proposed a design concept of an ultra high-performance lightweight cementitious composite (UHP-LCC), which had a compressive strength of higher than 120 MPa and an air-dried density down to around 1800 kg/m{sup 3}. The UHP-LCCs were innovatively developed by incorporating micro-sized hollow particles with a high strength shell (hollow glass microspheres, HGM) into an ultra-high performance cementitious composite (UHPC). The roles of HGM in the UHP-LCCs were investigated by evaluating the reactivity of the HGM and the mechanisms on achieving the excellent mechanical properties, low density and superior durability were revealed. The Chapelle test results showed that the HGM exhibited some pozzolanic reactivity, which facilitated the reaction between the shell of HGM and the alkali hydration products of the paste matrix. This chemical reaction was conducive to improving the HGM-paste interface and enhancing the mechanical properties. With the use of microspheres with a high stiff shell, the fundamental properties of the UHP-LCCs including thermal insulation, sound absorption, resistance to water ingress and electrical resistivity were improved significantly. The strategies for preparing the UHP-LCCs with high structural efficiency and great performance were proposed. The results of this study provide a new approach for designing and producing a lightweight UHPC, which would be a promising material for long-span structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Physical-chemical-mechanical quantitative assessment of the microstructural evolution in Portland-limestone cement pastes exposed to magnesium sulfate attack at low temperature

The changes in structural integrity and microstructure of Portland-limestone cement pastes were investigated in the course of magnesium sulfate attack at low temperature. A deterioration front, consisting of three distinct layers (brucite, gypsum, leached cement matrix), swelled in time due to the expansive nature of the deterioration products, generating cracks and subsequently detaching from the sound cement matrix, continuously promoting the process. Gypsum and thaumasite characterized the leached matrix, which experienced extensive cross-linking of the aluminosilicate structures, as a result of decalcification and dealumination of the calcium silicate hydrates (C(A)SH), impairing the overall mechanical performance. CSH of low packing density was most severely affected by the process, as confirmed by the significant drop in nano-mechanical properties. The increased rate of deterioration with limestone content was tentatively attributed to the prevalent morphology of the CSH phase. Results were validated by thermodynamic simulations, indicating that the real systems did not reach equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

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

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗
DOE OSTI2023

A neutron tomography study to visualize fluoride salt (FLiNaK) intrusion in nuclear-grade graphite

Manufactured graphite is a preferred material for in-core components of molten salt reactors and fluoride salt-cooled high-temperature reactors, which are in permanent contact with liquid salts. However, owing to the porous nature of nuclear graphite, under certain conditions, molten salts may intrude graphite's pores and affect graphite's properties and functionality. Therefore, a better understanding of molten salt intrusion (distribution across sample cross section and penetration depth) is needed to assess its effects. Here, in this work, we have demonstrated the use of neutron imaging (computed tomography) in the evaluation of salt penetration and distribution of a wide range of graphite grades with diverse microstructures that have been subjected to FLiNaK (LiF–NaF–KF) intrusion at 750 °C and 5 bar pressure for 12 h. Because of the great neutron attenuation contrast from scattering and adsorption between Li (from FLiNaK) and the graphite matrix, we have obtained direct visualization of FLiNaK salt distribution in the salt-impregnated graphites for the first time. Three-dimensional reconstructed images and cross-sectional concentration profiles demonstrate that salt penetration and density distribution are greatly dependent on the microstructural properties of the graphite grade.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗