Next Article in Journal
Study of Interactions between Amlodipine and Quercetin on Human Serum Albumin: Spectroscopic and Modeling Approaches
Previous Article in Journal
Antiradical and Antioxidative Activity of Azocalix[4]arene Derivatives: Combined Experimental and Theoretical Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 24
Issue 3
10.3390/molecules24030486
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polar Metallocenes

by
Haiwu Zhang
1,
B.Yu. Yavorsky
1 and
R.E. Cohen
1,2,*
1
Department of Earth and Environmental Sciences, Ludwig-Maximilians Universität München, Theresienstr, 41 80333 Munich, Germany
2
Extreme Materials Initiative, Geophysical Laboratory, Carnegie Institution for Science, Washington, DC 20015, USA
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(3), 486; https://doi.org/10.3390/molecules24030486
Submission received: 24 December 2018 / Revised: 23 January 2019 / Accepted: 26 January 2019 / Published: 29 January 2019
(This article belongs to the Section Physical Chemistry)
Browse Figures
Versions Notes

Abstract

:
Crystalline polar metallocenes are potentially useful active materials as piezoelectrics, ferroelectrics, and multiferroics. Within density functional theory (DFT), we computed structural properties, energy differences for various phases, molecular configurations, and magnetic states, computed polarizations for different polar crystal structures, and computed dipole moments for the constituent molecules with a Wannier function analysis. Of the systems studied, Mn2(C9H9N)2 is the most promising as a multiferroic material, since the ground state is both polar and ferromagnetic. We found that the predicted crystalline polarizations are 30–40% higher than the values that would be obtained from the dipole moments of the isolated constituent molecules, due to the local effects of the self-consistent internal electric field, indicating high polarizabilities.

Graphical Abstract

1. Introduction

Metallocenes have manifold uses ranging from fire retardants, catalysts, and cancer pharmaceuticals [1,2,3]. Almost all applications of metallocenes are for their molecular properties; however, here we consider crystalline metallocenes and their potential as active, multiferroic, materials.
The discovery of ferrocene, which has two cyclopentadienyl (Cp, or C5H5−) rings with a transition metal sandwiched between, was the beginning of the field of organometallic chemistry, with the 1973 Nobel prize awarded to Geoffrey Wilkinson and Ernst Otto Fisher [4,5]. There are many possible metallocenes, since the charge of the cyclic conjugated ligands -CnHnq (n = 3–8) can be positive, negative, or neutral [6]; the hydrogen atom(s) can be replaced by isovalent functional groups (-X, X = F, Cl, Br, I, OH, CH3, etc.) [7,8,9,10]; the -C5H5− ring can be replaced by other heterocycles; the Fe cation can be replaced by other transition metals with different spin states [11,12,13], and metallocene molecules can be doubled [14] or coordinated by other molecules or ligands. Here, we explored building polar crystals through the crystallization of polar molecules. Ferrocene itself has no (when the rings thermally rotate) or only a small possible dipole moment (as a static molecule) since it is a symmetric molecule. It crystallizes in a monoclinic structure (P21/a) at room temperature [15], which transforms metastably to triclinic at 164 K [16], or stably to orthorhombic [17]. The centrosymmetric staggered molecular structure within the P21/a phase transfers into two slightly distinct types of non-centrosymmetric, almost eclipsed, ferrocene molecules within the triclinic phase [16,18,19]. Density functional theory (DFT) calculations show that the ferrocene molecule with a staggered configuration has a very weak dipole moment (0.006 D), and the two types of molecules within the triclinic crystalline ferrocene increase to 0.08 D and 0.04 D, respectively [20]. Recently, using first-principle calculations, Wu et al. [21] predicted that crystallized M(C6H6) and M(C5H5) nanowires become ferroelectric or even multiferroic, depending on the decorating functional groups. Nanowires of M(Cp)∞ (M = Ti, Cr, Fe) are predicted to be ferromagnetic half-metals, and [M(FeCp2)]∞ (M = Sc, Ti, V, Mn) are ferromagnetic semiconductors [22,23].

2. Assessing van der Waals (vdW) Functionals

Density functional theory gives reasonable geometries for metallocene molecules [24,25,26]. Xu et al. [27] compared several gradient-corrected functionals, as well as the local density approximation (LDA), and found that the B3LYP functional gives the best agreement between the experimental and theoretical molecular geometry. However, for molecular crystals, accurate van der Waals forces must be included. We computed the equilibrium crystalline geometry starting from the experimental room temperature structure, the magnetic ground states, and magnetic moments of ferrocene, nickelocene, cobaltocene, and vanadocene, which are experimentally available, to assess the performance of the vdW functionals DFT-D, rVV10, and vdW-DF2 [28,29,30].
We found that the vdW-DF2 functional is most accurate of the three (Table 1), and provides quite good agreement with experiments, within ~3% of the observed cell volumes. In contrast, DFT-D and rVV10 methods significantly underestimate cell volumes. Therefore, we employed the vdW-DF2 functional for the subsequent calculations.

3. Dipole Moments for Isolated Molecules

Next, we considere maximizing the dipole moments for isolated molecules. We relaxed the isolated molecules in periodic cells of various sizes, to compare directly with our plane wave pseudopotential computations using the same basis set for crystals [32]. We computed the dipole moments using maximally localized Wannier functions [33]. We found that the dipole moment of an isolated ferrocene molecule with an eclipsed configuration is 0.08 D, with the direction being parallel to the two Cp rings (Figure 1, Table 2), consistent with simulation results of J. M. B.-García, et al. [20]. The molecular dipole moments can be increased by ligand decorations or by replacing the Cp with heterocyclic rings. The replacing of one C atom by Si increases the dipole moment substantially to 0.25 D for Fe(C9H10Si), with the positive polarization towards Si and the unsubstituted Cp ring. This is consistent with the weaker electronegativity of Si relative to C. When we replaced one H atom by F, we found more electrons on the -C5H4F− ring than on the Cp side, consistent with the stronger electronegativity of F with respect to H. Therefore, Fe(C10H9F) has a large dipole moment (1.70 D), pointing away from the unsubstituted Cp ring. Experimentally, metallocenes with 1–5 halogens per ring have been synthesized [10,34,35,36,37,38,39]. When we replaced F with heavier halogens, we found that the dipole moment increased to 1.77, 1.82, and 1.81 D in Fe(C10H9Cl), Fe(C10H9Br), and Fe(C10H9I), respectively (Table 3). We propose that the increase in dipole moment was simply due to the larger bond lengths of the positive and negative charges, opposite to what one would expect from electronegativity, which decreases with increasing atomic number.
Next, we considered replacing a -C5H5− ring by a charge neutral furan (-C4H4O). This sandwich structure was synthesized experimentally for the Ru-metallocene [40]. The iron atom donated one less electron to the neutral -C4H4O ring, giving rise to Fe+ (d7), and 19 electrons for the molecule. We obtained, for furanferrocene Fe(C9H9O), a magnetic moment of 1 μB (indicating a low-spin iron), and a large dipole moment (1.32 D). When we replaced a Cp ring with pyrrole (-C4H4N−), we found a larger dipole moment (2.08 D). The dipole moment of Fe(C9H9N) points from the N atom to the other side of the heterocycle and to the unsubstituted Cp ring. The electron-donating capacity of B, P, As, Sb, and Bi are weaker than that of N; however, the dipole moment decreases to 0.87, 1.71, 1.42, 0.70 and 0.52 D in Fe(C9H9B), Fe(C9H9P), Fe(C9H9As), Fe(C9H9Sb), and Fe(C9H9Bi), respectively (Table 4).
The magnitude of the dipole moment of azametallocene molecules shows a strong dependence on the atomic number of the cations (Table 2). Similar to M2(C10H10)2, the cations in azametallocene molecules donate two electrons (one to -C5H5− and the other to -C4H4N−) to form a stable aromatic configuration having 4k + 2 (k is an integer) electrons [41]. We found for vanadium in V(C9H9N) that the most stable 3d electronic configuration is the one with three parallel spins. Iron in Fe(C9H9N) is in a low spin state with zero magnetic moment. These results agree well with TM2(C10H10)2 (Table 1), as the substitution of -C5H5− by -C4H4N− does not change the number of electrons on the cations.
The relative stabilities and dipole moments of isolated metallocene molecules with two ring rotations show low rotation barriers (Erotation = Estaggered − Eeclipsed), which are lower than the rotation barrier within the crystalline ferrocene (4.4 ± 0.5 kJ/mol), [42] and comparable to ~kBT at room temperature (~2.6 kJ/mol, Table 2, Table 3 and Table 4). Therefore, the rotation of the Cp rings is unhindered at room temperature. The negative values of Erotation for Fe(C9H9O) and Ni(C9H9N) suggest that the staggered configuration, rather than the eclipsed configuration, is the ground state. Furthermore, we found that the dipole moments are dependent on the relative rotation of the rings. Upon rotation from the ground state to the transition state, the dipole moments vary by 10%, suggesting that ring rotation can be activated by applying time varying electric fields.

4. Crystalline Polar Metallocenes

Returning to crystalline properties, we considere polar metallocenes composed of molecules that display promising dipole moments. Azaferrocene Fe2(C9H9N)2 is isomorphous with ferrocene and crystallizes at room temperature in the monoclinic structure [43,44]. Thus, we first constructed the initial monoclinic azaferrocene using the room temperature [15] and low temperature [45] monoclinic ferrocene as prototype structures. They were found to be polar at their ground state, as shown in Figure 2a (M1) and b (M2), respectively. The net polarization depends on the magnitude as well as the arrangement of the dipole moments of the two constituent molecules in the unit cell. Therefore, we further consider two extreme cases, namely, the two constituent molecules antiparallel (AFE) and almost parallel (FE) to each other, relaxed from the triclinic ferrocene structure [10]. Within the AFE configuration, the dipole moments of the two molecules cancel each other, giving rise to a zero polarization (Figure 2c). The almost parallel arrangement of the dipole moments within the FE configuration yields the largest polarization among various structures investigated (Figure 2d).
The energy differences of the AFE configuration with respect to the FE configuration (ΔEpolar = EFE − EAFE, meV/per cell) are also shown in Table 5 to evaluate their relative stabilities. Mn2(C9H9N)2 has a negative ΔEpolar, which is highly desirable as the FE configuration is the ground state. All the other systems investigated have an AFE configuration. To investigate the stability of the crystals, we calculated the binding energy by Ebind = Emolecule − Eground/2, where Emolecule is the energy of the isolated molecules at their ground state, and Eground is the total energy per unit cell of triclinic polar metallocenes with the FE configuration at their magnetic ground state. The binding energies are large, ranging from 0.57 to 0.91 eV, and the relaxed cell volumes of azametallocenes decrease as the computed binding energies increase. The low energy differences between the AFE and FE configuration together with the large binding energies, suggest that such systems with ferroelectric ordering can be fabricated by applying an external bias field during crystallization and could be stable over a wide temperature range.
We further investigated the magnetic and electronic properties of polar metallocenes with the FE configuration. We computed the energy differences between the antiferromagnetic and ferromagnetic states (ΔEMS = EAFM − EFM) to estimate the exchange splitting, as shown in Table 6. Like ferrocene, Fe2(C10H9F)2 and Fe2(C9H9N)2 are non-magnetic. Fe2(C9H9O)2 has an AFM ground state. Similar to V2(C10H10)2 and Ni2(C10H10)2 (Table 1), we found that V2(C9H9N)2, Cr2(C9H9N)2, and Ni2(C9H9N)2 have AFM ground states, which are more stable than the FM states by 6.19, 1.60, and 4.93 meV per unit cell, respectively. Unlike Co2(C10H10)2, Co2(C9H9N)2 has a FM ground state, and has the largest ΔEMS (47.57 meV) among different azametallocenes investigated. Similarly, Mn2(C9H9N)2 has a FM ground state with a quite small energy difference (0.47 meV). Therefore, Mn2(C9H9N)2 is the most promising as ferroic material, since FM and FE ordering can be realized simultaneously.
With an FE configuration, the metallocene crystals formed from molecules with large dipole moments show large polarizations. The large polarizations observed for azametallocenes (e.g., >5 μC/cm2, Table 6) with respect to Fe2(C10H9F)2 and Fe2(C9H9O)2, are due to the large dipole moments of the component molecules. Similar to the phenomena observed for water [46,47], the polarization for the crystalline phase is about 30–40% higher than the values obtained based on the dipole moments for isolated molecules. This is due to the polarization of the molecules by the local self-consistent electrical field, and indicates high polarizabilities. Although the polarization of azametallocenes is smaller than inorganic perovskites (e.g., ~30 μC/cm2 for BaTiO3 [48]), a large piezoelectric response is possible, since metallocenes rae mechanically softer.

5. Materials and Methods

All the calculations were performed using the pseudopotential plane-wave method implemented in the Quantum Espresso code [32]. Three different van der Waals functionals, namely, DFT-D, rVV10, and vdW-DF2, were incorporated to describe the van der Waals interactions [49,50,51,52,53,54]. The van der Waals interaction is part of the correlation energy, and can be applied to the different exchange functionals. In particular, in Quantum Espresso, the non-local functionals vdW-DF2 and rVV10 are used together with the revised Perdew-Wang functional (rPW86) [55,56], and the DFT-D correction was applied to the generalized gradient approximation of Perdew, Burke, and Ernzerhof (PBE-GGA) [57].
The lattice constants as well as atomic positions are optimized until the forces were smaller than 0.01 eV/Å. The kinetic energy cutoffs for wavefunctions and charge density were set to 544.22 (40) and 2721.13 (200) eV(Ry), respectively. For crystalline metallocenes, the Brillouin zone was sampled by the Monkhorst-Pack scheme (4 × 4 × 4 k points for structure optimization, 6 × 6 × 6 k points for the subsequent self-consistent and property calculations). The ‘isolated molecules’ were simulated by putting the molecules into fixed unit cells with different sizes. We found unit cells with a size of 24 × 24 × 24 Å3 were large enough to well describe the internal stress and dipole moment of the isolated molecules. Maximally localized Wannier functions were generated by the Wannier90 code [33] to determine the negative charge center, based on which the dipole moments of the isolated molecules were calculated. GBRV pseudopotentials were employed for all the calculations [58].

6. Conclusions

We found polar metallocenes to be a worthwhile group to study as new active materials. Only recently has first-principles study of these systems become possible with new local van der Waals functionals, which we found to be indispensable. The vdW-DF2 functional gives accurate predictions for the structures and magnetic ground states of metallocene molecules and crystals. The dipole moment of isolated ferrocene molecules could be effectively tuned by ligand decorations or by replacing the Cp with heterocyclic rings. The magnetic moment of the ferrocene molecules could also be tuned by changing the charge of the cyclic conjugated ligands and/or by substitution of iron with other cations that had a different spin state. Our results suggest that crystalline polar metallocenes are worthy of study as potential ferroelectrics and multiferroics.

Author Contributions

R.E.C. conceived the idea; H.Z., B.Y.Y., and R.E.C. performed the simulations; H.Z. and R.E.C. designed the figures and analyzed the results; H.Z. drafted the manuscript; R.E.C. provided primary comments and editing; all authors approved the final manuscript.

Funding

This work is partly supported by the U. S. Office of Naval Research Grants No. N00014-12-1-1038, N00014-14-1-0516, and N00014-17-1-2768, the Carnegie Institution for Science, and the European Research Council Advanced Grant ToMCaT.

Acknowledgments

We thank Karlheinz Sünkel for helpful discussions. The authors gratefully acknowledge the Gauss Center for Supercomputing e.V. (www.gauss-centre.eu) for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (LRZ, www.lrz.de).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nicolaev, G.A.; Nicolaev, A. Transport in ferrocene single molecules for terahertz applications. Phys. Chem. Chem. Phys. 2014, 16, 5. [Google Scholar]
  2. Cunha, T.F.D.; Calderini, D.; Skouteris, D. Analysis of Partition Functions for Metallocenes: Ferrocene, Ruthenocene, and Osmocene. J. Phys. Chem. A. 2016, 120, 5282–5287. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, P.; Jiang, X.; Hu, J.; Huang, X.; Zhao, J.; Ahuja, R. Prediction of huge magnetic anisotropies in 5d transition metallocenes. J. Phys. Condens. Matter. 2017, 29, 435802. [Google Scholar] [CrossRef] [PubMed]
  4. Fischer, O.; Pfab, W. Cyclopentadien-Metallkomplexe, ein neuer Typ metallorganischer Verbindungen. Z. Nat. B 1952, 7, 3. [Google Scholar] [CrossRef]
  5. Wilkinson, G.R.M.; Whiting, M.C. The Structure of Iron Bis-Cyclopentadienyl. J. Am. Chem. Soc. 1952, 74, 2125–2126. [Google Scholar] [CrossRef]
  6. Frunzke, J.; Lein, M.; Frenking, G. Structures, Metal-Ligand Bond Strength, and Bonding Analysis of Ferrocene Derivatives with Group-15 Heteroligands Fe(η5-E5)2 and FeCp(η5-E5) (E=N, P, As, Sb). A Theoretical Study. Organometallics 2002, 21, 3351–3359. [Google Scholar] [CrossRef]
  7. Dombrowski, K.E.; Baldwin, W.; Sheats, J.E. Metallocenes in biochemistry, microbiology and medicine. J. Organomet. Chem. 1986, 302, 26. [Google Scholar] [CrossRef]
  8. Dunitz, J.D.; Orgel, L.E. Electronic Structure of Metal biscyclopentadienyls. J. Chem. Phys. 1955, 23, 954–958. [Google Scholar] [CrossRef]
  9. Vargas-Caamal, A.; Pan, S.; Ortiz-Chi, F.; Cabellos, J.L.; Boto, R.A.; Contreras-Garcia, J.; Restrepo, A.; Chattaraj, P.K.; Merino, G. How strong are the metallocene–metallocene interactions? Cases of ferrocene, ruthenocene, and osmocene. Phys. Chem. Chem. Phys. 2016, 18, 550–556. [Google Scholar] [CrossRef]
  10. Sünkel, K.; Weigand, S.; Hoffmann, A.; Blomeyer, S.; Reuter, C.G.; Vishnevskiy, Y.V.; Mitzel, N.W. Synthesis and Characterization of 1,2,3,4,5-Pentafluoroferrocene. J. Am. Chem. Soc. 2015, 137, 126–129. [Google Scholar] [CrossRef]
  11. Sohn, Y.S.; Hendrickson, D.N.; Gray, H.B. Electronic Structure of Metallocenes. J. Am. Chem. Soc. 1971, 93, 3603–3612. [Google Scholar]
  12. Haaland, A. Molecular Structure and Bonding in the 3d Metallocenes. Acc. Chem. Res. 1979, 12, 415–422. [Google Scholar] [CrossRef]
  13. Little, W.F. Metallocenes. Surv. Prog. Chem. 1963, 1, 78. [Google Scholar]
  14. Wu, X.; Zeng, X. Double Metallocene Nanowires. J. Am. Chem. Soc. 2009, 131, 14246–14248. [Google Scholar] [CrossRef] [PubMed]
  15. Seiler, P.; Dunitz, J.D. A New Interpretation of the Disordered Crystal Structure of Ferrocene. Acta Cryst. 1979, 35, 1068–1074. [Google Scholar] [CrossRef]
  16. Dunitz, J.D. Phase Changes and Chemical Reactions in Molecular Crystals. Acta Cryst. B 1995, 51, 619–631. [Google Scholar] [CrossRef]
  17. Seiler, P.; Dunitz, J.D. Low-temperature Crystallization of Orthorhombic Ferrocene: Structure Analysis at 98 K. Acta Cryst. Sect. B-Struct. Sci. 1982, 38, 1741–1745. [Google Scholar] [CrossRef]
  18. Seiler, P.; Dunitz, J.D. The Structure of Triclinic Ferrocene at 101, 123 and 148 K. Acta Cryst. 1979, B35, 2020–2032. [Google Scholar] [CrossRef]
  19. Edwards, J.W.; Kington, G.L.; Mason, R. The thermodynamic properties of ferrocene. Trans. Faraday Soc. 1960, 56, 660–667. [Google Scholar] [CrossRef]
  20. Bermúdez-García, J.M.; Yáñez-Vilar, S.; Castro-García, S.; Señarís-Rodríguez, M.A.; Sánchez-Andújar, M. New properties in old systems: Cooperative electric order in ferrocene and ammonia-borane. RSC Adv. 2015, 5, 83818–83824. [Google Scholar] [CrossRef]
  21. Wu, M.; Burton, J.D.; Tsymbal, E.Y.; Zeng, X.C.; Jena, P. Multiferroic Materials Based on Organic Transition-Metal Molecular Nanowires. J. Am. Chem. Soc. 2012, 134, 14423–14429. [Google Scholar] [CrossRef] [PubMed]
  22. Shen, L.; Yang, S.-W.; Ng, M.-F.; Ligatchev, V.; Zhou, L.; Feng, Y. Charge-Transfer-Based Mechanism for Half-Metallicity and Ferromagnetism in One-Dimensional Organometallic Sandwich Molecular Wires. J. Am. Chem. Soc. 2008, 130, 13956–13960. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, X.; Wang, J.; Gao, Y.; Zeng, X. Ab Initio Study of Structural and Magnetic Properties of TMn(ferrocene)n+1 (TM = Sc, Ti, V, Mn) Sandwich Clusters and Nanowires (n = ∞). ACS Nano 2009, 3, 537–545. [Google Scholar] [CrossRef]
  24. Fey, N. Organometallic molecular modelling-the computational chemistry of metallocenes: A review. J. Chem. Technol. Biotechnol. 1999, 74, 852–862. [Google Scholar] [CrossRef]
  25. Coriani, S.; Haaland, A.; Helgaker, T.; Jørgensen, P. The Equilibrium Structure of Ferrocene. ChemPhysChem 2006, 7, 245–249. [Google Scholar] [CrossRef] [PubMed]
  26. Aliabad, H.A.R.; Chahkandi, M. Comprehensive SPHYB and B3LYP-DFT Studies of Two Types of Ferrocene. Z. Anorg. Allg. Chem 2017, 643, 420–431. [Google Scholar] [CrossRef]
  27. Xu, Z.-F.; Xie, Y.; Feng, W.-L.; Schaefer, H.F. Systematic Investigation of Electronic and Molecular Structures for the First Transition Metal Series Metallocenes M(C5H5)2 (M=V, Cr, Mn, Fe, Co, and Ni). J. Phys. Chem. A 2003, 107, 2716–2729. [Google Scholar] [CrossRef]
  28. Seller, P.; Duntitz, J.D. The Structure of Nickelocene at Room Temperature and at 101 K. Acta Cryst. 1980, 36, 2255–2260. [Google Scholar] [CrossRef]
  29. Yu, M.; Antipin, R.B. Redetermination of the Cobaltocene Crystal Structure at 100 K and 297 K: Comparison with Ferrocene and Nickelocene. Struct. Chem. 1993, 4, 91–101. [Google Scholar]
  30. Yu, M.; Antipin, R.B. Structure of Vanadocene in the Temperature Interval 108–357 K and the Nature of the Ring Disorder. Acta Cryst. 1996, 52, 314–322. [Google Scholar]
  31. Robertson, R.E.; Mcconnell, H.M. The magnetic resonance properties of some sandwich compounds. J. Phys. Chem. 1960, 64, 70–77. [Google Scholar] [CrossRef]
  32. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 2009, 21, 395502. [Google Scholar] [CrossRef] [PubMed]
  33. Marzari, N.; Mostofi, A.A.; Yates, J.R.; Souza, I.; Vanderbilt, D. Maximally localized Wannier functions: Theory and applications. Rev. Mod. Phys. 2012, 84, 1419–1475. [Google Scholar] [CrossRef]
  34. Rosenblum, M.; Fish, R.W. A convenient synthesis of some haloferrocenes. J. Org. Chem. 1965, 30, 1253–1254. [Google Scholar]
  35. Bernhartzeder, S.; Sunkel, K. Coordination chemistry of perhalogenated cyclopentadienes and alkynes. Part 30. New high-yield syntheses of monochloroferrocene and 1,2,3,4,5-pentachloroferrocene. Molecular structures of 1,2-dichloroferrocene and 1,2,3-trichloroferrocene. J. Organomet. Chem. 2012, 716, 146–149. [Google Scholar] [CrossRef]
  36. Sunkel, K.; Bernhartzeder, S. Coordination chemistry of perhalogenated cyclopentadienes and alkynes. XXVIII [1] new high-yield synthesis of monobromoferrocene and simplified procedure for the synthesis of pentabromoferrocene. Molecular structures of 1,2,3-tribromoferrocene and 1,2,3,4,5-pentabromoferrocene. J. Organomet. Chem. 2011, 696, 1536–1540. [Google Scholar]
  37. Sünkel, K. Coordination chemistry of pentahalocyclopentadienyls. Chem. Ber. -Recl. 1997, 130, 1721–1730. [Google Scholar]
  38. Hnetinka, C.A.; Hunter, A.D.; Zeller, M.; Lesley, M.J.G. 1,1′-Dibromoferrocene. Acta Cryst. E 2004, 60, m1806–m1807. [Google Scholar] [CrossRef]
  39. Romanov, A.S.; Mulroy, J.M.; Khrustalev, V.N.; Antipin, M.Y.; Timofeeva, T.V. Monohalogenated ferrocenes C5H5FeC5H4X (X = Cl, Br and I) and a second polymorph of C5H5FeC5H4I. Acta Cryst. C 2009, 65, m426–m430. [Google Scholar] [CrossRef]
  40. Chaudret, B.; Jalon, F.A. Facile preparation of π-arene complexes of ruthenium [(η5-C5Me5)Ru(Arene)]X including a π-pyridine and the first π-furan derivatives. J. Chem. Soc. Chem. Commun. 1988, 711–713. [Google Scholar] [CrossRef]
  41. Da, H.; Jin, H.M.; Lim, K.H.; Ligatchev, V.; Ng, M.-F.; Khoo, K.H.; Yang, S.-W. Half metal semiconductor reversible switch in bimetallic sandwich molecular wire via redox reactions. Nano Energy 2012, 1, 297–302. [Google Scholar] [CrossRef]
  42. Gardner, A.B.; Howard, J.; Waddington, T.C. The dynamics of ring rotation in ferrocene, nickelocene and ruthenocene by incoherent quasi-elastic neutron scattering. Chem. Phys. 1981, 57, 453–460. [Google Scholar] [CrossRef]
  43. Brock, C.P.; Fu, Y. Rigid-Body Disorder Models for the High-Temperature Phase of Ferrocene. Acta Cryst. 1997, 928–938. [Google Scholar] [CrossRef]
  44. Clec’h, G.; Calvarin, G. Etude Structurale Des Phases Desordonnees De Co(C5H5)2, Mg(C5H5)2, Et Transitions De Phases De (C4H4)N)Fe(C5H5). Mol. Cryst. Liq. Cryst. 2011, 128, 305–320. [Google Scholar] [CrossRef]
  45. Chhor, K.; Pommier, C. Low temperature thermodynamic study of the stable and metastable phases of (C4H4N)Fe(C5H5). J. Chem. Thermodyn. 1984, 16, 503–517. [Google Scholar] [CrossRef]
  46. Gregory, J.K.; Clary, D.C.; Liu, K.; Brown, M.G.; Saykally, R.J. The Water Dipole Moment in Water Clusters. Science. 1997, 275, 814–817. [Google Scholar] [CrossRef] [PubMed]
  47. Silvestrelli, P.L.; Parrinello, M. Water Molecule Dipole in the Gas and in the Liquid Phase. Phys. Rev. Lett. 1999, 82, 3308–3311. [Google Scholar] [CrossRef]
  48. Shieh, J.; Yeh, J.H.; Shu, Y.C.; Yen, J.H. Hysteresis behaviors of barium titanate single crystals based on the operation of multiple 90° switching systems. Mater. Sci. Eng. B 2009, 161, 50–54. [Google Scholar] [CrossRef]
  49. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef]
  50. Sabatini, R.; Gorni, T.; de Gironcoli, S. Nonlocal van der Waals density functional made simple and efficient. Phys. Rev. B. 2013, 87, 4. [Google Scholar] [CrossRef]
  51. Vydrov, O.A.; Van Voorhis, T. Nonlocal van der Waals density functional: The simpler the better. J. Chem. Phys. 2010, 133, 10. [Google Scholar] [CrossRef] [PubMed]
  52. Thonhauser, T.; Zuluaga, S.; Arter, C.A.; Berland, K.; Schroder, E.; Hyldgaard, P. Spin Signature of Nonlocal Correlation Binding in Metal-Organic Frameworks. Phys. Rev. Lett. 2015, 115, 136402. [Google Scholar] [CrossRef] [PubMed]
  53. Thonhauser, T.; Cooper, V.R.; Li, S.; Puzder, A.; Hyldgaard, P.; Langreth, D.C. Van der Waals density functional: Self-consistent potential and the nature of the van der Waals bond. Phys. Rev. B 2007, 76, 11. [Google Scholar] [CrossRef]
  54. Langreth, D.C.; Dion, M.; Rydberg, H.; Schröder, E.; Hyldgaard, P.; Lundqvist, B.I. Van der Waals density functional theory with applications. Int. J. Quantum Chem. 2005, 101, 599–610. [Google Scholar] [CrossRef]
  55. Murray, E.D.; Lee, K.; Langreth, D.C. Investigation of Exchange Energy Density Functional Accuracy for Interacting Molecules. J. Chem. Theory Comput. 2009, 5, 2754–2762. [Google Scholar] [CrossRef] [PubMed]
  56. Lee, K.; Murray, E.D.; Kong, L.; Lundquist, B.I.; Langreth, D.C. Higher-accuracy van der Waals density functional. Phys. Rev. B. 2010, 82, 081101(R). [Google Scholar] [CrossRef]
  57. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef] [PubMed]
  58. Garrity, K.F.; Bennett, J.W.; Rabe, K.M.; Vanderbilt, D. Pseudopotentials for high-throughput DFT calculations. Comput. Mater. Sci. 2014, 81, 446–452. [Google Scholar] [CrossRef]
Molecules 24 00486 g001 550
Figure 1. The structure and dipole moments of isolated Fe(C10H10), Fe(C9H10Si), Fe(C10H9F), Fe(C9H9O), Fe(C9H9N), V(C9H9N), Cr(C9H9N), Mn(C9H9N), Co(C9H9N), and Ni(C9H9N) molecules with an eclipsed starting configuration. The green arrows denote the direction of the dipole moment for each molecule. Key: C: brown; Fe: shallow brown; H: pink; O: light blue; N: red; V: cyan; Cr: dark blue; Mn: purple; Co: blue; Ni: silver.
Figure 1. The structure and dipole moments of isolated Fe(C10H10), Fe(C9H10Si), Fe(C10H9F), Fe(C9H9O), Fe(C9H9N), V(C9H9N), Cr(C9H9N), Mn(C9H9N), Co(C9H9N), and Ni(C9H9N) molecules with an eclipsed starting configuration. The green arrows denote the direction of the dipole moment for each molecule. Key: C: brown; Fe: shallow brown; H: pink; O: light blue; N: red; V: cyan; Cr: dark blue; Mn: purple; Co: blue; Ni: silver.
Molecules 24 00486 g001
Molecules 24 00486 g002 550
Figure 2. Optimized crystal structure of azaferrocene constructed using the (a) room temperature and (b) low temperature monoclinic ferrocenes as prototype structures; triclinic ferrocene with the two constituent azaferrocene molecules (c) antiparallel (AFE) and (d) almost parallel (FE) to each other. Azaferrocene constructed using the two monoclinic ferrocene structures has a polar ground state. The two structures are denoted as M1 and M2, respectively. Key: C: dark brown; Fe: brown; H: pink; N: red.
Figure 2. Optimized crystal structure of azaferrocene constructed using the (a) room temperature and (b) low temperature monoclinic ferrocenes as prototype structures; triclinic ferrocene with the two constituent azaferrocene molecules (c) antiparallel (AFE) and (d) almost parallel (FE) to each other. Azaferrocene constructed using the two monoclinic ferrocene structures has a polar ground state. The two structures are denoted as M1 and M2, respectively. Key: C: dark brown; Fe: brown; H: pink; N: red.
Molecules 24 00486 g002
Table
Table 1. Static properties for ordered crystalline metallocenes. Results for magnetic order (FM, AFM and NM denote ferromagnetic, antiferromagnetic and non-spin-polarized, respectively), energy difference (ΔEMS, eV/molecule) relative to magnetic ground state, magnetic moment (M, μB/molecule), and cell volume (Vcell, Å3), of ferrocene, nickelocene, cobaltocene, and vanadocene, calculated using DFT-D, rVV10, and vdW-DF2 functionals. For ferrocene, a starting configuration with AFM or FM converges to the non-magnetic state for all three functionals. Therefore, the magnetic moment for the FM state is fixed to 8.0 μB per unit cell for comparison. Note that Ni2C20H20 with an AFM starting configuration also converges to non-magnetic with the rVV10 functional.
Table 1. Static properties for ordered crystalline metallocenes. Results for magnetic order (FM, AFM and NM denote ferromagnetic, antiferromagnetic and non-spin-polarized, respectively), energy difference (ΔEMS, eV/molecule) relative to magnetic ground state, magnetic moment (M, μB/molecule), and cell volume (Vcell, Å3), of ferrocene, nickelocene, cobaltocene, and vanadocene, calculated using DFT-D, rVV10, and vdW-DF2 functionals. For ferrocene, a starting configuration with AFM or FM converges to the non-magnetic state for all three functionals. Therefore, the magnetic moment for the FM state is fixed to 8.0 μB per unit cell for comparison. Note that Ni2C20H20 with an AFM starting configuration also converges to non-magnetic with the rVV10 functional.
DFT-DrVV10vdW-DF2Exp
FMAFMNMFMAFMNMFMAFMNM
Fe2(C10H10)2ΔEMS2.222-0.001.871-0.001.55-0.00-
M4.05-0.004.05-0.004.05-0.00.00 a
Vcell366.43-344.37387.12-363.45416.37-391.85395.00 b
Ni2(C10H10)2ΔEMS0.00.0010.3100.00-0.3650.000.0370.166-
M2.182.150.002.18-0.002.172.150.002.89 a
Vcell358.10358.23353.87373.25-370.79398.08384.73412.09403.67 b
Co2(C10H10)2ΔEMS1.1810.000.2250.9480.001.4320.7550.000.468-
M3.111.240.003.111.250.003.101.250.001.73 a
Vcell366.09351.22350.59386.45371.93450.32421.09405.40462.36418.57 b
V2(C10H10)2ΔEMS0.0030.001.2000.0030.001.2520.0010.0001.104-
M3.113.110.003.093.180.003.163.080.003.78 a
Vcell370.57369.95356.99387.93387.48376.42418.79418.54418.97414.77 b
a Reference [31]; b References [15,28,29,30], for ferrocene, nickelocene, cobaltocene, and vanadocene, respectively.
Table
Table 2. Number of 3d electrons on the transition metal (TM) ions (N); magnetic moment (M, μB), dipole moment (μ, D) of isolated metallocene molecules; rotation barriers (Erotation = Estaggered − Eeclipsed, kJ/mol) and dipole moment differences (Δμ = μstaggered − μeclipsed, D) between the two different configurations.
Table 2. Number of 3d electrons on the transition metal (TM) ions (N); magnetic moment (M, μB), dipole moment (μ, D) of isolated metallocene molecules; rotation barriers (Erotation = Estaggered − Eeclipsed, kJ/mol) and dipole moment differences (Δμ = μstaggered − μeclipsed, D) between the two different configurations.
Fe(C10H10)Fe(C9SiH10)Fe(C10H9F)Fe(C9H9O)Fe(C9H9N)V(C9H9N)Cr(C9H9N)Mn(C9H9N)Co(C9H9N)Ni(C9H9N)
N6667634578
M0001032112
μ0.080.251.701.322.081.752.371.842.202.80
Erotation2.613.382.41−3.762.220.290.972.121.06−1.74
Δμ−0.010.040.040.07−0.250.31−0.400.260.02−0.03
Table
Table 3. Dipole moment (μ, D) of isolated metallocene molecules with an eclipsed configuration; activation energy (Eact = Etransition state − Eground state, kJ/mol) and dipole moment difference (Δμ = μtransition state − μground state, D) between the transition state and the ground state. The number of 3d electrons on Fe is 6, and the magnetic moment of each isolated molecule is 0.
Table 3. Dipole moment (μ, D) of isolated metallocene molecules with an eclipsed configuration; activation energy (Eact = Etransition state − Eground state, kJ/mol) and dipole moment difference (Δμ = μtransition state − μground state, D) between the transition state and the ground state. The number of 3d electrons on Fe is 6, and the magnetic moment of each isolated molecule is 0.
Fe(C10H9Cl)Fe(C10H9Br)Fe(C10H9I)
μ1.771.821.81
Eact2.192.162.01
Δμ0.01−0.21−0.02
Table
Table 4. Dipole moment (μ, D) of isolated metallocene molecules with an eclipsed configuration; activation energy (Eact = Etransition state − Eground state, kJ/mol) and dipole moment difference (Δμ = μtransition state − μground state, D) between the transition state and the ground state. The number of the 3d electrons rest on Fe is 6, and the magnetic moment of the isolated molecules is 0.
Table 4. Dipole moment (μ, D) of isolated metallocene molecules with an eclipsed configuration; activation energy (Eact = Etransition state − Eground state, kJ/mol) and dipole moment difference (Δμ = μtransition state − μground state, D) between the transition state and the ground state. The number of the 3d electrons rest on Fe is 6, and the magnetic moment of the isolated molecules is 0.
Fe(C9H9B)Fe(C9H9P)Fe(C9H9As)Fe(C9H9Sb)Fe(C9H9Bi)
μ0.871.711.420.700.52
Eact3.412.603.603.083.19
Δμ0.100.01−0.110.160.30
Table
Table 5. Cell volume (Vcell, Å3) at magnetic ground state, energy difference of the triclinic polar metallocenes with the FE configuration with respect to the AFE configuration (ΔEpolar = EFE − EAFE, meV per cell), binding energy (Ebind = Emolecule − Eground/2, eV per molecule).
Table 5. Cell volume (Vcell, Å3) at magnetic ground state, energy difference of the triclinic polar metallocenes with the FE configuration with respect to the AFE configuration (ΔEpolar = EFE − EAFE, meV per cell), binding energy (Ebind = Emolecule − Eground/2, eV per molecule).
Fe2(C10H9F)2Fe2(C9H9O)2V2(C9H9N)2Cr2(C9H9N)2Mn2(C9H9N)2Fe2(C9H9N)2Co2(C9H9N)2Ni2(C9H9N)2
Vcell404.58382.58404.87392.73386.79378.17385.17382.98
ΔEpolar80.9674.9612.225.99−0.8717.6744.786.53
Ebind0.890.760.900.850.760.910.870.57
Table
Table 6. Magnetic ground state (MGS), energy difference between AFM and FM states (ΔEMS = EAFM − EFM, meV/per cell), magnetic moment (Mcell, μB), polarization (P, μC/cm2) by Berry’s phase calculations, and polarization (P*, μC/cm2) calculated by summing the dipole moment of isolated molecules for triclinic polar metallocenes with the FE configuration.
Table 6. Magnetic ground state (MGS), energy difference between AFM and FM states (ΔEMS = EAFM − EFM, meV/per cell), magnetic moment (Mcell, μB), polarization (P, μC/cm2) by Berry’s phase calculations, and polarization (P*, μC/cm2) calculated by summing the dipole moment of isolated molecules for triclinic polar metallocenes with the FE configuration.
Fe2(C10H9F)2Fe2(C9H9O)2V2(C9H9N)2Cr2(C9H9N)2Mn2(C9H9N)2Fe2(C9H9N)2Co2(C9H9N)2Ni2(C9H9N)2
MGS/AFMAFMAFMFM/FMAFM
ΔEMS/−7.92−6.20−1.600.47/47.57−4.93
Mcell00002.002.00
P4.403.165.036.065.435.565.796.56
P*2.721.912.773.923.073.553.694.80

Share and Cite

MDPI and ACS Style

Zhang, H.; Yavorsky, B.Y.; Cohen, R.E. Polar Metallocenes. Molecules 2019, 24, 486. https://doi.org/10.3390/molecules24030486

AMA Style

Zhang H, Yavorsky BY, Cohen RE. Polar Metallocenes. Molecules. 2019; 24(3):486. https://doi.org/10.3390/molecules24030486

Chicago/Turabian Style

Zhang, Haiwu, B.Yu. Yavorsky, and R.E. Cohen. 2019. "Polar Metallocenes" Molecules 24, no. 3: 486. https://doi.org/10.3390/molecules24030486

APA Style

Zhang, H., Yavorsky, B. Y., & Cohen, R. E. (2019). Polar Metallocenes. Molecules, 24(3), 486. https://doi.org/10.3390/molecules24030486

Article Metrics

Back to TopTop