Molecular Dynamics Simulation of Phosphorus Trichloride (PCl3) and Phosphorus Triiodide (PI3) Liquids by Using the Force-fields Derived from a Quantum Chemical Approach

In this study, phosphorus trichloride (PCl3) and phosphorus triiodide (PI3), as the condensed inorganic materials, were investigated based on the molecular dynamics simulation. To this purpose, molecular dynamics (MD) simulations were performed by applying force field parameters derived from a quantum chemistry approach. The potential energies were computed at the B3LYP/6-31+G(d) and B3LYP/dgdzvpd levels of theory for different configurations of PCl3 and PI3, respectively. To determine force field parameters, a four-site all-atom force field model was used to correlate the potential energy data. Therefore, the force field parameters were applied to perform the molecular dynamics simulations. The MD simulations were performed to obtain the atomic number density, enthalpy, heat capacity, and radial distribution function in the NPT and NVT ensembles for PCl3 and PI3 dimers. There is a good consistency between the experimental data and simulation results over a wide range of experimental conditions.