The present work addresses the atomic structure and electronic properties of a recently discovered mineral qatranaite (CaZn2(OH)6·2H2O). The present study was performed theoretically by means of density functional theory- (DFT-) based calculations within the frame of local density approximation (LDA) and general gradient approximation (GGA). To determine the energy band gap width, we carried out the ultraviolet-visible spectroscopy (UV-Vis) measurements. The structure relaxation performed with use of LDA and GGA provides results matching the experimentally determined crystal parameters. Interestingly, in contrast to existing interpretation of experimental data, our DFT calculations revealed energy gap of direct characteristics. Accordingly, our UV-Vis experiments yield the band gap width of 3.9 eV.
Infrastruktura PL-GridNarodowe Centrum Nauki2016/21/B/ST8/027371. Introduction
Qatranaite named after the village Al Qatrana (situated 15 km southeast from the Amman Aquaba Asia-Jordan Desert Highway) is considered as a natural analogue of synthetic calcium hexahydroxodizincate dehydrate [1, 2]. The mineral, discovered in an altered pyrometamorphic spurrite rocks in the Daba-Siwaqa region of the Hatrurim complex located in Jordan [3–6], was for the first time described in 2016 by Stasiak et al. as CaZn2(OH)6·OH2O [7]. Qatranaite crystallized in the temperature range 150 to 200°C, similarly as Se-bearing thaumasite or afwillite, in the form of flatten (010) crystals, growing up to 0.3 mm [7, 8].
Qatranaite unit cell consists of octahedrally coordinated Ca2+ ions linked to four hydroxy groups and two water molecules, whereas Zn2+ ions are tetrahedrally coordinated with four OH− groups and form pyroxene-like chains [Zn2(OH)6]2−. The mineral may display charge deficiency in case all oxygen atoms participate in the formation of OH groups, which certainly affects electronic properties of this interesting material [8].
The potential application of calcium hexahydroxodizincate dehydrate-like systems sparked hope for improved charge/discharge reversibility in Zn battery electrodes. Furthermore, they are expected to restrict Zn-electrode shape change and inhibit dendrite formation [1, 2, 9–11]. Additionally, electrochemical properties of CaZn2(OH)6·2H2O make it suitable for negative electrodes in rechargeable zinc-air batteries [12, 13].
The energy gap of calcium hexahydroxodizincate dehydrate has been a subject of experimental investigations [2] from which the direct energy gap width of 3.1 eV has been concluded. Our UV-Vis measurements of qatranaite revealed that the spectrum can be equally well interpreted, assuming the presence of direct or indirect energy gap. In order to disclose the nature of the energy gap in the mineral, we carried out the quantum, DFT-based, electronic structure investigations by means of the pseudopotential method employed in the Quantum Espresso package [14]. To verify reliability of the calculations, we have also performed the structural optimization and compared the results with experimental data. An ab initio prediction for the energy band gap was compared with experimental results provided by UV-Vis measurements.
2. Computational and Experimental Details
An ab initio calculations were carried out using the Quantum Espresso code [14] by applying the scalar relativistic ultrasoft pseudopotentials [15] for Ca(3s23p64s2), H(1s2), and O(2s22p6) elements within both LDA and GGA calculations. For Zn, we used the same type pseudopotentials with (3d9.74s24p0.3) and (3d104s2) valence state configuration for LDA and GGA, respectively. Calculations were performed for two approximations of exchange correlation (XC) energy functional. The LDA XC functional in the form parameterized by Perdew and Zunger [16] and the GGA XC functional proposed by Perdew et al. [17] were employed. The plane wave kinetic energy and charge density cutoff of 46 and 368 Ry, respectively, were applied within LDA, while in the case of GGA, the corresponding cutoff values equaled 53 and 368 Ry. The structure-optimization, following the Broyden–Fletcher–Goldfarb-Shanno (BFGS) scheme, was carried out within LDA and GGA approaches. The usage of the 12 × 12 × 12 Monkhorst–Pack mesh [18] in reciprocal space enabled us to achieve satisfactory convergence in the structural optimization. The total energy of the qatranaite unit cell was evaluated with an accuracy of 5 meV, while forces were calculated with uncertainty less than 0.08 and 5.21 meV/Å for LDA and GGA, respectively. To determine precisely the details of electronic density of states, the additional calculations using the 20 × 20 × 20 Monkhorst–Pack mesh were performed.
To study the energy gap of qatranaite, the ultraviolet-visible spectroscopy (UV-Vis) was performed with the use of microspectrophotometer (CRAIC Technologies) equipped with standard halogen lamp and Zeiss 15x objective. The sample was prepared in the form of petrographic polished plate with the qatranaite crystals (diameter of 0.3 mm) embedded in the mineral matrix. The experiments were performed at room temperature and ambient pressure. To estimate the band gap width, we use the formula proposed by Wood and Tauc [19]:(1)hν⋅α∼hν−Egn,where α is the absorbance, h stands for the Planck constant, ν defines the photon’s frequency, Eg denotes the optical band gap energy, and n is a constant related to different kinds of electronic transitions. The n parameter equals 0.5, 2, 1.5, and 3 for direct, indirect, allowed, and forbidden transitions, respectively.
3. Results and Discussion3.1. Structural Properties
According to the experimental results provided by Stasiak et al. [7], the CaZn2(OH)6·2H2O mineral crystallizes in the P21/c-type structure with the following lattice parameters: a = 6.3889(8) Å, b = 10.969(1) Å, c = 5.7588(8) Å, and β = 101.95(1)°. The coordination of each Zn atom to four oxygen atoms and Ca atom to six oxygen atoms prompted the authors to imagine the qatranaite unit cell as a particular set of tetrahedral ZnO4 and octahedral CaO6 clusters (Figure 1). Peculiarity of the qatranaite unit cell is related to the presence of hydroxyl groups [OH] attached to hydrogen bridges to H2O molecules [20], which ensures the electric charge balance.
Schematic of the crystalline structure of qatranaite showing the tetrahedrally coordinated zinc (ZnO4; yellow) and octahedrally coordinated calcium (CaO6; blue) clusters. The red and green spheres denote oxygen and hydrogen atoms, respectively.
Tables 1 and 2 summarize the results of ab initio structural optimization of qatranaite unit cell carried out with the use of two approximations for XC energy functional. The calculated results were compared with available experimental data. It is worth emphasizing that the output of DFT-LDA calculations is in well agreement with the previously reported diffraction data [7, 8] (Table 1). The computed (LDA) and measured a, b, c lattice parameters and β angle differ less than 0.2, 2.3, 3.5, and 2.1%, respectively. GGA approach give slightly worse results (corresponding differences are 4.1, 2.9, 5.6, and 0.6%). It can be concluded that LDA XC functional better describes the nature of bonding in qatranaite.
Lattice parameters of qatranaite estimated with the LDA and GGA and compared with the measured values.
Parameters
LDA
GGA
Exp. [7]
a (Å)
6.388
6.652
6.389
b (Å)
10.720
11.286
10.969
c (Å)
5.960
6.081
5.759
β (°)
104.05
102.52
101.95
Experimental and calculated coordinates of Wyckoff positions in qatranaite single crystal.
Atom
Exp. [8]
DFT-LDA
DFT-GGA
X
Y
Z
x
Y
z
x
y
z
Zn1-4e
0.5332
0.6650
0.6639
0.5300
0.6670
0.6385
0.5313
0.6669
0.6536
O1-4e
0.4687
0.6814
0.3122
0.4451
0.6933
0.3030
0.4540
0.6891
0.3108
H1-4e
0.5227
0.61710
0.2675
—
—
—
—
—
—
O2-4e
0.3380
0.5427
0.7455
0.3557
0.5309
0.7117
0.3440
0.5363
0.7084
H2-4e
0.3226
0.5543
0.8576
—
—
—
—
—
—
O3-4e
0.8349
0.6236
0.7394
0.8377
0.6388
0.6993
0.8311
0.6308
0.7010
H3-4e
0.8739
0.6108
0.85960
—
—
—
—
—
—
O4-4e
0.0373
0.6603
0.2367
0.0206
0.6422
0.2193
0.0262
0.6486
0.2274
H4-4e
0.1712
0.6704
0.2490
—
—
—
—
—
—
H5-4e
−0.01980
0.7304
0.2400
—
—
—
—
—
—
Ca1-2a
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
For hydrogen atoms, “—” stands for experimental values of coordinates.
In Table 2, we compare experimental and optimized coordinates of Wyckoff positions in the structure of CaZn2(OH)6·2H2O mineral. To avoid the issue of high mobility of hydrogen atoms [21], their experimentally obtained coordinates were frozen during structural optimization. We are satisfied with the overall agreement between calculated and measured coordinates of Wyckoff positions [8]. The noticeable difference in Wyckoff position can be observed for oxygen atoms, which may be attributed to their relative high mobility [22].
Based on structural data collected in Table 2, we also calculated an interatomic distances in tetrahedral ZnO4 and octahedral CaO6 clusters. The obtained results are presented in Table 3. The comparison of DFT structural results with experimental data confirms sufficient reliability of our calculations.
Interatomic distances in ZnO4 and CaO6 clusters.
Interatomic distances (Å)
Zn1
Ca1
O1
O1
O2
O3
O2
O2
O3
O3
O4
O4
LDA
1.96
1.94
1.95
1.93
2.34
2.34
2.29
2.30
2.30
2.29
GGA
2.05
2.00
2.00
1.99
2.40
2.40
2.39
2.35
2.35
2.39
Exp. [7]
1.91
1.99
1.95
1.94
2.37
2.37
2.37
2.34
2.34
2.34
The atomic symbols (e.g., O1 and O2) denote atoms occupying different Wyckoff positions.
3.2. Electronic Properties
Physical properties of a novel material are essential for its intended applications in electronics and photonics. Consequently, the energy band structure and the details of qatranaite’s partial density of states (PDOS) are of common interest. Indeed, our DFT-determined band diagram (Figure 2) indicates the insulating property of qatranaite with a direct band gap, at the Γ-point of reciprocal k-lattice, of 3.2 eV (LDA) and 3.3 eV (GGA).
The electronic band structure of crystalline qatranaite determined within LDA (green lines) and GGA (red lines). The zero of the energy scale is shifted to the top of the valence band.
Closer inspection of Figure 2 reveals that the selection of LDA and GGA of exchange-correlation energy does not significantly affect the shape of energy bands. Furthermore, the partial DOS analysis (Figure 3) reveals the electronic states near the top of the valence-band filled by O-p and Zn-d electrons with slight contribution of Ca electrons.
Partial density of states determined for qatranaite. The zero of the energy scale is shifted to the top of the valence band.
3.3. UV-Vis Absorption Spectroscopy
The results obtained from our DFT calculations demonstrated that qatranaite is the direct band gap material, and consequently, the n value of 0.5 have to be used in equation (1). The energy band gap, Eg = 3.9 eV, was estimated by fitting equation (1) to the experimental points at the linear portion of the curve (Figure 4). The latter value exceeds the one deduced from DFT calculations (Eg = 3.2 eV). This confirms the common perception that DFT calculations frequently lead to underestimation of the band gap [23]. Therefore, it may be concluded that the true energy gap in qatranaite should significantly exceed the calculated DFT value.
UV-Vis spectrum of qatranaite (black dots). Linear fit, according to equation (1), shows the energy gap width Eg.
Finally, the calculations and experimental results obtained for natural qatranaite were compared with previously published data for synthetic counterpart [2]. Thus, Xavier et al. measured the energy gap Eg of 3.1 eV, the value significantly lower than the one obtained by our UV-Vis absorption spectroscopy measurements. Such a high discrepancy is due to disputable assumptions imposed by Xavier et al. [2], who fitted their absorption results in such a way that stipulates n = 2 (equation (1)). The differences between experimentally derived energy gap for qatranaite (3.9 eV) and synthetic CaZn2(OH)6·2H2O (3.1 eV) [2] may stem from chemical and structure imperfections present in mineral qatranaite.
4. Conclusions
In sum, we have demonstrated the DFT calculations within LDA and GGA schemes accurately predict the atomic structure, which encourage the further study of qatranaite electronic properties. The calculated lattice parameters at the ground state conditions are in agreement with the experimental results. It is contend, therefore, that LDA is more appropriate for estimation of lattice parameters of qatranaite. The LDA (GGA) predicted direct band gap energy at Γ point of the first Brillouin zone provided the lower value 3.2 eV (3.3 eV) compared to our experimentally obtained 3.9 eV data. Moreover, it is worth noting that our theoretical and experimental results reveal the presence of direct gap in qatranaite mineral where available experimental data for its synthetic analogue pointed on indirect energy gap. We believe that the output of our calculations contains information being of importance to the modern electronic and photonic industries.
Data Availability
The raw data from all implemented techniques (UV‐VIS, DFT) used to support the findings of this study are available from the corresponding author upon request.
Disclosure
The authors would like to inform that a small part of their simulation results were previously published in a poster form.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was partially supported by PL-Grid Infrastructure. DCH acknowledges the support from the National Science Centre of Poland (Grant no. 2016/21/B/ST8/02737).
YangC.-C.ChienW.-C.ChenP.-W.WuC.-Y.Synthesis and characterization of nano-sized calcium zincate powder and its application to Ni-Zn batteries2009391394410.1007/s10800-008-9637-92-s2.0-58149515938XavierC. S.SczancoskiJ. C.CavalcanteL. S.A new processing method of CaZn2(OH)6·2H2O powders: photoluminescence and growth mechanism200911122173217910.1016/j.solidstatesciences.2009.09.0022-s2.0-70450240962BentorY.Israel19603Paris, FranceCentre National de la Recherche Scientifique80GrossS.The mineralogy of the Hatairim formation, Israel197770180SokolE.NovikovI.ZateevaS.VapnikY.ShagamR.KozmenkoO.Combustion metamorphism in the Nabi Musa dome: new implications for a mud volcanic origin of the Mottled Zone, Dead Sea area201022441443810.1111/j.1365-2117.2010.00462.x2-s2.0-77955133377VapnikY.SharyginV.SokolE.ShagamR.Paralavas in a combustion metamorphic complex: Hatrurim Basin, Israel200718121StasiakM.GaluskinE. V.KuszJ.IMA commission on new minerals nomenclature and classification (CNMNC), Qatranaite, IMA 2016-02420168032915922VapnikY.GaluskinE. V.GaluskinaI. O.“QatranaiteCaZn2(OH)6·2H2O—a new mineral from altered pyrometamorphic rocks of the Hatrurim Complex, Daba-Siwaqa, Jordan, in review,” European Journal of Mineralogy, 2019StahlR.JacobsH.AnorgZ.Zur kristallstruktur von CaZn2(OH)6·2H2O199762381287128910.1002/zaac.19976230819ZhuX. M.YangH. X.AiX. P.YuJ. X.CaoY. L.Structural and electrochemical characterization of mechanochemically synthesized calcium zincate as rechargeable anodic materials200333760761210.1023/a:10249992071782-s2.0-0042733320YuJ.YangH.AiX.ZhuX.A study of calcium zincate as negative electrode materials for secondary batteries20011031939710.1016/s0378-7753(01)00833-32-s2.0-0035977370MinY.-J.OhS.-J.KimM.-S.ChoiJ.-H.EomS.Calcium zincate as an efficient reversible negative electrode material for rechargeable zinc–air batteries20182417MainarA. R.IruinE.ColmenaresL. C.An overview of progress in electrolytes for secondary zinc-air batteries and other storage systems based on zinc201815130432810.1016/j.est.2017.12.0042-s2.0-85040337783GiannozziP.BaroniS.BoniniN.QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials2009213939550210.1088/0953-8984/21/39/3955022-s2.0-70349568754http://www.quantum-espresso.org/pseudopotentials/pslibrary/PerdewJ. P.ZungerA.Self-interaction correction to density-functional approximations for many-electron systems198123105048507910.1103/physrevb.23.50482-s2.0-26144450583PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple199677183865386810.1103/physrevlett.77.38652-s2.0-4243943295MonkhorstH. J.PackJ. D.Special points for Brillouin-zone integrations197613125188519210.1103/physrevb.13.51882-s2.0-1842816907WoodD. L.TaucJ.Weak absorption tails in amorphous semiconductors1972583144315110.1103/physrevb.5.31442-s2.0-0001134590LinT.-C.MollahM. Y. A.VempatiR. K.CockeD. L.Synthesis and characterization of calcium hydroxyzincate using X-ray diffraction, FT-IR spectroscopy, and scanning force microscopy19957101974197810.1021/cm00058a0312-s2.0-33751154999GasparD.PereiraL.GehrkeK.GallerB.FortunatoE.MartinsR.High mobility hydrogenated zinc oxide thin films201716325526210.1016/j.solmat.2017.01.0302-s2.0-85011320899BoccardM.RodkeyN.HolmanZ. C.High-mobility hydrogenated indium oxide without introducing water during sputtering20169229730310.1016/j.egypro.2016.07.0832-s2.0-85014429441PerdewJ. P.ParrR. G.LevyM.BalduzJ. L.Density-functional theory for fractional particle number: derivative discontinuities of the energy198249231691169410.1103/physrevlett.49.16912-s2.0-0001132752