Computational chemistry and molecular simulations of phosphoric acid

AbstractQuantum mechanics (QM) calculations, molecular dynamics (MD) simulations using the condensed‐phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field, and the atom‐centered density matrix propagation (ADMP) approach have been used to investigate properties of phosphoric acid (PA). QM using B3LYP/6‐31++G(d,p) density functional theory were used to calculate gas‐phase proton affinities and interaction energies of PA and its derivatives. Detailed single coordinate driving, followed by quadratic synchronous transit optimization was used to determine energy barriers for different proton transfer (PT) pathways. Determined energy barrier heights in ascending order are (unit: kJ/mol): H3O+→H3PO4 (0); H4P2O7→H3PO4 (2.61); H3PO4→H2PO (5.31); H4PO→H3PO4 (∼7.33); H3PO4→H4P2O7/H3PO4→H3PO4 (15.99); H4P2O7→H2O (28.61); H3PO4→H2O (47.14). The COMPASS force field was used to study condensed‐phase properties of PA. Good agreement between experimental data and MD results including density, radial distribution functions, and self‐diffusion coefficient at different temperatures provides validation of the COMPASS force field for PA. Finally, preliminary ADMP st