DFT Method Used for Prediction of Molecular and Electronic Structures of Mn(VI) Macrocyclic Complexes with Porhyrazine/Phthalocyanine and Two Oxo Ligands
1
Kazan Department of Joint Supercomputer Center of Russian Academy of Sciences—Branch of Federal Scientific Center “Scientific Research Institute for System Analysis of the RAS”, Lobachevskii Street 2/31, 420111 Kazan, Russia
2
Department of Analytical Chemistry, Certification and Quality Management, Kazan National Research Technological University, K. Marx Street 68, 420015 Kazan, Russia
*
Author to whom correspondence should be addressed.
Materials 2023, 16(6), 2394; https://doi.org/10.3390/ma16062394
Submission received: 6 February 2023
/
Revised: 4 March 2023
/
Accepted: 15 March 2023
/
Published: 16 March 2023
(This article belongs to the Special Issue Density Functional Theory Application on Chemical Calculation)
Abstract
:By using the data of the DFT quantum chemical calculation in the OPBE/TZVP and B3PW91/TZVP levels, the possibility of the existence of a manganese(VI) heteroligand complex containing porphyrazine or its tetra[benzo] derivative (phthalocyanine) and two oxygen (O2−) ligands, which is still unknown for this element, is shown. The parameters of the molecular structure, multiplicity of the ground state, NBO analysis data and standard thermodynamic parameters (enthalpy ΔH0f, entropy S0f and Gibbs’s energy ΔG0f of formation) of each of these metal macrocyclic compounds are presented. Additionally, it is noted that, based on the totality of structural data obtained by the above versions of the DFT method, the existence of a similar complex of manganese with di[benzo] derivative of porhyrazine and two oxygen (O2−) ligands seems doubtful.
1. Introduction
In studies [1,2,3,4,5,6], a quantum-chemical calculation of the parameters of molecular structures of coordination compounds having [M(P)(O)2] (I), [M(dbP)(O)2] (II) and [M(Pc)(O)2] (III) formulas [where M is Cr, Fe, Co, Ni, and Cu and P2−, dbP2− and Pc2− are double deprotonated forms of porphyrazine H2P, trans-di[benzo]porphyrazine H2dbP and phthalocyanine H2Pc, respectively] (Scheme 1) was carried out.
According to the data from these studies, obtained using two versions of the DFT method, namely DFT B3PW91/TZVP and DFT OPBE/TZVP, the possibility of the existence of all three types of complexes I–III, in the case of M=Cr, Fe and Co [1,2,3,4,5], but only of the same type, namely II in the case of M=Cu [6], was shown. In this regard, it should be noted that there are contradictory results for M=Cu: according to the DFT OPBE/TZVP data, complexes I and III can also exist, whereas, according to the DFT B3PW91/TZVP data, they cannot. At the same time, for Cr and Co, the oxidation state in each of these compounds is VI (at least formally), which is the maximum among the reliably established values of this parameter for these chemical elements. In this connection, it seems interesting to consider whether such complexes can be formed in the case of M=Mn, where the oxidation state VI is intermediate between the minimum (–III) and maximum (VII) reliably established for the given 3d element. What is remarkable is that any information about manganese compounds I, II and III in the literature devoted to macrocyclic ligands, such as porphyrins, porphyrazines and their various substituted compounds [7,8,9,10,11,12,13,14], and also anywhere else, is absent. Nevertheless, these objects are of enough considerable interest for preparative coordination chemistry and the chemistry of macrocyclic compounds. In addition, they may be helpful from a purely practical point of view, since in principle they might be employed at least as potential catalysts for diverse reactions of inorganic and organic synthesis. In this connection, our study was devoted to establishing the possibility of the existence of complexes having I, II and III formulas using modern methods of quantum-chemical calculation and the density functional theory (DFT) and, in the case of a positive answer to this question, the determination of the quantitative characteristics of their molecular and electronic structures.
2. Method
As in our studies [1,2,3,4,5,6] cited above, and also in the earlier ones [15,16,17], the DFT B3PW91/TZVP with a combination of B3PW91 functional [18,19] and TZVP basis set [20,21] was used in this study. According to data [22], this DFT version has a minimum of “normal error” compared with other DFT versions. Such conclusion was confirmed with a comparison of the results of determination of the parameters of molecular structure for 3d-element macrocyclic complexes with phthalocyanine, obtained using different DFT versions, with the experimental values of these parameters. Additionally, for comparison as in our studies [23,24,25], quantum-chemical calculations were done using the DFT OPBE/TZVP version combining TZVP and the OPBE functional [26,27]. For the 3d-element coordination compounds, the indicated DFT version more adequately predicted the relative energy stabilities of high-spin and low-spin states and also reliably described key geometric parameters of molecular structures of these compounds [27,28,29,30,31]. Calculations were performed using Gaussian09 software [32]. The data obtained as a result of calculations were visualized using the ChemCraft 1.8 program. The correspondence of the discovered stationary points to energy minima was proved in all cases using the calculation of energy-second derivatives with respect to atom coordinates, whereinto all equilibrium structures corresponding to minima of the potential energy surfaces had only real positive frequency values. In accordance with the theory of the structure of the atoms, Mn(VI), which is in the complexes I–III, must have 3d1 electronic configuration, and that is why spin multiplicities (MS) 2 and 4 were considered for the central ion indicated above. Among the structures optimized at such multiplicities, the one with the lowest total energy was selected. To calculate the parameters of molecular structures with multiplicity greater than 1, we used the unrestricted method (UB3PW91 and UOPBE). In addition, the energetically most favorable structure was always checked according to the STABLE = OPT procedure, while what was important, the wave function corresponding to this structure, was stable in all cases. Natural Bond Orbital (NBO) analysis was carried out using NBO 3.1 version in the framework of Gaussian09 software [32], according to the methodology submitted by [33]. The standard thermodynamic parameters of a formation (ΔH0f, S0f and ΔG0f) for the metal macrocyclic compounds indicated above were computed employing the methodology [34].
3. Results and Discussion
According to the data obtained as a result of our quantum-chemical calculation carried out using the B3PW91/TZVP as well as OPBE/TZVP methods, one can state with sufficient certainty the existence of two of the three complexes of types I–III indicated above, where M=Mn, namely [Mn(P)(O)2] and [Mn(Pc)(O)2]. In the case of complex II, the result turned out to be ambiguous, since according to the data obtained with the first of these two DFT methods, a complex with such a structure can exist as an isolated molecule, while according to the data obtained with the second method, it cannot. (In this regard, it is interesting that in the case of complexes similar in composition, where M=Cu, only complex II is capable of existing as an isolated molecule, while complexes I and III, according to the results of calculations with the DFT B3PW91/TZVP method, are not capable [6]). Taking into account this important fact, we will discuss further only the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes, since their existence is predicted using both of these DFT methods.
The chemical bond lengths between atoms and bond angles for complexes I and III, calculated using each of the DFT methods indicated above, are presented in Table 1. The images of the molecular structure of the macrocyclic compounds under examination, obtained with the DFT B3PW91/TZVP method, are shown in Figure 1. In this regard, we would like to note that the images of the structures of each of these two complexes obtained with the OPBE/TZVP method are very similar to images calculated with the B3PW91/TZVP method (see Supplementary Materials, Figure S1). As may be seen from the presented data, the key bond lengths, namely (Mn1N1), (Mn1N2), (Mn1N3) and (Mn1N4) in the [Mn(Pc)(O)2] macrocyclic complex, are equal to each other in the framework of the DFT B3PW91/TZVP method and within the DFT OPBE/TZVP method; however, in the [Mn(P)(O)2] macrocyclic complex, they are equal to each other only in pairs. This difference can be explained by the fact that, as is easy to see even with a cursory glance at the structural formula [Mn(P)(O)2], the nitrogen atoms in the MnN4 chelate node are not completely equivalent to each other, because two of them are bonded to neighboring atom carbons only with single bonds, while the other two have one single and one double bond. However, a very noticeable (more than 10 pm) difference in the lengths of the bonds (Mn1O1) and (Mn1O2), which takes place in the same complex according to the DFT B3PW91/TZVP data, while according to the DFT OPBE/TZVP data, the lengths of these bonds are almost the same, attracts attention (Table 1). In this regard, it should be noted that a similar inequality was noted earlier in other complexes of 3d elements formed by porphyrazine and two axial oxo ligands [1,5]. The reason for the inequality between the indicated bond lengths remains unclear. It is interesting that such an inequality is not observed if fluoride anions act as axial ligands [35,36].
The MnN4 chelate node, in the case of [Mn(P)(O)2], according to the data of both DFT methods used in this article, is almost flat (its deviation from coplanarity does not exceed 0.5°), and in the case of [Mn(Pc)(O)2], it is perfectly flat. At the same time, characteristically, all the angles in this chelate node still differ from 90° in the first of these complexes, although slightly, and then in the second, they are equal to 90° (Table 1). All four 6-membered metal chelate rings as well as all four 5-membered non-chelate rings that are adjacent to 6-membered metal chelate rings contain one nitrogen atom.
Carbon atoms also have either a perfectly flat structure or practically no different from it. Such a conclusion can be made if we take into account the sums of bond angles in each of these structural fragments, which are either equal to 360.0°, 720.0° and 540.0° or very close to these values. In addition, as can be seen from Table 1, 6-membered chelate rings as well as 5-membered non-chelate rings are almost identical to each other in terms of bond lengths between the corresponding atoms and in terms of the assortment of bond angles contained in them. According to the DFT B3PW91/TZVP method, the oxygen and manganese atoms in both [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes are under consideration.
When forming an angle (O1Mn1O2) of 180.0° between them, according to the DFT OPBE/TZVP method, the specified value of the angle (O1Mn1O2) takes place only in the case of [Mn(Pc)(O)2], while for [Mn(P)(O)2] it is more than 10° less (Table 1). As for the bond angles (OMnN) formed by the “axial” oxygen atoms, the Mn atom and the nitrogen atoms entering the MnN4 chelate node, in the case of [Mn(Pc)(O)2], they are equal to each other, and each of them is equal to 90°, while in the case of [Mn(P)(O)2], they are equal only in pairs, four of them are less than 90°, and four are more than 90° (Table 1). Despite this, as well as the above differences between the structural data obtained with these two independent versions of the DFT method, we can still state that, on the whole, these data are in fairly good agreement with each other. Concluding the discussion of the structural data, we noted that the electric moments of dipole (dipole moments) of the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes studied by us, according to the results of calculations using the DFT B3PW91/TZVP method, are 0.18 and 0.00 Debye units, respectively, and according to the DFT OPBE/TZVP method, they are 0.28 and 0.03 Debye units, respectively. As can be seen from these values, they are quite small, which is in rather good agreement with the key parameters of their structures presented in Table 1.
The ground state of both complexes considered by us in the framework of the DFT OPBE/TZVP method is a spin doublet. According to this method, the nearest excited state (spin quartet) is significantly higher in energy than the ground state (by 140.4 kJ/mol in the case of [Mn(P)(O)2] and by 96.5 kJ/mol in the case of [Mn(Pc)(O)2]). Such a situation seems quite natural for tetragonal (pseudo-octahedral) complexes with the 3d1 electron configuration and the coordination number of the central metal ion equal to 6, to which [Mn(P)(O)2] and [Mn(Pc)(O)2] belong. The B3PW91/TZVP method gives somewhat different results; according to it, the ground state of [Mn(Pc)(O)2] is also a spin doublet (the nearest excited state lies higher by 47.4 kJ/mol), while the ground state of [Mn(P)(O)2] is a spin quartet (although the energy distance between it and the nearest doublet state is only 5.2 kJ/mol). However, according to numerous statistical data for complexes with a configuration of d1, which have a coordination number of 6 and a symmetry group of D4h, C2h or C2v, in our opinion, the ground state is precisely the spin doublet, because the result obtained with the OPBE/TZVP method seems in our case to be more reliable. Images of higher occupied (HOMO) and lower unoccupied (LUMO) molecular orbitals of the studied Mn(VI) macrocyclic compounds are shown in Figure 2.
The key data of NBO analysis for complexes I and III are presented in Table 2, and full data can be found in the Supplementary Materials. It is noteworthy that the values of the effective charge on the Mn atoms in both these macrocyclic compounds, according to the NBO analysis, both in the framework of the DFT B3PW91/TZVP method and in the framework of the DFT OPBE/TZVP method, are much less than +1.00. Furthermore, the effective charges on the oxygen atoms are very different from the values (−2.00). The given fact, as well as the values of effective charges on other atoms present in these complexes (in particular, on donor nitrogen atoms N1, N2, N3 and N4), are a direct indication of a very high degree of electron density delocalization in these metal macrocyclic compounds. The distribution patterns of the spin density in the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes are shown in Figure 3.
The standard thermodynamic parameters of formation in the framework of the isobaric process for the given macrocyclic metal chelate under examination (ΔH0f, S0f and ΔG0f) are presented in Table 3. As can be seen from the data of this Table, all values of Gibbs’ energy ΔG0f are positive. It follows that neither of these two complexes can be obtained from simple substances formed by the chemical elements that make up these complexes (i.e., C, N, O and Mn). Nevertheless, judging by the data of our DFT quantum chemical calculation, the molecular structure of [Mn(P)(O)2] and [Mn(Pc)(O)2] and the full set of their geometric parameters can be realized as a whole, and therefore both of them are able to exist (at least in the gas phase in the form of separate isolated molecules).
4. Conclusions
According to data obtained as a result of our quantum-chemical calculation, both DFT versions indicated above, namely OPBE/TZVP and B3PW91/TZVP, unambiguously predict the possibility of the formation of heteroligand macrocyclic complexes [Mn(P)(O)2] (I) and [Mn(Pc)(O)2] (III) with doubly deprotonated forms of porphyrazine and its tetra[benzo] derivative (phthalocyanine); however, on the whole, they do not unambiguously confirm the possibility of the existence of an analogous complex with a doubly deprotonated form of trans-di[benzo]porphyrazine (II). In this connection, attention is drawn to the fact that, in the case of M=Cu, the reverse situation takes place: both of these methods confirm the existence of a complex of a similar type with trans-di[benzo]porphyrazine but do not unequivocally confirm the possibility of the formation of complexes with porphyrazine and phthalocyanine. With this in mind, it can be assumed that the relative stability of type II complexes in the Cr–Cu series increases, while the relative stability of complexes I and III decreases.
In the studied complexes, Mn is bonded to four nitrogen atoms and two oxygen atoms. In accordance with the generally accepted definition of the term “oxidation degree”, which was given in Refs. [1,6], we can rightly assume that the oxidation state of the central atom, namely Mn, in each of the complexes [Mn(P)(O)2] and [Mn(Pc)(O)2] is +6. The point is that in both of these complexes, the Mn atom forms six bonds by the exchange mechanism with atoms that have a higher electronegativity compared with its own, namely two double Mn=O bonds with two oxygen atoms and two single Mn–N bonds. Since the oxidation state of any chemical element is defined as the modulus of the oxidation degree and is displayed by the corresponding Roman numeral, the oxidation state Mn is exactly equal to VI in both [Mn(P)(O)2] and [Mn(Pc)(O)2]. The fact that, according to the NBO analysis data, the real charge on the manganese atom in each of these complexes differs significantly from the value of +6.00 Å, is not significant in determining the degree (state) of oxidation, since this parameter is in no way related to the above definition. Be that as it may, in connection with the foregoing, it should now confirm the possibility of the formation of these two compounds in a real experiment.
Supplementary Materials
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma16062394/s1, NBO analysis data for [Mn(P)O2]; NBO analysis data for [Mn(Pc)O2]; Figure S1: Molecular structures of the [Mn(P)(O)2] (a) and [Mn(Pc)(O)2] (b) complexes obtained as a result of DFT OPBE/TZVP quantum-chemical calculation.
Author Contributions
Conceptualization, O.V.M.; methodology, O.V.M. and D.V.C.; software, D.V.C.; validation, O.V.M. and D.V.C.; formal analysis, O.V.M. and D.V.C.; investigation, O.V.M. and D.V.C.; resources, D.V.C.; data curation, D.V.C.; writing—original draft preparation, O.V.M. and D.V.C.; writing—review and editing, O.V.M.; visualization, O.V.M. and D.V.C.; supervision, O.V.M.; project administration, O.V.M.; funding acquisition, D.V.C. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
This study did not require institutional approval.
Informed Consent Statement
Not applicable.
Data Availability Statement
No unpublished data were created or analyzed in this article.
Acknowledgments
All quantum-chemical calculations were performed at the Joint Supercomputer Center of Russian Academy of Sciences—Branch of Federal Scientific Center “Scientific Research Institute for System Analysis of the RAS”, and we acknowledge their technical support. The contribution of author Denis V. Chachkov was funded by the state assignment to the Federal Scientific Center “Scientific Research Institute for System Analysis of the RAS” for scientific research. This study was carried out using the equipment of the Center for Collective Use “Nanomaterials and Nanotechnology” at the Kazan National Research Technological University, with the financial support of the Ministry of Science and Higher Education of the Russian Federation under agreement No. 075-15-2021-699.
Conflicts of Interest
The authors declare that they have no conflicts of interest, financial or otherwise.
References
- Mikhailov, O.V.; Chachkov, D.V. M(VI) Oxidation State Stabilization in Iron, Cobalt and Nickel Heteroligand Metal Chelates Containing 3,7,11,15-Tetraazaporphine and Two Axial Oxo Ligands: Quantum-Chemical Simulation. Int. J. Mol. Sci. 2020, 21, 1494. [Google Scholar] [CrossRef] [Green Version]
- Mikhailov, O.V.; Chachkov, D.V. Stabilization of dioxochromium(VI) in the complex with tetra[benzo]porphyrazine and two oxo ligands: DFT quantum-chemical consideration. Eur. Chem. Bull. 2020, 9, 416–419. [Google Scholar] [CrossRef]
- Chachkov, D.V.; Mikhailov, O.V. Heteroligand complexes of chromium, manganese, and iron with trans-dibenzoporphyrazine and two oxo ligands: DFT calculations. Russ. Chem. Bull. 2022, 71, 656–664. [Google Scholar] [CrossRef]
- Chachkov, D.V.; Mikhailov, O.V. DFT Quantum-chemical prediction of molecular structure of iron(VI) macrocyclic complex with phthalocyanine and two oxo ligands. J. Porph. Phthalocyanines 2022, 26, 367–375. [Google Scholar] [CrossRef]
- Mikhailov, O.V.; Chachkov, D.V. DFT Quantum-Chemical Modeling Molecular Structures of Cobalt Macrocyclic Complexes with Porphyrazine or Its Benzo-Derivatives and Two Oxygen Acido Ligands. Int. J. Mol. Sci. 2020, 21, 9085. [Google Scholar] [CrossRef] [PubMed]
- Mikhailov, O.V.; Chachkov, D.V. Copper macrocyclic complex with trans-di[benzo]porphyrazine and two oxo ligands: DFT quantum-chemical design. J. Porph. Phthalocyanines 2022, 26, 180–185. [Google Scholar] [CrossRef]
- Kasuda, K.; Tsutsui, M. Some new developments in the chemistry of metallophthalocyanines. Coord. Chem. Rev. 1980, 32, 67–95. [Google Scholar] [CrossRef] [Green Version]
- Thomas, A.L. Phthalocyanines. Research & Applications; CRC Press: London, UK, 1990. [Google Scholar]
- Sliva, W.; Mianovska, B. Metalloporphyrin arrays. Transit. Met. Chem. 2000, 25, 491–504. [Google Scholar] [CrossRef]
- Spasojević, I.; Ines Batinić-Haberle, I. Manganese(III) complexes with porphyrins and related compounds as catalytic scavengers of superoxide. Inorg. Chim. Acta 2001, 317, 230–242. [Google Scholar] [CrossRef]
- Mamardashvili, G.M.; Mamardashvili, N.Z.; Koifman, O.I. Self-assembling systems based on porphyrins. Russ. Chem. Rev. 2008, 77, 59–75. [Google Scholar] [CrossRef]
- Donzello, M.P.; Ercolani, C.; Novakova, V.; Zimcik, P.; Stuzhin, P.A. Tetrapyrazinoporphyrazines and their metal derivatives. Part I: Synthesis and basic structural information. Coord. Chem. Rev. 2016, 309, 107–179. [Google Scholar] [CrossRef]
- Lomova, T.N. Axial Coordinated Metal Porphyrins in Science and Practice; URSS: Moscow, Russia, 2018; 700p. (In Russian) [Google Scholar]
- Khelevina, O.G.; Malyasova, A.S. 40 years with porphyrazines. J. Porph. Phthalocyanines 2019, 23, 1251–1264. [Google Scholar] [CrossRef] [Green Version]
- Mikhailov, O.V.; Chachkov, D.V. (H,H)-Isomerism of cis- and trans-di[benzo]porphyrazines: Quantum chemical modeling within the framework of the DFT method. J. Porph. Phthalocyanines 2021, 25, 858–865. [Google Scholar] [CrossRef]
- Chachkov, D.V.; Mikhailov, O.V. Nickel macrocyclic complexes with porphyrazine and some [benzo]substituted, oxo and fluoro ligands: DFT analysis. J. Porph. Phthalocyanines 2022, 26, 222–231. [Google Scholar] [CrossRef]
- Mikhailov, O.V.; Chachkov, D.V. New heteroligand complex of cobalt with phthalocyanine, oxo and fluoro ligands: DFT consideration. J. Porph. Phthalocyanines 2022, 26, 316–324. [Google Scholar] [CrossRef]
- Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
- Perdew, J.P.; Burke, K.; Wang, Y. Generalized gradient approximation for the exchange-correlation hole of a many-electron system. Phys. Rev. B 1996, 54, 16533–16539. [Google Scholar] [CrossRef] [Green Version]
- Schaefer, A.; Horn, H.; Ahlrichs, R. Fully optimized contracted Gaussian basis sets for atoms Li to Kr. J. Chem. Phys. 1992, 97, 2571–2577. [Google Scholar] [CrossRef] [Green Version]
- Schaefer, A.; Huber, C.; Ahlrichs, R. Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr. J. Chem. Phys. 1994, 100, 5829–5835. [Google Scholar] [CrossRef]
- Medvedev, M.G.; Bushmarinov, I.S.; Sun, J.; Perdew, J.P.; Lyssenko, K.A. Density functional theory is straying from the path toward the exact functional. Science 2017, 355, 49–52. [Google Scholar] [CrossRef]
- Chachkov, D.V.; Mikhailov, O.V. Quantum-chemical calculation of molecular structures of (5656)macrotetracyclic 3d-metal complexes “self-assembled” in quaternary systems M(II) ion-ethanedithioamide-formaldehyde-ammonia by the density functional theory method. Russ. J. Inorg. Chem. 2014, 59, 218–223. [Google Scholar] [CrossRef]
- Chachkov, D.V.; Mikhailov, O.V. Structure of (5656)macrotetracyclic chelates in the ternary systems M(II)-ethanedithioamide-acetone (M = Mn, Fe, Co, Ni, Cu, Zn) according to DFT calculations. Russ. J. Inorg. Chem. 2013, 58, 1073–1078. [Google Scholar] [CrossRef]
- Mikhailov, O.V.; Chachkov, D.V. Copper(IV) Stabilization in Macrocyclic Complexes with 3,7,11,15-Tetraazaporphine, Its Di[benzo]- or Tetra[benzo] Derivatives and Oxide Anion: Quantum-Chemical Research. Materials 2020, 13, 3162. [Google Scholar] [CrossRef] [PubMed]
- Hoe, W.M.; Cohen, A.; Handy, N.C. Assessment of a new local exchange functional OPTX. Chem. Phys. Lett. 2001, 341, 319–328. [Google Scholar] [CrossRef]
- Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1997, 78, 1396. [Google Scholar] [CrossRef] [Green Version]
- Paulsen, H.; Duelund, L.; Winkler, H.; Toftlund, H.; Trautwein, A.X. Free Energy of Spin-Crossover Complexes Calculated with Density Functional Methods. Inorg. Chem. 2001, 40, 2201–2203. [Google Scholar] [CrossRef] [PubMed]
- Swart, M.; Groenhof, A.R.; Ehlers, A.W.; Lammertsma, K. Validation of Exchange−Correlation Functionals for Spin States of Iron Complexes. J. Phys. Chem. A 2004, 108, 5479–5483. [Google Scholar] [CrossRef]
- Swart, M.; Ehlers, A.W.; Lammertsma, K. Performance of the OPBE exchange-correlation functional. Mol. Phys. 2004, 102, 2467–2474. [Google Scholar] [CrossRef]
- Swart, M. Metal–ligand bonding in metallocenes: Differentiation between spin state, electrostatic and covalent bonding. Inorg. Chim. Acta 2007, 360, 179–189. [Google Scholar] [CrossRef]
- Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision A.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
- Weinhold, F.; Landis, C.R.; Glendening, E.D. What is NBO analysis and how is it useful? Int. Rev. Phys. Chem. 2016, 35, 399–440. [Google Scholar] [CrossRef]
- Ochterski, J.W. Thermochemistry in Gaussian; Gaussian, Inc.: Wallingford, CT, USA, 2000. [Google Scholar]
- Mikhailov, O.V.; Chachkov, D.V. About possibility of stabilization of unusual copper(IV) oxidation state in complexes with porphyrazine and two fluorine ligands: Quantum-chemical design. Inorg. Chem. Commun. 2019, 106, 224–227. [Google Scholar] [CrossRef]
- Chachkov, D.V.; Mikhailov, O.V. Molecular structures of heteroligand ScIII complexes with porphyrazine, its dibenzo and tetrabenzo derivatives, and fluoride anion, as determined from DFT calculations. Russ. Chem. Bull. 2021, 70, 276–282. [Google Scholar] [CrossRef]
Scheme 1.
Structural formulas of complexes I, II and III.
Figure 1.
The pictures ofthe molecular structures of the [Mn(P)(O)2] (a) and [Mn(Pc)(O)2] (b) macrocyclic complexes obtained using the DFT B3PW91/TZVP method.
Figure 1.
The pictures ofthe molecular structures of the [Mn(P)(O)2] (a) and [Mn(Pc)(O)2] (b) macrocyclic complexes obtained using the DFT B3PW91/TZVP method.
Figure 2.
HOMO and LUMO images in the [Mn(P)(O)2] (left) and [Mn(Pc)(O)2] (right) complexes, according to the DFT B3PW91/TZVP and DFT OPBE/TZVP methods. The energies values of the given MOs (in brackets) are expressed in eV. The symbol “alpha” represents the electron with spin (+1/2), and the symbol “beta” represents the electron with spin (−l/2).
Figure 2.
HOMO and LUMO images in the [Mn(P)(O)2] (left) and [Mn(Pc)(O)2] (right) complexes, according to the DFT B3PW91/TZVP and DFT OPBE/TZVP methods. The energies values of the given MOs (in brackets) are expressed in eV. The symbol “alpha” represents the electron with spin (+1/2), and the symbol “beta” represents the electron with spin (−l/2).
Figure 3.
The pictures of spin density distribution in the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes obtained using the DFT B3PW91/TZVP (a,b) and DFT OPBE/TZVP methods (c,d).
Figure 3.
The pictures of spin density distribution in the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes obtained using the DFT B3PW91/TZVP (a,b) and DFT OPBE/TZVP methods (c,d).
Table 1.
Bond lengths and bond angles in the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes, calculated using the DFT B3PW91/TZVP and DFT OPBE/TZVP methods.
Table 1.
Bond lengths and bond angles in the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes, calculated using the DFT B3PW91/TZVP and DFT OPBE/TZVP methods.
| Macrocyclic Complex | [Mn(P)(O)2] | [Mn(Pc)(O)2] | ||
|---|---|---|---|---|
| Parameter of Molecular Structure | Calculated Using DFT B3PW91/TZVP | Calculated Using DFT OPBE/TZVP | Calculated Using DFT B3PW91/TZVP | Calculated Using DFT OPBE/TZVP |
| Mn-N bond lengths in the MnN4 chelate node, pm | ||||
| Mn1N1 | 193.2 | 195.0 | 198.2 | 198.0 |
| Mn1N2 | 196.1 | 195.0 | 198.2 | 198.0 |
| Mn1N3 | 193.2 | 196.5 | 198.2 | 198.0 |
| Mn1N4 | 196.1 | 196.5 | 198.2 | 198.0 |
| Axial Mn-O bond lengths, pm | ||||
| Mn1O1 | 187.9 | 165.1 | 164.0 | 164.9 |
| Mn1O2 | 177.7 | 165.3 | 164.0 | 164.9 |
| C-N bond lengths in 6-numbered chelate rings, pm | ||||
| N1C3 | 137.0 | 136.5 | 135.3 | 136.2 |
| N1C4 | 137.0 | 137.0 | 135.3 | 136.2 |
| N2C1 | 135.9 | 137.0 | 135.3 | 136.2 |
| N2C2 | 135.9 | 136.5 | 135.3 | 136.2 |
| N3C7 | 137.0 | 136.5 | 135.3 | 136.2 |
| N3C8 | 137.0 | 135.8 | 135.3 | 136.2 |
| N4C5 | 135.9 | 135.8 | 135.3 | 136.2 |
| N4C6 | 135.9 | 136.5 | 135.3 | 136.2 |
| N5C2 | 132.0 | 132.5 | 132.2 | 132.3 |
| N5C3 | 132.0 | 132.5 | 132.2 | 132.3 |
| N6C6 | 132.0 | 132.5 | 132.2 | 132.3 |
| N6C7 | 132.0 | 132.5 | 132.2 | 132.3 |
| N7C4 | 132.0 | 132.5 | 132.2 | 132.3 |
| N7C5 | 132.0 | 132.6 | 132.2 | 132.3 |
| N8C1 | 132.0 | 132.5 | 132.2 | 132.3 |
| N8C8 | 132.0 | 132.6 | 132.2 | 132.3 |
| C-C bond lengths in 5-numbered chelate ring (N1C4C9C10C3), pm | ||||
| C4C9 | 144.4 | 145.4 | 146.8 | 146.3 |
| C9C10 | 135.5 | 135.8 | 139.9 | 140.5 |
| C10C3 | 144.4 | 145.3 | 146.8 | 146.3 |
| Bond angles in the MnN4 chelate node, deg | ||||
| (N1Mn1N2) | 89.9 | 89.4 | 90.0 | 90.0 |
| (N2Mn1N3) | 89.9 | 89.8 | 90.0 | 90.0 |
| (N3Mn1N4) | 89.9 | 91.1 | 90.0 | 90.0 |
| (N4Fe1N1) | 89.9 | 89.8 | 90.0 | 90.0 |
| Bond angles sum (BAS), deg | 359.6 | 360.1 | 360.0 | 360.0 |
| Non-bond angles between N atoms in N4 grouping, deg | ||||
| (N1N2N3) | 89.1 | 90.7 | 90.0 | 90.0 |
| (N2N3N4) | 90.9 | 89.3 | 90.0 | 90.0 |
| (N3N4N1) | 89.1 | 89.3 | 90.0 | 90.0 |
| (N4N1N2) | 90.9 | 90.7 | 90.0 | 90.0 |
| Non-bond angles sum (NBAS), deg | 360.0 | 360.0 | 360.0 | 360.0 |
| Bond angles in 6-numbered chelate ring (Mn1N1C4N7C5N4), deg | ||||
| (Mn1N1C4) | 125.8 | 125.4 | 125.0 | 124.8 |
| (N1C4N7) | 127.9 | 128.2 | 128.2 | 127.7 |
| (C4N7C5) | 122.9 | 123.2 | 123.7 | 123.5 |
| (N7C5N4) | 127.4 | 126.7 | 128.1 | 127.7 |
| (C5N4Mn1) | 125.9 | 126.7 | 125.0 | 124.8 |
| (N4Mn1N1) | 89.9 | 89.8 | 90.0 | 90.0 |
| Bond angles sum (BAS61), deg | 719.8 | 720.0 | 720.0 | 718.5 |
| Bond angles in 6-numbered chelate ring (Mn1N4C6N6C7N3), deg | ||||
| (Mn1N4C6) | 125.9 | 124.7 | 125.0 | 124.8 |
| (N4C6N6) | 127.4 | 127.8 | 128.2 | 127.7 |
| (C6N6C7) | 122.9 | 124.1 | 123.7 | 123.5 |
| (N6C7N3) | 127.9 | 127.8 | 128.1 | 127.7 |
| (C7N3Mn1) | 125.8 | 124.7 | 125.0 | 124.8 |
| (N3Mn1N4) | 89.9 | 91.1 | 90.0 | 90.0 |
| Bond angles sum (BAS62), deg | 719.8 | 720.0 | 720.0 | 718.5 |
| Bond angles in 5-numbered ring (C3N1C4C9C10), deg | ||||
| (C3N1C4) | 108.2 | 108.0 | 110.0 | 109.8 |
| (N1C4C9) | 108.5 | 108.7 | 108.8 | 108.8 |
| (C4C9C10) | 107.4 | 107.2 | 106.2 | 106.3 |
| (C9C10C3) | 107.4 | 107.0 | 106.2 | 106.3 |
| (C10C3N1) | 108.5 | 109.1 | 108.8 | 108.8 |
| Bond angles sum (BAS51), deg | 540.0 | 540.0 | 540.0 | 540.0 |
| Bond angles in 5-numbered ring (C1N2C2C12C11), deg | ||||
| (C1N2C2) | 108.2 | 108.0 | 110.0 | 109.8 |
| (N2C2C12) | 108.9 | 109.1 | 108.8 | 108.8 |
| (C2C12C11) | 107.0 | 107.0 | 106.2 | 106.3 |
| (C12C11C1) | 107.0 | 107.2 | 106.2 | 106.3 |
| (C11C1N2) | 108.9 | 108.7 | 108.8 | 108.8 |
| Bond angles sum (BAS51), deg | 540.0 | 540.0 | 540.0 | 540.0 |
| Bond angles between O, Fe and N atoms, deg | ||||
| (O1Mn1N1) | 86.2 | 93.8 | 90.0 | 90.0 |
| (O1Mn1N2) | 88.9 | 93.8 | 90.0 | 90.0 |
| (O1Mn1N3) | 86.2 | 86.2 | 90.0 | 90.0 |
| (O1Mn1N4) | 88.9 | 86.2 | 90.0 | 90.4 |
| (O2Mn1N1) | 93.8 | 93.8 | 90.0 | 90.0 |
| (O2Mn1N2) | 91.1 | 93.8 | 90.0 | 89.6 |
| (O2Mn1N3) | 93.8 | 86.2 | 90.0 | 90.0 |
| (O2Mn1N4) | 91.1 | 86.2 | 90.0 | 90.0 |
| Bond angles between O and Mn atoms, deg | ||||
| (O1Mn1O2) | 180.0 | 169.2 | 180.0 | 180.0 |
Table 2.
NBO analysis data for the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes in the ground state according to DFT B3PW91/TZVP and DFT OPBE/TZVP.
Table 2.
NBO analysis data for the [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes in the ground state according to DFT B3PW91/TZVP and DFT OPBE/TZVP.
| Complex | Calculated Using DFT | Charges on the Atoms in Electron Charge Units (ē) | <S**2> a | ||||
|---|---|---|---|---|---|---|---|
| Mn1 | N1 (N3) | N2 (N4) | O1 | O2 | |||
| [Mn(P)(O)2] | B3PW91/TZVP | 0.2002 | −0.3617 (−0.3617) | −0.3561 (−0.3561) | −0.3698 | −0.4155 | 3.9600 |
| OPBE/TZVP | −0.0715 | −0.3051 (−0.3024) | −0.3051 (−0.3024) | −0.0999 | −0.0996 | 0.8218 | |
| [Mn(Pc)(O)2] | B3PW91/TZVP | −0.1560 | −0.3180 (−0.3180) | −0.3180 (−0.3180) | −0.2483 | −0.2483 | 0.7803 |
| OPBE/TZVP | −0.1752 | −0.2816 (−0.2816) | −0.2817 (−0.2817) | −0.1493 | −0.1500 | 0.7668 | |
a <S**2> are numerical values of the operator of the square of the intrinsic angular momentum of the total spin of the system.
Table 3.
Standard thermodynamic parameters of formation for [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes calculated using the OPBE/TZVP and B3PW91/TZVP methods.
Table 3.
Standard thermodynamic parameters of formation for [Mn(P)(O)2] and [Mn(Pc)(O)2] complexes calculated using the OPBE/TZVP and B3PW91/TZVP methods.
| Complex | DFT Version | ΔH0f, 298, kJ/mol | S0f, 298, J/mol ∙ K | ΔG0f, 298, kJ/mol |
|---|---|---|---|---|
| [Mn(P)(O)2] | OPBE/TZVP | 264.9 | 769.8 | 488.3 |
| B3PW91/TZVP | 689.7 | 759.4 | 916.2 | |
| [Mn(Pc)(O)2] | OPBE/TZVP | 151.1 | 1129.2 | 422.8 |
| B3PW91/TZVP | 756.8 | 1123.2 | 1030.3 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Chachkov, D.V.; Mikhailov, O.V. DFT Method Used for Prediction of Molecular and Electronic Structures of Mn(VI) Macrocyclic Complexes with Porhyrazine/Phthalocyanine and Two Oxo Ligands. Materials 2023, 16, 2394. https://doi.org/10.3390/ma16062394
AMA Style
Chachkov DV, Mikhailov OV. DFT Method Used for Prediction of Molecular and Electronic Structures of Mn(VI) Macrocyclic Complexes with Porhyrazine/Phthalocyanine and Two Oxo Ligands. Materials. 2023; 16(6):2394. https://doi.org/10.3390/ma16062394
Chicago/Turabian StyleChachkov, Denis V., and Oleg V. Mikhailov. 2023. "DFT Method Used for Prediction of Molecular and Electronic Structures of Mn(VI) Macrocyclic Complexes with Porhyrazine/Phthalocyanine and Two Oxo Ligands" Materials 16, no. 6: 2394. https://doi.org/10.3390/ma16062394
APA StyleChachkov, D. V., & Mikhailov, O. V. (2023). DFT Method Used for Prediction of Molecular and Electronic Structures of Mn(VI) Macrocyclic Complexes with Porhyrazine/Phthalocyanine and Two Oxo Ligands. Materials, 16(6), 2394. https://doi.org/10.3390/ma16062394
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
