SciELO - Scientific Electronic Library Online

 
vol.42 issue164Preliminary study of anti-insuline IgG antibodies and immune complexes in Colombian patients with type 2 diabetesBioeconomy: A sustainable future author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista de la Academia Colombiana de Ciencias Exactas, Físicas y Naturales

Print version ISSN 0370-3908

Rev. acad. colomb. cienc. exact. fis. nat. vol.42 no.164 Bogotá JulySept. 2018

https://doi.org/10.18257/raccefyn.653 

Physical Sciences

Thermodynamic properties, electronic and crystallographic structure, and magnetic response of the Sr 2 HoNbO 6 material

Propiedades termodinámicas, estructuras electrónica y cristalográfica y respuesta magnética del material Sr 2 HoNbO 6

Críspulo E. Deluque Toro 1  

Ariday S. Mosquera Polo 1  

Jorge I. Villa Hernández 2  

David A. Landínez Téllez 2  

Jairo Roa-Rojas 2   * 

1 Grupo de Nuevos Materiales, Facultad de Ingeniería, Universidad del Magdalena, Santa Marta, Colombia

2 Grupo de Física de Nuevos Materiales, Departamento de Física, Universidad Nacional de Colombia, Bogotá, Colombia


Abstract

In this work we used the Wien2k code, within the framework of the Kohn-Sham Density Functional Theory (DFT), applying the Full-Potential Linearized Augmented Plane Wave method (FP-LAPW) and adopting the Generalized Gradient approximation (GGA) for the exchange-correlation energy due to Perdew, Burke, and Ernzerhof, as well as the Local Density approximation (LDA) for the calculation of the Density of States and band structure of the Sr 2 HoNbO 6 double perovskite. For calculations, we considered the Fmm (#225) space group, which was experimentally obtained from X-ray diffraction measurements and Rietveld refinement. The experimental lattice parameter was 8.018 Å, which is 99.2% in agreement with the theoretical prediction from the minimization of energy through the Murnaghan state equation. From the measurements of magnetic susceptibility as a function of temperature and the adjustment with the Curie law, we obtained a value for the effective magnetic moment of 10.01 μB, which is close to the theoretical expected from Hund’s Rule (10.60 μB). An energy gap of 3.3 eV between the valence band and the conduction band revealed the insulator character of the Sr 2 HoNbO 6 complex perovskite for the spin up configuration, but a semiconductor feature was observed for the spin down polarization, with an energy gap of 0.77 eV. The thermodynamic properties were calculated from the state equation by using the Debye quasiharmonic model. A specific heat behavior of CV≈CP was found at temperatures below T = 500 K, with Dulong-Petit limit values doubling those reported for perovskite materials.

Key words: Perovskite material; electronic structure; crystalline structure; thermodynamic properties.

Resumen

En el presente trabajo se utilizó el código Wien2k, en el marco de la teoría del funcional de la densidad de Kohn- Sham, aplicando el método de ondas planas aumentadas y linealizadas (full-potential linearized augmented plane wave, FP-LAPW) y adoptando la aproximación de gradiente generalizado (GGA) de Perdew, Burke y Ernzerhof para la energía de intercambio y correlación, así como la aproximación de densidad local (local density approximation, LDA) para el cálculo de la densidad de estados y la estructura de bandas de la perovskita doble Sr 2 HoNbO 6 . Para los cálculos se consideró el grupo Fmm (#225), experimentalmente obtenido a partir de mediciones de difracción de rayos X y del método de refinamiento de Rietveld. El parámetro de red experimental fue de 8.018 Å, el cual concuerda en un 99,2 % con las predicciones teóricas efectuadas a partir de la minimización de la energía mediante la ecuación de estado de Murnaghan. A partir de mediciones de susceptibilidad magnética en función de la temperatura y del ajuste con la ley de Curie, se obtuvo el valor del momento magnético efectivo 10,01 μB. Este valor es muy cercano del esperado teóricamente a partir de las reglas de Hund (10,60 μB). La brecha de energía determinada entre las bandas de valencia y de conducción fue de 3,3 eV, lo que revela el carácter aislante de la perovskita compleja Sr 2 HoNbO 6 para la configuración de espín hacia arriba, en tanto que se observó el carácter semiconductor para la polarización de espín hacia abajo, con una brecha de energía de 0,77 eV. Las propiedades termodinámicas se calcularon a partir de la ecuación de estado usando el modelo cuasi-armónico de Debye. Un comportamiento del calor específico, con CV≈CP, se encontró a temperaturas inferiores a T = 500 K, con valores del límite de Dulong-Petit que doblaban los que se han reportado para materiales del tipo de la perovskita.

Palabras clave: material tipo perovskita; estructura electrónica; estructura cristalina; propiedades termodinámicas.

Introduction

Perovskite materials (SHNO) represent a great percentage of systems which are actually investigated by the physics and chemistry of solids. In general, perovskites is represented by the ideal formula ABX3, where A is an alkaline earth element, B represents a transition metal or rare earth element and X, usually is the oxygen (Hazen, 1988). The modifications of atomic radii of A and B introduce structural distortions and new crystalline phases while inclusions of rare earth elements give the possibility to produce materials with exotic electric and magnetic properties (Kruth & West, 2001;Llamosa, et al., 2009). Partial substitutions of the A and B cations give rise to complex materials such as the double perovskites with generic formula A2BB’O6 (Triana, et al. , 2012) and other more complex perovskites (Moreno, et al., 2008). Its chemical configuration offers multiple chances to combine different elements, as well as the possibility to synthesize new materials involving a great variety of physical properties. Depending on magnetic and electric characteristics of B and B’ it is relatively easy to create new perovskite systems with half-metallic properties (Erten, et al., 2011), magnetoelectric response (Molegraaf, et al., 2009) or magnetic ordering (Lee & Marianetti, 2018), which offer promissory perspectives in the new field of spintronics technology (Pilo, et al., 2017). In this paper, we propose the synthesis and characterization of the Sr 2 HoNbO 6 ceramic material. We describe the crystalline structure of this double perovskite and perform morphological, magnetic, and electronic analyses. Furthermore, we present the results of our measurements of the magnetic response as a function of temperature. Moreover, considering that in recent years the density functional theory (DFT) has become a good tool to study electronic properties in perovskite-like material (Bonilla, et al., 2008;Landínez-Téllez, et al., 2013), we conducted a study of the electronic properties of these materials in order to establish the type of hybridization between the orbitals of HoO 6 and NbO 6 octahedra present in the structure. This calculation is the greatest motivation of this work and the experimental corroboration of the theoretical predictions with regard to the energy gap.

Experiments

Samples of SHNO were prepared by the solid-state reaction technique. Stoichiometric ratios of precursor powders SrCO3, Ho2O3, and Nb2O5 (purity 99,90%) were finely ground and thoroughly mixed. The precursor powder was calcined at 1100 °C for 148 h in ambient atmosphere. Then, the calcined mixture was again crushed, finely ground and pressed at 6 ton/cm2 pressures to form a disc (10 mm diameter, 2 mm thickness). This disc was sintered at 1300 °C for 17 h and furnace-cooled to room temperature. The structural characterization was carried out by analyzing the X-ray diffraction (XRD) data obtained from a Panalytical® X-Pert PRO MPD diffractometer (CuKα=1,540598 Å). Refinements of the experimental data were performed through the GSAS code. The XRD pattern of the sample was recorded for crystallographic phase characterization. The morphological characterization of samples was systematically effectuated from scanning electron microscopy (SEM), by using VEGA 3 equipment. The field cooling measurement of the magnetic susceptibility as a function of temperature was carried out with an MPMS® Quantum Design SQUID.

Theoretical calculations. We used the Wien2k program (Blaha, et al., 2001) within the framework of the Kohn-Sham Density Functional Theory (DFT) (Hohenberg & Khon, 1964;Khon & Sham, 1965), applied the Full-Potential Linear Augmented Plane Wave method (FP-LAPW) and adopted the Generalized Gradient (GGA) approximation for the exchange-correlation energy (Perdew, et al., 1996). The self-consistent process in the numeric package Wien2k allows the total energy calculation for constituent elements and alloys. Taking the experimental unit cell data as input, all the structures studied in this work were fully relaxed with respect to their lattice parameters and the internal degrees of freedom compatible with the space group symmetry of the crystal structure. The resulting energies versus volume functions were fitted to the equation of state by Murnaghan (Murnaghan, 1944) in order to obtain de minimum energy value, the bulk modulus, its pressure derivative and the equilibrium lattice parameters and associated volume. The muffin-tin radius of the elements was 2,5; 2,19; 1,8 and 1,6 (in u.a.) Sr, Ho, Nb and O, respectively, angular momentum up to l = 10 inside the muffin-tin sphere, a maximum vector in the reciprocal space of Gmáx =12.0, RMKAX=8, and a mesh of 1500 points in the first Brillouin zone (equivalent to a maximum of 56 k points in the irreducible Brillouin zone). Finally, the convergence criteria for the selfconsistent calculation were 0.0001 Ry for the total energies, 0.0005 u.a. in the charge and 1.0 mRy/u.a. in the internal forces. Spin polarization was included in the calculations. The thermal expansion of the crystal by the dependence of volume and the temperature on equilibrium, the Debye temperature, the specific heat at constant volume and the pressure, as well as the Grüneisen parameter as a function of temperature, were predicted from the Debye quasi-harmonic model (Blanco, et al., 2004; Deluque-Toro, et al., 2018).

Results and discussion

Crystalline structure. Figure 1a shows the XRD pattern for the SHNO material, which consists of strong peaks characteristic of a primitive cubic perovskite plus a few weak line reflections arising from the super-lattice resulting in the crystallization of this material in a cubic double perovskite Fmm space group (#225) (Howard, et al., 2003; Madueño, et al., 2006). No evidence of distortions from the cubic symmetry was observed in the XRD pattern. The basic perovskite composition is ABO3, where A is a large ion suitable for the 12-coordinated cube-octahedral sites and B is a smaller ion suitable for the 6-coordinated octahedral site (Lufaso, et al., 2006). Cubic complex perovskite with mixed species on a site (particularly the B site) may be represented by multiples of this formula unit and a larger unit cell, e.g., A2BB’O6 (Corredor, et al., 2011). Thus, in the SHNO composition, the Sr+2 with the largest ionic radius (1,13 Å) occupies position A of the cubic complex perovskite, Ho+3 (ionic radius 0,97 Å) and Nb+5 (ionic radius 0,70 Å) cations occupy the B and B’ positions. The ionic radii were calculated by using the SPuDs software (Triana, et al., 2011). Due to the ordering of B and B’ on the octahedral site of the ABO3 unit cell, there is a doubling in the lattice parameter on the basic cubic-perovskite unit-cell. Thus, the whole XRD pattern of SHNO can be indexed in the A2BB’O6 cubic structure with the cell edge a=2ap, where ap represents the cell lattice of the cubic perovskite. The presence of the superstructure reflection lines (311), (531), (533) and (551) in the XRD pattern of SHNO is the signature of an ordered cubic complex perovskite structure. In a substitutional solid solution A2BB’O6, there is a random arrangement of B and B’ on equivalent lattice positions in the crystal structure. Upon stable heat treatment, the random solid solution rearranges into a structure in which B and B’ occupy the same set of positions but in a regular way; such a structure is described as a superstructure (Howard, et al., 2003).

Figure 1 a) XRD pattern for the complex perovskite Sr2HoNbO6. Symbols represent experimental diffraction data and the base line is the difference between experimental and simulated patterns (continuous line). B) Structure of the Sr2HoNbO6 cubic complex perovskite constructed from the Rietveld refinement of XRD experimental data. 

In the superstructure, the position occupied by B and B’ are no longer equivalent and this feature is exhibited in the XRD pattern of the material by the presence of superstructure reflection lines. For a double cubic perovskite of the formula A2BB’O6 the intensity, in particular, the (111) superstructure reflection is proportional to the difference in the scattering power of the B and B’ atoms, when all the atoms are situated in the ideal position (Howard, et al., 2003). A disordered arrangement of B and B’ should result in zero intensity. Therefore, Ho+3 and Nb+5 cation ordering in SHNO in B and B’ positions is clearly distinguished by the presence of the significant intensity of (311), (531), (533) and (551) superstructural reflection lines. The lattice parameter of SHNO, obtained from the XRD data, was a=8,018(0) Å.

In the study and description of crystallographic structures, the concept of “positions” is fundamental (Wondratschek, 2006). The so-called Wyckoff positions for this structure were experimentally determined to be Sr2+=8c, Ho3+=4a, Nb5+=4b, and O=24e. In this notation, a, b, c, and e are the Wyckoff letters that determine all the points x for which the site-symmetry groups are conjugate subgroups of the Fmm (Deluque-Toro, et al., 2014). These letters constitute just a coding frame for the Wyckoff positions starting with “a” at the bottom position and continuing in alphabetical order (Wondratschek, 2006). The number of equivalent points per unit cell, which accompanies the letter of Wyckoff, is known as the multiplicity of the Wyckoff position and is related to the possible coordinates. In the case of the Fmm (#225) space group, the Wyckoff letter “a” has multiplicity value of 4 and the coordinate is (0,0,0). In the same way, “b” coordinate is (½,0,0) and its multiplicity is 4. On the other hand, “c” has a multiplicity of 8 and its coordinate is (¼,¼,¼). Finally, “e” has a multiplicity of 24, with coordinate (0,2678, 0, 0).

The tolerance factor calculated from the experimental data is τ=0.9232, which is in accordance with expected regime values for the cubic perovskite structures (Triana, et al. , 2015). From the refinement analysis, the octahedral distribution in the structure of the Sr 2 HoNbO 6 material was constructed as shown in Figure 1b where the structure of the SHNO double perovskite constructed from the data of Rietveld refinement is shown. As expected from the discussion above, the cationic ordering of the Ho3+ and Nb5+ in the B and B’ sites of the cubic double perovskite can be observed clearly. On the other hand, no tilting or rotation of the HoO6 and NbO 6 octahedra is observed in the perfect cubic structure.

Microgranular morphology. The Secondary Electrons Image in Figure 2 for the SHNO reveals a qualitative approximation to the surface microstructure.

Figure 2 Surface microstructure of the SHNO material 

Strongly coupled particles are observed in the surface morphology, which is a characteristic of good sinterization in the polycrystalline sample. The picture evidences the formation of micrometric clusters of nanometric grains. It is clear that the grains could end up having a few tens of nm, forming groups with the appearance of clusters of several sizes in μm. The particle size D was calculated using the Scherrer equation, which is based on the widening of the maximum intensity peak of the XRD pattern due to the particle size exclusively, and defined by the ratio D=Kλ/β(θ)cosθ, where K is the Scherrer constant given by , λ=1,540598 Å is the radiation length and β(θ)=Δ(2θ) represents the halfwidening of the diffraction profile. The calculated particle size for the SHNO was D=83,2 nm. As D depends on the nucleation and growth rate, then, the value obtained for the SHNO material is attributed to the decomposition reaction toward oxides, which may be responsible of the refinement of the particle size D. Given that the solid-state reaction destroys the original morphologies of the precursor oxides, fine particles with similar size distributions can be obtained. These fine particles might be combined once and again to form agglomerates with size distributions as shown in Figure 2. Since the increase of the annealing temperature or the sintering time increases the particle size D of the material, the grain morphology changes from the common plate-like form observed in ceramics to a spherical-like form.

Magnetic response. The magnetic properties of SHNO have been investigated by measuring the DC magnetic susceptibility in the temperature range from 50 to 300 K and on the application of a 500 Oe magnetic field intensity. Figure 3 shows the temperature dependence of the inverse of the DC magnetic susceptibility as a function of temperature for SHNO.

Figure 3 Inverse DC magnetic susceptibility as a function of temperature for SHNO (data points) and Curie fitting of the paramagnetic characteristic (solid line) 

The magnetic susceptibility data in this figure can be well fitted with the Curie law: χ = χ0 + C/T, where is the Curie constant, N is Avogadros’s number, μeff is the effective magnetic moment ( μeff = Peff μB ), represents the effective Bohr magneton number, μB is the Bohr magneton, kB is the Boltzmann constant, and χ 0 is the temperature independent susceptibility term. The value for this term is χ 0 = 0.00945 emu/mol.Oe. The Curie constant, estimated from the fitting in Figure 3, is C = 12.55892 emu.K/ mol.Oe. The effective magnetic moment due to the Ho3+ ion, calculated from the Curie constant, is μeff = 10.01 μB . This value corresponds 94.4% with the expected Peff = g[J(J + 1)]1⁄2 for an isolated Ho+3 cation with a 4f115d06s2 configuration, calculated by Hund’s rule, where g represents the Landé factor and J=S+L, S and L being the spin and orbital angular momenta. This difference may be attributed to the crystal field effects of the trivalent Ho+3 cation, which explain the magnetic susceptibility in SHNO.

Electronic structure. In order to obtain the most accurate results for the theoretical calculations of the density of states, we determined the optimal lattice parameters corresponding to the minimal energy value. Figure 4a shows the energy as a function of volume and the respective fitting to the Murnaghan state equation (Murnaghan, 1944). Each one of the points constitutes an individual calculation and the line corresponds to a fitting with the Murnaghan state equation, carried out by using the concept of the least-square fitting method (Murnaghan, 1944). The minimum of the energy was obtained for -46521,912456 Ry, the equilibrium volume was 529,147244 Å3, and the volume modulus, Bo = 137,6 GPa. The theoretical value of the lattice parameter was 8,088 Å, which is 99,2% in agreement with the experimental result.

Figure 4 a) Minimization of the total energy as a function of volume for Sr2HoNbO6 considering the Fmm space group. b) Band structure calculated for the Sr2HoNbO6 by considering electronic polarization 

Figure 4b shows the band structure for both spin up and spin down configurations. E=0 eV corresponds to the Fermi level. It can be seen that these materials evidence insulator behavior with gap energy through the Fermi level from -1,3 eV up to 2,1 eV (width: 3,3 eV) for the spin-up orientation. On the other hand, for the spin-down configuration, this energy gap remains with a particular characteristic: Two isolated bands appear for E=0 eV and E=1,1 eV. In order to clarify the electronic structure of this material close to the Fermi level, we calculated the density of the states taking into account the angular momenta and spin polarization. As exemplified in Figure 5a, there is no substantial contribution of the Sr2+ orbitals in any of the bands, the valence or the conduction near the Fermi level. Likewise, their orbitals in the conduction band occur above E = 2,5 eV and reveal very low density values. It is clear from the results that the greatest contribution of the electronic states in the valence band is due to the 2p-O orbitals close to de Fermi level for both spin-up and spin-down polarizations. These states correspond to the bands located in the regime of energies between -5,4 eV and -1,3 eV and they define the value of the energy gap in the valence band. Moreover, as observed in Figure 5b, it is remarkable that the 4d-Nb orbitals are the biggest contributors to the states corresponding to the valence band above E=2,5 eV for the spin-up and spin-down configurations. As expected, for the characteristic of the nonmagnetic material, the 4d-Nb5+ orbitals remain in a low spin configuration, so no deg states in the conduction band next to the Fermi level nor any of the available states in this regime are of the dt2g type. These states delimit the gap energy in the conduction band and their distribution is the reason why the effective magnetic moment is due solely to the contribution of the electrons 4f-Ho3+.

Figure 5 a) Minimization of the total energy as a function of volume for Sr2HoNbO6 considering the Fmm space group. B) Band structure calculated for the Sr2HoNbO6 by considering electronic polarization 

Finally, as shown in Figure 5c, the greatest contribution of states to the valence band was attributed to the 4f-Ho orbitals for the spin-down polarization, while there were no states available in this band for the spin-up orientation. An interesting and unusual feature is observed near the Fermi level for the 4f-Ho orbitals in the spin-down configuration; just on the Fermi level, between -0.2 eV and 0.1 eV, and above the Fermi level, between 0.9 eV and 1.3 eV, there are two peaks with high densities of states and short energies amid the energy gap. The magnetic moment obtained from the asymmetry of the calculated spin-up and spin-down density of states around the Fermi level was calculated to be 4,0 μB. This value corresponds to 40% when compared to that expected from the Hund’s rules for the isolated Ho3+ cation (Kittel, 1997). This circumstance is attributed to the exchange and correlation potential applied, which do not contain corrections due to the existence of 4d-Nb orbitals. Our results are consistent with a recent report of GGA calculations for this material, with exactly the same magnetic moment value of 4,0 μB, attributed mostly to the contribution of the 4f-Ho3+ orbitals and an insulating behavior determined from calculations of the band structure and the density of states (Khandy & Gupta, 2018).

Thermodynamical properties. We analyzed the effects of temperature and pressure on the thermodynamic properties of the SHNO material from the state equation under the considerations of the quasi-harmonic approximation of the Debye model as presented below. Figure 6 shows the results of specific heat at constant volume, CV (a), and at constant pressures, CP (b), as functions of temperature. As can be seen in the figures, the temperature was varied between 0 K and 1500 K for eight applied pressure values, from 0 to 35 Gpa. Figure 7 shows that below T=500 K for all applied pressure values, CV≈CP. In Figure 6(a) the trend of specific heat towards the Dulong-Petit limit, which is the specific heat value independent of temperature, is clearer than in Figure 6(b). As the temperature increases from this Dulong-Petit limit value, each of the atoms in the material absorbs the same amount of energy proportional to this temperature increase.

Figure 6 Specific heat Cv (a) and Cp (b) calculated by using the quasi-harmonic Debye model for the SHNO material from the Murnaghan state equation. 

Figure 7 Debye temperature (a), thermal expansion coefficient (b), and Grüneisen coefficient (c) for the SHNO calculated by the application of the quasi-harmonic Debye model 

In Figure 7 we present the results of the calculations of the Debye temperature (a), the thermal expansion (b), and the Grüneisen coefficient (c) as functions of the temperature. We calculated the dependence of ΘD(T) (as a function of temperature) plotting isobar curves as shown in Figure 7(a). The result in Figure 7(a) shows that the Debye temperature increased substantially with increasing pressure, from 494.33 K for P = 0 GPa to 692.87 K for P=35 GPa (at T=0 K), because as the pressure increases, all the velocities of the elastic waves also increase gradually and their increase directly affects the Debye temperature. Likewise, a mild nonlinear decreasing behavior of ΘD(T) can be observed with increasing temperature for applied pressures. This decrease in the Debye temperature as a function of temperature is characteristic of perovskite-type materials (Anderson, 1998). Interpreting the Debye temperature as the highest temperature that can be achieved as a result of a single normal vibration, it can be argued that the effect of pressure is the increase of cation-anion vibration frequencies while the effect of temperature is the expansive distortion of the structure, increasing the wavelength of the vibrations, decreasing the frequency and, consequently, the Debye temperature.

Figure 7(b) shows the dependence of the thermal expansion coefficient, α, with temperature and pressure. In the figure, α(T) decreases drastically with increasing pressure. On the other hand, it can be seen that at low temperatures (between 0 K and 300 K) α(T) grows rapidly with temperature, and above 300 K it grows smoothly, following an approximately linear behavior. In spite of the SHNO cubic structure, its perovskite-type character confers it an anisometric structure and, therefore, very great differences of the thermal expansion must be present along the different crystallographic directions. This behavior can be associated with the structural distortions to which the perovskitetype materials are sensitive because both temperature and pressure can give rise to inclinations and/or rotations of the HoO 6 and NbO 6 octahedrons, elongation of the structure in certain crystallographic directions, and eventual contractions in other directions.

Thus, the type of response in α(T) is associated with the type of distortion or transition occurring because of the application of temperature and pressure to the material. This circumstance could slightly divert any experimental results with respect to the theoretical result of Figure 7(b). On the other hand, the relatively low value of α(T) is characteristic of ceramic materials (Dedova, et al., 2016).

There was an aspect that we did not consider in the calculations: as the temperature increases, structural phase transitions can occur and the material could lose its cubic structure. Finally, the Grüneisen parameter, shown in Figure 7(c), shows a gradual decrease with the increase in applied pressure (for example, from 2.10 for P=0 GPa to 1.80 for P=35 GPa, at T=0 K). This behavior was observed for all the temperatures considered in the calculations. On the other hand, the Grüneisen coefficient exhibited a smooth and non-linear increase as temperature increased. These characteristics, observed in the Grüneisen parameter as a function of pressure and temperature, are related to the alteration in the vibration frequency of a crystalline lattice according to the discussion presented above in relation to the coefficient of thermal expansion (Kim, et al., 2017).

Conclusions

We obtained Sr 2 HoNbO 6 samples with a double cubic structure (Fmm space group) using the solid state reaction technique. Rietveld refinement of the XRD experimental data showed the occurrence of a superstructure with the presence of Ho3+ and Nb5+ cations interspersed along the crystallographic axes of the cell. The theoretical lattice parameters calculated from the minimization of energy with respect to the cell volume were in agreement with the experimental values determined by the Rietveld refinement of the XRD patterns. SEM images revealed the micrometric character of the granular surface with evidence of clusters constituted by nanometric crystallites. The measurements of magnetic susceptibility as a function of temperature indicated the paramagnetic feature of this complex perovskite material. The effective magnetic moment obtained from the Curie analysis of the paramagnetic behavior in curves of magnetic susceptibility as a function of temperature was 94,4% in agreement with the expected effective magnetic moment calculated from Hund rules. Calculations of the band and electronic structures suggest that this material behaves as an insulator for both spin-up and spin-down polarizations, with a gap energy of E=3,3 eV. On the other hand, an anomaly due to the 4f-Ho orbitals took place for the spindown configuration, and two intra-bands with high states/eV and a localized energy regime were observed in the middle of the gap. For this spin polarization, the band structure suggests that this material behaves as a semiconductor. The calculations of the thermodynamic properties from the state equation by means of the quasi-harmonic approximation of the Debye model showed that the effects of the interatomic vibrations gave rise to a specific heat that decreases with the applied pressure and evidenced a Dulong-Petit limit relatively high (240,63 J/mol.K) compared to other results reported for simple perovskite-type materials. Likewise, we conclude that CV≈CP, mainly due to the low value of the coefficient of thermal expansion. The Debye temperature showed an increasing dependence on the pressure and a decreasing one with the temperature. By contrast, both the coefficient of thermal expansion and the Grüneisen parameter decreased with the pressure and increased slightly as a function of temperature. Except for the high specificity of the Dulong-Petit heat value, the behavior observed in the thermodynamic properties was, in general, characteristic of the perovskite-like ceramic materials. We emphasize that the thermodynamic properties depend fundamentally on the value of the Grüneisen parameter, which is very sensitive both to the numerical errors in the derivation and to the smoothness of the energy equation. For this reason, it would be highly advisable to obtain experimental results of the thermodynamic properties in order to establish more objectively the microscopic mechanisms that originate them.

Acknowledgments

This work was partially financed by División de Investigación, sede Bototá, Universidad Nacional de Colombia, and Fonciencias, Universidad del Magdalena.

References

Anderson O.L. (1998). Thermoelastic properties of MgSiO3 perovskite using the Debye approach. Amer. Mineral. 83: 23-35. [ Links ]

Blaha P., Schwarz K., Madsen G.K.H., Kvasnicka D., Luitz J. (2001). WIEN2k, An Aug-mented Plane Wave + Local Orbitals Program for Calculating Crystal Properties. Karlheinz Schwarz, Techn. Universität Wien, Austria. ISBN 3-9501031-1-2. [ Links ]

Blanco, M.A., Francisco, E., Luaña, V. (2004). GIBBS: Isothermal- isobaric thermodynamics of solids from energy curves using a quasi-harmonic Debye model. Comp. Phys. Commun. 158: 57-72. [ Links ]

Bonilla M., Landínez-Téllez D.A., Arbey-Rodríguez J., Albino- Aguiar J., Roa-Rojas J. (2008). Study of half-metallic behavior in Sr2CoWO6 perovskite by ab initio DFT calculations. J. Magn. Magn. Mater. 320: e397-e399. [ Links ]

Corredor L.T., Landínez-Téllez D.A., Pimentel Jr. J., Pureur P., Roa-Rojas J. (2011). Magnetic, Structural and Morphological Characterization of Sr2GdRuO6 Double Perovskite. J. Mod. Phys. 2: 154-157. [ Links ]

Dedova E.S., Shadrin V.S., Petrushina M.Y., Kulkov S.N. (2016). The Study on Thermal Expansion of Ceramic Composites with Addition of ZrW2O8. Mater. Sci. Engin. 116: 012030(1)- 012030(5). [ Links ]

Deluque-Toro C.E., Arbey-Rodríguez M.J., Landínez-Téllez D.A., Moreno-Salazar N.O., Roa-Rojas J. (2014). First principles study of the structural and electronic properties of double perovskite Ba2YTaO6 in cubic and tetragonal phases. Phys. B. 455: 18-21. [ Links ]

Deluque-Toro C.E., Landínez-Téllez D.A., Roa-Rojas J. (2018). Ab-initio analysis of magnetic, structural, electronic and thermodynamic properties of the Ba2TiMnO6 manganite. DYNA. 85: 27-36. [ Links ]

Erten O., Meetei O.N., Mukherjee A., Randeria M., Trivedi N., Woodward P. (2011). Theory of Half-Metallic Ferrimagnetism in Double Perovskites. Phys. Rev. Lett. 107: 257201(1)-257201(4). [ Links ]

Hazen R.M. (1988). Perovskites. Sci. Am. 258: 74-81. [ Links ]

Hohenberg P., Khon W. (1964). Inhomogeneous electron gas. Phys. Rev. 136: B864-B871. [ Links ]

Howard C.H., Kennedy B.J., Woodward P.M. (2003). Ordered double perovskites - a group-theoretical analysis. Acta Cryst. B. 59: 463-471. [ Links ]

Khandy S.A., Gupta D.C. (2018). Electronic structure, magnetism and thermoelectric properties of double perovskite Sr2HoNbO6. J. Magn. Magn. Mater. 458: 176-182 [ Links ]

Khon W., Sham L.S. (1965). Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 140: A1133- A1138. [ Links ]

Kim J., Chen X., Shih P-C., Yang H. (2017). Porous Perovskite- Type Lanthanum Cobaltite as Electrocatalysts toward Oxygen Evolution Reaction. Sustainable Chem. Eng. 5: 10910-10917. [ Links ]

Kittel C. (1997). Introduction to Solid State Physics, 6th Ed. John Wiley & Sons, Inc., New York. p. 482 [ Links ]

Kruth A., West A.R. (2001). Electrical properties of the oxygendeficient perovskites, Ca2Mn2−xNbxOγ: 0 ≤ x ≤ 1.2, with Mn valence varying from +2.0 to +4.0. J. Mater. Chem. 11: 153-159. [ Links ]

Landínez-Téllez D.A., Llamosa D.P., Deluque-Toro C.E., Gil- Rebaza A.V., Roa-Rojas J. (2013). Structural, magnetic, multiferroic and electronic properties of Sr2ZrMnO6 double perovskite. J. Mol. Struct. 1034: 233-237. [ Links ]

Lee A.T., Marianetti C.A. (2018). Structural and metal-insulator transitions in rhenium-based double perovskites via orbital ordering. Phys. Rev. B. 97: 045102(1)-045102(5). [ Links ]

Llamosa D.P., Landínez-Téllez D.A., Roa-Rojas J. (2009). Magnetic and structural behavior of Sr2ZrMnO6 double perovskite. Phys B. 404: 2726-2729. [ Links ]

Lufaso M.W., Barnes P.W., Woodward P.M. (2006). Structure prediction of ordered and disordered multiple octahedral cation perovskites using SPuDS. Acta Cryst. B . 62: 397-410. [ Links ]

Madueño Q., Landínez-Téllez D.A., Roa-Rojas J. (2006). Production and characterization of Ba2NdSbO6 complex perovskite as a substrate for YBa2Cu3O7-d superconducting films. Mod. Phys. Lett. B. 20: 427-437. [ Links ]

Molegraaf H.J.A., Hoffman J., Vaz C.A.F., Gariglio S., van der Marel D., Ahn C.H., Triscone J-M. (2009). Magnetoelectric Effects in Complex Oxides with Competing Ground States. Adv. Mater. 21: 3470-3474. [ Links ]

Moreno L.C., Valencia J.S., Landínez-Téllez D.A., Arbey- Rodríguez J., Martínez M.L., Roa-Rojas J., Fajardo F. (2008). Preparation and structural study of LaMnO3 magnetic material. J. Magn. Magn. Mater. 320: e19-e21. [ Links ]

Murnaghan, F.D. (1944). The Compressibility of Media under Extreme Pressures. Proceedings of National Academy of Sciences, USA. 30: 244-247. [ Links ]

Perdew J.P., Burke K., Ernzerhof M. (1996). Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77: 3865-3868. [ Links ]

Pilo J., Miranda A., Trejo A., Carvajal E., Cruz-Irisson M. (2017). Bidimensional perovskite systems for spintronic applications. J. Mol. Model. 23: 322-325. [ Links ]

Qiang L., Duo-Hui H., Qi-Long C., Fan-Hou W. (2013). Phase transition and thermodynamic properties of BiFeO3 from first-principles calculations. Chin. Phys. B . 22: 037101(1)- 037101(4). [ Links ]

Triana C.A., Corredor L.T., Landínez-Téllez D.A., Roa-Rojas J. (2011). High temperature-induced phase transitions in Sr2Gd RuO6 complex perovskite. Mater. Res. Bull. 46: 2478-2483. [ Links ]

Triana C.A., Landínez-Téllez D.A., Arbey-Rodríguez J., Fajardo F., Roa-Rojas J. (2012). Electronic, crystal structure and morphological properties of the Sr2DyRuO6 double perovskite. Mater. Lett. 82: 116-119. [ Links ]

Triana C.A., Landínez-Téllez D.A., Roa-Rojas J. (2015). General study on the crystal, electronic and band structures, the morphological characterization, and the magnetic properties of the Sr2DyRuO6 complex perovskite. Mater. Charact. 99: 128. [ Links ]

Wondratschek, W. 2006. International Tables for Crystallography A, Chapter 8.3, p. 732-740. [ Links ]

Editor: Ángela Camacho Beltrán

Received: April 28, 2018; Accepted: June 25, 2018

*Corresponding autor: Jairo Roa Rojas, jroar@unal.edu.co

Authors’ contributions

. C.E. Deluque Toro, A.S. Mosquera Polo, and J. Roa- Rojas carried out the DFT and thermodynamic properties calculations. D.A. Landínez Téllez and J.I. Villa Hernández performed the synthesis of samples and the experimental measurements. J. Roa-Rojas analyzed the theoretical and experimental results, drafted, proofread, and reviewed the final document.

Conflict of interests

. The authors declare that there is no conflict of interest of any kind that affects the publication of the results of our research work.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License