Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HoPS”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Defect Thermodynamics and Transport Properties of Proton Conducting Oxide BaZr 1–x Y x O 3–δ (x ≤ 0.1) Guided by Density Functional Theory Modeling

Density functional theory-based thermodynamic modeling was performed to determine the effect of humidity and H 2 /O 2 gas pressure on the defect chemistry and transport properties of proton conducting oxide BaZr 1-x Y x O 3-δ (x ≤ 0.1) in the temperature range of 800–1200K, relevant for solid oxide fuel/electrolysis cell applications. The first-principles charge defect analysis was carried out to obtain the defect energetic ensembles as a function of Fermi level as well as the dependence of the hydration and oxidation reaction energies on complex defect-dopant configurations. It is shown that oxygen vacancies introduced to compensate for the extra Ba vacancies may cause a slight increase in the concentration of hydroxyl and proton species upon hydration of Ba deficient BaZr 1-x Y x O 3-δ , while the overall proton diffusivity is predicted to decrease with increased proton hopping barriers upon trapping near the Ba vacancies. The developed defect model is further demonstrated to be able to describe the bulk defect chemistry and transport properties of BaZr 0.9 Y 0.1 O 3-δ under the solid oxide cell operating conditions. Finally, the theoretical framework developed in this work further allows the inclusion of configurational, energetic and electronic characteristics of the defects that can be used as a complementary tool to experimental measurements for testing various mechanistic pathways or to provide essential mechanistic data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The hierarchical energy landscape of screw dislocation motion in refractory high-entropy alloys

High-entropy alloys (HEAs) with concentrated solid solutions are conceived to possess a rugged atomic and energy landscape in which dislocation motion necessarily proceeds to accommodate mechanical deformation. Fundamental questions remain as to how rough the energy landscape is and to what extent it can be influenced by the local ordering of their constituent elements. Here, we construct and report the potential energy landscape (PEL) governing screw dislocation motion in refractory HEAs that reveals a hierarchical and multilevel structure with a collection of small basins nested in large metabasin. This striking feature pertaining to HEAs exerts a trapping force and back stress on saddle point activations, retarding dislocation movement. Further, by introducing chemical short-range order, the energy landscape is smoothed but skewed to different degrees that shifts the rate-liming process from kink-glide to kink-pair nucleation. The chemical disorder-roughened PEL in HEAs, analogous to structural disorder induced in metallic glasses, signifies the role of various barrier-hopping processes underlying the extraordinary mechanical behaviors of the two distinct groups of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Barrier-free predictions of short-range ordering/clustering kinetics in binary FCC solid solutions

We present comparisons of kinetic Monte Carlo (kMC) simulations of isothermal short-range ordering (SRO) and clustering (SRC) kinetics in binary FCC alloys with a mean-field concentration wave (CW) model. We find that the CW model is able to give order-of-magnitude agreement with kMC simulations for ordering/clustering relaxation times over a wide range of temperatures and compositions. The advantage of the CW model is that it does not require parameterization of vacancy hopping energy barriers, which, for a concentrated alloy, becomes prohibitive. We assess limits in the accuracy of the model and discuss the effect of cooling rates as well as the extension to multi-component systems. Ultimately, the simplicity and performance of the CW model compared to kMC simulations suggests that it is a useful tool to connect with models of properties dependent on SRO/SRC as well as for designing thermal treatments to control formation of SRO/SRC.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spectral properties of three-dimensional Anderson model

The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level statistics. While the former studies the sensitivity of Hamiltonian eigenlevels upon inserting a magnetic flux, the latter studies the properties of unperturbed eigenlevels. We define two versions of dimensionless conductance, the first corresponding to the width of the level curvature distribution relative to the mean level spacing, and the second corresponding to the ratio of the Heisenberg time and the Thouless time obtained from the spectral form factor. We show that both conductances look remarkably similar around the localization transition, in particular, they predict a nearly identical critical point consistent with other well-established measures of the transition. We then study some further properties of those quantities: for level curvatures, we discuss particular similarities and differences between the width of the level curvature distribution and the characteristic energy studied by Edwards and Thouless (1972) in their pioneering work, in which the hopping at one lattice edge is changed from periodic to antiperiodic boundary conditions. In the context of the spectral form factor, we show that at the critical point it enters a broad time-independent regime, in which its value is consistent with the level compressibility obtained from the level variance. Finally, we test the scaling solution of the average level spacing ratio in the crossover regime using the cost function minimization approach introduced recently in Šuntajs et al. (2020). The latter approach seeks for the optimal scaling solution in the vicinity of the crossing point, while at the same time it allows for the crossing point to drift due to finite-size corrections. We find that the extracted transition point and the scaling coefficient agree with those from the literature to high numerical accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Strong-coupling expansion of multi-band interacting models: Mapping onto the transverse-field J 1 - J 2 Ising model

We investigate a class of two-dimensional two-band microscopic models in which the inter-band repulsive interactions play the dominant role. We first demonstrate three different schemes of constraining the ratios between the three types of inter-band interactions – density-density, spin exchange, and pair-hopping – that render the model free of the fermionic sign-problem for any filling and, consequently, amenable to efficient Quantum Monte Carlo simulations. We then study the behavior of these sign-problem-free models in the strong-coupling regime. In the cases where spin-rotational invariance is preserved or lowered to a planar symmetry, the strong-coupling ground state is a quantum paramagnet. However, in the case where there is only a residual Ising symmetry, the strong-coupling expansion maps onto the transverse-field J 1 - J 2 Ising model, whose pseudospins are associated with local inter-band magnetic order. We show that by varying the band structure parameters within a reasonable range of values, a variety of ground states and quantum critical points can be accessed in the strong-coupling regime, some of which are not realized in the weak-coupling regime. Here, we compare these results with the case of the single-band Hubbard model, where only intra-band repulsion is present, and whose strong-coupling behavior is captured by a simple Heisenberg model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Band structure dependent electronic localization in macroscopic films of single-chirality single-wall carbon nanotubes

Significant understanding has been achieved over the last few decades regarding chirality-dependent properties of single-wall carbon nanotubes (SWCNTs), primarily through single-tube studies. However, macroscopic manifestations of chirality dependence have been limited, especially in electronic transport, despite the fact that such distinct behaviors are needed for many applications of SWCNT-based devices. In addition, developing reliable transport theory is challenging since a description of localization phenomena in an assembly of nanoobjects requires precise knowledge of disorder on multiple spatial scales, particularly if the ensemble is heterogeneous. Here, we report an observation of pronounced chirality-dependent electronic localization in temperature and magnetic field dependent conductivity measurements on macroscopic films of single-chirality SWCNTs. The samples included large-gap semiconducting (6,5) and (10,3) films, narrow-gap semiconducting (7,4) and (8,5) films, and armchair metallic (6,6) films. Experimental data and theoretical calculations revealed Mott variable-range-hopping dominated transport in all samples, while localization lengths fall into three distinct categories depending on their band gaps. Armchair films have the largest localization length. Here our detailed analyses on electronic transport properties of single-chirality SWCNT films provide significant new insight into electronic transport in ensembles of nanoobjects, offering foundations for designing and deploying macroscopic SWCNT solid-state devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermodynamics and kinetics of Pb intercalation under graphene on SiC(0001)

SiC-supported graphene intercalated by a two-dimensional Pb monolayer can provide an appealing platform for spintronic applications. Such a monolayer structure is thermodynamically ultrastable, as observed in recent experiments. However, important fundamentals such as the structure of intercalated phases, locations of intercalated atoms, thermodynamic preference for intercalation, and intercalation pathways for this system have not yet been understood conclusively. Here, in this work, extensive density functional theory calculations are performed to assess Pb intercalation thermodynamics and kinetics under graphene on SiC(0001). We find that intercalation of isolated Pb atoms is strongly disfavored over adsorption on top of graphene. However, intercalation of complete Pb layers in the gallery between SiC and graphene buffer layer is strongly favored over supported Pb monolayers and moderately favored over formation of supported large three-dimensional Pb islands. We also find that initiation of intercalation either by individual Pb atoms directly penetrating graphene or by hopping under a static graphene step edge is energetically prohibitive at experimental temperatures. Consequently, more complex intercalation pathways are operative and further analyzed. We demonstrated that once an intercalated Pb monolayer forms around a graphene step edge, a facile Pb mass transport by Pb vacancy-mediated diffusion can be triggered for continued growth of the intercalated monolayer.

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

Revisiting point defect thermodynamics in group IVB and VB transition metal carbides

We present a comprehensive re-examination of point defect thermodynamics in group IVB and VB transition metal carbides (TMCs) with the rocksalt structure using a combination of density functional theory (DFT) calculations and a statistical mechanical Wagner-Schottky model within the canonical ensemble. The most stable configurations of point defects were discovered using basin-hopping global optimization, driven by either a machine learning interatomic potential (MLIP) or DFT. A key finding is the identification of previously unreported dicarbon antisites—a C–C dimer occupying a metal site—as the structural (constitutional) defects on the carbon-rich side of stoichiometry in all group IVB and VB TMCs except TaC. Furthermore, dicarbon antisite-containing thermal defect complexes, such as quadruple and interbranch defects, can dominate in TMCs under specific stoichiometric and temperature conditions. In conclusion, by incorporating dicarbon antisites into the defect landscape, this work provides a revised understanding of the thermodynamics of point defects in TMCs.

Carbides↗
DOE OSTI2021

A critical dislocation velocity for serration mechanism transition in a nickel-chromium solid solution alloy

Here, the influence of strain rate across three orders of magnitude (1.70 × 10 –5 /s to 1.43 × 10 –2 /s) along with the effect of the plastic strain accumulation (up to 10%) on the serrated plastic flow were investigated in the nickel-chromium (Ni-Cr) solid solution alloy Nimonic 75 by performing constant-strain-rate tension testing at 600 °C. As the strain rate decreased, the critical strain for the onset of serrations transitioned from normal behavior to inverse behavior. The serrated flow was characterized as Type A+B serration at high strain rate (1.43 × 10 –2 /s). In the intermediate strain-rate regime (1.43 × 10 –3 /s and 1.45 × 10 –4 /s), Type B serrations were observed and followed by a transformation to Type C+B serrations. At the low strain rate (1.70 × 10 –5 /s), the plastic flow immediately displayed Type C serrations, which later evolved into Type C+B serrations. Regardless of the strain rate, plastic strain, or dislocation density, a critical dislocation velocity falling in the range of 1.2 × 10 –6 – 2.2 × 10 –6 m/s was identified to signify the onset of Type C serration, whereby the mobile dislocations break free from the solute cloud for short bursts of deformation. Finally, a novel model by solute rearrangement across dislocation cores was used to understand how the critical dislocation velocity is quantitatively determined by the rate at which solute atoms are able to hop across the glide plane as a partial dislocation core moves through the lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thin film growth effects on electrical conductivity in entropy stabilized oxides

Entropy stabilization has garnered significant attention as a new approach to designing novel materials. Much of the work in this area has focused on bulk ceramic processing, leaving entropy-stabilized thin films relatively under-explored. Following an extensive multi-variable investigation of polycrystalline (Mg 0.2 Co 0.2 Ni 0.2 Cu 0.2 Zn 0.2 )O thin films deposited via pulsed laser deposition (PLD), it is shown here that substrate temperature and deposition pressure have strong and repeatable effects on film texture and lattice parameter. Further analysis shows that films deposited at lower temperatures and under lower oxygen chamber pressure are 40x less electrically resistive than otherwise identical films grown at higher temperature and pressure. Annealing these films in an oxygen-rich environment increases their electrical resistivity to match that of the films grown at higher temperatures and pressures. Because of this, the electric conductivity is hypothesized to be the result of polaron hopping mediated by transition metal valence changes which compensate for oxygen off-stoichiometry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of the iodate sorption mechanism by CoAl LDH through experiments and ab initio molecular dynamics simulations

Radioiodine released during the nuclear-fuel cycle constitutes a persistent radiological hazard. In this study, the IO 3 – uptake mechanism of CoAl LDH was resolved by combining pH-controlled sorption experiments, synchrotron XAFS, and DFT-based AIMD simulations. At pH close to 6, approximately 90% of IO 3 – was removed, and the equilibrium distribution coefficient reached about 1.7 × 10 4 mL g –1 . EXAFS analysis indicated an average iodine–oxygen bond length of 1.81 Å and a coordination number near 3, with the fit R-factor equal to 0.002. The simulations faithfully reproduced the experimental spectrum and revealed transient proton hopping events that generated metastable I–O–H species inside the interlayer, thereby confirming nitrate-to-iodate exchange as the controlling capture pathway. Atomic density profiles and radial distribution functions further showed that IO 3 – adopt an end-on orientation perpendicular to the hydroxide sheets, while water molecules mediate proton migration without disturbing the host lattice. In conclusion, the integrated experimental–computational evidence demonstrates that CoAl LDH can rapidly and selectively sequester IO 3 – under near-neutral conditions, offering atomic scale guidance for the rational engineering of layered sorbents for advanced radioactive-waste treatment.

Kang, Jaehyuk [Jeju National Univ. (Korea, Republi↗
DOE OSTI2020

Comparison of continuum and cross-core theories of dynamic strain aging

Dynamic strain aging (DSA) is the process of solute atoms segregating around dislocations on the timescale of loading. Continuum theories of DSA derived from elasticity theory have been shown to severely overpredict both the timescale and strengthening of DSA. Recently, cross-core theory was developed to reconcile this gap, invoking a special single-atomic-hop diffusion mechanism across the core of an extended dislocation. In this work, we show that the classical continuum theory expression for the rate of solute segregation is in error. After correcting this error, we show that continuum theory predictions match cross-core theory when the elevated diffusivity near the dislocation core is accounted for. Finally, our findings indicate that continuum theory is still a useful tool for studying dislocation-solute interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2020

Theoretical study of the electronic, thermodynamic, and thermo-conductive properties of y-LiAlO2 with 6Li isotope substitutions for tritium production

Lithium aluminate has attracted researchers’ interests due to its wide applications. The structural, electronic, optical, lattice phonon thermodynamic, and thermo-conductive properties of ?-LiAlO2 with 6Li-isotope substitution are investigated by density functional theory and lattice phonon dynamics. The calculated results show that ?-LiAlO2 possesses a wide band-gap of 4.63 eV. The valence band is mainly formed by p orbitals of Li, O, and Al. The calculated lattice thermal conductivity of ?-LiAlO2 is found to be in good agreement over a wide range of temperatures with the available experimental data. The calculated dielectric matrix shows anisotropic behavior along z and x (or y) axes. Due to the wide band-gap, the optical conductivity of ?-LiAlO2 is mainly contributed from the electron hopping between valence-conduction bands. Substitution of 7Li with 6Li isotope in ?-LiAlO2 leads to observable differences in the lattice phonon frequencies at higher frequencies (>8THz) and to slight changes in corresponding thermal conductivity and the infrared and Raman spectra which are in good agreement with the measured data. The relevance of these properties for tritium production through 6Li absorbing neutron in nuclear reactors is discussed.

ab initio thermodynamics, ?-LiAlO2 and 6Li-isotope↗
DOE OSTI2026

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗
DOE OSTI2021

Vacancy surface migration mechanisms in dilute nickel-chromium alloys

We investigate the unitary mechanisms related to the surface migration of vacancies in dilute Ni-Cr alloys via first-principle calculations. We survey a complete set of surface and sub-surface migration paths for vacancies near the (100) free surface and calculate the corresponding migration barriers. Our results show that a vacancy migrating towards the free surface will face lower energy barriers to migrate via atomic exchange with a neighboring Cr atom rather than with a Ni atom. Once a vacancy reaches the free surface, it will be trapped there. Our results also reveal that, when a Cr atom sits in the atomic plane just below the free surface, any in-plane vacancy hopping jump that would result in the vacancy sitting on top a subsurface Cr atom is energetically unfavorable. Taken together, these fundamental unitary surface migration mechanisms offer insights into the complex interactions between surface segregation and vacancy migration phenomena in Ni-Cr-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unraveling small-scale defects in irradiated ThO 2 using kinetic Monte Carlo simulations

Point defects and their clusters generated through irradiation can have significant impact on the physical and mechanical properties of materials. However, direct experimental visualization of these small-scale defects using high-resolution scanning transmission electron microscopy remains a challenging task. In this study, using thorium dioxide (ThO 2 ) with the fluorite structure as a model system, we demonstrate the use of ab initio basin-hopping simulations in synergy with object kinetic Monte Carlo simulations as a powerful tool for identifying small defect complexes in irradiated materials. In addition to providing quantitative insights into defect evolution in ThO 2 under irradiation, our study reveals an unexpected role of bound anti-Schottky defect clusters in mediating defect transport. Remarkably, despite their poor thermal stability against dissociation at high temperatures, the transient formation of bound anti-Schottky defects under irradiation and their subsequent migration provide the dominant mechanism for the growth of large interstitial loops that have been experimentally observed in ThO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability and mobility of tungsten clusters on tungsten (110) surface: Ab initio and atomistic simulations

Quantifying the surface transport properties of tungsten (W) is of prime importance to understand the formation of nano-fuzz in fusion plasma-facing conditions. The stability and mobility of W adatom clusters (W n , n = 2-9) on the W(110) surface has been investigated by computer simulations, including ab initio calculations using density functional theory (DFT) and molecular statics (MS) simulations with multiple W interatomic potentials. The DFT results demonstrate that the sequential binding energy generally increases with number of W adatoms, except for the 5 th and 7 th W adatoms. The most common elemental migration steps of W n (n>2) clusters are observed to consist of monomer and dimer hops, while larger W n clusters can also diffuse by dissociation and recombination of smaller clusters. The threshold migration energy of W 9 is the highest, then followed by W 8 , W 4 , and W 6 , while W 3 , W 5 , and W 7 have similar migration energies. Compared to DFT, each interatomic potential evaluated overestimates the binding energies of W n clusters. Additionally, the embedded-atom potential developed by Juslin and Wirth adequately predicts the threshold migration energy of W n (n>2) clusters, although it predicts different underlying migration mechanisms. The results show that interaction mechanism between W adatoms controls the stability and mobility of W n clusters on the W(110) surface.

36 MATERIALS SCIENCE↗