I am trying to model the active sites of several MOF's using Gaussian 09 Unix package. One such MOF I have to model is MIL-100 and I am constantly finding difficulties in obtaining a converged structure due to its complex active site, which include three Ferric ions and an Oxide ion. I tried to start building up the active site with only three Ferric ions and the Oxide ion at first and then adding up the organic ligands but even the four atomic system failed to converge.
Is there any specific method I should proceed in obtaining the optimized active site? I tried using a small basis sets like 3-21G but the problem persists. I am using B3LYP method.
The active site of MIL-100 is attached herewith.