# Vacancy-induced Fano resonances in zigzag phosphorene nanoribbons

###### Abstract

Motivated by recent scanning tunneling microscopy/spectroscopy experiments on probing single vacancies in black phosphorus, we present a theory for Fano antiresonances induced by coupling between vacancy states and edge states of zigzag phosphorene nanoribbons (zPNRs). To this end, in the first step, using the tight-binding Hamiltonian, we obtain an analytic solution on the lattice for the state associated to a single vacancy located in the bulk phosphorene which shows a highly anisotropic localization in real space. For a finite zigzag ribbon, in the absence of particle-hole symmetry, the localized state induced by vacancies can couple to the wave functions of the edge states which results in the formation of a new bound state. The energy of vacancy bound state lies inside the quasi-flat band composed of edge states when the vacancy locates sufficiently far away from the edge. Then, we employ the T-matrix Lippmann-Schwinger approach to obtain an explicit analytical expression for the scattering amplitude of the edge electrons of a zPNR by the presence of a single vacancy which shows a Fano resonance profile with a tunable dip. We demonstrate that varying the position of the vacancy produces substantially different effects on the resonance width, resonance energy position, and the asymmetry parameter of Fano line shape. Furthermore, the validity of the theoretical descriptions is verified numerically by using the Landauer approach.

###### pacs:

## I Introduction

The existence of great transport properties and potential applications of a new semiconducting material, called black phosphorus (BP), has attracted a lot of attention recently Li2014 (); Li2016 (); Liu2014 (); Xia2014 (); Neto-2014 (); Ling2015 (). BP has a layered structure with van der Waals interaction between individual atomic layers and it is possible to exfoliate a single atomic layer from its bulk which is called phosphorene Li2014 (); Liu2014 (). In this new specimen of two-dimensional (2D) material, each phosphorus atom is covalently bonded with three nearest neighbors via hybridization to form a puckered 2D honeycomb structure. This unique property of anisotropy on one hand, and a tunable and large direct band gap, on the other hand, makes phosphorene a desirable candidate material for different applications with specific electronic, mechanical, thermal, and transport features Carvalho2016 (); Neto2014 (); Qiao2014 (); Peeters2014 (); Guinea2014 (); Yang2014 (); Katsnelson2015 (). Furthermore, confining one of the dimensions of phosphorene to a ribbon with the zigzag edge which is called zigzag phosphorene nanoribbon introduces degenerate quasi-flat bands in the middle of the gap detached entirely from the valence and conduction bands Ezawa2014 (); Peeters-skew (). Due to the presence of such edge modes, the quantum transport properties of zPNRs resembles a quasi-one-dimensional (Q1D) system in this energy window Ezawa2014 (); Asgari2018 (); Amini2018 ().

Despite the vast literature regarding different theoretical and experimental researches on perfect PNRs as mentioned above, the body of work dedicated to their defective counterparts is surprisingly small. For instance, the electronic band structure of defective phosphorene, as well as the formation energy of atomic vacancies, are studied numerically in Ref. Hu2015 (); Wang2015 (). But among these studies, there is an experiment carried out recently Katsnelson2017 () which is the observation of atomic vacancies on the surface of bulk BP. The spatial distribution of wave function amplitudes around a vacancy which is observed by scanning tunneling microscopy tomographic imaging exhibits strong anisotropy. Another significant aspect of the presence of vacancies is their effects on quantum transport through the edges of the zPNRs which is computed numerically in Ref. Peeters2018 () using the tight-binding approach. It is found that single vacancies can create quasi-localized states which can affect the transport properties of zPNRs by introducing antiresonant dips in the edge conductance of the ribbons.

However, the results outlined in the previous paragraph bring up an important question of whether the conductance dips computed for the electronic transmission through the edges of zPNR is due to the coupling of anisotropic vacancy state to the continuum states of the edges resulting in the so-called Fano antiresonance in the wave propagation. Indeed, Fano resonance Fano () is one of the most important resonance mechanisms in which a discrete level couples to an energetically close-by continuum of states giving rise to new scattering paths. It is quite predictable that a single vacancy can introduce a new localized state within the energy window of the edge modes. Therefore, a possible scenario to understand the corresponding conductance dips of vacant zPNRs can be due to the coupling of the edge states forming a quasi-flat band to the vacancy wave function. In this work, we present the theory of such resonance scattering due to the vacancy defect leading to Fano resonance in the transmission of edge channels in zPNRs. To this end, we first use a tight-binding model to obtain and analyze in detail the localized state induced by a vacancy defect in the energy interval of the edge band in zPNRs. Then, we employ the general scattering theory based on the Green’s function approach and the Lippmann-Schwinger equation to investigate the antiresonance dips in the electronic transmission through the edge states of zPNRs due to the presence of such a localized state.

The paper is organized as follows. In section II, we provide the rigorous derivation of localized wave function around a single vacancy in phosphorene using the tight-binding Hamiltonian. In this section, the energy shift of such a bound state due to the absence of electron-hole symmetry is explained. In section III, we generally introduce the theory of a Fano antiresonance induced by the coupling of an impurity state to a continuum of states using the T-matrix Lippmann-Schwinger approach. We employ this formalism in section IV to obtain the Fano antiresonance dips in the electronic transmission through the edge states of zPNR which is induced by a single vacancy near the edge. Section V is devoted to discussion on the obtained results and, finally, we wrap up the paper with the summary and concluding remarks in section VI.

## Ii Localized states around vacancies in bulk phosphorene

Phosphorene consists of a non-planar puckered honeycomb lattice of phosphorus atoms with two in-equivalent sublattices denoted by and as shown in Fig. 1. We choose the and axes along the zigzag and armchair chains respectively. Now, it is possible to describe each lattice site by where and are armchair and zigzag chain numbers respectively and . The energy band structure of phosphorene can be well described using a tight-binding model Ezawa2014 (); Rudenko2014 ()

(1) |

where represents summation over only considered neighbors. Here, is the hopping energy between sites and , is the electron creation (annihilation) operator in site , and stands for hermitian conjugate. It is shown that considering only five hopping integrals eV, eV, eV, eV, and eV in this model, which is shown in Fig. 1, provides a reasonable description of the phosphorene band structure Rudenko2014 (). However, for analytical calculations, it is reasonable to neglect hopping terms and which have a smaller contribution to the Hamiltonian in comparison to and terms. It is important to note that the presence of breaks both electron-hole and lattice inversion symmetries of phosphorene and we will take into account its effect perturbatively later. Therefore, we can rewrite the above Hamiltonian keeping only the important terms of first, second and fourth nearest neighbors as

(2) |

Without loss of any generality, let us now consider a single vacancy placed at site . This is equivalent to turning the hopping terms that connect this site to its neighbors off. In the presence of particle-hole symmetry, , this implies the appearance of a new zero energy state which has only non-zero amplitudes on the sublattice Pereira2008 () and we call it . Therefore, it should satisfy the condition

(3) |

Later, we will calculate the effect of perturbatively which shifts this energy toward the negative values.

Indeed, a single vacancy can divide the phosphorene plane into an upper semi-infinite ribbon with beard edge and a lower one with the zigzag edge which is shown by black dashed lines in Fig. 1. The upper semi-infinite ribbon with beard edge does not support localized edge states Ezawa2014 (). So, the vacancy wave function amplitudes are zero in the upper half-plane. On the other hand, due to the absence of hopping between site and the vacant site, a possible solution in the lower half-plane can start with non-zero amplitude on this site, namely . Also Eq. (3) implies that which results in . With the same kind of reasoning it is evident that . Next, one can apply the same procedure to obtain the amplitudes of the wave function on the next zigzag chains which we do not pursue it here. Instead, we focus on a more systematic way of constructing the vacancy wave function Pereira2006 (); Vozmediano () below.

Let us consider the lower semi-infinite zigzag ribbon in which the amplitudes of the vacancy wave function is non-zero. The typical edge states of this ribbon have the following wave functions Amini2018 ()

(4) |

where is the wave vector measured in units of (and is the lattice constant) along the horizontal axis. Here, is a constant value equal to 0 (0.5) for even (odd) , , and . Also, the superscript denotes the type of sublattice in which the wave function is localized and is the zigzag chain number measured from the edge (starting from zero). Therefore, the idea is to express the zero energy vacancy state as a superposition of these edge states with zero components everywhere along the zigzag edge except on site . The following combination of the edge states can satisfy such a condition

(5) |

in which is the normalization coefficient given by

(6) |

Now, we can obtain the amplitudes of the wave function in Eq. (5) on different lattice sites. Let us start from

(7) |

This ensures our desired form of the wave function on the zigzag chain. Clearly, we can go on and calculate the wave function amplitudes on the zigzag chain as

(8) |

This is equivalent to the results of the first method which we described before. Fig. 2 represents the corresponding wave function amplitudes around a vacancy located on the site which is denoted by symbol. As it is shown, the vacancy wave function is highly anisotropic and non-zero in a region like a triangle over the apposite sublattice in which the vacancy is located. This vacancy wave function resembles the topological boundary state emerges at the corners of a nanodisk Ezawa2018 ().

Now, It is straightforward to calculate the energy correction to this zero-energy mode by the presence of term in Eq. (2). In doing so, we employ the first order perturbation theory which leads to

(9) |

with in which the integration can be carried out using contour integration Amini2018 (). Eqs. (9) and (5) will be used in the next section to calculate the electronic transmission through the edge states of a zPNR in presence of a single vacancy.

However, before ending this section, let us consider a finite ribbon of phosphorene with zigzag boundary and locate a vacancy far away from the edges. Fig. 3, represents the corresponding LDOS of such a vacancy state which is obtained numerically in agreement with the analytical solution. We can compare our results with the LDOS of the edge states shown along the zigzag edges of the ribbon.

## Iii Fano resonances due to the coupling of a localized state with continuum states

In order to study the phenomenon of Fano resonance, we can consider the impurity state as a quantum dot (QD) coupled to a Q1D quantum system. The Hamiltonian describing the system is

(10) |

where the first term is the Hamiltonian of a Q1D lattice with translational invariant which can be written using the plane wave basis as

(11) |

in which is the lattice constant and is the energy dispersion relation for electrons moving on the Q1D lattice in the absence of impurity. The second term

(12) |

describes a QD with energy and, finally, the last term describes the coupling between QD and Q1D as follows

(13) |

in which is the coupling strength at momentum . From now on, we set the lattice spacing to one, .

Since we are interested in the transmission coefficient (where is the transmission amplitude), we can start with the Lippmann-Schwinger equation in the scattering theory

(14) |

where

(15) |

in which is the Green’s operator in the absence of coupling between Q1D and QD. Since , it is formally possible to write

(16) |

in which

(17) |

and

(18) |

are the Green’s functions of clean Q1D and QD respectively. With these expressions in hand, the first term of the series in Eq. (15) can be evaluated as

(19) |

where

(20) |

and

(21) |

Doing the same analysis for the second term in Eq. (15) and collecting these equations, we find,

(22) |

where we defined

(23) |

According to Eq. (14), the transmission amplitude at energy can be obtained as,

(24) |

for an incident state such that , hence,

(25) |

By having a quick look at the final expression for , we can see that the coupling of an extra state (due to the presence of QD) with the continuum states with the same energy can induce an antiresonance in the transmission amplitude. We will use this formalism in the next section to find this kind of resonant scattering due to the presence of vacancies in zPNRs.

### iii.1 Antiresonance scattering of a 1D chain coupled to a QD

There is, however, at least one task remaining, namely, to show that the above formalism works fine. In order to check our formalism, let us consider the trivial example of scattering of an electron traveling along a 1D chain by a QD (as shown in Fig. 4) which can be described by the FanoâAnderson tight-binding Hamiltonian FlachRMP (). The energy dispersion of this model is given by

(26) |

where is the hopping energy cost of an electron to the nearest neighbor sites. We also consider a QD which is attached to the chain with on-site energy and coupling constant . Therefore, the resulting transmission is obtained by substituting these parameters into Eq. (25) which yields the following expression

(27) |

where . Hence, we arrive at the exact solution Orellana () for the Fano antiresonance in the transmission amplitude of this system.

## Iv Fano resonances in vacant zPNRs

In this section, we use the above formalism to study the Fano resonance phenomenon which is induced by the presence of a single vacancy defect in ZPNRs. In so doing, our approach is very similar to what was followed in the previous section for a 1D chain coupled to a QD. Let us consider a single vacancy which is located at the th zigzag chain with respect to the edge of ribbon as shown in Fig. 5.

As we already discussed, this vacancy introduces an additional impurity state whose wave function is . Now, the significant idea is that we would replace the role of Q1D and QD systems by the edge modes of the zPNR and vacancy state respectively. Accordingly, we need to replace the required quantities , and in Eq. (25) with their counterparts of the vacant zPNR. This will be demonstrated shortly.

Let us start with first quantity which is the electronic energy dispersion relation of the edge states in zPNR which is obtained in Ref. Amini2018 () as

(28) |

Before proceeding to the other quantities, we need to solve a technical problem. It is important to note that according to Eqs. (4) and (5), the states and are not orthogonal, namely . Thus, we need to define an orthogonal state as

(29) |

such that the required condition is satisfied. Here,

(30) |

is a constant that can be used to normalize the wave function. It is now possible to obtain the energy of this bound state which resembles the quantity . Due to the presence of the term the corresponding bound-state energy of is changed from zero and can be obtained perturbatively as

(31) |

where s are given by

(32) |

and

(33) |

is given by Eq. (9) and . It is now important, however, to note that if the bound state energy lies within the energy interval of the edge modes, namely , then it is expected to observe the Fano antiresonance for the edge modes. For instance,

(34a) | |||

(34b) | |||

(34c) | |||

(34d) |

which shows that for one can not observe the Fano antiresonance because it’s corresponding bonding energy lies outside the bandwidth.

We are now in a position to calculate the only remaining quantity which is the coupling parameter between the vacancy state and edge states, namely

(35) |

This results in

(36) |

Collecting these results and inserting into Eq. (25), we can obtain the transmission amplitude through the edge states of a zPNR in presence of a single vacancy which is located on th zigzag chain. Thus,

(37) |

where and the integration can be carried out analytically using contour integration Amini2018 (); Amini2018-2 ().

## V Results and discussion

In this section we consider a zPNR in presence of a single vacancy which is located at site . As we already discussed, only when the Fano antiresonance takes place. Therefore, for simplicity, we first assume which is the first place in which the vacancy can induce antiresonance. In Fig. 6 we show the corresponding LDOS of the edge states and localized vacancy state obtained numerically. In order to obtain the transmission amplitude of such a system, we can insert Eqs. (31) and (36) into Eq. (37) which results in

(38) |

in which and should be evaluated for and is rather a complicated expression which we do not present here. Instead, we present a graphical representation of the transmission coefficient using the transmission amplitude in Eq. (38) as a function of energy in Fig. 7. It is clearly shown in Fig. 7 that there exists a Fano antiresonance around resonance energy (and close to it’s corresponding bound state energy ) for which the transmission through the edge channel vanishes. This is due to the destructive interference of the traveling wave which can either bypass the vacancy state or populate it and return back to the continuum. Let us now check the deviation of represented expression from the standard line shape of the Fano resonances FlachRMP () of the form

(39) |

where and , , and are the line width, peak position, and asymmetry parameter respectively. Fig. 7 shows the fitted curve (dashed lines) with this Fano formula which shows an excellent agreement to a Fano line shape with asymmetry parameter . This difference with the case of an exactly 1D chain which shows a symmetric transmission profile is reasonable since the zPNR is a Q1D system and the edge states have non-zero amplitudes within the localization radius near the edge. In order to check our analytical expression, we also use the numerical Landauer approach that works on the basis of the recursive Greenâs function technique and we have used it in Ref. Amini2018 (). Fig. 7 (discrete circles) shows the resulting transmission which is obtained from the numerical calculation and is in perfect agreement with the analytical solution.

Next, we study the influence of varying the vacancy position on the Fano lineshape properties. To this end, we change the vacancy position to different distances from the edge, namely increasing to . Fig. 8 presents the transmission probability versus the incident electron energy for various vacancy positions. It is clear that displacing the vacancy away from the edge of the zPNR causes several important differences: (i) a shift in the resonance energy toward the center of the edge band; (ii) suppression of the asymmetry parameter , and (iii) decreasing the line width of the peak, . The reason is that by increasing the distance of vacancy from the edge, the coupling parameter decreases due to the suppression of wave functions overlap and since the resonance width is directly proportional to the coupling strength, therefore decreases. Furthermore, far away from the edge, one can consider the edge states as an exactly 1D chain, therefore, the transmission profile should be symmetric and locate at the band center. Before ending this discussion, we should emphasize that our results are completely consistent with the results of Ref. Peeters2018 ().

Finally, it is worth pointing out that the presence of hopping parameters and does not change the overall picture of the transmission coefficient in general. To clarify this point, we present a plot of transmission coefficient versus energy for both cases, namely the case in which and the case with non-zero and , in Fig. 9. As presented in Fig. 9, the resonant peaks of transmission coefficient in presence of and (which is shown with dashed lines) takes place with a small energy shift. This is reasonable because the presence of and increases the energy bandwidth of the edge band.

## Vi Summary

In this paper, we have demonstrated the existence of a new type of impurity state associated with vacancies in monolayer phosphorene which exhibits a strongly anisotropic localization profile in real space. We found that the localized state induced by vacancies can overlap with the wave functions of the edge states in a confined ribbon of phosphorene which results in the formation of dips in the electronic transmission through the edge modes. The latter takes place when the energy of the vacancy bound state lies inside the quasi flat band composed of the edge modes. We have derived analytically and validated numerically the theory of such antiresonances in zPNRs with a vacancy using the T-matrix Lippmann-Schwinger approach. Modeling the continuum states of the edges of zPNRs as a Q1D system and the vacancy state as a QD, we found that it is exactly a Fano phenomenon due to the coupling of vacancy state to the continuum of the edge modes. We have tried to shed light on various effects of displacing the vacancy on such Fano antiresonances that occur in a zPNR with single vacancies.

###### Acknowledgements.

We gratefully acknowledge discussion with F. M. Peeters and L. L. Li. MA acknowledges the support of the Abdus Salam (ICTP) associateship program. MSh would also like to acknowledge the office of graduate studies at the University of Isfahan for their support and research facilities.## References

- (1) L. Li, Y. Yu, G. J. Ye, Q. Ge, X. Ou, H. Wu, D. Feng, X. H. Chen, and Y. Zhang, Nat. Nanotechnol. 9, 372 (2014).
- (2) L. Li, F. Yang, G. J. Ye, Z. Zhang, Z. Zhu, W. Lou, L. Li, X. Zhou, K. Watanabe, T. Taniguchi, K. Chang, Y. Wang, X. H. Chen, and Y. Zhang, Nat. Nanotechnol. 11, 593 (2016).
- (3) H. Liu, A. T. Neal, Z. Zhu, Z. Luo, X. Xu, D. Tomanek, and P. D. Ye, ACS Nano 8, 4033 (2014).
- (4) F. Xia, H. Wang, and Y. Jia, Nat. Commun. 5, 4458 (2014).
- (5) S. P. Koenig, R. A. Doganov, H. Schmidt, A. H. Castro Neto, and B. Zyilmaz, Appl. Phys. Lett. 104, 103106 (2014).
- (6) X. Ling, H. Wang, S. Huang, F. Xia, and M. S. Dresselhaus, Proc. Natl. Acad. Sci. (USA) 112, 4523 (2015).
- (7) Alexandra Carvalho, Min Wang, Xi Zhu, Aleksandr S. Rodin, Haibin Su and Antonio H. Castro Neto, Nat. Rev. Mater. 1, 16061 (2016).
- (8) J. Qiao, X. Kong, Z. X. Hu, F. Yang, and W. Ji, Nat. Commun. 5, 4475 (2014).
- (9) A. S. Rodin, A. Carvalho, and A. H. Castro Neto, Phys. Rev. Lett. 112, 176801 (2014).
- (10) D. Cakir, H. Sahin, and F. M. Peeters, Phys. Rev. B 90, 205421 (2014).
- (11) T. Low, R. Roldan, H. Wang, F. Xia, P. Avouris, L. M. Moreno, and F. Guinea, Phys. Rev. Lett. 113, 106802 (2014).
- (12) R. Fei and L. Yang, Nano Lett. 14, 2884 (2014).
- (13) S. Yuan, A. N. Rudenko, and M. I. Katsnelson, Phys. Rev. B 91, 115436 (2015).
- (14) M. Ezawa, New Journal of Physics 16, 115004 (2014).
- (15) M. M. Grujic, M. Ezawa, M. Z. Tadic, and F. M. Peeters, Phys. Rev. B 93, 245413 (2016).
- (16) Z. Nourbakhsh and R. Asgari, arXiv:1803.00751(2018).
- (17) M. Amini and M. Soltani, arXiv:1810.03042 (2018).
- (18) W. Hu and J. Yang, J. Phys. Chem. C 119, 20474 (2015).
- (19) V. Wang, Y. Kawazoe, and W. T. Geng, Phys. Rev. B 91, 045433 (2015).
- (20) B. Kiraly, N. Hauptmann, A. N. Rudenko, M. I. Katsnelson, and A. A. Khajetoorians, Nano Lett. 17, 3607 (2017).
- (21) L. L. Li, and F. M. Peeters, Phys. Rev. B 97, 075414 (2018).
- (22) U. Fano Phys. Rev. 124, 1866 (1961).
- (23) A. N. Rudenko, M. I. Katsnelson, Phys. Rev. B 89, 201408 (2014).
- (24) V. M. Pereira, J. M. B. Lopes dos Santos, and A. H. Castro Neto, Phys. Rev. B 77, 115109 (2008).
- (25) V. M. Pereira, F. Guinea, J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, Phys. Rev. Lett. 96, 036801 (2006).
- (26) E. V. Castro, M. P. LÃ³pez-Sancho, and M. A. H. Vozmediano, Phys. Rev. Lett. 104, 036802 (2010).
- (27) M. Ezawa, Phys. Rev. B 98, 045125 (2018).
- (28) A. E. Miroshnichenko, S. Flach, and Y. S. Kivshar Rev. Mod. Phys. 82 2257 (2010).
- (29) P. A. Orellana, F. DomÃnguez-Adame, I. GÃ³mez, and M. L. LadrÃ³n de Guevara, Phys. Rev. B 67, 085321 (2003).
- (30) M. Amini, M. Soltani, E. Ghanbari-Adivi, and M. Sharbafiun arXiv:1810.03486 (2018).