Ionic and Thermal Transport in Na-Ion-Conducting Ceramic Electrolytes

We have studied the ionic and thermal transport properties along with the thermodynamic key properties of a Na-ion-conducting phosphate ceramic. The system Na1+xAlxTi2−x(PO4)3 (NATP) with x = 0.3 was taken as a NASICON-structured model system which is a candidate material for solid electrolytes in post-Li energy storage. The commercially available powder (NEI Coorp., USA) was consolidated using cold isostatic pressing before sintering. In order to compare NATP with the “classical” NASICON system, Na1+xZr2(SiO4)x(PO4)3−x (NaZSiP) was synthesized with compositions of x = 1.7 and x = 2, respectively, and characterized with regard to their ionic and thermal transport behavior. While ionic conductivity of the NaZSiP compositions was about more than two orders of magnitude higher than in NATP, the thermal conductivity of the NASICON compound showed an opposite behavior. The room temperature value was about a factor two higher in NATP compared to NaZSiP. While the thermal conductivity decreases with increasing temperature in NATP, it increases with increasing temperature in NaZSiP. However, the overall change of this thermal transport parameter over the measured temperature range from room temperature up to 800 °C appeared to be relatively small.


Introduction
Solid ionic conductors are crucial components [1][2][3] of all solid-state batteries since they are transporting ions in electrochemical cells from the cathode to the anode and vice versa. The main advantage compared to conventional cells [4] is the lack of a liquid organic electrolyte, which is often considered as a safety risk [5] due to their high chemical reactivity causing thermal run-away events at higher temperatures.
Since the last decade Li-ion cells containing liquid electrolytes are dominating the market of portable and stationary devices for electrochemical storage. However, up to now all solid-state cells are only available as prototypes on the laboratory scale. Furthermore, due to the expected increasing number of Li-ion cells for the growing market of electrical vehicles and stationary storage for intermittent renewable energy sources, limited Lithium resources and higher costs may call for alternative elements in future decades. A promising alternative at least for stationary storage could be Sodium [6,7], which is cost effective without the limitation in the Li resources and exhibits comparable physicochemical and electrochemical properties.
Phosphate-based ceramic electrolytes are promising candidates [8][9][10] for allsolid-state Na-ion cells. Specifically, the NASICON (Na Super Ionic Conductor) structured electrolytes can be tailored to required properties by doping and structural tuning. The NASICON framework is composed of metal oxide MO 6 -octahedra which are connected to SiO 4 -and PO 4 -tetrahedra, respectively. Within this 3D structure, the Na + ions are migrating by hopping from occupied to empty lattice sites along the c-axis. In NATP (space group R 3 c), the Na ions are located at two different sites with Na1 (Wyckoff position 18e), which is nearly empty, and Na2 (Wyckoff position 6b), which is fully occupied. In NaZSiP (space group C2/c), the Na ions are distributed over three different sites: Na1 (Wyckoff position 4d) has an occupation level of 80 %, Na2 (Wyckoff position 4e) is completely occupied, and Na3 (Wyckoff position 8f) is occupied to 60 %.
Although a large number of studies is published on the relationship between the ionic transport and the structural and morphological properties [11][12][13][14][15][16][17][18] in phosphate-based Na-ion conductors, up to now the number of publications [19][20][21][22] dealing with the thermophysical or thermodynamic properties is significantly lower. Up to now, to the authors best knowledge, investigations which studied the ionic conductivity along with the thermal conductivity are still missing. Therefore, a possible correlation between the ionic and the thermal transport remains an open question due to the lack of combined ionic and thermal transport data, respectively. Therefore, our research work within this study is focused on the determination of the thermophysical properties like the thermal conductivity and the specific heat capacity and also on the evaluation of possible interrelations between mechanisms of ionic and thermal transport as well as thermodynamic properties of ceramic Na + ion conductors.

Experimental
Commercially available NATP powder (NEI Coorp., USA) with a stoichiometry of Na 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 was compacted to a green state using Cold Isostatic Pressing (CIP) at pressures of about 300 MPa. After a detailed in situ sintering study in order to optimize the heat treatment during the consolidation, NATP samples were sintered at 900 °C with a holding time of 11 h which results in a density of 2.86 g·cm −3 , which is 96 % of the theoretical density. X-ray powder diffraction (XRD) confirmed the rhombohedral lattice structure (R3 c) of NATP. By chemical analysis it could be shown that the NATP is nearly phase pure with only a small amount of AlPO 4 as impurity phase, which was less than 2.5 wt %, located in the grain boundaries.
Since there is no commercially available source for Na 1+x Zr 2 (SiO 4 ) x (PO 4 ) 3−x (NaZSiP), a synthetization route had to established for the production of the starting powder. Since a simple solid-state reaction process is known to be problematic [13,23,24] achieving the target stoichiometry due to the high-temperature processes involved, the so-called solution-assisted solid-state reaction (SA-SSR) [25] was adapted for the production of the NaZSiP powder. For the consolidation of the ceramic powder, a field-assisted sintering technology (FAST) at a temperature of 900 °C and also CIP followed by pressureless sintering at 1250 °C for 5 h was applied. Structural analysis by XRD confirmed the NASICON structure which shows a monoclinic distortion, and also traces of a ZrO 2 phase could be found.
The ionic conductivity was determined on the sintered samples by impedance spectroscopy in the frequency range from 100 Hz to 10 MHz as a function of temperature from RT up to 200 °C. The cylindrical sample with a typical thickness of 1 mm and a diameter of 10-12 mm is mounted in a heated ceramic sample holder which is built in a shielded chamber. Before the measurements, Au contacts with a thickness of about 150 nm are sputtered onto both faces of the sample. The voltage and current are applied by Au-coated spring-loaded contact pins in a two-point arrangement. The temperature is controlled by a K-type thermocouple located close to the sample. The temperature stability of the system is better than 1 °C.
The impedance is measured using an excitation AC voltage with a small amplitude (V A = 100 mV) at a given frequency f. The voltage V can be expressed as follows [26]: In a linear system, the response current signal I(t) has an amplitude I A and is shifted in phase: The complex impedance can be calculated as the ratio of the input voltage V(t) and the measured current I(t): The impedance can be, therefore, expressed in terms of a magnitude or absolute value Z A =|Z| and a phase shift . Alternatively, it can be written as a sum of a real Z′ and an imaginary part Z″, respectively: The impedance spectrum is recorded by a precision LCR meter (LCR 8110, GW Instek) over the frequency range given above. The results are plotted in the complex plane (cole-cole plot) which can be used for a first estimate of the bulk resistance by an inspection of the diameter of the depressed semicircle. For a more accurate evaluation, the experimental data are fitted to an equivalent circuit model (ECM). The real part of the frequency-dependent conductivity σ′, which is an important parameter for the ion migration in solid electrolytes [27], is calculated using the following expression: with d and A the thickness and surface area of the sample, respectively.
For each material, three different samples are taken from the same synthesis. The recording of the impedance spectra was repeated three times for each temperature level. The standard uncertainties of the bulk resistance and ionic conductivity are evaluated according to the guidelines of GUM (Guide to the Expression of Uncertainty in Measurements) [28].
In addition to the ionic transport also the thermal transport properties were measured by a Laser Flash Apparatus (LFA 427, Netzsch GmbH) on the sintered samples, which enables the direct determination of the thermal diffusivity. This particular property can be converted to the thermal conductivity if the specific heat and the density of a sample are known. Furthermore, changes of the thermal diffusivity within a given material system can be directly related to a variation in the phonon mean free path [29], which is affected by scattering processes at grain and phase boundaries, respectively, but also by defects. This is important for solid electrolytes since the number and distribution of point defects i.e., vacancies are crucial for the transport of ions within a crystal structure.
In a Laser Flash measurement [30], a short laser pulse (t pulse = 1 ms) is directed onto the front face of a cylindrical sample and the temperature rise is recorded as a function of time on the rear face by an infrared detector. Under ideal conditions, the thermal diffusivity α can be calculated from the time-temperature history by the following: where d is the thickness of the sample and t 0.5 the so-called half max time, i.e., the time in which the temperature rise on the rear face reaches half of the maximum value. Since this simple model expression is only valid under adiabatic conditions, a large number of more sophisticated models [31] have been developed which are able to consider deviations from ideal conditions by fitting the complete temperature-time curve. In this work a data reduction algorithm is used which is able to include radiation heat losses and finite size effects within the so-called combined model [32].
For the measurement the sample is placed in heating chamber which can be operated under vacuum or inert gas in a temperature range from RT to 1500 °C. Before the measurement the faces of the sample are coated with thin graphite film with a typical thickness of 10 µm-20 µm. This is necessary since the ceramic samples are translucent for the infrared radiation of the pulsed laser (Nd:YAG, 1064 nm) and for the InSb detector which records the temperature rise on the rear face of the sample.
The temperature stability of the furnace system is better than 1.5 °C. At each given temperature level the flash measurements were repeated for six times. Three different sample pieces from the same synthesis were included in a complete measurement run. The data were analyzed according to the GUM guidelines [28] in order to evaluate the standard uncertainty.
The specific heat capacity of the densified, sintered NASCION materials was measured using the standard three-step procedure [33] in a differential scanning calorimeter DSC 404 (Netzsch GmbH, Germany) for the high-temperature range from 220 °C to 900 °C and DSC 204 (Netzsch GmbH, Germany) for the low-temperature range from 25 °C to 320 °C.
The principle of differential scanning calorimetry (DSC) is based on the highprecision measurement of the heat flow during a predefined temperature cycle. During a "temperature scan", the temperature difference between a sample and a reference crucible placed symmetrically on a disk-type sample holder surrounded by a furnace is recorded. The heat flow rate between the crucibles can be determined as a function of temperature and time. For the measurement of the specific heat capacity [34], three "temperature scans" with identical heating rates, start, and end temperatures, respectively, are necessary. In the first step, a baseline is recorded with two empty crucibles. In the second and third run, the heat flow rates are recorded with the sample and a standard reference material placed in the sample crucible. Using data of the heat flow rate of the three scans, the specific heat c p of the sample can be calculated as follows: where Ф B , Ф R , and Ф S are the heat flow rates of the baseline, reference, and sample run, respectively, m R and m S mass of the reference and the sample, and C R is the specific heat capacity of the reference material. Synthetic sapphire (NIST, SRM 720 [35]) is used as a reference material.
In this work three different samples from the same batch were used for the heat capacity measurements. A typical sample size is 5 mm in diameter, with a thickness less than 1 mm and a typical sample mass of about 40 mg. The bottom side of the sample is polished in order to achieve a good contact to the crucible. For every sample, the temperature scan was repeated for three times resulting in three heating and three cooling ramps. Between each heating or cooling ramp, an isothermal phase with a duration of 30 min was introduced in order to establish thermal equilibrium before the next ramp. During the measurement, the furnace chamber and the sample holder are purged with high-purity Argon gas. Al crucibles are used in the lowtemperature DSC204 and Pt crucibles for the high-temperature part in the DSC404. The data from the three samples and the three heating ramps are used to evaluate the standard uncertainty of the specific heat according to GUM [28].
The density of the sintered samples was measured by the Archimedes method. This is also called hydrostatic weighing since the weight of the sample is measured in air and also immersed in a liquid. From the both measurements, the density ρ of a sample can be calculated as follows: where ρ liq is the density of the liquid, m air and m liq are the masses measured in air and in the liquid, respectively. In this work, high-purity ethanol was used a liquid with a density of 0.789 g·cm −3 . For the determination of the masses, a high-precision balance (Mettler Toledo MS-XSR105) is used which is equipped with an additional density measurement kit (MS-DNY-54) which allows the weighing in air and in the liquid without any additional sample preparation or changes in the balance setup. A typical sample mass varies between 600 mg and 1.3 g. Three different samples from the same batch are used in the density measurement. Each measurement was repeated for three times. The standard uncertainty of the density values was evaluated according to the guidelines of GUM [28].

Results and Discussion
Typical results extracted from the impedance spectroscopy on NATP are shown in Fig. 1. In the complex plane Z″ vs. Z′ a single depressed semicircle is visible. The bulk resistance can be estimated by measuring the diameter of the semicircle. This estimate is then further on used as an initial value for the fit procedure based on an equivalent circuit model (ECM). The ECM used for the fitting of the experimental data is shown in the inset of Fig. 1 (right). It consists of a combination of a resistance R SE representing the bulk resistance of the electrolyte in parallel to a constant phase element CPE SE for the bulk capacitance and an additional CPE BLE for the blocking electrodes [26,27]. The ECM fit shown as a continuous line in Fig. 1 yields a value for the bulk resistance R SE = 1.5 ± 0.05 × 10 -6 Ω and for the conductivity σ′ = 2.7 ± 0.095 × 10 -7 S·cm −1 . The uncertainties were calculated according to the guidelines of GUM [28] with a coverage factor of k = 2. The experimental data of the measured ionic conductivity in Fig. 2 are plotted in a log scale vs 1000/T in order to determine the activation energy E a for the ionic transport. The lines are the results of a calculation of the ionic conductivity σ vs. T −1 using the Arrhenius equation: where σ 0 is the pre-exponential factor and E a the activation energy, which have been used as fitting parameters in Eq. 9 to fit the experimental data. The uncertainties of the measured ionic conductivity were estimated to be about ± 3.5 %. For the activation energy, a standard deviation of ± 5 % was derived from the fitting procedure.
The two NASICON systems differ strongly in the ionic conductivity although the values of the activation energy E a are comparable. The room temperature value of NATP is σ Ion = 2.7 × 10 -7 (Ω·cm) −1 being significantly below the level of 10 -4 (Ω·cm) −1 in NaZSiP. The ionic conductivity increases by more than three orders of magnitude from NATP to NaZSiP. However, the activation barrier for ion transport, i.e., the activation energy, is only lowered by a factor of two.
This large difference in the ionic conductivity comparing the NATP (Na 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 ) with Na 2.7 ZSiP and further to Na 3.0 ZSiP (i.e., Na 1+x Zr 2 (SiO 4 ) x (PO 4 ) 3−x ), x = 1.7-2.0) is not surprising at a first glance since the Na content rises from 1.3 over 2.7 to 3.0 per formula unit, which increases also the number of Na + charge carriers within the structure. Together with the increased Na concentration additional sites appear for the Na ions which are only partially occupied. This effect is important for the ion migration through the structure since a balanced configuration of occupied and vacant lattice sites is crucial for ionic conduction. In NASICON systems the distance between two Na sites is about 3 Å depending on the Fig. 2 Ionic conductivity of the NASICON systems as a function of the inverse temperature specific composition. The ion transport occurs by thermal activation which explains the increasing conductivity within a given system with increasing temperature. Our conductivity data can be compared with the literature values of NASICON systems. Mouahid et. al. [36] synthesized a series of NATP ranging from 0 < x < 0.9 using a solid-state reaction. For the lowest Na concentration of x = 0.4, they determined a conductivity of σ 25 °C = 5.6 × 10 -8 S·cm −1 and an activation energy of E a = 0.54 eV. At a composition with x = 0.9, they found a room temperature conductivity of 1.3 × 10 -7 S·cm −1 and an activation energy of 0.51 eV. While the ionic conductivity was found to be lower than our values the activation energy is significantly lower. However, this can be explained by a sintering process at lower temperatures compared to our heat treatment. Lower sintering temperatures may lead to higher porosities which reduces the ionic conductivity.
Nieto-Munez and co-workers [12] synthesized glass-ceramics of NATP (Na 1+x Al x Ti 2−x (PO 4 ) 3 with 0 < x < 1.4) by a melt-quenching route. The glass transition as well as crystallization temperature were determined in a DSC and later heat treated the material at the onset of the crystallization temperature which was determined to be around 700 °C. The room temperature values of the ionic conductivity are comparable with the corresponding conductivity data of this work while the activation energy is slightly lower with a value of 0.54 eV for the composition of Na 1.4 Al 0.4 Ti 1.6 (PO 4 ) 3 , (x = 0.4) which is closest to our studied composition with x = 0.3.
In the "classical" NASICON system Na 3 Zr 2 (SiO 4 ) 2 (PO 4 ), the landscape of possible ionic conductivity values seems to be much wider [11,13,23]. Depending on the exact composition, processing route, and the heat treatment, room temperature values ranging from 4 × 10 -4 S·cm −1 up to about 1 mS/cm are found, while the activation energy is typically around 0.3 eV. The conductivity data of this work are within this range.
The measured thermal diffusivity values of NaZSiP and NATP are shown in Fig. 3. Standard uncertainties were calculated with a coverage factor of k = 2. It is obvious from the diagram that the uncertainties are higher in NATP than in the NaZ-SiP materials. This effect can be often observed with increasing thermal diffusivity values which leads to a lower signal to noise ratio in the laser Flash signal [31] inducing a higher scatter of the measured data. However, the relative uncertainty is about ± 3 % for NaZSiP and NATP which is typical in a Laser Flash apparatus. The thermal diffusivity values decrease with increasing temperature following a relationship of ∝ T −p with p = 0.17 for NATP and p = 0.04 for both NaZSiP compositions. Typical values for technical ceramics are in the range 0.5 < p < 1.0 [29,37] depending on the density of defects. This very low value of the exponent p reflects the strong influence of the defect scattering on the thermal diffusivity. Although the variation of the thermal diffusivity is relatively small within the temperature range of RT up to 900 °C, its thermal behavior is opposite to the ionic conductivity which shows a strong increase between RT and 200 °C by about two orders of magnitude. Furthermore, the NATP system with the lowest ionic conductivity exhibits the highest thermal diffusivity values compared to NaZSiP over the whole temperature range. Since the thermal diffusivity is directly related to the phonon mean free path, this gives some insights into the structure property relationship of ionic conductors, specifically in Na-conducting NASICON systems. Obviously, the Na concentration and the vacancy distribution are exposing an opposite effect on the ionic and thermal transport, respectively. While both are crucial for a high ionic conductivity and a low activation energy, they are perturbing the phonon transport with an effect not only on the level or absolute value of the thermal diffusivity but also on its temperature variation.
The heat capacity data of the studied NASICON systems are shown in Fig. 4. The relative uncertainties (coverage factor k = 2) are about ± 3.5 % in the low-temperature region, but increases up to ± 4.5 % at temperatures above 700 °C. In NATP, the apparent heat capacity is continuously growing with increasing temperature up to about 750 °C. Above this temperature, a peak centered around 800 °C is visible, which can be related to phase transformations in NATP. The specific nature of these transformations has to be evaluated in further studies, since no comparable data can be found in the literature. However, it should be noted here that a similar behavior was also observed in other phosphate-based Li + and Na +-conducting systems [38,39], which could be at least partially explained by a rotational distortion of the MO 6 -octahedra. This effect causes a rotation of the PO 4 -tetrahedra within the flexible polyhedral framework. The degree of distortion and the induced rotation depends obviously on the Na + doping level and as well as on the occupation of the different Na sites, respectively.
The two NaZSiP compositions show differences in the specific heat data below 300 °C and above 700 °C, while in between both systems can be characterized by nearly the same specific heat values. The enthalpy excursion in the apparent specific heat in the temperature range 100 °C < T < 200 °C of Na 3.0 ZSiP is due to structural phase transformation [19,20,40] from the monoclinic C2/c to the rhombohedral R 3 c structure. This transformation cannot be observed in the Na 2.7 ZSiP composition since it stays within the monoclinic C2/c structure over the whole temperature range. However, at temperature above 750 °C, the specific heat capacity values as a function of temperature of Na 2.7 ZSiP are increasing with a higher slope compared to Na 3.0 ZSiP. This additional anharmonic contribution could be due to an increased thermal expansion or further vacancy formation.
Studies of the specific heat capacity in Li-conducting NASICON systems, namely, LAGP Li 1+x Al x Ge 2−x (PO 4 ) 3 [41] with different compositions showed similar high-temperature behavior. However, the LAGP exhibited the anharmonic contributions for all studied compositions.
The values of the thermal conductivity λ of NATP and the NaZSiP compositions were calculated using the data of the thermal diffusivity α, the heat capacity Cp, and the density σ, respectively, by the following equation: For the calculation with Eq. 10, the specific heat data Cp(T) were taken from the experimentally determined data ignoring the peaks due to the phase transitions, where extrapolated baseline values were used. The measured density values at room temperature are given in Table 1. The density was not corrected due to thermal (10) = ⋅ C p ⋅ expansion, since the volume expansion in NATP and NaZSiP [42,43] is less than 1 % at 900 °C. The results of the thermal conductivity as a function of the temperature from RT to 900 °C are shown in Fig. 5. The combined standard uncertainty of the thermal conductivity was calculated considering the uncertainties of the thermal diffusivity, the heat capacity, and the density, respectively. The temperature variation of λ(T) in NATP is very small. It varies only by 15 % over the whole temperature range. However, λ(T) decreases with increasing temperature, but temperature variation of thermal conductivity is even lower than that of the thermal diffusivity, which is due to the effect that the specific heat increases toward higher temperatures and compensates the stronger decrease of the thermal diffusivity. The physical effect behind this behavior is the strong phonon scattering by point defects in the NATP, which leads to a low and nearly temperature-independent phonon mean free path.
In the NaZSiP compositions, a different temperature-dependent behavior can be observed. Although, the variation of λ(T) over the investigated temperature range is also small, the values tend to increase with increasing temperature. While the variation of the thermal diffusivity with temperature is much smaller in NaZSiP than in NATP, the monotonic growth of the specific heat data with temperatures determines the behavior of the thermal conductivity and leads to increasing thermal conductivity values with increasing temperatures.
It should be noted here that a 1/T behavior of the thermal conductivity [36,37], i.e., an exponent close to -1, is typically expected in non-conducting ceramic materials. Therefore, it can be concluded that a small temperature exponent tending to zero can obviously be taken as a kind of signature for superionic conductors, which results from a high concentration of point defects, which is needed for the thermally activated transport of the ions. For the further discussion, it might be helpful to take a closer look on the relationship of the thermal conductivity to the phonon mean where C is the heat capacity per unit volume, v is the average phonon group velocity, and l Ph is the mean free path of the phonons between scattering events. It is obvious that the thermal diffusivity α can be also related to the phonon mean free path as follows: without an explicit dependence on the heat capacity as in expression 11. Since the average phonon group velocity is almost independent of the temperature, it is straightforward to assume that the temperature dependence of the thermal diffusivity is related to a slowly decreasing mean free path with increasing temperature, i.e., it follows: as already discussed for the temperature dependence of the thermal diffusivity.
For the calculation of the mean free path using Eq. 11 the sound velocity v was calculated from density and elastic modulus data [45,46] which was measured by an indentation technique. The calculated values of sound velocity and phonon mean free path l Ph at room temperature are also shown in Table 1 together with data of the density and the elastic modulus. Under the assumption that the sound velocity is nearly independent of the temperature and the phonon mean free path (Eq. 13) follows the temperature dependence of the thermal diffusivity (Eq. 12), the thermal conductivity can be modeled according the Debye kinetic theory [44] given in Eq. 11. The results of this calculation are shown by the straight lines in Fig. 4.
It appears to be very interesting that the estimated room temperature values of the phonon mean free path derived from the thermal conductivity data can be related at least approximately to characteristic atomic distances within the NASICON framework. The Na-Na atomic distances in NATP are varying between 3.2 Å and 4.6 Å [35], respectively, and it is important to note that completely filled Na sites are surrounded by only partially filled sites. Therefore, an estimated phonon mean free path of 4.6 Å reflects the fact that the phonon propagation and consequently the thermal transport is also affected by this particular defect arrangement. The same effects can be observed in NaZSiP systems, in which the shortest hopping distance [23] between two Na sites varies between 3.5 Å and 3.7 Å. However, the Na-O distance [15] are showing variations from 2.4 Å to 3.1 Å. The estimations for the phonon mean free path is about 3.0 Å in Na 3.0 ZSiP and 2.6 Å in Na 2.7 ZSiP. Also, within this system the phonon mean free path seems to reflect structural details of the NASI-CON framework.
Although a direct comparison with other literature data of the NATP or the NaZSiP systems is not possible due the fact that thermal diffusivity or thermal conductivity data are not available up to now, the data of this work can be related to thermal transport studies in Na-or Li-ion-conducting phosphate-based systems. Suleiman and Lunden [47] studied the thermal conductivity and diffusivity in the ion conducting Na 3 PO 4 . They found low thermal conductivity values comparable with the results of this work. The thermal conductivity data between room temperature and 530 °C vary only slowly but give evidence of phase transitions in Na 3 PO 4 . However, since there are no data for the ionic conductivity are given, ionic and thermal transport could not be related. Petkov and Asabina [48] studied the thermophysical properties of NZP ceramics with a NaZr 2 (PO 4 ) 3 type structure. Within the temperature range from RT to 400 °C they found monotonical increasing values of the thermal conductivity for NaZr 2 (PO 4 ) 3 and Na 5 Zr(PO 4 ) 3 starting from 0.6 W·mK −1 and 0.4 W·mK −1 , respectively. This low level of the thermal conductivity values may be due to the relatively high amount (25 % to 30 %) of porosity in the samples. However, the ionic conductivity was not measured. The thermal conductivity and diffusivity, respectively, was studied [49] in the Li + conducting Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 analog system to NATP. The ionic conducting ceramic electrolyte was consolidated by a field-assisted sintering technique. The thermal transport parameter showed a strong dependence on the sintering temperature. Within the range of sintering temperatures from 800 °C to 1000 °C the values increased by a factor of more than two which could be related to growing grain sizes and increasing density. An investigation of the thermophysical properties [41] of Li + conducting Li 1+x Al x Ge 2−x (PO 4 ) 3 with x = 0.3, 0.5 and 0.7, respectively, came to similar results as in this work for the Na + conducting NaZSiP system. While the ionic conductivity increased with the Li + doping level, the thermal conductivity decreases accordingly.

Conclusion
In this work, we have studied the ionic and thermal transport properties of two different Na + ion conducting NASICON-structured systems. The thermal diffusivity and the thermal conductivity, respectively, were measured for the first time within these ion conductors together with the ionic conductivity. The ionic conductivity in NaZSiP is more than two orders of magnitude higher compared to NATP. However, the values of the thermal diffusivity and the thermal conductivity of NATP are about two times higher than in NaZSiP. Obviously, the defect concentration and the vacancy arrangement seem to play an important role for both the ionic and thermal transport, respectively. While a balanced vacancy configuration, i.e., a framework of partially occupied lattice sites is crucial for the migration of the ions, it has the opposite effect on the thermal transport since increasing defect concentrations are leading to a stronger scattering of phonons and reducing the thermal conductivity.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.