S. Rezaa,
M. Maazabc and
M. S. Islam
*a
aDepartment of Physics, University of Rajshahi, Rajshahi 6205, Bangladesh. E-mail: sislamru@gmail.com
bUNESCO-UNISA Africa Chair in Nanosciences-Nanotechnology, College of Graduate Studies, University of South Africa, Muckleneuk Ridge, PO Box 392, Pretoria, South Africa
cNanosciences African Network (NANOAFNET), Materials Research Dept., iThemba LABS-National Research Foundation of South Africa, 1 Old Faure Road, Somerset West, PO Box 722, Western Cape 7129, South Africa
First published on 8th June 2023
We performed a pressure-driven study of zinc pyrovanadate, Zn2V2O7, using the first-principles approach under the framework of density functional theory (DFT). Zn2V2O7 crystalizes in a monoclinic (α-phase) structure with the space group C2/c at ambient pressure. In comparison with the ambient phase, there are four different high-pressure phases, namely β, γ, κ and δ, found at 0.7, 3.8, 4.8 and 5.3 GPa, respectively. The detailed crystallographic analysis as well as their structures is consistent with the theory and experiment reported in the literature. All phases including the ambient phase are mechanically stable, elastically anisotropic and malleable. The compressibility of the studied pyrovanadate is higher than that of the other meta- and pyrovanadates. The energy dispersion of these studied phases reveals that they are indirect band gap semiconductors with wide band gap energies. The band gap energies follow a reduced trend with pressure except the κ-phase. The effective masses for all of these studied phases were computed from their corresponding band structures. The values of energy gaps obtained from the band structures are almost similar to the optical band gap obtained from the optical absorption spectra, as estimated by the Wood–Tauc theory.
Tremendous attention has recently been paid to study Zn2V2O7, isostructural with thortveitite mineral Sc2Si2O7,12 due to its above-mentioned fascinating behavior. First, Makarov et al. experimentally crystallized zinc pyrovanadate and observed two polymorphic structures within the transition temperature range of 608–620 °C.13 The low-temperature phase of α-Zn2V2O7 immediately crystallizes in a monoclinic system with the space group (SG) of C2/c.14 The other polymorph of Zn2V2O7, high-temperature β-phase, forms a thortveitite structure in the same crystal with the SG of C2/m.15 Besides structural transitions upon temperature stimulation, the behavior of Zn2V2O7 is interesting under high-pressure conditions. A much recent report on Zn2V2O7 based on XRD patterns at a pressure up to 12 GPa has surprisingly indicated at least three structural phase transitions but no pressure-driven transitions have been shown in metavanadate ZnV2O6 and orthovanadate Zn3V2O8.16 Interestingly, the studied pyrovanadate under ambient conditions is highly compressible in comparison with the meta- and ortho-vanadates. The first structural transition was found at 0.7 GPa, where the ambient phase of α-polymorph was transformed into β-polymorph although the structure was monoclinic.14,16 The second high-pressure phase transition of the γ-phase appeared at 3.8 GPa, and the proposed structure was triclinic with SG P. The third phase transition, namely, the post γ-phase appeared at 10.8 GPa, but the structure was unidentified.16
Very recently, first-principles calculations based on density functional theory (DFT) revealed several new additional phases of κ-Zn2V2O7 isomorphic to Ni2V2O7 (SG P21/c),17 δ-Zn2V2O7 structure corresponding to Hg2V2O7 (SG Pnma),18 ε-Zn2V2O7 similar to Pb2V2O7 (SG P21/c),19 and ω-Zn2V2O7 corresponding to Sr2V2O7 (SG P),20 in association with α-, β- and γ-phases on the basis of crystal-chemistry arguments.21,22 In this work, they predicted that the post γ-phase could be the κ-phase, although the structures are largely different. This comparison to define the high-pressure phase is still a subject of debate. The present work used first-principles DFT-encoded CASTEP to determine the post γ-phase in connection with the above-mentioned phases except ε- and ω-phases. In association with energy dispersion and optical absorption, the purpose is to show mechanical stability under compression by studying the elastic behavior of all these phases. In this work, the behavior of this pyrovanadate and their phase-related discussion will be carried out under pressure.
Phases | Crystal system | Space group | Z | Lattice parameters | Bond length (Å) | Ref. | ||||
---|---|---|---|---|---|---|---|---|---|---|
Axial length (Å) | Axial angle (°) | |||||||||
a | b | c | Zn–O | V–O | ||||||
α | Monoclinic | C2/c | 4 | 7.455 | 8.380 | 10.254 | α = γ = 90 | 2.157 | 1.789 | This work |
7.429 | 8.340 | 10.098 (ref. 14) | β = 111.069 | |||||||
β | C2/m | 2 | 6.975 | 8.503 | 5.046 | α = γ = 90 | 2.137 | 1.866 | ||
6.932 | 8.440 | 5.033 (ref. 15) | β = 108.846 | |||||||
γ | Triclinic | P![]() |
2 | 13.636 | 5.444 | 5.108 | α = 77.916 | 2.059 | 1.906 | |
13.621 | 5.235 | 4.923 (ref. 16) | β = 107.948 | |||||||
γ = 130.960 | ||||||||||
κ | Monoclinic | P21/c | 4 | 6.651 | 8.324 | 9.450 | α = γ = 90 | 2.113 | 1.722 | |
6.615 | 8.394 | 9.492 (ref. 21) | β = 100.133 | |||||||
δ | Orthorhombic | Pnma | 4 | 6.840 | 3.538 | 19.530 | α = β = γ = 90 | 2.323 | 1.926 | |
6.865 | 3.578 | 19.601 (ref. 21) |
As illustrated in Fig. 1a, the cations (Zn) in the ambient phase occupy sites with the nearest neighbors of five oxygen (O) atoms as compared with the cations in the thortveitite-like structure coordinated to be six-fold.14,16 The average distances of both of these polyhedral made of cations and anions are about 2 Å. However, the vanadium (V) atoms are tetrahedrally coordinated. Each pair of tetrahedra [VO4], slightly distorted, is connected by a common O atom located at a common corner to form the [V2O7]4− pyrovanadate anions present in the thortveitite-like minerals.30–32 The two cations per formula unit in this structure donate four electrons to the [V2O7] group, which is converted into [V2O7]4− anions based on the concept given by Zintl-Klemm.30 A similar structural pattern is seen in the first transition phase of Zn2V2O7 except the coordination of cations and anions, as shown in Fig. 1b. The Zn atoms in the β-phase form a polyhedra, ZnO6, but the vanadium is coordinated with five O atoms. In the case of γ-phase, the two cations, each being five-coordinated, are connected by two common oxygen atoms, while a single oxygen atom is common between the anionic VO5 and the two cationic polyhedra, as depicted in Fig. 1c. The κ-phase in Fig. 1d is structurally consistent with the thortveitite-like structure, where Zn and V atoms are bonded with six and four O atoms, respectively. Here, one oxygen atom connects two tetrahedra of V and a polyhedron of Zn. Higher coordination of cations (octahedra) and anions (hexahedra) is formed in the δ-phase, as shown in Fig. 1e. The three hexahedra in this structure are linked by common oxygen atoms.
In order to study the mechanical stability of a material, the mechanical behavior of a crystal lattice can be expressed by its matrix of second-order elastic/stiffness constants:
![]() | (1) |
→ | C11 | C12 | C13 | C15 | C22 | C23 | C25 | C33 | C35 | C44 | C46 | C55 | C66 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Phases | |||||||||||||
α | 135 | 61 | 61 | 5 | 145 | 64 | 7 | 120 | 11 | 44 | 8 | 50 | 32 |
β | 356 | 41 | 91 | 23 | 184 | 105 | −10 | 186 | 9 | 30 | −29 | 38 | −15 |
κ | 169 | 72 | 72 | 17 | 221 | 100 | 1 | 156 | 8 | 56 | 3 | 40 | 55 |
δ | 181 | 109 | 81 | — | 275 | 90 | — | 201 | — | 34 | — | 53 | 67 |
γ | 162 | 87 | 79 | 7 | 161 | 91 | 5 | 161 | 4 | 41 | 2 | 42 | 37 |
→ | C14 | C16 | C24 | C26 | C34 | C36 | C45 | C56 | |||||
γ | −0.5 | −12 | −3 | −6 | −12 | −10 | −6 | 4 |
Depending on the stiffness constants, we determine the mechanical stability of these studied phases of Zn2V2O7 using the well-known Born–Huang criteria.34 All these phases are mechanically stable except β, which has negative values for C66. However, the majority of the criteria based on stiffness values in the β-phase are fulfilled. It is noteworthy to mention that the stability of γ is unaffected by modest negative values of some elastic constants. This would rather suggest that the crystal in these phases may have minor internal stresses.35
As a polycrystalline mechanical property, we also estimated the modulus of elasticity, e.g. the bulk modulus BH = (BV + BR)/2, where BV is Voigt's bulk modulus and BR is Reuss's bulk modulus; and the shear modulus GH = (GV + GR)/2, where GV is Voigt's shear modulus and GR is Reuss's shear modulus using the average approximations given by Hill.36 Young's modulus (E) and Poisson's ratio (ν) were also computed using the relationships: E = 9BHGH/(3BH + GH) and ν = (3BH − E)/(6BH), respectively.37 The calculated values of BH, GH, E and ν, were compared with the previously reported DFT results, which are listed in Table 3. It is well known that the value of bulk modulus describes the compressibility of a material. In this context, β, γ, κ and δ phases are the least compressible in comparison with the α phase.21 Moreover, the compressibility can be enhanced by increasing the unit cell volume because of an inverse relationship between the bulk modulus and the unit cell volume under ambient conditions.21,38 The calculated values of bulk modulus were higher than the result published in the literature.21,38 Again the value of either bulk modulus or shear modulus can indirectly measure the hardness of a material. Using Hill's approximation, the shear modulus is much smaller than the bulk modulus, reflecting that the shear deformation is easier as shown in Table 3.36 Generally, the stiffness of an elastic material can be identified by the value of Young's modulus, so that the β phase is, in this case, highly stiffer than the others. The value of Pugh's ratio (B/G) distinguishes the ductile (>1.74) and brittle (<1.74) nature of materials.39 As listed in Table 3, the most malleable behavior is expected in these studied phases of Zn2V2O7. The malleability can also be correlated with the value of Poisson's ratio. A material will be ductile if ν > 0.26, otherwise it will be brittle.39 The calculated values of Poisson's ratio, presented in Table 3, confirm the malleable character.
The mechanical performance of a material can also be signified by the number of useful indicators tabulated in Table 4. The machinability index μM (=BH/C44) is such an indicator that can be useful for describing the plasticity and lubricating behavior of a material.40–42 It is observed that a lower value of C44 gives better dry lubricity. However, a high value of μM exhibits excellent lubricating properties, high plastic strain value, lower friction value and feed forces. All phases in Zn2V2O7 suggest a high degree of machinability. Within the studied phases, β and δ are more machinable which is consistent with the high value of BH. The Kleinman parameter ζ (= (C11 + 8C12)/(7C11 + 2C12)) typically has a value between 0 and 1. According to the theory given by Kleinman,43 it is reasonable to consider that the bond bending contribution instead of bond stretching is dominant for the mechanical strength in almost all phases of Zn2V2O7 except the β phase due to its high value of C11. Indeed, this parameter reaches the maximum value when C11 equals C12. A perfect isotropic crystal has a zero value of the universal anisotropy factor AU.44 The degree of anisotropy is defined from the deviation of AU = 0. In this case, all phases of Zn2V2O7 possess anisotropic nature. Understanding hardness is crucial for comprehending the elastic and flexible capabilities of a material. It is instructive to notice that the β-phase is reasonably harder than the other phases of Zn2V2O7, which is compatible with the calculated value of bulk modulus in Table 3.45
Phases | μM | ζ | AU | HV | ρ | νa | ΘD | Ref. |
---|---|---|---|---|---|---|---|---|
α | 1.98 | 0.58 | 0.40 | 3.13 | 3.83 | 3.42 | 431 | This work |
β | 4.33 | 0.27 | −2.64 | 8.14 | 4.04 | 4.51 | 573 | |
γ | 2.65 | 0.66 | 0.27 | 1.96 | 4.20 | 3.37 | 434 | |
κ | 1.99 | 0.56 | 0.37 | 4.36 | 4.25 | 3.67 | 482 | |
δ | 3.92 | 0.71 | 0.46 | 4.18 | 4.84 | 3.75 | 503 |
Using the computed polycrystalline elastic moduli, the average elastic wave velocity νa is defined as follows:46
![]() | (2) |
![]() | (3) |
In addition to the computed values regarding the elastic moduli, e.g., bulk modulus/linear compressibility, Young's modulus, shear modulus and Poisson's ratio, it is possible to visualize the nature of anisotropy in three dimensions using ELATE, which is an open-source Python module for the analysis of elastic tensors.27 The first two moduli can be represented by a single unit vector and characterized by two angles θ(0 ≤ θ ≤ π) and φ(0 ≤ φ ≤ 2π) in spherical coordinates but the latter two are described by two orthogonal unit vectors, which can be parameterized by three angles θ, φ and ξ(0 ≤ ξ ≤ 2π). In the case of linear compressibility, most of the materials show positive compressibility indicating compression axially upon isostatic pressure. There may be some materials whose linear compressibility can be both positive and negative.27,48,49 For a negative linear compressibility, the structure expands in a particular direction upon pressure but the volume as a whole can be reduced. Although the moduli as a function of two unit vectors are difficult to represent in 3D, we determine Poisson's ratio (shear modulus not shown), following the work by Gaillac et al.,27 in the spherical coordinates (θ, φ) with the requirement of minimal and maximal values over all possible values of an extra dimension ξ. The 3-dimensional sketch of the largest and smallest possible values of the linear compressibility and Poisson's ratio of anisotropic surfaces for all of these studied phases except β is shown in Fig. 2. All the moduli expressed by either one unit vector or two orthogonal unit vectors as depicted in 3D diagrams are shown in Fig. 2, and it is observed that these 3D shapes strongly deviate from the spherical one, indicating a strong anisotropic nature compatible with the calculation of the universal anisotropy factor. The minimum and maximum values of different elastic moduli estimated by ELATE are tabulated in Table 5.
Phases | Young's modulus (GPa) | Linear compressibility (T Pa−1) | Shear modulus (GPa) | Poisson's ratio | ||||
---|---|---|---|---|---|---|---|---|
Min | Max | Min | Max | Min | Max | Min | Max | |
α | 71 | 128 | 2.9 | 5.1 | 27 | 50 | 0.11 | 0.51 |
γ | 76 | 135 | 1.4 | 4.4 | 30 | 49 | 0.15 | 0.52 |
κ | 86 | 158 | 1.6 | 5.1 | 36 | 62 | 0.18 | 0.57 |
δ | 104 | 199 | 1.3 | 3.4 | 34 | 72 | 0.13 | 0.55 |
To analyze the electronic behavior of zinc pyrovanadate, we first computed the energy dispersion and the corresponding density of states (DOS) of the studied phases under high-pressure situations using first-principles DFT calculations with the PBE functional. Fig. 3a shows the electronic band structure and the DOSs of an ambient α-phase of Zn2V2O7. It is shown that an ambient phase of studied pyrovanadate is an indirect band gap semiconductor, indicated by the solid black line arrows, with a wide gap of 2.45 eV listed in Table 6. The recent DFT calculation21-encoded CRYSTAL14 program using hybrid functionals, e.g., B3LYP, HSE06, found the high value of energy band gap in comparison with the experimental value (∼3.5 eV) determined from the study of photoluminescence.50 The diffuse-reflectance spectra measurements51,52 reported a gap value 2.86 eV, which is almost consistent with the present calculations. Earlier estimation regarding the hetero-structured photocatalyst made by TiO2/vanadates presented a similar value of band gap energy of an ambient phase.53 Very recent computations21 using hybrid functionals have overestimated the band gap of vanadates inconsistent with the experimental value, but our calculations using the PBE functional are comparatively compatible with the experimental result given by diffuse-reflectance measurements.51,54
Phases | Functionals | Gap | Energy gap (eV) | Energy gap (eV) (W–T theory) | DOS (states per eV) | Effective mass (m0) | |
---|---|---|---|---|---|---|---|
Electron | Hole | ||||||
α | PBE | Indirect | 2.45 (this work) | 2.00 | 9.75 | 0.027 | 0.016 |
PBE | 2.29 (ref. 21) | ||||||
HSE06 | 3.98 (ref. 21) | ||||||
B3LYP | 4.20 (ref. 21) | ||||||
Expt. | — | 2.5, 2.86 (ref. 51–53) | |||||
3.50 (ref. 50) | |||||||
β | PBE | Indirect | 2.24 (this work) | 1.75 | 4.33 | 0.008 | 0.035 |
PBE | 2.00 (ref. 21) | ||||||
HSE06 | 3.63 (ref. 21) | ||||||
γ | PBE | Indirect | 2.14 (this work) | 1.65 | 3.72 | 0.015 | 0.022 |
PBE | Direct | 2.13 (ref. 21) | |||||
HSE06 | 3.78 (ref. 21) | ||||||
κ | PBE | Indirect | 2.77 (this work) | 2.10 | 14.10 | 0.527 | 0.014 |
PBE | 2.80 (ref. 21) | ||||||
HSE06 | 4.56 (ref. 21) | ||||||
δ | PBE | Indirect | 1.54 (this work) | 1.49 | 3.19 | 0.032 | 0.052 |
PBE | 1.46 (ref. 21) | ||||||
HSE06 | 2.91 (ref. 21) |
As can be presented in Fig. 3, the other phases of this pyrovanadate in our computations show the similar nature of semiconductivity, where the position of the top of the valence band (VB) and the bottom of the conduction band (CB) is different. According to the work given by Diaz-Anichtchenko et al., the γ-phase exhibits a direct band gap semiconductor.21 The estimated values of the energy gap of β, γ, κ and δ-phases are 2.24, 2.14, 2.77 and 1.54 eV, respectively, as presented in Table 6. Moreover, the energy band gap of the post-γ phase (not shown) is nearly the same as the gap value of the γ-phase, so the prior prediction regarding the post-γ phase21 is not consistent yet. The right panel of Fig. 3 illustrates the calculated density of states (DOSs) for these bulk phases of Zn2V2O7 with atom projected contributions (partial DOS). It is seen that the top of the VB is mainly dominated by the O 2p orbitals, while the states near the bottom of the CB are dominated by V 3d orbitals with partial hybridization of O 2p orbitals. However, the contribution from Zn orbital electrons has no significant influence near the Fermi level so the energy band gap of these studied phases of pyrovanadate is almost comparable with the other metavanadates and orthovanadates.21,55,56
Out of these band gap energies for the phases of Zn2V2O7, the highest energy gap is seen in the κ-phase, which was a similar finding of the Díaz-Anichtchenko group under the DFT study.21 In their calculations, they have shown the variation of the energy gap under pressure. All phases of the pyrovanadate show a reduced trend as a function of pressure, resulting in an increase in the hybridization between V 3d and O 2p electrons similar to the metavanadate like ZnV2O6, except the κ-phase. However, the opposite trend of the band gap energy with pressure indicates the lack of attraction between bonding and antibonding states.28 Moreover, the pressure coefficient at ambient pressure in the κ-phase is a completely opposite sign and a much smaller value in comparison with the other phases. Such repulsion as well as the value of the pressure coefficient is responsible for showing the highest band gap energy in the κ-phase.
The effective masses for the electron or hole
at the band extremes can also be calculated from the energy dispersion using the following equation:
![]() | (4) |
The absorption coefficient is an important characteristic to understand the electronic nature of a material. The kind of electronic transition (either indirect or direct) in the energy bands can also be determined by studying the optical absorption spectra as a function of energy studied for the unpolarized light. The absorption spectra until 40 eV for the studied phases are depicted in Fig. 4a. It is seen that, the absorption for all these phases starts at some energy, which confirms the semiconducting behaviour already mentioned in the energy dispersion results. The absorption spectra imply that the light photon can efficiently be absorbed within the range of ∼7–25 eV, equivalent to the ultraviolet region; while the energy exceeds 25 eV, the absorption falls sharply.
![]() | ||
Fig. 4 Energy dependence of (a) absorption spectra along with the Wood–Tauc plots (inset) and (b) conductivity for the studied phases of pyrovanadate. |
The optical band gap energy (Eg) can also be estimated by the Wood–Tauc (W–T) theory, as plotted in the insets of Fig. 4a, given by the relation αℏν ∝ (ℏν − Eg)k, where α is the absorbance, ℏ is the Planck constant, ν is the frequency and k is a constant associated with the different optical transitions.55 For an allowed indirect transition, the value of k is 2. By fitting with the best linear relation for a value of 2, the values of the band gap energy for different phases were also calculated, listed in Table 6. The spectra of photoconductivity, shown in Fig. 4b, start with some photon energy, which indicates that all the studied phases of Zn2V2O7 have a band gap as evident from the energy dispersion.
This journal is © The Royal Society of Chemistry 2023 |