Search NASA⌕ Search

SEARCH · Search NASA

Results for “structural disorder”

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 487 records · Page 27

Long-Range Ordered Amorphous Atomic Chains as Building Blocks of a Superconducting Quasi-One-Dimensional Crystal

Crystalline and amorphous structures are two of the most common solid-state phases. Crystals having orientational and periodic translation symmetries are usually both short-range and long-range ordered, while amorphous materials have no long-range order. Short-range ordered but long-range disordered materials are generally categorized into amorphous phases. In contrast to the extensively studied crystalline and amorphous phases, the combination of short-range disordered and long-range ordered structures at the atomic level is extremely rare and so far has only been reported for solvated fullerenes under compression. Here, a report on the creation and investigation of a superconducting quasi-1D material with long-range ordered amorphous building blocks is presented. Using a diamond anvil cell, monocrystalline (TaSe4)(2)I is compressed and a system is created where the TaSe(4)atomic chains are in amorphous state without breaking the orientational and periodic translation symmetries of the chain lattice. Strikingly, along with the amorphization of the atomic chains, the insulating (TaSe4)(2)I becomes a superconductor. The data provide critical insight into a new phase of solid-state materials. The findings demonstrate a first ever case where superconductivity is hosted by a lattice with periodic but amorphous constituent atomic chains.

combination of crystalline and amorphous structure↗

Vacancy-mediated complex phase selection in high entropy alloys

Phase selection in Ti-Zr-Hf-Al high-entropy alloys was investigated by in-situ high-energy x-ray diffraction, single-crystal x-ray diffraction, and KKR-CPA electronic-structure methods that address disorder and vacancies, predicting formation enthalpy and chemical short-range order (SRO). Samples with varying Al content were synthesized, characterized, and computationally assessed to ascertain the composition-dependent phase selection, as increased Al content often acts as a stabilizer of a body-centered-cubic structure. Equiatomic TiZrHfAl was especially interesting due to its observed bcc superstructure – a variant of γ-brass with 4 vacancies per cell (not 2 as in γ-brass). We highlight how vacancy ordering mediates selection of this variant of γ-brass, which is driven by vacancy-atom SRO that dramatically suppress all atomic SRO. As vacancies are inherent in processing refractory systems, we expect that similar discoveries await in other high entropy alloys or in revisiting older experimental data.

36 MATERIALS SCIENCE↗

Effect of solvothermal synthesis parameters on the crystallite size and atomic structure of cobalt iron oxide nanoparticles

We here investigate how the synthesis method affects the crystallite size and atomic structure of cobalt iron oxide nanoparticles. By using a simple solvothermal method, we first synthesized cobalt ferrite nanoparticles of ca. 2 and 7 nm, characterized by Transmission Electron Microscopy (TEM), Small Angle X-ray scattering (SAXS), X-ray and neutron total scattering. The smallest particle size corresponds to only a few spinel unit cells. Nevertheless, Pair Distribution Function (PDF) analysis of X-ray and neutron total scattering data shows that the atomic structure, even in the smallest nanoparticles, is well described by the spinel structure, although with significant disorder and a contraction of the unit cell parameter. These effects can be explained by the surface oxidation of the small nanoparticles, which is confirmed by X-ray near edge absorption spectroscopy (XANES). Neutron total scattering data and PDF analysis reveal a higher degree of inversion in the spinel structure of the smallest nanoparticles. Neutron total scattering data also allow magnetic PDF (mPDF) analysis, which shows that the ferrimagnetic domains correspond to ca. 80% of the crystallite size in the larger particles. A similar but less well-defined magnetic ordering was observed for the smallest nanoparticles. Finally, we used a co-precipitation synthesis method at room temperature to synthesize ferrite nanoparticles similar in size to the smallest crystallites synthesized by the solvothermal method. Structural analysis with PDF demonstrates that the ferrite nanoparticles synthesized via this method exhibit a significantly more defective structure compared to those synthesized via a solvothermal method.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

High electron mobility of Al x Ga1− x N evaluated by unfolding the DFT band structure

We calculate the alloy-disorder-limited electron mobility of AlxGa1−xN from first principles. AlxGa1−xN is a technologically important ultra-wide-bandgap alloy with promise in light emitting diodes and high-power transistors. Alloying introduces statistical disorder, which causes electrons to scatter between different crystal-momentum states, leading to a reduction in mobility for intermediate alloy compositions. The corresponding lifetime, which appears as an energy broadening in the band structure, can be evaluated by unfolding the band structure from the supercell basis to the primitive-cell basis. We fit the first-principles band broadening with a model scattering potential and evaluate the low-field electron mobility using the semiclassical Boltzmann transport equation in the relaxation-time approximation. Our calculated mobility is in agreement with experimental values. We also find the lowest alloy-scattering electron mobility (total electron mobility) across the entire composition range to be 186 cm2/V s (136 cm2/V s), which is comparable to the highest electron mobility predicted in the competitor system, β-(AlxGa1−x)2O3. Our results elucidate the intrinsic limits imposed by alloy disorder on electron transport in AlxGa1−xN.

Pant, Nick (ORCID:0000000254506045)↗

Reduced Graphene Oxide Aerogels with Functionalization-Mediated Disordered Stacking for Sodium-Ion Batteries

Surface modified reduced graphene oxide (rGO) aerogels were synthesized using the hydrothermal method. Ethylene diamine (EDA) and α-cyclodextrin (CD) were used to functionalize the surface of the graphene oxide layers. The oxygen reduction and surface modification occurred in-situ during the hydrothermal self-assembly process. The chemical functionality and structure of the resulting ethylene diamine modified (rGO-EDA) and cyclodextrin modified (rGO-CD) aerogels as well as of the pristine unmodified rGO aerogel were studied using XPS, SEM, XRD, and SANS techniques. The overall surface composition showed a significant decrease in the oxygen content for all synthesized aerogels. The surface modified aerogels were characterized by a disordered stacking of the assembled rGO layers. The surface functionalities resulted in a broad distribution of the interlayer spacing and introduced structural heterogeneities. Such disordered structures can enable a better adsorption mechanism of the sodium ions. Coin cells based on the synthesized aerogels and sodium metal were assembled and tested at several charge and discharge rates. The correlation between the surface functionality of the rGO, the induced structural heterogeneities due to the disordered stacking, and the electrochemical performance of sodium-ion batteries were investigated. Operando XRD measurements were carried out during the battery cycling to investigate the adsorption or intercalation nature of the sodiation mechanism.

25 ENERGY STORAGE↗

Multi-scale investigation of heterogeneous swift heavy ion tracks in stannate pyrochlore

Er 2 Sn 2 O 7 pyrochlore was irradiated with swift heavy Au ions (2.2 GeV), and the induced structural modifications were systematically examined using complementary characterization techniques including transmission electron microscopy (TEM), X-ray diffraction (XRD), and neutron total scattering with pair distribution function (PDF) analysis. Each technique probes different aspects and length scales of the transformed material regions. TEM revealed a core–shell ion track structure—an amorphous core surrounded by a disordered, anion-deficient fluorite shell—which was confirmed by XRD. Neutron total scattering, with sensitivity to the oxygen sublattice, provided relative fractions of amorphous and disordered fluorite phases and confirmed the presence of a defective pyrochlore phase, which largely maintains its structural ordering but is clearly distinct from the pristine pyrochlore matrix. This defect-rich pyrochlore phase forms a halo extending radially beyond the well-characterized core–shell track morphology observed in electron micrographs. Despite their differing long-range periodicity, the short-range structures of the amorphous, disordered, and defective pyrochlore phases are all modeled well with a weberite-type configuration. Evolution of the phase fractions with increasing ion fluence was examined to ascertain the phase-to-phase pathways that occur during primary and secondary ion impact. Furthermore, this approach extends knowledge about the multi-scale response of stannate pyrochlores to swift heavy ion irradiation in the electronic energy loss regime and improves existing track-overlap models.

36 MATERIALS SCIENCE↗

Quasiclassical Theory of C4-Symmetric Magnetic Order in Disordered Multiband Metals

Recent experimental studies performed in the normal state of iron-based superconductors have discovered the existence of the C 4 -symmetric (tetragonal) itinerant magnetic state. This state can be described as a spin density wave with two distinct magnetic vectors Q → 1 and Q → 2 . Given an itinerant nature of magnetism in iron-pnictides, we develop a quasiclassical theory of tetragonal magnetic order in disordered three-band metal with anisotropic band structure. Within our model we find that the C 4 -symmetric magnetism competes with the C 2 -symmetric state with a single Q → magnetic structure vector. Our main results is that disorder promotes tetragonal magnetic state which is in agreement with earlier theoretical studies.

74.45.c, 74.50. r, 74.20.R↗

MS25: Materials Science-Focused Benchmark Data Set for Machine Learning Interatomic Potentials

Here, we present MS25, a benchmark data set for evaluating machine learning interatomic potentials (MLIPs) across diverse materials-relevant systems including MgO surfaces, liquid water, zeolites, a catalytic Pt surface reaction, high-entropy alloys (HEAs), and disordered Zr-oxides. Five MLIP architectures (MACE, NequIP, Allegro, MTP, and Torch-ANI) are trained and tested, focusing not only on traditional metrics (energies, forces, and stresses) but also explicitly validating derived physical observables such as lattice constants, volumes, and reaction barriers. We find that most models reach comparable accuracy on standard error metrics across the simple systems, although equivariant MLIPs offer 1.5–2× improvements over nonequivariant MLIPs in energy and force error for structurally complex or compositionally disordered environments such as HEAs and Zr–O systems. Our analysis highlights that low errors in energy and force predictions do not guarantee reliable observables, emphasizing the necessity of explicit validation. We demonstrate limitations in cross-framework transferability, as models trained on one zeolite framework (CHA) fail to reliably generalize to predictions of structurally distinct frameworks (e.g., MFI). Size-extensive tests show some dependence on system size for MgO, resulting from forced periodicity. The HEA and Zr–O data sets are identified as challenging tests for future benchmarks and MLIP model architecture developments as they show significant differentiation in error between MLIP architectures and are still relatively difficult at 1000 training images. Moving forward, we recommend that benchmarking efforts shift their focus from marginal accuracy improvements in energy and force errors toward identifying and understanding model failure modes, rigorously assessing transferability, and evaluating how their errors affect observable predictions. For researchers looking to choose an MLIP architecture, we suggest selecting equivariant MLIP architectures if the complexity of the system is a challenge. For simple materials problems, auxiliary features such as integration with molecular dynamics engines, trade-offs between computational data set generation cost vs MLIP inference speed, and framework integration may play a more important decision factor than small differences in error metrics that are unlikely to matter for production-level research.

chemical structure↗

Enhancing the Electrocatalysis of LiNi 0.5 Co 0.2 Mn 0.3 O 2 by Introducing Lithium Deficiency for Oxygen Evolution Reaction

LiNi 0.5 Co 0.2 Mn 0.3 O 2 (NCM523), as a cathode material for rechargeable lithium-ion batteries, has attracted considerable attention and been successfully commercialized for decades. NCM is also a promising electrocatalyst for the oxygen evolution reaction (OER), and the catalytic activity is highly correlated to its structure. In this paper, we successfully obtain NCM523 with three different structures: spinel NCM synthesized at low temperature (LT-NCM), disordered NCM (DO-NCM) with lithium deficiency obtained at high temperature, and layered hexagonal NCM at high temperature (HT-NCM). By introducing lithium deficiency to tune the valence state of transition metals in NCM from Ni 2+ to Ni 3+ , DO-NCM exhibits the best catalytic activity with the lowest onset potential (~1.48 V) and Tafel slope (~85.6 mV dec -1 ), whereas HT-NCM exhibits the worst catalytic activity with the highest onset potential (~1.63 V) and Tafel slope (~241.8 mV dec -1 ).

36 MATERIALS SCIENCE↗

Layered polymorph of titanium triiodide

A previously unreported layered spin ½ triangular lattice polymorph of TiI 3 is described, synthesized under 6 GPa of applied pressure at 900 °C, but stable at atmospheric pressure. This air-sensitive material has a CdI 2 -type layered structure (P-3m1 (no. 164), a = 4.012 Å and c = 6.641 Å at 120 K, Z = 1 of Ti 0.667 I 2 ) with an in-plane triangular lattice, related to that of TiI 4 (Ti 0.5 I 2 ). Although the TiI 3 formula is consistent with expectations for a layered honeycomb lattice of spin ½ Ti(III), there appears to be disorder in the crystal structure. Magnetic susceptibility and heat capacity measurements suggest that the material undergoes several low temperature phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cation-Induced Disruption of the Local Structure of Water in Faujasite

Here, the changes in the structure of water confined in faujasite were resolved by a series of ab initio molecular dynamic simulations of explicit water molecules in faujasite with varying Si/Al ratio and counterions present. Radial distribution functions of water oxygen atom pairs, as well as counterions to Al pairs, were calculated to quantify the structure of water in faujasite and the distribution of counterions in the water-filled faujasite pores. Radial distribution functions reveal that confinement without counterions results in more water molecules hydrogen bound to each other. In hydrogen-exchanged faujasite, the increase in aluminum content led to the disruption of the confined water structure as charge-compensating protons dissociated from Al tetrahedra and entered the liquid water in the pores. In addition, a potential of mean force of proton transfer shows that protons are thermodynamically favored to be in the liquid water rather than bound on Al tetrahedra. Exchanged Na + and La 3+ counterions, which remain closer to Al tetrahedra in liquid water than H + , also disordered the confined water structure. Two La 3+ ions disrupted water less than did six Na + ions, suggesting that, for the alkali and lanthanide ions, the number of counterions in FAU has a larger effect on the structure of water than the charge of the ion itself.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-assembly of co-continuous nanostructured copolymer templates with compositional and architectural dispersity

In our previous DOE-supported project (DE-SC0016208), we demonstrated that a particular architecture dubbed ‘randomly end-linked copolymer networks’ (RECNs) led to especially robust self-assembly of disordered cocontinuous nanostructures. Such structures, wherein both microphases are fully percolating in three dimensions, are attractive for a wide-variety of energy-relevant applications including purification membranes, separators in batteries and fuel cells, catalyst supports, and high surface area electrodes. In the current project we sought to substantially extend our fundamental understanding of how different types of randomness in blocky copolymer architecture contribute to the robust formation of disordered cocontinuous nanostructures. To do so, we filled in the sizable gap between the previously studied cases of randomly linked linear multiblock copolymers and RECNs by studying two intermediate cases: (i) random miktoarm star polymers, which contain dispersity in the number of the two chemically-incompatible strands located at each junction and therefore the preferred interfacial curvature of each polymer, and (ii) randomly branched copolymers, which also possess bridging strands that stretch across domains and apply local ‘elastic’ forces to the interfaces. In addition, we studied a different network architecture, statistically crosslinked copolymer networks, that is more straightforward to prepare using chain growth polymerization methods, enabling more facile formation of cocontinuous nanostructures. Finally, we extended our approaches from model systems based on polystyrene and poly(lactic acid) to copolymer architectures containing either (i) polysulfone, an engineering polymer that provides excellent mechanical properties for membrane applications, or (ii) polyacrylonitrile, providing a convenient route to nanoporous carbon materials.

36 MATERIALS SCIENCE↗

Pair distribution function analysis of nano-object assemblies

The atomic pair distribution function (aPDF) analysis technique, also known as the total scattering method, which considers both Bragg and diffuse scattering, has been used extensively to probe local atomic arrangements in crystalline and disordered materials. In contrast, there have been limited applications of the PDF in self-assembled nanomaterials, which represent a class of materials built from nanoscale objects, such as nano-colloids, micelles and proteins. As distinguished from atoms, nano-objects have polydispersity in size and shape, and such form-factor effects complicate the application of PDF analysis to nano-systems. Here, in this paper, the application of the PDF is extended to spherical nano-object assemblies and the formulae for the nano-PDF (nPDF) are derived, showing some differences from the aPDF. By numerical simulations, the properties of the nPDF (peak broadening and pattern profile) are studied systematically as a function of structural features, such as nano-object parameters (size and size polydispersity) and assembly structural features (size, shape, structure type and lattice disorder), and of data processing parameters ( q cut-off and `missing' data in ultra-small-angle regions). The nPDF analysis method is found to provide an effective route to revealing not only nanoscale but also mesoscale structural properties, for example the morphology of a nano-assembly.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Lone-Electron-Pair Micelles Strengthen Bond Anharmonicity in MnPb16Sb14S38 Complex Sulfosalt Leading to Ultra-Low Thermal Conductivity

Designing crystalline solids in which intrinsic extremely low lattice thermal conductivity mainly arises from its unique bonding nature rather than structure complexity and/or atomic disorder could promote thermal energy manipulation and utilization for applications ranging from thermoelectric energy conversion to thermal barrier coatings. Here, we report an extremely low lattice thermal conductivity of ~0.34 W m-1 K-1 at 300 K in the new complex sulfosalt MnPb16Sb14S38. We attribute the ultra-low lattice thermal conductivity to a synergistic combination of scattering mechanisms involving (1) strong bond anharmonicity in various structural building units, owing to the presence of stereoactive lone-electron-pair (LEP) micelles, and (2) phonon scattering at the interfaces between building units of increasing size and complexity. Remarkably, low-temperature heat capacity measurement revealed a Cp value of 0.206 J g-1 K-1 at T > 300 K, which is 22% lower than the Dulong-Petit value (0.274 J g-1 K-1). Further analysis of the Cp data and sound velocity (ν = 1834 m/s) measurement yielded Debye temperature values of 161 K and 187 K, respectively. Here, the resulting Grüneisen parameter, γ = 1.65, further supports strong bond anharmonicity as the dominant mechanism responsible for the observed extremely low lattice thermal conductivity.

36 MATERIALS SCIENCE↗

Epitaxial intergrowths and local oxide relaxations in natural bixbyite Fe 2– x Mn x O 3

The scattering pattern of a crystal obeys the symmetry of the crystal structure through the corresponding Laue group. This is usually also true for the diffuse scattering, containing information about disorder, but here a case is reported where the diffuse scattering is of lower symmetry than the parent crystal structure. The mineral bixbyite has been studied by X-ray and neutron scattering techniques since 1928 with some of the most recent studies characterizing the low-temperature transition to a magnetically disordered spin-glass state. However, bixbyite also exhibits structural disorder, and here single-crystal X-ray and neutron scattering is used to characterize the different modes of disorder present. One-dimensional rods of diffuse scattering are observed in the cubic mineral bixbyite, which break the expected symmetry of the scattering pattern. It is shown that this scattering arises from epitaxial intergrowths of the related mineral, braunite. The presence of this disorder mode is found to be directly observable as well-defined residuals in the average structure refined against the Bragg diffraction. An additional three-dimensional diffuse scattering component is observed in neutron scattering data, which is shown to originate from the substitutional disorder on the Fe/Mn sites. This occupational disorder gives rise to local relaxations of the oxide sublattice, and the pattern of oxide displacements can be rationalized based on crystal-field theory. The combined use of neutron and X-ray single-crystal scattering techniques highlights their great complementarity. In particular, the large sample requirements for neutron scattering experiments prove to be an obstacle in solving the intergrowth disorder due to several growth orientations, whereas for X-ray scattering the one-dimensional nature of the intergrowth disorder renders solving this a more tractable task. On the other hand, the oxide relaxations cannot be resolved using X-rays due to the low Mn/Fe contrast. By combining the two approaches both types of disorder have been characterized.

36 MATERIALS SCIENCE↗

Recent Results from Epitaxial Growth on Step Free 4H-SiC Mesas

This paper updates recent progress made in growth, characterization, and understanding of high quality homoepitaxial and heteroepitaxial films grown on step-free 4H-SiC mesas. First, we report initial achievement of step-free 4H-SiC surfaces with carbon-face surface polarity. Next, we will describe further observations of how step-free 4H-SiC thin lateral cantilever evolution is significantly impacted by crystal faceting behavior that imposes non-uniform film thickness on cantilever undersides. Finally, recent investigations of in-plane lattice constant mismatch strain relief mechanisms observed for heteroepitaxial growth of 3C-SiC as well as 2H-AlN/GaN heterofilms on step-free 4H-SiC mesas will be reviewed. In both cases, the complete elimination of atomic heterointerface steps on the mesa structure enables uniquely well-ordered misfit dislocation arrays to form near the heterointerfaces with remarkable lack of dislocations threading vertically into the heteroepilayers. In the case of 3C-SiC heterofilms, it has been proposed that dislocation half-loops nucleate at mesa edges and glide laterally along the step-free 3C/4H interfaces. In contrast, 3C-SiC and 2H-AlN/GaN heterofilms grown on 4H-SiC mesas with steps exhibit highly disordered interface misfit dislocation structure coupled with 100X greater density of dislocations threading through the thickness of the heteroepilayers. These results indicate that the presence of steps at the heteroepitaxial interface (i.e., on the initial heteroepitaxial nucleation surface) plays a highly important role in the defect structure, quality, and relaxation mechanisms of single-crystal heteroepitaxial films.

Neudeck, Philip G.↗

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

The transformation of diamond to graphite: Experiments reveal the presence of an intermediate linear carbon phase

Natural diamonds that have been partially replaced by graphite have been observed to occur in natural rocks. While the graphite-to-diamond phase transition has been extensively studied the opposite of this (diamond to graphite) remains poorly understood. We performed high-pressure and temperature hydrous and anhydrous experiments up to 1.0 GPa and 1300 °C using Amplex premium virgin synthetic diamonds (20–40 μm size) as the starting material mixed with Mg(OH) 2 as a source of H 2 O for the hydrous experiments. The experiments revealed that the diamond-to-graphite transformation at P = 1.0 GPa and T = 1300 °C was triggered by the presence of H 2 O and was accomplished through a three-stage process. Stage 1: diamond reacts with a supercritical H 2 O producing an intermediate 200–500 nm size “globular carbon” phase. This phase is a linear carbon chain; i.e. a polyyne or carbyne. Stage 2: the linear carbon chains are unstable and highly reactive, and they decompose by zigzagging and cross-linking to form sp 2 -hybridized structures. Stage 3: normal, disordered, and onion-like graphite is produced by the decomposition of the sp-hybridized carbon chains which are re-organized into sp 2 bonds. Our experiments show that there is no direct transformation from sp 3 C-bonds into sp 2 C-bonds. Our hydrous high-pressure and high-temperature experiments show that the diamond-to-graphite transformation requires an intermediate metastable phase of a linear hydrocarbon. This process also provides a simple mechanism for the substitution of other elements into the graphite structure (e.g. H, S, O).

36 MATERIALS SCIENCE↗