Search

The properties of a new method of performing molecular dynamic simulations of complex chemical processes are presented. The method is formulated to give a time-dependent, multilevel representation of the total potential that is derived from spatially resolved quantum mechanical regions. An illustrative simulation is performed on a 110 atom system to demonstrate the continuity and energy conserving properties of the method. The effect of a discontinuous total potential upon the kinetic energy of the system is examined. The discontinuities in the magnitude of atomic force vectors due to changing the electronic structure during the simulation are examined as well as the effect that these discontinuities have upon the atomic kinetic energies. The method, while not conserving total energy, does yield canonical (NVT) simulations. The time reversibility property of the simulation with an extremely discontinuous total potential is discussed. The computational scaling associated with the formation of the spatially resolved, time-dependent groups is also investigated.

Density Functional calculations have been performed to calculate the one-electron oxidation potential for ferrocene and the redox couples for a series of small transition metal compounds of the first-, second-, and third-row elements. The solvation effects are incorporated via a self-consistent reaction field (SCRF), using the polarized continuum model (PCM). From our study of seven different density functionals combined with three different basis sets for ferrocene, we find that no density functional method can reproduce the redox trends from experiment when referencing our results to the experimental absolute standard hydrogen electrode (SHE) potential. In addition, including additional necessary assumptions such as solvation effects does not lead to any conclusion regarding the appropriate functional. However, we propose that if one references their transition metal compounds results to the calculated absolute half-cell potential of ferrocene, they can circumvent the additional assumptions necessary to predict a redox couple. Upon employing this method on several organometallic and inorganic complexes, we obtained very good correlation between calculated and experimental values (R(2) = 0.97), making it possible to predict trends with a high level of confidence. The hybrid functional B3LYP systematically underestimates the redox potential; however, the linear correlation between DFT and experiment is good (R(2) = 0.96) when including a baseline shift. This protocol is a powerful tool that allows theoretical chemists to predict the redox potential in solution of several transition metal complexes a priori and aids in the rational design of redox-active catalysts.