Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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

Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization

Abstract A fundamental understanding of the interaction between microstructure and underlying physical mechanisms is essential, especially for developing more accurate multi-physics models for heterogeneous materials. Effects of microstructure on the material response at the macroscale are modeled by using the generalized thermomechanics. In this study, strain gradient theory is employed as a higher-order theory on the macroscale with thermodynamics modeled as a first-order theory on the microscale. Hence, energy depends only on the temperature such that we circumvent an extension of Fourier’s law and analyze the “simplest” thermo-mechanical model in strain gradient elasticity. Developing multiphysics models for heterogeneous materials is indeed a challenge and even this “simplest” model in generalized thermomechanics creates dozens of parameters to be determined. We develop a thermo-mechanical framework, in which microstructure is modeled as a periodic structure and through asymptotic homogenization approach, higher-order parameters at macroscopic scale are calculated. To illustrate the importance of higher-order parameters in overall thermo-mechanical response of a heterogeneous materials, finite element method (FEM) is employed with the aid of open-source codes (FEniCS). Verification example of a bulk system and several case studies of porous structures demonstrate how such numerical framework can be beneficial in the design of materials with tailored microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Cascade electrocatalysis via AgCu single-atom alloy and Ag nanoparticles in CO 2 electroreduction toward multicarbon products

Electrocatalytic CO 2 reduction into value-added multicarbon products offers a means to close the anthropogenic carbon cycle using renewable electricity. However, the unsatisfactory catalytic selectivity for multicarbon products severely hinders the practical application of this technology. In this paper, we report a cascade AgCu single-atom and nanoparticle electrocatalyst, in which Ag nanoparticles produce CO and AgCu single-atom alloys promote C-C coupling kinetics. As a result, a Faradaic efficiency (FE) of 94 ± 4% toward multicarbon products is achieved with the as-prepared AgCu single-atom and nanoparticle catalyst under ~720 mA cm -2 working current density at -0.65 V in a flow cell with alkaline electrolyte. Density functional theory calculations further demonstrate that the high multicarbon product selectivity results from cooperation between AgCu single-atom alloys and Ag nanoparticles, wherein the Ag single-atom doping of Cu nanoparticles increases the adsorption energy of *CO on Cu sites due to the asymmetric bonding of the Cu atom to the adjacent Ag atom with a compressive strain.

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

High-Burnup LOCA Burst Susceptibility BISON Analysis in PWRs and BWRs

Accurately assessing high-burnup fuel behavior during loss-of-coolant accidents (LOCAs) is essential for understanding fuel fragmentation, relocation, and dispersal (FFRD) risks across the US light-water reactor fleet. This work updates previous Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program multiphysics LOCA analyses for a pressurized water reactor (PWR) and a boiling water reactor (BWR) by incorporating recent model and material property advancements in the BISON fuel performance code, including a high-burnup structure (HBS) model, revised cladding burst criteria, and updated thermal–mechanical correlations. This update was needed to support ongoing industry initiatives and upcoming regulatory changes. Full-core, rod-resolved operating histories generated using Virtual Environment for Reactor Analysis (VERA) and system-level LOCA conditions obtained from TRACE were applied to statistically representative rod samples in BISON to evaluate burst behavior and FFRD susceptibility. These calculations used two cladding burst correlations and three fuel pulverization models so that the predictions of these models could be compared. The updated PWR simulations show markedly improved numerical stability as the number of crashed simulations decreased by 95% compared to the previous study, and hence higher confidence in results. The updated PWR simulations predicted cladding bursts exclusively among once-burned, high-power rods, with two different cladding burst models identifying the same burst-susceptible population. Resulting FFRD susceptibility estimates are significantly reduced compared with earlier studies, driven by cooler predicted fuel and plenum temperatures, lower hoop strains, and reduced fission gas release in the updated models. In contrast, none of the BWR rods were predicted to burst under either burst criterion, reaffirming minimal BWR FFRD susceptibility even with updated HBS and material models. Comparisons between the PWR and BWR end-of-cycle predictions are made. Comparison with prior work highlights significant shifts in PWR fuel performance metrics and confirmation of earlier BWR conclusions. Overall, the updated results underscore the importance of having high-resolution detailed modeling capability and continuously integrating evolving material models and physics into high-resolution multiphysics simulations. The unified assessment presented here strengthens confidence in predicting high-burnup LOCA behavior by improving agreement between different cladding burst correlations. These results also provide an improved foundation for future BISON model development, FFRD susceptibility calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Ultrafast photo-induced phonon hardening due to Pauli blocking in MAPbI 3 single-crystal and polycrystalline perovskites

Metal-halide perovskite semiconductors have attracted intense interest over the past decade, particularly for applications in photovoltaics. Low-energy optical phonons combined with significant crystal anharmonicity play an important role in charge-carrier cooling and scattering in these materials, strongly affecting their optoelectronic properties. We have observed optical phonons associated with Pb–I stretching in both MAPbI 3 single crystals and polycrystalline thin films as a function of temperature by measuring their terahertz conductivity spectra with and without photoexcitation. An anomalous bond hardening was observed under above-bandgap illumination for both single-crystal and polycrystalline MAPbI 3 . First-principles calculations reproduced this photo-induced bond hardening and identified a related lattice contraction (photostriction), with the mechanism revealed as Pauli blocking. For single-crystal MAPbI 3 , phonon lifetimes were significantly longer and phonon frequencies shifted less with temperature, compared with polycrystalline MAPbI 3 . We attribute these differences to increased crystalline disorder, associated with grain boundaries and strain in the polycrystalline MAPbI 3 . Thus we provide fundamental insight into the photoexcitation and electron–phonon coupling in MAPbI 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Instability of Rock-Salt Cubic NbN in Density Functional Calculations

All-nitride semiconductor/superconductor heterojunctions utilizing cubic niobium nitride (NbN) are a promising approach to superconducting quantum circuits for next-generation quantum-information systems. However, there are fundamental open questions on the atomic structure of NbN. In calculations with several levels of density functional theory (DFT), we find that the cubic rocksalt structure NbN (Fm-3m, 225) is energetically very unstable against the ground state hexagonal NbN in a tungsten carbide (P-6m2, 187) type lattice. To better understand the appearance of a cubic phase in numerous experiments, we perform a DFT study on possible NbN structures, determining the energy ordering between different polymorphs from databases and structure prediction. We perform supercell calculations of disordered NbN, finding that the rocksalt structure is dynamically unstable and relaxes to a lower energy monoclinic phase (C2/m, 12), which retains an approximate average cubic symmetry. However, the associated energy gain is not substantial enough for a plausible explanation for cubic NbN. We further investigate the role of external factors such as in-plane strain during epitaxial growth as well as the presence of off-stoichiometry and impurity doping on the energy ordering.

CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SU↗
DOE OSTI2025

Electronic Excitation‐Driven β‐Ga 2 O 3 Metastability Transformation and Self‐Organization Mechanism: β→κ/γ/δ Phases

Here, irradiation-driven multiphase self-organization presents emergent opportunities for the customization of nanoscale engineering properties, dynamically tuning strain-field distributions and interfacial electronic structures. Responding to intense electronic excitation-induced energy deposition, the dominant phase transformations, with varying Gibbs free energy $Δ,G^o_f$ are confirmed as β → κ → γ → δ that are located in specific microregions for Gallium (III) oxide (Ga 2 O 3 ), as follows: (i) Surface-localized interstitial accumulation under compressive stress triggers β → δ via semi-coherent interface formation. (ii) Tensile stress within latent tracks drives vacancy-mediated oxygen layer truncation (4/12 periodicity along ⟨0001⟩), stabilizing coherent 4H (ABCB) κ and 3C (ABC) β (ABC) interfaces through strain-compensated octahedral distortion. (iii) Screw dislocation-mediated lattice relaxation induces β → γ via cation disordering (Ga 3+ occupancy at β-interstitial sites), forming metastable spinel γ with mixed occupancy across 16d/8a Wyckoff sites. Irradiation-driven β-Ga 2 O 3 →κ/γ/δ transitions, as mechanistically revealed via inelastic thermal spike (i-TS) calculations and molecular dynamics simulations, induce defect-mediated nonlinear photoresponse, critical for optoelectronic engineering.

Han, Xinqing [Shandong Univ., Jinan (China)]↗
DOE OSTI2021

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chapter 7.1: Microalgae

Microalgae is a unique biomass resource that does not need to compete for land and water with other biomass feedstocks because it can be cultivated on low-quality unencumbered land using noncompetitive water types including saline and wastewater. It is included as a complementary resource alongside other biomass feedstocks reported in this study, albeit at higher biomass production costs reflective of more capital-intensive farming operations than typical for terrestrial biomass. Higher biomass costs can be offset by the potential to produce value-added coproducts unique to compositional constituents of microalgae. Relative to BT16, this chapter reflects the latest analysis from the 2022 Algae Harmonization Update, which uses the latest parameterized and high-performing saline algal strain, second-generation carbon capture of point-source waste CO2 and high-pressure pipeline transport resolved to specific point-source types, saline water sourcing up to 40,000 mg/L total dissolved solids for source and makeup water salinity, blowdown water treatment and recycle, and brine disposal handling. National-scale algal biomass availability potential was calculated at 152 million tons/yr ash-free dry weight (AFDW) (191 million tons/yr dry weight) at an average biomass productivity of 26.2 g/m2 -day AFDW (about 50 tons/acre/yr dry weight).1 The algal biomass was cultivated on 3.9 million acres of multi-criteria screened and potentially available land for CONUS and fixed 268 million tons of waste CO2 based on biomass uptake. The algae biomass can be produced at an average MBSP of $674/ton AFDW ($536/ton dry weight) in 2020 dollars,2 corresponding to a total energy potential of 3.3 quads/yr at an average MBSP of $31.2/MMBtu (higher-heating-value [HHV] basis).

algae↗
DOE OSTI2021

X-ray diffraction and in situ pressurization of dentine apatite reveals nanocrystal modulus stiffening upon carbonate removal

Bone-like materials comprise carbonated-hydroxyapatite nanocrystals (c-Ap) embedding a fibrillar collagen matrix. The mineral particles stiffen the nanocomposite by tight attachment to the protein fibrils creating a high strength and toughness material. The nanometer dimensions of c-Ap crystals make it very challenging to measure their mechanical properties. Mineral in bony tissues such as dentine contains 2~6 wt.% carbonate with possibly different elastic properties as compared with crystalline hydroxyapatite. Here we determine strain in biogenic apatite nanocrystals by directly measuring atomic deformation in pig dentine before and after removing carbonate. Transmission electron microscopy revealed the platy 3D morphology while atom probe tomography revealed carbon inside the calcium rich domains. High-energy X-ray diffraction in combination with in situ hydrostatic pressurization quantified reversible c-Ap deformations. Additionally, crystal strains differed between annealed and ashed (decarbonated) samples, following 1 or 10 h heating at 250 °C or 550 °C respectively. Measured bulk moduli (K) and a-/c-lattice deformation ratios (η) were used to generate synthetic K syn and ηsyn identifying the most likely elastic constants C 33 and C 13 for c-Ap. These were then used to calculate the nanoparticle elastic moduli. For ashed samples, we find an average E 11 =107 GPa and E 33 =128 GPa corresponding to ~5% and ~17% stiffening of the a-/c-axes of the nanocrystals as compared with the biogenic nanocrystals in annealed samples. Ashed samples exhibit ~10% lower Poisson's ratios as compared with the 0.25~0.36 range of carbonated apatite. Carbonate in c-Ap may therefore serve for tuning local deformability within bony tissues.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Barrier Heights for Diels–Alder Transition States Leading to Pentacyclic Adducts: A Benchmark Study of Crowded, Strained Transition States of Large Molecules

Theoretical characterization of reactions of complex molecules depends on providing consistent accuracy for the relative energies of intermediates and transition states. Here we employ the DLPNO-CCSD(T) method with core–valence correlation, large basis sets, and extrapolation to the CBS limit to provide benchmark values for Diels–Alder transition states leading to competitive strained pentacyclic adducts. We then used those benchmarks to test a diverse set of wave function and density functional methods for the absolute and relative barrier heights of these transition states. Our results show that only a few of the tested density functionals can predict the absolute barrier heights satisfactorily, although relative barrier heights are more accurate. The most accurate functionals tested are ωB97M-V, M11plus, ωB97X-V, PBE-D3(0), M11, and MN15 with MUDs from best estimates less than 3.0 kcal. Furthermore, these findings can guide selection of density functionals for future studies of crowded, strained transition states of large molecules.

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

Small Molecule Inhibitor-Modulated Al 2 O 3 Atomic Layer Deposition on Monolayer MoS 2 for Controlled Nucleation

The integration of ultrathin dielectrics on two-dimensional (2D) semiconductors is essential for advancing beyond-Si electronics. However, the intrinsic inertness of van der Waals 2D basal planes remains a primary bottleneck to achieving uniform dielectric nucleation and growth. Here, in this study, we introduce a small molecule inhibitor (SMI)-modulated thermal atomic layer deposition (ALD) strategy, exemplified by aluminum oxide (Al 2 O 3 ) ALD on monolayer molybdenum disulfide (1L MoS 2 ) with acetic acid (HAc) SMI. The ABC-type sequence comprises HAc inhibitor (A), trimethylaluminum (TMA) precursor (B), and deionized H 2 O coreactant (C). In situ quartz crystal microbalance (QCM) studies reveal robust HAc adsorption on Al 2 O 3 and suppression of subsequent oxide growth on HAc-passivated surfaces. When applied to 1L MoS 2 , this inhibitory pathway enables HAc to selectively passivate nascent Al 2 O 3 nuclei formed on the MoS 2 surface, limiting their three-dimensional (3D) island coarsening and redirecting precursor adsorption toward the uncovered basal plane. Consequently, nearly continuous ultrathin (∼1.5 nm) Al 2 O 3 films are achieved on 1L MoS 2 with markedly improved uniformity compared to standard Al 2 O 3 ALD using TMA and H 2 O, as validated by atomic force microscopy (AFM), cross-sectional scanning transmission electron microscopy (STEM), and energy-dispersive X-ray spectroscopy (EDS). Density functional theory (DFT) calculations further provide atomistic insight into HAc-modulated Al 2 O 3 nucleation, corroborating the energetic preference of HAc for Al 2 O 3 over MoS 2 and attenuated TMA adsorption on HAc-passivated surfaces. Spatially resolved Raman spectroscopy also confirms that the HAc-modulated process preserves the structural integrity of 1L MoS 2 , with only minimal strain and doping perturbations observed after dielectric deposition. This SMI-modulated approach offers a broadly applicable framework for controlling ALD nucleation across various inhibitors, ALD chemistries, and 2D materials, opening opportunities for reliable dielectric integration in next-generation nanoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Semi-Dirac and Weyl fermions in transition metal oxides

We illustrate that a class of compounds with I4/mcm crystalline symmetry hosts three-dimensional semi-Dirac fermions. Unlike the known two-dimensional semi-Dirac points, the degeneracy of these three-dimensional semi-Dirac points is not lifted by spin-orbit coupling due to the protection by a nonsymmorphic symmetry—mirror reflection in the a–b plane and a translation along the c axis. This crystalline symmetry is found in tetragonal perovskite oxides, realizable in thin films by epitaxial strain that results in a 0 a 0 c - -type octahedral rotation. Interestingly, with broken time-reversal symmetry, two pairs of Weyl points emerge from the semi-Dirac points within the Brillouin zone, and an additional lattice distortion leads to an enhanced intrinsic anomalous Hall effect. The ability to tune the Berry phase by epitaxial strain can be useful in novel oxide-based electronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Controllable topological insulator phases in litharge-phase InBi monolayer

Despite recent advances of layered square-net topological material models that possess ideal semimetallic electronic structures and promising potential in material applications, the identification of experimentally accessible two-dimensional square-net materials with related topological properties has proven challenging. Due to the highly tunable physical and topological properties of III-V semiconductors, we revisit the class of III-V materials and observe that the litharge-phase InBi is a layered square-net material and can be exfoliated into the InBi monolayer. We present a comprehensive first-principles study of the energy landscape of the InBi monolayer. We identify a paraelastic phase and three ferroelastic phases and study their topological properties. Specifically, we show that the paraelastic InBi monolayer is a trivial insulator due to the orbital-ordering-induced band inversion occurring between states with the same parity. Substituting one Bi atom per cell with another V-group element (N, P, As) or applying an electric field that breaks the inversion symmetry and changes the orbital onsite energy, the paraelastic InBi monolayer can be driven into the topological insulator phase. Furthermore, one of the ferroelastic phases of pure InBi, which can be obtained by gently straining the paraelastic phase, also possesses such topological insulating properties. Furthermore, these results provide several experimentally accessible routes to tune the nontrivial topology in the InBi monolayer, including creating heterostructures with piezoelectric or ferroelectric substrates and applying mechanical strain, making the InBi monolayer an ideal platform to study the interplay of reduced dimensionality, square-net chemical bonding networks, and band topology.

2-dimensional systems↗
DOE OSTI2024

2023 Billion-Ton Report: An Assessment of U.S. Renewable Carbon Resources - Microalgae

Microalgae is a unique biomass resource that does not need to compete for land and water with other biomass feedstocks because it can be cultivated on low-quality unencumbered land using noncompetitive water types including saline and wastewater. It is included as a complementary resource alongside other biomass feedstocks reported in this study, albeit at higher biomass production costs reflective of more capital-intensive farming operations than typical for terrestrial biomass. Higher biomass costs can be offset by the potential to produce value-added coproducts unique to compositional constituents of microalgae. Relative to the 2016 Billion-Ton Report, this chapter reflects the latest analysis from the 2022 Algae Harmonization Update, which uses the latest parameterized and high-performing saline algal strain, second-generation carbon capture of point-source waste CO 2 , and high-pressure pipeline transport resolved to specific point-source types, saline water sourcing up to 40,000 mg/L total dissolved solids for source and makeup water salinity, blowdown water treatment and recycle, and brine disposal handling. National-scale algal biomass availability potential was calculated at 152 million tons/yr ash-free dry weight (AFDW) (191 million tons/yr dry weight) at an average biomass

09 BIOMASS FUELS↗
DOE OSTI2020

An algorithm for resolving intragranular orientation fields using coupled far-field and near-field high energy $\mathrm{X}$-ray diffraction microscopy

Here, we present a novel algorithm for reconstructing spatial intragranular lattice orientation fields using both far-field and near-field high energy X-ray diffraction microscopy (HEDM) measurements. An established far-field indexing algorithm is modified to include lattice orientation distribution information (grain orientation envelopes) in addition to average grain orientations. The near-field data reduction algorithm utilizes this enriched far-field orientation data as a seed for reconstructing intragranular spatial maps of lattice orientation from the diffraction images. The primary benefit of the new algorithm is a significant decrease in the number of trial calculations that must be performed in the reconstruction process compared to existing methodologies while maintaining scalability. The resulting gains in efficiency facilitate the use of relatively modest computational resources and improve throughput at the point of measurement. We provide two example applications: volumetric orientation field reconstructions for a Ti-Al alloy both before and after the application of 3% uniaxial strain. The results showcase the efficiency of the new method and the ability to resolve subtle changes in microstructure, which are associated with incipient plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural and optical characterization of dilute Bi-doped GaN nanostructures grown by molecular beam epitaxy

We have carried out detailed studies on the epitaxy and characterization of dilute Bi-doped GaN nanostructures. A comprehensive investigation of Bi-doped GaN nanowires and quasi-film epitaxial growth conditions has been performed. Scanning electron microscopy studies show that lowering the GaBiN growth temperature causes gradual changes in top c-plane nanowire morphology due to the incremental incorporation of foreign Bi atoms. This trend is further substantiated by the secondary ion mass spectroscopy analysis of a multi-layer Bi-doped GaN quasi-film. However, it is also found that the amount of Bi incorporation into the GaN lattice is relatively independent of the N2 flow rate variation under the growth conditions investigated. Furthermore, room-temperature micro-Raman spectra show that there are additional peaks near 530, 650, and 729 cm−1 wave numbers in the Bi-doped GaN samples, which can primarily be attributed to Bi local vibrational modes, indicative of a small amount of Bi incorporation in the GaN lattice. Moreover, phonon calculations with density functional theory indicate that Bi replacing the N sites is the likely origin of the experimentally measured Raman modes. X-ray photoelectron spectroscopy measurements have also been obtained to deduce the electronic interaction between the Bi dopant atom and the GaN nanostructure. Such one-dimensional nanowires permit the synthesis of dislocation-free highly mismatched alloys due to strain relaxation, allowing efficient light absorption and charge carrier extraction that is relevant for solar energy harvesting and artificial photosynthesis.

Materials Science↗
DOE OSTI2020

The Role of Local Inhomogeneities on Dendrite Growth in LLZO-Based Solid Electrolytes

The majority of the ceramic solid electrolytes (LLZO, LATP) demonstrate polycrystalline grain/grain-boundary (G/GB) microstructure. Higher lithium (Li) concentration and lower mechanical stiffness result in current focusing at the GBs. Growth of Li dendrites through local inhomogeneities and subsequent short circuit of the cell is a major concern. Recent studies have revealed that bulk Li metal is a viscoplastic material that has low (~0.3 MPa) and high (~1.0 MPa) yield strength during deformation at smaller and larger rates of strain, respectively. It has been argued that during deposition at smaller current densities, due to its lower yield strength, Li metal should demonstrate plastic flow against stiff ceramic electrolytes, and Li dendrites will be prevented from penetrating through solid electrolytes. In this manuscript, a multiscale modeling framework has been developed for predicting properties of GBs and the bulk of ceramic electrolytes using atomistic calculations for input to mesoscale models. Using the parameters obtained from the atomistic simulations, the mesoscale model reveals that, given enough time, even at low charge rates, lithium dendrites can grow through the GBs of LLZO. The present multiscale model results also provide information regarding the dendrite growth velocity through LLZO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗