Parrydeep Kaur
Sachdeva
abc,
Shuchi
Gupta
b and
Chandan
Bera
*a
aInstitute of Nano Science and Technology, Knowledge City, Sector-81, S. A. S Nagar, Mohali, Punjab 140306, India. E-mail: chandan@inst.ac.in
bUniversity Institute of Engineering and Technology, Panjab University, Sector-25, Chandigarh 160014, India
cDepartment of Physics, Panjab University, Sector-14, Chandigarh 160014, India
First published on 7th April 2021
Two dimensional (2D) chalcogenide monolayers have diversified applications in optoelectronics, piezotronics, sensors and energy harvesting. The group-IV tellurene monolayer is one such emerging material in the 2D family owing to its piezoelectric, thermoelectric and optoelectronic properties. In this paper, the mechanical and piezoelectric properties of 2D tellurene in centrosymmetric β and non-centrosymmetric β′ phases are investigated using density functional theory. β′-Te has shown a negative Poisson's ratio of −0.024 along the zigzag direction. Giant in-plane piezoelectric coefficients of −83.89 × 10−10 C m−1 and −42.58 × 10−10 C m−1 are observed for β′-Te under biaxial and uniaxial strains, respectively. The predicted values are remarkably higher, that is 23 and 12 times the piezoelectric coefficient of a MoS2 monolayer with biaxial and uniaxial strain in the zigzag direction, respectively. A large thermal expansion coefficient of tellurene is also estimated using quasi harmonic approximation. High piezoelectricity combined with exotic mechanical and thermal properties makes tellurene a very promising candidate in nanoelectronics.
Another emergent characteristic observed in monoelemental 2D materials such as graphene20 and phosphorene21,22 is negative Poisson's ratio (NPR). In-plane NPR emerges in graphene while out-of-plane NPR exists in black phosphorus. The difference lies in the mechanism for the structural compensation of the strain. NPR has been observed in a few 2D materials having honeycomb lattices such as silicene, h-BN, h-GaN, h-SiC, and h-BAs.23 A few 2D heterostructures constructed with graphene/MoS2 and graphene/h-BN24 also show negative Poisson's ratio. Initially, the NPR behaviour was explored in systems such as crystalline SiO2,25 metallic nanoplates,26 2D isotropic systems27 and chiral honeycomb systems.28 These materials with a negative Poisson's ratio have many promising applications in defense, biomedical, etc.21,29 Mechanical metamaterials have superior physical properties compared to conventional materials with a positive Poisson's ratio.
Preferable thermal properties when combined with their mechanical properties become very suitable for device applications. Among 2D chalcogenides, graphene displays negative thermal expansion in the temperature range 0–1000 K.30,31 Doping graphene with nitrogen and boron makes the thermal expansion coefficient more negative. Negative thermal expansion is caused by the transverse acoustic modes present in graphene, which result in highly negative Grüneisen parameter. Simulations by Sevik30 revealed that h-BN exhibits a high negative thermal expansion (−6.5 × 10−6 K−1) compared to graphene (−3.6 × 10−6 K−1). Huang et al.32 reported negative to high positive variation of thermal expansion with temperature in MoS2. A similar variation in thermal expansion has been observed for MoTe2.33 Another TMDC, 2H-WS2 also shows negative thermal expansion up to 1100 K.34 A negative volumetric expansion exists in GeTe depending upon the thermal evolution of its structure. Herein, the coupling of acoustic and soft transverse optical modes causes the negative thermal expansion.35 Black phosphorus possesses anisotropic thermal expansion with negative values in the zigzag direction.36
Since the mechanical exfoliation of graphene,37 researchers from all over the globe have come to terms with the possibility of extracting 2D materials from their bulk counterparts. Ultrathin Te films have already been synthesized using different approaches such as by extracting monolayers from layered bulk tellurium via the liquid-phase exfoliation method or by growing Te monolayers on a substrate using physical vapor deposition, epitaxial formation by molecular beam epitaxy and synthesis of free standing 2D Te by the solution method.38–43 Wu et al.44 constructed a 2D FET device with a large drain current of 300 mA mm−1, high mobility and high on/off ratio. Te has superior transport properties that make it suitable for thermoelectric applications.41,43,45,46 It has high carrier mobilities leading to good electrical conductivity and the heavy Te atoms result in flat bands towards the conduction region which increase the Seebeck coefficient.47 The unusually low thermal conductivity in tellurene is a major contributing factor towards increasing its power factor and efficiency in a thermoelectric device.48,49 Also, theoretical reports have shown that Te has a direct band gap,50,51 and hence, optical to electrical energy can be harnessed from it.41,42,52 Te is fundamentally centrosymmetric but becomes non-centrosymmetric upon strain engineering. Also, strain is an unavoidable entity during the fabrication of thin films. Nowadays, strain is also intentionally used to modulate the properties.
In this article, we have investigated the mechanical and piezoelectric properties of single layer tellurene under strain along with its temperature dependent expansion. The motivation of this work lies in a few previous reports indicating the origin of ferroelectricity and polarizations in monoelemental tellurene.53,54 Two stable structures of 2D Te are obtained in this study; where one structure is centrosymmetric and the other one is non-centrosymmetric. Anomalous piezoelectric coefficients and Poisson's ratio are observed in the non-centrosymmetric 2D Te.
β-Te (Fig. 1(b)) shows an indirect 1.26 eV band gap in the Γ → BA direction with the PBE functional. The inclusion of spin–orbit coupling changes the band gap to 0.96 eV and makes it direct. The band gap calculated using the HSE functional is 1.81 eV and is indirect. In a similar manner, the incorporation of SOC makes the band gap direct and reduces it to 1.48 eV. The calculated band gaps for β-Te using different functionals are in reasonable agreement with previous theoretical reports.50,52,65 On the other hand, the β′-Te monolayer (Fig. 1(d)) displays indirect band gaps of 1.18 eV and 1.13 eV with PBE and PBE + SOC, respectively. The HSE and HSE + SOC functionals also result in indirect band gaps of 1.89 eV and 1.74 eV, larger than the ones obtained using the PBE functional. The band structures obtained using HSE with and without spin–orbit coupling are shown explicitly in the ESI, Fig. S5.†
The mechanical strength of 2D structures to withstand strain can be measured by computing their Poisson's ratio. In the case of a monolayer, the in-plane Poisson's ratio is the most significant component. We investigated this quantity for single layer β-Te and β′-Te by deforming them uniaxially in the zigzag and armchair directions with the strain in the range of −10% to +10%. Once the structures get relaxed and optimized, energy–strain relations are obtained (ESI, Fig. S1 and S2†).
Poisson's ratio is computed as , where εa is the strain generated along the a-axis in response to the applied strain εb along the b-axis or vice versa. β-Te displays positive Poisson's ratios (PPR) of 0.49 and 0.29 for uniaxial loading in the zigzag and armchair directions as shown in the ESI (Fig. S3).† The Poisson's ratio values for β-Te perfectly align with Gao et al.'s reported values.48 However, there is a small deviation in the Poisson's ratio in the zigzag direction (0.58) calculated by Liu et al.66
The Poisson's ratio of the β′-Te monolayer is shown in Fig. 2(a) as a function of uniaxial deformation along the zigzag direction (b-axis) and Fig. 2(b) presents the Poisson's ratio for deformation along the armchair direction (a-axis). In Fig. 2(a), the data are fitted to the 7th order function, y = −v1x + v2x2 + v3x3 + v4x4 + v5x5 + v6x6 + v7x7 + c, which gives the linear Poisson's ratio v1 = 0.01 and constant, c = −5.21. In Fig. 2(b), the calculated data are fitted to the function, y = −v1x + v2x2 + v3x3 + v4x4 + v5x5 + c. The linear Poisson's ratio obtained is v1 = 0.29 and c = −7.57.
The εavs. εb plot (Fig. 2(a)) shows that tellurene expands (or compressive strain decreases) in the a-axis when compressive strain along the b-axis is applied. This is referred to as positive Poisson's ratio. Tellurene also expands in the a-axis under applied tensile strain along the b-axis. When a material elongates in an axis upon stretching its other axis or vice versa, it is known as negative Poisson's ratio behavior. Therefore, uniaxial deformation in tellurene leads to positive Poisson's ratio upon applying negative strain and negative Poisson's ratio upon applying positive strain. Materials exhibiting negative Poisson's ratio are called auxetic materials due to their unusual mechanical behaviour.67 NPR in tellurene is attributed to the existence of a six-chair like arrangement of Te atoms in the puckered monolayer. The structural evolution for the in-plane NPR behavior is discussed in the ESI, Table S1.†
Quantitatively, a high PPR of 0.48 is obtained in the negative strain range which is close to the critical value, 0.5, for most of the materials.5,68 For the applied positive strain, a linear negative Poisson's ratio equivalent to −0.024 is achieved. The linear NPR is quite small and this indicates that the tensile stress leads to a small lateral strain only. 1 T-MX2 (M = Mo, W, Tc, Re; X = S, Se, Te) monolayers also display small in-plane negative Poisson's ratios ranging from −0.03 to −0.37.67 Borophene also shows a NPR of −0.04 along a and −0.02 along b.69 δ-phosphorene has a NPR of −0.267.70 Black phosphorus exhibits a theoretical out-of-plane NPR of −0.027,21 which has been experimentally validated.71 Silica in its different 2D phases is reported to have a negative Poisson's ratio in the range of −0.02 to −0.12.72 Other puckered 2D materials having a NPR include group V-enes (BP, As, and Sb) and group IV monochalcogenides (SiS, SiSe, GeS, GeSe, SnS, SnSe, and SnTe).73
Now, the effect of applied strain in the a-axis of β′-Te is observed on the b-axis (Fig. 2(b)). The application of strain in the armchair direction (a-axis) results in a strain in the zigzag direction (b-axis) in a way that only a positive Poisson's ratio is obtained for both the positive and negative applied strains. A linear positive Poisson's ratio of 0.29 is obtained from the εb vs. εa plot. β′-Te exhibits different Poisson's ratios for strain applied to the armchair and zigzag directions. This substantial anisotropy in the in-plane mechanical properties is related to the geometric anisotropy in the monolayer (Fig. 1(a)). The obtained PPR is comparable to that of 2D materials such as graphene (v = 0.18), BN (v = 0.22), BP (v = 0.27),74 MoS2 (v = 0.30),4 ZrS2 (v = 0.19), ZrSe2 (v = 0.20), HfS2 (v = 0.18) and HfSe2 (v = 0.21).5
Mechanical strength of a nanomaterial can also be determined by its 2D elastic constants. The following conditions govern the mechanical stability of tellurene monolayers.75
C11C22 − C122 > 0 and C11, C22, C66 > 0 | (1) |
![]() | (2) |
The above conditions (eqn. (1)) are satisfied by both β-Te and β′-Te. By definition, elastic constant is , where σij and εkl are the stress and strain tensors, respectively. In Voigt notation, the fourth order tensor has been reduced to the second order tensor as seen in eqn (2). For a highly stretchable material, elastic constants should be small and hence, it should have a small Young's modulus. This is preferable for easy displacement of the charges in 2D sheets to generate piezoelectricity. The anisotropic 2D Young's modulus for monoclinic tellurene is given by66
![]() | (3) |
Table 1 presents the elastic constants and Young's modulus of β-Te and β′-Te. The calculated elastic constants and Young's modulus of β-Te are close to the ones reported by Liu et al. at 0 K temperature.66 Elastic constants C22 and C12 of β′-Te are smaller than the constants observed for β-Te. This implies that β′-Te is softer than β-Te in the b-direction whereas β-Te is more flexible in the a-direction. Also, the Young's modulus of β′-Te along the b-axis (Y22) is much smaller than that along the a-axis (Y11), indicating that β′-Te is softer in the zigzag direction. Tellurene is a low stiffness material as it has elastic constants and Young's modulus smaller than those of many 2D materials such as graphene (C11 = 351.9 N m−1, Y = 341 N m−1),74 BN (C11 = 290 N m−1, Y = 275.9 N m−1),74 BP (C11 = 145 N m−1, Y = 135.6 N m−1),74 MoS2 (C11 = 134.9 N m−1, Y = 125 N m−1),4 MoTe2 (C11 = 88 N m−1, Y = 79.9 N m−1),9,76 ZrS2 (C11 = 75.57 N m−1, Y = 72.55 N m−1), ZrSe2 (C11 = 68.86 N m−1, Y = 65.83 N m−1), HfS2 (C11 = 80.77 N m−1, Y = 77.99 N m−1) and HfSe2 (C11 = 75.07 N m−1, Y = 71.85 N m−1).5 Elastic constants of SnS, SnSe, GeS and GeSe12 are also larger than those of tellurene.
Monolayer | C 11 | C 22 | C 12 | Y 11 | Y 22 |
---|---|---|---|---|---|
β-Te | 13.52 | 23.68 | 6.69 | 11.63 | 20.37 |
β′-Te | 16.42 | 3.27 | 2.01 | 15.18 | 3.02 |
Further, we have investigated the piezoelectric properties of β-Te and β′-Te. Also, the effect of strain is observed on their piezoelectric coefficients. The stress tensor is defined as
![]() | (4) |
![]() | (5) |
In two dimension, the number of independent elements reduces further due to increased symmetry. Therefore, elements e34 and e35 become zero. As the stress coefficients are evaluated directly using DFPT, the piezoelectric strain coefficients dij can be determined from the stress coefficients and elastic constants Cij by the following relation,
![]() | (6) |
On solving the equations obtained from this relation (ESI†), the values of dij are obtained, which are experimentally relevant.
β-tellurene, being centrosymmetric, does not exhibit piezoelectricity. However, it shows piezoelectricity under the effect of strain. The phase transition to a P2 structure is obtained at 4.5% tensile biaxial strain resulting in high piezoelectric coefficients (Table 2). The complete variation of piezoelectricity under strain is shown in the ESI (Fig. S4†). Uniaxial strain has no effect on its inversion symmetry, thus leading to zero piezoelectricity. The piezoelectric stress coefficients of single layer non-centrosymmetric β′-tellurene have also been calculated and are tabulated in Table 2. It displays quite large coefficients than β-Te. Fig. 3(a)–(c) show the variation of all the elements of the stress tensor under the applied biaxial strain. Strain led to an enhancement in all the components of the piezoelectric tensor resulting in peaks for e21, e22, e23 and e26 at a compressive biaxial strain of −5.5% (Table 2). The shear strain components e14 and e15 increase with the increase in tensile strain. The highest element obtained is the in-plane polarization (e22) having a value of −83.89 × 10−10 C m−1 (d22 = −333.14 pm V−1). For all elements, the ionic contribution is higher than the electronic contribution, highlighting the ionic nature of bonds between Te atoms.
Monolayer | Strain | Ionic | Electronic | Net contribution | Net contribution | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
e 21 | e 22 | e 23 | e 26 | e 21 | e 22 | e 23 | e 26 | e 21 | e 22 | e 23 | e 26 | d 21 | d 22 | d 23 | d 26 | ||
β-Te | Biaxial strain (4.5%) | −3.88 | −10.98 | −0.90 | 0.29 | 1.36 | −1.78 | 0.81 | −0.52 | −2.52 | −12.76 | −0.09 | −0.23 | 36.69 | −310.32 | −15.18 | −5.14 |
β′-Te | Strain (0%) | −3.12 | −7.44 | −0.78 | 0.31 | 1.35 | −1.81 | 0.73 | 0.47 | −1.77 | −9.25 | −0.05 | 0.78 | 25.79 | −298.78 | 12.12 | 15.54 |
Biaxial strain (−5.5%) | −36.68 | −83.72 | −3.81 | −6.30 | 0.12 | −0.17 | 0.13 | −0.04 | −36.56 | −83.89 | −3.68 | −6.34 | −106.94 | −333.14 | −95.15 | −97.52 | |
Uniaxial strain along zigzag (−4.9%) | −17.79 | −41.82 | −1.25 | −0.71 | 0.56 | −0.76 | 0.59 | −0.23 | −17.23 | −42.58 | −0.66 | −0.94 | −112.48 | −302.22 | −30.37 | −13.77 |
For uniaxial strain in the zigzag direction (Fig. 3(d)–(f)), the variation of piezoelectric coefficients is similar to the variation under biaxial strain although the magnitude is different. Again, in-plane polarization (−42.58 × 10−10 C m−1) (d22 = −302.22 pm V−1) is the highest for uniaxial deformation in the b-axis. Also, we observed a phase change from monoclinic P2 to P2/m (denotes the centrosymmetry) below biaxial and uniaxial compressive strains −5.5% and −4.9%, respectively. We explore the variation in piezoelectric coefficients by deforming the structures under the application of uniaxial strain along the armchair direction also. However, it is observed that these structures attain the centrosymmetric nature in the whole range of −10% to +10% and hence, no piezoelectricity is obtained along the a-axis.
The in-plane piezoelectric coefficient of tellurene under biaxial and uniaxial strain in the zigzag direction is higher than that of MoS2 by a factor of 23 and 12, respectively. This in-plane piezoelectric coefficient is remarkably huge and hence tellurene proves to be a great replacement for the existing traditional piezoelectric materials such as AlN (5.6 pm V−1),78 GaN (3.1 pm V−1),78 wurtzite-BN (0.60 pm V−1), MoS2 (3.73 pm V−1) and quartz (2.3 pm V−1).1 The coefficients are found to be larger than the reported giant piezoelectricity in SnS (144.76 pm V−1), SnSe (250.58 pm V−1), GeS (75.43 pm V−1), and GeSe (212.13).12
Spontaneous polarization is further analyzed using Born effective charges occurring in the 2D sheet. BEC determine the change in polarization (Pi) caused by applied strain or the force (Fi) experienced by atoms due to the applied electric field (εj) as
, where Ω is the volume and uj is the displacement in the system.
The diagonal dynamic charges on the β′-Te atoms are given in Table S2.† We find that for biaxial strain, charge on Te5 increases by 53% along the armchair direction and by 105% along the zigzag direction while Te2 and Te3 show a reduction in charges by the same amounts in the respective directions. Under a uniaxial strain of −4.9% in the zigzag direction, the increase and decrease of dynamic charges are nearly 48% and 93% for respective Te atoms. Dynamic charges developed along the z-direction are too small indicating no out-of-plane polarization. Typically, the materials having only one element are non-polar. But since Te exhibits multivalent nature and due to the different coordination of Te atoms (Fig. 1(a)) in this particular monolayer, there are different charges on the atoms. BECs indicate that Te5 acts as a donor while Te2 and Te3 atoms act as acceptors. Hence, the origin of the huge piezoelectricity is attributed to the relatively large BECs of Te atoms.
As the mechanical properties of tellurene ascertain that it is highly stretchable, it becomes important to determine the temperature dependent behavior of sheets. We calculated the thermal expansion of the Te monolayer using the quasi harmonic approximation. Both β and β′-Te are strained from −5% to +5% and phonon frequencies are calculated for all 11 structures. The equilibrium volume V(T) is determined from the free energies F(V, T) calculated at each temperature up to 400 K. The volume expansion and variation of coefficient of thermal expansion (α) with temperature are shown in Fig. 4. It is expressed as
![]() | (7) |
![]() | ||
Fig. 4 Temperature dependence of (a) volume and (b) coefficient of thermal expansion of β and β′ tellurene. |
Liu et al.66 have shown anisotropic thermal expansion along the armchair and zigzag directions in β-Te where the negative expansion occurs particularly along the a-axis in the low temperature range. They have also observed large thermal expansion along the b-axis. We have obtained a similar thermal expansion of β-Te lying in the same range. The thermal expansion attained in β′-Te is higher than that in β-Te. At low temperatures, the negative thermal expansion in β-Te and β′-Te is possibly because of the negative values of Grüneisen parameter.66 Grüneisen parameter describes the change in vibrational properties with the change in volume. The variation of Grüneisen parameter with respect to temperature is shown in Fig. S6.†
In comparison to MoS2,32 tellurene is highly expandable and has thermal expansion higher by nearly 4 times for β-Te and 7 times for β′-Te at room temperature, 300 K. Also, expansion in tellurene is larger than in black phosphorene (9 ∼10 × 10−6 K−1)79 and MoSe2 (8 × 10−6 K−1).30 Along with the large thermal expansion at room temperature, tellurene shows negative thermal expansion from 0 to 32 K. There are a few 2D materials that show negative thermal expansion in a broad range of temperature. This includes graphene,30,31 doped graphene31 and h-BN30 exhibiting negative thermal expansion up to 1000 K and black phosphorene displaying negative expansion up to 150 K.22 GeTe also falls into the category of materials having negative thermal expansion.35
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0na00930j |
This journal is © The Royal Society of Chemistry 2021 |