R. N. Wang†
*a,
X. H. Zhenga,
J. Lana,
X. Q. Shia and
Z. Zengab
aKey Laboratory of Materials Physics, Institute of Solid State Physics, Chinese Academy of Sciences, Hefei 230031, P. R. China. E-mail: rnwang@qq.com
bDepartment of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
First published on 5th December 2013
We investigate the structural, magnetic properties and spin polarized transport of monatomic carbon chains covalently interconnected between zigzag graphene nanoribbon (ZGNR) leads with first-principles and non-equilibrium Green's function methods. It is revealed that an even numbered carbon chain is nonmagnetic when connected at the center of the armchair side of ZGNRs, but it becomes spin polarized and a net magnetic moment (about 0.1 μB) is developed when it moves to the edge. In contrast, an odd-numbered carbon chain is magnetic when it lies at the center of the armchair side and the magnetic moment decreases when the carbon chain moves to the edge. The total magnetic moment of the ZGNR/Cn/ZGNR junction depends on the magnetic alignment between two zigzag edges. Compared with that of ZGNR/C9/ZGNR with C9 chain anchoring at the center, the magnetic moment of the junction increases by 2.2 μB in the antiferromagnetic configuration and decreases by 1.5 μB in the ferromagnetic alignment when the carbon chains migrate from the center to one edge. Thus, the ZGNR/Cn/ZGNR junctions show rich magnetic varieties by modulating the carbon chain anchoring sites and the zigzag edge magnetic coupling. Further, spin-dependent transport calculations indicate that the ZGNR/C9/ZGNR junction behaves as a perfect spin-filter in spintronic circuits.
With the above-mentioned technique,9,10 it is very easy to create junctions in which monoatomic carbon chains are interconnected between two graphene nanoribbons (GNR) via sp2 bonding, herein referred to as GNR/Cn/GNR. Such all-carbon junctions have aroused much research interest.12–16 Experimentally, the conductance of such junctions could be ballistic, independent of structural imperfections such as distortion and absorbent except for oxygen impurities.17 Moreover, the equilibrium conductance parity oscillation differs from that of carbon chains directly interconnected to metallic electrodes in that the former is much more pronounced and could appear as the difference between on- and off-states, while the latter is not.18 In addition, in such a pure carbon-based junction, spin transport is feasible due to its long spin-coherence distance and time in the presence of weak spin–orbit coupling.19 Fürst et al. found a completely spin-polarized transmission due to the intrinsic spin-polarization of zigzag graphene nanoribbon (ZGNR).20
Yet in addition, it is well known that ZGNRs possess electron states localized near the edges with energies very close to Fermi energy (Ef),21 which may greatly affect the magnetic and transport behaviors of the sandwiched carbon chains. However, until now there have been few concerns22 on the effects of ZGNR's edge states on the magnetic and spin polarized transport properties of carbon chains. Therefore, in this work, we pay special attention to the magnetic and spin transport properties of the ZGNR/Cn/ZGNR junctions when the carbon chain migrates from the center to one zigzag edge along the armchair side of the cutting ZGNR. It is found that the magnetic moments of the carbon chain and whole junction and the spin-polarized transport of such junctions are dramatically altered as the carbon chain moves from the center to the edge, due to the change of spatial symmetry.23 The ZGNR/Cn/ZGNR complex systems can be modulated to show magnetic varieties by tailoring the magnetic coupling of ZGNR edges, the parity of carbon chains, and the anchoring sites to achieve various characteristics, such as semiconducting non-spin-polarized, metallic spin-polarized, and even semiconducting spin-polarized.
The paper is organized as follows. In Section II, we briefly describe the simulation model and computational method. The results and discussions are presented in detail in Section III. The paper ends with Section IV where a short summary is given.
The vacuum layers between two sheets are 20 Å so that each junction is laterally isolated from its periodic images. Firstly, we need to obtain the suitable electrode/electrode distance along the z direction. As a first step, we assume that the carbon atomic chains are in the ideal cumulene configuration and the initial distance between carbon atoms in the chain is set to 1.283 Å which is the C–C bond length of an infinite carbon chain.24 And then, by specifying several electrode/electrode distances varying, all the junctions are fully optimized by using SIESTA code25 until the absolute value of forces acting on each atom is less than 0.03 eV Å−1. In f
Further, the total energy as a function of z spacing is obtained. At last, a least-squares fit of the energy as the distance yields the optimal lattice constants, which corresponds to the minimum total energy.
And then, the spin-polarized transport calculations are carried out with the SMEAGOL program26 which combines density functional theory with non-equilibrium Green's function method. In ZGNR/Cn/ZGNR junctions, the valence electrons of carbon atoms are 2s22p2 and don't contain localized d or f electrons, so density functional theory is shown to be very suitable in predicting the ground states.27,28 However, the DFT + U method should be used when the interaction of graphene with 3d transition-metal adatoms is investigated. Because it is found that the electronic configuration and the adsorption geometries are very sensitive to the effects of local on-site Coulomb interactions U in the transition-metal d orbitals. All the calculations were performed within the generalized gradient approximation parameterized by Perdew, Burke, and Ernzerhof29 and with standard norm-conserving pseudopotentials.30 The real-space grid cutoff energy is 300 Ry and a Fermi–Dirac distribution function with an electronic temperature of 20 meV is used to populate the energy levels. The k-point sampling mesh for the leads is 1 × 1 × 20. Basis set quality, k-point sampling, kinetic energy cutoff, and electronic temperature have been tested with convergence studies on the structural, magnetic, and spin-polarized transport properties of ZGNR/Cn/ZGNR junctions.
![]() | ||
Fig. 2 The variation of carbon–carbon bond length with respect to C–C bond index for the carbon chain in all ZGNR/Cn/ZGNR systems. |
For cZGNR/Cn/cZGNR junctions, there is a reversal symmetry for the carbon–carbon bonds. For even n, the inversion symmetry point is on a bond and makes no restrictions on the character of the central bond. Due to a Peierls-type distortion, the single and triple bond alternation has a much lower energy. So the bond length for the seven C–C bonds (di) inside the chain of the cZGNR/C8/cZGNR junction oscillates in the range of 1.301 to 1.422 Å. However, if n is odd, the symmetrical point is at the center carbon atom and this carbon atom forms double bonds with its neighbors and all the other bonds thus have to be double bonds. Therefore, the bond length alternation for the cZGNR/C9/cZGNR junctions is less than that of the cZGNR/C8/cZGNR junctions and is in the range of 1.317 to 1.391 Å. But the d0s for the cZGNR/Cn/cGNR junctions are nearly the same and 1.498 and 1.478 Å for the cZGNR/C8/cZGNR and cZGNR/C9/cZGNR junctions, respectively. When the carbon chain moves to the edge site, the dis of eZGNR/C8/eZGNR junctions are similar to those of cZGNR/C8/cZGNR junctions. But for the eZGNR/C9/eZGNR junctions, dis change a lot in comparison with those of cZGNR/C9/cZGNR junctions due to the mutual effects and actions between the edge states and the states localized on the carbon chains, which will be discussed in detail in the following section. In general, the length of a single σ bond is 1.47 Å, whereas the double (σ + π) and triple (σ + 2π) bonds shrink to 1.380 and 1.210 Å, respectively. In conclusion, the bond length d0s are nearly close to the single σ bonds. The bonds in the inner Cn chains with odd n are close to a double (σ + π) bond. Whereas for the even n carbon chains, the bond lengths, dis, show single–triple bond alternation.
Extensive calculations were carried out for all the junctions using a ferromagnetic (FM, ZGNR/Cn/ZGNR↑↑) and an antiferromagnetic (AFM, ZGNR/Cn/ZGNR↑↓) coupling between the opposite edges of ZGNRs. It is shown that the total energies of the ZGNR/Cn/ZGNR↑↑ junctions are larger than those of the ZGNR/Cn/ZGNR↑↓ junctions, which is similar to the previous results.31 However, the total energies of ZGNR/C9/ZGNR↑↓ and ZGNR/C9/ZGNR↑↑ systems only differ at most by 7.76 meV when the carbon chain migrates from the center to the edge along the armchair side of ZGNR, suggesting that the alignment of edge magnetization in eZGNR/C9/eZGNR junctions can be easily influenced by thermal fluctuation.32
In the following, we investigate how the magnetic properties of all ZGNR/Cn/ZGNR systems change when carbon chains migrate from the center to one edge. It is found that the considered ZGNR/Cn/ZGNR junctions display an interesting magnetic behavior during the migration. The spin density (ρ↑ − ρ↓) is computed and presented in Fig. 3 and 4 when carbon chains anchored at the center and edge site, respectively. The magnetization of the ZGNR/Cn/ZGNR junctions is predicted to reside mainly on zigzag edges (shown in Fig. 3 and 4), in agreement with the predicted magnetism of zigzag edges in GNRs.33 However, the presence of a net spin density localized on the chain atoms depends exclusively on the parity of the chain and not the edge magnetic coupling of ZGNR. The same dependence of magnetization on the chain parity has also been reported by Ravagnan et al.34 for Cn chains with n from 4 to 12 with different kinds of termination.
This is because that the spin polarization of linear carbon chains is in fact an intrinsic property of the chains themselves, only depending on the resulting presence of unpaired electrons which is responsible for magnetism in all-carbon nanostructures.35–38 In a free-standing carbon chain with n atoms, 2s and 2pz valence electrons are hybridized, and then form two σ bonds with adjacent atoms making the carbon chain linear and stable. Then, the remaining two-fourth valence electrons (2px and 2py) make π bonds, which are 2-fold degenerate. Thus, there are (n − 1)/2 fully occupied π orbitals for odd n and (n/2 − 1) fully occupied π orbitals plus one half-filled π orbital for even n. Due to Hund's rule, a free odd-numbered carbon chain has 0 μB magnetic moment while a free even-numbered carbon chain has a 2 μB magnetic moment. When a carbon chain is covalently connected to ZGNR electrodes, it will accept two electrons from the electrodes. So, a new unoccupied π orbital is added for odd n or the partially occupied π orbital is filled for even n. This leads to the completely oppositive magnetic moment of carbon chain. Hence, for ZGNR/C8/ZGNR junctions, all the bonds inside the carbon chain are occupied with paired electrons and no magnetism is expected (shown in Fig. 3(a) and (b)). While, for ZGNR/C9/ZGNR junctions, a net spin density is observed (shown in Fig. 3(c) and (d) and Fig. 4(c) and (d)). As is well known that half-filled orbitals contribute to the electrical conduction while fully occupied orbitals do not. Obviously, the above-mentioned magnetic features will also be reflected in the transport properties of ZGNR/Cn/ZGNR junctions which will be discussed later.
In addition, comparing the spin density of eZGNR/C8/eZGNR junctions (Fig. 3(a) and (b)) with that of cZGNR/C8/cZGNR junctions (Fig. 4(a) and (b)), a small spin density is excited by edge states of ZGNRs when C8 chain moves from the center to one edge. The concrete magnetic moment per carbon atom inside the chain is analyzed by Mulliken population and shown in Fig. 5 and Table. 1. The magnetic moment of the C8 chain changes very little and increases from 0 to ∼0.1 μB. For ZGNR/C9/ZGNR junctions, an oscillating spin density is obviously observed on C9 chain as illustrated in Fig. 3(c) and (d) and 4(c) and (d), which is consistent with the analysis of Mulliken population (Fig. 5(c) and (d)). When the C9 chain anchors at the center site of the armchair side of ZGNRs, the magnetic moment on the chain is about 1.6 μB, as reported in Table 1. When the C9 chain migrates from the center site to the edge site, the magnetic moment on the C9 chain decreases from 1.653 to 1.207 μB for the AFM coupling and from 1.620 to 1.221 μB for the FM coupling, respectively. In conclusion, the magnetic moment of C8 chain increases from 0.0 to 0.1 μB, while that of the C9 chain decreases by 0.4 μB when the chain migrates from the center to the edge site along the armchair side of ZGNRs. It seems that the changes in the magnetic moment of carbon chains also depend on the parity during the migration.
![]() | ||
Fig. 5 Magnetic moment (μB) per atom as a function of carbon atom index (shown in Fig. 1) when the C8 (up panel) or C9 (down panel) chains are contacted by ZGNR↑↓ (left panel) or ZGNR↑↑ (right panel) from the center to the edge. |
Mchain (μB) | Mtotal (μB) | |
---|---|---|
cZGNR/C8/cZGNR↑↓ | 0.000 | 0.001 |
eZGNR/C8/eZGNR↑↓ | 0.093 | 0.001 |
cZGNR/C8/cZGNR↑↑ | 0.008 | 11.707 |
eZGNR/C8/eZGNR↑↑ | 0.087 | 11.640 |
cZGNR/C9/cZGNR↑↓ | 1.653 | 1.933 |
eZGNR/C9/eZGNR↑↓ | 1.207 | 4.163 |
cZGNR/C9/cZGNR↑↑ | 1.620 | 13.448 |
eZGNR/C9/eZGNR↑↑ | 1.221 | 11.927 |
The total magnetic moment on the whole ZGNR/C8/ZGNR junction is nearly zero for AFM coupling and 11.7 μB for FM coupling, respectively. In addition, the magnetic moment of ZGNR/C8/ZGNR junctions is hardly affected by the magnetic coupling of ZGNR. For ZGNR/C9/ZGNR junctions, the magnetic moments are more complex and all of them are spin-polarized regardless of the carbon chain anchoring sites and the magnetic coupling between the opposite edges of ZGNR. The total magnetic moment on the whole ZGNR/C9/ZGNR junctions has completely reversed charge and increases from 1.933 to 4.163 μB for the AFM coupling and decreases from 13.448 to 11.927 μB for the FM coupling, respectively.
For the FM alignment, one zigzag edge magnetic moment is suppressed and the other zigzag edge magnetic moment is hardly affected (shown in Fig. 4(d)), which makes the total magnetic moment decrease for the ZGNR/C9/ZGNR↑↑ junctions. This is consistent with the detailed magnetic moment per carbon atom as a function of the carbon index (shown in the Fig. 1) along the armchair side of ZGNRs in the Fig. 6. Nevertheless, the understanding of the magnetic moment changes for ZGNR/C9/ZGNR↑↓ junctions may be somewhat complicated. As is well known, graphene is composed of two different kinds of sublattice and the 1 and −1 (as denoted in Fig. 1) sites are the same. For odd-numbered carbon chains in contact with ZGNR leads, if they connect to the same sublattice, the sublattice symmetry is conserved and the electrodes have slight effects on the magnetic moments of the carbon chains. Otherwise, the symmetry is broken and the magnetic moment of carbon chains will decrease when the odd carbon chains are in contact with two different sublattices. So, the carbon chain in the ZGNR/C9/ZGNR junction has less than 2 μB magnetic moments. When the odd carbon chains migrate from the center to the edge, the magnetic moment will further decrease because edge states of ZGNR exponentially decay from the edge to the center. In contrast, for even carbon chains, no matter whether it is connected to the same sublattice or different sublattices, it always obtains a small finite magnetic moment which will increase with the carbon chain migrating from the center to the edge.
![]() | ||
Fig. 6 Magnetic moment (μB) per atom as a function of carbon atom index (shown in Fig. 1) when the C9 chain is contacted by ZGNR↑↓ ((a) and (b)) or ZGNR↑↑ ((c) and (d)) from the center ((a) and (c)) to the edge ((b) and (d)). |
Further, with n increasing, the odd–even conductivity oscillation magnitude becomes smaller and smaller and finally saturates to a constant, namely, 1 G0. This trend is easy to understand since the effect of the electrodes on the electronic structure of a long carbon chain is smaller than that on a short carbon chain.17 However, the above-mentioned conclusions are still applicable. The magnetic moment of ZGNR will also increase for AFM and decrease for FM whilst carbon chains move to the edge, because the symmetry is broken.
The transmission spectra of ZGNR/Cn/ZGNR junctions for the C8 and C9 chains are shown in Fig. 7 and 8, respectively. As mentioned above, the ZGNR/C8/ZGNR junctions are indeed semiconducting, regardless of the anchoring site and the magnetic ordering. From Fig. 7(a), we see that the transmission spectra of cZGNR/C8/cZGNR↑↓ junctions do not show spin-polarization. In contrast, when the carbon chain deviates from the center, the removal of degeneracy for β- and α-spin makes the junction a semiconducting spin-polarized structure (as shown in Fig. 7(b)). All ZGNR/C8/ZGNR↑↑ junctions also exhibit semiconducting spin-polarized features.
![]() | ||
Fig. 7 Transmission spectra of the cZGNR/C8/cZGNR systems: (a) cZGNR/C8/cZGNR↑↓; (b) eZGNR/C8/cZGNR↑↓; (c) cZGNR/C8/cZGNR↑↑; (d) eZGNR/C8/cZGNR↑↑. |
![]() | ||
Fig. 8 Transmission spectra of the cZGNR/C9/cZGNR systems: (a) cZGNR/C9/cZGNR↑↓; (b) eZGNR/C9/cZGNR↑↓; (c) cZGNR/C9/cZGNR↑↑; (d) eZGNR/C9/cZGNR↑↑. |
However, it is quite different for the ZGNR/C9/ZGNR junction. It is metallic in the FM state and semiconducting in the AFM state with an energy gap of 350 meV which is ∼295 meV in previous calculations.28 Outside the gap region, the conductance is no larger than 1 G0 (G0 denotes the quantum conductance and equals e2/h for each spin) in the presence of the scattering since the carbon atomic chain could provide at most one channel. The transmission spectra of ZGNR/C9/ZGNR↑↓ junctions show that it is always semiconducting and spin-polarized, no matter where the C9 chain anchors at the armchair profile of ZGNR↑↓ (as shown in panels (a) and (b) of Fig. 8). As far as we know, the FM alignment of ZGNR makes it conducting and the transmission is nonzero around Ef. However, the conductance also never exceeds 1 G0, as mentioned above. All the transmission spectra of ZGNR/C9/ZGNR↑↑ junctions are spin-polarized. The conductance around Ef becomes larger and larger when the carbon chain deviates from the edge, which shows that the contact between the carbon chain and ZGNRs is increasingly strong. As seen in panel (c) of Fig. 8, the obtained spin-resolved conductance at Ef is 0.019 G0 or 0.608 G0 for α- or β-spin channels, when the C9 chain locates at the center. When the C9 chain moves to the edge, the obtained spin-resolved conductance at Ef is 0.051 G0 and 0.003 G0 for α- and β-spin channels, respectively. To quantify the transmission spin polarization, the spin filter efficiency (SFE)39 in the linear response regime is then calculated from the transmission by its definition
Footnote |
† Current address: Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, P. R. China. |
This journal is © The Royal Society of Chemistry 2014 |