# High Pressure Study of Structural and Electronic Properties of PbSe

Authors: P. Bhambhani; K. Kabra; B. K. Sharma; G. Sharma
Journal: Journal of Solid State Physics
 (2014)

---

## Abstract

High pressure structural phase transition and electronic properties have been investigated using the linear combination of atomic orbitals (LCAO) method with two exchange-correlation approximations, the generalized gradient approximation (GGA) and local density approximation (LDA). The present study shows phase transitions from B1 to B27 and B27 to B2 at 6.24 GPa and 16.39 GPa, respectively. Lattice constant, bulk modulus, and energy gap of pressure-induced PbSe are found to be in good agreement with previous theoretical and experimental results. Variation of electronic band structure with pressure shows direct band gap alongL point of the Brillouin zone.


---

## Body

## 1. Introduction

Systematic studies of IV–VI semiconductor compounds on theoretical as well as experimental background have given overwhelming properties which not only are important for the scientific inquest but have wide scale technological applications [1–9]. These compounds possess intriguing array of narrow energy gap, high carrier mobilities, high dielectric constant, positive temperature coefficient of the band gap and negative pressure dependence of the band gap, and high anisotropy, which are the basic requirements for the application in the branch of the thermoelectronics, optoelectronics, and spintronics [10–15].Such fascinating properties of the IV–VI semiconductor compounds become a subject of extensive research on the lead chalcogenides PbX (X = S, Se, Te). On account of narrow forbidden gap of PbSe ~0.26 eV (300 K), anomalous order of band gap, and low thermal conductivity, it appeared to be a potential candidate for the technological applications in the field of thermoelectric devices [10, 16–18]. Further, low transition pressure and narrow energy gap of the intermediate phase of this compound make it a suitable candidate for the fabrication of other electronic devices [18, 19]. The anomalous properties of PbSe and IV–VI semiconductor compounds as compared to III–V and II–VI semiconductors are due to interaction of Pbs-electron valence bands with the Sep-electrons valence band [12].Excellent properties of this material were aggravated to investigate pressure-induced structural and electronic properties. In the last few decades, a number of investigations on its pressure-induced phase transition and structural properties have been reported; see, for example, [9–12, 18–23]. At ambient conditions, PbSe crystallizes in sixfold-coordinated rock-salt type B1 structure, with space group Fm3
-m. At pressure between 3 and 5 GPa, it transforms from B1 phase to an orthorhombic intermediate phase [11]. The intermediate phase may be GeS (B16, Pbnm), FeB (B27, Pnma), and CrB/TlI (B33, Cmcm) [19, 21, 24, 25]. Further, with increasing pressure between 15 and 16 GPa, structural phase transition takes place from the orthorhombic phase to eightfold coordinated CsCl B2-type phase (space group Pm3
-m) [8–11, 19, 21]. The phase transition from low to high pressure, that is, from intermediate to B2-type phase, is reported as semiconducting to metallic transition [11, 25–27]. Chattopadhyay et al. [8, 9] from X-ray diffraction studies using synchrotron radiation reported TlI, B33-type structure for the intermediate phase of PbSe, while Ovsyannikov et al. reported it to be of GeS-type [22]. This intermediate phase has also been experimentally proposed to be orthorhombic Pnma structure [28]. The theoretical investigations on the intermediate phase of PbSe have less attention in the literature. Thus, the present paper deals with the pressure-induced structural phase transition and electronic properties of PbSe at different pressures using LCAO method.


## 2. Computational Details

Quantum mechanical periodic LCAO calculations for the ground state total crystal energy are performed using the computer code CRYSTAL06 [29]. There are two fundamental schemes to compute the electronic structure and related properties of solids: one is the Hartee-Fock (HF) approximation and the other is the DFT [30, 31]. In HF, the exchange potential is included exactly and additional terms describe the correlation effects approximately but, in the DFT, both exchange and correlation effects are included. To gain the quantitative validity of results, we have used both the GGA and LDA approximations. The Kohn-Sham Hamiltonian for GGA is constructed by considering the exchange and correlation scheme of Becke [32] and Perdew-Burke-Ernzerhof (PBE) [33], respectively, while LDA approximation is based on Becke [32] and von Barth-Hedin (VBH) [34] schemes, respectively. The basis sets for the S were taken from [35] and due to nonavailability of all electron basis sets of Pb, only its valence part was considered in the present computations [36]. The self-consistent calculations are performed using 29, 35, and 343k-points in the irreducible Brillouin zone (BZ) for B1, B2, and B27 phases, respectively, with sufficient tolerances. The level of numerical approximation in evaluating the Coulomb and exchange series appearing in the SCF equations for periodic systems is controlled by five tolerances [29]. In the present calculations, tolerances of the order of 10−6, 10−6, 10−6, 10−6, and 10−10 were considered. To achieve self-consistency, 70% mixing of successive cycles is considered and self-consistency is achieved within 10, 12, and 18 cycles for B1, B2, and B27 phases, respectively.The cohesive energy [37] from the Compton profile data can be derived as follows:
(1)
E
coh
=
∫
0
∞
p
z
2
J
s
p
z
-
J
f
a
p
z
d
p
z
,
where p
z is component of electron momentum along the z-axis and J
s and J
f
a represent Compton profiles of solid and free atoms, respectively. The solid Compton profile (J
s) values are calculated from DFT-LCAO method and the free-atom Compton profile (J
f
a) values are taken from Biggs et al. [38–40].


## 3. Results and Discussion

### 3.1. Phase Transition and Structural Parameters

The structural parameters were theoretically evaluated first by calculating the total lattice energies over a set of different volumes around the equilibrium values and then were fitted to the third order Birch-Murnaghan equation of state (EOS) using the following relations:(2)
E
V
=
E
0
+
9
V
0
B
0
16
V
0
V
2
/
3
-
1
3
B
0
′
+
V
0
V
2
/
3
-
1
2
6
-
4
V
0
V
2
/
3
,
where E
0 is minimum energy, V
0 is corresponding volume, and B
0 is the bulk modulus at zero pressure and B
0
′
=
d
B
0
/
d
P. The third order Birch-Murnaghan isothermal equation of state [41, 42] relating pressure and volume is given by
(3)
P
V
=
3
B
0
2
V
0
V
7
/
3
-
V
0
V
5
/
3
×
1
+
3
4
B
0
′
-
4
V
0
V
2
/
3
-
1
.
Calculations were performed in both GGA and LDA frameworks and the respective parameters, that is, lattice constant (a), bulk modulus (B
0), and pressure derivative of bulk modulus (B
0
′) for B1, B27, and B2 phases of PbSe, are given in Table 1. As per the device application, high bulk modulus and its pressure derivative are essential parameters for photovoltaic applications in optoelectronic devices [43, 44]. The bulk modulus also defines resistance to change in volume against any mechanical deformation [44]. Therefore, we computed the values of bulk modulus and its derivative using Birch-Murnaghan EOS and summarized them in Table 1. The value of bulk modulus obtained from Birch-Murnaghan EOS agrees well with the value reported by Rached et al. [45]. This variation of B
0 is in accordance with the theory of Cohen [46] who obtained an empirical expression for the bulk modulus based on the nearest-neighbor distance. As seen from Table 1, the pressure derivative of bulk modulus (B
0
′) is almost the same for all the three phases with average value of 4.0. On comparing the structural parameters obtained from two different schemes, that is, GGA and LDA, we observed that lattice parameter obtained from LDA approach is larger than GGA while the bulk modulus obtained from GGA is higher as compared to the one obtained from LDA.Table 1
Calculated equilibrium lattice parameter (a), bulk modulus (B
0), pressure derivatives of bulk modulus (B
0
′), and cohesive energy (E
coh) for PbSe.
Parameters
Present calculations
Experimental results
Other calculations
GGA
LDA
NaCl-type, B1 phase
a
0 (Å)
6.08
6.17
6.117 [29], 6.133 [9], 6.124 [56], 6.117 [57]
6.211 [10], 6.098 [23], 6.226 [12], 6.214 [12]
B
0
(GPa)
66.96
64.72
28.8 [8], 54.1 [56], 54.1 [58]
48.27 [10], 60.8 [23], 44.5 [12], 49.5 [12], 66.0 [42]
B
0
′
4
3.7
4.1 [8]
4.35 [10], 4.56 [23], 3.64 [12], 3.8 [12]
E
coh (eV)
6.95
—
—
7.39 [44]
Orthorhombic Pnma phase (B27, FeB type)
a
0 (Å)
a = 11.38
a = 11.48,
a = 11.61,
a = 8.07,
b = 4.374
b = 4.474,
b = 4.00,
b = 4.49,
c = 4.259
c = 4.359
c = 4.39 [29]
c = 6.29 [10]
B
0 (GPa)
68.94
65.5
—
—
B
0
′
4.9
4.9
—
—
CsCl-type, B2 phase
a
0 (Å)
3.709
3.74
3.379 [9]
—
B
0 (GPa)
70
64.72
—
62.0 [12]
B
0
′
3.9
3.7
—
2.40 [12]Figures1(a) and 1(b) show the variation of the total energy versus volume computed for PbSe in B1, B27, and B2 phases obtained from both GGA and LDA approximations together with EOS data. As evident from Figures 1(a) and 1(b), the energy state of B1 phase is lower than that of the other two phases. This phase could be regarded as a stable phase of PbSe system [26].The calculated energy versus volume for B1, B27, and B2 phases of PbSe: using the PBE (a) and VBH (b) correlation functional. The scattered points show calculated energies and the solid lines show the fittedE
(
V
) curves according to Birch-Murnaghan equation of state (EOS).
(a)
(b)In order to investigate the phase transition pressure (P
t) for intermediate and CsCl-type phase, we computed variation in enthalpy with pressure using GGA and LDA approximations considering certain constraint of temperature for all the three phases of PbSe compound. Since a stable phase of any crystal is coupled with its minimum free energy, we have constrained the temperature to zero in order to obtain a thermodynamically stable phase which has lowest enthalpy,H =E + PV, under given pressure [26, 47]. At P
t enthalpy of the two consecutive phases is equal [16, 48]. Difference in enthalpy of high pressure phases (B27 and B2) and that in stable phase (B1) for both GGA and LDA approximations are plotted in Figures 2(a) and 2(b), respectively. From GGA approach, the transition from the stable B1 phase to Pnma intermediate phase is predicted to be 6.24 GPa and the transition from intermediate to CsCl-type phase is 16.39 GPa, while the P
t evaluated from LDA approach is 7.97 GPa for transition from NaCl phase to an orthorhombic Pnma (B27, FeB type) phase and 18.25 GPa for transition from Pnma (B27) phase to the CsCl-type (B2) phase.Enthalpy versus pressure calculated from first-principles total energy for the B1, B27, and B2 structures of PbSe using the PBE (a) and VBH (b) correlation functional.
(a)
(b)These investigated transition pressures are reported in Table2. Since the transition pressures derived from GGA approximation are closer to the experimental results as compared with the transition pressures derived from LDA approximation, for this reason further investigation is done with GGA approximation. The recent investigations of Li et al. [49] on the pressure-induced structure transitions of PbSe using density functional theory are in agreement with the present data. The calculation of cohesive energy from the experimental Compton profile is a difficult task due to weighting of p
2 in the high momentum region [50, 51]. Therefore, we have used the DFT-PBE profile for the calculation of cohesive energy of PbSe. The calculated cohesive energy along with available data are given in Table 1. From this table, it is evident that our calculated value of E
coh shows fairly good agreement with earlier data [44].Table 2
Transition pressure (in GPa) for PbSe.
Transition
Present
Experimental
Other calculations
GGA
LDA
B1 → B27
6.24
7.97
4.5 [9], 4.5 [8]
7.45 [10], 10.0 [11]
B27 → B2
16.39
18.25
16.0 [9], 16.0 [8]
18.76 [10], 20.0 [11]The pressure versus volume data for B1, B27, and B2 phases were computed, using the calculated values ofB
0 and B
0
′, using (3), and plotted in Figure 3. The curve shows that at P
t there is sudden collapse in volume.Figure 3
Variation of volume with increasing pressure for PbSe using GGA.


### 3.2. Pressure Dependent Electronic Properties

Knowledge of energy band gaps at different pressure is important to provide the information about their potential efficiency in fabricating various optoelectronic and thermoelectronic devices. In order to compute the energy band gap at different pressures, the equilibrium lattice constants are used and different symmetry points likeW (0.5, 0.25, 0.75), L (0.5, 0.5, 0.5), Γ (0.0, 0.0, 0.0), X (0.5, 0.0, 0.0), W (0.5, 0.25, 0.75), and K (0.375, 0.375, 0.75) are considered in the BZ. The energy band gaps of PbSe in B1 phase at different hydrostatics pressures are reported in Table 3.Table 3
Band gap of PbSe in B1 phase at different pressures.
Pressure (GPa)
Present calculation
Other calculations
−4
0.338
—
−2
0.252
—
0
0.185
0.21 [13], 0.15 [12], 0.16 [52]
2
0.118
—
4
0.05
—As seen from Figure4, the valence band maximum and the conduction band minimum lie at L point of the BZ. The occupied 6s-band of Pb atom lies close to the valence band of Se atom which is of p-type. Since the symmetry of Pb 6s band is the same as that of the highest valence band at the L point, there exists strong level repulsion which determines the feature of band structure and gives direct band gap along L point of the BZ at all considered pressures [12].Figure 4
Energy band structure of PbSe in B1 phase.Since the band structures at all considered pressures are similar with only a small change around the Fermi level, it is desirable to present band structure of B1 phase only at zero pressure (Figure4). The variation of band gap with pressure for B1 phase is plotted in Figure 5. It shows that with increasing pressure minima of conduction band come closer to the Fermi energy level and thus reduce the band gap. A similar feature of energy gap with pressure is reported by Nabi et al. [52]. A number of groups have observed the zero pressure band gap of PbSe ~0.26–0.28 eV computationally and ~0.27–0.31 eV at ambient temperature, experimentally [53–55]. Thus, the present zero pressure value, that is, 0.185 eV of PbSe, is lower than these values.Figure 5
Variation of the energy band gap of PbSe under pressure. The solid line shows the linear fit and the dots show the calculated values.


## 3.1. Phase Transition and Structural Parameters

The structural parameters were theoretically evaluated first by calculating the total lattice energies over a set of different volumes around the equilibrium values and then were fitted to the third order Birch-Murnaghan equation of state (EOS) using the following relations:(2)
E
V
=
E
0
+
9
V
0
B
0
16
V
0
V
2
/
3
-
1
3
B
0
′
+
V
0
V
2
/
3
-
1
2
6
-
4
V
0
V
2
/
3
,
where E
0 is minimum energy, V
0 is corresponding volume, and B
0 is the bulk modulus at zero pressure and B
0
′
=
d
B
0
/
d
P. The third order Birch-Murnaghan isothermal equation of state [41, 42] relating pressure and volume is given by
(3)
P
V
=
3
B
0
2
V
0
V
7
/
3
-
V
0
V
5
/
3
×
1
+
3
4
B
0
′
-
4
V
0
V
2
/
3
-
1
.
Calculations were performed in both GGA and LDA frameworks and the respective parameters, that is, lattice constant (a), bulk modulus (B
0), and pressure derivative of bulk modulus (B
0
′) for B1, B27, and B2 phases of PbSe, are given in Table 1. As per the device application, high bulk modulus and its pressure derivative are essential parameters for photovoltaic applications in optoelectronic devices [43, 44]. The bulk modulus also defines resistance to change in volume against any mechanical deformation [44]. Therefore, we computed the values of bulk modulus and its derivative using Birch-Murnaghan EOS and summarized them in Table 1. The value of bulk modulus obtained from Birch-Murnaghan EOS agrees well with the value reported by Rached et al. [45]. This variation of B
0 is in accordance with the theory of Cohen [46] who obtained an empirical expression for the bulk modulus based on the nearest-neighbor distance. As seen from Table 1, the pressure derivative of bulk modulus (B
0
′) is almost the same for all the three phases with average value of 4.0. On comparing the structural parameters obtained from two different schemes, that is, GGA and LDA, we observed that lattice parameter obtained from LDA approach is larger than GGA while the bulk modulus obtained from GGA is higher as compared to the one obtained from LDA.Table 1
Calculated equilibrium lattice parameter (a), bulk modulus (B
0), pressure derivatives of bulk modulus (B
0
′), and cohesive energy (E
coh) for PbSe.
Parameters
Present calculations
Experimental results
Other calculations
GGA
LDA
NaCl-type, B1 phase
a
0 (Å)
6.08
6.17
6.117 [29], 6.133 [9], 6.124 [56], 6.117 [57]
6.211 [10], 6.098 [23], 6.226 [12], 6.214 [12]
B
0
(GPa)
66.96
64.72
28.8 [8], 54.1 [56], 54.1 [58]
48.27 [10], 60.8 [23], 44.5 [12], 49.5 [12], 66.0 [42]
B
0
′
4
3.7
4.1 [8]
4.35 [10], 4.56 [23], 3.64 [12], 3.8 [12]
E
coh (eV)
6.95
—
—
7.39 [44]
Orthorhombic Pnma phase (B27, FeB type)
a
0 (Å)
a = 11.38
a = 11.48,
a = 11.61,
a = 8.07,
b = 4.374
b = 4.474,
b = 4.00,
b = 4.49,
c = 4.259
c = 4.359
c = 4.39 [29]
c = 6.29 [10]
B
0 (GPa)
68.94
65.5
—
—
B
0
′
4.9
4.9
—
—
CsCl-type, B2 phase
a
0 (Å)
3.709
3.74
3.379 [9]
—
B
0 (GPa)
70
64.72
—
62.0 [12]
B
0
′
3.9
3.7
—
2.40 [12]Figures1(a) and 1(b) show the variation of the total energy versus volume computed for PbSe in B1, B27, and B2 phases obtained from both GGA and LDA approximations together with EOS data. As evident from Figures 1(a) and 1(b), the energy state of B1 phase is lower than that of the other two phases. This phase could be regarded as a stable phase of PbSe system [26].The calculated energy versus volume for B1, B27, and B2 phases of PbSe: using the PBE (a) and VBH (b) correlation functional. The scattered points show calculated energies and the solid lines show the fittedE
(
V
) curves according to Birch-Murnaghan equation of state (EOS).
(a)
(b)In order to investigate the phase transition pressure (P
t) for intermediate and CsCl-type phase, we computed variation in enthalpy with pressure using GGA and LDA approximations considering certain constraint of temperature for all the three phases of PbSe compound. Since a stable phase of any crystal is coupled with its minimum free energy, we have constrained the temperature to zero in order to obtain a thermodynamically stable phase which has lowest enthalpy,H =E + PV, under given pressure [26, 47]. At P
t enthalpy of the two consecutive phases is equal [16, 48]. Difference in enthalpy of high pressure phases (B27 and B2) and that in stable phase (B1) for both GGA and LDA approximations are plotted in Figures 2(a) and 2(b), respectively. From GGA approach, the transition from the stable B1 phase to Pnma intermediate phase is predicted to be 6.24 GPa and the transition from intermediate to CsCl-type phase is 16.39 GPa, while the P
t evaluated from LDA approach is 7.97 GPa for transition from NaCl phase to an orthorhombic Pnma (B27, FeB type) phase and 18.25 GPa for transition from Pnma (B27) phase to the CsCl-type (B2) phase.Enthalpy versus pressure calculated from first-principles total energy for the B1, B27, and B2 structures of PbSe using the PBE (a) and VBH (b) correlation functional.
(a)
(b)These investigated transition pressures are reported in Table2. Since the transition pressures derived from GGA approximation are closer to the experimental results as compared with the transition pressures derived from LDA approximation, for this reason further investigation is done with GGA approximation. The recent investigations of Li et al. [49] on the pressure-induced structure transitions of PbSe using density functional theory are in agreement with the present data. The calculation of cohesive energy from the experimental Compton profile is a difficult task due to weighting of p
2 in the high momentum region [50, 51]. Therefore, we have used the DFT-PBE profile for the calculation of cohesive energy of PbSe. The calculated cohesive energy along with available data are given in Table 1. From this table, it is evident that our calculated value of E
coh shows fairly good agreement with earlier data [44].Table 2
Transition pressure (in GPa) for PbSe.
Transition
Present
Experimental
Other calculations
GGA
LDA
B1 → B27
6.24
7.97
4.5 [9], 4.5 [8]
7.45 [10], 10.0 [11]
B27 → B2
16.39
18.25
16.0 [9], 16.0 [8]
18.76 [10], 20.0 [11]The pressure versus volume data for B1, B27, and B2 phases were computed, using the calculated values ofB
0 and B
0
′, using (3), and plotted in Figure 3. The curve shows that at P
t there is sudden collapse in volume.Figure 3
Variation of volume with increasing pressure for PbSe using GGA.


## 3.2. Pressure Dependent Electronic Properties

Knowledge of energy band gaps at different pressure is important to provide the information about their potential efficiency in fabricating various optoelectronic and thermoelectronic devices. In order to compute the energy band gap at different pressures, the equilibrium lattice constants are used and different symmetry points likeW (0.5, 0.25, 0.75), L (0.5, 0.5, 0.5), Γ (0.0, 0.0, 0.0), X (0.5, 0.0, 0.0), W (0.5, 0.25, 0.75), and K (0.375, 0.375, 0.75) are considered in the BZ. The energy band gaps of PbSe in B1 phase at different hydrostatics pressures are reported in Table 3.Table 3
Band gap of PbSe in B1 phase at different pressures.
Pressure (GPa)
Present calculation
Other calculations
−4
0.338
—
−2
0.252
—
0
0.185
0.21 [13], 0.15 [12], 0.16 [52]
2
0.118
—
4
0.05
—As seen from Figure4, the valence band maximum and the conduction band minimum lie at L point of the BZ. The occupied 6s-band of Pb atom lies close to the valence band of Se atom which is of p-type. Since the symmetry of Pb 6s band is the same as that of the highest valence band at the L point, there exists strong level repulsion which determines the feature of band structure and gives direct band gap along L point of the BZ at all considered pressures [12].Figure 4
Energy band structure of PbSe in B1 phase.Since the band structures at all considered pressures are similar with only a small change around the Fermi level, it is desirable to present band structure of B1 phase only at zero pressure (Figure4). The variation of band gap with pressure for B1 phase is plotted in Figure 5. It shows that with increasing pressure minima of conduction band come closer to the Fermi energy level and thus reduce the band gap. A similar feature of energy gap with pressure is reported by Nabi et al. [52]. A number of groups have observed the zero pressure band gap of PbSe ~0.26–0.28 eV computationally and ~0.27–0.31 eV at ambient temperature, experimentally [53–55]. Thus, the present zero pressure value, that is, 0.185 eV of PbSe, is lower than these values.Figure 5
Variation of the energy band gap of PbSe under pressure. The solid line shows the linear fit and the dots show the calculated values.


## 4. Conclusions

The pressure-induced phase transition and electronic properties of PbSe investigated using LCAO method were found to be in good agreement with previous theoretical and experimental results. The intermediate phase of PbSe between B1 and B2 phases is found to be B27 phase. The variation of energy band gap of PbSe under pressure is reported and it is found that, unlike II–VI and III–V compounds, PbSe shows direct band gap atL point of BZ.



---

*Source: 921092-2014-12-30.xml*