Chinese Physics Letters, 2017, Vol. 34, No. 9, Article code 097303 Fano Resonance Effect in CO-Adsorbed Zigzag Graphene Nanoribbons * Gao Wang(王皋)1, Meng-Qiu Long(龙孟秋)2**, Dan Zhang(张丹)2,3 Affiliations 1Advanced Research Center, Central South University, Changsha 410083 2Institute of Super Microstructure and Ultrafast Process, Central South University, Changsha 410083 3School of Science, Hunan University of Technology, Zhuzhou 412007 Received 5 June 2017 *Supported by the National Natural Science Foundation of China under Grant No 21673296, the Science and technology Plan of Hunan Province under Grant No 2015RS4002, and the Hunan Provincial Natural Science Foundation under Grant No 2017JJ3063.
**Corresponding author. Email: mqlong@csu.edu.cn
Citation Text: Wang G, Long M Q and Zhang D 2017 Chin. Phys. Lett. 34 097303 Abstract Quantum interference plays an important role in tuning the transport property of nano-devices. Using the non-equilibrium Green's Function method in combination with density functional theory, we investigate the influence to the transport property of a CO molecule adsorbed on one edge of a zigzag graphene nanoribbon device. Our results show that the CO molecule-adsorbed zigzag graphene nanoribbon devices can exhibit the Fano resonance phenomenon. Moreover, the distance between CO molecules and zigzag graphene nanoribbons is closely related to the energy sites of the Fano resonance. Our theoretical analyses indicate that the Fano resonance would be attributed to the interaction between CO molecules and the edge of the zigzag graphene nanoribbon device, which results in the localization of electrons and significantly changes the transmission spectrum. DOI:10.1088/0256-307X/34/9/097303 PACS:73.20.Fz, 73.22.Pr, 73.50.Bk © 2017 Chinese Physics Society Article Text The 20th century witnessed the development of science and technology, and researchers have entered into a Nano era. Many new nano-materials have been produced, especially for graphene. It is a single atomic layer of graphite. Novoselov et al.[1,2] have successfully synthesized it experimentally and demonstrated that it is stable at room temperature. In addition, numerous efforts have been focusing on the researches of graphene-based materials due to their excellent physical properties, such as 2D massless Dirac fermions, unique quantum Hall effect, and giant mobility of charge carriers.[3-5] Moreover, it can also be tailored into armchair graphene nanoribbons (AGNRs) and zigzag graphene nanoribbons (ZGNRs) according to the tailoring direction. Till now, many methods have been proposed to tune their electronic properties. Among them, doping,[6-8] adding electric field[9,10] and building defects[11-13] are the most common ones. Since the band gaps of GNRs are sensitive to their size and geometries, they are excellent candidates for semiconductor devices. To realize different functions of GNRs devices, it is important to manipulate the transport properties with a controllable manner. To date, many methods have been reported.[14-19] As a nano-electronic device, there is no denying that the particle property plays a leading role in transport properties, while its wave property cannot be ignored. The interference of waves propagated along different paths is one of the distinctive features. Thus quantum interference also plays an important role in the transport properties of devices.[20,21] Many works[22-25] have been devoted to investigate the quantum interference effect of graphene nanoribbons, and significant conclusions have been obtained. Hong et al.[21] designed a reversible device consisting of a carbon nanotube decorated with magnetic molecules, and obtained that the coupling between discrete spin-states of a magnetic molecule and continue bands of a carbon nanotube induces the spin-dependent Fano resonance, which can significantly change the conductance property. In this work, the quantum interference between CO molecules and the ZGNR device is investigated. The Fano resonance effect can be found in the transport property. Moreover, it is sensitive to the distance between CO molecules and the ZGNR device. A detailed theoretical explanation is also given. Researchers have shown that the electronic properties of ZGNRs are sensitive to the adsorption sites,[26,27] and adsorbing on both edges plays an important function for the electronic properties. Our research is based on the configuration of CO molecules adsorbed on one of the edges of the ZGNR device, as shown in Fig. 1, and each edged carbon atom is terminated by one hydrogen atom to avoid the dangling bonds. The whole system can be mainly divided into three parts: left electrode, central scattering region and right electrode. The left (right) electrode is semi-infinite 6-ZGNRs (here 6 is the ribbon width), and the scattering region is 6-ZGNRs adsorbed one CO molecule at one of the edges. The distance between the CO molecule and ZGNRs is labeled by $d$, as shown in Fig. 1. In the following we take $d=1.3$ Å as an example (short for M1) to have a detailed study.
cpl-34-9-097303-fig1.png
Fig. 1. Geometry structure of the CO molecule adsorbed on the ZGNR device.
Our calculations are performed with the Atomistix ToolKit (ATK) package, which is based on the density functional theory (DFT) combined with the non-equilibrium Green function (NEGF) method.[28,29] The structural relaxations of the systems are performed by fixing the CO molecule, then relaxing the device till the residual force on each atom is smaller than 0.01 eV/Å, and the total energy is converged to 10$^{-5}$ eV. The local-density approximation (LDA) is used for the exchange-correlation potential. The real space grid techniques are used with the energy cutoff of 150 Ry as a required cutoff energy in numerical integrations. The $k$-point sampling is 1, 1 and 100 in the $x$, $y$ and $z$ directions, respectively, where $z$ is the direction of electron transport. Open boundary conditions are used to describe the electronic and the transport properties of the device. A vacuum layer of 12 Å is added in $x$ and $y$ directions to avoid the interaction between adjacent ribbons. The wave functions of all the atoms are expanded by double-zeta polarized basis set. The NEGF-DFT self-consistency is controlled by a numerical tolerance of $4\times10^{-5}$ eV. The electron temperature is 300 K in the transport calculations. The transmission coefficient $T(E,V_{\rm b})$ can be obtained by $$\begin{alignat}{1} \!\!\!\!\!\!T(E,V_{\rm b})={\rm Tr}[{\it \Gamma} _{\rm L} (E,V_{\rm b})G^{\rm R}(E){\it \Gamma} _{\rm R} (E,V_{\rm b})G^{\rm A}(E)],~~ \tag {1} \end{alignat} $$ where $G^{\rm R(A)}$ is the retarded (advanced) Green's function of the scattering region, and ${\it \Gamma} _{\rm L(R)}$ is the contact broadening function associated with the left (right) electrode. The transmission spectrum $T(E,V_{\rm b})$ of the M1 device at zero bias has been calculated, as presented in Fig. 2(c). To have a better understanding about the transmission spectrum of the M1 device, we also give the transmission spectrum of the pristine 6-ZGNR device (short for M0) as a reference, which can be seen in Fig. 2(a). For the transmission spectrum of M0, there are perfect transmission steps. When a CO molecule is adsorbed on one edge of the ZGNR with $d =1.3$ Å, the transmission spectrum changes a lot, and one can find that the transmission steps are destroyed, instead two transmission dips appear at $-$0.12 eV and 1.12 eV energy sites, respectively. The transmission peak at 0 eV has decreased from three to two. All of these indicate that the CO molecule adsorption can effectively tune the transport properties of the ZGNR device. To have a further understanding about the transmission dips of the M1 device, we calculate the density of states (DOS) of M0 and M1, as shown in Figs. 2(b) and 2(d), respectively. Comparing Fig. 2(b) with Fig. 2(d), it can be seen that two more peaks appear in the M1 DOS, corresponding to the energies of $-$0.12 eV and 1.12 eV. More interesting is that the transmission dips of M1 appear just right in the energies of $-$0.12 eV and 1.12 eV, which is the same as that of the two more DOS peaks that appear. Thus we can predict that the DOS peaks appearing at the energies of $-$0.12 eV and 1.12 eV have a close relationship with the transmission dips of the M1 device.
cpl-34-9-097303-fig2.png
Fig. 2. [(a), (c)] The transmission coefficient of M0 and M1, [(b), (d)] the DOS of M0 and M1.
To continue probing the origin of the transmission dips in the M1 device, we give the project density of state (PDOS) and the local density of state (LDOS) of the M1 device in the following. Figures 3(a)–3(c) are the transmission spectrum of the whole device, PDOS of the CO molecule and PDOS of the carbon atoms in the up side (the green ball inside Fig. 3(c)), respectively. We can find that the peak of the DOS at 1.12 eV mainly comes from the CO molecule, while the peak of DOS at $-$0.12 eV mainly comes from the edge near the CO molecule of ZGNR. Figures 3(d)–3(f) present the LDOS of the device at different energy sites. Figure 3(d) is the LDOS of the device at 0 eV. We can see that the DOS is nonlocal along the down edge of the device, but is suppressed at the CO molecule adsorbed site along the up edge, thus the corresponding transmission peak decreases compared with that of the pristine ZGNR device. However, for the LDOS at 1.12 eV, shown in Fig. 3(e), there are quasi-bound states, the electrons in these states are mainly located on the CO molecule. Interestingly, the LDOS is mainly located at the carbon atoms on the up edge at the energy of $-$0.12 eV, as shown in Fig. 3(f). From the above, we obtain that the interaction between the CO molecule and the up edge of the ZGNR device results in the localization of electrons and then significantly changes the transmission spectrum.
cpl-34-9-097303-fig3.png
Fig. 3. (a)–(c) The transmission spectrum of M1, PDOS of the CO molecule and PDOS of green carbons shown inside in (c); (d)–(f) the LDOS for M1 at 0 eV, 1.12 eV and $-$0.12 eV, respectively. The isovalue is 2 Å$^{-3}$.
The Fano–Anderson model is used to explain this phenomenon.[30] As shown in Fig. 4, the ZGNRs adsorbed with CO molecule can be modeled as a 1D quantum wire (QW) with a side quantum dot (QD). The system Hamiltonian can be written as $$\begin{align} H=H_{\rm QW} +H_{\rm QD} +H_{\rm QD-QW},~~ \tag {2} \end{align} $$ where $H_{\rm QW}$ corresponds to the free-particle Hamiltonian on a lattice with spacing $d$. Their eigen-function is expressed as Bloch solutions $$\begin{align} |k\rangle =\sum\limits_{j=-\infty}^\infty {e^{ikdj}} |j\rangle,~~ \tag {3} \end{align} $$ where $|k\rangle$ is the momentum eigenstate, and $|j\rangle$ is a Wanier state localized at site $j$. The dispersion relation associated with these Bloch states is $$\begin{align} \varepsilon =2t\cos (kd),~~ \tag {4} \end{align} $$ where $t$ is the hopping coefficient. The energy level of the QD is $\varepsilon _{\rm L}$, and the coupling between QD and QW is $t$. The transmission coefficients can be obtained through a series of calculations[20,23] as $$\begin{align} T(\varepsilon)=\frac{({\varepsilon -\varepsilon _{\rm L}})^2}{({\varepsilon -\varepsilon _{\rm L}})^2+{\it \Gamma} ^2}.~~ \tag {5} \end{align} $$ From Eq. (5), we can obtain that the transmission coefficient is zero when $\varepsilon =\varepsilon _{\rm L}$. To the system of CO molecules adsorbed on ZGNRs, we can see the ZGNRs as QW, and the CO molecule as QD. When the energy eigenvalues of the CO molecule are equal to the energy eigenvalues of the ZGNRs, there will be a transmission dip in the corresponding energy site, which gives a better explanation for the former calculation. As can be seen from Eq. (5), if the parameter $\varepsilon _{\rm L}$ is changed, due to the fact that the energy spectrum of ZGNRs is continued, there still exist equal values for $\varepsilon$. Thus the transmission coefficient can also become zero.
cpl-34-9-097303-fig4.png
Fig. 4. Model of side-coupled quantum dot attached to a perfect quantum wire.
cpl-34-9-097303-fig5.png
Fig. 5. The transmission spectrum (a) and DOS of the CO adsorbed ZGNR device (b) with different $d$ in the zero bias region.
Next we study the influence of different distances between CO molecules and the ZGNR device (controlled by $d$) to the system's transmission properties. Figure 5(a) shows the transmission spectra of different systems at zero bias, where $d$ changes from 1.6 Å to 1.0 Å. From Figs. 5(a) and 2(c), we can see that when the distance is increased from 1.3 Å to 1.5 Å, there are still two transmission dips of the transmission spectrum at zero bias. Meanwhile, the transmission dips in the positive energy region are moving to a higher energy site. If the distance $d$ keeps on increasing, the transmission dip will disappear, as well as the transmission spectrum of M0, which indicates that the interaction between the CO molecule and ZGNR is too small to change the transmission spectrum of the device. However, when the distance $d$ decreases from 1.3 Å to 1.2 Å, the transmission dip still exists and the energy site is changed. When $d$ continues decreasing, the transmission dip in the positive energy region moves to a lower energy site, meanwhile the number of transmission dips is increased. Then we give the corresponding DOSs of different systems in Fig. 5(b). It is clear that the energy sites of transmission dips coincide with the energy sites of the DOS peaks. To have a clear picture of this phenomenon, we choose the positive energy region as an example; Fig. 6(a) gives the DOS of different devices and Fig. 6(b) shows the PDOS of the CO molecule from different devices. From them we can see that the energy sites of DOS peaks are completely consistent with that of PDOS of the CO molecule, thus the peaks appearing in the DOS mainly come from that of the CO molecule. Thus the variation of $d$ would cause the changes of the DOS distribution of the CO molecule in different energy sites, and results in the Fano resonance in different energy sites. From the above, we can see that the energy positions of the Fano resonance can be controlled by tuning $d$, and as $d$ decreases, the CO molecule is closer to the ZGNR device, the overlappings of wave functions are increased, thus the Fano dip shifts to a lower energy to stabilize the system.
cpl-34-9-097303-fig6.png
Fig. 6. The DOS (a) and PDOS (b) of the CO molecule in the positive energy area in the zero bias region with different $d$ for the devices.
In summary, we have investigated the influence to the transport property of the CO molecule adsorbed on one edge of the zigzag graphene nanoribbons device. It is shown that, when $d$ is larger than 1.6 Å, the CO molecule adsorption has no influence on the transport property of the ZGNR device. However, when $d$ is smaller than 1.6 Å, the mutual interaction between the CO molecule and the ZGNR edge can result in the Fano resonance effect. Moreover, with the decrease of $d$, the transmission dip moves to lower energy. Our further analysis indicates that the Fano resonance would be attributed to the interaction between the CO molecule and the ZGNR edge, which results in the localization of the electrons and significantly changes the transmission spectrum. These results would be useful for designing nano-devices.
References Electric Field Effect in Atomically Thin Carbon FilmsTwo-dimensional gas of massless Dirac fermions in grapheneGraphene: Status and ProspectsThe electronic properties of grapheneElectronic transport in two-dimensional grapheneDesigning of spin-filtering devices in zigzag graphene nanoribbons heterojunctions by asymmetric hydrogenation and B-N dopingElectronic structures and transport properties of armchair graphene nanoribbons by ordered dopingModulating the spin transport behaviors in ZBNCNRs by edge hydrogenation and position of BN chainHalf-metallic graphene nanoribbonsSpatially Separated Spin Carriers in Spin-Semiconducting Graphene NanoribbonsEffect of edge reconstruction and electron-electron interactions on quantum transport in graphene nanoribbonsTuning spin polarization and spin transport of zigzag graphene nanoribbons by line defectsElectronic transport properties on V-shaped-notched zigzag graphene nanoribbons junctionsGraphene-based bipolar spin diode and spin transistor: Rectification and amplification of spin-polarized currentElectronic transport properties on transition-metal terminated zigzag graphene nanoribbonsControl of electronic transport in nanohole defective zigzag graphene nanoribbon by means of side alkene chainCombined effect of strain and defects on the conductance of graphene nanoribbonsImproving spin-filtering efficiency in graphene and boron nitride nanoribbon heterostructure decorated with chromium-ligandPerfect spin filtering, rectifying and negative differential resistance effects in armchair graphene nanoribbonsTransport through a quantum wire with a side quantum-dot arrayFano-Resonance-Driven Spin-Valve Effect Using Single-Molecule MagnetsControl of electron transport through Fano resonances in molecular wiresQuantum conductance of graphene nanoribbons with edge defectsSpin interference and the Fano effect in electron transport through a mesoscopic ring side-coupled with a quantum dotQuasi-bound states and Fano effect in T-shaped graphene nanoribbonsEdge states and flat bands in graphene nanoribbons with arbitrary geometriesEdge-functionalized and substitutionally doped graphene nanoribbons: Electronic and spin propertiesAb initio modeling of quantum transport properties of molecular electronic devicesDensity-functional method for nonequilibrium electron transportFano resonances in nanoscale structures
[1] Novoselov K S, Geim A K, Morozov S V, Jiang D, Zhang Y, Dubonos S V, Grigorieva I V and Firsov A A 2004 Science 306 666
[2] Novoselov K S, Geim A K, Morozov S V, Jiang D, Katsnelson M I, Grigorieva I V, Dubonos S V and Firsov A A 2005 Nature 438 197
[3] Geim A K 2009 Science 324 1530
[4] Castro A H, Guinea F, Peres N M R K, Novoselov K S and Geim A K 2009 Rev. Mod. Phys. 81 109
[5] Sarma S D, Adam S, Hwang E H and Rossi E 2011 Rev. Mod. Phys. 83 407
[6] Zhang D, Long M Q, Zhang X J, Ouyang F P, Li M J and Xu H 2015 J. Appl. Phys. 117 014311
[7] Liu J, Zhang Z H, Deng X Q, Fan Z Q and Tang G P 2015 Org. Electron. 18 135
[8] Ouyang J, Long M Q, Zhang X J, Zhang D, He J and Gao Y L 2016 AIP Adv. 6 035116
[9] Son Y W, Cohen M L and Louie S G 2006 Nature 444 347
[10] Wang Z F, Jin S and Liu F 2013 Phys. Rev. Lett. 111 096803
[11] Ihnatsenka S and Kirczenow G 2013 Phys. Rev. B 88 125430
[12] Tang G P, Zhang Z H, Deng X Q, Fan Z Q and Zhu H L 2015 Phys. Chem. Chem. Phys. 17 638
[13] Zhang X J, Chen K Q, Tang L M and Long M Q 2011 Phys. Lett. A 375 3319
[14] Zeng M G, Shen L, Zhou M, Zhang C and Feng Y P 2011 Phys. Rev. B 83 115427
[15] Cao C, Chen L N, Long M Q, Huang W R and Xu H 2012 J. Appl. Phys. 111 113708
[16] Zou Y, Long M Q, Li M J, Zhang X J, Zhang Q T and Xu H 2015 RSC Adv. 5 19152
[17] Lehmann T, Ryndyk D A and Cuniberti G 2013 Phys. Rev. B 88 125420
[18] Zeng J, Chen L Z and Chen K Q 2014 Org. Electron. 15 1012
[19] Zhang D, Long M Q, Zhang X J, Cui L L, Li X M and Xu H 2017 J. Appl. Phys. 121 093903
[20] Orellana P A, Dominguez-Adame F, Gómez I and de Guevara M L L 2003 Phys. Rev. B 67 085321
[21] Hong K and Kim A W Y 2013 Chem. Int. Ed. 52 3389
[22] Papadopoulos T A, Grace I M and Lambert C J 2006 Phys. Rev. B 74 193306
[23] Li T C and Lu S P 2008 Phys. Rev. B 77 085408
[24] Ding G H and Dong B 2010 J. Phys.: Condens. Matter 22 135301
[25] Xu J G, Wang L and Weng M Q 2013 J. Appl. Phys. 114 153701
[26] Jaskólski W, Ayuela A, Pelc M, Santos H and Chico L 2011 Phys. Rev. B 83 235424
[27] Cervantes-Sodi F, Csányi G, Piscanec S and Ferrari A C 2008 Phys. Rev. B 77 165427
[28] Taylor J, Guo H and Wang J 2001 Phys. Rev. B 63 245407
[29] Brandbyge M, Mozos J, Ordejon P, Taylor J and Stokbro K 2002 Phys. Rev. B 65 165401
[30] Miroshnichenko A E, Flach S and Kivshar Y S 2010 Rev. Mod. Phys. 82 2257