ABSTRACT
A multiscale density functional theory-quantum mechanics/molecular mechanics (DFT-QM/MM) scheme is presented, based on an efficient electrostatic coupling between the electronic density obtained from a grid-based projector augmented wave (GPAW) implementation of density functional theory and a classical potential energy function. The scheme is implemented in a general fashion and can be used with various choices for the descriptions of the QM or MM regions. Tests on H2O clusters, ranging from dimer to decamer show that no systematic energy errors are introduced by the coupling that exceeds the differences in the QM and MM descriptions. Over 1 ns of liquid water, Born-Oppenheimer QM/MM molecular dynamics (MD) are sampled combining 10 parallel simulations, showing consistent liquid water structure over the QM/MM border. The method is applied in extensive parallel MD simulations of an aqueous solution of the diplatinum [Pt2(P2O5H2)4]4- complex (PtPOP), spanning a total time period of roughly half a nanosecond. An average Pt-Pt distance deviating only 0.01 Å from experimental results, and a ground-state Pt-Pt oscillation frequency deviating by <2% from experimental results were obtained. The simulations highlight a remarkable harmonicity of the Pt-Pt oscillation, while also showing clear signs of Pt-H hydrogen bonding and directional coordination of water molecules along the Pt-Pt axis of the complex.
ABSTRACT
Os(II)/(III) and Co(II)/(III) polypyridine complexes in aqueous solution are robust molecular entities both in freely solute state and adsorbed on Au(111)- and Pt(111)-electrode surfaces. This class of robust coordination chemical compounds have recently been characterized by electrochemical scanning tunneling microscopy (in situ STM). The Os-complexes were found to display strong tunneling spectroscopic (STS) features at the level of resolution of the single molecule while STS features of the Co complexes, although clear, were much weaker. The data was framed by concise but phenomenological theory of interfacial electrochemical electron transfer extended to the electrochemical in situ STM configuration. With a view on first-principle insight into the in situ STM behavior of robust redox (as opposed to nonredox) molecules, we present in this report a density functional theory (DFT) study of the complexes in both free and adsorbate state, in either state exposed to both stoichiometric counterions and a large assembly of solvent water molecules. The oxidation states of the complexes were controlled, first by introducing chlorine counter atoms followed by spontaneous attraction of electrons from the complexes, also at first in electrostatically neutral form. Second, the solvent is found to provide strong dielectric screening of this charge transfer process and to be crucial for achieving the full chemically meaningful charge separated ionic oxidation states. The molecular charge and structure of the complexes in the presence of the solvent, are conserved upon adsorption, whereas the structural features of the different oxidation states are completely lost upon adsorption under vacuum conditions. Detailed microscopic insight such as offered by the present study will be important in molecular-based approaches to "smart" redox molecules enclosed in in situ STM or other nanoscale and single-molecules scale configurations in condensed matter environments.