Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Machine Learning Approach for the Prediction of Formability and Thermodynamic Stability of Single and Double Perovskite Oxides

Perovskite oxides continue to attract huge interest due to their fascinating and wide-ranging properties for diverse applications. The tunability of these properties may be further enhanced by increasing their compositional complexity via double perovskite-ordered configurations containing multiple cations. In this work, we focus on an exhaustive chemical space of single and double oxide perovskites and optimally explore this space to identify novel compositions that are likely to form stable compounds. Critically, we examine the relationship between formability, the practical ability to synthesize a compound, and stability, the thermodynamic preference to form the structure. Our formability and stability training data sets were enumerated from the available experimental literature and in-house density functional theory computations and contained 1505 and 3469 examples, respectively, representing state-of-the-art in the current open literature in perovskite and double perovskite compounds. Subsequently, cross-validated and highly accurate machine learning classification models are built using these training data sets and employed to screen for novel stable oxide perovskites. The study identifies (1) atomic features relevant to prediction of formability and stability in perovskite and double perovskite compounds, (2) the importance of including energy contributions due to local structural relaxations going beyond the high symmetry perovskite phase, and (3) 437,828 double perovskite compounds that are likely to be stable and 891,188 compounds that are likely to be formable. From the intersection of this large chemical space of formable and stable oxide perovskites, 414 compositions are identified as the most promising candidates for future experimental synthesis of novel oxide perovskites. The developed models may be generalized and have implications beyond perovskite discovery if applied to other families of compounds.

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

Unified Approach to Polarons and Phonon-Induced Band Structure Renormalization

Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green’s function theory of band structure renormalization that incorporates localization and self-trapping. Here, we show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Soft fragmentation on the celestial sphere

We develop two approaches to the problem of soft fragmentation of hadrons in a gauge theory for high energy processes. The first approach directly adapts the standard resummation of the parton distribution function’s anomalous dimension (that of twist-two local operators) in the forward scattering regime, using k T -factorization and BFKL theory, to the case of the fragmentation function by exploiting the mapping between the dynamics of eikonal lines on transverse-plane to the celestial-sphere. Critically, to correctly resum the anomalous dimension of the fragmentation function under this mapping, one must pay careful attention to the role of regularization, despite the manifest collinear or infra- red finiteness of the BFKL equation. The anomalous dependence on energy in the celestial case, arising due to the mismatch of dimensionality between positions and angles, drives the differences between the space-like and time-like anomalous dimension of parton densities, even in a conformal theory. The second approach adapts an angular-ordered evolution equation, but working in 4 – 2ϵ dimensions at all angles. The two approaches are united by demanding that the anomalous dimension in 4 – 2ϵ dimensions for the parton distribution function determines the kernel for the angular-ordered evolution to all orders.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. Here, we finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

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

Fundamental Studies of the Vibrational, Electronic, and Photophysical Properties of Tetrapyrrolic Architectures

The ability to capture and utilize light in the near-ultraviolet (NUV), visible and near-infrared (NIR-I and NIR-II) spectral regions (i.e., 320–400, 400–700, 700–1000, 1000–1700 nm) is essential for any solar-energy conversion scheme. Nature employs chlorophylls and bacteriochlorophylls in light-harvesting architectures to absorb light in the blue and red/NIR regions. Accessory pigments (carotenoids, bilins) augment absorption of the (bacterio)chlorophylls in the green region. The harvested energy is funneled to a reaction center protein, where charge separation occurs. Subsequent migration of the electron and the hole stabilizes and stores the energy from light via redox chemistry. The long-term objective of the Bocian/Holten&Kirmaier/Lindsey research program under this DOE grant has been to develop tetrapyrrole-based molecular architectures that absorb sunlight, funnel energy and separate charge with high efficiency. Integral to the program has been iterative cycles of design, synthesis and characterization that provided deep insights into the relationships between chemical composition, electronic structure, and key static and dynamic properties (vibrational, redox, photophysical, energy/charge transfer) of tetrapyrrolic systems. Such architectures included monomers, dyads, larger arrays, and complexes with accessory components. The objective was to develop molecular designs and guiding principles to enhance current and future energy-conversion schemes. Molecular arrays targeted to address one or more fundamental questions concerning light harvesting and energy/charge transfer were constructed from analogues of the naturally occurring hemes, chlorophylls and bacteriochlorophylls. Diverse, tunable synthetic building blocks were prepared that spanned the three respective tetrapyrrole families, which are the porphyrins, chlorins and bacteriochlorins. Thus, the research focused on porphyrins as well as synthetic surrogates for chlorophylls (chlorins, 13 1 -oxophorbines and chlorin-imides) and bacteriochlorophylls (bacteriochlorins, bacterio-13 1 -oxophorbines and bacteriochlorin-imides), generically termed hydroporphyrins. Although the three tetrapyrrole classes (porphyrins, chlorins and bacteriochlorins) absorb light strongly in the violet-blue spectral region, the long-wavelength absorption band typically lies in the green-orange, red, and NIR regions, respectively, with increasing intensity. Understanding the spectra, electronic structure, and energy/charge-transfer properties of such tetrapyrrolic macrocycles is of central importance for the rational design of molecular architectures for solar-energy conversion. Our integrated program of molecular design and synthesis coupled with a variety of spectroscopic, electrochemical, and computational studies have probed from first principles how structural and electronic properties of tetrapyrrolic macrocycles dictate spectral properties as well as the rates of ground-state hole/electron transfer and excited-state energy flow in multicomponent architectures. Individual molecules and multicomponent architectures were designed to test ideas of fundamental importance, often requiring the development of new synthetic methodology. The members of the collaborative team had almost daily discussions by phone and/or e-mail concerning design of molecules, flow of compounds between the labs, planning of physical characterization studies, discussing results and analysis and integrating into design of next generation architectures, and the preparation of manuscripts. Furthermore, students and postdocs in the different labs routinely communicated with one another to facilitate the advancement of the research activities. In short, a highly integrated and collaborative research program was well established among the groups. The research effort involved molecular design and synthesis of synthetic molecular architectures by the Lindsey group integrated with physicochemical and photophysical characterization by the Bocian group and the Holten&Kirmaier group (Figure 2). The Bocian group carried out electrochemical, electron paramagnetic resonance (EPR), resonance Raman (RR), and Fourier-transform infrared (FT-IR) studies, as well as density functional theory (DFT) calculations and the time-dependent extension (TDDFT) to gain insight into excited-state properties. The Holten&Kirmaier group carried out static and time-resolved absorption and fluorescence spectroscopy studies and simulated absorption spectra using molecular orbital (MO) energies from DFT as input to the four-orbital model to complement the TDDFT calculations. The combined measurements provided understanding of the vibrational/electronic properties of the individual molecules and the changes that occur upon incorporation into multicomponent architectures. This information underpinned elucidating the mechanisms and timescales of ground-state hole/electron transfer and excited-state energy and charge transfer.

14 SOLAR ENERGY↗
NASA NTRS1983

Numerical models for cometary dust environments

The development of a numerical model for the dust environment near a comet is described. The major components of the model are a simple treatment of the acceleration of the dust by the gas outflow from an isothermal icy nucleus, a detailed treatment of the dust particle trajectories controlled by solar gravity and radiation pressure, including the formation of the dust envelope and tail, and a simple treatment of the radiation from the dust in the visible and thermal infrared. The results of recent comet observations are used to estimate parameter values for the model, including dust-to-gas ratios by light and by mass, the distribution of dust and gas emission on the nuclear surface, the distributions of dust particle size, density, mass, and radiative properties, and other quantities. Gas and dust production rates as functions of time are derived from observed light curves using the simple photometric theory. Sample results are presented for the dust concentration, flux and fluence along the trajectory of a spacecraft which flies by Comet P/Halley in 1986.

Divine, N.↗
NASA NTRS1992

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗
DOE OSTI2024

Zirconium‐Based Metal–Organic Frameworks with Free Hydroxy Groups for Enhanced Perfluorooctanoic Acid Uptake in Water

Perfluorooctanoic acid (PFOA) is a highly recalcitrant organic pollutant, and its bioaccumulation severely endangers human health. While various methods are developed for PFOA removal, the targeted design of adsorbents with high efficiency and reusability remains largely unexplored. Here, in this work, the rational design and synthesis of two novel zirconium‐based metal‒organic frameworks (MOFs) bearing free ortho ‐hydroxy sites, namely noninterpenetrated PCN‐1001 and twofold interpenetrated PCN‐1002, are presented. Single crystal analysis of the pure ligand reveals that intramolecular hydrogen bonding plays a pivotal role in directing the formation of MOFs with free hydroxy groups. Furthermore, the transformation from PCN‐1001 to PCN‐1002 is realized. Compared to PCN‐1001, PCN‐1002 displays higher chemical stability due to interpenetration, thereby demonstrating an exceptional PFOA adsorption capacity of up to 632 mg g −1 (1.53 mmol g −1 ), which is comparable to the reported record values. Moreover, PCN‐1002 shows rapid kinetics, high selectivity, and long‐life cycles in PFOA removal tests. Solid‐state nuclear magnetic resonance results and density functional theory calculations reveal that multiple hydrogen bonds between the free ortho ‐hydroxy sites and PFOA, along with Lewis acid‐base interaction, work collaboratively to enhance PFOA adsorption.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Single-Point Extrapolation to the Complete Basis Set Limit through Deep Learning

Machine learning (ML) offers an attractive method for making predictions about molecular systems while circumventing the need to run expensive electronic structure calculations. Once trained on ab initio data, the promise of ML is to deliver accurate predictions of molecular properties that were previously computationally infeasible. In this work, we develop and train a graph neural network model to correct the basis set incompleteness error (BSIE) between a small and large basis set at the RHF and B3LYP levels of theory. Our results show that, when compared to fitting to the total potential, an ML model fitted to correct the BSIE is better at generalizing to systems not seen during training. We test this ability by training on single molecules while evaluating on molecular complexes. We also show that ensemble models yield better behaved potentials in situations where the training data is insufficient. However, even when only fitting to the BSIE, acceptable performance is only achieved when the training data sufficiently resemble the systems one wants to make predictions on. As a result, the test error of the final model trained to predict the difference between the cc-pVDZ and cc-pV5Z potential is 0.184 kcal/mol for the B3LYP density functional, and the ensemble model accurately reproduces the large basis set interaction energy curves on the S66x8 dataset.

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

Unraveling the Mechanism of Interfacial Charge Transfer and Photoresponsivity of WS 2 Quantum Dots/MoS 2 (0D–2D) Heterostructure-Based Transistors

To assess the efficacy of a mixed-dimensional van der Waals (vdW) heterostructure in modulating the optoelectronic responses of nanodevices, the charge transport properties of the transition-metal dichalcogenide (TMD)-based heterostructure comprising zero-dimensional (0D) WS 2 quantum dots (QDs) and two-dimensional (2D) MoS 2 flakes are critically analyzed. Herein, a facile strategy was materialized in developing an atomically thin phototransistor assembled from mechanically exfoliated MoS 2 and WS 2 QDs synthesized using a one-pot hydrothermal route. The amalgamated photodetectors exhibited a high responsivity of ∼8000 A/W at an incident power of 0.05 nW of white light, surpassing that of the pristine MoS 2 devices. Furthermore, the detectivity of pristine MoS 2 , which was on the order of 10 10 , increased to 10 12 Jones for the WS 2 QDs/MoS 2 heterostructure photodetector, outperforming other WS 2 -based materials. The quasiparticle band gap and density of states (DOS) are further analyzed to elucidate the photophysics of the WS 2 QD/MoS 2 hybrid assembly. The difference in the work function between MoS 2 and WS 2 QDs gives rise to an electric field across the 0D–2D interface, facilitating effective charge separation and migration and contributing to the enhancement of photoresponsivity. The analysis of optical responses using density functional theory (DFT) revealed stronger absorption and less reflection over a broader spectrum of wavelengths for the heterostructure compared to the pristine materials. The estimated optical conductivity aligns well with the experimentally predicted maximum photoresponsivity under visible light, which is attributed to the high absorbance of 2D MoS 2 . In conclusion, combining diverse spectroscopic and imaging techniques with quantum simulation provides insights that clarify the pertinence of 0D–2D TMDs in designing phototransistors.

WS2 QDs/MoS2 heterostructure↗
DOE OSTI2025

Regioisomeric Engineering for Multicharge and Spin Stabilization in Two-Electron Organic Catholytes

Developing multicharge and spin stabilization strategies is fundamental to enhancing the lifetime of functional organic materials, particularly for long-term energy storage in multiredox organic redox flow batteries. Current approaches are limited to the incorporation of electronic substituents to increase or decrease the overall electron density or bulky substituents to sterically shield reactive sites. With the aim to further expand the molecular toolbox for charge and spin stabilization, we introduce regioisomerism as a scaffold-diversifying design element that considers the collective and cumulative electronic and steric contributions from all of the substituents based on their relative regioisomeric arrangements. Through a systematic study of regioisomers of near-planar aromatic cyclic triindoles and nonplanar nonaromatic cyclic tetraindoles, we demonstrate that this regioisomeric engineering strategy significantly enhances the H-cell cycling stability in the above two new classes of 2e – catholytes, even when current strategies failed to stabilize the multicharged species. Density functional theory calculations reveal that the strategy operates by redistributing the charge and spin densities while highlighting the role of aromaticity in charge stabilization. The most stable 2e – catholyte candidate was paired with a viologen derivative anolyte to achieve a proof-of-concept all-organic flow battery with 1.26–1.49 V, 98% capacity retention, and only 0.0117% fade/h and 0.00563% fade/cycle over 400 cycles (192 h), which is the highest capacity retention ever reported over 400 cycles in a multielectron all-organic flow battery setup. We anticipate regioisomeric engineering to be a promising strategy complementary to conventional electronic and steric approaches for multicharge and spin stabilization in other functional organic materials.

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

DC Hall coefficient of the strongly correlated Hubbard model

The Hall coefficient is related to the effective carrier density and Fermi surface topology in non-interacting and weakly interacting systems. In strongly correlated systems, the relation between the Hall coefficient and single-particle properties is less clear. Clarifying this relation would give insight into the nature of transport in strongly correlated materials that lack well-formed quasiparticles. In this work, we investigate the DC Hall coefficient of the Hubbard model using determinant quantum Monte Carlo in conjunction with a recently developed expansion of magneto-transport coefficients in terms of thermodynamic susceptibilities. At leading order in the expansion, we observe a change of sign in the Hall coefficient as a function of temperature and interaction strength, which we relate to a change in the topology of the apparent Fermi surface. We also combine our Hall coefficient results with optical conductivity values to evaluate the Hall angle, as well as effective mobility and effective mass based on Drude theory of metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

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

Evaluation of Gamma Ray Production - Quarterly Report FY22Q2

The LLNL team has performed the first systematic study of the effect of two-body weak currents in β decay rates across the entire mass table. In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. We examined the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. We extended the charge-changing finiteamplitude pnFAM code to incorporate the contributions of two-body currents to the usual onebody linear response in the Gamow-Teller channel. We found that the two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows; see Fig. 1. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. This new physics capability is unique and provides the highest-fidelity description of β decay on the market. The underlying computational implementation can be reused to calculate γ-strength functions for finite-range nuclear potentials such as the Gogny force, which will allow us to benchmark some of our future developments against published results. The paper is in print in Phys. Rev. C.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multi-modal characterization of the B2 phase in the Ta-Re binary system

The energy and transportation industries demand materials that retain their mechanical property at high temperatures. Refractory complex concentrated alloys (RCCAs) with a BCC + B2 microstructure offer a potential solution, where maintaining the high temperature mechanical properties can be achieved by precipitation strengthening. This depends on the B2 phase in RCCAs being thermodynamically stable with a high solvus temperature. Recently, we predicted the high temperature stability of the B2 structure in the Ta-Re binary system, using density functional theory. Here, in this study, we provide experimental evidence for the existence of this phase for the first time, using a Ta 65 Re 35 alloy. Despite Ta-Re binary phase diagrams predicting a single-phase BCC microstructure for Ta 65 Re 35 , we show that a high Z nanoscale secondary phase appears after heat treatment at 1550 °C and 1100 °C. Scanning transmission electron microscopy (STEM) revealed that this phase has a cubic structure and is equiatomic TaRe though B2 superlattice reflections were absent in fast Fourier transforms (FFT) and diffraction patterns (DPs). DP simulations indicate that the B2 TaRe superlattice reflections are up to two orders of magnitude weaker than their fundamental reflections, making their detection challenging via electron microscopy. Neutron diffraction confirmed the second phase had a B2 structure. This study identified a previously unobserved high temperature stable B2 phase in the Ta-Re system, enabling the development of new high temperature BCC + B2 RCCAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗