Elastic Bending Modulus for Single-Layer Black Phosphorus

Hao-Yu Zhang,Jin-Wu Jiang
DOI: https://doi.org/10.1088/0022-3727/48/45/455305
2015-04-17
Abstract:We derive an analytic formula for the elastic bending modulus of single-layer black phosphorus (SLBP) based on the valence force field model. The obtained elastic bending modulus is 4.8028 eV and 7.9905 eV along the armchair and zigzag directions in the SLBP, respectively. These values are obviously larger than the bending modulus of 1.4 eV in graphene due to the intrinsic finite thickness for SLBP. Our derivation analytically illustrates that the elastic bending modulus of the SLBP is proportional to the square of the intrinsic thickness of the SLBP.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to calculate the elastic bending modulus of single - layer black phosphorus (SLBP). Specifically, the author derives the elastic bending modulus of SLBP in the armchair and zigzag directions by analytical methods and analyzes its anisotropic characteristics. ### Main problems: 1. **Lack of research on the elastic bending modulus of SLBP**: Although the elastic bending moduli of two - dimensional materials such as graphene and molybdenum disulfide have been studied, the elastic bending modulus of single - layer black phosphorus has not been predicted or calculated yet. 2. **Understanding the unique properties of SLBP**: Due to the unique wrinkled structure of SLBP, it exhibits significant anisotropy in different directions. Therefore, understanding its bending modulus is crucial for further exploring its applications. ### Solutions: - **Using the valence force field model (VFFM)**: The author describes the inter - atomic interactions based on the VFFM model and derives the elastic bending modulus of SLBP. - **Considering the influence of finite thickness**: Analytical derivations show that the elastic bending modulus of SLBP is proportional to the square of its finite thickness, which explains why the bending modulus of SLBP is significantly larger than that of single - layer graphene (1.4 eV), being 4.8028 eV (armchair direction) and 7.9905 eV (zigzag direction) respectively. - **Verifying the analytical results**: By comparing with the numerical simulation results, the correctness of the analytical expressions is verified. ### Formula summary: - The definition of the elastic bending modulus \(D\) is: \[ D=\frac{\partial^{2}W}{\partial\kappa^{2}} \] where \(W\) is the bending energy density and \(\kappa\) is the curvature. - Definition of geometric coefficients: \[ \alpha_{i}^{\sigma}=\left.\frac{\partial d_{i}^{\sigma}}{\partial\kappa}\right|_{\kappa = 0},\quad\beta_{i}^{\sigma}=\left.\frac{\partial\theta_{i}^{\sigma}}{\partial\kappa}\right|_{\kappa = 0} \] These coefficients are used to describe the changes in bond length and angle caused by bending. Through these formulas and models, the author successfully analyzes the elastic bending modulus of SLBP and reveals the physical mechanism of its anisotropy.