Abstract
The terahertz band is an increasingly important spectrum in a wide range of applications from bioimaging and medical diagnostics to security and wireless communications. We propose a tunable terahertz coherent radiation source based on graphene plasmon-induced transition radiation. The transition radiation in terahertz regime arises from the graphene plasmons, which are excited by a normally incident bunched electron beam. We analyze the field-intensities and spectral-angular distributions of the transition radiation with respect to Fermi energy, substrate dielectric permittivity, and electron bunch energy for both the coherent and incoherent radiation. The effect of electron bunching on the radiation pattern is discussed. The mechanism of plasmon frequency-selective transition radiation is discovered.
© 2017 Optical Society of America
1. Introduction
Plasmons are the collective excitations of conduction electrons at the metal’s surface [1]. They can be excited either by an incident electromagnetic wave or by an electron beam on the metal-dielectric interface, to form coupled oscillations called surface plasmon polaritons (SPPs) [2,3]. The confinement of SPPs on the metal-dielectric interface not only allows propagation of light below the diffraction limit, it also enhances radiative emission from the surface as well. One of the plasmonic materials available today is graphene. Different from metals, the electrons in graphene are uniquely two-dimensional massless fermions that exhibit exceptional electronic properties [4], offering unduly performance enhancement in photonics, optoelectronics and plamonics [5–9]. Also, unlike in noble metals, plasmonic excitations in graphene can be tuned by shifting graphene’s Fermi energy through chemical doping or electrostatic bias. Within the carrier density of lowly-doped graphene, the graphene plasmon frequency lies right at the terahertz (THz) frequency range. As such, electrons moving uniformly and parallel to the graphene sheet on a substrate are predicted to generate SPP-enhanced THz radiation [10–16]. Further, inelastic electron tunneling has been proposed as a low-energy pathway for the highly-efficient excitation of surface plasmons at the THz frequencies [17]. On the other hand, when electrons are in normal or oblique incidence to the graphene surface, both SPP and transition radiation (TR) will be generated [18]. The TR is determined by the polarization current of the graphene’s surface, which is influenced by the plasmons, and therefore are dynamically tunable through electrostatic gating. The ability to effectively generate tunable coherent THz waves from graphene has recently became one of the most attractive subjects, as they have wide applications in bio-imaging, biomolecule identification, security inspection and ultrafast wireless communications.
In this work, we propose a tunable THz coherent radiation source based on graphene SPP-induce TR. From classical electromagnetic theory, we take into account the polarization current of the graphene’s surface to study the TR from graphene. The properties of THz TR are analyzed for the case when a travelling electron-bunch is at normal incidence to the graphene boundary. In Section 2, we theoretically study the field distribution and spectral-angular distribution of TR excited by a bunch. Then, in Section 3, we present the calculation results and discussion.
2. Transition radiation from graphene by a bunch
With the remarkable development of ultrafast technology, emission from an electron beam can now be precisely controlled. Recently, a TR-based THz source with energies in the order of μJ/pulse has been reported [19]. This coherent TR is produced by laser-driven electron bunches across a vacuum-solid boundary, with a reasonable number of electrons (Ne ~1.32 × 1010) in the bunch.
When a bunched beam transmits through a graphene layer on a dielectric substrate, it excites the in-plane graphene SPPs. Then, the polarization current density induced by SPPs will result in out-of-plane TR. The induced two-dimensional polarization current density J// is given by σE// according to Ohm’s law where σ is the conductivity of graphene and E// is the in-plane electric field.
A. Radiation field
Let us consider a graphene layer as an infinitely thin conducting sheet placed at the interface (z = 0) between vacuum (z>0) and a dielectric substrate (z<0) with the relative permittivity ε1 and ε2, respectively, as shown in Fig. 1. When an electron bunch normally passes through the graphene monolayer, it induces a polarization current density, which serves as source of TR.
The Maxwell’s wave equation can be written for Fourier images of electric field and polarization current density as follows
For an electron bunch, the charge density and current density in time domain can be found fromwhere e is charge of electron, the subscript j denotes the jth electron, and r0j and vj are the position vector and the velocity of the jth electron, respectively.Each electron in the bunch can be described as an object in a six-dimensional position-momentum space, which consists of a three-dimensional position distribution and a three-dimensional momentum distribution. For simplicity, we will assume that there is no correlation between the position and the momentum of electrons and all electrons move with the same velocity. Following this assumption, the bunch distribution can be reduced to a three-dimensional position distribution. Substituting the Fourier transform of the charge density and the beam current density in Eq. (2) into the Fourier transform of the wave equation in Eq. (1), the Fourier image of electric field accompanied by a traveling bunch with velocity v is now written as
where is the wave vector, , and , as shown in Fig. 1. The subscripts // and z represent the direction parallel with the plane of graphene layer and the direction of the bunch trajectory, respectively.Subsequently, the magnetic field of a travelling electron bunch is derived as
where the subscript φ denotes the angle from graphene plane.The TR fields can be calculated by solving the wave equations in both the vacuum and substrate and by applying the continuity condition of the tangential electric and magnetic fields across the boundary of the graphene layer. The complete solutions of the wave equations contain both the charged particle fields and the homogeneous TR fields. From the continuity equations across the interface (z = 0), according to the fact that the homogeneous fields are convergent,, the homogeneous radiation fields induced by graphene in vacuum (z>0, zone I) and substrate (z<0, zone II) are written as:
where is the propagation wave vector. Factors A and B are determined by the boundary conditions on the vacuum-substrate interface:and , and and are the total electric and magnetic field in the zone I(II), respectively.The spatial distribution of the electron beam can be considered as a Gaussian distribution through, where σz and σr are the root-mean-square longitudinal length of beam duration and transverse width, respectively. By solving Eq. (2-5) with the boundary conditions, we get the factors A and B of the radiation field as
whereσg is the intraband conductivity of graphene monolayer, A0 and B0 are the amplitudes of the radiation field excited by a single electron, N is the number of electrons, and F = F//Fz is the spatial factor which is closely related to the Fourier transform of the electron beam’s spatial distribution. For a Gaussian spatial beam distribution, the spatial factor is expressed byWhen the denominators of A0 and B0 are equal to zero, the corresponding k// will represent the wave vector of the graphene SPP. One can see that when the evanescent waves from an incident electron bunch excites the SPPs on graphene layer, it will spontaneously induce a polarization current density according to the Drude conductivity and the graphene plasmon frequency, which serves as the source of the TR at the same frequency. After straightforward mathematical manipulations, the dispersion of the graphene SPPs on the interface can be written as:
Within linear response theory, graphene’s dynamic conductivity is the sum of interband electronic transitions and intraband Drude-like processes. However, in the THz range, the intraband contribution plays a dominant role in the optical conductivity. Thus a Drude-like intraband conductivity of graphene layer can be obtained [20]
where τ is the carrier relaxation time, and is the Fermi energy, which is determined by carrier concentration ns and Fermi velocity vF.Finally, the electric field of TR in the graphene plane can be expressed as
where is the radiation field caused by a single electron, and J0 is the Bessel function. When the condition NF>1 is satisfied, the radiation field from the bunched electrons must be larger than that from an individual electron. In this way we can take into account the beam bunching effect on the TR field. Assuming the number of electrons in the bunch to be ~1010, the bunch parameters are σr = 30μm and σz = 3μm, which corresponds to a moderate beam energy with a fs pulse.B. Spectral-angular distribution
The total far-field TR energy from an electron beam through the graphene layer can be obtained from the integral of the Poynting vector
where and are the radiation fields in parallel with the graphene plane and the direction of the bunch trajectory, respectively.The spectral angular radiation distribution in the frequency domain is given by the expression
where Ω is the solid angle. One can get the distribution of the TR from graphene over the angles and frequencies byThe summations over electrons in Eq. (9) and Eq. (12) are already taken into account by the three-dimensional electron position distributions. The radiation distribution yieldswhere, E0// and E0z are the average radiation fields from a single electron parallel to xy-plane and z-plane, respectively. The first term on the right-hand side of Eq. (13) is the contribution from incoherent radiation, which scales as N, and latter term is from coherent radiation, which scales as N2. One can see that the coherent radiation is related to the bunch size, while the incoherent radiation is not.In our calculations, we ignore contributions from Cherenkov radiation in the substrate, Bremsstrahlung as well as characteristic X-rays from graphene or the substrate, and assume graphene is not damaged for electron bunch energies up to 100keV.
3. Results
Given reasonable parameters for graphene’s conductivity, τ = 0.25ps, EF = 0.1eV, and dielectric constants ε1 = 1 (z>0), ε2 = 4 (z<0), we obtained numerical results for the TR field distributions in the r-z plane at 5THz. The results for an individual electron and an electron bunch (N = 1010, σr = 30μm and σz = 3μm), both with a beam energy of 3keV, are shown in Fig. 2(a) and 2(b) respectively. From the figures, it is observed that the distribution angles of the radiation fields in vacuum are larger than that in the substrate for both cases. Meanwhile, the amplitude of the radiation fields excited by an electron bunch is up to 107 higher than that by a single electron. Moreover, the radiation angles for bunched electrons are much smaller than that for a single electron. The inset can clearly identify the field contour in the middle part of the figure where the contour is overlapped.
Next, we perform numerical calculations on the spectral-angular distribution for incoherent and coherent radiations with a beam energy 3keV and r = 500μm, while the other parameters remain the same. The results are shown in Fig. 3. The radiation energy concentrates in a small angular range near θ = 0 (in the backward direction) and θ = π (in forward direction). The intensity of coherent radiation at the lower THz band is far stronger than that of the incoherent one, and the angular distribution of the coherent radiation has two individual frequency peaks near θ = 0 and θ = π respectively. In contrast, incoherent radiation still exists in the higher THz band beyond the graphene plasmon operating frequency.
The Poynting fluxes with bunch energies of 1keV, 10keV and 100keV are shown in Fig. 4. The intensity of TR is seen to increase with bunch energy. We also find that the radiation in backward direction increases linearly with the bunch energy, while radiation in forward direction increases nonlinearly. Figure 4(d) shows the TR from 10keV electron bunches with respect to angle for various frequencies, having independent peaks in backward (BW) and forward (FW) directions. At 0.3-0.5THz, the spectral-angular distribution of the radiation intensity mainly lies in the BW direction while at 1.1-1.3THz it lies in the FW direction. In all cases, TR in direction parallel with the graphene layer (θ = π/2) is negligible. The inset in Fig. 4(d) clearly shows the sharp change of the Poynting fluxes at all frequencies near the angle θ = π/2 due to the different dielectric property on two sides of the graphene layer.
We proceed to calculate the TR Poynting flux with varying Fermi energy. In Fig. 5, the three grouped curves, from bottom to top, correspond to the bunch energy of 1keV, 10keV and 100keV respectively. In each group, curves are obtained with varying Fermi energy from 0.1eV to 0.3eV. The peak intensity in the BW direction shown in Fig. 5(a) increases with the Fermi energy and is proportional to the electron bunch energy. The shape of the TR spectrum remains the same throughout the varying Fermi energy and electron bunch energy. In contrast, the shape of the radiation spectrum changes with the electron bunch energy for TR in the FW direction, as shown in Fig. 5(b). The peak intensity decreases with increase of Fermi energy, and the peak frequency shifts upwards with the increase of electron bunch energy. For the electron bunch energy of 100keV, the radiation power in the FW direction is about 400 times higher than that in the BW direction. Overall, these results imply that the radiation intensity and peak frequency can be tuned by external doping (via chemical or electrostatic gating) of graphene. Additionally, as well-known, Bremsstrahlung and characteristic X-rays can be generated from graphene or the substrate by the kilo-voltage energy of electron beam. These radiations mainly distribute in the angle near θ = π/2, while TR mainly in forward and backward, thus they could be detected individually in different azimuth angle.
The graphene SPP dispersion is also influenced by the dielectric permittivity of substrate. Here we discuss the effect of dielectric substrate to the TR. For the case when the electron bunch energy is 10keV, the Poynting flux with varying dielectric permittivity in the BW and FW directions is shown in Fig. 6. One can see that, the radiation intensity in the BW direction increases with increase of dielectric constant, while in the FW direction it first increases and then decreases. This is due to the discrepancy between the two boundary conditions for the graphene SPP dispersion for the BW and FW radiation respectively. For the radiation in FW direction, the intensity maximum can be obtained near ε2 = 7.5 and there are no significant changes for ε2 from 6 to 10.
4. Summary
In this paper, we proposed a tunable coherent THz radiation source based on graphene SPP-induced TR. We present a detailed theoretical analysis of TR in the THz frequency range, which is excited by bunched electrons uniformly passing through a graphene layer with normal incidence. The radiation fields, spectral-angular distribution and Ponyting flux of the TR are derived and numerically analyzed for coherent radiation in the backward and forward directions. The results show that the TR from graphene has unique characteristics arising from the high tunable graphene plasmon dispersion. The frequency of TR can be tuned externally by Fermi energy, dielectric constant of substrate, electron bunch size and energy in the frequency range of 0.4-2THz. The THz TR will find potential uses in THz pulsed imaging (TPI) [21,22], biological imaging [23,24], security and wireless communications applications.
Funding
Fundamental Research Funds for the Central Universities (ZYGX2015J038); National Key Basic Research Program of China (2013CB933603, 2014CB339801); China Scholarship Council (201506075086).
References and links
1. D. Bohm and D. Pines, “A Collective Description of Electron Interactions: III. Coulomb Interactions in a Degenerate Electron Gas,” Phys. Rev. 92(3), 609–625 (1953). [CrossRef]
2. R. H. Ritchie, “Plasma Losses by Fast Electrons in Thin Films,” Phys. Rev. 106(5), 874–881 (1957). [CrossRef]
3. C. J. Powell and J. B. Swan, “The Origin of the Characteristic Electron Energy Losses in Aluminum,” Phys. Rev. 115(4), 869–875 (1959). [CrossRef]
4. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, “Two-dimensional gas of massless Dirac fermions in graphene,” Nature 438(7065), 197–200 (2005). [CrossRef] [PubMed]
5. M. Liu, X. Yin, E. Ulin-Avila, B. Geng, T. Zentgraf, L. Ju, F. Wang, and X. Zhang, “A graphene-based broadband optical modulator,” Nature 474(7349), 64–67 (2011). [CrossRef] [PubMed]
6. L. Ju, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, H. A. Bechtel, X. Liang, A. Zettl, Y. R. Shen, and F. Wang, “Graphene plasmonics for tunable terahertz metamaterials,” Nat. Nanotechnol. 6(10), 630–634 (2011). [CrossRef] [PubMed]
7. A. Vakil and N. Engheta, “Transformation optics using graphene,” Science 332(6035), 1291–1294 (2011). [CrossRef] [PubMed]
8. A. N. Grigorenko, M. Polini, and K. S. Novoselov, “Graphene plasmonics,” Nat. Photonics 6(11), 749–758 (2012). [CrossRef]
9. Y. D. Kim, H. Kim, Y. Cho, J. H. Ryoo, C. H. Park, P. Kim, Y. S. Kim, S. Lee, Y. Li, S. N. Park, Y. S. Yoo, D. Yoon, V. E. Dorgan, E. Pop, T. F. Heinz, J. Hone, S. H. Chun, H. Cheong, S. W. Lee, M. H. Bae, and Y. D. Park, “Bright visible light emission from graphene,” Nat. Nanotechnol. 10(8), 676–681 (2015). [CrossRef] [PubMed]
10. S. G. Liu, C. Zhang, M. Hu, X. X. Chen, P. Zhang, S. Gong, T. Zhao, and R. B. Zhong, “Coherent and tunable terahertz radiation from graphene surface plasmon polaritons excited by an electron beam,” Appl. Phys. Lett. 104(20), 201104 (2014). [CrossRef]
11. S. Gong, T. Zhao, M. Sanderson, M. Hu, R. B. Zhong, X. X. Chen, P. Zhang, C. Zhang, and S. G. Liu, “Transformation of surface plasmon polaritons to radiation in graphene in terahertz regime,” Appl. Phys. Lett. 106(22), 223107 (2015). [CrossRef]
12. K. J. A. Ooi, W. S. Koh, H. S. Chu, D. T. H. Tan, and L. K. Ang, “Efficiencies of Aloof-Scattered Electron Beam Excitation of Metal and Graphene Plasmons,” IEEE Trans. Plasma Sci. 43(4), 951–956 (2015). [CrossRef]
13. F. J. Garcıía de Abajo, “Multiple excitation of confined graphene plasmons by single free electrons,” ACS Nano 7(12), 11409–11419 (2013). [CrossRef] [PubMed]
14. K. J. A. Ooi, Y. S. Ang, J. L. Cheng, L. K. Ang, and D. T. H. Tan, “Electronic scattering of graphene plasmons in the nonlinear regime,” IEEE J. Sel. Top. Quantum Electron. 27(4), 5100206 (2017).
15. S. Liu, P. Zhang, W. Liu, S. Gong, R. Zhong, Y. Zhang, and M. Hu, “Surface Polariton Cherenkov Light Radiation Source,” Phys. Rev. Lett. 109(15), 153902 (2012). [CrossRef] [PubMed]
16. L. J. Wong, I. Kaminer, O. Ilic, J. D. Joannopoulos, and M. Soljačić, “Towards graphene plasmon-based free-electron infrared to X-ray sources,” Nat. Photonics 10(1), 46–52 (2015). [CrossRef]
17. K. J. A. Ooi, H. S. Chu, C. Y. Hsieh, D. T. H. Tan, and L. K. Ang, “Highly efficient mid-infrared on-chip electrical generation of graphene plasmons by inelastic electron tunnelling excitation,” Phys. Rev. Appl. 3(5), 054001 (2015). [CrossRef]
18. V. Ginsburg and I. Frank, “Radiation of a uniformly moving electron due to its transition from one medium into another,” J. Phys. 9, 353 (1945).
19. G. Q. Liao, Y. T. Li, Y. H. Zhang, H. Liu, X. L. Ge, S. Yang, W. Q. Wei, X. H. Yuan, Y. Q. Deng, B. J. Zhu, Z. Zhang, W. M. Wang, Z. M. Sheng, L. M. Chen, X. Lu, J. L. Ma, X. Wang, and J. Zhang, “Demonstration of Coherent Terahertz Transition Radiation from Relativistic Laser-Solid Interactions,” Phys. Rev. Lett. 116(20), 205003 (2016). [CrossRef] [PubMed]
20. M. Jablan, H. Buljan, and M. Soljacic, “Plasmonics in Graphene at Infrared Frequencies,” Phys. Rev. B 80(24), 245435 (2009). [CrossRef]
21. A. J. Fitzgerald, V. P. Wallace, M. Jimenez-Linan, L. Bobrow, R. J. Pye, A. D. Purushotham, and D. D. Arnone, “Terahertz Pulsed Imaging of Human Breast Tumors,” Radiology 239(2), 533–540 (2006). [CrossRef] [PubMed]
22. P. Tewari, Z. D. Taylor, D. Bennett, R. S. Singh, M. O. Culjat, C. P. Kealey, J. P. Hubschman, S. White, A. Cochran, E. R. Brown, and W. S. Grundfest, “Terahertz imaging of biological tissues,” Stud. Health Technol. Inform. 163, 653–657 (2011). [PubMed]
23. R. H. Clothier and N. Bourne, “Effects of THz exposure on human primary keratinocyte differentiation and viability,” J. Biol. Phys. 29(2-3), 179–185 (2003). [CrossRef] [PubMed]
24. G. J. Wilmink, B. D. Rivest, C. C. Roth, B. L. Ibey, J. A. Payne, L. X. Cundin, J. E. Grundt, X. Peralta, D. G. Mixon, and W. P. Roach, “In vitro investigation of the biological effects associated with human dermal fibroblasts exposed to 2.52 THz radiation,” Lasers Surg. Med. 43(2), 152–163 (2011). [CrossRef] [PubMed]