PART IV

THE WURTZITE-ROCKSALT TRANSITION IN ZnO

By Katsuhiko Inoue

1975


I. ABSTRACT

Pressure induced phase transition from wurtzite to rocksalt structure of Zn0 has been investigated up to 85 kbar and 1200Ž. Phase identification at high pressure and temperature is based on the x-ray diffraction method using energy dispersive technique instead of the usual quenching method. The transition curve is approximately expressed by a linear relation, P(kbar)= -0.023 T(Ž) + 74. The lattice parameters of both phases at various temperatures along the phase boundary were also measured. The volume change and thus calculated transition entropy were found to be almost temperature independent. The compressibility measurements at room temperature indicated the reduction of the axial ratio c/a with increasing pressure for the wurtzite type phase. These distortions of the hexagonal lattice, which are likely to be related to the phase transition in ZnO, are discussed with special reference to the theory of Keffer and Portis in term of the partial polar binding in wurtzite type compounds.


II. INTRODUCTION
For the group II-VI compounds (e.g. BeO, ZnO, CdS), two structural types, the zincblende and wurtzite are well known with four-fold coordination of both cation and anion. The zincblende structure is face centered cubic and can be considered as composed of cubic close packed tetrahedra. On the other hand, the wurtzite structure belongs to the space group C6mc and consists of two hexagonal close packed sublattices which are displaced against each other along c axis by a distance of about 3c/8.
Many of these II-VI Compounds with zincblende and wurtzite structure have been studied under high pressure and found to undergo phase transitions to more dense structures generally to rocksalt or closely related structures (for a review see C.J.M.Rooymans, 1969) , which are common structure in transition metal mono-oxides and the coordination of both cation and anion is 6 in this structure. These phase transitions are characterized in general by a large drop in electrical resistance as well as the increase in coordination number (Samara and Drickamer, 1962).
The phase transitions of wurtzite or zincblende structure to rocksalt one are rapid and reversible, so that the quenching method often used for the determination of phase stability in oxides can't be used. P-T phase diagrams for some II-VI compounds have been established by means of DTA, piston-displacement technique and electrical resistance measurement.(Jayaraman et al., 1963 and Cline and Stephens, 1965). It is needless to say that in situ x-ray diffraction technique at high pressure and high temperature has intrinsic advantages for the establishment of P-T phase diagram. The use of this technique, however, has been limited to structure determination of high pressure polymorphs or compressibility measurement only near room temperature because of various experimental difficulties to achieve high temperatures. As reported in Part Ill, the application of SSD have much improved the high temperature capability in high pressure x-ray diffraction.
In this paper, phase stability between wurtzite and rocksalt phase in zincite (ZnO) up to 1200Ž and 85 kbar is determined by means of in situ x-ray diffraction method. The energy dispersive method was adopted for this purpose and pyrophyllite was used as pressure transmitting medium. The lattice parameters of both phases are also reported at various temperatures along the phase boundary. For zincite, Bate et al.(1962) have reported that the high pressure rocksalt can be retained to 1 atm after pressure and temperature treatment in contrast with the behavior of the majority of other II-VI compounds and they established P-T phase diagram of this compound by quenching method. Akimoto et al. (unpublished data) reported that high pressure phase of ZnO could not be quenched successfully to the ambient condition and pointed out the ambiguity in the P-T phase diagram for ZnO reported by Bates et al.
On the other hand, Soga and Anderson (1967) have been pointed out anomalous elastic behavior in polycrystalline ZnO i.e. shear wave velocity decreases with pressure. Similar phenomenon was also reported in ƒ¿-Quartz (Soga,1968). These results are significant to seismology and tectonophysics because a low value of dVs/dp allows a low velocity zone for shear waves to be compatible with a small thermal gradient. Anderson (1968) proposed a relationship of this anomaly with coordination number for pairs in the lattice and the ratio of the shear modulus to the bulk modulus by extending the Blackman's (1958) lattice dynamic theory. The negative pressure derivatives of shear modulus are also observed in other wurtzite type compounds such as CdS (Corll, 1967) and ZnS (Chang and Barsch, 1973).
Considering these anomalous elastic properties in wurtzite type Compounds, Corll (1967) proposed a mechanism for the pressure-induced structural transition in CdS based upon the possibility of microscopic instability developing under high applied stress. The proposed mechanism is a continual deformation from the wurtzite structure to the rocksalt one due to a shear of the (001) basal planes of wurtzite structure plus a linear compression of the hexagonal axis. Chang and Barsch (1973) also discussed a transformation mechanism for ZnS based on Born's (1939) mechanical stability criterion by using the second order elastic constants measured as a function of pressure and temperature.
The compression measurement by x-ray diffraction technique cannot give extensive information about shear instability directly as in the ultrasonic measurement. However, the experimental determination of the pressure effect on the lattice parameters of wurtzite type ZnO at room temperature is expected to give us some information about transformation mechanism from the pressure dependence of the axial ratio, c/a.


lIl. EXPERIMENTAL
A cubic anvil apparatus with energy dispersive x-ray diffraction system was used throughout this study. A detailed description of the apparatus and techniques are given in Part I and II. Pyrophyllite and boron nitride cells were used for the experiment for phase identification at high pressure and high temperature. The SSD was fixed at Bragg angles near 12‹. ZnO powder samples were intimately mixed with boron nitride powder to avoid the grain growth at high temperatures.
About thirty diffraction patterns were recorded at various P-T conditions extending to 85 kbar and 1200Ž. The pressure was not measured directly by NaCl internal standard, but was determined from the oil pressure of the 250 ton press. The relation between internal pressure and oil pressure was calibrated in advance for the same sample assembly by means of NaCl and Decker's eq. of state (1971). Temperature was measured by Chromel-Alumel thermocouple (0.2 mm in dia.) lead out through a pair of gaskets and the effect of pressure on its emf was not corrected.
In the compressibility measurement of wurtzite phase at room temperature, pressure medium made of amorphous boron powder filled plastics was used and the sample powder was intimately mixed with fine polyethylene powder. The reflections 101, 102, 110, 103 and 112 were used to calculate the lattice parameters by means of least squares method.


IV. RESULTS AND DISCUSSION
A typical energy dispersive powder diffraction pattern for rocksalt type ZnO just as brought to 85 kbar and 800Ž is shown in Fig, 1.

Fig. 1 Energy dispersive x-ray diffraction pattern at 85 kbar and 800Ž for ZnO. Indices are for high pressure rocksalt structure.

Reflections 200, 220, 311 and 222 of rocksalt structure are clearly resolved. Reflections of boron nitride sample holder, with which reflections from sample are often interfered in cases of low Z sample materials, are weak and negligible in this case because of the contrast of absorption in ZnO and boron nitride. Searching the wurtzite-rocksalt phase boundary, diffraction patterns were recorded in the interval of 100 deg. at various constant press loads. A series of patterns at 45 ton press load (corresponding to about 50 kbar) are shown in Fig. 2.

Fig. 2 Series of x-ray diffraction patterns to search the wurtzite-rocksalt phase boundary at 45 ton constant press load. patterns were recorded for 400 sec accumulating time sequentially from bottom to top. At 900Ž existence of both phases are clearly identified.

Exposure time for each pattern was only 400 sec, that seems enough to distinguish which phase is stable at the condition. Only Wurtzite phases are shown in the patterns from room temperature before heating up to 800Ž, however, at 900Ž reflection 220 of rocksalt type ZnO appears suddenly between the reflections 110 and 103 of wurtzite type phase. Further heating to 1000Ž causes the increase of intensities of peaks from the new phase relative to those from low temperature phase. After subsequent heating up to 1200Ž there appears to be no trace of wurtzite phase at last. After rapidly quenched from 1200Ž to room temperature, only metastable rocksalt phase can be seen and there seems to be no reversion to wurtzite phase. Finally after releasing the pressure, rocksalt phase disappears and only wurtzite phase with broad reflections is observed. The results of phase identification at various P-T conditions are summarized in Fig. 3.

 

Fig. 3 P-T phase diagram for ZnO. The broken line shows the results of Bates et a1. by quenching method. At P-T conditions the arrows indicate, the lattice parameters of both phases, therefore volume changes were measured.

The results of Bates et al. (1962) are also plotted together in the figure. Since the old pressure scale was adopted in their original report, the correction was given based on the new fixed-point pressure scale calibrated by NaCl internal standard (Jeffery et al., 1966). A remarkable discrepancy between the Bates' results and the present ones is shown not only in the position in P-T space but in the sign of the slope of the phase boundary. In the quenching experiment the rocksalt phase could be retained only with ammonium chloride as a catalyst and furthermore, new phase has not been prepared free from contamination by Zincite even with long runs. The disagreement with our in situ phase identification is presumably coming from the failure of quenching in pure ZnO. The phase boundary determined in the present study can be expressed;
P(kbar) = -0.023 T(Ž) + 74. (1)
The slope of the boundary is relatively small but definitely negative in comparison with the positive value of + 0.0425 kbar/deg in the quenching experiment. The lattice parameters of both phases along the equilibrium line are listed in Table 1.

The transition enthalpies and entropies are calculated by using Clausius-Clapeyron eq. to be almost temperature independent as shown in the table. Figure 4 shows a typical pattern used for the compression measurement of wurtzite type phase at room temperature.

Fig.4. Energy dispersive x-ray diffraction pattern of wurtzite type ZnO at 75 kbar and room temperature.

The effect of pressure on the volume, linear compressions and axial ratio, c/a for wurtzite phase of ZnO is shown in Fig.5 and 6,respectively.

Fig. 5 Pressure dependence of volume compression for wurzite type ZnO at room temperature. The solid line shows the ultrasonic equation of state.

Fig. 6 The effect of pressure on the lattice parameters and the axial ratio of wurtzite type ZnO at room temperature.

The solid line shows the ultrasonic eq. of state based on 2 parameter Birch eq. using the isothermal bulk modulus and its pressure derivative measured by Soga and Anderson (1967). The present results gives appreciably smaller bulk modulus than that determined in the ultrasonic experiment. Similar discrepancies are frequently encountered in many other high pressure x-ray diffraction experiments and have been partly accounted in terms of the stress inhomonogeneity in the aggregate of fine grains having different elastic properties (Sato, 1973). Accordingly it can be said that the present results on volume and linear compressions do not represent the intrinsic elastic properties of ZnO and probably there appears no correlation with the anomaly in shear property in polycrystalline ZnO. 0n the other hand, it should be noted that the pressure dependence of axial ratio, however, still remains valid in a qualitative sense.
Even at 1 atm, the axial ratio for ZnO is 1.602, far from the ideal value (1.633) for the hexagonal close packing. Such a small value in c/a in ZnO among the other wurtzite type compounds is considered to be coming from the large difference in ionic electronegativity with increasing pressure, c/a decreases monotonically to about 1.591 at 80 kbar in a further leaving sense from the ideal value. In the following discussion this anomalous variation in c/a will be related only qualitatively to the phase transition to rocksalt structure. Keffer and Portis(1957) have accounted the slight departure from ideal c/a ratio and from ideal sublattice displacement, u in wurtzite type compounds by the forces arising from partial polar binding and calculated these long range forces in terms of a postulated charge }fe on ion and certain lattice sums in order to balance against the elastic, piezoelectric and dielectric contributions caused by the distortion.
For the departure of c/a from the ideal value, they derived the following equation;

where N is the number of ions of one sign per unit volume, fe, the assumed effective charge on each ion, sij and dij the elastic compliance and piezoelectric strain constants, respectively. In the case of ZnO, we evaluated f to be a reasonable value, about 0.78 at 1 atm by using the x-ray and single-crystal elastic data (Bateman, 1968). The accompanying piezoelectric effect have been assumed to play only a minor role to the lattice contraction and the effect of charge transfer have been also neglected in the calculation. It can be concluded, however, that ZnO is considerably ionic and the anisotropy in the Coulomb interaction is mainly responsible to the distortion from the ideal hcp lattice at 1 atm in wurtzite type ZnO. There is no direct measurement on the pressure derivatives of the second order compliance constants for ZnO, but it is probable to assume that the quantity in the square parenthesis in eq. (2) is roughly proportional to s11 because s12 and s13 are small as compared to s11 and s33 and also the order of magnitude of Ýs11 /Ýp is about -1 x 10-24 cm4/dyne2.
Only less than 10% increase in the right side of eq. (2) was calculated at 50 kbar, for instance. The observed increase of about 25% at this pressure in the left side of eq. (2) i.e. the departure of c/a from the ideal value, cannot be explained without taking account of the change in ionic contribution with increasing pressure. Consequently the effective charge on each ion is expected to become larger at the higher pressures. In the ZnX (X=0, S, Se, Te) series, the large drops in the room temperature resistivity ZnS, ZnSe and ZnTe were found near 240-245, 165, 140-145 kbar, respectively (Samara et a1,1962) and confirmed to be the first order phase transitions from wurtzite structure to rocksalt one by means of high pressure x-ray diffraction technique (Smith and Martin, 1965). The sequence of the transition pressures in these three compounds is consistent with that of the transitions from diamond structure to ƒÀSn one in the group IV elements, where the element with lower atomic number undergoes phase transition at lower pressure. Relatively low transition pressure in ZnO observed in this study is violating the tendency mentioned above.
Jayaraman et al. (1963a,b) have pointed out in the generalization of the phase transitions in group IV elements, Ill-V and II-IV compounds that the wurtzite type compound with larger electronegativity difference tend to transform to rocksalt structure at lower pressure and to become metallic at higher pressure. It may be concluded that the relatively ionic nature of binding even at 1 atm among other covalent ZnX (X=S, Se, Te) and presumably the tendency of becoming more ionic at higher pressures are responsible to the transition to 6-coordinated rocksalt structure at such a low pressure.


REFERENCES
Anderson,0.L. (1968) J. Geophys. Res., 73, 7707.
Anderson,0.L. and R.C.Liebermann (1970) Phys. Earth Planet. Interiors, 3, 61.
Bates,C.H., W.B.White and A.Roy (1962) Science, 131, 993.
Bateman,T.B. (1962) J. Appl. Phys., 33, 3309.
Blackman,M. (1957) Proc. Phys. Soc. London Ser. B, 70, 827.
Born,M. (19391 J. Chem. Phys., 7, 591.
Chang,E. and G.R.Barsch (1973) J. Phys. Chem. Solids, 34, 1543.
Cline,C.F. and D.A.Stephens L1965) J. Appl. Phys., 36, 2869.
Corll, J..A. (1967) Phys. Rev., 157, 623. .
Decker,D.L. (1971) J. Appl. Phys., 42, 3239.
Jayaraman,A., W,.Element,Jr. and G.C.Kennedy (1963) Phys. Rev., 130, 540.
Jayaraman,A., W.Klement,Jr. and G.C.Kennedy (1963b) Phys. Rev., 130, 2277. .
Jeffery,R.N., I.D.Barnett, H.A.Vanfleet and H.T.Hall (1966) J. Appl. Phys., 37, 3172.
Keffer,F. and Portis,A.M. (1957) J. Chem. Phys., 27, 675.
Rooymans,C.J.M. (1969) Advances in high pressure research, vol.2 ed. R.S.Bradley, 1-100, Academic Press, London, New York.
Samara,G.A. and H.G.Drickamer (1962) J. Phys. Chem. Solids, 23, 457. .
Sato,Y. (1974) read at IV AIRAPT Internat. Conf., Kyoto.
Smith,P.L. and J.E,Martin (1965) Phys. Letters, 19, 541.
Soga,N. and O.L.Anderson (1967) J. Appl. Phys., 38, 2985.
Soga,N. (1968) J. Geophys. Res., 73, 827.