Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion processes”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Peak radial growth of diffuse-porous species occurs during periods of lower water availability than for ring-porous and coniferous trees

Climate models project warmer summer temperatures will increase the frequency and heat severity of droughts in temperate forests of Eastern North America. Hotter droughts are increasingly documented to affect tree growth and forest dynamics, with critical impacts on tree mortality, carbon sequestration and timber provision. The growing acknowledgement of the dominant role of drought timing on tree vulnerability to water deficit raises the issue of our limited understanding of radial growth phenology for most temperate tree species. Here, we use well-replicated dendrometer band data sampled frequently during the growing season to assess the growth phenology of 610 trees from 15 temperate species over 6 years. Patterns of diameter growth follow a typical logistic shape, with growth rates reaching a maximum in June, and then decreasing until process termination. On average, we find that diffuse-porous species take 16–18 days less than other wood-structure types to put on 50% of their annual diameter growth. However, their peak growth rate occurs almost a full month later than ring-porous and conifer species (ca. 24 ± 4 days; mean ± 95% credible interval). Unlike other species, the growth phenology of diffuse-porous species in our dataset is highly correlated with their spring foliar phenology. We also find that the later window of growth in diffuse-porous species, coinciding with peak evapotranspiration and lower water availability, exposes them to a higher water deficit of 88 ± 19 mm (mean ± SE) during their peak growth than ring-porous and coniferous species (15 ± 35 mm and 30 ± 30 mm, respectively). Given the high climatic sensitivity of wood formation, our findings highlight the importance of wood porosity as one predictor of species climatic sensitivity to the projected intensification of the drought regime in the coming decades.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Water Sorption Study of Various Polymers and Polymeric Composites

In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗
DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2020

Diffuser Braze Development

At the request of SRTE, through the PDRD project, SRNL is developing a vacuum brazing process to facilitate the fabrication of hydrogen purification systems. In particular, this effort is directed at tube-style hydrogen diffusers. The diffusers are comprised of numerous palladium-silver tubes that are brazed to a nickel fitting, which is subsequently welded to a stainless steel tube, fitting, or tube sheet. Historically, torch brazing has been used as a heat source; for this work vacuum furnace brazing was suggested to improve the uniformity of the process. This report describes the braze process development for braze selection and simple nickel fitting to tube brazing. Five braze alloys were selected and sessile drop testing was conducted over a range of temperatures from 25 to 75° C over the liquidus temperature. Samples were characterized for wetting angle, erosion, and brazability. The cup-to-tube braze joints were inspected using digital radiography, ultrasonic testing, and metallographic examination. Two of the five alloys exhibited excellent fabrication properties and are considered for further development.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Dynamic response of a pulsed Burke-Schumann diffusion flame

Turbulent flames are often envisioned as an ensemble of random vortices interacting with the combustion process. A better understanding of the vortex-flame interactions therefore would be useful in improving the modeling of turbulent diffusion flames. Substantial simplification may be made by investigating controlled interactions in a laminar flame, as opposed to random interactions in a turbulent flame. The general goals of the research project are to improve our understanding of (1) the influence of buoyancy on co-flow diffusion flames and (2) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. As a first step toward objective (2), we conducted a joint experimental and numerical investigation of the vortex-flame interaction. Vortices were produced by mechanically pulsing the fuel flow at a low frequency, e.g., 10 Hz. Experiments were conducted using a nonflickering Burke-Schumann flame in both microgravity (mu-g) and normal gravity (1g) as a means of varying the buoyant force without modification of the pressure (i.e., density). The effects of buoyant convection may then be determined by a comparison of the mu-g and 1g results. The mu-g results may also reveal the important mechanisms which are masked or overwhelmed by buoyant convection in 1g. A numerical investigation was conducted using a validated, time-accurate numerical code to study the underlying physics during the flame interaction and to assist the interpretation of the experimental results.

Sheu, Jyh-Cherng↗
NASA NTRS1989

Soliton and strong Langmuir turbulence in solar flare processes

The occurrence of modulational instability in the current sheet of a solar flare is investigated. Special attention is given to the plasma microinstability in this sheet and its relation to the flare process. It is found that solitons or strong Langmuir turbulence are likely to occur in the diffusion region under solar flare conditions in which the electric resistivity could be enhanced by several orders of magnitude in the region, resulting in significant heating and stochastic acceleration of particles. A numerical example is used to demonstrate the transition of the magnetic field velocity and plasma density from the outer MHD region into the diffusive region and then back out again with the completion of the energy conversion process. This is all made possible by an increase in resistivity of four to five orders of magnitude over the classical value.

Song, M. T.↗
DOE OSTI2023

Experimental and computational design tools for industrial drying processes: Challenges in process‐limit prediction

Abstract The coating and drying of inks and slurries are important steps to manufacture a plethora of products. Drying processes, particularly, comprise energy‐intensive steps that affect product cost and quality. Prior work has highlighted failures of various multicomponent diffusivity models to conserve mass in dryer modeling and challenges in predicting process limits given variability in published values of key thermodynamic parameters. Herein, we develop a computational model and benchtop drying experiments to investigate these concerns for drying polymer‐laden coatings. Model predictions of process limits in a single‐zone drying oven demonstrate that published variability in Flory–Huggins parameter yields large variations in predicted operating temperatures above which blistering occurs. This indicates that caution should be exercised when choosing approaches to obtain or predict the Flory–Huggins parameter, and that both benchtop drying experiments and a set of additional experiments, such as sorption experiments, are needed to fully characterize and optimize a given drying process.

Parrish, Chance↗
NASA NTRS1999

Experimental Observations on a Low Strain Counter-Flow Diffusion Flame: Flow and Bouyancy Effects

Diffusion flames are of great interest in fire safety and many industrial processes. The counter-flow configuration provides a constant strain flow, and therefore is ideal to study the structure of diffusion flames. Most studies have concentrated on the high velocity, high strain limit, since buoyantly induced instabilities will disintegrate the planar flame as the velocity decreases. Only recently, experimental studies in microgravity conditions have begun to explore the low strain regimes. Numerical work has shown the coupling between gas phase reaction rates, soot reaction rates, and radiation. For these programs, size, geometry and experimental conditions have been chosen to keep the flame unaffected by the physical boundaries. When the physical boundaries can not be considered infinitely far from the reaction zone discrepancies arise. A computational study that includes boundary effects and accounts for the deviations occurring when the major potential flow assumptions are relaxed was presented by Borlik et al. This development properly incorporates all heat loss terms and shows the possibility of extinction in the low strain regime. A major constraint of studying the low strain regime is buoyancy. Buoyant instabilities have been shown to have a significant effect on the nature of reactants and heat transport, and can introduce instabilities on the flow that result in phenomena such as flickering or fingering. The counter-flow configuration has been shown to provide a flame with no symmetry disrupting instabilities for inlet velocities greater than 50 mm/s. As the velocity approaches this limit, the characteristic length of the experiment has to be reduced to a few millimetres so as to keep the Rayleigh number (Ra L = (β g 0 L 3 ∆T)/(𝑎v)) below 2000. In this work, a rectangular counter-flow burner was used to study a two-dimensional counter-flow diffusion flame. Flow visualisation and Particle Image Velocimetry served to describe the nature of the stagnation plane for strain rates smaller than 100 (1/s). These experiments were conducted with a non-reacting flow. Video images of a propane air diffusion flame were used to describe the behaviour of a diffusion flame in this regime. Flame geometry and pulsation frequency are described.

J. A. Sutula↗
DOE OSTI2016

Performance Evaluation of Spectroscopic Detectors for LEU Hold-up Measurements

The hold-up measurement of low-enriched uranium materials may require use of alternate detector types relative to the measurement of highly enriched uranium. This is in part due to the difference in process scale (i.e., the components are generally larger for low-enriched uranium systems), but also because the characteristic gamma-ray lines from 235 U used for assay of highly enriched uranium will be present at a much reduced intensity (on a per gram of uranium basis) at lower enrichments. Researchers at Oak Ridge National Laboratory examined the performance of several standard detector types, e.g., NaI(Tl), LaBr3(Ce), and HPGe, to select a suitable candidate for measuring and quantifying low-enriched uranium hold-up in process pipes and equipment at the Portsmouth gaseous diffusion plant. Detector characteristics, such as energy resolution (full width at half maximum) and net peak count rates at gamma ray energies spanning a range of 60–1332 keV, were measured for the above-mentioned detector types using the same sources and in the same geometry. Uranium enrichment standards (Certified Reference Material no. 969 and Certified Reference Material no. 146) were measured using each of the detector candidates in the same geometry. The net count rates recorded by each detector at 186 keV and 1,001 keV were plotted as a function of enrichment (atom percentage). Background measurements were made in unshielded and shielded configurations under both ambient and elevated conditions of 238 U activity. The highly enriched uranium hold-up measurement campaign at the Portsmouth plant was performed on process equipment that had been cleaned out. Therefore, in most cases, the thickness of the uranium deposits was less than the “infinite thickness” for the 186 keV gamma rays to be completely self-attenuated. Because of this, in addition to measuring the 186 keV gamma, the 1,001 keV gamma ray from 234m Pa—a daughter of 238 U in secular equilibrium with its parent—will also need to be measured. Based on the performance criteria of detection efficiency, energy resolution, peak-to-continuum ratios, minimum detectable limits, and the weight of the shielded probe, a shielded (0.5 in. thick lead shield) 2 × 2 in. NaI(Tl) detector is recommended for use. The recommended approach is to carry out analysis using data from both 186 keV and 1,001 keV gamma rays, and select a best result based on propagated uncertainty estimates. It is also highly recommended that a two-point gain stabilization scheme based on an 241 Am seed embedded in the probe be implemented. Shielding configurations to reduce the impact of background interference on the measurement of 1,001 keV gamma-ray are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Stochastic modeling of catalytic processes in nanoporous materials: Beyond mean-field approach

Transport and reaction in zeolites and other porous materials, such as mesoporous silica particles, has been a focus of interest in recent years. This is in part due to the possibility of anomalous transport effects (e.g. single-file diffusion) and its impact in the reaction yield in catalytic processes. Computational simulations are often used to study these complex nonequilibrium systems. Computer simulations using Molecular Dynamics (MD) techniques are prohibitive, so instead coarse grained one-dimensional models with the aid of Kinetic Monte Carlo (KMC) simulations are used. Both techniques can be computationally expensive, both time and resource wise. These coarse-grained systems can be exactly described by a set of coupled stochastic master equations, that describe the reaction-diffusion kinetics of the system. The equations can be written exactly, however, coupling between the equations and terms within the equations make it impossible to solve them exactly; approximations must be made. One of the most common methods to obtain approximate solutions is to use Mean Field (MF) theory. MF treatments yield reasonable results at high ratios of reaction rate k to hop rate h of the particles, but fail completely at low k=h due to the over-estimation of fluxes of particles within the pore. We develop a method to estimate fluxes and intrapore diffusivity in simple one- dimensional reaction-diffusion models at high and low k=h, where the pores are coupled to an equilibrated three-dimensional fluid. We thus successfully describe analytically these simple reaction-diffusion one-dimensional systems. Extensions to models considering behavior with long range steric interactions and wider pores require determination of multiple boundary conditions. We give a prescription to estimate the required parameters for these simulations. For one dimensional systems, if single-file diffusion is relaxed, additional parameters to describe particle exchange have to be introduced. We use Langevin Molecular Dynamics (MD) simulations to assess these parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Nanoscale insight into the p-n junction of alkali-incorporated Cu(In,Ga)Se 2 solar cells

The effects of alkali diffusion and post-deposition treatment in three-stage processed Cu(In,Ga)Se 2 solar cells are examined by using atom probe tomography and electrical property measurements. Cells, for which the substrate was treated at 650 °C to induce alkali diffusion from the substrate prior to absorber deposition, exhibited high open-circuit voltage (758 mV) and efficiency (18.2%) and also exhibited a 50 to 100-nm-thick ordered vacancy compound layer at the metallurgical junction. Surprisingly, these high-temperature samples exhibited higher concentrations of K at the junction (1.8 at.%) than post-deposition treatment samples (0.4 at.%). A model that uses Ga/(Ga + In) and Cu/(Ga + In) profiles to predict bandgaps (+/-17.9 meV) of 22 Cu(In,Ga)Se2 solar cells reported in literature was discussed and ultimately used to predict band properties at the nanoscale by using atom probe tomography data. The high-temperature samples exhibited a greater drop in the valence band maximum (200 meV) due to a lower Cu/(Ga + In) ratio than the post-deposition treatment samples. There was an anticorrelation of K concentrations and Cu/(Ga + In) ratios for all samples, regardless of processing conditions. In conclusion, changes in elemental profiles at the active junctions correlate well with the electrical behaviour of these devices.

14 SOLAR ENERGY↗
DOE OSTI2024

Spatiotemporal Reaction Dynamics Control in Two‐Photon Polymerization for Enhancing Writing Characteristics

Since 2001, 3D microfabrication based on two‐photon polymerization (TPP) has drawn extensive attention and interest in biology, optics, photonics, material science, and high‐energy physics. The in‐volume fabrication capability due to the threshold behavior of two‐photon absorption enables TPP higher flexibility compared with other nanofabrication techniques. However, as determined by the in‐volume fabrication feature as well as various reaction dynamics, the writing characteristics of TPP, such as throughput, accuracy, surface quality, and fabrication capability, are still limited. Herein, a comprehensive study is performed on the spatiotemporal behavior of reaction dynamics during TPP fabrication, mainly focusing on spatiotemporal characteristics of radical diffusion, photothermal effect, microscale mechanics, and voxel stacking process. Based on the study, a nonsequential fabrication method is established to simultaneously improve key writing characteristics of TPP and realize sharp features, high speeds, large overhang structure, and smooth surfaces. The method established in this work can be applied to improve the performance of functional devices for various fields.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Cyclic Oxidation Testing and Modelling: A NASA Lewis Perspective

The Materials Division of the NASA Lewis Research Center has been heavily involved in the cyclic oxidation of high temperature materials for 30 years. Cyclic furnace and burner rig apparati have been developed, refined, and replicated to provide a large scale facility capable of evaluating many materials by a standard technique. Material behavior is characterized by weight change data obtained throughout the test, which has been modelled in a step-wise process of scale growth and spallation. This model and a coupled diffusion model have successfully described cyclic behavior for a number of systems and have provided insights regarding life prediction and variations in the spalling process. Performance ranking and mechanistic studies are discussed primarily for superalloys and coating alloys. Similar cyclic oxidation studies have been performed on steels, intermetallic compounds, thermal barrier coatings, ceramics, and ceramic composites. The most common oxidation test was performed in air at temperatures ranging from 800 deg. to 1600 C, for times up to 10000 h, and for cycle durations of 0.1 to 1000 h. Less controlled, but important, test parameters are the cooling temperature and humidity level. Heating and cooling rates are not likely to affect scale spallation. Broad experience has usually allowed for considerable focus and simplification of these test parameters, while still revealing the principal aspects of material behavior and performance. Extensive testing has been performed to statistically model the compositional effects of experimental alloys and to construct a comprehensive database of complex commercial alloys.

Smialek, J. L.↗
DOE OSTI2025

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗
DOE OSTI2020

A phase field study of the thermal migration of gas bubbles in UO 2 nuclear fuel under temperature gradient

Phase field models are developed to study the gas bubble migration in uranium dioxide nuclear fuel in which a large temperature gradient exists during the operation. In this work, thermal diffusion mechanism for nanosized gas bubbles and vapor transport process for micron-sized gas bubbles are considered, respectively. In both cases, gas bubbles migrate to the high-temperature area. Due to the velocity difference between leading and trailing edges of the gas bubbles, nanosized gas bubbles are elongated along the temperature gradient direction when thermal diffusion is dominated. Micron-sized gas bubbles are either compressed along temperature gradient direction to form lenticular shape bubbles or elongated along temperature gradient direction, depending on the location of the gas bubbles within the fuel pellet. Initial gas bubble radius has no significant effect on the gas bubble migration velocity for both thermal diffusion and vapor transport mechanisms. We notice that the shape change of the gas bubble due to vapor transport mechanism has no significant effect on the migration velocity. Furthermore, the center cavity formation is also captured by our model which is due to the migration and accumulation of lenticular gas bubbles at the center of the fuel pellet. The modeling results compare well with experimental observations and theoretical analysis in the literature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Electrochemical properties and lithium ion diffusion in Li{sub 4}FeSbO{sub 6} studied by first principle

Due to the high capacity, Li-rich materials Li{sub 2}MO{sub 3} (M = transition metal) have attracted considerable attention as the next generation of Li-ion batteries. Li{sub 4}FeSbO{sub 6} is a new Li-rich layered oxide material with antiferromagnet honeycomb structure. In this work, the electrochemical behavior, charging process and oxygen stability of Li{sub x}FeSbO{sub 6} (0 ≤ x<4) material are investigated by using the spin-polarized density functional theory. Moreover, climbing image nudged elastic band (CI-NEB) method is used to calculate the Li ion migration properties. The results show that Li{sub 4}FeSbO{sub 6} exhibits poor inherent conductivity with the band gap of 2.45 eV. During charging process, the structure presents the excellent volume stability. When Li is extracted, the charge compensation is mainly contributed by the oxygen atoms through analyzing the Bader charges of each element. In addition, oxygen evolution reactions will occur in Li{sub x}FeSbO{sub 6} (x ≤ 1.5), which will decay the capacities during cycling process. Finally, we calculated that the lithium ion can diffuse in a three-dimensional pathway with the activation barriers from 0.36 eV to 0.67 eV. - Graphical abstract: Li{sub 4}FeSbO{sub 6} is a new Li-rich layered oxide material with antiferromagnet honeycomb structure. We calculated that the lithium ion can diffuse in a three-dimensional pathway Display Omitted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗