Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Study”

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

Photoelectron Spectroscopic and ab Initio Computational Studies of the Anion, HThO –

The synergetic combination of anion photoelectron spectroscopy and high-level relativistic coupled-cluster calculations was employed to study the anion, HThO-. The atomic connectivity of this anion was found to be HThO- and not ThOH-. Vibrational and electronic energy spacings in the HThO- photoelectron spectrum were measured and calculated, with good agreement between them being found. Computations yielded electronic energies and equilibrium structures as well as enabling orbital analyses. The adiabatic electron affinity of HThO was determined to be 1.297 ± 0.035 eV.

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

Ligand effects on electronic structure and bonding in U(iii) coordination complexes: a combined MCD, EPR and computational study

The trivalent oxidation state of uranium has been shown to undergo unique reactivity, from its ability to activate a variety of small molecules to its role in the catalytic reduction of ethene to ethane amongst others. Central to this unique reactivity and ability to rationally design ligands for isotope separation is the underlying uranium electronic structure. While electronic structure studies of U(IV), U(V), and U(VI) have been extensive, by comparison, analogous studies of more reduced oxidation states such as U(III) remains underdeveloped. Herein we report a combined MCD and EPR spectroscopic approach along with density functional theory and multireference wavefunction calculations to elucidate the effects of ligand perturbation in three uranium(III) Tp* complexes. Overall, the experimental and computational insight suggests that the change in ligand environment across this series of U(III) complexes resulted in only minor perturbations in the uranium electronic structure. This combined approach was also used to redefine the electronic ground state of a U(III) complex with a redox non-innocent Bipy - ligand. Overall, these studies demonstrate the efficacy of the combined experimental and theoretical approach towards evaluating electronic structure and bonding in U(III) complexes and provide important insight into the challenges in altering ligand environments to modify bonding and reactivity in uranium coordination chemistry.

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

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A molecular beam and computational study on the barrierless gas phase formation of (iso)quinoline in low temperature extraterrestrial environments

Despite remarkable progress toward the understanding of the formation pathways leading to polycyclic aromatic hydrocarbons (PAHs) in combustion systems and in deep space, the complex reaction pathways leading to nitrogen-substituted PAHs (NPAHs) at low temperatures of molecular clouds and hydrocarbon-rich, nitrogen-containing atmospheres of planets and their moons like Titan have remained largely obscure. Here, we demonstrate through laboratory experiments and computations that the simplest prototype of NPAHs – quinoline and isoquinoline (C 9 H 7 N) – can be synthesized via rapid and de-facto barrier-less reactions involving o-, m- and p-pyridinyl radicals (C 5 H 4 N˙) with vinylacetylene (C 4 H 4 ) under low-temperature conditions.

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

Vibrational Spectra of HNIW and its Isotopologues: A Combined Experimental and Computational Study

Incorporating isotopically labelled materials in degradation experiments could help unravel the mechanism(s) of decomposition through use of the kinetic isotope effect. Characterizing synthesized isotopologues however requires an understanding of what observable signals are affected by the isotopic substitution. As vibrational spectroscopy can distinguish between isotopologues, it is an ideal characterization technique to evaluate isotopic variants. To this end, the vibrational spectra of HNIW and its deuterated ( 2 H), 13 C, 15 N (all), 15 N (nitro), and 18 O isotopologues have been computationally predicted in the gas phase using density functional theory. These results are compared to experimentally measured FTIR/ATR and Raman spectra of both unsubstituted HNIW and 15 N-labeled HNIW in which the six nitro groups were synthetically tagged with 15 N atoms ( 15 N nitro -HNIW). The experimental isotopic frequency shift for the -NO 2 asymmetric stretching frequencies agrees with that theoretically calculated (~35.7 cm-1 vs. 36.5cm-1, respectively). Furthermore, analysis of the theoretically predicted frequency shifts for all isotopologues suggest the -NO 2 bending modes are lower in frequency than previously reported. This assignment is supported by the experimentally measured isotopic shift of ~10.1 cm -1 for these features (consistent with the predicted shift of ~13.1 cm -1 ). This work expands our current understanding of the vibrational modes in HNIW as well as provides a method for future work on similar systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Why Are 5-Thioglycopyranosyl Donors More Axially Selective than their Glycopyranosyl Counterparts? A Low and Variable Temperature NMR Spectroscopy and Computational Study

5-Thioglycopyranosyl donors differ in reactivity and selectivity from simple glycopyranosyl donors. An extensive study has been conducted on the nature and stability of the reactive intermediates generated on the activation of per-O-acetyl and per-O-methyl 5-thioglucopyranosyl donors and the corresponding glucopyranosyl donors. Variable temperature nuclear magnetic resonance (NMR) studies with per-O-methylated or per-O-acetyl glycosyl sulfoxides and trichloroacetimidates on activation with trifluoromethanesulfonic anhydride or trimethylsilyl triflate are reported. These show that following initial adduct formation with the promoter conversion of the 5-thioglucopyranosyl donors to the 5-thioglucopyranosyl triflates requires higher temperatures than that of the glucopyranosyl donors to the glucopyranosyl triflates. It is demonstrated that neighboring group participation is a less important phenomenon for the peracetylated thioglucosyl donors than for the peracetylated glucosyl donors. A simple thiocarbenium ion was generated by protonation of 2,3-dihydro-4H-thiopyran at low temperature and characterized by NMR spectroscopy. However, the corresponding 5-thioglucopyranosyl thenium ions were not observed in any of the NMR studies of the 5-thiopyranosyl donors: the electron-withdrawing C–O bonds around the thiopyranoside core discourage thiocarbenium ion formation, just as they discourage oxocarbenium ion formation. Density functional theory (DFT) calculations reveal the tetrahydrothiopyranyl thiocarbenium ion to be approximately 2.5 kcal/mol lower in energy than the corresponding tetrahydropyranyl oxocarbenium ion relative to the corresponding covalent triflates. However, the computations reveal a 5.8 kcal/mol activation barrier for conversion of the tetrahydrothiopyranyl triflate to the thiocarbenium ion, while formation of the oxocarbenium ion–triflate ion pair from tetrahydropyranyl triflate requires only 2.6 kcal·mol –1 . Overall, the greater axial selectivity of 5-thioglycopyranosyl donors compared to analogous glycopyranosyl donors derives from (i) the lower kinetic reactivity necessitating higher reaction temperatures, (ii) the greater stability of the thiocarbenium ion over the oxocarbenium ion facilitating equilibration under thermodynamic conditions, (iii) the greater magnitude of the anomeric effect in the 5-thiosugars, and (iv) decreased neighboring group participation in the per-esterified 5-thiosugars.

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

A crossed molecular beams and computational study on the unusual reactivity of banana bonds of cyclopropane (c-C 3 H 6 ; X 1 A$^{'}_{1}$ ) through insertion by ground state carbon atoms (C( 3 P j ))

The mechanism and chemical dynamics of the reaction of ground electronic state atomic carbon C( 3 P j ) with cyclopropane c-C 3 H 6 (X 1 A$^{'}_{1}$) have been explored by combining crossed molecular beams experiments with electronic structure calculations of the pertinent triplet C 4 H 6 potential energy surface and statistical computations of product branching ratios under single-collision conditions. The experimental findings suggest that the reaction proceeds via indirect scattering dynamics through triplet C 4 H 6 reaction intermediate(s) leading to C 4 H 5 product(s) plus atomic hydrogen via a tight exit transition state, with the overall reaction exoergicity evaluated as 231 ± 52 kJ mol –1 . The calculations indicate that C( 3 P j ) can easily insert into one of the three equivalent C–C ‘banana’ bonds of cyclopropane overcoming a low barrier of only 2 kJ mol –1 following the formation of a van der Waals reactant complex stabilized by 15 kJ mol –1 . The carbon atom insertion into one of the six C–H bonds is also feasible via a slightly higher barrier of 5 kJ mol –1 . Here, these results highlight an unusual reactivity of cyclopropane's banana C–C bonds, which behave more like unsaturated C–C bonds with a π-character than saturated σ C–C bonds, which are known to be generally unreactive toward the ground electronic state atomic carbon such as in ethane (C 2 H 6 ). The statistical theory predicts the overall product branching ratios at the experimental collision energy as 50% for 1-butyn-4-yl, 33% for 1,3-butadien-2-yl, i-C 4 H 5 , and 11% for 1,3-butadien-1-yl, n-C 4 H 5 , with i-C 4 H 5 (230 kJ mol –1 below the reactants) favored by the C–C insertion providing the best match with the experimentally observed reaction exoergicity. The C( 3 P j ) + c-C 3 H 6 reaction is predicted to be a source of C 4 H 5 radicals under the conditions where its low entrance barriers can be overcome, such as in planetary atmospheres or in circumstellar envelopes but not in cold molecular clouds. Both i- and n-C 4 H 5 can further react with acetylene eventually producing the first aromatic ring and hence, the reaction of the atomic carbon with c-C 3 H 6 can be considered as an initial step toward the formation of benzene.

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

A state-of-the-art review of experimental and computational studies of granular materials: Properties, advances, challenges, and future directions

Modeling of heterogeneous materials and media is a problem of fundamental importance to a wide class of phenomena and systems, ranging from condensed matter physics, soft materials, and composite media to porous media, biological systems, geosystems, ceramic engineering, pharmaceutical science and even in space discoveries. Among the most important materials are granular systems, which have received intense interest from the engineering, physics, and mathematics communities. In this review paper, the recent developments and new advances in experimental, and computational methods on a variety of scales and physics that extend understanding to a wide range of materials and phenomena are reviewed. Experimental advances include computed neutron and nanometer-scale tomography, magnetic resonance imaging, refractive index matching, digital image correlation, acoustic emission analysis, and the most recent 4D techniques. Furthermore, a tremendous shift has occurred from the continuum scale to micro-scale and developing multiscale approaches. As such, various computational methods, including, constitutive modeling, discrete modeling, and multiscale approaches, have been developed. In conclusion, aside from all these evolutions, more complicated modeling called coupled, or multiphysics, systems representing a simultaneous presence of heat, fluid, chemical variation, and mechanical effect are also explored.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Different bonding type along each crystallographic axis: Computational study of poly($\textit{p}$-phenylene terephthalamide)

Poly($\textit{p}$-phenylene terephthalamide) (PPTA) exhibits van der Waals (vdW) bonding along its $\textit{a}$ axis, hydrogen bonding along its $\textit{b}$ axis, and covalent bonding along its $\textit{c}$ axis. In this work, we explore the structural and mechanical properties of PPTA using density functional theory with various functionals including LDA, PBE, PBE+rVV10L, SCAN, and SCAN+rVV10, compared with available experiments. The hierarchy of nonempirical semilocal functionals LDA, PBE, and SCAN (not fitted to any multicenter bonded system) includes differing amounts of intermediate-range vdW interaction. rVV10 is the long-range vdW correction. (rVV10L differs from rVV10 only in the value of a range parameter.) Among the tested functionals, SCAN shows the best performance for the lattice parameters of PPTA along the two crystal directions involving vdW or hydrogen-bond interaction. The equilibrium lattice constants obtained by SCAN and PBE+rVV10L are closest to experimental data, while SCAN+rVV10 slightly overbinds the system. We study the mechanical response of PPTA by applying strain along three lattice directions. Due to the inclusion of vdW interaction, SCAN, PBE+rVV10L, and SCAN+rVV10 all exhibit correct bonding strain-energy curves. On the contrary, PBE strongly underestimates the vdW interaction needed to resist uniaxial stretching along the $\textit{a}$ axis. The Young's modulus and yield strength of PPTA are computed and compared with previous results. The experimental values are much smaller than the computed ones, mainly due to the fact that the PPTA fiber samples used for measurements are mechanically weaker than the perfect molecular crystals considered in the simulations. Interestingly, when a compressive uniaxial stress of 25 Kbar is applied along the $\textit{b}$ axis, a structural phase transition, in which the hydrogen bonds reform along one diagonal of the $\textit{ab}$ rectangle, is predicted by SCAN.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

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

Sequential Bond Dissociation Energies of Th + (CO) x , x = 3–6: Guided Ion Beam Collision-Induced Dissociation and Quantum Computational Studies

Collision-induced dissociation (CID) of [Th,xC,xO] + , x = 3–6, with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). Products are formed exclusively by the loss of CO ligands. Analyses of the kinetic energy-dependent CID product cross sections yield bond dissociation energies (BDEs) of (CO) x–1 Th + –CO at 0 K as 1.09 ± 0.05, 0.82 ± 0.07, 0.63 ± 0.05, and 0.70 ± 0.05 eV, respectively. Different structures of [Th,xC,xO] + were explored using various electronic structure methods, and BDEs for CO ligand loss from precursor [Th,xC,xO] + complexes were computed. Both experimental and theoretical results corroborate that the structures of [Th,xC,xO] + , x = 3–6, formed experimentally are homoleptic thorium cation carbonyl complexes, Th + (CO) x . The nonmonotonic trend in experimental BDEs is reproduced theoretically, although ambiguities in the spin states of the x = 4–6 complexes (doublet or quartet) remain. BDEs calculated at the coupled cluster with single, double, and perturbative triple excitations (CCSD(T))/cc-pVXZ//B3LYP/cc-PVXZ (X = T and Q) level and a complete basis set (CBS) extrapolation agree reasonably well with the experimental values for all complexes. Thorium oxide ketenylidene carbonyl cations, OTh + CCO(CO) y , y = 1–4, were calculated to be the most stable structures of [Th,xC,xO] + , x = 3–6, respectively; however, these are not observed in our experiment. Potential energy profiles (PEPs) having either quartet or doublet spin calculated at the B3LYP/cc-pVQZ level suggest that the failure to observe OTh + CCO(CO) y , y = 1–4, is the result of a barrier corresponding to the C–C bond formation, making the formation of OTh + CCO(CO) y inaccessible kinetically under the present experimental conditions.

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

Experimental and computational studies of load ratio effect on fatigue crack growth rates in Alloy 709 at elevated temperatures

Abstract Experimental testing and finite element simulations of fatigue crack growth were performed in austenitic 20Cr‐25Ni (Alloy 709) steel at different load ratios and elevated temperatures. The experimental tests were performed using compact tension specimens, and crack growth rates were measured at stress intensity factors ranging between 5 and 35 MPa√m. Fractographic analysis using scanning electron microscopy indicated crack surface roughness and secondary cracking depending on testing temperature and load ratio. Finite element simulations of fatigue crack growth were performed to compute plasticity‐induced crack opening loads and predict crack growth rates. Predictions of fatigue crack growth rates using finite element simulations were performed using the computed crack‐tip opening loads, and they are shown to match well the experimental measurements.

Engineering↗
DOE OSTI2024

Activation of H 2 O by ThO 2 – Experimental and Computational Studies

Here, a synergetic study that utilized anion photoelectron spectroscopy and high-level abinitio calculations has explored the activation of H 2 O molecules by ThO 2 – molecular anions. Both experiment and theory found conclusive evidence for said activation. In the experiments, this appeared as a tell-tale directional shift in the spectral profile of the anionic complex that ruled out physisorption,i.e., ThO 2 – (H 2 O), and implied chemisorption. In the computations, good agreement was found between the calculated and measured vertical detachment energies, and the atomic connectivity (the structure) of the resulting anionic complex was found to be [OTh(OH) 2 ] – .

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

Effect of isotropic and anisotropic porous microstructure on electrochemical performance of Li ion battery cathodes: An experimental and computational study

Liquid electrolyte mass transport is a major limitation affecting high-power Li ion batteries. Fast discharging causes Li salt depletion in the current collector region of the cathode which produces overpotential in the electrolyte and consequently a drop of cell voltage to below the cut-off voltage, especially at higher electrode thickness and discharge rate. In this study, through experiment and simulation, we have investigated the effect of electrode thickness, mass loading, discharge rate and tortuosity on electrolyte mass transport and final derived areal capacity and specific energy for electrodes having isotropic (normal tape casting) and anisotropic (freeze tape casting) porous microstructure. Furthermore, the macroporous channels in freeze tape cast electrodes facilitate Li salt transport and reduce the Li salt mass transport limitations even at high electrode thickness and discharge rates, and high electrode tortuosity. In our results, computer simulations show that freeze tape cast electrodes may be fully discharged up to 750 μm thickness at 1 C rate compared to 300 μm for normal tape cast electrodes with the same mass loading. Freeze tape cast electrodes also show stable maximum areal capacity for C rates about double the maximum C rates of their normal tape cast electrode counterparts with the same mass loading.

25 ENERGY STORAGE↗
DOE OSTI2023

State-of-the-art and review of condensation heat transfer for small modular reactor passive safety: Computational studies

The small modular reactor (SMR) is a promising option with added safety features, economical manufacturing, reliable parts, portability, and scalable energy capacity that emits no greenhouse gas during its operating lifespan. The SMR safety systems, however, need to be evaluated for design and licensing. Thus, they require the verification and validation of the physics models and correlations. This study focuses on state-of-the-art condensation heat transfer analysis and a review of previous studies related to the passive containment cooling system (PCCS) of a SMR. In the PCCS of a SMR, due to its smaller size containment, filmwise condensation is dominant and therefore emphasized in this study. Furthermore, previous condensation heat transfer studies for PCCSs did not make the SMR the primary focus, so a critical review for formulating the state-of-the-art is necessary. A previous review covered experimental condensation heat transfer studies with a brief overview of associated test facilities and empirical correlations. This review covers the empirical, resistance-layer and theoretical (numerical and commercial CFD) approaches.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Surface Functionalized Barium Titanate Nanoparticles: A Combined Experimental and Computational Study

Barium titanate (BTO) nanoparticles show great potential for use in electrostatic capacitors with high energy density. This includes both polymer composite and sintered capacitors. However, questions about the nanoparticles' size distribution, amount of agglomeration, and surface ligand effect on performance properties remain. Reducing particle agglomeration is a crucial step to understanding the properties of nanoscale particles, as agglomeration has significant effects on the composite dielectric constant. BTO surface functionalization using phosphonic acids is known reduce BTO nanoparticle agglomeration. We explore solution synthesized 10 nm BTO particles with tert-butylphosphonic acid ligands. Recent methods to quantifying agglomeration using an epoxy matrix before imaging shows that tert-butylphosphonic acid ligands reduce BTO agglomeration by 33%. Thermometric, spectroscopic, and computational methods provide confirmation of ligand binding and provide evidence of multiple ligand binding modes on the BTO particle surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational study of Li3BO3 and Li3BN2 I: Electrolyte properties of pure and doped crystals

Both Li3BO3 and Li3BN2 materials have promising properties for use in all-solid-state batteries and other technologies dependent on electrolytes with significant ionic conductivity. As the first of a two-part study, this paper reports the analysis of detailed simulations of Li ion diffusion in the monoclinic forms of these materials. Using both NEB and MD methods, it is clear that Li ion migration via vacancy mechanisms provides the most efficient ion transport in each material. While the results suggest that interstitial defects in these materials do not play a direct role in Li ion migration, their relative stability seems to enhance vacancy production via the formation of Frenkel-type defects. This may partially explain why the Li ion conductivities computed from MD simulations of samples initially containing a single Li ion vacancy are in reasonable agreement with measured values of this work for Li3BO3 and those reported in the literature for poorly crystalline samples of both materials. The possibility of increasing vacancy concentrations by substitutional doping (F for O in Li3BO3 and C for B in Li3BN2) is also examined, finding simulated conductivities comparable to those of the ideal vacancy model.

Li, Yan↗