Rupture Mechanism of Aromatic Systems from Graphite Probed with Molecular Dynamics Simulations

Intermolecular interactions involving aromatic rings are of pivotal importance in many areas of chemistry, biology and materials science. Mimicking recent atomic force microscopy (AFM) experiments that measured the adhesion forces of single pi-pi complexes, here interactions between pyrene/coronene and graphite have been probed by force-probe molecular dynamics (FPMD) simulations. The pyrene or coronene molecule was connected to a virtual spring through a flexible poly(ethylene glycol) (PEG) linker and was pulled away from graphite in water under constant velocity. Pyrene and coronene showed similar unbinding pathways featuring four states, with a transition and an intermediate state connecting the bound and unbound states in terms of distance and interplanar angles. Transient conformations with tilted orientations (approximately 40 degrees) and with one side of the aromatic structure still in contact with the graphite surface (approximately 70 degrees) were identified as the transition and intermediate states, respectively, similar to previously observed perpendicularly stacked benzene dimers. The distance to transition state x(tr) was determined to be 0.23 +/- 0.03 nm both for pyrene/graphite and coronene/graphite. The complexes share similar unbinding pathways, but coronene binds to graphite more strongly than to pyrene.