Journal of Theoretical Chemistry

Volume 2013 (2013), Article ID 734354, 4 pages

http://dx.doi.org/10.1155/2013/734354

## A Computational Determination of the Lowest Energy Electronic and Geometric States of First Row Transition Metal Dioxygen Dications

^{1}Environmental Resource Management, 1 Beacon Street, 5th Floor, Boston, MA 02108, USA^{2}Department of Chemistry, Biochemistry, and Physics, Marist College, 3399 North Road, Poughkeepsie, NY 12601, USA

Received 28 March 2013; Accepted 13 August 2013

Academic Editors: G. Borosky, P. Derosa, G. Pacchioni, A. Stavrakoudis, A. Tilocca, and B. M. Wong

Copyright © 2013 Jillian Lennartz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The lowest energy geometric structures and electronic spin states of first row transition metal (TM) dioxygen dication molecules ([TM–O_{2}]^{2+}; TM = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn) have been determined at the B3LYP/LANL2DZ level of theory (along with an extra -type polarization function added to the O atoms). In order to further verify the spin states, CASSCF(6 + , 9) energy points were determined ( = number of TM electrons). It has been found that with the exception of [Sc–O_{2}]^{2+}, [V–O_{2}]^{2+}, [Co–O_{2}]^{2+}, and [Ni–O_{2}]^{2+}, all [TM–O_{2}]^{2+} molecules take on a high-spin state. [Sc–O_{2}]^{2+} adopts a trigonal structure, while [Ti–O_{2}]^{2+}-[Mn–O_{2}]^{2+} are essentially linear and [Fe–O_{2}]^{2+}-[Zn–O_{2}]^{2+} are bent. It is further noted that the O–O bond decreases from 130.0 pm to 118.1 pm as the TM changes from Sc to Zn. However, the TM–O_{2} bond lengths fluctuate between values of 182.2 pm for [Ni–O_{2}]^{2+} and 232.2 pm for [Zn–O_{2}]^{2+}.

#### 1. Introduction

Transition metals (TMs) bound to an O_{2} molecule are the chemically active site in many industrial [1–5] and biological [1, 2, 6–10] molecules. Considering the TM–O_{2} bond to lie along the -axis of a Cartesian coordinate system, in TM–O_{2} complexes, bonds can form by the interactions of TM , , and atomic orbitals with suitable orbitals on the O_{2} ligand in order to make () or ( and ) bonds (Scheme 1). While there have been numerous studies of neutral TM–O_{2} [11, 12] complexes, in many TM–O_{2} containing molecules of biological and industrial importance, the TM is bound to a porphyrin ring in the 2^{−} oxidation state. Therefore, [TM–O_{2}]^{2+} molecules are the simplest approximations to this important class of TM–O_{2} containing molecules.

Probably the most studied TM–O_{2} bond is the Fe–O_{2} bond due to its role in heme-containing proteins [13–16]. This bond has been described as low-spin Fe(II) accepting an electron pair from excited singlet-state O_{2} forming the bond, while donating an electron pair to oxygen to form a bond. An alternative view involves Fe(III) interacting with O^{2−}. In a third possibility, triplet Fe(II) couples with triplet O_{2} to form a closed shell singlet. Computational results vary depending on the extent of electron correlation and basis set size. It is clear, however, that the overall bonding scheme involves some degree of both O_{2}Fe and FeO_{2} and electron transfer.

Further complicating matters, the number of near degenerate TM -orbitals leads to the possibility of many low lying spin states. In addition, there are a number of different ways that the O_{2} unit can be arranged relative to the TM. Herein, we seek to determine the lowest energy electronic and geometric states of [TM–O_{2}]^{2+} molecules (TM = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn) as a precursor to a more in-depth study of the TM–O_{2} bond.

#### 2. Theoretical Methods

In order to locate the lowest energy geometrical configuration, B3LYP [17–19] geometry optimizations were performed using the LANL2DZ [20–22] effective core potential basis set with an additional -type polarization function () on the O atoms [23]. Optimizations began with trigonal, linear, and bent starting structures as shown in Scheme 2. Optimizations at all three starting geometries were performed for singlet, triplet, and quintet spin states for [TM–O_{2}]^{2+} molecules with an even number of electrons (TM = Ti, Cr, Fe, Ni, and Zn) and doublet, quartet, and sextet spin states for molecules with an odd number of electrons (TM = Sc, V, Mn, Co, and Cu).

In order to further verify the spin state of each molecule, CASSCF [24–29] energy point calculations were performed for the various different spin states at the lowest energy B3LYP conformation. The CASSCF active space consisted of nine orbitals: the 5 TM -type orbitals along with the and O_{2} combinations (Scheme 3). Six electrons were placed into this active space from the O_{2} and orbitals, along with the TM electrons (i.e., 1 for Sc, 2 for Ti, etc.).

All calculations were carried out with the Gaussian 03 [30] suite of programs.

#### 3. Results and Discussion

The geometric parameters and spin states of the lowest energy structures of all [TM–O_{2}]^{2+} molecules studied herein are reported in Table 1.

[Sc–O_{2}]^{2+} adopts a trigonal structure as in **1** of Scheme 2. This can be thought of as the Sc^{2+} atom bonding to the system of O_{2}. The single electron of Sc^{2+} spin pairs with the single electron in the MO of O_{2} as depicted in Scheme 3 (MO , **1**), while one electron remains in the MO of O_{2} (MO , **1**) resulting in a doublet ground state.

[Ti–O_{2}]^{2+}, [V–O_{2}]^{2+}, [Cr–O_{2}]^{2+}, and [Mn–O_{2}]^{2+} are essentially linear with the O_{2} binding end-on to the TM^{2+}. The slight bending in [Ti–O_{2}]^{2+} breaks the degeneracy of the and MOs. The observed triplet is the result of the subtle balance between orbital energies and favorable electron exchange. Doublet [V–O_{2}]^{2+} has one electron in the MO, and triplet [CrO_{2}]^{2+} has one electron in each of the degenerate MOs (Scheme 3, **2**). In order for the high spin quartet to be the lowest energy state of [Mn–O_{2}]^{2+}, one electron must remain in each of the MOs, while another single electron occupies the MO (Scheme 3, **2**). This situation is only possible if the and MOs are close enough in energy that exchange stabilization outweighs the destabilization introduced by placing an electron in the higher energy MO.

Upon examination of Scheme 3, it would appear that [Fe–O_{2}]^{2+} should continue the trend of linearity by completely filling the MO thus forming a stable singlet. Indeed, a stable linear singlet was found; however, it was 18.2 kJ/mol above the bent triplet ground state. Upon bending, the linear MO (Scheme 3, **2**) splits into the bent and MOs (Scheme 3, **3**). This results in a stabilization of the MO due to increased bonding interaction, while the MO experiences little change energetically. The triplet is then formed by placing a single electron in the and MOs. The bent triplet is favored over the bent singlet due to favorable exchange between electrons in these energetically close orbitals.

Adding more electrons causes [Co–O_{2}]^{2+}, [Ni–O_{2}]^{2+}, [Cu–O_{2}]^{2+}, and [Zn–O_{2}]^{2+} to become increasingly more bent. For [Co–O_{2}]^{2+} and [Ni–O_{2}]^{2+}, the and MOs are filled one at a time resulting in doublet and singlet ground states, respectively. On the other hand, [Cu–O_{2}]^{2+} and [Zn–O_{2}]^{2+} adopt high spin quartet and triplet ground states, respectively, indicating that the , , and MOs (Scheme 3, **3**) are close in energy.

It is interesting to note that the O–O bond length decreases from a value of 130.0 pm in [Sc–O_{2}]^{2+} to 118.1 pm in [Zn–O_{2}]^{2+} (Table 1). Such a decrease could be caused by increasing the electron density in O–O bonding MOs or decreasing electron density in O–O antibonding MOs. This range is both above and below the O–O bond distance in O_{2} (121 pm [31]), but within the limits of O^{2−} (135 pm [31], an extra electron in the antibonding MO) and O^{2+} (112 pm [31], one less electron in the antibonding MO). However, the TM–O_{2} bond lengths are not so consistent and fluctuate between values of 182.2 pm in [Ni–O_{2}]^{2+} and 232.2 pm in [Zn–O_{2}]^{2+} (Table 1).

#### 4. Conclusions

In the current study, we present the lowest energy geometrical conformations and spin states of first row [TM–O_{2}]^{2+} molecules. [Sc–O_{2}]^{2+} is trigonal, while [Ti–O_{2}]^{2+}-[Mn–O_{2}]^{2+} are linear and [Fe–O_{2}]^{2+}-[Zn–O_{2}]^{2+} are bent. The O–O bond length decreases as the TM goes from Sc to Zn, whereas the TM–O_{2} bond length fluctuates. For [Mn–O_{2}]^{2+}, [Fe–O_{2}]^{2+}, [Cu–O_{2}]^{2+}, and [Zn–O_{2}]^{2+} molecules, high spin states are lowest in energy indicating degenerate or near-degenerate orbitals.

Although it is clear that there is some level of electron transfer from the O_{2} to the TM as well as some degree of TM to O_{2} back donation, a more detailed analysis is required in order to fully understand the bonding in these unique molecules and how that extends to larger molecules of biological and industrial importance. Furthermore, we have shown that while the current results can be used as a starting point for studies of TM–O_{2} containing molecules with donor and acceptor ligands, caution must be exercised as both high- and low-spin ground states are possible.

#### Acknowledgments

This work was supported by the Research Corporation Cottrell College Science Award CC6519. The authors thank the School of Science and the office of the Academic Vice President of Academic Affairs at Marist College for continued support of research.

#### References

- A. L. Feig and S. J. Lippard, “Reactions of non-heme iron(II) centers with dioxygen in biology and chemistry,”
*Chemical Reviews*, vol. 94, no. 3, pp. 759–805, 1994. View at Google Scholar · View at Scopus - N. Kitajima and Y. Moro-oka, “Copper-dioxygen complexes. Inorganic and bioinorganic perspectives,”
*Chemical Reviews*, vol. 94, no. 3, pp. 737–757, 1994. View at Google Scholar · View at Scopus - J. Muzart, “Chromium-catalyzed oxidations in organic synthesis,”
*Journal of Chemical Reviews*, vol. 92, no. 1, pp. 113–140, 1992. View at Publisher · View at Google Scholar - D. Dolphin, T. G. Traylor, and L. Y. Xie, “Polyhaloporphyrins: unusual ligands for metals and metal-catalyzed oxidations,”
*Accounts of Chemical Research*, vol. 30, no. 6, pp. 251–259, 1997. View at Google Scholar · View at Scopus - D. H. Busch and N. W. Alcock, “Iron and cobalt “lacunar ” complexes as dioxygen carriers,”
*Chemical Reviews*, vol. 94, no. 3, pp. 585–623, 1994. View at Google Scholar · View at Scopus - G. H. Loew and D. Harris, “Role of the heme active site and protein environment in structure, spectra, and function of the cytochrome P450s,”
*Chemical Reviews*, vol. 100, no. 2, pp. 407–420, 2000. View at Publisher · View at Google Scholar - A. Butler, M. J. Clague, and G. E. Meister, “Vanadium peroxide complexes,”
*Chemical Reviews*, vol. 94, no. 3, pp. 625–638, 1994. View at Google Scholar · View at Scopus - K. Shikama, “The molecular mechanism of autoxidation for myoglobin and hemoglobin: A venerable puzzle,”
*Chemical Reviews*, vol. 98, no. 4, pp. 1357–1374, 1998. View at Publisher · View at Google Scholar - M. H. Dickman and M. T. Pope, “Peroxo and superoxo complexes of chromium, molybdenum, and tungsten,”
*Chemical Reviews*, vol. 94, no. 3, pp. 569–584, 1994. View at Google Scholar · View at Scopus - J. P. Collman and L. Fu, “Synthetic models for hemoglobin and myoglobin,”
*Accounts of Chemical Research*, vol. 32, no. 6, pp. 455–463, 1992. View at Publisher · View at Google Scholar - E. L. Uzunova, H. Mikosh, and G. St. Nikolov, “Electronic structure of oxide, peroxide, and superoxide clusters of the 3d elements: a comparative density functional study,”
*Journal of Chemical Physics*, vol. 128, Article ID 094307, 2008. View at Publisher · View at Google Scholar - Y. Gong, M. Zhou, and L. Andrews, “Spectroscopic and theoretical studies of transition metal oxides and dioxygen complexes,”
*Chemical Reviews*, vol. 109, no. 12, pp. 6765–6808, 2009. View at Publisher · View at Google Scholar · View at Scopus - I. Bytheway and M. B. Hall, “Theoretical calculations of metal-dioxygen complexes,”
*Chemical Reviews*, vol. 94, no. 3, pp. 639–658, 1994. View at Google Scholar · View at Scopus - B. D. Olafson and W. A. Goddard III, “Molecular description of dioxygen bonding in hemoglobin,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 74, no. 4, pp. 1315–1319, 1977. View at Google Scholar - W. A. Goddard and B. D. Olafson, “Ozone model for bonding of an O
_{2}to heme in oxyhemoglobin,”*Proceedings of the National Academy of Sciences of the United States of America*, vol. 72, no. 6, pp. 2335–2339, 1975. View at Google Scholar - H. Chen, M. Ikeda-Saito, and S. Shaik, “Nature of the Fe-O
_{2}bonding in oxy-myoglobin: effect of the protein,”*Journal of the American Chemical Society*, vol. 130, no. 44, pp. 14778–14790, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,”
*Journal of Chemical Physics*, vol. 98, pp. 5648–5653, 1993. View at Publisher · View at Google Scholar - B. Miehlich, A. S. Savin, and H. Preuss, “Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr,”
*Chemical Physics Letters*, vol. 157, no. 3, pp. 200–206, 1989. View at Google Scholar - C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation energy formula into a functional of the electron density,”
*Physical Review B*, vol. 37, pp. 785–789, 1988. View at Google Scholar - P. J. Hay and W. R. Wadt, “Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg,”
*Journal of Chemical Physics*, vol. 82, p. 270, 1985. View at Publisher · View at Google Scholar - W. R. Wadt and P. J. Hay, “Ab initio effective core potentials for molecular calculations. Potentials for main group elements Na to Bi,”
*The Journal of Chemical Physics*, vol. 82, no. 1, pp. 284–298, 1985. View at Google Scholar · View at Scopus - P. J. Hay and W. R. Wadt, “Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals,”
*Journal of Chemical Physics*, vol. 82, pp. 299–311, 1985. View at Publisher · View at Google Scholar - C. E. Check, T. O. Faust, T. M. Bailey, B. J. Wright, T. M. Gilberrt, and L. S. Snderlin, “Addition of polarization and diffuse functions to the LANL2DZ basis set for P-block elements,”
*The Journal of Physical Chemistry A*, vol. 105, no. 34, pp. 8111–8116, 2001. View at Publisher · View at Google Scholar - D. Hegarty and M. A. Robb, “Application of unitary group methods to configuration interaction calculations,”
*Molecular Physics*, vol. 38, no. 6, pp. 1795–1812, 1979. View at Publisher · View at Google Scholar - R. H. E. Eade and M. A. Robb, “Direct minimization in mc scf theory. The quasi-newton method,”
*Chemical Physics Letters*, vol. 83, no. 2, pp. 362–368, 1981. View at Google Scholar - H. B. Schlegel and M. A. Robb, “MC SCF gradient optimization of the ${\text{H}}_{2}\text{CO}\to {\text{H}}_{2}+\text{CO}$ transition structure,”
*Chemical Physics Letters*, vol. 93, no. 1, pp. 43–46, 1982. View at Google Scholar - F. Bernardi, A. Bottini, J. J. W. McDougall, M. A. Robb, and H. B. Schlegel, “MCSCF gradient calculation of transition structures in organic reactions,”
*Faraday Symposia of the Chemical Society*, vol. 19, pp. 137–147, 1984. View at Publisher · View at Google Scholar - N. Yamamoto, T. Vreven, M. A. Robb, M. J. Frisch, and H. B. Schlegel, “A direct derivative MC-SCF procedure,”
*Chemical Physics Letters*, vol. 250, no. 3-4, pp. 373–378, 1996. View at Google Scholar - M. J. Frisch, I. N. Ragazos, M. A. Robb, and H. B. Schlegel, “An evaluation of three direct MC-SCF procedures,”
*Chemical Physics Letters*, vol. 189, no. 6, pp. 524–528, 1992. View at Google Scholar - Gaussian 03, Revision D. 02, M. J. Frisch et al., Wallingford, CT, 2003.
- NIST Chemistry WebBook, 2013, http://webbook.nist.gov/chemistry/.