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.

arXiv comment10 pages, 10 figures, code and data available

References in corpus (21)

• Gaussian Approximation Potentials: the accuracy of quantum mechanics, without the electrons
• Tunable Band Gap and Anisotropic Optical Response in Few-layer Black Phosphorus
• Semiconducting layered blue phosphorus: A computational study
• Anisotropic in-plane thermal conductivity observed in few-layer black phosphorus
• Long-range correlation energy calculated from coupled atomic response functions
• Unified theory of thermal transport in crystals and disordered solids
• 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
• Phononic thermal properties of two-dimensional materials
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials
• Anisotropic Thermal Conductivity of Exfoliated Black Phosphorus
• 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
• Unexpectedly Large Tunability of Lattice Thermal Conductivity of Monolayer Silicene via Mechanical Strain
• Thermal conductivity decomposition in two-dimensional materials: Application to graphene
• Optimizing many-body atomic descriptors for enhanced computational performance of machine learning based interatomic potentials
• Divergence of the Thermal Conductivity in Uniaxially Strained Graphene
• 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
• Phonon Localization in One-Dimensional Quasiperiodic Chains

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