---
abstract: 'The properties of organometallic wires \[TM$_2$(Ant)\]$_{\infty}$ constructed with transitional metals (TM = Sc, Ti, V, Cr, Mn and Fe) and anthracene (Ant) are investigated by first-principles calculations. As the gap between HOMO (Highest Occupied Molecular Orbital) and LUMO (Lowest Unoccupied Molecular Orbital) of Ant is much smaller than that of benzene (Bz), much larger charge transfer (CT) occurs between TMs and Ant, which results in much more diverse magnetic properties in \[TM$_2$(Ant)\]$_{\infty}$ than in \[TM(Bz)\]$_{\infty}$. Particularly, \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ are found to be half-metallic ferromagnets. As a result of this and the better structural stability, compared with \[TM(Bz)\]$_{\infty}$, \[TM$_2$(Ant)\]$_{\infty}$ (like \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$) may be better candidates of spintronic devices. Furthermore, as the HOMO-LUMO gap of small pieces of graphene (SPG), such as pentacene and coronene, decreases with the increase of polycyclic number, the CT effects may also fit for the TM-SPG sandwich polymers which can also act as good spintronic materials.'
author:
- 'Xianlong Wang, Xiaohong Zheng and Zhi Zeng$\footnote{Correspondence author: zzeng@theory.issp.ac.cn}$'
title: '**Organometallic Wires Constructed from Transitional Metals and Anthracene: A Theoretical Study**'
---

**Introduction**
================

One-dimensional organometallic sandwich polymers (OSPs) have attracted wide interests in recent years due to their potential applications in spintronic devices and information storage [@ref0MiyajimaK; @ref0LyonJT; @ref0MiyajimaK1; @ref0WangJL; @ref0MaslyukVV; @ref0KoleiniM; @ref0XiangHJ; @ref0WengHM; @ref0WengHM1; @ref0NagaoS; @ref0ZhouLP; @ref0ShenL; @ref0KurikawaT; @ref0KatzTJ; @ref0WuXJ; @ref0WangL]. One of the extensively studied systems is TM$_n$Bz$_{n+1}$ which demonstrates versatile properties. The property diversity arises from two facts. On one hand, the properties of \[TM(Bz)\]$_{\infty}$ are highly dependent on the TM ions. For example, \[V(Bz)\]$_{\infty}$ and \[Mn(Bz)\]$_{\infty}$ are ferromagnetic half-metals, while \[Cr(Bz)\]$_{\infty}$ is a nonmagnetic insulator. On the other hand, the type of organic molecules can also significantly influence the properties of the one-dimensional OSPs. For example, no magnetic moment appears in \[Mn(Cp)\]$_{\infty}$ (Cp = cyclopentadienyl), while \[Cr(Cp)\]$_{\infty}$ has a ferromagnetic ground state.[@ref0NagaoS; @ref0ZhouLP; @ref0ShenL] It is believed that the ferromagnetism obsvered in such OSPs (like V$_n$Bz$_{n+1}$[@ref0MiyajimaK]) may qualify them as good materials for spin polarized transport.[@ref0WangJL; @ref0MaslyukVV; @ref0KoleiniM; @ref0XiangHJ]

Because of the abundant magnetic properties and the possibility as candidates for spintronic devices and information storage, recently, the study on the one-dimensional OSPs has been extended to the cases of TM-pentalene (Pn), TM-naphthalene (Np) and TM-Ant.[@ref0KurikawaT; @ref0KatzTJ; @ref0WuXJ; @ref0WangL] The sandwich like clusters of TM-Np and TM-Ant were firstly synthesized in 1999.[@ref0KurikawaT] Further theoretical investigations showed that \[V$_2$(Ant)\]$_{\infty}$ is a half-metallic ferromagnet. Since in \[TM$_2$(Ant)\]$_{\infty}$, the Ant are connected by two TM ions, as shown in Fig 1(a), compared with \[TM(Bz)\]$_{\infty}$ and \[TM(Cp)\]$_{\infty}$, it should have better stability and conductivity and it may be more suitable for spintronic devices and information storage. Note that, previous theoretical investigation pointed out a possible structural transition from *D$_{6h}$* to *D$_2$* symmetry appearing on V$_n$Bz$_{n+1}$.[@ref0WangJL] Consequently, systematic investigate on this kind of polymers is quite interesting and necessary. Therefore, in this work, the structural, electronic and magnetic properties of \[TM$_2$(Ant)\]$_{\infty}$ (TM=Sc, Ti, V, Cr and Mn) are investigated by first principles calculations. We find that magnetic properties in \[TM$_2$(Ant)\]$_{\infty}$ are much different from those in TM(Bz)$_{\infty}$ and the difference arises from the different HOMO-LUMO gap of the organic molecules and the subsquent different CT between the TMs and the organic molecules. Especially, we find that \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ are half metals with the 100% spin polarized density of states (DOS) around the Fermi level, which is of great importance for spin transport.

Compupational Details
=====================

Full-electron density functional theory (DFT) implemented in the Dmol$^3$ package [@ref0DellyB; @ref0DellyB1] is used in our calculations. The Perdew, Burke, and Ernzerhof (PBE) [@ref0PerdewJP] functional is used for exchange and correlation interaction in the generalized gradient approximation(GGA). The basis set is composed of double numerical basis and a polarized function (DNP) is chosen. The supercell is chosen as 20$\times$20$\times$*c* Å$^3$, where $c$ is the lattaice constant along the periodic direction and is varied for different TM cases \[see Table I\]. The corresponding Brillouin zone is sampled by a 1$\times$1$\times$50 Monkhorst and Pack grid.[@ref0MonkhorstHJ] All structures are fully relaxed.

In order to investigate the different magnetic configurations and Peierls transition in \[TM$_2$(Ant)\]$_{\infty}$, a supercell containing four TMs and two Ants are used. For each kind of \[TM$_2$(Ant)\]$_{\infty}$, as shown in Figs. 1(b)-(f), five magnetic configurations are considered. The combination of two characters *XY* (*X* and *Y* correspond to *N*(non-magnetic), *F*(ferromagnetic) and *A*(anti-ferromagnetic)) are adopted for representing the related magnetic configurations, with *X* and *Y* describing the intralayer and interlayer magnetic coupling of TM atoms, respectively. For example, as shown in Fig. 1(d), *FA* means the ferromagnetic coupling between the inner layer TM atoms and the anti-ferromagnetic coupling between two nearest neighbor layers.

**Results and Discussions**
===========================

After full relaxation, the calculated structural, electronic and magnetic properties of \[TM$_2$(Ant)\]$_{\infty}$ are summarized in Table I. Peierls transition does not occur in \[TM$_2$(Ant)\]$_{\infty}$, and no magnetic moment is observed in the cases of \[Sc$_2$(Ant)\]$_{\infty}$ and \[Mn$_2$(Ant)\]$_{\infty}$. In contrast, Ti and Fe are spin polarized, but the magnetic ground states of \[Ti$_2$(Ant)\]$_{\infty}$ and \[Fe$_2$(Ant)\]$_{\infty}$ are *FA* and *AF*, respectively. Hence, the total magnetic moments of \[Ti$_2$(Ant)\]$_{\infty}$ and \[Fe$_2$(Ant)\]$_{\infty}$ supercells are zero. In addition, our calculations predict that *FF* is the ground state of both \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$, and that the supercells of \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ respectively possess a magnetic moment of 2.00$\mu_{B}$ and 0.96$\mu_{B}$ per TM ion, with a finite negative magnetic moment formed on the Ant part. Further analysis shows that \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ are half-metallic ferromagnets which can be seen from the band structures of \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ shown in Figs. 2(a) and 2(b). The \[V$_2$(Ant)\]$_{\infty}$ total DOS and the partial DOS of V 3*d* are shown in Fig. 2(c). It is found that the 100% spin-polarized DOS around Fermi level mainly come from the V 3*d* orbitals, with small contribution from the C 2*d*. A gap of 0.76 eV and 1.25 eV for \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ respectively is observed in the spin majority channel. Thus, \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ may be good candidates for spintronic devices and information storage. Note that our relsults on \[V$_2$(Ant)\]$_{\infty}$ presented here agree well with former calculations.[@ref0WangL]

In both \[TM$_2$(Ant)\]$_{\infty}$ and \[TM(Bz)\]$_{\infty}$, TMs are located in a crystal field with *D$_{6h}$* symmetry. Under *D$_{6h}$* symmetry crystal field, as shown in Fig. 3(a), TMs *d* orbital splits into a a$_1$ (d$_{3z^2-r^2}$), 2-fold e$_2$ (d$_{xy}$ and d$_{x^2-y^2}$) and e$_1$ (d$_{xz}$ and d$_{yz}$) orbitals. In the meantime, both Ant and Bz containing 14 and 6 C atoms obey the Hükckel rule, which means that a stable aromatic molecule should have 4*m*+2 carbon atoms (where *m* is an integer). However, the magnetic behaviors of \[TM$_2$(Ant)\]$_{\infty}$ are totally different from those of \[TM(Bz)\]$_{\infty}$.[@ref0XiangHJ] For example, \[Mn(Bz)\]$_{\infty}$ is a ferromagnet with 1.0$\mu_{B}$ magnetic moment per Mn ions, and non-magnetic state is the ground state of \[Cr(Bz)\]$_{\infty}$. Whereas, this situation is reversed in \[TM$_2$(Ant)\]$_{\infty}$, since \[Mn$_2$(Ant)\]$_{\infty}$ do not show any magnetic moment while Cr ions of \[Cr$_2$(Ant)\]$_{\infty}$ are ferromagnetically spin-polarized with 1$\mu_{B}$ magnetic moment. Both \[V(Bz)\]$_{\infty}$ and \[V$_2$(Ant)\]$_{\infty}$ are a spin-polarized sandwich polymers, but the magnetic moments of V ions in \[V(Bz)\]$_{\infty}$ and \[V$_2$(Ant)\]$_{\infty}$ are approximately 1.0$\mu_{B}$[@ref0XiangHJ] and 2.0$\mu_{B}$, respectively. An interesting question arises from these results: Why is there so big difference in the magnetic properties of \[TM$_2$(Ant)\]$_{\infty}$ and \[TM(Bz)\]$_{\infty}$?

Further analysis shows that, compared with \[TM(Bz)\]$_{\infty}$, there is a larger CT between TMs and Ant. We find that the magnetic properties of \[TM$_2$(Ant)\]$_{\infty}$ are attributed to the CT effect which was once proved to be responsible for the magnetic behaviors of \[TM(Cp)\]$_{\infty}$.[@ref0ShenL] For example, in \[V(Bz)\]$_{\infty}$, V (3d$^3$4s$^2$) has five valence electrons occupy the *d* orbitals with two and three electrons respectively occupying the *a$_1$* and doubly degenerate *e$_2$* orbitals, because of two *s* electrons fill 3*d* shell[@ref0MaslyukVV]. Hence, one unpaired electron exists in *e$_2$* orbitals, and \[V(Bz)\]$_{\infty}$ shows a magnetic moment of 1.0$\mu_{B}$ on each V ion. However, in \[V$_2$(Ant)\]$_{\infty}$, one electron at V *d* orbital is transferred to Ant, and *e$_2$* has two unpaired electrons. Thus, V ions in \[V$_2$(Ant)\]$_{\infty}$ possess a 2.0$\mu_{B}$ magnetic moment. In order to illustrate this effect more clearly, calculations are performed, with the unit cell containing two TMs and one Ant and the initial magnetic coupling of two TMs arranged as ferromagnetic. For each \[TM$_2$(Ant)\]$_{\infty}$, the calculated total unit-cell magnetic moments per TM are shown in Fig. 3(b), in which the magnetic moments and schematic occupation of valence electrons obtained by CT model are also presented. We can find that our calculated results agree well with the CT model. Furthermore, \[TM$_2$(Ant)\]$_{\infty}$ and \[TM(Cp)\]$_{\infty}$ are quite similar in magnetic properties, since Cp containing 5 C atoms does not obey the Hückel rule, Cp tends to get one additional electron from TMs. Therefore, CT occurs in \[TM(Cp)\]$_{\infty}$.[@ref0ShenL] However, according to Hückel rule, Ant is a stable aromatic molecule. We want to inspect why there is so big CT in \[TM$_2$(Ant)\]$_{\infty}$. Generally, the key factors determining the CT in this kind of polymers are: (1) the average number of carbon atom (ANC) per TM; (2) the strength of frontier orbital hybridization. From \[TM(Bz)\]$_{\infty}$ to \[TM$_2$(Ant)\]$_{\infty}$, the ANC increases from six to seven. This small change of ANC should not induce large CT difference between \[TM(Bz)\]$_{\infty}$ and \[TM$_2$(Ant)\]$_{\infty}$. Then the large CT in \[TM$_2$(Ant)\]$_{\infty}$ might arise from the second factor.

Our calculations show that the HOMO-LUMO gap of Ant is 2.7 eV, which is smaller than that of Bz(5.5 eV). Comparing with Bz, the LUMO of Ant is much closer to the HOMO of TM ions. Thus, the frontier orbital hybridization between TMs and Ant is much stronger in \[TM$_2$(Ant)\]$_{\infty}$ than in \[TM(Bz)\]$_{\infty}$ and the CT appears on \[TM$_2$(Ant)\]$_{\infty}$. This is also illustrated clearly in Figs. 3(c) and 3(d), where the partial DOS of TM 3*d* and Ant C 2*p* are shown. Taking Fermi level as the reference, the HOMO of both Bz and Ant are localized around -5.5 eV. However, as the gap of Ant is smaller than that of Bz, the minimum energy of hybdridization of V 3*d* and C 2*p* electrons in \[V$_2$(Ant)\]$_{\infty}$ is -4 eV, which is lower than that of \[V(Bz)\]$_{\infty}$ (-3 eV). Thus, as shown in Fig 3(c) and (d), compared with \[V(Bz)\]$_{\infty}$, the minority electrons of V in \[V$_2$(Ant)\]$_{\infty}$ is transferred to Ant, which can also be found in \[V(Cp)\]$_{\infty}$.[@ref0ShenL] Consequently, from \[V(Bz)\]$_{\infty}$ to \[V$_2$(Ant)\]$_{\infty}$, the magnetic moment of V ion increases from 1$\mu_{B}$ to 2$\mu_{B}$. The schematic view of this opinion is also shown in Fig. 3, where the frontier orbitals (HOMO and LUMO) of Bz, V and Ant are approximately shown in the left, middle and right, respectively. The bigger arrow means the larger CT in \[TM(Bz)\]$_{\infty}$ than that in \[TM(Bz)\]$_{\infty}$.

Double exchange (DE) was attributed as one of the ferromagnetic origins of the in \[TM(Bz)\]$_{\infty}$ and \[TM$_2$(Pn)\]$_{\infty}$.[@ref0XiangHJ; @ref0WuXJ] However, in \[V(Bz)\]$_{\infty}$, the ferromagnetic state was believed to arise from the negative polarization of carbon atoms since magnetic moments anti-polarized by TMs appears on benzene.[@ref0WengHM] In addition, in \[V(Cp)\]$_{\infty}$, the ferromagnetism originates from the CT effects since there is CT between TMs and Cp. [@ref0ShenL] In the cases of \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ that are studied in this work, strong coupling between TMs and Ant within the energy range of -4 to 0 eV can be found \[see Figs. 2(c) and 3(d)\]. Furthermore, the wave functions at the $\Gamma$ point of \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ near the Fermi level also confirm the strong TMs-Ant coupling. Therefore, by comparing CT and negative polarization effect, it can conclude that the DE may be mainly responsible for the ferromagnetism.

Furthermore, Ant can be considered as a small piece of graphene (SPG) with multi-Bzs. SPG with more benzene rings, such as coronene and pentacene (Pe), can also be used for constructing this kind of sandwich like nanostructures and nanowires. As a matter of fact, the sandwich like TM-corroene polymers have already been synthesized.[@ref0PozniakBP; @ref0FosterNR] Such kind of polymers built with SPG and TMs may also be suitable materials for spintronic devices and information recording. Except for the properties of better stability and conductivity mentioned at the beginning, as shown in Figs. 4(a) and 4(b), TM-SPG polymers can be built into warped and branched strutures. In Fig. 4(a), the relaxed structure of Sc$_4$V$_4$(Ant)$_5$ cluster are shown. As the radius of Sc is larger than the that of V, Sc$_4$V$_4$(Ant)$_5$ cluster shows a warp structure. By co-doping and introducing TM defects, warped structure can be realized in this kind of sandwich polymers. Fig. 4(b) shows the schematic branched structure constructed with Pe, Bz and TMs. In fact, warped and branched structures already achieved in carbon nanotubes and nanowires[@ref0LiJ; @ref0MengGW] are central to the design of nanodevices and nanocircuits. While the different shapes and versatile properties exhibited by the sandwich like TM-SPG nanostructures and nanowires may qualify them as promising building blocks for spintronics devices and information storage. As the HOMO-LUMO gap of SPG decrease with the size increase of SPG, the CT effects and magnetic behaviors in the \[V$_2$(Ant)\]$_{\infty}$ may also exists in TM-SPG. Thus, further experimental and theoretical investigations are highly suggested in this interesting field.

**Conclusion**
==============

In conclusion, the properties of \[TM$_2$(Ant)\]$_{\infty}$ are investigated by first principles calculations. As the HOMO-LUMO gap of Ant is smaller than that of Bz, there is a large CT in \[TM$_2$(Ant)\]$_{\infty}$, which results in different magnetic behaviors in \[TM$_2$(Ant)\]$_{\infty}$ from TM(Bz)$_{\infty}$. Especially, \[V$_2$(Ant)\]$_{\infty}$ and \[Cr$_2$(Ant)\]$_{\infty}$ exhibit 100% spin polarized DOS around Fermi level which is of great importance for spin transport. Furthermore, we indicate that CT effect and magnetic behaviors shown in \[V$_2$(Ant)\]$_{\infty}$ may also exist in TM-SPG. Finally, due to the excellent stability, conductivity and unique shape (like branch and warp) that can be possibly realized, TM-SPG may be prospective materials for spintronics and information storage.

**ACKNOWLEDGEMENTS** {#acknowledgements .unnumbered}
====================

This work was supported by the National Science Foundation of China under Grant No 10774148, No 10904148 and No 10847162, the special Funds for Major State Basic Research Project of China(973), Knowledge Innovation Program of Chinese Academy of Sciences, and Director Grants of CASHIPS. Part of the calculations were performed in Center for Computational Science of CASHIPS and the Shanghai Supercomputer Center.

[00]{} K. Miyajima, A. Nakajima, S. Yabushita, M.B. Knickelbein and K. kaya, J. AM. Chem. Soc. [**126**]{}, 13202 (2004)

J.T. Lyon and L. Andrews, J. Phys. Chem. A [**109**]{}. 431 (2005)

K. Miyajima, S. Yabushita, M.B. Knickelbein and A. Nakajima, J. AM. Chem. Soc. [**129**]{}, 8473 (2007)

J.L. Wang, P.H. Acioli and J. Jellinek, J. AM. Chem. Soc. [**127**]{}, 2812 (2005)

V.V. Maslyuk, A. Bagrets, V. Meded, A. Arnold, F. Ever, M. Brandbyge, T. Bredow and I. Mertig, Phys. Rev. Lett.[**97**]{}, 097201 (2006).

M. Koleini, M. Paulsson and M. Brandbyge, Phys. Rev. Lett.[**98**]{}, 197202 (2007).

H.J. Xiang, J.L. Yang, J.G. Hou and Q.S. Zhu, J. AM. Chem. Soc. [**128**]{}, 2310 (2006)

H.M. Weng, T. Ozake and K. Terakura, J. Phys. Society of Japan [**77**]{}, 014301 (2008)

H. Weng, T. Ozaki and K. Terakura, Phys. Rev.B [**79**]{}, 235118 (2009). S. Nagao, A. Kata and A. Nakajima, J. AM. Chem. Soc. [**122**]{}, 4221 (2000)

L.P. Zhou, S.W. Yang, M.F. Ng, M.B. Sullivan, V.B.C. Tan and L. Shen, J. AM. Chem. Soc. [**130**]{}, 4023 (2008)

L. Shen, S.W. Yang, M.F. Ng, V. Ligatchev, L.P. Zhou and Y.P. Feng, J. AM. Chem. Soc. [**130**]{}, 13956 (2008) T. Kurikawa, H. Tkeda, M. Hirano, K. Judai, T. Arita, S. Nagao, A. Nakajima and K. Kaya, Organometallics. [**18**]{}, 1430 (1999)

T.J. Katz and N.J. Acton, J. AM. Chem. Soc. [**94**]{}, 3281 (1972)

X.J. Wu and X.C. Zeng, J. AM. Chem. Soc. [**131**]{}, 14246 (2009)

L. Wang, Z.X. Cai, J.Y. Wang, J. Lu, G.F Luo, L. Lai, J. Zhou, R. Qin, Z.X. Gao, D.P. Yu, G.P. Li, W.N. Mei, and S. Sanvito, Nano Letters [**8**]{}, 3640 (2008) B. Delly, J. Chem. Phys. [**92**]{}, 508 (1990)

B. Delly, J. Chem. Phys. [**113**]{}, 7756 (2000)

J.P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. [**77**]{}, 3865 (1996).

H.J. Monkhorst and J.D. Pack, Phys. Rev. B [**13**]{}, 5188 (1976).

B.P. Pozniak and R.C. Dunbar, J. AM. Chem. Soc. [**119**]{}, 10439 (1997)

N.R. Foster, G.A. Grieves, J.W. Buchanan, N.D. Flynn and M.A. Duncan, J. Phys. Chem. A [**104**]{}. 11055 (2000)

J. Li, C. Papadopoulos and J. Xu, nature [**402**]{}. 253 (1999)

G.W. Meng, Y.J. Jung, A. Cao, R. Vajtai and P.M. Ajayan, Proc. Natl. Acad. Sci. USA [**102**]{}. 7074 (2005)

   meV   *c*(Å)          GS          band gap (eV)   M($\mu_{B}$)       VED          
  ----- -------- ------------------ --------------- -------------- -------------- -- --
   Sc     7.89      *NN*(metal)          0.00            0.00       3d$^1$4s$^2$     
   Ti     7.37      *FA*(metal)          0.00            0.00       3d$^2$4s$^2$     
    V     7.01    *FF*(half metal)   0.76(direct)        2.00       3d$^3$4s$^2$     
   Cr     6.79    *FF*(half metal)   1.25(direct)        0.96       3d$^4$4s$^2$     
   Mn     6.64     *NN*(semicond)    0.66(direct)        0.00       3d$^5$4s$^2$     
   Fe     6.87      *AF*(metal)          0.00            0.00       3d$^6$4s$^2$     

  :  The lattice constant in the *c* axial direction, the electronic ground state (BS), the band gap for semiconductor and half metal cases, the total supercell magnetic moment (M) per TM and valence electron distribution (VED) of \[TM$_2$(Ant)\]$_{\infty}$ (TM= Sc, Ti, V, Cr, Mn and Fe).

1.0cm

3.0cm ![image](Fig1-xlw.eps){width="11.5cm" height="6.5cm"}

3.0cm ![image](Fig2-xlw.eps){width="8cm" height="11cm"}

3.0cm ![image](Fig3-xlw.eps){width="10cm" height="8cm"}

3.0cm ![image](Fig4-xlw.eps)
