Abstract

We have studied the structural and thermophysical properties of cadmium oxide (CdO), using the Three-Body Potential (TBP) model. Phase transition pressures are associated with a sudden collapse in volume. The phase transition pressures and related volume collapses obtained from this model show a generally good agreement with available experimental others data. The thermophysical properties like molecular force constant, Debye temperature, and so forth, of CdO are also reported.

1. Introduction

The group of IIB-VIA oxides have presented a great deal of interest because of their applications in various technologies [1]. The semiconducting compounds of this group crystallize mostly in the zincblende (B3), wurtzite (B4), or both structures. Cadmium oxide (CdO) is one of the binary oxides having important electronic, structural, and optical properties. Cadmium oxide occurs naturally as the rare mineral monteponite. CdO is a semiconductor with a band gap of 2.16 eV at room temperature. It normally crystallizes in a cubic sodium chloride (NaCl) rock-salt structure, with octahedral cation and anion centers. However under pressure it shows a first-order structural phase transition from NaCl (B1) to CsCl (B2) structure [2]. First-principles calculations of the crystal structures, and phase transition, and elastic properties of cadmium oxide (CdO) have been carried out with the plane-wave pseudopotential density functional theory method by Peng et al. [3]. Liu et al. studied the B1 to B2 phase transition pressure at about 90.6 GPa for CdO [4]. Guerrero-Moreno et al. [5] observed the ground-state properties of CdO with B1 to B2 structure, using the first 2 principles calculations.

We have applied the Three-Body Potential (TBP) model to the present compound to study the high pressure phase transition and other properties. The need of inclusion of three-body interaction forces was emphasized by many workers for the betterment of results [68]. Earlier calculations for B1-B2 transitions were based on two-body potential mainly. They concluded that possible reasons for disagreements include the failure of the two-body potential model. Since these studies were based on two-body potentials and could not explain Cauchy violations (C12C44). They remarked that results could be improved by including the effect of nonrigidity of ions in the model. This Three-Body Potential (TBP) model consists of long-range Coulomb energy, three body interactions corresponding to the nearest neighbour separation, vdW (van der Waal) interaction, and energy due to the overlap repulsion represented by Hafemiester and Flygare (HF) [9] type potential and extended up to the second-neighbour ions. The purpose of this work is to investigate the structural and thermophysical properties of CdO.

2. Potential Model and Method

Application of pressure directly results in compression leading to the increased charge transfer (or three body interaction effect [10]) due to the deformation of the overlapping electron shell of the adjacent ions (or nonrigidity of ions) in solids.

These effects have been incorporated in the Gibbs free energy (𝐺=𝑈+𝑃𝑉𝑇𝑆) as a function of pressure and three-body interactions (TBI) [10], which are the most dominant among the many body interactions. Here, 𝑈 is the internal energy of the system equivalent to the lattice energy at temperature near zero and 𝑆 is the entropy. At temperature 𝑇=0 K and pressure (𝑃) the Gibbs free energies for rock salt (B1, real) and CsCl (B2, hypothetical) structures are given by 𝐺B1(𝑟)=𝑈B1(𝑟)+𝑃𝑉B1(𝑟),(1)𝐺B2𝑟=𝑈B2𝑟+𝑃𝑉B2𝑟.(2) With 𝑉𝐵1(=2.00𝑟3) and 𝑉𝐵2(=1.54𝑟3) as unit cell volumes for 𝐵1 and 𝐵2 phases, respectively. The first terms in (1) and (2) are lattice energies for 𝐵1 and 𝐵2 structures and they are expressed as 𝑈𝐵1(𝑟)=𝛼𝑚𝑧2𝑒2𝑟12𝛼𝑚𝑧𝑒2𝑓(𝑟)𝑟𝐶𝑟6+𝐷𝑟8+6𝑏𝛽𝑖𝑗𝑟exp𝑖+𝑟𝑗𝑟𝜌+6𝑏𝛽𝑖𝑖exp2𝑟𝑖1.414𝑟𝜌+6𝑏𝛽𝑗𝑗exp2𝑟𝑗1.414𝑟𝜌,𝑈(3)𝐵2𝑟=𝛼𝑚𝑧2𝑒2𝑟16𝛼𝑚𝑧𝑒2𝑓𝑟𝑟𝐶𝑟6+𝐷𝑟8+8𝑏𝛽𝑖𝑗𝑟exp𝑖+𝑟𝑗𝑟𝜌+3𝑏𝛽𝑖𝑖exp2𝑟𝑖1.154𝑟𝜌+3𝑏𝛽𝑗𝑗exp2𝑟𝑗1.154𝑟𝜌(4) with 𝛼𝑚 and 𝛼𝑚 as the Madelung constants for NaCl and CsCl structure, respectively. C(C′) and D(D′) are the overall vander der Waal coefficients of B1 (B2) phases, 𝛽𝑖𝑗(𝑖,𝑗=1,2) are the Pauling coefficients. Ze is the ionic charge and 𝑏(𝜌) are the hardness (range) parameters, 𝑟(𝑟) are the nearest neighbour separations for NaCl (CsCl) structure and 𝑓(𝑟) is the three body force parameter.

These lattice energies consists of long-range Coulomb energy (first term), three-body interactions corresponding to the nearest neighbour separation 𝑟(𝑟) (second term), vdW (van der Waal) interaction (third term), and energy due to the overlap repulsion represented by Hafemeister and Flygare (HF) type potential and extended up to the second neighbour ions (remaining terms).

3. Results and Discussion

The Gibb’s free energies contain three model parameters [𝑏,𝜌,𝑓(𝑟)]. The values of these parameters have been evaluated using the first- and second-order space derivatives of the cohesive energy (𝑈) expressed as 𝑑𝑈𝑑𝑟𝑟=𝑟0𝑑=0,2𝑈𝑑𝑟2𝑟=𝑟0=9𝑘𝑟0𝐵,(5) and following method adopted earlier [11]. Using these model parameters and the minimization technique, phase transition pressures of CdO have been computed. The input data of the crystal and calculated model parameters are listed in Table 1. We have followed the technique of minimization of Gibbs free energies of real and hypothetical phases. We have minimized 𝐺𝐵1(𝑟) and 𝐺𝐵2(𝑟) given by (3) and (4) at different pressures in order to obtain the interionic separations 𝑟 and 𝑟 corresponding to 𝐵1 and 𝐵2 phases associated with minimum energies. The factor Δ𝐺[𝐺𝐵1(𝑟)𝐺𝐵2(𝑟)] plays an important role in stability of structures. The phase transition occurs when Δ𝐺 approaches zero (Δ𝐺0). The phase transition pressure (𝑃𝑡) is the pressure at which Δ𝐺 approaches zero. At 𝑃𝑡 these compounds undergo a (𝐵1-𝐵2) transition associated with a sudden collapse in volume showing a first-order phase transition. Figure 1 shows our present computed phase transition pressure for NaCl-type (𝐵1) to CsCl-type (𝐵2) structures in CdO at 90 GPa. The present phase transition pressure is illustrated by arrow in Figure 1. The calculated values of phase transition pressure have been listed in Table 2 and compared with their experimental and other theoretical results. It is interesting to note from Table 2 and Figure 1 that the phase transition pressures (𝑃𝑡), obtained from our model, are in general in closer agreement with experimental data [4] and match equally well with other theoretical results [5].

The compression curves are plotted in Figure 2. The values of the volume collapses (Δ𝑉(𝑝)/𝑉(0)) are depicted in Table 2. The experimental and theoretical values of volume collapses are not available for the present compounds. It is clear that during the phase transition from NaCl to CsCl, the volume discontinuity in pressure volume phase diagram identifies the occurrence of first-order phase transition and the same trend as the other theoretical approach. In Figure 2 the pressure versus volume graph has been plotted.

In additional to know the behaviour of the interionic distance with pressure for the present oxide, we present the variation of nearest-neighbor (nn) and next-nearest neighbor (nnn) distances for both the 𝐵1 and 𝐵2 phases with pressure in Figure 3. The interionic distances of the present oxide decrease on increasing the pressure. The open circles represent the nearest-neighbor (nn) and solid circles represent next nearest neighbor (nnn) distance in Figure 3 for CdO.

To further increase the applicability of our model, we have calculated the molecular force constant (𝑓), infrared absorption frequency (𝜐0), Debye temperature (𝜃𝐷), and Grunneisen parameter (𝛾) which are directly derived from the cohesive energy, 𝑈(𝑟).

The compressibility is well known to be given by 𝛽=3𝐾𝑟0𝑓(6) in terms of molecular force constants 1𝑓=3𝑈𝑆𝑅𝑘𝑘2(𝑟)+𝑟𝑈𝑆𝑅𝑘𝑘(𝑟)𝑟=𝑟0.(7) With 𝑈𝑆𝑅𝑘𝑘(𝑟) as the short range nearest neighbour (𝑘𝑘) part of 𝑈(𝑟) given by the last three terms in (3) and (4). This force constant 𝑓 leads to the infrared absorption frequency with the knowledge of the reduced mass (𝜇) of the crystals. The thermal expansion coefficient (𝛼𝑣) can be calculated with the knowledge of specific heat (𝐶𝑣). The expressions have been given in our earlier paper [12].

We have calculated the thermophysical properties of CdO. The thermophysical properties provide us the interesting information about the substance. The Debye characteristic temperature 𝜃𝐷 reflects its structure stability, the strength of bonds between its separate elements, structure defects availability, and its density. The calculated thermophysical properties have been listed in Table 3. Due to the lack of experimental and theoretical data, we could not compare them. We have compared the value of Debye temperature with theoretical results provided by Peng et al. [3]. Our result shows the same trend as reported by others. As to the best of our knowledge the value of the thermal properties for present compounds have not yet been measured or calculated, hence our results can serve as a prediction for future investigations.

In view of the overall attainments, it may be concluded that there is generally a good agreement of Three-Body Potential (TBP) model with the available experimental and theoretical values. Finally, it may be concluded that the present model has successfully predicted the compression curves and phase diagrams giving the phase transition pressures, associated volume collapses, and elastic properties correctly for cadmium oxide.