Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-fidelity modeling of breakdown in helium: initiation processes and secondary electron emission

Understanding the role of physical processes contributing to breakdown is critical for many applications in which breakdown is undesirable, such as capacitors, and applications in which controlled breakdown is intended, such as plasma medicine, lightning protection, and materials processing. The electron emission from the cathode is a critical source of electrons which then undergo impact ionization to produce electrical breakdown. In this study, the role of secondary electron yields due to photons (γ ph ) and ions (γ i ) in direct current breakdown is investigated using a particle-in-cell direct simulation Monte Carlo model. The plasma studied is a one-dimensional discharge in 50 Torr of pure helium with a platinum cathode, gap size of 1.15 cm, and voltages of 1.2–1.8 kV. The current traces are compared with experimental measurements. Larger values of γ ph generally result in a faster breakdown, while larger values of γ i result in a larger maximum current. The 58.4 nm photons emitted from He(2 1 P) are the primary source of electrons at the cathode before the cathode fall is developed. Of the values of γ ph and γ i investigated, those which provide the best agreement with the experimental current measurements are γ ph = 0.005 and γ i = 0.01. In this work, these values are significantly lower than those in the literature for pristine platinum or for a graphitic carbon film which we speculate may cover the platinum. This difference is in part due to the limitations of a one-dimensional model but may also indicate surface conditions and exposure to a plasma can have a significant effect on the secondary electron yields. The effects of applied voltage and the current produced by a UV diode which was used to initiate the discharge, are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

A kinetic study of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas

We present particle-in-cell simulations with Monte Carlo collisions of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas. We study the energy balance in the one-dimensional expansion of a hot-spot by simulating Coulomb collisions, fusion reactions, and bremsstrahlung emission with a Monte Carlo model and inverse bremsstrahlung absorption using a new PIC model. This allows us to self-consistently capture the alpha particle heating and radiative losses in the expanding hot-spot and surrounding cold fuel. After verifying our model in a code-to-code comparison with both kinetic and fluid codes for the case of a deuterium–tritium hot-spot, we simulate the expansion of a proton–boron hot-spot initialized at 200 keV and 1,000 g/cm 3 . Our model predicts that energy radiated by the hot-spot is recaptured by the surrounding high-density opaque fuel reducing the expansion work done by the propagating burn wave. As a result, we find the net fusion energy produced over the course of $20$~ps is twice the initial hot-spot energy independent of whether radiation physics is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗
DOE OSTI2025

Mechanism of Mesoscale Woodpile Development via Photoelectrochemical Deposition of Se–Te

A combination of experiments and optical modeling provided insight into the mechanism of mesoscale woodpile formation in response to an orthogonal shift in polarization during photoelectrochemical deposition of Se–Te. Cathodic deposition of semiconducting Se–Te using spatially uniform, linearly polarized illumination produced arrays of lamellae that were aligned parallel to the optical E-field oscillation. Continued deposition in conjunction with an orthogonal shift in the polarization direction then produced aligned bridging features that spanned the void space between, and were orthogonal to, the preexisting lamellae. The height and pitch, respectively, in each layer of the woodpile were a function of the charge density and illumination wavelength during deposition. A Monte Carlo model, in which material addition was scaled by the absorption magnitude obtained from electromagnetic simulations, produced morphologies that were nominally identical to those observed experimentally. Here, the formation of mesoscale woodpiles is consistent with a mechanism that involves a series of spontaneously initiated, concerted light–matter interactions during the photoelectrochemical deposition process.

absorption↗
DOE OSTI2022

Intra-Specific Variability in Plant Hydraulic Parameters Inferred From Model Inversion of Sap Flux Data

Understanding plant hydraulic regulation is critical for predicting plant and ecosystem responses to projected increases in drought stress. Plant hydraulic regulation is controlled by observable, diverse plant hydraulic traits that can vary as much across individuals of the same species as they do across different species. Direct measurements of plant hydraulic traits from a range of ecosystems remain limited in comparison to other, more readily measured traits (e.g., specific leaf area). Furthermore, plant hydraulic trait measurements, often made at leaf or branch levels, are not easily scaled to whole-plant values that are typically used to predict plant and ecosystem fluxes. In this study, multiple whole-plant hydraulic parameters are inferred from observations of plant water use (i.e., sap flow), soil properties, and meteorological data. We use a Markov Chain Monte Carlo model inversion approach to obtain the best estimates and uncertainty of plant hydraulic parameters that capture whole-plant effective embolism resistance and stomatal sensitivity to decreasing plant water potential. We then use the inferred values in the model to estimate whole-tree water use and isohydricity. This approach reliably infers whole-plant parameter values with enough specificity to resolve inter- and intra-specific differences, and thus supplements time- and labor-intensive direct measurements of traits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Cometary Activity Begins at Kuiper Belt Distances: Evidence from C/2017 K2

We study the development of activity in the incoming long-period comet C/2017 K2 over the heliocentric distance range 9 ≲ r {sub H} ≲ 16 au. The comet continues to be characterized by a coma of submillimeter-sized and larger particles ejected at low velocity. In a fixed co-moving volume around the nucleus we find that the scattering cross section of the coma, C, is related to the heliocentric distance by a power law, C∝r{sub H}{sup -s}, with heliocentric index s = 1.14 ± 0.05. This dependence is significantly weaker than the r {sub H} {sup -2} variation of the insolation as a result of two effects. These are, first, the heliocentric dependence of the dust velocity and, second, a lag effect due to very slow-moving particles ejected long before the observations were taken. A Monte Carlo model of the photometry shows that dust production beginning at r {sub H} ~ 35 au is needed to match the measured heliocentric index, with only a slight dependence on the particle size distribution. Mass-loss rates in dust at 10 au are of order 10{sup 3} kg s{sup -1}, while loss rates in gas may be much smaller, depending on the unknown dust to gas ratio. Consequently, the ratio of the nongravitational acceleration to the local solar gravity, α', may, depending on the nucleus size, attain values of ~10{sup -7} ≲ α' ≲ 10{sup -5}, comparable to values found in short-period comets at much smaller distances. Nongravitational acceleration in C/2017 K2 and similarly distant comets, while presently unmeasured, may limit the accuracy with which we can infer the properties of the Oort cloud from the orbits of long-period comets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Solar Neutrino Detection with a Pixelated Liquid-Argon Time Projection Chamber

This thesis presents a study of low-energy solar neutrino detection using large-scale LArTPCs, focussing on novel pixelated readout technologies. Solar neutrinos offer a unique probe of fundamental neutrino properties and solar physics, but their detection in the MeV range is challenged by backgrounds. We investigate two complementary technologies: SoLAr, which integrates LArPix-based pixelated charge collection with Silicon Photomultipliers (SiPMs) in a hybrid anode design for simultaneous charge and light detection; and Q-Pix, a triggerless pixelated readout architecture based on charge integrate-reset circuits with local clocks, where Reset Time Differences encode ionisation waveforms via time-to-charge conversion. Two SoLAr prototypes were developed and operated, demonstrating VUV-sensitive SiPM performance in liquid argon and accurate charge-light signal matching with a charge detection threshold of $\sim 100 \mathrm{keV}$. We also implement a complete simulation and reconstruction framework, incorporating realistic detector geometry, electron transport, readout response, and detailed signal and background models, including intrinsic argon and radon progeny, as well as site-specific $\gamma$-ray and neutron fluxes. For Q-Pix, we demonstrate that with a pixel size of $4\times 4$ mm$^2$ and a reset threshold of 1 fC ($\sim 0.1475$ MeV), full-scale operation produces data volumes below 1 PB per 10 ktonne-year. For SoLAr, assuming a shielded DUNE-like detector and 100 ktonne-year exposure, we project uncertainties of $0.90\times 10^{-5}$ eV$^2$ on $\Delta m^2_{21}$ and 0.033 on $\sin^2\theta_{12}$, improving to $0.46\times 10^{-5}$ eV$^2$ and 0.025 with 400 kilotonne-year. At this higher exposure, we also obtain a day–night flux asymmetry at the level of $( -5.6 \pm 3.6 ) \%$. Combining Monte Carlo modelling, hardware validation, and advanced reconstruction techniques, this work establishes a path toward next-generation ktonne-scale LArTPCs as observatories for precision solar neutrino physics.

Ruiz Ferreira, Guilherme [Manchester U.] (ORCID:00↗
DOE OSTI2021

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

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

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Study of dijet events with large rapidity separation in proton-proton collisions at $\sqrt{s}$ = 2.76 TeV

The cross sections for inclusive and Mueller-Navelet dijet production are measured as a function of the rapidity separation between the jets in proton-proton collisions at $\sqrt{s}$ = 2.76 TeV for jets with transverse momentum p T > 35 GeV and rapidity |y| < 4.7. Various dijet production cross section ratios are also measured. A veto on additional jets with p T > 20 GeV is introduced to improve the sensitivity to the effects of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution. The measurement is compared with the predictions of various Monte Carlo models based on leading-order and next-to-leading-order calculations including the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi leading-logarithm (LL) parton shower as well as the LL BFKL resummation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Characterization of Two Extraordinary AmLi Neutron Sources

Two legacy AmLi sources, MRC-101 and MRC-103, are routinely used for neutron detector characterization and testing campaigns that Pacific Northwest National Laboratory manages for the Defense Threat Reduction Agency within the U.S. Department of Defense. However, applicability of source measurement data has been substantially limited due to the lack of original manufacturer documentation and cross comparison to other sources. This paper summarizes results of the measurement campaign conducted to better quantify source radionuclide content, neutron energy spectrum and emission rate of the AmLi sources. The source characterization study consisted of three major consecutive measurement exercises. First, gamma spectrometry measurements using a high purity germanium detector were used to estimate the radionuclide content and potential impurities of the AmLi sources. Then neutron spectral measurements were taken to determine the unique neutron energy distribution. Finally, assessments of the total neutron yields of the sources were made with measurements in a neutron well counter. The neutron detection efficiency for these measurements was determined via detailed Monte Carlo modeling based on the neutron energy distribution measured with the rotating neutron spectrometer. The most exciting and unexpected finding of the study was a significant contamination of the source material with beryllium, contributing up to 17% to the total neutron emission rate and correspondingly changing the neutron energy spectrum. Also, the obtained lithium-only portion of the measured energy distribution of the sources was compared against several AmLi spectra previously reported by different researchers.

AWCC, AmLi, beryllium contamination, neutron spect↗
DOE OSTI2024

Measurements of Lund subjet multiplicities in 13 TeV proton-proton collisions with the ATLAS detector

This Letter presents a differential cross-section measurement of Lund subjet multiplicities, suitable for testing current and future parton shower Monte Carlo algorithms. This measurement is made in dijet events in 140 fb -1 of $\sqrt{s}$ =13 TeV proton–proton collision data collected with the ATLAS detector at CERN's Large Hadron Collider. The data are unfolded to account for acceptance and detector-related effects, and are then compared with several Monte Carlo models and to recent resummed analytical calculations. The experimental precision achieved in the measurement allows tests of higher-order effects in QCD predictions. Most predictions fail to accurately describe the measured data, particularly at large values of jet transverse momentum accessible at the Large Hadron Collider, indicating the measurement's utility as an input to future parton shower developments and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Unraveling the Early-Stage Ordering of Krypton Solid Bubbles in Molybdenum: A Multimodal Study

Self-organization of defects such as fission gas bubbles in materials can lead to high inventory capacity for fission gas storage and help mitigate swelling caused by fission gases in nuclear fuel materials under radiation in nuclear reactors. Here, we report the physical mechanism of self-organization of krypton (Kr) gas bubbles in molybdenum (Mo) under ion implantation. The ion fluence and temperature-dependent formation of Kr solid bubble superlattice (SBS) in Mo were investigated by using both synchrotron-based small-angle X-ray scattering (SAXS) and transmission electron microscopy (TEM). Early stage self-organization of gas bubbles is observed at a fluence of 2.5 × 10 16 ions/cm 2 at temperatures of 300–400 °C. The bubble lattice constant increases with increasing implantation temperature from 300 to 400 °C. Both experiments and atomic kinetic Monte Carlo modeling indicate that the Kr solid bubbles are weakly ordered in comparison to previously studied helium (He) gas bubble superlattice (GBS) while the lattice constant are relatively smaller for Kr SBS compared to that of He GBS. The irradiation conditions suggest that spinodal decomposition, which is a form of phase separation, probably precedes gas bubble ordering in Mo. Altogether, our work sheds light on the formation mechanism of noble gas superlattice toward the development of radiation-tolerant materials which are important for the design of advanced nuclear reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Dissociation Pathways of BCl 3 on the Si(100)-2 × 1 Surface

BCl 3 is a promising acceptor precursor for atomic-precision δ-doping of silicon, as it has been observed to rapidly dissociate into boron doped into the silicon surface and surface chlorine, which can be removed upon annealing. The chemical pathway and the resulting kinetics, through which BCl 3 adsorbs and dissociates on silicon, however, have only been partially explained. Here, in this work, we use density functional theory to expand the dissociation reactions of BCl 3 to include reactions that take place across multiple silicon dimer rows and reactions which end in a bare B atom either at the surface, substituted for a surface silicon, or in a subsurface position. We further simulate the resulting scanning tunneling microscopy images for each of these BCl x dissociation fragments, demonstrating that they often display distinct features that may allow for relatively confident experimental identification. Finally, we input the full dissociation pathway for BCl 3 into a kinetic Monte Carlo model, which simulates realistic reaction pathways as a function of environmental conditions, such as the pressure and temperature of dosing. We find that BCl 2 is broadly dominant at low temperatures, while high temperatures and ample space on the silicon surface for dissociation encourage the formation of bridging BCl fragments and B substitutions on the surface. This work provides the chemical mechanisms for understanding atomic-precision doping of Si with B, enabling a number of relevant quantum applications, such as bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.

Campbell, Quinn T. [Sandia National Laboratories (↗
DOE OSTI2023

Unexpected Thermomechanical Behavior of Off-Stoichiometry Epoxy/Amine Materials

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure–property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate T g . Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in T g , suggesting that spatiotemporal contributions-not captured by our model-were significant determinants of material properties compared to polymer architecture for three-component systems. Here these findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Overcoming the Entropy Penalty of Direct Air Capture for Efficient Gigatonne Removal of Carbon Dioxide

Atmospheric carbon poses an existential threat to civilization via global climate change. Hundreds of gigatonnes of carbon dioxide must be removed from earth’s atmosphere in the next three decades, necessitating a low-cost, energy-efficient process to extract low concentrations of carbon dioxide for conversion to a stable material permanently stored for thousands of years. In this work, the challenge of removing gigatonnes of CO 2 is described via the scale of effort and the thermodynamics of collecting and reducing this diffuse chemical, the accumulation of which imparts a substantial entropy penalty on any atmospheric carbon capture process. The methods of CO 2 reduction combined with upstream direct air capture (DAC) including absorption, membrane separation, and adsorption are compared with biomass torrefaction and permanent burial (BTB). A Monte Carlo model assesses the mass, energy, and economics of the full process of biomass torrefaction from biomass collection and transport to stable carbon burial to determine that 95% of scenarios could remove carbon for less than $200 per CO 2 -tonne-equivalent. Torrefied carbon is further discussed for its long-term stability and availability at the scale required to substantially mitigate the threat of climate change.

biomass↗
DOE OSTI2024

Amorphous Oxyhalide Matters for Achieving Lithium Superionic Conduction

The recently surged halide-based solid electrolytes (SEs) are great candidates for high-performance all-solid-state batteries (ASSBs), due to their decent ionic conductivity, wide electrochemical stability window, and good compatibility with high-voltage oxide cathodes. In contrast to the crystalline phases in halide SEs, amorphous components are rarely understood but play an important role in Li-ion conduction. Here, we reveal that the presence of amorphous component is common in halide-based SEs that are prepared via mechanochemical method. The fast Li-ion migration is found to be associated with the local chemistry of the amorphous proportion. Taking Zr-based halide SEs as an example, the amorphization process can be regulated by incorporating O, resulting in the formation of corner-sharing Zr–O/Cl polyhedrons. This structural configuration has been confirmed through X-ray absorption spectroscopy, pair distribution function analyses, and Reverse Monte Carlo modeling. The unique structure significantly reduces the energy barriers for Li-ion transport. As a result, an enhanced ionic conductivity of (1.35 ± 0.07) × 10 –3 S cm –1 at 25 °C can be achieved for amorphous Li 3 ZrCl 4 O 1.5 . In addition to the improved ionic conductivity, amorphization of Zr-based halide SEs via incorporation of O leads to good mechanical deformability and promising electrochemical performance. In conclusion, these findings provide deep insights into the rational design of desirable halide SEs for high-performance ASSBs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗