Abstract:Thermal conductivity of single layer black phosphorus (BP) is investigated by combining density functional calculations and Peierls-Boltzmann transport equation. Differing from isotropic and divergent thermal conductivities in two-dimensional graphene and MoS$_2$, an compelling co-existence of size-dependent and size-independent thermal conductivities are discovered for single layer BP along zigzag (ZZ) and armchair (AM) direction, respectively. Besides, thermal conductivities in single layer BP are found to be highly anisotropic because of orientation dependent group velocities, e.g., thermal conductivities at 300 K are 83.5 and 24.3 W/m-K along ZZ and AM directions for single layer BP with a size of 10 $\mu m$, respectively.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **Thermal conductivity characteristics of monolayer black phosphorus (BP) and its relationship with size**.
Specifically, by combining density functional theory (DFT) calculations and the Peierls - Boltzmann transport equation (PBTE), the researchers explored the thermal conductivity behavior of monolayer BP in different directions. Different from the isotropic and divergent thermal conductivity observed in two - dimensional materials such as graphene and molybdenum disulfide (MoS₂), monolayer BP exhibits a unique phenomenon: its thermal conductivity **depends on size in the zigzag (ZZ) direction**, while it **does not depend on size in the armchair (AM) direction**.
In addition, the study also found that the thermal conductivity of monolayer BP has significant anisotropy, that is, the thermal conductivity along the ZZ direction is much higher than that along the AM direction. For example, at 300 K, the thermal conductivity of a 10 - micron - sized monolayer BP along the ZZ direction is 83.5 W/m - K, while it is only 24.3 W/m - K along the AM direction.
These findings are crucial for understanding the thermal transport properties of monolayer BP and its potential in nanoelectronics and thermoelectric applications. Specifically, this size - dependence and anisotropy provide a theoretical basis for designing efficient thermal management schemes and contribute to the development of high - performance thermoelectric materials.
### Key Formulas
1. **Calculation formula for thermal conductivity**:
\[
\kappa_\alpha=\frac{1}{N_q V}\sum_{q,j}C_{q,j}v^2_{q,j,\alpha}\tau_{q,j}
\]
where \( N_q \) is the total number of q - points sampled in the Brillouin zone, \( V \) is the volume of the monolayer BP unit cell, \( C_{q,j} \) is the specific heat capacity, \( v_{q,j,\alpha} \) is the group velocity along the transport direction \( \alpha \), and \( \tau_{q,j} \) is the relaxation time of the phonon mode.
2. **Expression for ZA - mode frequency**:
\[
\omega^2(q)=A|q|^4 + B|q|^2
\]
where \( A \) is determined by the bending stiffness \( t \) and the two - dimensional mass density \( \rho_{2D} \), that is, \( A=\frac{t}{\rho_{2D}} \).
3. **Integral form of thermal conductivity with respect to wave vector**:
\[
\kappa\sim\sum_j\int q\int\theta C(v\cos\theta)^2\tau q dq d\theta\sim\sum_j\int q C v^2\tau q dq
\]
These formulas and results reveal the unique behavior of thermal conductivity in monolayer BP and provide theoretical support for further research and applications.