Research Article
Lattice Instability of 2H-TaS2
Department of Physics, St. Joseph�s College, Tiruchirappalli, 620 002, India
N. Lawrence
Department of Physics, St. Joseph�s College, Tiruchirappalli, 620 002, India
Charge Density Wave (CDW) is a dramatic phenomenon where coupling of the electron-phonon works with the Fermi surface nesting results in the spatial modulation of the crystalline lattice and conduction-electron density. The formation of CDWs in one-dimensional systems was first suggested by Peierls (1995) and Frohlich (1954). However, at low temperatures the lattice will not be stable. The instability is due to the periodic potential presence which will break the distribution of Fermi surface which induces charge density wave. The 2H-TaS2 has hexagonal structure and has interesting phase change behaviour. It is in the ICCDW phase below 75 K and in normal phase above this transition temperature (Sugai, 1985). In ICCDW phase, it has unit cell of size 3a0x3a0xco, where a0 and c0 are the lattice constants.
Many experimental studies (Brouet et al., 2004; Shin et al., 2005; Lavagnini et al., 2008) are available to study the characteristics of CDW system such as 2H-TaS2. The available theoretical model describes the electron-phonon interaction in the CDW formation. But a few experimental studies (Harper et al., 1977; Sugai et al., 1981; Wilson et al., 2001) are available to investigate the variation of thermal properties with temperature in the context of CDW. Hence, it has been planned to investigate the dynamical behaviour of the lattice and the contribution of phonons towards the thermal properties of this system. In this study, phonon frequency distribution, thermal properties includes Debye Waller factor and specific heat capacity are calculated and in normal and ICCDW phases the thermal conductivity for 2H-TaS2 are reported (Balaguru et al., 2002).
CALCULATION PROCEDURE
The space group of 2H-TaS2 is D46h. It has two layers and six atoms per unit cell namely two Ta atoms (Type 1 and 2) and four S atoms (Type 3, 4, 5 and 6). In Cartesian frame, the basic vectors of unit are given by:
and:
The unit vectors of the reciprocal Lattice are given by:
and:
In lower layer, the Ta1, two S3 and S4 are denoted by:
and:
respectively. In upper layer, the coordinates of Ta2, two S5 and S6 are:
and:
respectively, where a, c and u are lattice parameters having values 3.35 Å, 12.10 Å and 0.118, respectively.
Born-Von Karman model has been used to illustrate the force constant tensor of the Ta-S interaction. The lower layer force constant is represented by:
and the upper layer is represented by six parameters:
using the symmetry property the force constant tensor can be created for the remaining S neighbors. For the S-S interactions of same type, the force constant parameters are given by A5 and B5. A3 and B3 denotes the force constant parameter for Ta-Ta interactions of same type and A4 and B4 represents the Ta-Ta interactions of different type. Pawleys potential (Rinaldi and Pawley, 1975) will fix the values for these parameters which will be in units of 104 dyne cm-1 as:
The dynamical matrix represents Hermitian matrix of (18x18) (as shown in Appendix 1) will be diffused into (36x36) matrix. The diagonalization of this real symmetric matrix results in phonon frequencies and polarization vectors.
In normal phase, by uniformly dividing the Brillouin zone, phonon frequencies for a set of 84 wave-vector points (Di Salvo et al., 1976). In ICCDW phase, to calculate phonon frequencies for a set of 336 wave vector points a folding technique has been used. Using values and vectors, Debye Waller factor and thermal conductivity of 2H-TaS2 t in the variation of specific heat capacity with reference to temperature in both the phases are reported by following the previous study (Berman, 1976; Nunez-Regueiro et al., 1985).
Figure 1 represents phonon frequency distribution curves. The blank circled C line and the black circled B curve represent the normal phase and ICCDW phase, respectively. From the figures, most of the frequency modes are distributed in low frequency of the ICCDW as compared to normal one which results in less electron scattering by phonons in the ICCDW phase. In addition to the beginning of CDW and because of the less electron phonon scattering, the resistivity is expected to decrease in the ICCDW phase. Hence, the present result supports the superconducting behavior of this material at low temperature as observed by Smith et al. (1975).
Figure 2 represents the 0-265 k variation in specific heat. Since, the phase transition occurs at 75 K, in ICCDW phase the frequency distribution data has been used below this transition and in normal phase above the transition phase. From the graph, discontinuity occurs at the transition temperature.
Fig. 1: | Phonon frequency distribution |
Fig. 2: | Specific heat variation with temperature |
Fig. 3: | Debye waller factor of tantalum atoms |
At the transition temperature, this discontinuity may be attributed to the change of the distributed phonon frequency. The observed results are comparable with Craven and Meyer (1977).
Using the respective Eigen values and vectors for both the phases, Ta and S Debye Waller factor values at different temperatures have been calculated. The curves for both the types of atoms are shown in Fig. 3 and 4, respectively. At the transition temperature, an anomalous behavior is observed from these curves that may be due to the electron-phonon interaction change in the ICCDW phase and due to the increase in the size of the unit cell. From these figures, Ta atoms are displaced lesser than chalcogen atoms due to the heavier nature of Ta atom and is comparable with the results of Moncton et al. (1977).
For the range of 25-175 K, the thermal conductivity and temperature curve are shown in Fig. 5 that shows a discontinuity at 75 K, the transition temperature as expected and also results in a drastic decrease that occurs due to the heavy change in the concentration of electrons and mainly due to phonon frequency modes in larger number in ICCDW phase.
Fig. 4: | Debye waller factor of sulfur atoms |
Fig. 5: | Thermal conductivity of 2H-TaS2 |
Thus, the lattice vibration forms the contributions to the thermal conductivity and hence, the phonon frequency modes. This curve is similar to the curve observed for 2H-NbSe2.
The phonon frequencies of the normal phase for 2H-TaS2 have been computed using 84 wave vector points obtained by uniformly dividing the Brillouin zone. For ICCDW phase, they were computed using 336 wave vector points generated using folding technique. The phonon frequency distribution curves have been drawn and using them, the thermal properties such as the specific heat capacity, the Debye Waller factor for both atoms (Ta and S) and the thermal conductivity for both the phases have been computed. The present results are in good agreement with the experimental results.
Appendix 1: | From the diagonalization of the dynamical matrix, the phonon frequencies and their polarization vectors can be obtained. Maradudin et al. (1971) having elements as follows: |
where, Q1, Q2 and Q3 are the phonon wave vector components along b1, b2 and b3. Fig. 1: Phonon frequency distribution.