Ab initio simulation of the structure and dynamics of white phosphorus (P4) at low temperatures
Ab initio simulations of low-temperature white phosphorus (β-P4) were performed using Car–Parrinello density-functional-theory molecular dynamics and a plane-wave basis set. Gradient corrections to the local density approximation were used to obtain a better description of the weak intermolecular forces. Starting with unit-cell parameters and atomic coordinates recently verified using neutron scattering, molecular dynamics trajectories over 9–10 ps were obtained for temperatures of 50 and 145 K. The computations show that the distortions from ideal tetrahedral P4 molecules found in the reduction of the neutron data may be a finite-temperature effect and not an artifact of the data analysis. Examination of the molecular reorientational dynamics shows that orientational order persists throughout the simulations, which suggests that β-P4 is an orientational crystal.