SEARCH · Engineering Papers
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.
Research projects on defect structure, mechanical behavior, electrical and optical properties, and thermodynamics of solid materials at university
Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.
When kinetic conditions dominate and dictate the growth in a nanoparticle synthesis, properties of the synthesis environment can have considerable effects on the properties of the products. Such effects were studied here, where the solution environment was changed via the addition of KF to the hydrothermal syntheses of KTa 1– x Nb x O 3 and KTaO 3 . One result demonstrated the straightforward cause-and-effect relationship between the solution and reaction kinetics: KF directly increased the solution stability of the Ta species and therefore decreased its reaction rate, resulting in a change in composition heterogeneity of Ta and Nb in the KTa 1– x Nb x O 3 particles. However, not all effects are so simple; changing the chemical potential of the solution with KF can also promote the formation of particles with anisotropic defect enhanced kinetic Wulff shapes instead of cuboidal shapes. The increased F chemical potential in the solution enabled the formation of planar defects in the bulk, which accelerated growth in-plane to form particles characterized by flat rectangular flake geometries. In this work, thermodynamic modeling with density functional theory calculations confirmed that sufficient KF concentrations can drive the formation of a defect phase K n +1 Ta n O 3 n F.
Thermoelectric materials enable the energy conversion of waste heat into electricity, helpful to relieve global energy crisis. Here, we report a systematic investigation on high-temperature thermoelectric gadolinium selenides, cubic Gd 3-x Se 4 (x = 0.16, 0.21 and 0.25) and orthorhombic Gd 2 Se 3-y (y = 0.02, 0.06 and 0.08). High energy synchrotron x-ray diffraction and total scattering have been used to investigate the crystallographic and local structures. Atomic-scale clusters of Gd vacancy in the cubic phase are observed by employing the reverse Monte Carlo simulation. For cubic Gd 3-x Se 4 , adjusting Gd vacancy triggers the effect of multiple conduction bands, confirmed by the increase in effective masses and theoretical calculations. A reasonable peak zT value of 0.27 is achieved at 850 K for Gd 3-x Se 4 (x = 0.16). On the other hand, tuning Se vacancy enables the optimization of electron concentration for the orthorhombic Gd 2 Se 3-y . More significantly, its low deformation potential (Ξ = 12eV) gives rise to enhanced electron mobility and higher peak zT value of 0.54 at 850 K for Gd 2 Se 3-y (y = 0.02). Intriguingly, a higher zT of 1.2 at 1200 K is reasonably predicted by quality factor analysis. Finally, this work extends the scope of high-temperature thermoelectric materials and facilitates the exploration of novel high-temperature thermoelectric materials
Oxyhydrides are an intriguing class of materials in which there is partial replacement of the oxide ion with hydride ions and oxygen vacancies. Conventional synthesis relies on vacuum sealed ampules, using long reaction times at high temperatures, limiting accessibility. Here, we demonstrate a rapid, ambient-pressure route to oxyhydride formation using NaBH 4 under flowing argon in just 1 h. This approach significantly lowers experimental barriers, enabling broader exploration of these materials. The maximum hydride incorporation, obtained using a reaction temperature of 400 °C, is given by the formula SrTiO 2.945 H 0.049 , where the hydride, oxygen vacancy, and unpaired electron concentrations are determined through thermogravimetric analysis, quantitative solid-state nuclear magnetic resonance (NMR) spectroscopy, and electron paramagnetic resonance (EPR) spectroscopy. Density functional theory simulations of the 1 H NMR shifts for candidate point defects support the assignment of the observed hydride peak, validating the efficacy of the synthetic approach. A combination of in situ and ex situ studies of the reaction pathway reveal that hydride incorporation into the perovskite occurs via direct solid-solid reaction, with higher reaction temperatures favoring NaBH 4 decomposition and H 2 (g) release over oxyhydride formation. The electronic defect structure established from the EPR and NMR studies indicate that, at ambient temperature, anion vacant sites are occupied by single electrons, whereas hydride sites do not trap electrons. As a result, this work establishes a scalable synthesis strategy and provides a computational-experimental framework for understanding defect chemistry in oxyhydrides, opening pathways for their integration into energy and electronic applications.
Proposed viscoplastic material models include in their formulation observed material response but do not generally incorporate principles from thermodynamics, statistical mechanics, and quantum mechanics. Numerous hypotheses were made for material response based on first principles. Many of these hypotheses were tested experimentally. The proposed viscoplastic theories and the experimental basis of these hypotheses must be checked against the hypotheses. The physics of thermodynamics, statistical mechanics and quantum mechanics, and the effects of defects, are reviewed for their application to the development of constitutive laws.
The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.
In conventional processing, metals go through multiple manufacturing steps including casting, plastic deformation, and heat treatment to achieve the desired property. In additive manufacturing (AM) the same target must be reached in one fabrication process, involving solidification and cyclic remelting. The thermodynamic and kinetic differences between the solid and liquid phases lead to constitutional undercooling, local variations in the solidification interval, and unexpected precipitation of secondary phases. These features may cause many undesired defects, one of which is the so-called hot cracking. The response of the thermodynamic and kinetic nature of these phenomena to high cooling rates provides access to the knowledge-based and tailored design of alloys for AM. Here, we illustrate such an approach by solving the hot cracking problem, using the commercially important IN738LC superalloy as a model material. The same approach could also be applied to adapt other hot-cracking susceptible alloy systems for AM.
The defect fluorite yttrium niobate Y 3 NbO 7 and pyrochlore yttrium titanate Y 2 Ti 2 O 7 solid solutions have been synthesized via a solid state synthesis route. The resulting stoichiometry of the oxides is Y 2+x Ti 2-2x Nb x O 7 , where x = 0 to x = 1. All of the samples were single-phase; however, for those with a predominant fluorite phase, a small amount of additional pyrochlore phase was detected. The volume of the solid solution unit cells linearly increases with increase in yttrium niobate content. The water uptake increases with (x) and the protonic defect concentration reaches almost 4.5 × 10 -3 mol mol -1 at 300 °C. The calculated enthalpy of formation from oxides suggests strong stability for all of the compositions, with the values of enthalpy ranging from -84.6 to -114.3 kJ mol -1 . The total conductivity does not have a visible dependence on Y 3 NbO 7 content. For each compound, the total conductivity is higher in wet air. Interestingly, for samples where x < 0.5, the ratio of conductivity in hydrogen to air increases with increasing temperature, while for x > 0.5, the trend is the opposite.
Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.
Surface segregation is a phenomenon that depends on the delicate interplay between thermodynamic driving forces and kinetic obstacles, for which elevated temperature is often needed to enhance the atom mobility and reach equilibrium. Using the classic system of Cu 3 Au ( 100 ) under nonisothermal conditions, in this study, we show an adatom process underlying transient surface segregation dynamics through the temperature-change-driven creation and annihilation of thermal vacancies in the bulk and the resulting bulk-surface mass exchanges. This is demonstrated by monitoring the surface composition evolution of Cu 3 Au ( 100 ) with temperature changes between 250 and 500 ° C , showing that the increase in temperature decreases monotonically the surface Au concentration as a result of the transfer of more Cu than Au from the bulk to the surface to form Cu-rich clusters of adatoms. Such a bulk thermal defect effect is expected to be universal in inducing the disparity in the bulk-surface mass exchanges of dissimilar atoms in multicomponent materials because of the inherent differences in the vacancy formation energies of the constituent atoms.
Undoped mercury cadmium telluride crystals were subjected to high temperature equilibration at temperatures ranging from 400 C to 655 C in various Hg atmospheres. Hall effect and mobility measurements were carried out on the crystals quenched to room temperature subsequent to the high temperature equilibration. The variation of the hole concentration in the cooled crystals at 77 K as a function of the partial pressure of Hg at the equlibration temperatures, together with a comparison of the hole mobility in the undoped samples with that in the copper and phosphorous doped samples yielded a defect model for the undoped crystals, according to which, the undoped crystals are essentially intrinsic at the equilibration temperatures and the native acceptor defects are doubly ionized. Native donor defects appear to be negligible in concentration, implying that the p to n conversion in these alloys is mainly due to residual foreign donor impurities. The thermodynamic constants for the intrinsic excitation process as well as for the incorporation of the doubly ionized native acceptor defects in the undoped crystals were obtained.
Multicomponent alloys undergoing irradiation damage develop radiation-induced composition redistribution at point defect sinks such as grain boundaries (GBs) and dislocations. Such redistribution results in undesired changes to their mechanical behavior and corrosion resistance. Additively manufactured alloys proposed for future nuclear applications are expected to demonstrate a distinct response to irradiation owing to their unique microstructure with as-solidified dislocation density and chemical microsegregation. To capture the composition redistribution in such systems, we develop a mesoscale model with coupled evolution of atomic and point defect components in the presence of dislocation density, dislocation heterogeneity, and thermodynamic interactions at the GB. The model is parameterized for an FCC Fe–Cr–Ni alloy as a representative system for austenitic stainless steels, and simulations are performed in 1D and 2D as a function of irradiation temperature, dose, dislocation density, and grain size. Radiation-induced segregation (RIS) characterized by Cr depletion and Ni enrichment is predicted at both the GB and the dislocation cell wall, with RIS being lower in magnitude but wider at the cell wall. Strongly biased absorption of self-interstitials by dislocations is found to suppress Ni enrichment but slightly enhance Cr depletion under certain conditions. Thermodynamic segregation at the GB is predicted to be narrower and opposite in sign to RIS for both Cr and Ni. Importantly, non-monotonic segregation is found to occur when both thermodynamic and RIS mechanisms are considered, providing a novel physical interpretation of experimental observations. The model is expected to serve as a key tool in accelerated qualification of irradiated materials.
Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.
The past decade has seen a significant increase in research efforts aimed at understanding the thermodynamics of low-dimensional phases existing in many materials systems, ranging from two-dimensional materials to core regions of extended defects in crystalline solids. We review the current status of theoretical, computational, and experimental research on the “defect phases,” focusing on grain boundaries (GBs) in elemental and multicomponent polycrystalline materials. After reviewing the generalized concept of a phase of any dimensionality, we discuss recent progress in atomistic computer simulations of GB phase transformations and phase coexistences, including the observation of one-dimensional defects separating GB phases (defects in defects). Computational predictions compare well with experimental observations of multiple GB phases and segregation-induced phase transformations. An intriguing open question of GB thermodynamics is whether the GB free energy can be driven to a zero value by increasing solute segregation. We review recent efforts to understand this ultimate thermodynamic stabilization of GB phases and the possible polycrystalline microstructures that may arise. An outlook for future research in the field is discussed.
This proposed study aims at elevating Idaho National Laboratory (INL)’s portfolio in nuclear material research on two major fronts: 1) expanding INL’s material characterization capabilities to provide holistic microstructural and microchemical characterization ranging from atomic scale to millimeter scale, and 2) filling the knowledge gap on the effect of microstructural evolution on material phase stability at extreme conditions resembling those inside nuclear reactors, using Yttrium-Titanium-Oxide (YTO) as a model material. Previous research indicated the kinetics of phase transition of YTO materials varies between bulk samples and powder samples, suggesting that microstructural characteristics such as defect, micro-strain, and crystallite size could play a significant role in phase transition. The effect of crystallite size on the thermodynamics of phase transition has been well documented for other ceramic materials at nanometer crystallite size, with potential impacts from surface and strain energy and defects. In this proposed study, we will perform microstructural (micro-strain, crystallite size, defects, phase) and microchemical (elemental mapping) analyses through in-situ high temperature measurements using X-ray diffraction (XRD; at millimeter scale for powder and bulk samples) and transmission electron microscopy (TEM, at nanometer scale and atomic scale for powder and thin foil samples), to provide a comprehensive view of the effect of microstructural and microchemical evolution on phase transition and phase stability of YTO materials. These findings will also establish a roadmap and provide new insights on the phase stability of other materials that are critical to nuclear energy and other renewable energy infrastructure. Besides advancing our understanding of this fundamental material science question, this project will also enable the development of state-of-the-art material characterization capabilities at INL/MFC for microstructural and microchemical analysis at high temperature conditions. The new collaboration between the High-Resolution Materials Characterization group at INL and the Electron Microscopy group at the University of Florida (UF) in this project will open the possibilities for more fruitful and synergetic cutting-edge nuclear material research in the future.
Elucidating the factors limiting quantum coherence in real materials is essential to the development of quantum technologies. Here we report a strategic approach to determine the effect of lattice dynamics on spin coherence lifetimes using oxygen deficient double perovskites as host materials. In addition to obtaining millisecond T 1 spin-lattice lifetimes at T ~ 10 K, measurable quantum superpositions were observed up to room temperature. We determine that T 2 enhancement in Sr 2 CaWO 6-δ over previously studied Ba 2 CaWO 6-δ is caused by a dynamically-driven increase in effective site symmetry around the dominant paramagnetic site, assigned as W 5+ via electron paramagnetic resonance spectroscopy. Further, a combination of experimental and computational techniques enabled quantification of the relative strength of spin-phonon coupling of each phonon mode. This analysis demonstrates the effect of thermodynamics and site symmetry on the spin lifetimes of W 5+ paramagnetic defects, an important step in the process of reducing decoherence to produce longer-lived qubits.
Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.