Variable thermal transport in black, blue, and violet phosphorene from extensive atomistic simulations with a neuroevolution potential

Abstract

Phosphorus has diverse chemical bonds and even in its two-dimensional form there are three stable allotropes: black phosphorene (Black-P), blue phosphorene (Blue-P), and violet phosphorene (Violet-P). Due to the complexity of these structures, no efficient and accurate classical interatomic potential has been developed for them. In this paper, we develop an efficient machine-learned neuroevolution potential model for these allotropes and apply it to study thermal transport in them via extensive molecular dynamics (MD) simulations. Based on the homogeneous nonequilibrium MD method, the thermal conductivities are predicted to be $12.5\pm 0.2$ (Black-P in armchair direction), $78.4\pm 0.4$ (Black-P in zigzag direction), $128\pm 3$ (Blue-P), and $2.36\pm 0.05$ (Violet-P) $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. The underlying reasons for the significantly different thermal conductivity values in these allotropes are unraveled through spectral decomposition, phonon eigenmodes, and phonon participation ratio. Under external tensile strain, the thermal conductivity in black-P and violet-P are finite, while that in blue-P appears unbounded due to the linearization of the flexural phonon dispersion that increases the phonon mean free paths in the zero-frequency limit.

10 pages, 10 figures, code and data available

References in corpus (7)

• Semiconducting layered blue phosphorus: A computational study
• Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• First-Principles Prediction of Phononic Thermal Conductivity of Silicene: a Comparison with Graphene
• Co-existence of size-dependent and size-independent thermal conductivities in single layer black phosphorus
• Improving the accuracy of the neuroevolution machine learning potential for multi-component systems
• Spectral Decomposition of Thermal Conductivity: Comparing Velocity Decomposition Methods in Homogeneous Molecular Dynamics Simulations

Cited by in corpus (5)

• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• Accurate prediction of heat conductivity of water by a neuroevolution potential
• Mechanisms of temperature-dependent thermal transport in amorphous silica from machine-learning molecular dynamics
• Dissimilar thermal transport properties in $\kappa$-Ga${}_2$O${}_3$ and $\beta$-Ga${}_2$O${}_3$ revealed by machine-learning homogeneous nonequilibrium molecular dynamics simulations
• Thermal conductivities of monolayer graphene oxide from machine learning molecular dynamics simulations