The ligand field molecular mechanics (LFMM) model, which incorporates the ligand field stabilization energy (LFSE) directly into the potential energy expression of molecular mechanics (MM), has been implemented in the "chemically aware" molecular operating environment (MOE) software package. The new program, christened DommiMOE, is derived from our original in-house code that has been linked to MOE via its applications programming interface and a number of other routines written in MOE's native scientific vector language (SVL). DommiMOE automates the assignment of atom types and their associated parameters and popular force fields available in MOE such as MMFF94, AMBER, and CHARMM can be easily extended to provide a transition metal simulation capability. Some of the unique features of the LFMM are illustrated using MMFF94 and some simple [MCl)]2- and [Ni(NH3)n]2+ species. These studies also demonstrate how density functional theory calculations, especially on experimentally inaccessible systems, provide important data for designing improved LFMM parameters. DommiMOE treats Jahn-Teller distortions automatically, and can compute the relative energies of different spin states for Ni(II) complexes using a single set of LFMM parameters.

To account for the distortion of the coordination sphere that takes place in complexes containing open-shell metal cations such as Cu(II), we implemented, in sum of interactions between fragments ab initio computed (SIBFA) molecular mechanics, an additional contribution to take into account the ligand field splitting of the metal d orbitals. This term, based on the angular overlap model, has been parameterized for Cu(II) coordinated to oxygen and nitrogen ligands. The comparison of the results obtained from density functional theory computations on the one hand and SIBFA or SIBFA-LF on the other shows that SIBFA-LF gives geometric arrangements similar to those obtained from quantum mechanical computations. Moreover, the geometric improvement takes place without downgrading the energetic agreement obtained from SIBFA. The systems considered are Cu(II) interacting with six water molecules, four ammonia or four imidazoles, and four water plus two formate anions.

No abstract

No abstract