Structural Features and Defect Equilibrium in Cubic PrBa1−xSrxFe2O6−δ

The structure, oxygen non-stoichiometry, and defect equilibrium in perovskite-type PrBa1−xSrxFe2O6−δ (x = 0, 0.25, 0.50) synthesized at 1350 °C were studied. For all compositions, X-ray diffraction testifies to the formation of a cubic structure (S.G. Pm3¯m), but an electron diffraction study reveals additional diffuse satellites around each Bragg spot, indicating the primary incommensurate modulation with wave vectors about ±0.43a*. The results were interpreted as a sign of the short-order in both A-cation and anion sublattices in the areas of a few nanometers in size, and of an intermediate state before the formation of an ordered superstructure. An increase in oxygen deficiency was found to promote the ordering, whereas partial substitution of barium by strontium caused the opposite effect. The oxygen content in oxides as a function of oxygen partial pressure and temperature was measured by coulometric titration, and the data were used for the modeling of defect equilibrium in oxides. The simulation results implied oxygen vacancy ordering in PrBa1−xSrxFe2O6−δ that is in agreement with the electron diffraction study. Besides oxidation and charge disproportionation reactions, the reactions of oxygen vacancy distribution between non-equivalent anion positions, and their trapping in clusters with Pr3+ ions were taken into account by the model. It was demonstrated that an increase in the strontium content in Pr0.5Ba0.5−xSrxFeO3−δ suppressed ordering of oxygen vacancies, increased the binding energy of oxygen ions in the oxides, and resulted in an increase in the concentration of p-type carriers.


Introduction
The perovskite-type ferrites R x A 1−x FeO 3−δ , where R is a rare-earth element and A is an alkaline earth one, are known to possess mixed oxygen-ion and electron conductivity (MIEC), which makes them promising functional materials for electrochemical devices for energy conversion and storage, such as solid oxide fuel cells (SOFCs), membrane reactors for oxygen separation from air and synthesis gas production [1][2][3][4][5][6].
Ion and electron conductivity of p-and n-type are known to be strongly affected by the oxygen content in the oxide (3−δ), temperature, type and fraction of the R cation [7][8][9][10]. Although in the most MIEC ferrites R x A 1−x FeO 3−δ lanthanum is used as R, materials based on perovskite-type ferrites with other R cations attracted much attention in recent years [6]. For example, the PrBaFe 2-x Co x O 5+δ oxides were successfully tested as cathode materials for intermediate-temperature SOFCs [11]. RBaFe 2 O 5+δ ferrites were reported to exhibit high stability under reducing atmosphere at elevated temperatures [12]. PrBaFe 2 O 5+δ was shown to demonstrate good electro-catalytic activity and chemical stability in the H 2 Materials 2022, 15, 4390 2 of 20 atmosphere [13]. This allows considering PrBaFe 2 O 5+δ based oxides as promising SOFCs anode materials [14][15][16].
A specific feature of Pr 0.5 Ba 0.5 FeO 3−δ , in comparison with La 0.5 Ba 0.5 FeO 3−δ , is a greater difference in the size of A-site residing cations, which should favor their ordering. Although a number of studies reported about Pr 0.5 Ba 0.5 FeO 3−δ (or PrBaFe 2 O 6−δ ) with a cubic perovskite structure [17,18], it was demonstrated that perfectly ordered double perovskite PrBaFe 2 O 5 with an orthorhombic (SG Pmmm) structure could be synthesized via a special route under reducing atmosphere and appropriate temperature conditions [19]. Moreover, a literature analysis demonstrates that simple synthesis in air can also result in the formation of PrBaFe 2 O 6−δ with different lattice symmetry. Besides the cubic structure [15,17,20], a tetragonal one with a P4/mmm space group [14,21], and orthorhombic one with Pmmm [22,23] are also reported for this composition prepared in air. A comparison of the synthesis conditions used in the above mentioned studies demonstrates that the oxides with low symmetry were obtained at relatively low temperatures-950 or 1000 • C [14,[21][22][23], whereas synthesis at 1330 or 1350 • C provides the formation of PrBaFe 2 O 6−δ with a cubic structure [15,17,20]. This feature shows that although Pr 3+ and Ba 2+ cations in ferrites tend to be ordered due to a considerable size difference, synthesis at high temperature ensures the formation of a substantially disordered state in the A sublattice. It is natural to expect that the degree of ordering can affect the characteristics of charge transport in PrBaFe 2 O 6−δ . This should be taken into account in order to achieve reproducible operational properties of materials in electrochemical devices.
This work is aimed at studying the defect formation in PrBa 1−x Sr x Fe 2 O 6−δ . The interest in strontium substitution is attributed to its reported favorable impact on the electrical conductivity under both reducing and oxidizing conditions [15]. In order to minimize the uncertainties caused by possible structural ordering, the synthesis of oxides was performed at 1350 • C, which was demonstrated to provide formation of the cubic PrBaFe 2 O 6−δ structure. Since ordering on a scale below the XRD sensitivity limit but affecting the thermodynamic and transport properties of materials cannot be excluded, electron diffraction was used for the local structure examination. In the results of the study, it was found that the PrBaFe 2 O 6−δ oxide consists of the ordered domains with the size of a few nanometers. Due to a strong disposition to oxygen coordination less than 12-fold, Pr 3+ ions in the ordered domains trap oxygen vacancies in steady clusters [Pr 3+ -V O3 ] that favors reduced oxygen content in oxide. The substitution of barium with strontium was found to promote the disordering of cation sublattice that in turn causes the decomposition of [Pr 3+ -V O3 ] clusters and release of oxygen vacancies, an increase in oxygen content and concomitant increase in p-type carrier concentration. All of these effects are expected to be favorable for the enhanced electrode performance.

Experimental Section
The PrBa 1−x Sr x Fe 2 O 6−δ (x = 0, 0.25, 0.50) ferrites were prepared by the glycine-nitrate combustion synthesis. Metal iron, praseodymium oxide, high purity barium and strontium carbonates taken in the necessary proportion were dissolved in nitric acid. Glycine was added to the obtained solution in a 50% excess to nitrates, after which the mixture was stirred and heated until self-ignition occurred. The combustion product was annealed at 900 • C for 10 h in order to remove organic and carbon residues. Then, the product was milled, pressed into discs 20 mm in diameter and about 2 mm thick, and sintered at 1350 • C for 10 h. The obtained samples were milled into powder to be used in the further study. In order to evaluate material stability under reducing conditions, the samples were held in a gas mixture of 89%Ar-1%CO-10%CO 2 at 950 • C (p O 2~1 0 −15 atm) for 5 h and then cooled down at a rate of 5 • /min for structural analysis.
The phase purity of the materials and their structural characteristics were examined by X-ray diffraction (XRD) on a Shimadzu XRD-7000 diffractometer with CuK α radiation. The analysis of XRD data was performed by the Rietveld refinement method using the GSAS-II software [24]. The local structure of PrBa 1−x Sr x Fe 2 O 6−δ was studied by electron diffraction with a JEM-200 CX microscope. Powder specimens of oxides were treated in an isobutyl alcohol bath with ultrasonic stirring, and then used to form sediments upon a carbon film. The deposited film was supported with a copper grid.
The oxygen content in oxides was measured in relation to temperature by a Setaram TG-92 thermal analyzer in air. The measurements were performed in a cooling mode at a rate of 1 • /min after preliminary sample exposure at 950 • C for 5 h. The oxygen content in PrBa 1−x Sr x Fe 2 O 6−δ as a function of oxygen partial pressure (p O 2 ) at different temperatures was measured by coulometric titration. Measurements at oxygen partial pressure above 10 −3 atm were carried out in the isothermal mode by a stepwise decrease in p O 2 and subsequent holding for sample equilibration with the gas phase. Below 10 −3 atm, a step temperature mode was used for measurements because it provides a faster equilibrium state achievement [25]. The experimental data were recorded only when a variation in the logarithm of the oxygen partial pressure over the sample became less than 0.01 per hour. This equilibrium criterion provided good reversibility of measurement results. More details related to the experimental setup and measurement conditions can be found elsewhere [26,27].

Model of a Point Defect Equilibrium
In the studied ranges of temperature and partial pressure of oxygen, the iron in PrBa 1−x Sr x Fe 2 O 6−δ can be in 2+ , 3+ and 4+ oxidation states. Therefore, the composition of a solid solution can be presented as Pr 3+ Ba 2+ 1 Taking into account the disposition of barium and praseodymium to ordering, it is assumed that irrespective of the results of the XRD analysis, the oxides can have an ordered double perovskite structure at a low scale, as shown in Figure 1. There are three oxygen positions in double perovskite: O1, O2 and O3 with different bond energy. As in the case of PrBaCo 2 O 6−δ , the O1 position is considered to be completely occupied with oxygen ions, unlike O3, where the probability of vacancy appearance is the highest. Besides, the formation of oxygen vacancies in the O2 position is also possible [28].
The oxygen partial pressure in the gas phase affects the oxygen content in the oxides via the oxidation reaction: where V O3 designates oxygen vacancy at the O3 site. The charge disproportionation reaction on iron ions can also contribute to defect formation: The redistribution of oxygen ions between the O2 and O3 sites is assumed to be possible: The oxygen partial pressure in the gas phase affects the oxygen content in the oxides via the oxidation reaction: where O3 V designates oxygen vacancy at the O3 site. The charge disproportionation reaction on iron ions can also contribute to defect formation: Fe Fe The redistribution of oxygen ions between the O2 and O3 sites is assumed to be possible: It should be noted that in contrast to La 3+ , Ba 2+ or Sr 2+ , which can have a 12-fold oxygen coordination necessary for the cubic perovskite structure, for Pr 3+ ion, the maximum coordination number is known to be 9 [29]. This should restrict the number of oxygen ions around praseodymium and cause the trapping of oxygen vacancies in clusters: The temperature dependences of the equilibrium constants of these reactions are governed by the well-known thermodynamic expression: where R is the gas constant, and ∆G 0 i , ∆H 0 i and ∆S 0 i are the standard Gibbs free energy, enthalpy and entropy changes for the defect formation reactions, respectively. Additionally, the requirement of site conservation in the iron sublattice and the electroneutrality condition should be taken into account: Preliminary calculations for the PrBaFe 2 O 6−δ demonstrated that the vacancy concentration in the O2 position does not exceed 10 −9 over the studied ranges of temperature and partial pressure of oxygen. Therefore, Equations (3) and (8) were excluded from further consideration for this composition.
The joint solution of Equations (1)-(11) allows one to combine the partial pressure of oxygen, the oxygen content in oxides, and the thermodynamic parameters of reactions (1)- (4) in the single equation omitted here due to complexity and space limitation.

Partial Molar Enthalpy and Entropy of Oxygen in Pr
The thermodynamic equilibrium between the oxide and the surrounding gas phase implies that the respective chemical potentials µ O and µ O 2 are equal: The chemical potential of oxygen in the gas phase is determined as follows: where µ o O 2 is the chemical potential of oxygen in the gas phase under standard conditions. Thus, the chemical potential of oxygen in oxides relative to the standard state in the gas phase ∆µ O can be expressed as: In turn, this chemical potential of oxygen in oxides is related to the respective partial molar enthalpy ∆h O and entropy ∆s O : The Gibbs-Helmholtz relations can be obtained from Equations (9) and (10):

Statistical Thermodynamic Modeling
The total Gibbs energy for the PrBa 1−x Sr x Fe 2 O 6−δ solid solution is determined, as follows: where G 0 is the standard Gibbs energy, µ 0 is the chemical potential of the respective component and S conf is the total configurational entropy of the solid solution.
The configurational entropy in the Stirling approximation is defined by the formula [30]: . The chemical potential of oxygen can be determined as: Using Equations (1)- (11), (20) and (21), one can obtain: The detailed derivation of Equations (22)- (25) and expressions for are omitted because of their complexity and space limitation.

Structure and Oxygen Content
The XRD patterns of the PrBa 1−x Sr x Fe 2 O 6−δ samples, as-prepared in air and subjected to reductive treatment at 950 • C and presented in Figure 2, indicate that all oxides are single-phase with a cubic perovskite-type structure (S.G. Pm3m). This demonstrates that the disordered state of the A sublattice provided by synthesis at high temperature is steady and cannot be notably modified by reductive treatment at 950 • C. The unit cell parameter of 3.9387(1) Å for the as-prepared Pr 0.5 Ba 0.5 FeO 3−δ is in reasonable agreement with 3.930 and 3.934 Å, which were reported earlier [17,20]. A decrease in the lattice parameter with increasing strontium content observed in Figure 3 for both as-prepared and reduced series of oxides is naturally attributed to the smaller size of Sr 2+ ions than that of Ba 2+ . The larger lattice parameter of the reduced oxides is explained by higher fractions of Fe 3+ (RCN6 = 0.645 Å ) and Fe 2+ (RCN6 = 0.75 Å ) ions with a larger radius than of Fe 4+ (RCN6 = 0.585 Å ) [29]. Some deviation of the lattice parameter dependence on the strontium content from linear behavior should be attributed to the non-trivial relation of 3−δ vs. x. A decrease in the lattice parameter with increasing strontium content observed in Figure 3 for both as-prepared and reduced series of oxides is naturally attributed to the smaller size of Sr 2+ ions than that of Ba 2+ . The larger lattice parameter of the reduced oxides is explained by higher fractions of Fe 3+ (R CN6 = 0.645 Å) and Fe 2+ (R CN6 = 0.75 Å) ions with a larger radius than of Fe 4+ (R CN6 = 0.585 Å) [29]. Some deviation of the lattice parameter dependence on the strontium content from linear behavior should be attributed to the non-trivial relation of 3−δ vs. x.  An electron-diffraction (ED) study was carried out to investigate the local structure features of PrBaFe2O6−δ and the strontium-substituted derivative. For all SAED (selected area electron diffraction) patterns, the observed Bragg reflections correspond to a primitive cubic structure (  (Figures 4 and 5a). The average wave vectors are equal to ± 0.43 * ( * is the unit cell vector of reciprocal cubic lattice) and hence, are incommensurate with the unit cell vectors of perovskite cubic structure. Some electron-diffraction patterns of individual crystals reveal sharper incommensurate satellites, with their wave vectors being slightly larger and close to 0.44 * . For the as-prepared PrBaFe2O6−δ samples, in the SAED patterns, the incommensurate satellites are weak and more blurred ( Figure  5b). In addition, for the [001] SAED pattern, the positions of these satellites are given by the wave vectors q  0.4 * . For as-prepared PrBa1−xSrxFe2O6−δ, the incommensurate satellites are barely distinguishable (Figure 5c).  (Figures 4 and 5a). The average wave vectors are equal to ±0.43a * (a * is the unit cell vector of reciprocal cubic lattice) and hence, are incommensurate with the unit cell vectors of perovskite cubic structure. Some electron-diffraction patterns of individual crystals reveal sharper incommensurate satellites, with their wave vectors being slightly larger and close to 0.44a * . For the as-prepared PrBaFe 2 O 6−δ samples, in the SAED patterns, the incommensurate satellites are weak and more blurred (Figure 5b). In addition, for the [001] SAED pattern, the positions of these satellites are given by the wave vectors q ≈ ±0.4a * . For as-prepared PrBa 1−x Sr x Fe 2 O 6−δ , the incommensurate satellites are barely distinguishable (Figure 5c).
Such incommensurate modulations are believed to be created by atomic displacement and/or an occupancy (composition) probability wave with wave vector q for some atomic sites in the lattice, and the modulated state can be intermediate before the formation of ordered superstructures with changes in the composition or temperature [31]. For a range of the A-sublattice-doped ferrites R x Sr 1−x FeO 3−δ (R is Y, Ce), which are known to be prone to the formation of ordered superstructures, similar incommensurate modulations have been observed in previous TEM studies [32][33][34][35]. The modulated states and ordered phases were found mainly in the reduced samples. In the RBaFe 2 O 5+w ferrites, reduced at high temperatures, a few ordered phases with varying oxygen nonstoichiometry (w) were also recently obtained [12,36]. The oxygen-poor BaPrFe 2 O 5.01 was reported to have an orthorhombic double-perovskite structure (S.G. Pmmm) at room temperature. For compositions with 0.1 < w < 0.48, the orthorhombic distortion is unstable, and a tetragonal superstructure P4/mmm is formed at temperatures above 1000 • C. At temperatures between 600 and 1000 • C and w ≥ 0.5 the presence of the Pr 2 Ba 2 Fe 4 O 11 orthorhombic structure, Cmmm was confirmed [12]. Therefore, the observed incommensurate modulations in the reduced PrBaFe 2 O 6−δ can be considered as a sign of an intermediate state related to the order-disorder transformations in ferrite. A number of our experimental findings support the above assumption. The intensity of incommensurate satellites substantially decreases in the as-prepared PrBaFe 2 O 6−δ sample. In addition, Sr-doping of the as-prepared PrBaFe 2 O 6−δ sample leads to the almost complete disappearance of these satellites.    Similar to R x Sr 1−x FeO 3−δ (R is Y, Ce), all PrBaFe 2 O 6−δ superstructures are characterized by simultaneous ordering in the oxygen-vacancy and A-metal sublattices. In the ordered Pr 2 Ba 2 Fe 4 O 11 phase, the ordering of both sublattices (Ba/Pr and oxygen/vacancy) is accompanied by atom displacements from the ideal perovskite-like positions in the <100> directions of the basic cubic lattice [12]. According to the "modulation wave approach" [37], composition and displacive modulation waves with wave vectors q can be suggested in both sublattices of the reduced PrBaFe 2 O 6−δ , with modulation waves being poorly correlated because only the first order (primary) but not the higher order (±nq) satellites are observed. Such a state with only primary incommensurate diffuse satellites was proposed to be called a disordered modulated structure, which does not exhibit a long-range order, only a short-range (local) order [38]. The correlation length, which is inversely proportional to the finite width (or smearing) of diffuse incommensurate peaks, can be estimated at about a few nanometers. Hence, the size of the region with the modulation wave is close to this value. This estimation is generally consistent with HRTEM data on the nonuniform nanosized structure of Y 0.25 Sr 0.75 FeO 3−δ powder demonstrating similar primary incommensurate diffuse peaks in the SAED pattern [32]. It can be observed in Figure 4 that the correlation length and ordering tendency decreases in the series of samples: reduced PrBaFe 2 O 6−δ , as-prepared PrBaFe 2 O 6−δ , and as-prepared PrBa 0.5 Sr 0.5 Fe 2 O 6−δ . A decrease in value of the wave vector q is also observed with a decrease in the correlation length.
Thus, the electron diffraction study has found that in the reduced PrBaFe 2 O 6−δ , a disordered modulated (or short-range) structure but not an ordered superstructure is formed. The short-range is characterized by primary incommensurate modulations with wave vectors q = ±0.43a * (a is perovskite-lattice vector), which are considered to be related to displacement and composition waves in both Ba/Pr and oxygen/vacancy sublattices of perovskite structure. Oxygen nonstoichiometry, as well as Sr-doping, are crucial parameters of incommensurate modulations, correlation length and ordering tendency in Pr 0.5 Ba 0.5−x Sr x FeO 3−δ . Figure 6 presents isothermal plots of the oxygen content in PrBa 1−x Sr x Fe 2 O 6−δ versus the oxygen partial pressure at different temperatures. All isotherms contain a wide plateau with an inflection point, which corresponds to the nominal oxidation state 3+ of iron according to the formula Pr 3+ Ba 2+ 1−x Sr 2+ x Fe 3+ 2 O 2− 5.5 . These features were used to adjust the experimental arrays p O 2 -T-(6−δ) to the reference value of oxygen content (6−δ) = 5.5. The cross-like symbols on the 950 • C isotherms in Figure 6 designate the oxygen content in the oxides in air. These values indicated by the same symbols in Figure 7 were used as references to calculate oxygen content data from thermogravimetric measurement results.
As can be observed in Figure 7, all compositions at room temperature are oxygen deficient. This deficiency increases with increasing barium content and is presumably attributed to the ordering of oxygen vacancies accompanied by the formation of stable Pr 3+ − V O3 clusters. The maximum oxygen content, which PrBaFe 2 O 6−δ attains at slow cooling in air, is 5.78, which is in good agreement with the doubled value of~2.885 determined for Pr 0.5 Ba 0.5 FeO 3−δ in [14]. Another effect of strontium substitution, which

Defect Equilibrium Parameters and Thermodynamic Quantities
The experimental data on the oxygen content in PrBa1−xSrxFe2O6−δ, versus the oxygen partial pressure at different temperatures presented in Figure 6 by dots, were employed for the modeling of defect equilibrium in the oxides. The results of calculations presented in Figure 6 by solid lines correspond reasonably well to the experimental data. The obtained thermodynamic parameters of the defect formation reactions are summarized in Table 1. Table 1. Thermodynamic parameters of defect formation reactions in PrBa1−xSrxFe2O6−δ.
The enthalpy of the charge disproportionation reaction (2) is similar for all compositions of PrBa1−xSrxFe2O6−δ. In contrast, the enthalpy of the oxidation reaction (1) decreases with increasing strontium content in PrBa1−xSrxFe2O6−δ. The same trend in ΔHox variation was observed for La0.49Ba0.5−xSrxFeO3−δ [39]. A decrease in ΔHox upon an increase in strontium content indicates an increase in the binding energy of oxygen in the oxide, which to some degree can result from lattice shrinkage and the respective contraction of the Fe-O bond lengths. At the same time, a big difference in the values of ΔHox for La0.49Ba0.5FeO3−δ (−95 kJmol −1 ) and PrBaFe2O6−δ (−80 kJmol −1 ) ferrites clearly indicates the presence of other factors affecting the enthalpy of the oxidation reaction. The origin of this difference can be related to lower available oxygen coordination for Pr 3+ ions in

Defect Equilibrium Parameters and Thermodynamic Quantities
The experimental data on the oxygen content in PrBa 1−x Sr x Fe 2 O 6−δ , versus the oxygen partial pressure at different temperatures presented in Figure 6 by dots, were employed for the modeling of defect equilibrium in the oxides. The results of calculations presented in Figure 6 by solid lines correspond reasonably well to the experimental data. The obtained thermodynamic parameters of the defect formation reactions are summarized in Table 1. Table 1. Thermodynamic parameters of defect formation reactions in PrBa 1−x Sr x Fe 2 O 6−δ .
The enthalpy of the charge disproportionation reaction (2) is similar for all compositions of PrBa 1−x Sr x Fe 2 O 6−δ . In contrast, the enthalpy of the oxidation reaction (1) decreases with increasing strontium content in PrBa 1−x Sr x Fe 2 O 6−δ . The same trend in ∆H ox variation was observed for La 0.49 Ba 0.5−x Sr x FeO 3−δ [39]. A decrease in ∆H ox upon an increase in strontium content indicates an increase in the binding energy of oxygen in the oxide, which to some degree can result from lattice shrinkage and the respective contraction of the Fe-O bond lengths. At the same time, a big difference in the values of ∆H ox for La 0.49 Ba 0.5 FeO 3−δ (−95 kJmol −1 ) and PrBaFe 2 O 6−δ (−80 kJmol −1 ) ferrites clearly indicates the presence of other factors affecting the enthalpy of the oxidation reaction. The origin of this difference can be related to lower available oxygen coordination for Pr 3+ ions in comparison to La 3+ . It is easy to estimate that the average 9-fold oxygen coordination of Pr 3+ ions in PrBaFe 2 O 6−δ at the full 12-fold coordination of Ba 2+ ions corresponds to the oxygen content of 5.25, which is less than the experimental results in Figure 6. One can assume that excessive oxygen around the Pr 3+ ion is less favorable for this oxide and its incorporation results in a decrease in the binding energy of oxygen ions reflecting a higher ∆H o ox value, in comparison to La 0.49 Ba 0.5 FeO 3−δ with 12-fold oxygen coordination for La 3+ ions.
Another interesting trend that the data in Table 1 demonstrate is an increase in the ∆H c with an increase in strontium content in PrBa 1−x Sr x Fe 2 O 6−δ indicating that strontium substitution favors the release of V O3 oxygen vacancies from Pr 3+ − V O3 clusters.
In order to verify the adequacy of the used point defect model and the obtained thermodynamic parameters, the partial molar enthalpy and entropy of oxygen versus the oxygen content in the oxides are calculated and shown in Figures 8 and 9, respectively. The dots in the figures present the results obtained directly from the experimental data via Gibbs-Helmholtz relations (16) and (17), whereas the lines show the results of statistical thermodynamic calculations using Equations (23) and (24). comparison to La 3+ . It is easy to estimate that the average 9-fold oxygen coordination of Pr 3+ ions in PrBaFe2O6−δ at the full 12-fold coordination of Ba 2+ ions corresponds to the oxygen content of 5.25, which is less than the experimental results in Figure 6. One can assume that excessive oxygen around the Pr 3+ ion is less favorable for this oxide and its incorporation results in a decrease in the binding energy of oxygen ions reflecting a higher ox H   value, in comparison to La0.49Ba0.5FeO3−δ with 12-fold oxygen coordination for La 3+ ions.
Another interesting trend that the data in Table 1 demonstrate is an increase in the ΔHc with an increase in strontium content in PrBa1−xSrxFe2O6−δ indicating that strontium substitution favors the release of VO3 oxygen vacancies from In order to verify the adequacy of the used point defect model and the obtained thermodynamic parameters, the partial molar enthalpy and entropy of oxygen versus the oxygen content in the oxides are calculated and shown in Figure 8 and Figure 9, respectively. The dots in the figures present the results obtained directly from the experimental data via Gibbs-Helmholtz relations (16) and (17), whereas the lines show the results of statistical thermodynamic calculations using Equations (23) and (24).   According to Equation (23) ] O 2− O2 + 2Fe 3+ oxidation reaction. As for the partial molar entropy of oxygen, the behavior of this function, according to Equation (24), is conditioned mainly by the contribution of the configurational entropy. Since the partial molar enthalpy and entropy of oxygen are very sensitive to deviations in the oxygen content data, the general agreement of the experimental and calculated results indicates the correctness of the employed model that allows for considering the obtained defect concentrations reliable.
The concentrations of iron ions calculated with the obtained thermodynamic parameters of defect formation reactions are shown in Figure 10. It can be observed that in wide ranges of p O 2 the average oxidation state of iron ions is close to 3+, therefore the concentration of Fe 2+ and Fe 4+ is proportional to p -1/4 O 2 and p +1/4 O 2 , respectively. The dependence of these concentrations on the ratio of barium and strontium is determined by the values of ∆H o ox and ∆H o ox (2−3) . Comparison of the data in Table 1 and the curves in Figure 10 shows that the higher the values of ∆H o ox and ∆H o ox(2−3) , the higher the concentration of Fe 4+ and Fe 2+ ions. In general, strontium substitution favors increasing the p-type carrier concentration.
A significant impact of strontium content on the defect formation in PrBa 1−x Sr x Fe 2 O 6−δ is revealed in Figure 11, where oxygen vacancy distribution between available positions is shown. It can be observed that the main number of vacancies in all oxides are bound in clusters. A distinctive feature of the strontium-free composition is the absence of oxygen vacancies in the O2 positions and a very low concentration of free V O3 vacancies. Strontium substitution to the degree x = 0.25 results in two orders of magnitude increase in the concentration of free V O3 vacancies, and more importantly, in the appearance of oxygen vacancies in the O2 positions. At x = 0.5, the concentration of V O2 vacancies becomes higher than that of V O3 . All of these indicate that the substitution of barium by strontium results in a progressive disordering of the anion sublattice in PrBa 1−x Sr x Fe 2 O 6−δ that is in good agreement with the results of electron diffraction studies. The mechanism of strontium action can be assumed as follows. For small Pr 3+ ions, the 12-fold oxygen coordination required for a complete anion sublattice is unfavorable, but a separate Pr 3+ ion can acquire it under the force of the rest of the lattice. The ordering in A-sublattice caused by a big difference in size of Ba 2+ (R CN12 = 1.610 Å) and Pr 3+ (R CN9 = 1.179 Å) ions [29] collects Pr 3+ ions into layers and enhances their ability to lower coordination. This promotes the ordering of the anion sublattice via the formation of steady clusters [Pr 3+ -V O3 ], and decreases the overall oxygen content in the oxide. Partial substitution of Ba 2+ ions by smaller Sr 2+ (R CN12 = 1.440 Å) [29], weakens the tendency to cation ordering, dissipates Pr 3+ layers, destroys [Pr 3+ -V O3 ] clusters and results in the anion sublattice disordering ( Figure 11) and in an increase in the oxygen content in the oxide (Figure 7). The higher oxygen content makes a higher concentration of p-type carriers according to reaction (1), as can be observed in Figure 10.
An increase in the concentration of free vacancies should favor oxygen transport, therefore partial substitution of barium by strontium can be considered as a promising approach to improve the operational characteristics of PrBaFe 2 O 6−δ -based materials in electrochemical devices. A significant impact of strontium content on the defect formation in PrBa1−xSrxFe2O6−δ is revealed in Figure 11, where oxygen vacancy distribution between available positions is shown. It can be observed that the main number of vacancies in all oxides are bound in clusters. A distinctive feature of the strontium-free composition is the absence of oxygen vacancies in the O2 positions and a very low concentration of free VO3 vacancies. Strontium substitution to the degree x = 0.25 results in two orders of magnitude increase in the concentration of free VO3 vacancies, and more importantly, in the appearance of oxygen vacancies in the O2 positions. At x = 0.5, the concentration of VO2 vacancies becomes higher than that of VO3. All of these indicate that the substitution of barium by strontium results in a progressive disordering of the anion sublattice in PrBa1−xSrxFe2O6−δ that is in good agreement with the results of electron diffraction studies. The mechanism of strontium action can be assumed as follows. For small Pr 3+ ions, the 12-fold oxygen coordina- and in an increase in the oxygen content in the oxide (Figure 7). The higher oxygen content makes a higher concentration of p-type carriers according to reaction (1), as can be observed in Figure 10.
An increase in the concentration of free vacancies should favor oxygen transport, therefore partial substitution of barium by strontium can be considered as a promising approach to improve the operational characteristics of PrBaFe2O6−δ-based materials in electrochemical devices.

Conclusions
A series of PrBa 1−x Sr x Fe 2 O 6−δ perovskite-type oxides were synthesized at 1350 • C. All oxides, both as-prepared and reduced at 950 • C and p O 2~1 0 −15 atm were characterized by XRD to have a cubic structure, whereas an electron diffraction study detected primary incommensurate modulations with wave vectors q = ±0.43a * , which manifests a tendency to ordering in A-or/and anion sublattices.
The oxygen content in PrBa 1−x Sr x Fe 2 O 6−δ was measured versus oxygen partial pressure in the interval between 10 −15 and 0.5 atm at 750-950 • C by coulometric titration. The obtained data were used for the modeling of defect equilibrium in oxides. The point defect model assumed that the oxides consisted of chaotically oriented domains with an ordered double perovskite structure. Taking into account the reactions of oxygen vacancy distribution between non-equivalent positions, and their trapping in clusters with Pr 3+ ions, in addition to the usual oxidation and charge disproportionation reactions, allowed for obtaining an appropriate agreement with the experimental data.
The modeling results have demonstrated that the partial substitution of strontium for barium in PrBa 1−x Sr x Fe 2 O 6−δ leads to a progressive disordering of anion sublattice, an increase in the binding energy of oxygen in the oxides and an increase in the concentration of p-type carriers.
The correctness of the point defect model and the obtained thermodynamic parameters has been confirmed by a reasonably good agreement between the partial molar enthalpy and entropy of oxygen obtained from the experimental data and by statistical thermodynamic calculations.