Structure of crystalline and amorphous materials in the NASICON system Na

The structure of crystalline and amorphous materials in the sodium (Na) super-ionic conductor system Na 1 + x Al x Ge 2 − x (PO 4 ) 3 with x = 0, 0.4, and 0.8 was investigated by combining (i) neutron and x-ray powder diffraction and pair-distribution function analysis with (ii) 27 Al and 31 P magic angle spinning (MAS) and 31 P/ 23 Na double-resonance nuclear magnetic resonance (NMR) spectroscopy. A Rietveld analysis of the powder diffraction patterns shows that the x = 0 and x = 0.4 compositions crystallize into space group-type R¯3, whereas the x = 0.8 composition crystallizes into space group-type R¯3c. For the as-prepared glass, the pair-distribution functions and 27 Al MAS NMR spectra show the formation of sub-octahedral Ge and Al centered units, which leads to the creation of non-bridging oxygen (NBO) atoms. The influence of these atoms on the ion mobility is discussed. When the as-prepared glass is relaxed by thermal annealing, there is an increase in the Ge and Al coordination numbers that leads to a decrease in the fraction of NBO atoms. A model is proposed for the x = 0 glass in which super-structural units containing octahedral Ge (6) and tetrahedral P (3) motifs are embedded in a matrix of tetrahedral Ge (4) units, where superscripts denote the number of bridging oxygen atoms. The super-structural units can grow in size by a reaction in which NBO atoms on the P (3) motifs are used to convert Ge (4) to Ge (6) units. The resultant P (4) motifs thereby provide the nucleation sites for crystal growth via a homogeneous nucleation mechanism.


I. INTRODUCTION
Homogeneous (or internal) nucleation in bulk glass can be used to create a uniform distribution of crystallites of controllable size and shape throughout a material.This process is crucial in the production of many glass-ceramics that include photosensitive and photochromic materials, along with tough fracture-resistant ceramics that have uses ranging from cookware to missile nose-cones. 1 The creation of a uniform distribution of crystallites also makes homogeneous nucleation an attractive option for optimizing the properties of superionic conducting materials for battery and sensor applications.A case example is provided by the technologically relevant Na 1+x AlxGe 2−x (PO 4 ) 3 (0 ≤ x ≤ 1) (NAGP) system in which the structure of the crystalline phase is based on that of the sodium (Na) super-ionic conductor (NASICON) NaGe 2 (PO 4 ) 3 . 2nitial work on these highly conducting materials, which are candidates for the electrolyte in all-solid-state sodium ion The Journal of Chemical Physics ARTICLE scitation.org/journal/jcpbatteries, 3,4 focused on pressed powders.However, it was discovered that they can also be prepared via the glass-ceramic route, giving the significant advantage of moldable bulk materials. 5It is therefore desirable to know the glass structure and how this structure changes during thermal annealing and the early stages of crystal nucleation and growth.The latter is of general scientific interest for any glass undergoing homogeneous nucleation but has not been explored in non-siliceous systems. 6n this paper, we investigate the structure of NAGP materials with x = 0, 0.4, and 0.8.Here, the structure of the crystalline phase with x = 0 is based on a negatively charged three-dimensional framework in which tetrahedral PO 4 units are linked via bridging oxygen (BO) atoms to four octahedral GeO 6 units.This framework of corner-sharing polyhedra contains interstitial sites in which the Na + ions reside 2 and where the size of a connection between sites determines the extent to which it acts as a bottleneck for ion motion. 7The ionic conductivity increases when a proportion of the Ge 4+ ions are aliovalently substituted by Al 3+ ions, and additional Na + ions enter interstitial sites to compensate for the charge deficit. 8For glassy NAGP materials, the ionic conductivity also increases with x, and the addition of Al improves the glass-forming ability. 8The lithium analogs Li 1+x AlxGe 2−x (PO 4 ) 3 (LAGP) of these materials nucleate homogeneously 9 and have already found commercial usage as separator membranes in lithium-air batteries and related systems. 3,10omparatively little is known, however, about the Na-containing materials, although solid-state nuclear magnetic resonance (NMR) experiments on the LAGP 11 and NAGP 12 systems indicate that all structural aspects of the glass-to-crystal transition are the same, suggesting the same crystallization mechanisms.Homogeneous nucleation in NAGP glasses is supported by the results obtained from differential scanning calorimetry experiments on monolithic vs powdered samples. 8he structure of the NAGP materials was probed by combining neutron diffraction (ND) and x-ray diffraction (XRD) with 27 Al and 31 P magic angle spinning (MAS) and 31 P/ 23 Na double-resonance NMR spectroscopy.The structure of the glass was investigated in its as-prepared condition and after thermal annealing near the glasstransition temperature Tg in order to investigate the initial stages of thermal relaxation.The structure of the fully crystallized glass was also investigated by applying (i) the Rietveld method 13,14 to refine powder diffraction patterns and (ii) the pair-distribution function method.The results for the crystalline materials provided a starting point for interpreting the pair-distribution functions measured for the vitreous materials.Here, the analysis was helped by the provision of site-specific information on the coordination environments of P and Al from the solid-state NMR experiments.In turn, the pair-distribution functions provided information on the Ge coordination environment, which is largely inaccessible from NMR spectroscopy.
This paper is organized as follows.The diffraction theory for pair-distribution function analysis is outlined in Sec.II, and the experimental methods are described in Sec.III.The results are presented in Sec.IV and are discussed in Sec.V where a model is proposed for the glass structure.Conclusions are drawn in Sec.VI.

II. THEORY FOR PAIR-DISTRIBUTION FUNCTIONS
In a diffraction experiment, the total structure factor 15

S(k)
is measured, where k denotes the magnitude of the scattering vector; cα ≡ Nα/N denotes the atomic fraction of chemical species α, where Nα is the number of atoms of type α and N = ∑ α Nα is the total number of atoms; wα(k) represents the k-dependent x-ray atomic form-factor f α (k) or the k-independent coherent neutron scattering length bα of chemical species α; the mean value ⟨w(k)⟩ = ∑ α cαwα(k); and S αβ (k) is a Faber-Ziman partial structure factor.Neutral atom form-factors were used in the XRD data analysis. 16The neutron scattering lengths of the elements in NAGP materials are b Na = 3.63 (2) fm, b Al = 3.449 (5) fm, b Ge = 8.185 (20) fm, b P = 5.13(1) fm, and b O = 5.803(4) fm. 17he corresponding real-space information can be represented by the total pair-distribution function where r is a distance in real space and ⊗ denotes the onedimensional convolution operator.M(k) is a window function that originates from the finite k-range that is accessible to a diffractometer, which extends to a maximum value kmax.This accessibility has the effect of applying to S(k) − 1 the step modification function defined by M(k) = 1 for k ≤ kmax and M(k) = 0 for k > kmax.M(r) is the real-space manifestation of M(k) and is a symmetrical function.
In a neutron diffraction experiment, where the scattering lengths are k independent, where ρ is the atomic number density and g αβ (r) is a partial pairdistribution function.At r-values smaller than the distance of the closest approach between two atoms, g αβ (r) = 0.It follows that, in the absence of Fourier transform artifacts, such as those originating from M(k), D ′ (r) = −4πρr, which is often called the density line.
In order to distinguish between features in D ′ (r) that describe the glass structure from those that are an artifact of M(r), each peak i in rg αβ (r) can be represented by the Gaussian function where r i αβ , σ i αβ , and nβ α (i) are the peak position, standard deviation, and coordination number of chemical species β around α, respectively. 18The contribution of each peak p i αβ (r) to S(k) − 1 is given by where

ARTICLE scitation.org/journal/jcp
In the case of the ND data, the method of least squares was used to fit D ′ (r) to a sum of Gaussian functions convoluted with M(r) such that In the case of the XRD data, each p i αβ (k) function was Fourier transformed to real-space using the same M(k) function as used for the experimental data.For a selected region in r-space, the method of least squares was then used to fit an appropriate sum of these Fourier transforms to T ′ (r) = T(r) ⊗ M(r) using the program PXFIT (Hannon, private communication), where T(r) ≡ D(r) + T 0 (r) and T 0 (r) = 4πρr.For ease of comparison with the neutron diffraction results, the fitted functions are presented as The goodness-offit was assessed by the parameter 19

Rχ
calculated over the fitted range, where "exp" and "fit" denote the measured and fitted functions, respectively.The use of a step modification function in Eq. ( 2) not only leads to sharp peaks in D ′ (r), which is the reason for its application, but can also lead to Fourier transform artifacts that include large excursions below the density line D ′ (r) = −4πρr on each side of the first peak in D ′ (r).The latter can be suppressed by using the Lorch 20,21 modification function M(k) = sin(πk/kmax)/(πk/kmax) for k ≤ kmax and M(k) = 0 for k > kmax, but at the expense of broadened peaks.

A. Sample preparation and characterization
The glasses were made using the method described in Ref. 8. For each composition, part of the glass was kept in its as-prepared condition, part was annealed at a thermal treatment temperature T TT near to Tg for a time chosen to relax the glass structure while minimizing the formation of crystallites, and part was annealed at T TT for a time chosen to fully crystallize the material.For the x = 0 composition, two additional parts of the as-prepared glass were taken for further investigation.One part was annealed at T TT = 876 K for 0.25 h, and the other part was annealed at the same temperature for 0.5 h.
The glass transition temperature Tg and crystallization temperature Tx were measured for ∼20 mg monolithic samples using a differential scanning calorimeter (Netzsch DSC404) with a heating rate of 10 K/min.The mass density of the glasses was measured using a Quantachrome 1200e pycnometer operated with He gas.The density of the crystals was taken from the powder diffraction work.Parameters describing the samples and their thermal history are listed in Table I.

B. Powder diffraction experiments
The neutron powder diffraction patterns were measured using the fine resolution powder diffractometer E9 (FIREPOD) 22 at the reactor source BER II, Helmholtz-Zentrum Berlin (HZB), Germany.The datasets were collected at ambient temperature (295 K) using an incident wavelength λ = 1.7982(1)Å.The samples were handled in air and held in vanadium cans of 5 mm diameter.The x-ray powder diffraction patterns were measured using a Rigaku Smart-Lab system with Bragg-Brentano diffraction geometry, monochromatic copper Kα 1 radiation, and a fixed 0.66 ○ divergence slit.The instrument was equipped with a Johansson-type Ge-crystal monochromator and a HyPix-3000 2D single photon counting detector.
Rietveld refinement 13,14 of the crystal structures was performed using the FullProf suite. 23,24The ND and XRD patterns for each sample were refined simultaneously.The crystal structure data were standardized using STRUCTURE TIDY. 25 Crystal structures were drawn using the three-dimensional visualization system for electronic and structural analysis (VESTA) package, version 3.3.1. 26

C. Pair-distribution function experiments
Two sets of ND experiments were performed on the NAGP materials to measure the total pair-distribution functions.

TABLE I.
Parameters describing the Na 1+x AlxGe 2−x (PO 4 ) 3 samples and their thermal history.Glass transition temperature T g,as-prepared and number density ρ as-prepared of the as-prepared glass, crystallization temperature Tx and thermal treatment temperature T TT , time t anneal used to prepare the annealed glass from the as-prepared glass, glass transition temperature T g,annealed and number density ρ annealed of the annealed glass, thermal treatment time t cryst used to crystallize the as-prepared glass, number density ρ crystal of the crystalline phase, and melting point temperature T melt .
T g,as-prepared ρ as-prepared Tx T TT t anneal T g,annealed The Journal of Chemical Physics

ARTICLE scitation.org/journal/jcp
This instrument was chosen because of its stability and high count-rate over a wide k-range, which enables it to probe small differences in structure. 28The powdered samples were held in a cylindrical vanadium container of inner diameter 6.8 mm and wall thickness 0.1 mm.Diffraction patterns were measured at room temperature (≃298 K) for each of the samples in its container, the empty container, the empty instrument, and a cylindrical vanadium rod of diameter 6.08 mm for normalization purposes.The diffraction pattern was also measured for a slab of neutron absorbing 10 B 4 C of dimensions similar to the sample in order to estimate the effect of the sample's attenuation on the background count-rate at small scattering angles. 29The data analysis followed the procedure described in Ref. 30.
The second set of experiments used the diffractometer GEM at the ISIS pulsed neutron source. 31This instrument was chosen because it accesses a much larger kmax value than D4c, thus improving the resolution of the peaks in D ′ (r).The powdered samples were held in a cylindrical vanadium container of inner diameter 6 mm and wall thickness 40 μm.Diffraction patterns were measured at room temperature (≃298 K) for each of the samples in its container, the empty container, the empty instrument, and a cylindrical rod of the alloy V 0.9486 Nb 0.0514 (diameter 7.95 mm) for normalization purposes.The datasets were processed using the Gudrun analysis program 32 with inelasticity corrections calculated according to Ref. 33.
In each set of experiments, the relative counting times for the sample-in-container and empty container measurements were optimized in order to minimize the statistical error on the containercorrected intensity. 34Self-consistency checks were also made to assess the reliability of the measured functions.For instance, it is necessary that (i) each of the measured total structure factors satisfies the sum-rule relation ∫ (r) = −4πρr at r-values smaller than the distance of the closest approach between two atoms; and (iii) when these oscillations in D ′ (r) are set to the density-line values, the back Fourier transform is in good overall agreement with the original total structure factor.
The XRD experiment was performed on beamline 6-ID-D at the Advanced Photon Source with an incident photon energy of 80.019 keV.The scattered x-rays were counted using a Varex 4343CT amorphous silicon flat panel detector, which was placed at a distance of 271.932 mm from the sample position as deduced from the diffraction pattern measured for crystalline CeO 2 .The powdered glass samples were held in cylindrical Kapton polyimide tubes of 1.00(1) mm internal diameter and 0.051(6) mm wall thickness.Diffraction patterns were collected for each sample in its container, an empty container, and the empty instrument.The data were converted to one-dimensional diffraction patterns using FIT2D, 35 and the corrections for background scattering, beam polarization, attenuation, and Compton scattering were made using PDFgetX2. 36

D. Solid-state NMR experiments
Single-pulse 27 Al MAS NMR spectra were measured using a Bruker AVANCE 600 spectrometer with a 2.5 mm probe operated at 20.0 kHz.A pulse length of 0.34 μs, corresponding to a flip angle of π/6 and relaxation delays in the range 0.1-2 s, was used.Chemical shifts are reported against a 1M solution of Al(NO 3 ) 3 using Al(PO 3 ) 3 (isotropic chemical shift δ iso = −20.8ppm) as a secondary reference.Line shapes were fitted to the Czjzek model 37 simulating a wide distribution of quadrupolar coupling constants, as implemented in DMFit. 38olid-state 31 P NMR experiments were performed using an Agilent DD2 spectrometer interfaced with a 5.64 T magnet and a commercial 3.2 mm probe operated at a spinning speed of 10.0 kHz.Chemical shifts were referenced using powdered BPO 4 (δ iso = −29.27ppm against H 3 PO 4 85%) as a secondary reference.Experimental spectra were simulated using the DMFit software. 38 23Na{ 31 P} rotational echo double resonance (REDOR) experiments 39 were performed using a Bruker AVANCE 600 spectrometer with a 2.5 mm triple resonance probe and spinning frequency of 10 kHz, operated with the following conditions: 23 Na detection using a rotor synchronized Hahn spin-echo sequence with 180 ○ refocusing pulses of 3.25 μs length and a 31 P dipolar re-coupling π-pulse of 3.5 μs length, at a spinning speed of 15.0 kHz and a relaxation delay of 1 s. 31 P{ 23 Na} rotational echo adiabatic passage double resonance (REAPDOR) spectra were measured with the same apparatus used for the REDOR experiments and the pulse sequence proposed by Gullion, 40 applying all the 180 ○ pulses on the 31 P detection channel and the adiabatic passage pulse (1/3 of the length of the rotor period) to the 23 Na nuclei.Spinning and nutation frequencies were 10.0 and 80 kHz, respectively.

A. Powder diffraction
The NASICON aristotype crystallizes into space group-type R 3c, as originally determined by Hagman and Kierkegaard 41 for NaZr 2 (PO 4 ) 3 .In the solid solution series Na 1+x TixGe 2−x (PO 4 ) 3 at room temperature, the boundary phases NaTi 2 (PO 4 ) 3 (NTP) and NaGe 2 (PO 4 ) 3 (NGP) crystallize into this structure-type and the lower symmetry form R 3, respectively. 2,42The R 3c and R 3 space groups have a group-subgroup relation that leads to a splitting of crystallographic sites.According to the literature, the solid solutions crystallize into the NGP-type phase for x ≤ 1 or the NTP-type phase for x > 1. 2,42,43 The crystal structures of the NAGP compounds were refined using the Rietveld method.The x = 0 and x = 0.4 compositions were assigned to the NGP-type structure, whereas the x = 0.8 composition was assigned to the higher symmetry NTP-type structure.The Rietveld refinements are shown in Figs.1-3 where, for each experiment, the intensity is plotted vs sin θ/λ = k/4π and 2θ is the scattering angle. 15The corresponding crystal structure data are summarized in Tables II-IV where SOF represents a site occupation factor; B iso represents a Debye-Waller factor; and Rexp, Rp, and Rwp represent the expected, profile, and weighted profile R-factors, 44 respectively.
The crystalline NAGP compounds were found to be almost phase pure.For the x = 0 sample, the XRD pattern shows some peaks from a non-determined impurity phase, but according to the peak intensities, the phase fraction is <1%.For the x = 0.4 sample, no impurities could be identified.For the x = 0.8 sample, some impurity peaks and shoulders are present and the intensities suggest a phase fraction of a few percent.The diffraction background for the  II.

The Journal of Chemical Physics
x = 0 and x = 0.4 samples is flat, which suggests full crystallization.The background for the x = 0.8 sample is increased slightly, suggesting some amorphous residue.Maldonado-Manso et al. 43 reported differences between the nominal and actual compositions of NAGP and related compounds prepared via a solid-state reaction.To check for non-stoichiometry, different refinement strategies were tested.In the initial model, only the phosphorus SOF was set to unity.Refinement did not show any indication for non-stoichiometry at the oxygen sites, so the occupancies of these sites were also set to unity.The germanium site occupancies for the x = 0 sample were refined and showed negligible deviation from full occupancy.Hence, for the aluminum substituted samples, the combined occupation factor SOF(Ge) + SOF(Al) for each site was set to unity, but the Al/Ge-ratio was allowed to vary.The diffraction patterns for the x = 0 sample can be reproduced with two fully occupied sodium sites Na(1) and Na(2).Nevertheless, the aluminum substituted samples show occupancy of an additional Na(3) site.Accordingly, this site was included in the refinement of the x = 0 data with coordinates fixed at those obtained for the x = 0.4 sample.The refinement showed the Na(3) site to be unoccupied within the estimated standard deviation.For the NAGP samples with nominal compositions x = 0, 0.4, and 0.8,  III.
The crystal structures obtained from powder diffraction are illustrated in Fig. 4.They serve as a starting point for interpreting the pair-distribution functions.A summary of the coordination numbers and bond lengths is given in Table V.
The structures have a single phosphorus site in which a phosphorus atom is coordinated tetrahedrally by four symmetry independent oxygen atoms.The tetrahedra are quite regular with a narrow distribution of P-O distances ranging from 1.50 to 1.54 Å.The phosphorus site occupation was set to unity in all the Rietveld refinements.
In the lower symmetry R 3 structure, there are two sixfold coordinated germanium sites.The Ge  IV.
bond lengths is expected as Ge is substituted by Al in the NAGP compounds.
In the lower symmetry R 3 structure, there are two sixfold coordinated sodium sites Na(1) and Na(2).The Na(1) site is at the center of a trigonal prism formed by six symmetry equivalent oxygen atoms O(3), and the Na(2) site is similarly coordinated by O(4) oxygen atoms.In the higher symmetry R 3c structure, the Na(1) and Na(2) sites are equivalent.In the x = 0 and x = 0.4 compounds, the Na(1) site is fully occupied and the SOF was set to unity during the last refinement cycle.The Na(2) site is almost fully occupied and was refined.In the x = 0.8 compound, the single Na(1)/Na(2) site is fully occupied and its SOF was set to unity.The Na(1)-O and Na(2)-O distances are similar in all three compounds, ranging from 2.39 to 2.52 Å.The crystal structures also have a low-symmetry Na(3) site at the center of an irregular polyhedron of eight oxygen atoms with a broad distribution of Na(3)-O distances ranging from 2.04 to 2.78 Å (Table V).This Na(3) site also has two phosphorus atoms at a similar distance (2.6-2.8Å) to several of the oxygen atoms, and there is a third Na(3)-P distance of about 3 Å, which is slightly shorter than the Na-Na and the Na-Ge/Al distances.In the x = 0 compound, the Na(3) sites are unoccupied.In the x = 0.4 and x = 0.8 compounds, however, the Na(3) site occupancy is 14% vs 24%, respectively.According to the review by Guin and Tietz, 7 the sodium ions move from Na(1) to Na(2) sites via Na(3) sites.The high mobility leads to large Debye-Waller factors for the sodium atoms.

B. NMR spectroscopy
The 31 P MAS NMR spectra taken for both the glassy and crystallized NAGP samples with x = 0, 0.4, and 0.8 confirm the results obtained in previous work. 12The datasets for the x = 0 composition are shown in Fig. 5.The spectrum for the as-prepared glass is about 25 ppm wide and somewhat asymmetric.Its center of gravity [−28.9 ppm, Fig. 5(a)] suggests that P (3) units are the majority phosphate species, 12 where the superscript denotes the number of BO atoms per P atom.The spectrum differs significantly from that measured for the crystalline state [−37.4ppm, Fig. 5(d

Atom
Wyckoff-site   There are, however, many other possibilities for deconvoluting the 31 P MAS NMR spectrum into two Gaussians.Independent evidence in favor of this particular deconvolution is therefore required.Such evidence is provided by a REAPDOR experiment 40 where, in the 31 P{ 23 Na} double resonance experiment conducted here, the magnetic dipole-dipole coupling strength between the observed 31 P nuclei and the surrounding 23 Na nuclei was probed.This magnetic dipole coupling, which is averaged out by magicangle spinning, can be re-introduced by irradiating the 23 Na nuclei on resonance during the rotor period, resulting in a decrease in signal intensity.A comparison is then made between the intensities S 0 of the regular MAS NMR signal and S of the MAS NMR signal observed with simultaneous irradiation of the 23 Na nuclei for a fixed dipolar re-coupling time.The corresponding difference in signal intensity ΔS = S 0 − S depends on the strength of the dipolar coupling between the 31 P and 23 Na sites and thus has the potential to distinguish between different local environments in a glass.
Figure 6 shows the Fourier transforms of the echoes obtained after a 1.4 ms mixing time in NAGP glass with x = 0.The data indicate that the dipolar dephasing is indeed non-uniform as a function of the isotropic chemical shift: The 31 P nuclei on the highfrequency side of the resonance display stronger dipolar dephasing than those on the low-frequency side.Based on these data, it was possible to achieve a consistent deconvolution of the echo signal with dipolar dephasing (which emphasizes the less-strongly coupled 31 P nuclei) and the difference signal ΔS (which emphasizes the more-strongly coupled 31 P nuclei), in terms of two line shape components near −27(2) and −37(1) ppm (see Fig. 6).Table VI summarizes the line shape parameters for the two-component fit for the single-pulse spectra and the REAPDOR data.In the latter experiments, the ratio of areas of the high-frequency component to the low-frequency component was found to be lower than 2:1 for the S signal and higher than 2:1 for the ΔS signal, consistent with the difference in dipolar coupling strengths between the two sites.The component near −27 ppm is identified with P (3) units, and the component near −37 ppm is identified with P (4) units.The chemical shift of the latter is close to that found for the P (4) units in the crystalline state.VI.In (c), an additional small peak (violet curve) is included to account for a small amount of crystalline material.The sum of fitted functions is given by the red curve.In (d), the spectrum is fitted to a Gauss-Lorentz curve.
Figures 5(b) and 5(c) show the 31 P MAS NMR spectra for the annealed glasses with composition x = 0.In the case of the sample annealed at 873 K for 0.25 h, the data clearly show that the material remains completely amorphous.Some subtle line shape changes are, however, evident (Fig. S1), which may be attributed to a structural re-arrangement in the glass caused by the onset of the nucleation process.For example, the full-width at half-maximum (FWHM) of the line shape increases from 21.7(5) ppm for the as-prepared glass TABLE V. Element-oxygen distances for the first coordination spheres in crystalline NAGP with x = 0, 0.4, or 0.8.

Atom pair
Multiplicity x = 0 Distance (Å) x = 0.   P MAS-NMR spectra for glassy (as-prepared vs annealed) and crystalline NAGP with x = 0.The parameters describe the fractional area, δ iso , and FWHM of the peaks fitted to the spectra.For one of the annealed glasses, the crystalline contribution to the line shape is indicated by a subscript.The errors on the fitted areas are ±1% for the single-pulse spectra vs ±5% for the REAPDOR line shapes.The errors on δ iso and the FWHM are ±0.5 ppm for the single-pulse spectra vs ±2 ppm for the REAPDOR line shapes.5) ppm for the annealed glass.Table VI indicates that the change to the principal line shape involves a change to the ratio of the relative areas of the two line shape components, indicating a moderate increase in the fraction of P (4) units having the same local coordination environment as the crystalline phase.This finding was confirmed by fixing the chemical shifts and linewidths for the P (3) and P (4) sites at the average of the values listed in Table VI for the as-prepared and annealed glasses and fitting the spectra for both glassy samples with the peak areas as the only adjustable parameters.

Fitted
The peak areas showed a change in the ratio of P (4) :P (3) sites from 33.5(1.0):66.5(1.0) for the as-prepared glass to 36.0(1.0):64.0(1.0) for the annealed glass.In the case of the sample annealed at 876 K for 0.5 h, the change in peak areas is more dramatic, with the P (4) :P (3) ratio increasing to ∼1:1.For this sample, a small amount of crystalline material is also detected.Finally, the 23 Na/ 31 P double resonance experiments also allow for a comparison between the 31 P/ 23 Na dipolar coupling strengths of glassy and crystalline materials by plotting the normalized difference signal as a function of the dipolar mixing time NTr (= number of rotor cycles × rotor period).The corresponding 31 P{ 23 Na} REAP-DOR curves are compared in Fig. 7(a) and qualitatively illustrate that the internuclear 31 P/ 23 Na dipole-dipole interactions are significantly stronger in the crystalline rather than in the glassy material.A more quantitative estimate can be found on the basis of the 23 Na{ 31 P} REDOR results 39 shown in Fig. 7(b).In this experiment, the dephasing of the 23 Na nuclei in the dipolar field created by the 31 P spins is being observed.The simpler spin dynamics arising from the spin-1/2 character of the recoupled 31 P nuclei permits an approximate analysis in terms of dipolar second moments M 2(Na-P) , using the formula in the limit of short mixing times when ΔS/S 0 ≤ 0.2-0.3, 46where f is a calibration factor on the order of unity.The quantity M 2(Na-P) can be calculated from the internuclear distances using Van Vleck theory; 47 it is proportional to the sum of the inverse sixth powers of all Na-P internuclear distances present and can thus be computed from It follows that the mean Na-P (4) distance back-calculated from M 2(Na-P) in the NAGP crystal is shorter than the average of the Na-P (3) and Na-P (4) distances back-calculated from M 2(Na-P) in the glass.
The 27 Al MAS NMR spectra taken for the glassy and crystallized NAGP samples with x = 0.4 and 0.8 are shown in Fig. 8.The line shapes were fitted using a Czjzek model 37 that accounts for a distribution of quadrupolar coupling constants originating from a distribution of local field gradients.The fitted parameters are listed in Table S1.For the glasses, the spectra feature broad resonances at isotropic chemical shifts around 46(1), 15(1), and −12(1) ppm that can be attributed to aluminum in fourfold, fivefold, and sixfold coordination environments, respectively.For some of the glassy samples, a sharp peak near −20 ppm is observed, which indicates a small amount of a yet unidentified crystalline impurity with aluminum in sixfold coordination.For the x = 0.8 composition, the presence of this site in the as-prepared glass, but its absence in the annealed material, suggests a non-uniform distribution of more ordered regions when the glass is first formed.
For the crystallized materials prepared by thermally annealing the glass for an extended period of time, the 27 Al MAS NMR spectra show that crystalline NASICON is the dominant phase.The dataset for the x = 0.8 composition shows a small amount of fourfold coordinated aluminum in an AlPO 4 impurity phase.The results for both compositions indicate that a small fraction of aluminum atoms remain in a glass-like domain where they occur in all three coordination environments, which contributes toward an average Al-O coordination number less than six (Table VII).
In Table VII, the results for the aluminum speciation are compared to those obtained for a different set of NAGP samples in previous work. 12For the as-prepared glass, there is some variation in the speciation that may be influenced by the sample preparation procedure although the NMR spectra in Ref. 12 were measured at a lower field strength of 5.7 T (cf.(8)   discrepancies between the intensities of Bragg peaks measured using GEM vs D4c that originate from (i) differences between the kspace resolution function of the diffractometers and (ii) the possibility of some preferred orientation of the crystallites (Fig. 9).For the glassy materials, the GEM and D4c datasets for the asprepared glasses are in overall agreement (Fig. 10).The XRD patterns indicate a small amount of crystallinity for the annealed x = 0 and x = 0.8 samples (Fig. 12) that is not observed using ND with D4c, which may reflect the enhanced k-space resolution of the x-ray diffractometer.For all the glassy materials, the S(k) functions display a shoulder at k 1 ∼ 0.85 Å −1 that is indicative of ordering on a real-space length scale of 2π/k 1 ∼7.4 Å.In Na 2 O-GeO 2 -P 2 O5 glasses, the appearance of a similar feature has been linked to an inhomogeneous distribution of Na and Ge atoms. 48 Structure of NAGP with x = 0 For crystalline NAGP with x = 0, the D ′ (r) functions measured using GEM and XRD are shown in Fig. 13.By comparison with the powder diffraction work (Table V), the first peak at 1.52 Å and the second peak at 1.86 Å originate from P-O and Ge-O correlations, respectively, and the third peak will have a contribution from both Na-O and O-O correlations.The D ′ (r) functions were therefore fitted by using the following protocol.
A single peak was used to represent the P-O correlations.A single peak was also used to represent the Ge-O correlations in the ND work, whereas two peaks were used to represent these correlations in the XRD work, where the large atomic number of Ge (Z = 32) makes the results particularly sensitive to the germanium coordination environment.A single peak was used to represent the Na-O correlations with a coordination number fixed at the value found from the Rietveld refinement.Two peaks were used to represent the O-O correlations because each BO atom is shared between a PO 4 tetrahedron and a GeO 6 octahedron, i.e., it has three nearest-neighbor O atoms within a tetrahedron and four nearest-neighbor O atoms at a longer distance within an octahedron (the fifth O atom within an octahedron is at an even longer next nearest-neighbor distance).The P-O coordination number found from both ND and XRD is smaller than the value nO P = 4 expected from Rietveld refinement and from 31 P MAS NMR experiments (Secs.IV A and IV B).In the latter, the measured chemical shifts are typical of phosphorus in a fourfold coordination environment and are quite different to the chemical shift of phosphorus in a PO 3 coordination environment, which is observed near 100 ppm in the crystalline molecular systems P 4 O7 and P 4 O 6 S, as referenced to an aqueous solution of H 3 PO 4 . 49A measured P-O coordination number nO P < 4 is typical of diffraction work on phosphate materials, 48,50,51 where the shortfall is often attributed to the k-space resolution of the diffractometer. 19,52everal steps were taken in order to investigate the origin of the reduced nO P values.In this analysis, the GEM dataset for crystalline NAGP with x = 0 was used as an exemplar because (i) the large kmax value of ≃39 Å −1 leads to well-resolved peaks in D ′ (r), (ii) the analysis is simplified by the absence of Al-O correlations, and (iii) the crystal structure is known.(1) The coherent neutron scattering length of P is reported to be in the range 5.0-5.3fm [53][54][55] with a recommended value b P = 5.13(1) fm. 17The data analysis was repeated using b P values that extend over this range, but a notable increase in nO P was not observed.(2) The dataset was reanalyzed using the density ρ = 0.0843(1) Å −3 measured by He pycnometry.A fit to the D ′ (r) function using the strategy reported above led, however, to a reduced coordination number nO P = 3.44(5) with nO Ge = 5.71( 5).

The Journal of Chemical Physics
(3) The Lorch 20 modification function was used in Eq. ( 2) to suppress Fourier transform artifacts, which facilitates the extraction of coordination numbers by direct integration over the peaks in D ′ (r).Here, the mean coordination number of atoms of type β, contained in a volume defined by two concentric spheres of radii ri and rj centered on an atom of type α, is given by The resultant D ′ (r) function is shown in Fig. 13 V).
The fitted GEM and XRD D ′ (r) functions are shown in Figs.14(a) and 15(a), respectively, and the peak parameters are summarized in Tables S2 and S3.Inclusion of the second P-O peak increases the P-O coordination number to nO P ≃ 3.9 and leads to a Ge-O coordination number nO Ge = 5.90(4)-5.95(5), which compares to an expected value of nO Ge = 6.The Ge-O distance of 1.857(2) Å obtained from ND is in agreement with the mean Rietveld value of 1.857 Å (Table V), and the Ge-O distances of 1.847(1) and 1.871(1) Å obtained from XRD compare to average short and long Rietveld Ge-O distances of 1.841 and 1.873 Å, respectively (Table V).
The GEM results for the crystalline material (Table S2) provided starting parameters for fitting the D4c dataset.The fitted D ′ (r) function is shown in Fig. 16(a), where the peaks are broadened relative to GEM because of the smaller value kmax = 22.25 Å −1 of D4c and the peak parameters are summarized in Table S4.The fit gives nO P = 4.00(4) with nO Ge = 5.97 (5).For the as-prepared glass, the GEM and XRD datasets were fitted using two Gaussian peaks representing the nearest-neighbor P-O correlations and one Gaussian peak representing the nearestneighbor Ge-O correlations.Additional peaks at larger-r values were used in order to constrain the fitted functions at lower r.S5 and S6.The GEM results provided starting parameters for fitting the D4c dataset for the as-prepared glass, where the fitted D ′ (r) function is shown in Fig. 16(b) and the P-O and Ge-O peak parameters are summarized in Table S7.The results give a P-O coordination number of nO P = 4 within the experimental error and an average Ge-O coordination number nO Ge ∼ 5 that is significantly smaller than found for the crystalline phase.
For the annealed glass, the fits were initiated using starting parameters obtained for the as-prepared glass.The fitted D ′ (r) functions are shown in Fig. 15(c) (XRD) and in Fig. 16(c) and   S8 and S9.The results do not indicate a change to nO P but show that, relative to the as-prepared glass, nO Ge increases with the time of thermal annealing.

Structure of NAGP with x > 0
The GEM and XRD D ′ (r) functions for crystalline and glassy NAGP with x = 0.4 and x = 0.8 are shown in Figs.17-20.By comparison with the powder diffraction results (Table V), the first peak in each D ′ (r) function at ≃1.52 Å is attributed to P-O correlations, the second peak at ≃1.86 Å will have contributions from both Ge-O and Al-O correlations, and the third peak will have contributions from both Na-O and O-O correlations.In view of the results obtained in Sec.IV C 1, the nearest-neighbor P-O correlations in D ′ (r) were represented by two Gaussian peaks.The parameters describing the Gaussian peaks for the nearest-neighbor Al-O correlations were constrained using the information obtained from other techniques.
For the crystalline materials, the Al-O coordination number was fixed at nO Al = 6 and the position of the Gaussian peak compositions are shown in Figs.17(a) and 18(a), respectively, and the peak parameters are summarized in Table S2.The fitted XRD D ′ (r) functions for the x = 0.4 and x = 0.8 compositions are shown in Figs.19(a) and 20(a), respectively, and the peak parameters are summarized in Table S3.The P-O coordination numbers are in the range 3.85(3)-3.99(3),and within the experimental error, the Ge-O coordination number is consistent with the powder diffraction value nO Ge = 6.The GEM results (Table S2) provided starting parameters for fitting the D4c datasets for the crystalline materials.The fitted D ′ (r) functions for the x = 0.4 and x = 0.8 compositions are shown in Figs.21(a) and 22(a), respectively, and the fitted P-O and Ge-O parameters are summarized in Table S4.The fitted coordination numbers are nO P = 4 and nO Ge = 6 within the experimental error.
For the glassy materials, three Gaussian peaks were used to represent the nearest-neighbor Al-O correlations for fourfold, fivefold, and sixfold coordinated aluminum atoms.The relative abundance of these species was fixed in accordance with the aluminum speciation found from 27 Al MAS NMR experiments (Table VII).The peak positions for fourfold and fivefold coordinated aluminum atoms were set at 1.76 and 1.84 Å, respectively, which are the Al-O bond lengths calculated from the bond valence method. 56The position of the peak for sixfold coordinated aluminum atoms was set at the average Al-O bond length of 1.86 Å found from the powder diffraction work (Table V).The width of the Gaussian peak describing AlO 4 units was set at 0.05 Å, which is a typical value for AlO 4 units in alumino-phosphate glasses, 50 and the width of the Gaussian peak describing AlO 6 units was set at 0.07 Å, which is representative of the value found in the crystalline NAGP materials with x = 0.4 and x = 0.8 (Table S2).The width of the Gaussian peak describing AlO5 units was set in between these values at 0.06 Å.Additional peaks at larger-r values were used in order to constrain the fitted functions at lower r.

The Journal of Chemical Physics
For the as-prepared glasses with x = 0.  fitted peak parameters are summarized in Table S7.The fitted XRD D ′ (r) functions for the as-prepared glasses with x = 0.4 and x = 0.8 are shown in Figs.19(b) and 20(b), respectively, and the fitted P-O and Ge-O peak parameters are summarized in Table S6.The results give a P-O coordination number of nO P = 4 within the experimental error and indicate a reduction in the average Ge-O coordination number as more aluminum is added to the glass structure.
The XRD and D4c results obtained for the as-prepared glasses provided starting parameters for fitting the datasets for the annealed glasses, where the aluminum speciation was taken from 27 Al MAS NMR experiments (Table VII).The fitted XRD D ′ (r) functions for the x = 0.4 and x = 0.8 compositions are shown in Figs.19(c) and 20(c), respectively, and the fitted P-O and Ge-O peak parameters are summarized in Table S8.The fitted D4c D ′ (r) functions for the x = 0.4 and x = 0.8 compositions are shown in Figs.21(c) and 22(c), respectively, and the fitted P-O and Ge-O peak parameters are summarized in Table S9.The results do not show a change to nO P but, for each composition, indicate an increase in nO Ge relative to the as-prepared glass.V. DISCUSSION A. Crystal vs glass structure

The Journal of Chemical Physics
The P-O and Ge-O bond distances and coordination numbers obtained from the ND and XRD pair-distribution function analyses are summarized in Table VIII.The diffraction results for the glasses show Ge-O coordination numbers that are markedly smaller than nO Ge = 6, which indicates the presence of polyhedral GeO5 and/or GeO 4 units.This observation of a substantial difference between the crystal and glass structures is supported by solid-state NMR experiments. 12For example, the 31 P chemical shifts differ by more than 10 ppm; P-O-P linkages are absent in the crystal but may be present at low concentrations in the glass; and P (3) phosphate units dominate the structure of the glasses but are absent in the isochemical crystals where all the phosphate units are P (4) .In the case of NAGP with x = 0, 23 Na spin-echo decay spectroscopy finds that the second moment quantifying the strength of Na-Na interactions is about twice as large in the glass as compared to the isochemical crystal, which suggests that some Na-Na distances are shorter in the glass than in the crystal.Furthermore, Al is exclusively six-coordinated in the aliovalently (Al + Na) substituted crystalline materials but is found in a mixture of fourfold, fivefold, and sixfold coordination environments in the glass (Table VII).In order to investigate the network connectivity via the BO vs non-bridging oxygen (NBO) content of the NAGP materials, it is instructive to consider a simple model based on the oxygen coordination number. 48Let us assume that all oxygen atoms are part of a network and occupy either NBO or BO sites with coordination numbers of one or two, respectively.Then, the fraction of NBO atoms is given by where N NBO is the number of NBO atoms and N O is the total number of oxygen atoms.
In the following, it will be assumed that all the nearest-neighbor oxygen atoms to Ge and Al are BO atoms, in accordance with the crystal structures where each of these oxygen atoms is shared with a phosphorus atom.The number of NBO atoms per phosphorus atom is given by N NBO /N P = (c O /c P ) f NBO , where N P is the number of phosphorus atoms.Hence, if the P-O coordination number nO P = 4 and NBO atoms reside on phosphorus atoms alone, the number of BO atoms per phosphorus atom In the case of glassy (GeO 2 )x(NaPO 3 ) 1−x , for example, O-1s x-ray photoemission spectroscopy (XPS) experiments indicate that most of the NBO atoms are bound to P and not Ge, e.g., 100% for x = 0 vs 82.5% for x = 0.5. 65For M 2 O-P 2 O5 materials, Eq. ( 12) will also give the nearest-neighbor P-P coordination number nP P .The number of NBO atoms per M atom N NBO /N M = (c O /c M ) f NBO , where N M is the number of M atoms.The number of BO atoms per M atom The nearest-neighbor O-O coordination number is given by where the fraction of BO atoms f BO = 1 − f NBO .IX) match the measured values.For the cases of crystalline NAGP 2,43 and LAGP, 63 the Al 3+ ions substitute at the sites of the Ge 4+ ions such that nO Al = nO Ge = 6, so there is no change to N NBO /N P or N NBO /N M and f NBO = 0.By contrast, in glassy NAGP, the Ge-O and Al-O coordination numbers take values smaller than 6, so NBO atoms appear.The fraction of NBO atoms increases as the composition is changed from x = 0 to x = 0.8, leading to an increase in N NBO /N P .The N NBO /N P ratios are consistent with a mixture of P (3) and P (4) species, as found from 31 P solid-state NMR experiments (Sec.IV B).If all the phosphate groups are either P (3) or P (4) , it follows that the fraction of P (3) species for the as-prepared glass is 0.618(44) at x = 0 vs 0.770(52) at x = 0.8.The annealing procedure leads to a reduction in f NBO for each of the glasses, corresponding to a small increase in the Ge-O and Al-O coordination numbers, which leads to a decrease in both N NBO /N P and N NBO /N Na .

C. Super-structural units
A challenge is to provide a structural model for NAGP glasses in which the phosphate species are constrained to be either P (3) or P (4) and the mean Ge-O and Al-O coordination numbers are both substantially larger than four.A starting point is provided by the model of Ren and Eckert 66 for sodium phosphosilicate glasses, which invokes super-structural units containing sixfold coordinated silicon atoms (Fig. 24).For simplicity, the focus in the following will be on the NAGP glass with x = 0. BO atoms will be denoted by Ø, and as mentioned before, the number of BO atoms per Ge or P atom will be denoted by a superscript.
The NAGP glass composition for x = 0 can be written as (Na 2 P 6 GeO 18 ) 1/2 (GeO 2 ) 3/2 .In the model of Ren and Eckert, the Na 2 P 6 GeO 18 part of the glass is built from [Ge (6) super-structural units, where the negative charge on a unit is compensated by two Na + ions, and the GeO 2 part of the glass is built from charge-neutral corner-sharing tetrahedral Ge (4) units.The model therefore contains only P (3) , Ge (4) , and Ge (6) species.All NBO atoms reside within the P (3) motifs of a super-structural unit, where charge-shuttling ensures that they carry a partial charge of −(2/6)e = −0.33eand e is the elementary charge.The model delivers a coordination number nO Ge = 4.5 for NAGP glass with x = 0, which is smaller than the value found from the diffraction experiments (Table VIII).comprises at least three different subtypes, depending on whether the BO atoms are connected to three Ge, two Ge and one Al, or one Ge and two Al atoms (the environment with three Al atoms will have a negligible statistical probability).In principle, one could try to identify these neighbors with the help of 31 P{ 23 Na} REAP-DOR in combination with 31 P{ 27 Al} REAPDOR experiments, but the analysis will be complicated.

E. Ion transport
The electrical conductivity in electrolytes with the NASICON structure is considered to be predominantly ionic, with an ion transport number approaching unity. 68For these materials, in which the electrical conduction originates from the migration of a single charge carrier species, the ionic conductivity is given by σ = neμ, where n is the effective carrier concentration, e is the charge on the mobile Na + ions, and μ is the ion mobility.In the case of polycrystalline materials, there will be contributions to this conductivity from ion transport both in grains and at grain boundaries.A contribution from grain boundaries will not feature in conductivity measurements on glass, provided that the samples are unpowdered.
In crystalline NAGP, the electrical conductivity increases by almost two orders of magnitude and the activation energy Ea for the conductivity decreases as the composition is changed from x = 0 to x = 0.8 [Figs.23(a) and 23(b)], i.e., as Ge 4+ ions are aliovalently substituted by Al 3+ ions and additional Na + ions are introduced to compensate for the charge deficit.This substitution does not generate NBO atoms, i.e., in the absence of coordination defects, f NBO = 0 across the composition range.The conductivity is between 3.2 and 40 times larger than that of the isochemical as-prepared glass. 8or the glass, the composition dependence of σ for the crystalline materials can be roughly matched by scaling [Fig.23(a)].The conductivity of the glass also increases with the concentration of Na + ions, but the mobility of these ions is smaller in the more disordered glass structure.NBO atoms are present in the glass as part of the P (3) motifs on the super-structural units (Sec.V C).They are negatively polarized and presumably influence the Na + ion mobility, although the role that they play is unclear.While it has been suggested that NBO atoms at the P (2) sites in phosphate glasses can act as local traps for Na + ions, 69 the traps generated by NBO atoms at the P (3) sites of the super-structural units carry only a fractional negative charge and are thus expected to be significantly more shallow.

VI. CONCLUSIONS
The structure of crystalline NAGP with x = 0, 0.4, and 0.8 was investigated using neutron and x-ray powder diffraction and pair-distribution function analysis.The latter indicates deformable tetrahedral PO 4 units in that the distribution of P-O distances is broader than can be accommodated by using a single Gaussian peak to describe the nearest-neighbor P-O correlations.The ionic conductivity increases with x as more Na + ions are introduced.The ion transport does not involve NBO atoms because they are absent from the crystal structures.
The structure of the corresponding glasses was investigated by neutron and x-ray pair-distribution function analysis and solidstate NMR spectroscopy.There is a substantial reduction in the Ge-O and Al-O coordination numbers from the value nO Ge = nO Al = 6 found in the crystalline phases, which leads to the creation of NBO atoms.The fraction of these atoms f NBO is quantified via a structural model based on the oxygen coordination number.The ionic conductivity of the glass increases with the concentration of Na + ions, an enhancement that is accompanied by an increase in both f NBO and the ratio NBO:P vs a decrease in the ratio NBO:Na.The effect of NBO atoms on the Na + ion mobility is in need of further investigation.
A structural model is proposed for the x = 0 glass by adapting the model of Ren and Eckert 66 for phosphosilicate glasses.Here, the super-structural units of this model can be modified by using the donor function of a doubly bonded NBO atom on a P (3) motif to convert Ge (4) atoms to higher-coordinated germanium atoms, thereby converting P (3) to P (4) motifs and increasing the size of a super-structural unit.The P (4) species have a formal positive charge so that they do not attract sodium ions.Accordingly, they interact more weakly with Na + ions than the P (3) species.Justification for the structural model is provided by deconvolution of the 31 P MAS NMR spectra for the as-prepared glass that is consistent with the 31 P{ 23 Na} REAPDOR results.
When the as-prepared glass is thermally annealed, the 31 P MAS NMR spectra show an increase in the fraction of P (4) species, indicating a further progression to the super-structural modification process.This modification is also seen by a moderate increase in the Ge-O coordination number found from the diffraction experiments and a concomitant reduction in f NBO .The P (4) species generated within the modified super-structural units by the reaction between P (3) and Ge (4) motifs therefore provide the nucleation sites for crystal growth via a homogeneous nucleation mechanism.See the supplementary material for Fig. S1 that gives a more detailed comparison of the 31 P MAS NMR spectra shown in Fig. 5 for the as-prepared glass and the glass annealed at T TT = 873 K for 0.25 h. Figure S2 shows the fitted D ′ (r) functions measured using D4c for the NAGP samples with x = 0 annealed at T TT = 876 K for either 0.25 or 0.5 h.Table S1 lists the parameters obtained from the fits to the 27 Al MAS NMR spectra shown in Fig. 8. Tables S2-S9 list the parameters obtained from the fits to the real-space functions measured by ND or XRD that are shown in Figs.14-22 and Fig. S2.

FIG. 1 .
FIG. 1. Rietveld refinement results for the crystal structure of NAGP with x = 0.Each panel shows the measured intensity I obs (red curve), calculated intensity I calc (black curve), and difference of intensities I obs − I calc (blue curve), along with the positions of the Bragg reflections allowed by the space group for the structure (vertical black bars).The refined composition is Na 0.99(2) Ge 1.96(2) (PO 4 ) 3 , and the refinement details are given in TableII.

FIG. 2 .
FIG. 2. Rietveld refinement results for the crystal structure of NAGP with x = 0.4.Each panel shows the measured intensity I obs (red curve), calculated intensity I calc (black curve), and difference of intensities I obs − I calc (blue curve), along with the positions of the Bragg reflections allowed by the space group for the structure (vertical black bars).The refined composition is Na 1.42(2) Al 0.43 Ge 1.57(2) (PO 4 ) 3 , and the refinement details are given in TableIII.

FIG. 3 .
FIG. 3. Rietveld refinement results for the crystal structure of NAGP with x = 0.8.Each panel shows the measured intensity I obs (red curve), calculated intensity I calc (black curve), and difference of intensities I obs − I calc (blue curve), along with the positions of the Bragg reflections allowed by the space group for the structure (vertical black bars).The refined composition is Na 1.72(2) Al 0.804 Ge 1.196(12) (PO 4 ) 3 , and the refinement details are given in TableIV.

FIG. 4 .FIG. 5 .
FIG. 4. The R 3-type NAGP crystal structure.The unit cell is shown top left, the GeO 6 -octahedra are shown middle left, and the PO 4 -tetrahedron is shown bottom left.The different sodium centered polyhedra are shown on the right.The positions of the different sodium centered polyhedra in the unit cell are shown in the center.In the higher symmetry R 3c structure, several of the sites are equivalent: Ge(1) = Ge(2), Na(1) = Na(2), O(1) = O(2), and O(3) = O(4).

FIG. 6 .
FIG. 6.The 31 P{ 23 Na} REAPDOR results for NAGP glass with x = 0.The Fourier transforms of S (top) and ΔS = S 0 − S (bottom) are plotted, where the spectra were obtained for a dipolar mixing time of 1.4 ms.Each spectrum is fitted to two Gaussian peaks (blue and green curves), and the sum of fitted functions is given by the red curve.The weaker signals are attributed to spinning sidebands.

FIG. 7 .
FIG. 7.The (a)31 P{23 Na} REAPDOR and (b)23 Na{ 31 P} REDOR datasets for asprepared glassy vs crystalline NAGP with x = 0.In (b), the solid curves show fits to the datasets using Eq.(8) for the data range ΔS/S 0 ≤ 0.3.The results indicate consistently that the magnetic dipole-dipole coupling between the nuclei23 Na and 31 P is stronger in the crystal as compared to the glass.

FIG. 8 .
FIG. 8.The single-pulse27 Al MAS NMR spectra (black curves) measured for (a) as-prepared glassy, (b) annealed glassy, and (c) crystalline NAGP with (I) x = 0.4 or (II) x = 0.8.The Czjzek fits to the spectral components are given by the blue or cyan curves for Al(IV), the green curve for Al(V), and the violet or magenta curves for Al(VI).The sum of fitted functions is given by the red curve.
14.1 T in the present work) where the resolution of the individual 27 Al resonances is substantially lower.The results show an increase in the mean Al-O coordination number when the as-prepared glass is thermally annealed.C.Pair-distribution functionsThe S(k) functions measured by ND (using GEM or D4c) and XRD are shown in Figs.9-12.For the crystalline materials, there are

FIG. 9 .
FIG. 9.The total structure factors S(k) for crystalline NAGP with x = 0, 0.4, and 0.8 measured using GEM (red curves) or D4c (black curves).The vertical error bars are smaller than the line thickness at most k values.For clarity of presentation, several of the curves are displaced vertically and the GEM datasets are shown only to 30 Å −1 .

FIG. 10 .FIG. 11 .
FIG. 10.The total structure factors S(k) for as-prepared vs annealed glassy NAGP with x = 0, 0.4, and 0.8.The red curves give the GEM datasets, and the black or blue curves give the D4c datasets, where the latter correspond to the x = 0 glass annealed at T TT = 876 K for either 0.25 or 0.5 h.The vertical error bars are smaller than the line thickness at most k values.For clarity of presentation, several of the curves are displaced vertically and the GEM datasets are shown only to 30 Å −1 .

FIG. 12 .
FIG.12.The total structure factors S(k) for the as-prepared (black curves) vs annealed (red curves) glassy NAGP with x = 0, 0.4, and 0.8 measured using XRD.Several of the curves are displaced vertically for clarity of presentation.

FIG. 13 .
FIG. 13.The D ′ (r) functions for crystalline NAGP with x = 0 as obtained from the GEM data by the application of either (a) a Lorch or (b) a step modification function with kmax = 39.1 Å −1 and from (c) the XRD data by the application of a step modification function with kmax = 24.99Å −1 .In (a), the vertical downward arrow indicates the minimum after the second peak.In (b) and (c), the measured functions (filled circles) are fitted using a single Gaussian peak to describe the P-O correlations (blue broken curve).The other peaks show the fitted contributions from the Ge-O [red solid curve(s)], Na-O (magenta solid curve), and O-O (green broken curves) correlations.The black solid curves give the overall fits.The displaced green solid curves show the residuals where (b) Rχ = 5.06% or (c) Rχ = 2.97% for the fitted range 1.30-2.70Å.In (a) and (b), the curves are displaced vertically for clarity of presentation.
(a).Integration to the first minimum at 1.67 Å gives nO P = 3.7(1), whereas integration to the second minimum at 2.19 Å gives an increased value nO P = 4.3(2) if the Ge-O coordination number is set at nO Ge = 6.The measured D ′ (r) function therefore has sufficient area to accommodate a P-O coordination number of 4, which suggests a broader distribution of P-O distances than found by fitting a single peak.(4) In view of the findings in step (3), the fitting procedure was adapted to include a second P-O Gaussian peak.This requirement may originate from different coordination environments for the BO atoms in a phosphate tetrahedron.All these BO atoms are shared between P and Ge atoms but only two [O(3) and O(4)] have nearest-neighbor sodium atoms, where the Na-O distance is 2.39-2.52Å (Table

FIG. 14 .FIG. 15 .
FIG. 14.The fitted D ′ (r) functions for (a) crystalline and (b) as-prepared glassy NAGP with x = 0 measured using GEM.The filled circles give the measured function; the black solid curve gives the fitted function; and the other curves show the contributions from the P-O (blue broken curves), Ge-O (red solid curve), Na-O (magenta solid curve), and O-O (green broken curves) correlations.The displaced green solid curve shows the residual.For the glass, the Na-O and O-O correlations are introduced to constrain the peaks fitted at smaller r-values.
The fitted D ′ (r) functions are shown in Fig. 14(b) (ND) and Fig. 15(b) (XRD), and the P-O and Ge-O peak parameters are summarized in Tables

FIG. 16 .
FIG. 16.The fitted D ′ (r) functions for (a) crystalline, (b) as-prepared glassy, and (c) annealed glassy (T TT = 873 K) NAGP with x = 0 measured using D4c.The filled circles give the measured function; the black solid curve gives the fitted function; and the other curves show the contributions from the P-O (blue broken curves), Ge-O (red solid curve), Na-O (magenta solid curve), and O-O (green broken curves) correlations.The displaced green solid curve shows the residual.For the glass, the Na-O and O-O correlations are introduced to constrain the peaks fitted at smaller r-values.

Fig
Fig. S2 (ND).The P-O and Ge-O peak parameters are summarized in TablesS8 and S9.The results do not indicate a change to nO P but show that, relative to the as-prepared glass, nO Ge increases with the time of thermal annealing.

FIG. 17 . 8 FIG. 18 .FIG. 19 .
FIG. 17.The fitted D ′ (r) functions for (a) crystalline and (b) as-prepared glassy NAGP with x = 0.4 measured using GEM.The filled circles give the measured function; the black solid curve gives the fitted function; and the other curves show the contributions from the P-O (blue broken curves), Al-O (cyan solid curves), Ge-O (red solid curve), Na-O (magenta solid curve), and O-O (green broken curves) correlations.The displaced green solid curve shows the residual.For the glass, the Na-O and O-O correlations are introduced to constrain the peaks fitted at smaller r-values.

FIG. 20 .
FIG. 20.The fitted D ′ (r) functions for (a) crystalline, (b) as-prepared glassy, and (c) annealed glassy NAGP with x = 0.8 measured using XRD.The filled circles give the measured function; the black solid curve gives the fitted function; and the other curves show the contributions from the P-O (blue broken curves), Al-O (cyan solid curves), Ge-O [red solid curve(s)], Na-O (magenta solid curve), and O-O (green broken curves) correlations.The displaced green solid curve shows the residual.For the glass, the Na-O and O-O correlations are introduced to constrain the peaks fitted at smaller r-values.
4 and x = 0.8, the fitted GEM D ′ (r) functions are shown in Figs.17(b) and 18(b), respectively, and the fitted P-O and Ge-O peak parameters are summarized in Table S5.The results provided starting parameters for fitting the D4c datasets for the as-prepared glasses.The fitted D4c D ′ (r) functions for the x = 0.4 and x = 0.8 compositions are shown in Figs.21(b) and 22(b), respectively, and the

FIG. 21 .
FIG. 21.The fitted D ′ (r) functions for (a) crystalline, (b) as-prepared glassy, and (c) annealed glassy NAGP with x = 0.4 measured using D4c.The filled circles give the measured function; the black solid curve gives the fitted function; and the other curves show the contributions from the P-O (blue broken curves), Al-O (cyan solid curves), Ge-O (red solid curve), Na-O (magenta solid curve), and O-O (green broken curves) correlations.The displaced green solid curve shows the residual.For the glass, the Na-O and O-O correlations are introduced to constrain the peaks fitted at smaller r-values.

FIG. 22 .
FIG. 22.The fitted D ′ (r) functions for (a) crystalline, (b) as-prepared glassy, and (c) annealed glassy NAGP with x = 0.8 measured using D4c.The filled circles give the measured function; the black solid curve gives the fitted function; and the other curves show the contributions from the P-O (blue broken curves), Al-O (cyan solid curves), Ge-O (red solid curve), Na-O (magenta solid curve), and O-O (green broken curves) correlations.The displaced green solid curve shows the residual.For the glass, the Na-O and O-O correlations are introduced to constrain the peaks fitted at smaller r-values.

B
. Structural model for M 2 O-Al 2 O 3 -GeO 2 -P 2 O 5 materials For a fully polymerized network, nP O + nGe O + nAl O = 2 and f NBO = 0.For a fully depolymerized network, nP O + nGe O + nAl O = 1 and f NBO = 1.The model does not allow for the formation of threefold coordinated oxygen atoms as supported, for example, by 17 O NMR experiments on GeO 2 -P 2 O5 glass. 64The number of O-β bonds between oxygen and atoms of chemical species β is equal to the number of β-O bonds, i.e., N O nβ O = N β nO β , where N β is the number of β atoms, such that c O nβ O = c β nO β .Hence, the fraction of NBO atoms can be re-written as f NBO = (2c O − c P nO P − c Ge nO Ge − c Al nO Al )/c O .
.M.R., J.F.O.-M., and A.M.N.-M.tailored the heattreatments and discussed the ionic conductivity results.J.F.O.-M.and A.M.N.-M.prepared the samples.H.A. designed the powder diffraction work and analyzed the data.P.S.S. and A.Z. designed the pair-distribution function work; L.V.D.G., A.Z., P.S.S., R.M.D.S., and C.J.B. performed the diffraction experiments; and L.V.D.G. and R.M.D.S. analyzed the data.I.D.A.S. and H.E. performed the NMR experiments and analyzed the data.L.V.D.G. and P.S.S. devised the model in Sec.V C. P.S.S., H.A., and H.E. wrote this paper with input from all the co-authors.

TABLE VI .
Deconvolution of the 31

TABLE VII .
27e Al speciation in the NAGP materials with x = 0.4 and x = 0.8 found from the27Al MAS NMR experiments.The mean Al-O coordination number is also given.The error on the fraction of Al in a fourfold, fivefold, or sixfold coordination environment is ±2% in the present work vs ±3% in Ref. 12. Crystalline coordination environments are identified by a subscript.

TABLE VIII .
The mean P-O and Ge-O bond distances and coordination numbers obtained from Gaussian peak fits to the r-space functions measured using ND and XRD for the as-prepared and annealed NAGP glasses and corresponding crystalline phases.The shorter and longer P-O bond distances, obtained by representing the nearest-neighbor P-O correlations by two Gaussian functions, are distinguished.network connectivity.For example, the structure of crystalline NAGP with x = 0 is based on a network of corner-sharing tetrahedral PO 4 and octahedral GeO 6 units in which each BO atom is shared between a tetrahedron and an octahedron.Each BO atom therefore has three nearest-neighbor O atoms within a tetrahedron and four nearest-neighbor O atoms within an octahedron such that nO O = nO BO = 7.The parameters predicted by the structural model for a variety of crystalline M 2 O-Al 2 O 3 -GeO 2 -P 2 O5 materials (Table a TTT = 873 K for 0.25 h.b TTT = 876 K for 0.25 h.c TTT = 876 K for 0.5 h. the

TABLE IX .
The parameters predicted by the structural model of Sec.V B using nO P = 4 and the reported values of nO Ge and nO Al for several crystalline M 2 O-Al 2 O 3 -GeO 2 -P 2 O 5 materials.N NBO / N BO / N NBO / a x = 0 or 0.127.b x = 0, 0.24, 0.40, 0.44, 0.47, or 0.48.and a decrease in N NBO /N Na [Figs.23(c)-23(e)]

TABLE XI .
27e Al-O and Ge-O coordination numbers for the NAGP glasses with x = 0.4 and x = 0.8 found from the27Al MAS NMR (TableVII) and diffraction (TableVIII) experiments.