Phosphorene, the single layer counterpart of black phosphorus, is a novel two-dimensional semiconductor with high carrier mobility and a large fundamental direct band gap, which has attracted tremendous interest recently. Its potential applications in nano-electronics and thermoelectrics call for fundamental study of the phonon transport. Here, we calculate the intrinsic lattice thermal conductivity of phosphorene by solving the phonon Boltzmann transport equation (BTE) based on first-principles calculations. The thermal conductivity of phosphorene at 300 K is 30.15 W m(-1) K(-1) (zigzag) and 13.65 W m(-1) K(-1) (armchair), showing an obvious anisotropy along different directions. The calculated thermal conductivity fits perfectly to the inverse relationship with temperature when the temperature is higher than Debye temperature (ΘD = 278.66 K). In comparison to graphene, the minor contribution around 5% of the ZA mode is responsible for the low thermal conductivity of phosphorene. In addition, the representative mean free path (MFP), a critical size for phonon transport, is also obtained.

Silicene, the silicon-based counterpart of graphene with a two-dimensional honeycomb lattice, has attracted tremendous interest both theoretically and experimentally due to its significant potential industrial applications. From the aspect of theoretical study, the widely used classical molecular dynamics simulation is an appropriate way to investigate the transport phenomena and mechanisms in nanostructures such as silicene. Unfortunately, no available interatomic potential can precisely characterize the unique features of silicene. Here, we optimized the Stillinger-Weber potential parameters specifically for a single-layer Si sheet, which can accurately reproduce the low buckling structure of silicene and the full phonon dispersion curves obtained from ab initio calculations. By performing equilibrium and nonequilibrium molecular dynamics simulations and anharmonic lattice dynamics calculations with the new potential, we reveal that the three methods consistently yield an extremely low thermal conductivity of silicene and a short phonon mean-free path, suggesting silicene as a potential candidate for high-efficiency thermoelectric materials. Moreover, by qualifying the relative contributions of lattice vibrations in different directions, we found that the longitudinal phonon modes dominate the thermal transport in silicene, which is fundamentally different from graphene, despite the similarity of their two-dimensional honeycomb lattices.

We systematically investigated the geometric, electronic and thermoelectric (TE) properties of bulk black phosphorus (BP) under strain. The hinge-like structure of BP brings unusual mechanical responses such as anisotropic Young's modulus and negative Poisson's ratio. A sensitive electronic structure of BP makes it transform among metal, direct and indirect semiconductors under strain. The maximal figure of merit ZT of BP is found to be 0.72 at 800 K that could be enhanced to 0.87 by exerting an appropriate strain, revealing BP could be a potential medium-high temperature TE material. Such strain-induced enhancements of TE performance are often observed to occur at the boundary of the direct-indirect band gap transition, which can be attributed to the increase of degeneracy of energy valleys at the transition point. By comparing the structure of BP with SnSe, a family of potential TE materials with hinge-like structure are suggested. This study not only exposes various novel properties of BP under strain, but also proposes effective strategies to seek for better TE materials.

New classes two-dimensional (2D) materials beyond graphene, including layered and non-layered, and their heterostructures, are currently attracting increasing interest due to their promising applications in nanoelectronics, optoelectronics and clean energy, where thermal transport property is one of the fundamental physical parameters. In this paper, we systematically investigated the phonon transport properties of 2D orthorhombic group IV-VI compounds of GeS, GeSe, SnS and SnSe by solving the Boltzmann transport equation (BTE) based on first-principles calculations. Despite the similar puckered (hinge-like) structure along the armchair direction as phosphorene, the four monolayer compounds possess diverse anisotropic properties in many aspects, such as phonon group velocity, Young's modulus and lattice thermal conductivity (κ), etc. Especially, the κ along the zigzag and armchair directions of monolayer GeS shows the strongest anisotropy while monolayer SnS and SnSe shows an almost isotropy in phonon transport. The origin of the diverse anisotropy is fully studied and the underlying mechanism is discussed in detail. With limited size, the κ could be effectively lowered, and the anisotropy could be effectively modulated by nanostructuring, which would extend the applications in nanoscale thermoelectrics and thermal management. Our study offers fundamental understanding of the anisotropic phonon transport properties of 2D materials, and would be of significance for further study, modulation and applications in emerging technologies.

Strain engineering is one of the most promising and effective routes toward continuously tuning the electronic and optic properties of materials, while thermal properties are generally believed to be insensitive to mechanical strain. In this paper, the strain-dependent thermal conductivity of monolayer silicene under uniform bi-axial tension is computed by solving the phonon Boltzmann transport equation with force constants extracted from first-principles calculations. Unlike the commonly believed understanding that thermal conductivity only slightly decreases with increased tensile strain for bulk materials, it is found that the thermal conductivity of silicene first increases dramatically with strain and then slightly decreases when the applied strain increases further. At a tensile strain of 4%, the highest thermal conductivity is found to be about 7.5 times that of unstrained one. Such an unusual strain dependence is mainly attributed to the dramatic enhancement in the acoustic phonon lifetime. Such enhancement plausibly originates from the flattening of the buckling of the silicene structure upon stretching, which is unique for silicene as compared with other common two-dimensional materials. Our findings offer perspectives of modulating the thermal properties of low-dimensional structures for applications such as thermoelectrics, thermal circuits, and nanoelectronics.

Two-dimensional (2D) carbon allotrope called penta-graphene was recently proposed from first-principles calculations and various similar penta-structures emerged. Despite significant effort having been dedicated to electronic structures and mechanical properties, little research has been focused on thermal transport in penta-structures. Motivated by this, we performed a comparative study of thermal transport properties of three representative pentagonal structures, namely penta-graphene, penta-SiC2, and penta-SiN2, by solving the phonon Boltzmann transport equation with interatomic force constants extracted from first-principles calculations. Unexpectedly, the thermal conductivity of the three penta-structures exhibits diverse strain dependence, despite their very similar geometry structures. While the thermal conductivity of penta-graphene exhibits standard monotonic reduction by stretching, penta-SiC2 possesses an unusual nonmonotonic up-and-down behavior. More interestingly, the thermal conductivity of penta-SiN2 has 1 order of magnitude enhancement due to the strain induced buckled to planar structure transition. The mechanism governing the diverse strain dependence is identified as the competition between the change of phonon group velocity and phonon lifetime of acoustic phonon modes with combined effect from the unique structure transition for penta-SiN2. The disparate thermal transport behavior is further correlated to the fundamentally different bonding nature in the atomic structures with solid evidence from the distribution of deformation charge density and more in-depth molecular orbital analysis. The reported giant and robust tunability of thermal conductivity may inspire intensive research on other derivatives of penta-structures as potential materials for emerging nanoelectronic devices. The fundamental physics understood from this study also solidifies the strategy to engineer thermal transport properties of broad 2D materials by simple mechanical strain.

No abstract

Two-dimensional (2D) materials with graphene as a representative have been intensively studied for a long time. Recently, monolayer gallium nitride (ML GaN) with honeycomb structure was successfully fabricated in experiments, generating enormous research interest for its promising applications in nano- and opto-electronics. Considering all these applications are inevitably involved with thermal transport, systematic investigation of the phonon transport properties of 2D GaN is in demand. In this paper, by solving the Boltzmann transport equation (BTE) based on first-principles calculations, we performed a comprehensive study of the phonon transport properties of ML GaN, with detailed comparison to bulk GaN, 2D graphene, silicene and ML BN with similar honeycomb structure. Considering the similar planar structure of ML GaN to graphene, it is quite intriguing to find that the thermal conductivity (κ) of ML GaN (14.93 W mK) is more than two orders of magnitude lower than that of graphene and is even lower than that of silicene with a buckled structure. Systematic analysis is performed based on the study of the contribution from phonon branches, comparison among the mode level phonon group velocity and lifetime, the detailed process and channels of phonon-phonon scattering, and phonon anharmonicity with potential energy well. We found that, different from graphene and ML BN, the phonon-phonon scattering selection rule in 2D GaN is slightly broken by the lowered symmetry due to the large difference in the atomic radius and mass between Ga and N atoms. Further deep insight is gained from the electronic structure. Resulting from the special sp orbital hybridization mediated by the Ga-d orbital in ML GaN, the strongly polarized Ga-N bond, localized charge density, and its inhomogeneous distribution induce large phonon anharmonicity and lead to the intrinsic low κ of ML GaN. The orbitally driven low κ of ML GaN unraveled in this work would make 2D GaN prospective for applications in energy conversion such as thermoelectrics. Our study offers fundamental understanding of phonon transport in ML GaN within the framework of BTE and further electronic structure, which will enrich the studies of nanoscale phonon transport in 2D materials and shed light on further studies.