Anisotropic Thermoelectric Response in Two-Dimensional Puckered

Aug 8, 2016 - †Institute for Materials Science and Max Bergmann Center of Biomaterials, ‡Institut für Physikalische Chemie und Elektrochemie, §D...
5 downloads 14 Views 2MB Size
Article pubs.acs.org/JPCC

Anisotropic Thermoelectric Response in Two-Dimensional Puckered Structures Leonardo Medrano Sandonas,†,⊥ David Teich,‡ Rafael Gutierrez,*,† Tommy Lorenz,‡ Alessandro Pecchia,# Gotthard Seifert,‡ and Gianaurelio Cuniberti†,§,∥ †

Institute for Materials Science and Max Bergmann Center of Biomaterials, ‡Institut für Physikalische Chemie und Elektrochemie, Dresden Center for Computational Materials Science (DCMS), and ∥Center for Advancing Electronics Dresden, TU Dresden, 01062 Dresden, Germany ⊥ Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany # Consiglio Nazionale delle Ricerche, ISMN, Via Salaria km 29.6, 00017 Monterotondo, Rome, Italy §

S Supporting Information *

ABSTRACT: Two-dimensional semiconductor materials with puckered structure offer a novel playground to implement nanoscale thermoelectric, electronic, and optoelectronic devices with improved functionality. Using a combination of approaches to compute the electronic and phonon band structures with Green’s function based transport techniques, we address the thermoelectric performance of phosphorene, arsenene, and SnS monolayers. In particular, we study the influence of anisotropy in the electronic and phononic transport properties and its impact on the thermoelectric figure of merit ZT. Our results show no strong electronic anisotropy, but a strong thermal one, the effect being most pronounced in the case of SnS monolayers. This material also displays the largest figure of merit at room temperature for both transport directions, zigzag (ZT ∼ 0.95) and armchair (ZT ∼ 1.6), thus hinting at the high potential of these new materials in thermoelectric applications.



black phosphorus.19,20 The main feature of phosphorene is its puckered honeycomb structure, which leads to electronically relevant properties such as a nonzero bandgap of ∼1.1 eV, large current on−off ratios, high charge mobility, and fast photoresponse.21,22 Moreover, anisotropic behavior of its physical properties has been demonstrated.23−29 In this respect, Luo et al. have reported anisotropic in-plane thermal conductivity of suspended few-layer black phosphorus measured by microRaman spectroscopy.23 Here, it was also found that, independent of the number of layers, the thermal conductivity for the zigzag (ZZ) direction is larger than for the armchair (AC) direction; for black phosphorus films thicker than 15 nm values of ∼40 and ∼20 W/mK for ZZ and AC directions, respectively, were reported. This effect has been theoretically verified by using the Boltzmann transport equation (BTE)24,25 and nonequilibrium Green’s function (NEGF) techniques.26 Electronic anisotropy has also been experimentally observed in layered black phosphorus.27,28 The results indicate that the AC direction is the most preferable electron transport direction. Taking advantage of these anisotropic properties, Fei et al. have shown, by performing first-principles band structure calcu-

INTRODUCTION Thermoelectric (TE) materials interconvert heat to electricity and are of great interest in environmentally relevant issues such as waste heat recovery and solid-state cooling.1 The efficiency of a thermoelectric material is measured by its figure of merit, ZT = S2GT/(κel + κph), where G is the electrical conductance, S is the Seebeck coefficient, T is the absolute temperature, and κel(ph) is the electronic (phononic) thermal conductance. The larger ZT is the better is the performance of a TE material. Thus, promising TE materials should simultaneously have a high power factor (defined as the product GS2) and low thermal conductance, ideally implementing an electron crystal and a phonon glass. Hence, in order for TE materials to become competitive with conventional refrigerators and generators, routes for tuning both transport quantities, electrical and thermal, need to be developed. Nowadays, the ongoing quest for highly efficient TE materials has been centered on nanostructured materials.2−10 In particular, two-dimensional (2D) materials have been pointed out to provide a new direction for designing highperformance TE devices due to the high sensitivity of their thermal and electronic properties to finite size quantum effects.11−18 Recently, phosphorene (or black phosphorene), one of the newest member from the family of 2D materials, has been successfully produced by mechanical exfoliation from bulk © 2016 American Chemical Society

Received: May 17, 2016 Revised: August 5, 2016 Published: August 8, 2016 18841

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C G = e 2L 0

lations, that the thermoelectric performance of phosphorene can indeed be enhanced for transport in the AC direction, with values of ZT ∼ 1.0 at room temperature.29 In the past few years, a new family of 2D materials such as arsenene (a single layer of arsenic) and monolayer tin sulfide (SnS) has been discovered. They all may exist, like in phosphorene, orthorhombic puckered structures and semiconducting bandgaps,30−38 and, more interestingly, they all display anisotropic features in their physical properties. This has been reported e.g. in arsenene, where a different behavior for the electronic bandgap was found when uniaxial strain was applied along the ZZ and AC directions.31 A recent theoretical study by Zeraati et al. in orthorhombic arsenene has also shown, by solving the BTE, a highly anisotropic thermal conductivity at room temperature where the ZZ direction is the preferred direction for thermal transport (∼30.4 W/mK).34 Unlike arsenene, SnS nanosheets have been successfully synthesized by several techniques.35,39,40 It has also been found that SnS thin films display a strong angle-dependent Raman response, indicative of anisotropy in their optical properties.39,40 The large Seebeck coefficient and the low thermal conductivity theoretically found in SnS bulk produce comparative thermoelectric performance compared to SnSe bulk,37,41,42 which has been proclaimed to be an excellent candidate for thermoelectric devices.43 It has also been reported that the electronic bandgap of SnS increases as the number of layers decreases,32,35 but it is not clear how this tendency will affect the thermoelectric response. Moreover, Guo et al. have shown that the thermal conductivity for SnS bulk is very low compared to other 2D materials and is also anisotropic.42 Based on the studies performed in other phosphorene analogues,29,34,38 it is expected that the anisotropic effect remains in its monolayer form. However, to the best of our knowledge, this statement has not been yet verified. Motivated by the results reported above, we analyze here the influence of anisotropy on the electron and phonon transport properties of homo and heteroatomic two-dimensional materials with puckered structuresphosphorene, arsenene, and SnS monolayerswith the goal of identifying strategies to achieve a good thermoelectric performance. To do this, we employ a density-functional based tight-binding (DFTB) approach combined with Green’s function techniques. Our results show that the thermoelectric response of the studied materials, quantified by the figure of merit ZT, is larger along the AC direction and that SnS monolayers display the highest ZT values (e.g., ∼1.6 at 300 K and ∼4.5 at 800 K) followed by arsenene. In the next section, we describe the computational methodology and, then, proceed to discuss the results of the anisotropic thermoelectric properties of these materials.

S=

(1)

1 L1 qT L0

κel =

(2)

L 2⎤ 1⎡ ⎢L 2 − 1 ⎥ T⎣ L0 ⎦

(3)

where q is the electric charge of carriers, which is positive for holes and negative for electrons. The functions Lm are given in terms of the moments of the electronic transmission function: Lm =

2 h

⎡ ∂f (E , μ , T ) ⎤ ⎥τel(E) dE ⎦ ∂E

∫ (E − μ)m ⎢⎣−

(4)

with f(E,μ,T) being the Fermi−Dirac distribution function, μ the chemical potential, and h Planck’s constant. For phonon transport, the corresponding linear thermal conductance is given by κ ph =

ℏ2 2πkBT 2

∫0



ω2

eℏω / kBT τph(ω) dω (eℏω / kBT − 1)2

(5)

with kB being the Boltzmann constant. The prefactor in front of the phonon transmission function is obtained by a linear expansion in the applied temperature difference ΔT of the quantity NB(T + ΔT) − NB(T), where NB is the Bose−Einstein distribution. Notice that while for electrons the relevant integration range is usually located around the Fermi energy of the system, computing the thermal conductance of the vibrational system requires an integration over the whole phonon spectrum. Both transmission coefficients τel(E) and τph(ω) can be calculated in a similar way by using Green’s function techniques.46,47 We use a common partitioning scheme for both electron and phonons by dividing the system into three regions, namely left electrode, right electrode, and a central region. Notice however that the atomic composition of all three regions is the same for each given material; i.e., there are no effects related to interface scattering, neither for electronic nor for phononic degrees of freedom. The electron transmission coefficient τel(E) is then calculated as τel(E) = Trace(Gr ΓLGa ΓR)

(6) r

where the retarded Green’s function G is obtained from the Hamiltonian H, Gr = (EI − H − ΣrL − ΣrR)−1, and the broadening functions ΓL/R = i[ΣrL/R − ΣaL/R] with ΣL/R being the electrode self-energies. The phonon transmission coefficient τph(ω) is calculated in a similar way by the substitution H → K and EI → ω2I, where K is the dynamical matrix. Here, both the electronic Hamiltonian matrix and the interatomic force constants are obtained using the DFTB method as implemented in the DFTB+ software package.48,49 The DFTB method is based on density-functional based parametrization of a general tight-binding Hamiltonian. Hereby, the system eigenstates are expanded into a linear combination of atomic orbitals (LCAO) minimal valence basis set, which are determined from self-consistent atomic LDA calculations. Additionally, the effective Kohn−Sham potential is approximated as a superposition of localized potentials of neutral atoms; see the Supporting Information for additional details. A valence basis set has been used including the 3s, 3p, and 3d orbitals for phosphorus51 and the 4s, 4p, and 4d orbitals for arsenic. For SnS we used the 5s, 5p, and 4d orbitals for tin and



COMPUTATIONAL METHODS Neglecting electron−phonon interactions, which would need a separate self-consistent nonequilibrium study, both electronic and phonon transport channels can be considered as independent from each other.6 Then, the electronic and phononic (ballistic) transport properties can be computed within the Landauer approach in terms of the corresponding electronic τel(E) and phononic τph(ω) transmission coefficients. Using nonequilibrium thermodynamics and series expansions of the electrical current in the low bias, linear response regime, closed expressions for the different electronic transport coefficients appearing in the definition of the figure of merit can be written as6,12,44,45 18842

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C

Figure 1. Electron and phonon band structures for homo- and heteroatomic puckered materials: (a) phosphorene, (b) arsenene, and (c) tin sulfide (SnS) monolayer. It is also shown the atomistic view of the 2D nanostructures, highlighting the zigzag (ZZ) and armchair (AC) transport direction. On the electron dispersion, we indicate the type and magnitude of the bandgap.



the 3s, 3p, and 3d orbitals for sulfur.35,50 All crystal structures have been optimized by using the conjugate gradient method with a force-threshold criterion of 10−4 hartree per Bohr radius. Finally, the thickness of the monolayers are taken to be 5.24 Å for phosphorene, 5.5 Å for arsenene, and 5.6 Å for SnS. To express the thermoelectric parameters using an experimentally accessible quantity, one can easely relate the chemical potential to the carrier concentration by using the following expression:52 n=

∫ De(E) e(E−μ)/1k T + 1 dE B

RESULTS AND DISCUSSION Black phosphorus, black arsenic, and SnS are layered materials with orthorhombic structure in which individual layers are stacked together through van der Waals interactions. Their monolayer form has four atoms per unit cell, as it is depicted in Figure 1. Each atom is pyramidally bonded to three neighboring atoms of the same (phosphorene and arsenene, for homoatomic) or different (tin sulfide (SnS), for heteroatomic) species forming a puckered-like honeycomb lattice. As shown in Table 1, the lattice constants computed with the DFTB approach quantitatively agree (error ≤5%) with those obtained at the full DFT level by other authors for all three materials. Electronic and Phononic Band Structures. The electronic band structure of the puckered materials is displayed in the three central panels of Figure 1. For phosphorene, the valence-band maximum (VBM) occurs along the Γ−Y path while the conduction-band minimum (CBM) is located directly at the Γ point. This results in a bandgap of about 1.94 eV,

(7)

Here, n is the carried density, De is the electronic density of states obtained from the electronic structure calculations, kB is Boltzmann’s constant, and T is the temperature. Additional details of the temperature dependence of the chemical potential for the 2D puckered systems are shown in the Supporting Information. 18843

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C

phosphorene24,29 and arsenene,31,34 but not for SnS monolayers. On the basis of these results, we expect that SnS monolayers will display the largest anisotropy in the thermal conductance due to the predominance of acoustic modes in the thermal transport, and this may have an impact on its thermoelectric properties. Transport Properties. All transport calculations have been carried out using periodic boundary conditions in the perpendicular direction to transport and with the same number of unit cells. Moreover, to avoid additional effects arising from different material compositions of the electronic and phonon reservoirs and of the scattering region, we have chosen the reservoirs to consist of the same material as the scattering region. With this, we can reveal intrinsic transport features of the different systems. The electron and phonon transmission spectra are calculated for each material along the zigzag (ZZ) and armchair (AC) direction at equilibrium conditions. As it can be seen in Figure 2a,b for both phosphorene and arsenene,

Table 1. Calculated Lattice Constants of 2D Puckered Materials along Zigzag (ZZ) and Armchair (AC) Directionsa transport direction systems

ZZ [Å]

AC [Å]

other works (ZZ, AC) [Å]

phosphorene arsenene SnS monolayer

3.49 3.81 3.93

4.34 4.75 4.51

(3.28, 4.43),24 (3.32, 4.58)53 (3.68, 4.77),31 (3.69, 4.77)34 (4.03, 4.26),36 (4.01, 4.35)54

a

For comparison, the lattice constants from other published theoretical studies are given. In general, the difference between the DFTB lattice parameters agrees quite well (error ≤5%) with those at the full DFT level.

which is in agreement with other theoretical studies.19,29 In contrast to what is reported by Fei et al.,29 we observe an indirect bandgap, which is due to the fact that a second valenceband edge at the Γ point is slightly lower than the VBM. Puckered arsenene also shows an indirect bandgap of ∼1.08 eV, in agreement with previous works.31,55 For heteroatomic SnS monolayer, we obtained a direct gap of about 1.96 eV with VBM and CBM located at the Γ point. Despite showing a different type of bandgap to that observed in recent firstprinciple studies,32,35,36,54 our bandgap magnitude is nevertheless close to the reported values at a higher level of theory, e.g., using hybrid functionals for the exchange-correlation potential (see the Supporting Information for a more detailed discussion about the electronic band structure). Considering the underestimation of bandgaps in conventional DFT calculations, we believe that our results are within an adequate range to study the thermoelectric properties of these materials. We also remark at this place that the ballistic transport features of all the studied systems is less sensitive to the type of bandgap, since for each energy an integration over k vectors is carried out. Based on the phonon dispersion analysis (see the three right panels of Figure 1), all studied systems are mechanically stable and do not show imaginary modes. As one can see in the phonon spectrum, the acoustic branches display the typical dispersion of 2D materials: the longitudinal (LA) and transversal (TA) acoustic branches have linear dispersion as the wave vector approaches the Γ point, while out-of-plane ZA branches exhibit quadratic dispersion due to the rapid decay of transversal forces. Moreover, it also can be seen that the distribution and dispersion of the phonon branches for homoatomic puckered materials is almost identical, except for their maximum frequency value, which is a consequence of the mass difference between As (∼75 uma) and P (∼31 uma). In general, the DFTB phonon band structures for P and As agree quite well with those computed with DFT (the results are included in the Supporting Information). Only for SnS the high-frequency optical modes are shifted upward. Comparing the slope of the acoustic branches (i.e., the group velocity) for ZZ (Γ → X) and AC (Γ → Y) transport directions, it is expected that these materials will display strong thermal anisotropy. For instance, the longitudinal acoustic (LA) branch in phosphorene shows a group velocity of 8.35 and 4.74 km/s along the Γ−X (ZZ) and Γ−Y (AC) directions, respectively, which are very close to previous DFT results.25,29,53 Our values for arsenene, ZZ-5.01 km/s and AC-2.71 km/s, are also in agreement with those reported by Zeraati et al.34 SnS monolayer gives group velocities of ZZ-6.48 km/s and AC-2.14 km/s. It is worth mentioning that, up to now, thermal anisotropy has only been reported for

Figure 2. Transport properties. Electron and phonon transmission functions along the ZZ (blue line) and AC (red line with shadow) transport direction: (a, d) phosphorene, (b, e) arsenene, and (c, f) SnS monolayer, respectively.

the electron transmission functions for ZZ and AC directions are slightly different at regions near the band edges but become increasingly different around the band centers. Hence, we find a weak electron transport anisotropy, in agreement with previously reported DFT calculations for phosphorene.56 Interestingly, DFT as well as the DFTB calculations give distinctive different electron effective masses (DFTB: m* = 0.47me along AC direction and 1.43me along ZZ direction, DFT:57 m* = 0.17me along AC direction and 1.12me along ZZ direction) which would suggest within a classical transport model a distinct anisotropy in the electron transport. The situation is different, however, for SnS, where the transmission of the conduction band is 2 times larger for the ZZ direction than for the AC direction, showing significant anisotropic electronic properties (see Figure 2c). In the next step, we computed the phonon transmission function in the ballistic regime (see Figure 2d−f). As it is expected from the phonon dispersion analysis presented before, phosphorene and arsenene have relatively similar phonon 18844

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C

appreciably depend on the AC or ZZ directions (see Figure 3a−c). However, when combining the corresponding electronic conductances and Seebeck coefficients into the power factor (PF = S2G), anisotropic features become more evident, as shown in Figure 3d−f. The anisotropy is clearer seen for positive μ values, with larger PF in the ZZ direction. Notice also that the PF exhibits some sharp peaks at certain μ values, corresponding to the positions of the subband edges. Unlike graphene, the maximum value of PF for puckered materials is not located on the vicinity of the edge of the first conduction band.58 Some higher peaks were also found at doping concentrations where the Seebeck coefficient is small; for phosphorene and arsenene such high values are obtained in the μ > 0 region, while for SnS they lie in the μ < 0 domain. Despite the nonuniform temperature dependence of the Seebeck coefficient with respect to μ, the PF increases with increasing temperature for the three systems. Thus, PF peaks widen and can overlap with nearest-neighbor peaks. It is also seen in Figure 3 that increasing the temperature decreases the anisotropy in the PF, which could be further exploited to tune the thermoelectric performance. In fact, the high Seebeck coefficient combined with the remarkable anisotropic behavior displayed by the PF for arsenene and SnS monolayers, compared to those for phosphorene, could suggest better thermoelectric properties if only the electronic contributions are taken into account. These values are also higher than those theoretically obtained for novel two-dimensional materials such as monolayer transitionmetal dichalcogenides15,17,59−61 and graphyne.14,58 However, to have a more accurate understanding of the thermoelectric response, we need to consider the influence of phononic contributions to the thermal conductance. The contributions of electrons, κel, and phonons, κph, to the thermal conductance are plotted in Figure 4. Here, to have a

transmission functions (Figure 2d,e). The only differences are the phonon bandgaps (∼67 and ∼45 cm−1, respectively) and the maximum frequency value in the phonon dispersion. Likewise, we have found that thermal transport in the ZZ direction is larger than in the AC direction for almost the whole frequency range, in agreement with previous studies.24,29,31,34 In particular, SnS monolayer turns out to display the strongest thermal anisotropy. Anisotropic Thermoelectric Response. To gain additional insights, we plot in Figure 3a−c the Seebeck coefficient S

Figure 3. Thermoelectric properties. Chemical potential dependence of the Seebeck coefficient (S) and power factor (PF): (a, d) phosphorene, (b, e) arsenene, and (c, f) SnS monolayer, respectively. We plot the results for both zigzag and armchair direction at two different temperatures: 300 K (solid lines) and 600 K (dashed lines).

as a function of the chemical potential μ for the three systems at 300 K for both transport directions. We remark that μ is positive for n-type doping and negative for p-type doping. S is symmetric with respect to μ = 0 for homoatomic structures, but in the case of SnS monolayer, it displays an asymmetric behavior for |μ| > 2.0 eV (∼1.36 × 1019 m−2), related to the breaking of electron−hole symmetry in the electronic structure. The SnS monolayer displays the maximum Seebeck coefficient at 300 K, ∼1.45 × 10−3 eV/K, followed by phosphorene (∼1.3 × 10−3 eV/K) and arsenene (∼1.18 × 10−3 eV/K). Smax for SnS monolayer is comparable to that found for its bulk phase41 but larger than that computed for SnSe systems.38,41 For phosphorene, Smax is smaller than the value reported by Fei et al. using BTE.29 Moreover, we have found that S(μ) has a different behavior as a function of temperature for each material. For phosphorene and SnS monolayer, the μ value at which the Seebeck coefficient becomes largest lies rather close to the charge neutrality point (CNP), its magnitude first increases and then decreases, and at 600 K, it has almost the same value to that obtained at room temperature. Whereas, for arsenene, the Seebeck coefficient reaches its maximum at even lower μ values than in the two other systems, and Smax is almost halved at 600 K. It can be also seen that independently of the temperature, the Seebeck coefficient is nearly isotropic; i.e., it does not

Figure 4. Thermoelectric properties. Chemical potential dependence of the electronic (κel) and phononic (κph) contribution to the thermal conductance: (a, d) phosphorene, (b, e) arsenene, and (c, f) SnS monolayer, respectively. We plot the results of κel for both zigzag and armchair direction at two different temperatures: 300 K (solid lines) and 600 K (dashed lines). Here, we can easily see the thermal anisotropy in these 2D puckered materials, κph−ZZ > κph−AC. 18845

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C better comparison with other works, we have considered a length of 100 nm to convert thermal conductance [in units of W/K] into thermal conductivity [in units of W/mK]. For electrons (see Figure 4a−c), κel strongly depends on μ and on the temperature. After increasing the temperature, we can see an increase in the electronic contribution, although the differences between ZZ and AC directions discussed in Figure 3 remain. Moreover, all three systems have similar values of kel. This behavior has also been found in other heteroatomic puckered materials (bulk).41 Regarding the phonon contribution κph (see Figure 4d−f), we have found that for very low temperatures (T < 20 K) the thermal conductance is independent of the transport direction, whereas when the temperature increases, κph increases (because of the excitation of higher modes) and anisotropic effects become more pronounced due to the increasing contribution of acoustic modes to thermal transport. For the three materials it turns out that the thermal anisotropy is stronger than the electronic one, the latter being almost negligible. The preferred direction for phonon thermal conductance is along the ZZ direction, which is in good agreement with previous works.23,29,34,42,53 As it is expected due to the strong anisotropy in its phonon transmission function, SnS monolayer displays the highest thermal anisotropy at 300 K, PAZ = κph−ZZ/κph−AC ∼ 2.5, followed by phosphorene (PAZ ∼ 1.6), which is slightly more thermal anisotropic than arsenene (PAZ ∼ 1.35). The PAZ value for phosphorene is slightly different from that obtained by employing DFT calculations,26,53 while for arsenene there is more discrepancy from that reported by Zeraati et al.,34 which is mainly due to the methodology used. In the latter work they employed self-consistent calculations to solve the Boltzmann transport equation. Unlike the electronic contributions, κph remains constant above 200 K, independently of the transport direction, since the whole phonon spectrum has been already covered by the frequency integration in eq 5. Comparing both contributions, it can be seen that phonons are dominant in the thermal conductance only for μ ranges nearby the subband edges. It is worth mentioning that as previously discussed, we have found sharp peaks in the PF in those μ ranges (see Figure 3). Thus, the influence of anisotropy on the thermoelectric response will be governed by the thermal one. Since SnS monolayers display the most optimal ration of high power factor to small thermal conductance when compared with arsenene and phosphorene, we may expect that its thermoelectric performance may also be the best one. Combining now all the transport quantities, we show the results for the figure of merit, ZT, at a fixed temperature of 300 K for ZZ and AC direction in Figures 5a and 5b, respectively. We can see that in all studied materials ZT displays sharp peaks which are located at μ values inside the region where κel < κph. For phosphorene and SnS, the closest peaks to the CNP (μ = 0.0 eV) are located almost at the same μ values because of the similar bandgap energy (∼1.95 eV, see also Figure 1), whereas for arsenene, the main ZT peaks appear around |μ| ∼ 0.5 eV (∼8.0 × 1016 m−2). ZT peaks for transport in the AC direction are higher than those corresponding to the ZZ direction. This is a consequence of the competition between the weak electrical and strong thermal anisotropy found in all three systems. Another important factor is their extremely low phonon thermal conductance compared to conventional two-dimensional materials, which helps to increase ZT. Thus, SnS monolayer, which has the lowest thermal conductivity and the

Figure 5. Thermoelectric properties. Figure of merit as a function of the carrier concentration at 300 K for (a) zigzag and (b) armchair direction in phosphorene, arsenene, and SnS monolayer. In their inset, we show the main peaks of ZT for SnS monolayer. Temperature dependence of the maximal figure of merit, ZTmax, for the 2D puckered materials when transport is along the (c) zigzag and (d) armchair direction.

highest power factor, displays the best thermoelectric performance in both transport directions, reaching ZT = 0.95 and ZT = 1.6 at 300 K for ZZ and AC directions, respectively. These values are larger than those found for another heteroatomic puckered material, SnSe monolayer,38 and its bulk phase.37,42 We have also found that arsenene has better thermoelectric performance than phosphorene, showing a maximum figure of merit ZTmax ∼ 0.85 at 300 K. The ZT for phosphorene is close to the reported by Fei et al. employing the BTE29 and higher than bulk black phosphorus.62 Unlike the homoatomic materials, the main sharp peaks of the SnS monolayer are located at μ ∼ −2.1 and −3.4 eV (1.445 × 1019 and 1.452 × 1019 m−2, respectively; p-type doping), which are related to the subband edges of the electron dispersion (see insets in Figures 5a,b). In Figures 5c and 5d we then show the temperature dependence of ZTmax in the ZZ and AC directions. SnS monolayer has the highest ZT values for the whole temperature range, and at 800 K values as high as ZTZZ = 2.9 and ZTAC = 4.5 are obtained. These values are comparable to those found for tailored graphene nanoribbons.6,63,64 It is worth nothing that, despite having quite similar electronic transport properties, arsenene shows larger ZT values than phosphorene due to its lower thermal conductance. Moreover, ZTmax values are μ dependent, and they display a different behavior for each material after increasing the temperature. For instance, phosphorene shows a better thermoelectric performance for p-type doping, and μ corresponding to ZTmax changes from −0.98 to −0.93 eV for transport in the ZZ direction and from −0.98 to −0.90 eV for the AC direction, whereas for arsenene, n-type doping turns out to be more efficient in maximizing ZT, and the chemical potential varies in the range of 0.72−0.50 and 0.55−0.43 eV for ZZ and AC directions, respectively. In the case of SnS monolayer, finally, p-type doping is more efficient to optimize ZT; the chemical potential changes from −2.12 to −2.22 eV for the ZZ direction and from −3.42 to −3.28 eV for the AC direction. These effects are closely related to the temperature dependence of the chemical potential at a given 18846

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C

effect decreases when the temperature is raised. This happens mainly due to the temperature dependence of the electron transport properties and quantities related because the phonon thermal conductance remains almost constant for T > 200 K as shown in Figure 4. We have also observed that by increasing the temperature, besides improving the thermoelectric performance, the ZT peaks widen and, hence, reduce the gap around μ = 0.0 eV. This effect has also been reported for other nanomaterials.6,29,44,58

carrier concentration (see Supporting Information for details). Hence, to get the large thermoelectric response displayed by the SnS monolayer, we need a strong p-type carrier concentration. To further quantify the anisotropy of the thermoelectric response, we have defined the quantity MAZ = ZTAC/ZTZZ and plot it in Figure 6. The MAZ(μ,T) features are not symmetric



CONCLUSIONS In conclusion, we have performed a detailed theoretical study of the thermoelectric response of phosphorene, arsenene, and SnS monolayers. As we have shown, these systems show no strong electronic anisotropy, but a strong thermal one, which consequently leads to a strong anisotropic thermoelectric response. SnS monolayers show the highest thermal anisotropy (κph−ZZ > κph−AC) as well as the largest figure of merit ZT in both transport directions, attaining values of up to ZTZZ = 2.9 and ZTAC = 4.5 at 800 K, which can be reached by tuning the chemical potential to have a strong p-type carrier system. ZTmax(T) values for SnS monolayers are also larger than those reported for SnSe monolayers, which is accepted to be a material with a high potential for thermoelectric applications. We have also found that arsenene may be a good candidate for designing thermoelectric devices because of its high ZT ∼ 1.0 along the AC direction at moderate temperatures. Puckered arsenene displays the strongest anisotropic thermoelectric response, despite the fact that it has the weakest thermal anisotropy. Hence, these new puckered materials can be considered to offer a new playground to fine-tunee.g., through directional strain engineering or by producing nanoribbonsthe thermoelectric response.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.6b04969. Detailed description of DFTB method and comparison with DFT calculations (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail [email protected]; Ph +49 (0)351 4633 1419; Fax +49 (0)351 4633 1422 (R.G.). Notes

The authors declare no competing financial interest.



Figure 6. Quantitative analysis of the anisotropic thermoelectric performance of 2D puckered materials. We show the patterns of the temperature and chemical potential dependences of MAZ = ZTAC/ZTZZ for (a) phosphorene, (b) arsenene, and (c) SnS monolayer. The region in black means that ZTZZ < 0.01.

ACKNOWLEDGMENTS L.M.S. thanks the International Max Planck Research School Dynamical Processes in Atoms, Molecules and Solids and the Deutscher Akademischer Austauschdienst (DAAD) for the financial support. D.T. acknowledges financial support from the Junior Research Group Computer-Simulationen für das Materialdesign (CoSiMa). This work has also been partly supported by the German Research Foundation (DFG) within the Cluster of Excellence “Center for Advancing Electronics Dresden”. We acknowledge the Center for Information Services and High Performance Computing (ZIH) at TU Dresden for computational resources.

with respect to the charge neutrality point, which indicates that the anisotropic features are dependent on the doping type. Thus, the maximum MAZ is found for phosphorene and arsenene at n-type doping, while for SnS monolayer at p-type doping. Despite having the weaker thermal anisotropy, arsenene displays the maximal anisotropic response, MAZ = 4.0 at 200 K (Figure 6b), followed by SnS monolayer with MAZ = 3.5 (Figure 6c). Moreover, we have found this anisotropic 18847

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C



Conductivity Observed in Few-Layer Black Phosphorus. Nat. Commun. 2015, 6, 8572. (24) Jain, A.; McGaughey, A. J. H. Strongly Anisotropic In-Plane Thermal Transport in Single-Layer Black Phosphorene. Sci. Rep. 2015, 5, 8501. (25) Liu, T.-H.; Chang, C.-C. Anisotropic Thermal Transport in Phosphorene: Effects of Crystal Orientation. Nanoscale 2015, 7, 10648−10654. (26) Ong, Z.-Y.; Cai, Y.; Zhang, G.; Zhang, Y.-W. Strong Thermal Transport Anisotropy and Strain Modulation in Single-Layer Phosphorene. J. Phys. Chem. C 2014, 118, 25272−25277. (27) Xia, F.; Wang, H.; Jia, Y. Rediscovering Black Phosphorus as an Anisotropic Layered Material for Optoelectronics and Electronics. Nat. Commun. 2014, 5, 4458. (28) Sofer, Z.; Sedmidubský, D.; Huber, S.; Luxa, J.; Bousa, D.; Boothroyd, C.; Pumera, M. Layered Black Phosphorus: Strongly Anisotropic Magnetic, Electronic, and Electron-Transfer Properties. Angew. Chem., Int. Ed. 2016, 55, 3382−3386. (29) Fei, R.; Faghaninia, A.; Soklaski, R.; Yan, J.-A.; Lo, C.; Yang, L. Enhanced Thermoelectric Efficiency via Orthogonal Electrical and Thermal Conductances in Phosphorene. Nano Lett. 2014, 14, 6393− 6399. (30) Zhang, S.; Yan, Z.; Li, Y.; Chen, Z.; Zeng, H. Atomically Thin Arsenene and Antimonene: Semimetal-Semiconductor and IndirectDirect Band-Gap Transitions. Angew. Chem., Int. Ed. 2015, 54, 3112− 3115. (31) Kamal, C.; Ezawa, M. Arsenene: Two-dimensional Buckled and Puckered Honeycomb Arsenic Systems. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 085423. (32) Gomes, L. C.; Carvalho, A. Phosphorene analogues: Isoelectronic Two-Dimensional Group-IV Monochalcogenides with Orthorhombic Structure. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 085406. (33) Zhu, Z.; Guan, J.; Liu, D.; Tománek, D. Designing Isoelectronic Counterparts to Layered Group V Semiconductors. ACS Nano 2015, 9, 8284−8290. (34) Zeraati, M.; Vaez Allaei, S. M.; Abdolhosseini Sarsari, I.; Pourfath, M.; Donadio, D. Highly Anisotropic Thermal Conductivity of Arsenene: An ab-initio Study. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 085424. (35) Brent, J. R.; Lewis, D. J.; Lorenz, T.; Lewis, E. A.; Savjani, N.; Haigh, S. J.; Seifert, G.; Derby, B.; O’Brien, P. Tin(II) Sulfide (SnS) Nanosheets by Liquid-Phase Exfoliation of Herzenbergite: IV-VI Main Group Two-Dimensional Atomic Crystals. J. Am. Chem. Soc. 2015, 137, 12689−12696. (36) Fei, R.; Li, W.; Li, J.; Yang, L. Giant Piezoelectricity of Monolayer Group IV Monochalcogenides: SnSe, SnS, GeSe, and GeS. Appl. Phys. Lett. 2015, 107, 173104. (37) Hao, S.; Shi, F.; Dravid, V. P.; Kanatzidis, M. G.; Wolverton, C. Computational Prediction of High Thermoelectric Performance in Hole Doped Layered GeSe. Chem. Mater. 2016, 28, 3218−3226. (38) Wang, F. Q.; Zhang, S.; Yu, J.; Wang, Q. Thermoelectric Properties of Single-Layered SnSe Sheet. Nanoscale 2015, 7, 15962− 15970. (39) Xia, J.; Li, X.-Z.; Huang, X.; Mao, N.; Zhu, D.-D.; Wang, L.; Xu, H.; Meng, X.-M. Physical Vapor Deposition Synthesis of TwoDimensional Orthorhombic SnS Flakes with Strong Angle/Temperature-Dependent Raman Responses. Nanoscale 2016, 8, 2063−2070. (40) Patel, M.; Chavda, A.; Mukhopadhyay, I.; Kim, J.; Ray, A. Nanostructured SnS with Inherent Anisotropic Optical Properties for High Photoactivity. Nanoscale 2016, 8, 2293−2303. (41) Ding, G.; Gao, G.; Yao, K. High-Efficient Thermoelectric Materials: The Case of Orthorhombic IV-VI Compounds. Sci. Rep. 2015, 5, 9567. (42) Guo, R.; Wang, X.; Kuang, Y.; Huang, B. First-Principles Study of Anisotropic Thermoelectric Transport Properties of IV-VI Semiconductor Compounds SnSe and SnS. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 115202.

REFERENCES

(1) Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Nat. Mater. 2008, 7, 105−114. (2) Xu, Y.; Li, Z.; Duan, W. Thermal and Thermoelectric Properties of Graphene. Small 2014, 10, 2182−2199. (3) Heremans, J. P.; Dresselhaus, M. S.; Bell, L. E.; Morelli, D. T. When Thermoelectrics Reached the Nanoscale. Nat. Nanotechnol. 2013, 8, 471−473. (4) Dresselhaus, M.; Chen, G.; Tang, M.; Yang, R.; Lee, H.; Wang, D.; Ren, Z.; Fleurial, J.-P.; Gogna, P. New Directions for LowDimensional Thermoelectric Materials. Adv. Mater. 2007, 19, 1043− 1053. (5) Nozaki, D.; Avdoshenko, S. M.; Sevincli, H.; Cuniberti, G. Quantum Interference in Thermoelectric Molecular Junctions: A Toy Model Perspective. J. Appl. Phys. 2014, 116, 074308. (6) Sevincli, H.; Sevik, C.; Cagin, T.; Cuniberti, G. A Bottom-up Route to Enhance Thermoelectric Figures of Merit in Graphene Nanoribbons. Sci. Rep. 2013, 3, 1228. (7) Haskins, J.; Kinaci, A.; Sevik, C.; Sevincli, H.; Cuniberti, G.; Cagin, T. Control of Thermal and Electronic Transport in DefectEngineered Graphene Nanoribbons. ACS Nano 2011, 5, 3779−3787. (8) Carrete, J.; Mingo, N.; Wang, S.; Curtarolo, S. Nanograined HalfHeusler Semiconductors as Advanced Thermoelectrics: An ab-initio High-Throughput Statistical Study. Adv. Funct. Mater. 2014, 24, 7427− 7432. (9) Mingo, N.; Hauser, D.; Kobayashi, N. P.; Plissonnier, M.; Shakouri, A. Nanoparticle-in-Alloy Approach to Efficient Thermoelectrics: Silicides in SiGe. Nano Lett. 2009, 9, 711−715. (10) Kanatzidis, M. G. Nanostructured Thermoelectrics: The New Paradigm? Chem. Mater. 2010, 22, 648−659. (11) Miro, P.; Audiffred, M.; Heine, T. An Atlas of Two-dimensional Materials. Chem. Soc. Rev. 2014, 43, 6537−6554. (12) Chang, P.-H.; Bahramy, M. S.; Nagaosa, N.; Nikolic, B. K. Giant Thermoelectric Effect in Graphene-Based Topological Insulators with Heavy Adatoms and Nanopores. Nano Lett. 2014, 14, 3779−3784. (13) Zhang, G.; Zhang, Y.-W. Strain Effects on Thermoelectric Properties of Two-dimensional Materials. Mech. Mater. 2015, 91 (Part 2), 382−398. (14) Tan, X.; Shao, H.; Hu, T.; Liu, G.; Jiang, J.; Jiang, H. High Thermoelectric Performance in Two-dimensional Graphyne Sheets Predicted by First Principles Calculations. Phys. Chem. Chem. Phys. 2015, 17, 22872−22881. (15) Lee, C.; Hong, J.; Whangbo, M.-H.; Shim, J. H. Enhancing the Thermoelectric Properties of Layered Transition-Metal Dichalcogenides 2H-MQ2 (M = Mo, W; Q = S, Se, Te) by Layer Mixing: Density Functional Investigation. Chem. Mater. 2013, 25, 3745−3752. (16) Kannan, P. K.; Late, D. J.; Morgan, H.; Rout, C. S. Recent Developments in 2D Layered Inorganic Nanomaterials for Sensing. Nanoscale 2015, 7, 13293−13312. (17) Kumar, S.; Schwingenschloegl, U. Thermoelectric Response of Bulk and Monolayer MoSe2 and WSe2. Chem. Mater. 2015, 27, 1278− 1284. (18) Butler, S. Z.; et al. Progress, Challenges, and Opportunities in Two-Dimensional Materials Beyond Graphene. ACS Nano 2013, 7, 2898−2926. (19) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Black Phosphorus Field-Effect Transistors. Nat. Nanotechnol. 2014, 9, 372−377. (20) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tománek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033−4041. (21) Du, H.; Lin, X.; Xu, Z.; Chu, D. Recent Developments in Black Phosphorus Transistors. J. Mater. Chem. C 2015, 3, 8760−8775. (22) Pan, Y.; Wang, Y.; Ye, M.; Quhe, R.; Zhong, H.; Song, Z.; Peng, X.; Yu, D.; Yang, J.; Shi, J.; Lu, J. Monolayer Phosphorene-Metal Contacts. Chem. Mater. 2016, 28, 2100−2109. (23) Luo, Z.; Maassen, J.; Deng, Y.; Du, Y.; Garrelts, R. P.; Lundstrom, M. S.; Ye, P. D.; Xu, X. Anisotropic In-Plane Thermal 18848

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849

Article

The Journal of Physical Chemistry C (43) Zhao, L.-D.; Lo, S.-H.; Zhang, Y.; Sun, H.; Tan, G.; Uher, C.; Wolverton, C.; Dravid, V. P.; Kanatzidis, M. G. Ultralow Thermal Conductivity and High Thermoelectric Figure of Merit in SnSe Crystals. Nature 2015, 508, 373−377. (44) Medrano-Sandonas, L.; Gutierrez, R.; Pecchia, A.; Dianat, A.; Cuniberti, G. Thermoelectric Properties of Functionalized Graphene Grain Boundaries. J. Self-Assem. Mol. Electron. 2015, 2015, 1−20. (45) Gandi, A. N.; Schwingenschloegl, U. WS2 as an Excellent HighTemperature Thermoelectric Material. Chem. Mater. 2014, 26, 6628− 6637. (46) Pecchia, A.; Penazzi, G.; Salvucci, L.; Di Carlo, A. Nonequilibrium Green’s Functions in Density Functional Tight Binding: Method and Applications. New J. Phys. 2008, 10, 065022. (47) Zhang, W.; Fisher, T. S.; Mingo, N. The Atomistic Green’s Function Method: An Efficient Simulation Approach for Nanoscale Phonon Transport. Numer. Heat Transfer, Part B 2007, 51, 333−349. (48) Frauenheim, T.; Seifert, G.; Elsterner, M.; Hajnal, Z.; Jungnickel, G.; Porezag, D.; Suhai, S.; Scholz, R. A Self-Consistent Charge Density-Functional Based Tight-Binding Method for Predictive Materials Simulations in Physics, Chemistry and Biology. Phys. Status Solidi B 2000, 217, 41−62. (49) Seifert, G. Tight-Binding Density Functional Theory: An Approximate Kohn-Sham DFT Scheme. J. Phys. Chem. A 2007, 111, 5609−5613. (50) Lorenz, T.; Joswig, J.-O.; Seifert, G. Combined SnS@SnS2 Double Layers: Charge Transfer and Electronic Structure. Semicond. Sci. Technol. 2014, 29, 064006. (51) Seifert, G.; Hernández, E. Theoretical Prediction of Phosphorus Nanotubes. Chem. Phys. Lett. 2000, 318, 355−360. (52) Ohkubo, I.; Mori, T. Anisotropic Anomalies of Thermoelectric Transport Properties and Electronic Structures in Layered Complex Nitrides AMN2 (A = Na, Cu; M = Ta, Nb). Chem. Mater. 2015, 27, 7265−7275. (53) Qin, G.; Yan, Q.-B.; Qin, Z.; Yue, S.-Y.; Hu, M.; Su, G. Anisotropic Intrinsic Lattice Thermal Conductivity of Phosphorene from First Principles. Phys. Chem. Chem. Phys. 2015, 17, 4854−4858. (54) Kamal, C.; Chakrabarti, A.; Ezawa, M. Direct Band Gaps in Group IV-VI Monolayer Materials: Binary Counterparts of Phosphorene. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 93, 125428. (55) Zhang, Z.; Xie, J.; Yang, D.; Wang, Y.; Si, M.; Xue, D. Manifestation of Unexpected Semiconducting Properties in Few-layer Orthorhombic Arsenene. Appl. Phys. Express 2015, 8, 055201. (56) Ç akır, D.; Sahin, H.; Peeters, F. M. Tuning of the Electronic and Optical Properties of Single-layer Black Phosphorus by Strain. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 205421. (57) Qiao, J.; Kong, X.; Hu, Z.-X.; Yang, F.; Ji, W. High-Mobility Transport Anisotropy and Linear Dichroism in Few-Layer Black Phosphorus. Nat. Commun. 2014, 5, 4475. (58) Wang, X.-M.; Mo, D.-C.; Lu, S.-S. On the Thermoelectric Transport Properties of Graphyne by the First-principles Method. J. Chem. Phys. 2013, 138, 204704. (59) Jin, Z.; Liao, Q.; Fang, H.; Liu, Z.; Liu, W.; Ding, Z.; Luo, T.; Yang, N. A Revisit to High Thermoelectric Performance of SingleLayer MoS2. Sci. Rep. 2015, 5, 18342. (60) Huang, W.; Da, H.; Liang, G. Thermoelectric Performance of MX2 (M = Mo,W; X = S,Se) Monolayers. J. Appl. Phys. 2013, 113, 104304. (61) Huang, W.; Luo, X.; Gan, C. K.; Quek, S. Y.; Liang, G. Theoretical Study of Thermoelectric Properties of Few-Layer MoS2 and WSe2. Phys. Chem. Chem. Phys. 2014, 16, 10866−10874. (62) Qin, G.; Yan, Q.-B.; Qin, Z.; Yue, S.-Y.; Cui, H.-J.; Zheng, Q.-R.; Su, G. Hinge-Like Structure Induced Unusual Properties of Black Phosphorus and New Strategies to Improve the Thermoelectric Performance. Sci. Rep. 2014, 4, 6946. (63) Mazzamuto, F.; Hung Nguyen, V.; Apertet, Y.; Caër, C.; Chassat, C.; Saint-Martin, J.; Dollfus, P. Enhanced Thermoelectric Properties in Graphene Nanoribbons by Resonant Tunneling of Electrons. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 83, 235426.

(64) Saha, K. K.; Markussen, T.; Thygesen, K. S.; Nikolić, B. K. Multiterminal Single-Molecule Graphene-Nanoribbon Junctions with the Thermoelectric Figure of Merit Optimized Via Evanescent Mode Transport and Gate Voltage. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 041412.

18849

DOI: 10.1021/acs.jpcc.6b04969 J. Phys. Chem. C 2016, 120, 18841−18849