Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Exploiting black phosphorus based-Tamm plasmons in the terahertz region

Open Access Open Access

Abstract

Polarization-sensitive Tamm plasmons are investigated in a multi-layer photonic configuration where a monolayer black phosphorus (BP) is coated on a Bragg mirror separated by a dielectric. Owing to the in-plane anisotropy of BP, the Tamm plasmon can be excited selectively by tuning the BP carrier density. Cross-polarization conversion occurs when the armchair direction of BP makes an angle with the incident plan, i.e., ϕ≠0 or 90°. The BP-based Tamm device can be used as an intensity modulator with a modulation depth up to ∼100% and an insertion loss smaller than -0.55 dB. By analyzing the polarization evolution carefully, a multichannel polarization division multiplexing scheme is proposed and discussed. These findings open a new avenue for exploiting versatile tunable THz devices based on the monolayer of BP.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Polarization is an intrinsic property of electromagnetic waves. Active control and modulation of the polarization state of electromagnetic wave is of practical importance in many areas such as imaging, sensing and communications [13]. Many active materials with different tuning mechanisms have been proposed including electrically-tuned liquid crystals [4], mechanically-deformed elastic materials [5], and electrically-modulated graphene metamaterials [6,7]. Among them, graphene metamaterials show advantages owing to the its atomic thickness, high electrical modulation speed, and wide operation spectrum from sub-terahertz to near-infrared [8]. In the THz region, graphene metamaterials with subwavelength patterned structure have been designed to functioned as quarter wave plate [9], asymmetric transmission [6,10], polarization converter [11], etc. However, the electrical modulation of the patterned graphene is difficult. To address this problem, researchers assemble a complete large area graphene with anisotropic/chiral subwavelength structure made of metal or dielectric [12,13]. In 2017, Kim and co-workers combined a chiral metamaterial with a gated single-layer sheet of graphene [13]. They demonstrated that transmission of a terahertz wave with one circular polarization can be electrically controlled with large-intensity modulation depths (>99%) [13].

Black phosphorus (BP) has recently emerged as a strong competitor to graphene, due to its high carrier mobility and uniquely tunable band gap ranging from 0.3 to 2 eV [1416]. The striking in-plane anisotropy of BP promises novel polarization-dependent and angle-resolved optoelectronic devices without the need of patterned structure or additional subwavelength metal/dielectric structure [17,18]. In 2016, the anisotropic property of BP has been visualized by a polarized optical microscopy [19]. The crystallographic orientation of BP can be precisely determined by a polarized reflectance measurement [17,1921]. In 2017, Yang et al. compared the optical birefringence of different anisotropic 2D layered materials [22]. They found that birefringence in BP (∼0.245) is larger than those of anisotropic transition metal dichalcogenides such as ReS2 (∼0.037) and ReSe2 (∼0.047) and is larger than that of the current state of the art bulk materials (e.g., CaCO3) [22]. The layered BP can rotate polarization state with of ∼0.005° per atomic layer [22]. BP has found applications ranging from ultrafast photonic absorbers [23], multidisciplinary biomedical applications [24], vector soliton fiber lasers [25], to all-optical signal processing [26]. Recently, BP has expanded its applications to the THz region. In 2020, Guo et al. fabricated a BP-based THz photon detector with excellent sensitivity of 297 VW−1 and low noise equivalent power of 58 pW/Hz0.5, which can be used for THz imaging [27]. Taking the advantage of the anisotropic property of BP, unusual polarization-dependent perfect absorption has been demonstrated [28]. Using BP as saturable absorber, a switchable dual-wavelength Q-switched fiber laser was fabricated, which might be applied to effectively generate THz signal [29]. However, the instability of BP hinders the application of black phosphorus, and a breakthrough has been made in this field in recent years. Some recent work has proved how to improve the quality of monolayer black phosphorus to solve the problem of oxidation and instability [26].

Although the BP has strong anisotropy, the monolayer of BP is far from being able to manipulate the polarization state efficiently. Here, a BP-based Tamm (BPT) structure is designed with a monolayer BP siting on the top of the structure. The excitation of the Tamm plasmon enhances the interaction between BP and electromagnetic field, thus the polarization state of the reflected field can be modulated flexibly. Tamm plasmons are electromagnetic modes confined between a distributed Bragg reflector (DBR) and a metal (e.g. gold in the visible region, graphene in the THz region) [30]. Comparing to the surface-plasmons, the Tamm plasmons may be excited by both the TE and TM modes and present a relatively high quality factor [31]. Tamm plasmons therefore finds applications in refractive-index sensing [32,33], perfect absorption [34], and thermal emission [35].

2. Model and theory

As shown in Fig. 1(a), the BPT structure consists of a monolayer of BP, a defect layer, and a distributed Bragg reflector (DBR). BP shows metallic property in the THz region. The DBR is composed of 15 alternating dielectric layers of silicon (Si) and silicon dioxide (SiO2). The refractive indices of Si and SiO2 are respectively 3.42 and 1.9 with quarter-wave layer thickness for a center frequency of 1 THz, namely, dSi=21.93 µm, dSiO2=39.47 µm. A Si defect layer of thickness ds is sandwiched by BP and DBR, which can tune the resonance dip.

 figure: Fig. 1.

Fig. 1. (a) Schematic diagram of BP-based Tamm structure, where ds is the thickness of the Si defect layer between BP and the DBR. (b) Modulus of reflection coefficients |rBP| (for a Si-air interface containing BP), |rDBR| (for a Si-DBR interface), and |rpp| (for THz wave incident normally onto a BP-based Tamm structure) and phase ${\mathop{\rm Arg}\nolimits} [{r_{BP}}{r_{DBR}}\exp (2i\delta )]$ changing with frequency. (c) A single BP layer surrounding by two dielectrics with refractive index of nj and nj+1. (d) A stack of dielectrics containing a BP layer. (e) A dissymmetric cavity, enclosed by a sheet of BP and a DBR.

Download Full Size | PDF

Let us consider a monochromatic light beam illuminating onto BP-based Tamm structure. The photonic property of BP can be described by a surface conductivity [33]. Due to the in-plane anisotropic property, the conductivity is different along the armchair or zigzag crystalline directions. Under the Drude model, it is [36,37]

$${\sigma _{arm,zig}} = (i{D_{arm,zig}})/[\pi (\omega + i\eta /\hbar )],$$
where Darm,zig=πe2n/marm,zig is the Drude weight. marm=ћ2/(2γ2c), mzig=ћ2/2vc are respectively the electron mass along the armchair and zigzag directions, where η=10meV, Δ=2 eV, ηc2/0.4m0, vc= ћ2/1.4m0, γ=4a/π eVm, with the scale length of BP being a = 0.223 nm. The electron carrier density ρ can be changed by electric doping via a bias voltage [38,39].

By rotating the BP, the armchair axis of BP crystal makes an angle of ϕ with the incident plane. Therefore, the conductance matrix connecting the surface current and electric light field can be given by σ=[σpp, σps; σsp, σss], where σpp=σarmcos2ϕ+σzigsin2ϕ, σss=σarmsin2ϕ+σzigcos2ϕ, and σsp=σps=(σzig-σarm)sinϕcosϕ [15,34]. The cross conductivity σsp vanishes for isotropic 2D materials such as graphene [40]. For BP, however, σsp is generally nonzero except when the BP angle ϕ=0 and 90°.

Modified transfer matrix has been derived in our previous work to describe the light propagation through layered media with BP layers [41]. In this work, a 4×4 matrix is employed since the cross conductivity of BP will induce a coupling between p and s waves. According to Ref. [41], the transmission matrix for a dielectric interface containing BP layer (see Fig. 1(b)) is

$$T_{j,j + 1}^{BP} = \left[ {\begin{array}{cccc} {{k_{j + 1}}/{k_j} + {\eta_p} + {\xi_p}{\sigma_{xx}}}&{{k_{j + 1}}/{k_j} - {\eta_p} - {\xi_p}{\sigma_{xx}}}&{{\zeta_p}{\sigma_{xy}}}&{{\zeta_p}{\sigma_{xy}}}\\ {{k_{j + 1}}/{k_j} - {\eta_p} + {\xi_p}{\sigma_{xx}}}&{{k_{j + 1}}/{k_j} + {\eta_p} + {\xi_p}{\sigma_{xx}}}&{{\zeta_p}{\sigma_{xy}}}&{{\zeta_p}{\sigma_{xy}}}\\ {{\zeta_s}{\sigma_{yx}}}&{ - {\zeta_s}{\sigma_{yx}}}&{1 + {\eta_s} + {\xi_s}{\sigma_{yy}}}&{1 - {\eta_s} + {\xi_s}{\sigma_{yy}}}\\ { - {\zeta_s}{\sigma_{yx}}}&{{\zeta_s}{\sigma_{yx}}}&{1 - {\eta_s} - {\xi_s}{\sigma_{yy}}}&{1 + {\eta_s} - {\xi_s}{\sigma_{yy}}} \end{array}} \right],$$
where ${\eta _p} = {k_j}{k_{j + 1,z}}/{k_{j + 1}}{k_{jz}}$, ${\eta _s} = {k_{j + 1,z}}/{k_{jz}}$, ${\xi _p} = \omega \mu {k_{j + 1,z}}/{k_{j + 1}}{k_j}$, ${\xi _p} = \omega \mu /{k_{jz}}$, ${\zeta _p} = \omega \mu /{k_j}$, ${\zeta _s} = \omega \mu {k_{j + 1,z}}/{k_{jz}}{k_{j + 1}}$, µ are circular frequency and permeability. kl (l = j,j + 1) is the wavenumber of dielectric l, and klz is their z component. Without the BP layer, the upper-right and lower-left parts of the transmission matrix vanish, the upper-left and lower-right 2×2 matrices describe the transmission of p and s waves, independently. Thus, the 4×4 transmission matrix can be rewritten into two conventional 2×2 matrices.

The effect of BP layer is excluded in the propagation matrix as the thicknesses of BP layers are neglected. The propagation matrix in a homogenous medium j with a thickness of dj is

$${P_j} = \left[ {\begin{array}{cccc} {{e^{ - i{k_{jz}}{d_j}}}}&0&0&0\\ 0&{{e^{i{k_{jz}}{d_j}}}}&0&0\\ 0&0&{{e^{ - i{k_{jz}}{d_j}}}}&0\\ 0&0&0&{{e^{i{k_{jz}}{d_j}}}} \end{array}} \right].$$
As shown by Fig. 1(c), the BP-based Tamm structure is composing of a BP layer, a defect layer, and a DBR. Denote the electric field coefficients on left side of the first interface by A0, B0, C0, D0 and those on the right side of the end interface by AN+1, BN+1, CN+1, DN+1. These two sets of field coefficients are then related by a 4×4 transfer matrix M, namely,
$$\left[ {\begin{array}{c} {{A_0}}\\ {{B_0}}\\ {{C_0}}\\ {{D_0}} \end{array}} \right] = M\left[ {\begin{array}{c} {{A_{N + 1}}}\\ {{B_{N + 1}}}\\ {{C_{N + 1}}}\\ {{D_{N + 1}}} \end{array}} \right],$$
with $M = T_{0,s}^{BP}{P_s}{T_{a,b}}{P_b}{T_{b,a}}{P_a}\ldots {T_{b,a}}{P_a}{T_{b,a}}$. $T_{0,s}^{BP}$ describes an air-defect layer interface containing BP layer, while ${T_{a,b}}$ and ${T_{b,a}}$ describe the dielectric a-b and dielectric b-a interface, respectively. Ps, Pa, Pb are the propagation matrix in defect layer and media a and b, respectively.

According to Eq. (4), the Fresnel reflection coefficients can be obtained. By setting C0=BN+1=DN+1=0, we have

$${r_{pp}} = \frac{{{B_0}}}{{{A_0}}} = \frac{{{M_{21}}{M_{33}} - {M_{23}}{M_{31}}}}{{{M_{11}}{M_{33}} - {M_{13}}{M_{31}}}},$$
$${r_{sp}} = \frac{{{D_0}}}{{{A_0}}} = \frac{{{M_{41}}{M_{33}} - {M_{43}}{M_{31}}}}{{{M_{11}}{M_{33}} - {M_{13}}{M_{31}}}}.$$
Similarly, by setting A0=BN+1=DN+1=0, we obtain
$${r_{ss}} = \frac{{{D_0}}}{{{C_0}}} = \frac{{{M_{43}}{M_{11}} - {M_{41}}{M_{13}}}}{{{M_{11}}{M_{33}} - {M_{13}}{M_{31}}}},$$
$${r_{ps}} = \frac{{{A_0}}}{{{C_0}}} = \frac{{{M_{23}}{M_{11}} - {M_{21}}{M_{13}}}}{{{M_{11}}{M_{33}} - {M_{13}}{M_{31}}}}.$$
It should be noted that rps=rsp. In the same way, the transfer matrix describing the electromagnetic wave incident from the defect layer onto the BP layer and DBR can be constructed i.e., $M = T_{\textrm{s},0}^{BP}$ and $M = {T_{a,b}}{P_b}{T_{b,a}}{P_a}\ldots {T_{b,a}}{P_a}{T_{b,a}}$, respectively. Thus, the corresponding reflection coefficients rBP and rDBR can be obatined, as shown in Fig. 1(e).

The BPT structure can be regarded as a dissymmetric cavity, enclosed by a sheet of BP on one side and a DBR on the other [34]. For the case of ϕ=0 and 90°, the resonance condition for the cavity with incident p (TE) and s (TM) plane wave can be given respectively by [31,32]

$$r_{\textrm{BP}}^{pp,ss}{r_{DBR}}\exp (2i\delta ) \approx 1.$$
δ is the propagation phase within the deflection layer. Since $r_{\textrm{BP}}^{pp}$ and $r_{\textrm{BP}}^{\textrm{s}p}$ are different, the excitation conditions of Tamm plasmons for TE and TM modes are different, which is different from the Tamm plasmons based on graphene and metals. When ϕ≠0 or 90°, the TE and TM plane waves couple each other. In this situation, however, the cross polarization conversion effect can be enhanced.

3. Results and discussions

Figure 1(b) shows the modulus and phase of reflection coefficients, when a p-polarized 1 THz wave is illuminated normally onto the BPT structure with BP angle ϕ=0, i.e., the polarization direction is along the armchair axis of BP. The characteristic stop band of DBR is clearly observed in the red line of Fig. 1(b), which is centered at 1 THz with a bandwidth of ∼0.36 THz. The reflection coefficient for a Si–air interface containing BP layer is shown in light yellow line, with a relatively constant modulus of ∼0.633 over the spectral range. This relatively small modulus results in low quality factor of asymmetric cavity [42]. The green curve presents for comparison the full reflectance of BPT structure |rpp|, the Tamm plasmon appears as a dip at 1 THz. The resonance condition expressed in Eq. (2) requires ${\mathop{\rm Arg}\nolimits} [{r_{BP}}{r_{DBR}}\exp (2i\delta )] = 0$. As shown by the green line in Fig. 1(b), this zero-crossing point is observed at 1 THz, and coincides with the reflectance dip for |rpp|.

The Bragg structure acts as a perfect reflector. Si-SiO2 period numbers will affect the performance of Bragg structure. We have calculated the reflection coefficients of |rDBR| and |rpp| for different period numbers. As shown by Fig. 2, the |rDBR| increases with the period number, resulting in the decrease of the dip of |rpp|. Although the side lobes of |rpp| vary with the period number T, the main lobes are almost identical when T > 7. Throughout the paper, the period number is fixed to T = 15.

 figure: Fig. 2.

Fig. 2. The reflection coefficients of |rDBR| (a) and |rpp| (b) vs frequency for different Si-SiO2 period T in Bragg structure.

Download Full Size | PDF

The reflectance dip varies with the thickness of defect layer ds. The reflection coefficients |rpp| and |rss| as functions of frequency f at normal incidence for ds = 130, 133 and 136 µm are shown respectively in Fig. 3(a), where the vanishing |rps| is not shown since the BP angle is ϕ=0. The dip positions for |rpp| and |rss| have a little different owing to the anisotropic of BP layer. The dip position for |rpp| moves from 0.98 to 1.02 THz when ds increases from 130 to 136 µm. Further increase of ds can produce multiple reflection dips within the stop band between 0.81 ∼ 1.18THz (see Fig. 6(b) and 6(c)).

 figure: Fig. 3.

Fig. 3. (a) Reflection coefficients as functions of frequency f at normal incidence for different ds (a) and different carrier density (b). In (a), θ=0, ϕ=0°; θ=30°, ϕ=45° in (b).

Download Full Size | PDF

When the BP angle is set as ϕ=45°, the p and s plane waves couple each other, resulting in nonzero non-diagonal reflection coefficients, rps=rsp≠0. As shown by Fig. 3(b), the |rps| is up to 0.24 at f = 1 THz, when ρ=8.55×1017 m-2. At this carrier density, |rpp| has a dip at 1 THz, which is down to 4.55 × 10−3, while the dip of |rss| is 0.148. However, without BP, ρ=0, there is not reflectance dip within the stopband, and rps vanishes.

To achieve a change in carrier density, it is usually necessary to apply a bias voltage. To this end, a several tenths of nanometer (∼10−4 wavelength) insulating layer is inserted between the BP and top silicon layers, as shown by Fig. 4(a). The inserted insulating layer will hardly affect the Fresnel reflection coefficients of the BP-based Tamm structure. According to Ref. [18,43], the carrier density of BP can be electrically modulated from a minimum of 0.17×1017 m-2 to a maximum of 12×1017 m-2.

 figure: Fig. 4.

Fig. 4. (a) The schematic of modulation models for reflection coefficients. (b) The Fresnel reflection coefficient as function of carrier density ρ. (c,d) The reflection coefficient as function of frequency for ρ=4.68×1017 m-2 (blue star in (b)) and 11.68×1017 m-2 (red srar in (b)).

Download Full Size | PDF

A scheme for THz wave modulation is shown by Fig. 4(a), where both the incidence angle and BP angle are zero. By adjusting the BP carrier density, the amplitude changes dominate over the phase changes in the reflection coefficients because the monolayer BP is too thin for effective phase accumulation. When the carrier density ρ increases from 0, both |rpp| and |rss| decreases from 1 gradually. |rpp| reaches a minimum of 1.35×10−3 when ρ=4.71×1017 m-2, while the minimum of |rss| is 1.28×10−3 at ρ=11.73×1017 m-2, as shown by Fig. 4(b). Figure 4(c) shows the reflection coefficient as function of frequency for ρ=4.68×1017 m-2 (blue star in (b)) and 11.68×1017 m-2 (red star in (b)). One can conclude that, the Tamm plasmons for TE and TM wave can be selectively excited by controlling the carrier density of BP. This unique phenomenon originates from the anisotropic property of BP, cannot be found in isotropic materials such as metals and graphene. For isotropic materials, the reflection of TE and TM waves are identical at normal incidence. The intensity modulation depth for TE and TM waves can be defined as [44]

$$M{D_{p,s}} = [1 - \frac{{|{r_{pp,ss}}{|^2}}}{{|r_{pp,ss}^{\max }{|^2}}}] \times 100\%,$$
where $r_{pp,ss}^{\max }$ are the baselines defined as the maximum reachable reflection values of ${r_{pp,ss}}$. Assume that 0.17×1017 m-2 is the minimum carrier density of BP [18,43]. Thus, the maximum reachable reflection values of ${r_{pp,ss}}$ are 0.939 and 0.975, respectively. The intensity modulation depth for TE and TM waves are both 100%, with the insertion loss 10lg$\textrm{|}r_{pp,ss}^{\max }{\textrm{|}^\textrm{2}}$ of -0.55dB, -0.22dB, respectively. Comparing to the graphene-based THz modulator demonstrated by Chen et al. [44], our BP-based device can operate for both TE and TM modes with much lower insertion loss.

Figure 5(a) shows the reflection coefficients when THz wave illuminated obliquely with θ=60°. As shown by Fig. 5(a), the |rss| can be effectively modulated by the BP carrier density effectively. However, the |rpp| can be hardly changed by the carrier density. |rpp| increases a litter from 0.3 to 0.33, when ρ increases from 0 to 6×1017 m-2. In Fig. 5(b), |rps| is no longer zero, when the BP angle is 45°. The incident angle and thickness of defected layer are optimized respectively to be θ=60° and ds=354.9 µm to maximize the polarization rotation angle. |rps| increases gradually with ρ, and approach to an asymptotic value of 0.10. Therefore, the reflected THz wave can be flexibly modulated by tuning the BP carrier density.

 figure: Fig. 5.

Fig. 5. (a) Reflection coefficients (a,b) and polarization angles (c,d) as functions of BP carrier density for (a,c) θ=60°, ϕ=90°; (b,d) θ=60°, ϕ=45°.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. (a,b) Reflection coefficients |rps| (b) and |rss| (c) as functions of frequency and carrier density. (b) Stokes parameters of the reflected beam changing with ρ; (c) The evolution of polarization states in Poincaré sphere. The red, blue and green lines in (b,c) represent f = 0.83, 1.00, and 1.17 THz, respectively.

Download Full Size | PDF

In the following, we focus on attention on the polarization modulation by the carrier density. For the case of ϕ=90°, we assume that the incident polarization is diagonal polarization, i.e. [1,1]T. The polarization of the reflected wave is therefore [rpp,rss]T. As shown by Fig. 5(c), the reflected wave is in elliptical polarization, which varies with the carrier density. The polarization angle is defined as

$$\varphi = 0.5\arctan [S2/S1].$$
S1, S2, and S3 are the Stokes parameters of the reflected waves. When the carrier density increases from 0 to 6×1017 m-2, the polarization angle changes from 75° to -58.5° gradually.

In the case of ϕ=45°, the incident wave is assumed to be TE mode, namely, [0,1]T. The polarization of the reflected wave is [rps,rss]T. In this case, the reflected wave is linearly polarized. The polarization angle, increases from -90° with the carrier densy and approaches gradually to 81° with an angle change up to 171°. The insets show four special polarization states at ρ=0, 0.90×1017, 1.01×1017, 1.15×1017 and 2×1017 m-2, corresponding polarization angle corresponding polarization angles are -90°, -45°, 0°, 45°, 76°, respectively.

Multi-channel modulations can be achieved by increasing the thickness of the defect layer ds. Figures 5(a) and (b) shows the reflection coefficients for ϕ=0 at normal incidence, when ds=220.8 µm. As shown by Fig. 6(a), three dips can be found with the frequency in |rpp|, which are respectively f = 0.83, 1.00 and 1.17 THz at ρ=4.7×1017m-2. If the incident beam is circularly polarized, i.e., [1,i]T. The polarization state of the reflected wave is [rpp,irss]T. The Stokes parameters of the reflected beam S1, S2, and S3 can be calculated numerically, and the results are shown in Fig. 6(c). With the increase of the carrier density, parameter S3 increases gradually from -1 to + 1, and then decreases back to ∼0.5. Parameter |S2| is smaller than 0.2 for an arbitrary density. The evolution of polarization state on Poincaré sphere with the change of carrier density are shown in Fig. 6(d). The curves for three channels are almost overlapped. Therefore, one can control the polarizations of the reflected waves at three different frequency channels simultaneously.

As shown by Fig. 6(d), the reflected wave is in horizontal [1,0]T and vertical [0,1]T polarization states when ρ=4.71×1017 and 11.73×1017 m-2, respectively. When ρ=7.44×1017m-2, the polarization state is [1,ei(-π/4)]T. This elliptical polarization state can be considered as a superposition of horizontal and vertical polarization states with equal amplitude but a phase difference.

As demonstrated above, the BPT device can flexibly modulate the polarization state of the reflected light into two orthogonal polarization states as well as their superposition state by tuning the BP carrier density via bias voltage. Therefore, the proposed device can be employed for polarization-division multiplexing (PDM). PDM is a crucial technique that can significantly increase the transmission capacity of a single physical channel [13].

Figure 7 shows the multi-channel PDM scheme based on the BP-based Tamm structure. Three frequency channels are modulated simultaneously by the carrier density. As shown in Fig. 7, two independent information signals can be, respectively, encoded into two orthogonally polarized light beams with a binary format. Three different carrier densities of ρ1, ρ2, and ρ3 correspond to the horizontal, elliptical, and vertical polarization states, which associate to “10”, “11”, “01”, respectively. By tuning the carrier density, the reflected beams with a particular sequence of polarization states can be generated directly from the proposed structure without the need to couple multiple sources.

 figure: Fig. 7.

Fig. 7. The schematic of the multi-channel PDM technology, according to the superposition principle of the polarization state.

Download Full Size | PDF

The proposed BPT device can be fabricated by stacking the 21.93 µm Si and 39.47 µm SiO2 slices alternatively [45]. And then the Si defected layer with thickness ds>100 µm is stacked on the top. Before transferring BP layer, a nanoscale thickness insulator layer (SiO2 or Al2O3) is deposited on the defect layer by atomic layer deposition. Finally, an Au electrode is thermally deposited on the BP layer. A voltage is applied via the Si defected layer and Au electrode to tune the carrier density of BP. The Bragg reflector can be also made of alternating air spacers and SiO2 layers, which is simpler in fabrication [45].

Before experimental measurement, it is necessary to determine crystalline axes of BP. As illustrated by Ref. [21], the armchair and zigzag axes of BP are associated respectively with maximum and minimum brightness of R channel, when analyzing the reflection of normal incident visible light with a polarizer and a CCD.

4. Conclusions

A novel multichannel THz polarization modulation scheme was proposed based on the BP-based Tamm structure. The excitation of Tamm plasmons enables an efficient and flexible modulation of THz waves. When the BP angle ϕ=0 or 90°, a perfect intensity modulator for TM and/or TE modes with a modulation depth up to ∼100% and an insertion loss smaller than -0.55 dB is demonstrated. When ϕ=45°, the polarization of the reflected THz wave rotates when the BP carrier density changes. The maximum polarization rotation angle is up to 171°. The polarization state changing with BP carrier density are visually shown in Poincare sphere for different frequency channels, based on which the multi-channel PDM is achieved. We believe the proposed versatile BPT device has potential applications in THz imaging, sensing and communications.

Funding

Science & Technology Project of Guangzhou (201604040005, 201605030002, 201704030105, 201707010396); Natural Science Foundation of Guangdong Province (2016A030311019, 2016A030313079, 2016TQ03X962, 2017A010102006, 2017A030313375); National Natural Science Foundation of China (61475066, 61505069, 61705086); College Students' Innovative Entrepreneurial Training Plan Program (201910559054).

Disclosures

The authors declare no conflicts of interest.

References

1. L. Ju, B. S. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, H. A. Bechtel, X. G. Liang, A. Zettl, Y. R. Shen, and F. Wang, “Graphene plasmonics for tunable terahertz metamaterials,” Nat. Nanotechnol. 6(10), 630–634 (2011). [CrossRef]  

2. N. K. Grady, J. E. Heyes, D. R. Chowdhury, Y. Zeng, M. T. Reiten, A. K. Azad, A. J. Taylor, D. A. R. Dalvit, and H.-T. Chen, “Terahertz Metamaterials for Linear Polarization Conversion and Anomalous Refraction,” Science 340(6138), 1304–1307 (2013). [CrossRef]  

3. M. Tamagnone, C. Moldovan, J.-M. Poumirol, A. B. Kuzmenko, A. M. Ionescu, J. R. Mosig, and J. Perruisseau-Carrier, “Near optimal graphene terahertz non-reciprocal isolator,” Nat. Commun. 7(1), 11216 (2016). [CrossRef]  

4. L. Wang, X.-W. Lin, W. Hu, G.-H. Shao, P. Chen, L.-J. Liang, B.-B. Jin, P.-H. Wu, H. Qian, Y.-N. Lu, X. Liang, Z.-G. Zheng, and Y.-Q. Lu, “Broadband tunable liquid crystal terahertz waveplates driven with porous graphene electrodes,” Light: Sci. Appl. 4(2), e253 (2015). [CrossRef]  

5. L. Cong, P. Pitchappa, Y. Wu, L. Ke, C. Lee, N. Singh, H. Yang, and R. Singh, “Active Multifunctional Microelectromechanical System Metadevices: Applications in Polarization Control, Wavefront Deflection, and Holograms,” Adv. Opt. Mater. 5(2), 1600716 (2017). [CrossRef]  

6. Y. Huang, Z. Yao, F. Hu, C. Liu, L. Yu, Y. Jin, and X. Xu, “Tunable circular polarization conversion and asymmetric transmission of planar chiral graphene-metamaterial in terahertz region,” Carbon 119, 305–313 (2017). [CrossRef]  

7. P. Q. Liu, I. J. Luxmoore, S. A. Mikhailov, N. A. Savostianova, F. Valmorra, J. Faist, and G. R. Nash, “Highly tunable hybrid metamaterials employing split-ring resonators strongly coupled to graphene surface plasmons,” Nat. Commun. 6(1), 8969 (2015). [CrossRef]  

8. J. Cheng, F. Fan, and S. Chang, “Recent Progress on Graphene-Functionalized Metasurfaces for Tunable Phase and Polarization Control,” Nanomaterials 9(3), 398 (2019). [CrossRef]  

9. Y. Zhang, Y. Feng, B. Zhu, J. Zhao, and T. Jiang, “Switchable quarter-wave plate with graphene based metamaterial for broadband terahertz wave manipulation,” Opt. Express 23(21), 27230–27239 (2015). [CrossRef]  

10. J. Zhao, J. Song, T. Xu, T. Yang, and J. Zhou, “Controllable linear asymmetric transmission and perfect polarization conversion in a terahertz hybrid metal-graphene metasurface,” Opt. Express 27(7), 9773–9781 (2019). [CrossRef]  

11. Y. Zhang, Y. Feng, T. Jiang, J. Cao, J. Zhao, and B. Zhu, “Tunable broadband polarization rotator in terahertz frequency based on graphene metamaterial,” Carbon 133, 170–175 (2018). [CrossRef]  

12. T.-T. Kim, S. S. Oh, H.-D. Kim, H. S. Park, O. Hess, B. Min, and S. Zhang, “Electrical access to critical coupling of circularly polarized waves in graphene chiral metamaterials,” Sci. Adv. 3(9), e1701377 (2017). [CrossRef]  .

13. J. Li, P. Yu, H. Cheng, W. Liu, Z. Li, B. Xie, S. Chen, and J. Tian, “Optical Polarization Encoding Using Graphene-Loaded Plasmonic Metasurfaces,” Adv. Opt. Mater. 4(1), 91–98 (2016). [CrossRef]  

14. B. Deng, R. Frisenda, C. Li, X. Chen, A. Castellanos-Gomez, and F. Xia, “Progress on Black Phosphorus Photonics,” Adv. Opt. Mater. 6(19), 1800365 (2018). [CrossRef]  

15. H. Lin, B. Chen, S. Yang, W. Zhu, J. Yu, H. Guan, H. Lu, Y. Luo, and Z. Chen, “Photonic spin Hall effect of monolayer black phosphorus in the Terahertz region,” Nanophotonics 7(12), 1929–1937 (2018). [CrossRef]  

16. X. Wang and S. Lan, “Optical properties of black phosphorus,” Adv. Opt. Photonics 8(4), 618–655 (2016). [CrossRef]  

17. M. C. Sherrott, W. S. Whitney, D. Jariwala, S. Biswas, C. M. Went, J. Wong, G. R. Rossman, and H. A. Atwater, “Anisotropic Quantum Well Electro-Optics in Few-Layer Black Phosphorus,” Nano Lett. 19(1), 269–276 (2019). [CrossRef]  

18. H. Yuan, X. Liu, F. Afshinmanesh, W. Li, G. Xu, J. Sun, B. Lian, A. G. Curto, G. Ye, Y. Hikita, Z. Shen, S.-C. Zhang, X. Chen, M. Brongersma, H. Y. Hwang, and Y. Cui, “Polarization-sensitive broadband photodetector using a black phosphorus vertical p-n junction,” Nat. Nanotechnol. 10(8), 707–713 (2015). [CrossRef]  

19. S. Lan, S. Rodrigues, L. Kang, and W. Cai, “Visualizing Optical Phase Anisotropy in Black Phosphorus,” ACS Photonics 3(7), 1176–1181 (2016). [CrossRef]  

20. A. Islam, W. Du, V. Pashaei, H. Jia, Z. Wang, J. Lee, G. J. Ye, X. H. Chen, and P. X. L. Feng, “Discerning Black Phosphorus Crystal Orientation and Anisotropy by Polarized Reflectance Measurement,” ACS Appl. Mater. Interfaces 10(30), 25629–25637 (2018). [CrossRef]  

21. H. Jiang, H. Shi, X. Sun, and B. Gao, “Optical Anisotropy of Few-Layer Black Phosphorus Visualized by Scanning Polarization Modulation Microscopy,” ACS Photonics 5(6), 2509–2515 (2018). [CrossRef]  

22. H. Yang, H. Jussila, A. Autere, H.-P. Komsa, G. Ye, X. Chen, T. Hasan, and Z. Sun, “Optical Waveplates Based on Birefringence of Anisotropic Two-Dimensional Layered Materials,” ACS Photonics 4(12), 3023–3030 (2017). [CrossRef]  

23. M. Zhang, Q. Wu, F. Zhang, L. Chen, X. Jin, Y. Hu, Z. Zheng, and H. Zhang, “2D Black Phosphorus Saturable Absorbers for Ultrafast Photonics,” Adv. Opt. Mater. 7(1), 1800224 (2019). [CrossRef]  

24. M. Qiu, W. X. Ren, T. Jeong, M. Won, G. Y. Park, D. K. Sang, L.-P. Liu, H. Zhang, and J. S. Kim, “Omnipotent phosphorene: a next-generation, two-dimensional nanoplatform for multidisciplinary biomedical applications,” Chem. Soc. Rev. 47(15), 5588–5601 (2018). [CrossRef]  

25. Y. Song, S. Chen, Q. Zhang, L. Li, L. Zhao, H. Zhang, and D. Tang, “Vector soliton fiber laser passively mode locked by few layer black phosphorus-based optical saturable absorber,” Opt. Express 24(23), 25933–25942 (2016). [CrossRef]  

26. J. Zheng, Z. Yang, C. Si, Z. Liang, X. Chen, R. Cao, Z. Guo, K. Wang, Y. Zhang, J. Ji, M. Zhang, D. Fan, and H. Zhang, “Black Phosphorus Based All-Optical-Signal-Processing: Toward High Performances and Enhanced Stability,” ACS Photonics 4(6), 1466–1476 (2017). [CrossRef]  

27. W. Guo, Z. Dong, Y. Xu, C. Liu, D. Wei, L. Zhang, X. Shi, C. Guo, H. Xu, G. Chen, L. Wang, K. Zhang, X. Chen, and W. Lu, “Sensitive Terahertz Detection and Imaging Driven by the Photothermoelectric Effect in Ultrashort-Channel Black Phosphorus Devices,” Adv. Sci. 7(5), 1902699 (2020). [CrossRef]  

28. Y. M. Qing, H. F. Ma, and T. J. Cui, “Tailoring anisotropic perfect absorption in monolayer black phosphorus by critical coupling at terahertz frequencies,” Opt. Express 26(25), 32442–32450 (2018). [CrossRef]  

29. J. Liu, Y. Chen, Y. Li, H. Zhang, S. Zheng, and S. Xu, “Switchable dual-wavelength Q-switched fiber laser using multilayer black phosphorus as a saturable absorber,” Photonics Res. 6(3), 198–203 (2018). [CrossRef]  

30. M. Kaliteevski, I. Iorsh, S. Brand, R. A. Abram, J. M. Chamberlain, A. V. Kavokin, and I. A. Shelykh, “Tamm plasmon-polaritons: Possible electromagnetic states at the interface of a metal and a dielectric Bragg mirror,” Phys. Rev. B 76(16), 165415 (2007). [CrossRef]  

31. B. Auguie, A. Bruchhausen, and A. Fainstein, “Critical coupling to Tamm plasmons,” J. Opt. 17(3), 035003 (2015). [CrossRef]  

32. N. Li, T. Tang, J. Li, L. Luo, P. Sun, and J. Yao, “Highly sensitive sensors of fluid detection based on magneto-optical optical Tamm state,” Sens. Actuators, B 265, 644–651 (2018). [CrossRef]  

33. Y. Tsurimaki, J. K. Tong, V. N. Boriskin, A. Semenov, M. I. Ayzatsky, Y. P. Machekhin, G. Chen, and S. V. Boriskinalb, “Topological Engineering of Interfacial Optical Tamm States for Highly Sensitive Near-Singular-Phase Optical Detection,” ACS Photonics 5(3), 929–938 (2018). [CrossRef]  

34. X. Wang, X. Jiang, Q. You, J. Guo, X. Dai, and Y. Xiang, “Tunable and multichannel terahertz perfect absorber due to Tamm surface plasmons with graphene,” Photonics Res. 5(6), 536–542 (2017). [CrossRef]  

35. D. Geng, E. Cabello-Olmo, G. Lozano, and H. Miguez, “Tamm Plasmons Directionally Enhance Rare-Earth Nanophosphor Emission,” ACS Photonics 6(3), 634–641 (2019). [CrossRef]  

36. Z. Liu and K. Aydin, “Localized Surface Plasmons in Nanostructured Monolayer Black Phosphorus,” Nano Lett. 16(6), 3457–3462 (2016). [CrossRef]  

37. T. Low, R. Roldan, H. Wang, F. Xia, P. Avouris, L. Martin Moreno, and F. Guinea, “Plasmons and Screening in Monolayer and Multilayer Black Phosphorus,” Phys. Rev. Lett. 113(10), 106802 (2014). [CrossRef]  

38. A. Avsar, J. Y. Tan, M. Kurpas, M. Gmitra, K. Watanabe, T. Taniguchi, J. Fabian, and B. Ozyilmaz, “Gate-tunable black phosphorus spin valve with nanosecond spin lifetimes,” Nat. Phys. 13(9), 888–893 (2017). [CrossRef]  

39. M. Huang, M. Wang, C. Chen, Z. Ma, X. Li, J. Han, and Y. Wu, “Broadband Black-Phosphorus Photodetectors with High Responsivity,” Adv. Mater. 28(18), 3481–3485 (2016). [CrossRef]  

40. W. Zhu, M. Jiang, H. Guan, J. Yu, H. Lu, J. Zhang, and Z. Chen, “Tunable spin splitting of Laguerre-Gaussian beams in graphene metamaterials,” Photonics Res. 5(6), 684–688 (2017). [CrossRef]  

41. Y. Lin, X. Liu, H. Chen, X. Guo, J. Pan, J. Yu, H. Zheng, H. Guan, H. Lu, Y. Zhong, Y. Chen, Y. Luo, W. Zhu, and Z. Chen, “Tunable asymmetric spin splitting by black phosphorus sandwiched epsilon-near-zero-metamaterial in the terahertz region,” Opt. Express 27(11), 15868–15879 (2019). [CrossRef]  

42. H. M. Zhu, Y. P. Fu, F. Meng, X. X. Wu, Z. Z. Gong, Q. Ding, M. V. Gustafsson, M. T. Trinh, S. Jin, and X. Y. Zhu, “Lead halide perovskite nanowire lasers with low lasing thresholds and high quality factors,” Nat. Mater. 14(6), 636–642 (2015). [CrossRef]  

43. L. Li, F. Yang, G. J. Ye, Z. Zhang, Z. Zhu, W. Lou, X. Zhou, L. Li, K. Watanabe, T. Taniguchi, K. Chang, Y. Wang, X. H. Chen, and Y. Zhang, “Quantum Hall effect in black phosphorus two-dimensional electron system,” Nat. Nanotechnol. 11, 593–597 (2016). [CrossRef]  

44. Z. Chen, X. Chen, L. Tao, K. Chen, M. Long, X. Liu, K. Yan, R. I. Stantchev, E. Pickwell-MacPherson, and J.-B. Xu, “Graphene controlled Brewster angle device for ultra broadband terahertz modulation,” Nat. Commun. 9(1), 4909 (2018). [CrossRef]  

45. Y. Yu, J. H. Cai, J. D. Su, Z. P. Zhang, and H. Qin, “Fabrication and characterization of a wide-bandgap and high-Q terahertz distributed-Bragg-reflector micro cavities,” Opt. Commun. 426, 84–88 (2018). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. (a) Schematic diagram of BP-based Tamm structure, where ds is the thickness of the Si defect layer between BP and the DBR. (b) Modulus of reflection coefficients |rBP| (for a Si-air interface containing BP), |rDBR| (for a Si-DBR interface), and |rpp| (for THz wave incident normally onto a BP-based Tamm structure) and phase ${\mathop{\rm Arg}\nolimits} [{r_{BP}}{r_{DBR}}\exp (2i\delta )]$ changing with frequency. (c) A single BP layer surrounding by two dielectrics with refractive index of nj and nj+1. (d) A stack of dielectrics containing a BP layer. (e) A dissymmetric cavity, enclosed by a sheet of BP and a DBR.
Fig. 2.
Fig. 2. The reflection coefficients of |rDBR| (a) and |rpp| (b) vs frequency for different Si-SiO2 period T in Bragg structure.
Fig. 3.
Fig. 3. (a) Reflection coefficients as functions of frequency f at normal incidence for different ds (a) and different carrier density (b). In (a), θ=0, ϕ=0°; θ=30°, ϕ=45° in (b).
Fig. 4.
Fig. 4. (a) The schematic of modulation models for reflection coefficients. (b) The Fresnel reflection coefficient as function of carrier density ρ. (c,d) The reflection coefficient as function of frequency for ρ=4.68×1017 m-2 (blue star in (b)) and 11.68×1017 m-2 (red srar in (b)).
Fig. 5.
Fig. 5. (a) Reflection coefficients (a,b) and polarization angles (c,d) as functions of BP carrier density for (a,c) θ=60°, ϕ=90°; (b,d) θ=60°, ϕ=45°.
Fig. 6.
Fig. 6. (a,b) Reflection coefficients |rps| (b) and |rss| (c) as functions of frequency and carrier density. (b) Stokes parameters of the reflected beam changing with ρ; (c) The evolution of polarization states in Poincaré sphere. The red, blue and green lines in (b,c) represent f = 0.83, 1.00, and 1.17 THz, respectively.
Fig. 7.
Fig. 7. The schematic of the multi-channel PDM technology, according to the superposition principle of the polarization state.

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

σ a r m , z i g = ( i D a r m , z i g ) / [ π ( ω + i η / ) ] ,
T j , j + 1 B P = [ k j + 1 / k j + η p + ξ p σ x x k j + 1 / k j η p ξ p σ x x ζ p σ x y ζ p σ x y k j + 1 / k j η p + ξ p σ x x k j + 1 / k j + η p + ξ p σ x x ζ p σ x y ζ p σ x y ζ s σ y x ζ s σ y x 1 + η s + ξ s σ y y 1 η s + ξ s σ y y ζ s σ y x ζ s σ y x 1 η s ξ s σ y y 1 + η s ξ s σ y y ] ,
P j = [ e i k j z d j 0 0 0 0 e i k j z d j 0 0 0 0 e i k j z d j 0 0 0 0 e i k j z d j ] .
[ A 0 B 0 C 0 D 0 ] = M [ A N + 1 B N + 1 C N + 1 D N + 1 ] ,
r p p = B 0 A 0 = M 21 M 33 M 23 M 31 M 11 M 33 M 13 M 31 ,
r s p = D 0 A 0 = M 41 M 33 M 43 M 31 M 11 M 33 M 13 M 31 .
r s s = D 0 C 0 = M 43 M 11 M 41 M 13 M 11 M 33 M 13 M 31 ,
r p s = A 0 C 0 = M 23 M 11 M 21 M 13 M 11 M 33 M 13 M 31 .
r BP p p , s s r D B R exp ( 2 i δ ) 1.
M D p , s = [ 1 | r p p , s s | 2 | r p p , s s max | 2 ] × 100 % ,
φ = 0.5 arctan [ S 2 / S 1 ] .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.