Diisopropyl Methylphosphonate and Sarin Decomposition on Pristine vs. Hydroxylated Alumina Surfaces: Mechanistic Predictions from Ab Initio Molecular Dynamics

23 May 2024, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Understanding the reactivity of chemical warfare agents on material surfaces is essential for mitigating and controlling their chemical/biological activity. Using \textit{ab initio} molecular dynamics, we compare the adsorption and decomposition of sarin and its simulant, diisopropyl methylphosphonate (DIMP), on both pristine and hydroxylated alumina surfaces at various temperatures. Our extensive calculations confirm that DIMP is a suitable simulant for sarin experimental studies on pristine alumina surfaces due to their similar adsorption configurations and reaction dynamics. Moreover, our ab initio molecular dynamics simulations predict C--O bond cleavage to be the initial decomposition step for both DIMP and sarin on the pristine surface. We also identify an additional potential decomposition route for sarin via cleavage of the P--F bond. While similarities exist on the pristine surface, differences between sarin and DIMP emerge on the hydroxylated surface at elevated temperatures. This work provides detailed insights into these atomistic decomposition mechanisms and highlights the utility of predictive quantum dynamics simulations for evaluating the suitability of simulants to guide future experimental investigations of these hazardous compounds.

Keywords

DIMP
Sarin
ab initio molecular dynamics
metadynamics
chemical warfare agents
metal oxide
catalysis
Born-Oppenheimer molecular dynamics
surfaces
quantum calculations
decomposition
diisopropyl methylphosphonate
alumina

Supplementary materials

Title
Description
Actions
Title
Supplementary Material
Description
Additional details on metadynamics simulations, optimized structures of gamma-Al2O3, collective variables used to calculate free energies, time-dependent fluctuations of the C--O bond of DIMP and sarin, additional sarin decomposition mechanism, free energy profiles for DIMP and sarin, Cartesian coordinates for structures discussed in the main text, and sample CP2K input files.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.