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

Tuning the transmission lineshape of a photonic crystal slab guided-resonance mode by polarization control

Open Access Open Access

Abstract

We demonstrate a system consisting of a two-dimensional photonic crystal slab and two polarizers which has a tunable transmission lineshape. The lineshape can be tuned from a symmetric Lorentzian to a highly asymmetric Fano lineshape by rotating the output polarizer. We use temporal coupled mode theory to explain the measurement results. The theory also predicts tunable phase shift and group delay.

© 2013 Optical Society of America

1. Introduction

Controlling the lineshape of photonic resonators is of great interest for their applications in various fields [1]. For photonic crystal microcavities, lineshape tuning has been demonstrated by tuning the incident spot size [2] or polarization [3]. Similar effects have been investigated in plasmonic systems [4]. In this paper, we focus on guided resonance modes in photonic crystal slabs. For normally incident light, the transmission spectrum of a photonic crystal slab exhibits Fano lineshape features [5], which result from interference between discrete, guided resonance modes and a continuous Fabry-Pérot background [6, 7]. Past work has demonstrated that the lineshape of a guided resonance mode can be tuned by changing device dimensions [8] or incidence angle [9]. In this paper, we demonstrate that the lineshape of a guided resonance mode can be very simply tuned by rotating an external polarizer. The lineshape tuning method we present will enable better control of photonic resonance characteristics in such applications as optical filtering [10], sensing [11, 12] and signal delay.

Below, we first describe theoretically how rotating the angle of an external polarizer tunes the transmission lineshape of a photonic crystal slab from a symmetric (Lorentzian) shape to a highly asymmetric (Fano) shape. We then fabricate and characterize a Suzuki-phase lattice photonic crystal to validate the model and demonstrate lineshape tuning experimentally. We infer that the phase shift, and thus the group delay through the device, also depend on polarization. In particular, the group delay can be tuned from positive to negative by rotating the output polarizer. These effects all originate from changing the ratio between the guided-resonance transmission (indirect path) and the background, Fabry-Pérot transmission (direct path) of the photonic crystal.

2. Experimental system

Figure 1 shows a schematic of the system. Light propagates in the +z-direction. The first polarizer (P1) sets the polarization angle θi of light incident on the photonic crystal (PhC). The photonic crystal (PhC) is assumed to support a guided resonance mode (GRM) that is singly degenerate and couples to normally incident radiation with well-defined polarization. That is, there is a polarization direction for a linearly-polarized, normally-incident plane wave for which the overlap integral with the mode is finite, whereas the overlap integral for a plane wave polarized in the orthogonal direction is zero. We use a Suzuki-phase lattice photonic crystal for demonstration [13, 14]; other 1-D or 2-D periodic photonic crystal structures with the required mode characteristics may also be used. Due to coupling of incident light to the GRM, light transmitted through the PhC will in general not be polarized at angle θi; a similar effect has previously been used to demonstrate polarization conversion in a reflective geometry [15]. We place a second polarizer (P2) after the photonic crystal with polarization angle θ′p = θi + θp, where θp is the relative angle between the two polarizers. We assume that the component of transmitted light that is polarized orthogonal to θ′p is lost from the system (e.g. absorbed by a polarizer, or redirected elsewhere by a polarizing beam splitter).

 figure: Fig. 1

Fig. 1 Schematic of system for tuning transmission lineshape.

Download Full Size | PDF

A fully-suspended silicon PhC membrane device is fabricated in the silicon-on-insulator (SOI) platform using standard e-beam lithography and ICP-RIE dry etching. We then use wet under-etching to release the membrane and transfer it to a host carrier wafer [16]. Backside etching is used to tune the GRM to the wavelength of interest [17]. Fig. 2(a) shows the SEM picture of the fabricated device, which has a lattice constant a = 510 nm, hole diameter d = 346 nm and slab thickness t = 200 nm.

 figure: Fig. 2

Fig. 2 Fabricated structure (a) and simulated transmission spectra of x- and y-polarized incident light (b). The scale bar indicates 1μm.

Download Full Size | PDF

The simulated transmission spectrum is shown in Fig. 2(b). The spectrum is different for x- and y- polarized light; due to the rectangular unit cell of the Suzuki-phase lattice, the two directions are not equivalent. For the x-polarization, there is a resonance mode at 1570 nm with a quality factor of approximately 70. Simulations indicate that this mode corresponds to the fifth TE photonic band. The relatively low quality factor corresponds to a relatively wide peak, which simplifies the characterization of lineshape in experiments. For y-polarization, there is no resonance mode visible in the wavelength range shown.

3. Analytical model for transmission lineshape

We can write an analyical expression for the transmission lineshape by adding the transmission coefficients of the direct path (Fabry-Pérot resonance) (td) and guided resonance (tg).

The transmission coefficients of the photonic crystal slab are [6]:

td1=td0(cosθix^+sinθiy^)
tg1=tg0cosθix^
where td0 is the transmission coefficient through a homogeneous dielectric slab, and
tg0=f11+jε
is the transmission coefficient of the guided resonance. ε = (ωω0)/γ is the normalized detuning, ω0 is the resonance frequency, γ is the damping coefficient of the resonance, and j is the imaginary unit. The prefactor f is equal to −(td0 ± rd0) to ensure energy conservation. The plus and minus signs correspond to resonant modes that are even and odd with respect to the xy-mirror plane. The mode we consider in this paper is even, so we keep only the plus sign in the following derivation. Moreover, since the bandwidth of the guided resonance is much narrower than the background Fabry-Pérot features in the spectrum, we can neglect the dispersion of td0 and assume it to be a complex constant.

After the second polarizer,

td2=td0[cosθicos(θi+θp)+sinθisin(θi+θp)]p^=td0cosθpp^
tg2=tg0cosθicos(θi+θp)p^
where is the unit vector aligned with the second polarizer. The total transmission coefficient is thus:
t=td2+tg2=jp^cosθptd0ε+q1+jε

The transmission spectrum T is described by the Fano resonance formula:

T=|t|2=|td2+tg2|2=cos2θp|td0|2|ε+q|21+ε2
where
q=jrd0td0cosθicos(θi+θp)cosθpjcosθpcosθicos(θi+θp)cosθp

Because rd0/td0 is always purely imaginary [6], the first term of Eq. (8) is purely real, the second term is purely imaginary, and therefore, the asymmetry factor q is a complex number. This is in contrast to the system without the second polarizer, for which the asymmetry factor is real.

4. Comparison of theory and experiment

We measure the transmission spectrum while rotating the second polarizer at 4° intervals over 180°. We plot the measured data as a 2D function of frequency and the relative angle between the polarizers in Fig. 3(a).

 figure: Fig. 3

Fig. 3 Measured (a) and fitted (b) transmission spectra as a function of the relative angle between polarizers.

Download Full Size | PDF

We use Eq. (7) to fit the measured data. There are four unknown fitting parameters corresponding to the model: the transmission coefficient |td0|, the guided resonance frequency ω0, the resonance damping coefficient γ, and the incident polarization θi (i.e. the angle between the first polarizer and the mode polarization). The reflection coefficient can be calculated by |rd0|=1|td0|2. An additional fitting parameter F is used to take into account measurement normalization. Figure 3(b) shows the fitted function of Fig. 3(a), with fitting parameters F = 1.90, |td0| = 0.633, ω0a/2πc = 0.319, γa/2πc = 0.0238 and θi = 122°. Visually, the fitted spectrum reproduces the features of the experimental spectrum quite well, while smoothing out experimental noise.

To clearly illustrate lineshape tuning, we select several different angles from Fig. 3 and plot the corresponding spectra. Figure 4(a) and Fig. 4(b) show measured and fitted spectra, respectively. As the relative angle between the polarizers increases, the transmission spectrum changes from an asymmetric Fano lineshape at 10° to a symmetric Lorentzian lineshape at 90°. Further rotation of the second polarizer again yields a Fano lineshape. At θp = 150°, the transmission spectrum is a straight line. This occurs because the guided resonance path is perfectly cancelled by the second polarizer. The fitting result shows that θi is approximately 120°, and θ′p = θi + θp ≈ 270°, indicating that the second polarizer is nearly aligned with the y-axis, perpendicular to the mode polarization. The non-zero transmission arises from the Fabry-Pérot direct path.

 figure: Fig. 4

Fig. 4 Measured (a) and fitted (b) spectra for different relative angles between the two polarizers.

Download Full Size | PDF

5. Discussion

To gain better physical insight into the tunable lineshape, we may recast the transmission in the form of a transfer function for a first order linear system [18]. In particular, the transmission coefficient of Eq. (6) can be written as

t=p^Tεεzεεp
where T is a complex factor which represents the transmission coefficient far from resonance, and εp=εpR+jεpI and εz=εzR+jεzI are the complex pole and zero of the system. By comparing Eq. (6) and Eq. (9),
T=cosθptd0
εz=q
εp=j

As define above, ε is real. However, we are free to consider Eq. (9) as a complex function of a complex variable ε = εR + I. The response along the real axis then gives the physical behavior of the system.

The θp dependence of the zero point (εz) of the system shown in Eq. (11) plays an important role in tuning the line shape. For simplicity, we choose θi to be 45° and |td0|=0.5. Figure 5 shows the pole and zeros of the transmission coefficient for different θp. For all cases of θp, the pole point is fixed at (0, 1) in the complex plane while the zero point varies. In Fig. 6, we choose three special cases when θp is equal to 0°, 60° and 90°, respectively, to illustrate how zero position affects the system response. Different θp angles are shown in different rows in Fig. 6. The first and second columns show the system amplitude and phase response in the complex plane.

 figure: Fig. 5

Fig. 5 Pole (red cross) and zeros (blue circles) of transmission coefficients for different θp.

Download Full Size | PDF

Figure 6(a) shows the amplitude response for θp = 0°, when the polarizers are aligned. The zero is at (−0.5, 0.5) in the complex plane and is visible as an amplitude minimum (black region). The pole is at (0, 1) and is visible as an amplitude maximum (white region). The amplitude response along the real axis is shown as a superimposed red line (see red scale bar on right). The transmission first decreases and then increases as a function of frequency. The phase response is shown in Fig. 6(b). The zero and pole are connected by a line of phase discontinuity. The phase response on the real axis is again shown as a superimposed red line. Both the amplitude and phase response change with polarizer angle.

 figure: Fig. 6

Fig. 6 Complex amplitude (first column) and phase(second column) and modeled transmission amplitude, phase shift and group delay through second polarizer when θp = 0° ((a),(b),(c)), 60° ((d),(e),(f)), 90° ((g),(h),(i)).

Download Full Size | PDF

Figure 6(d) and Fig. 6(e) show the case of θp = 60°. The pole remains fixed at (0,1), while the zero has moved to (0.36,1.36). The real transmission and phase response reflect this change, showing inverted behavior relative to Fig. 6(a) and Fig. 6(b).

Figure 6(g) and Fig. 6(h) show the case when θp is equal to 90°, and the polarizers are perpendicular to one another. In this case, the zero is not visible on the plot. The amplitude lineshape is a symmetric Lorentzian, and the phase response decreases by π over the resonance.

We thus see that by tuning the angle of the output polarizer, we not only can change the lineshape of the transmission spectrum, but also the phase change through the structure. The tunable phase change can be used to obtain a tunable group delay, calculated as Δτ = −dϕ/dω. Figure 6(c), Fig. 6(f) and Fig. 6(i) show the dimensionless group delay, (Δτ)c/a. For a fixed frequency (fixed εR ≈ −0.5), the group delay can be tuned from negative to positive by changing θp.

The physical value of the delay can be obtained from the graph by multiplying by a/c = 1.7 fs. For the relatively low quality factor mode studied here, the delays are in the 1 to 2 fs range. The delay increases linearly with quality factor, Q. Recently, a photonic crystal lattice has been designed that supports a coupled, linearly-polarized guided resonance mode with theoretical quality factor as high as 105[19]. Such a Q value would increase the delay to the picosecond range.

In this paper, we have demonstrated lineshape and time delay tuning by fixing the incident polarization and tuning the collection polarization. We could also have fixed the collection polarization and tuned the incident polarization. This can be inferred from Eq. (8). Substituting θi with θ′pθp, fixing θ′p to be 45°, and tuning θp from 0° to 180°, the zero point of the system has the same trace as shown in Fig. 5, yielding a similar tuning mechanism for the lineshape.

6. Conclusion

In summary, we have proposed a method for tuning the transmission lineshape through a photonic crystal slab membrane using polarization control. The principle is to tune the transmission ratio between the direct, Fabry-Pérot path and the indirect, guided resonance path. We first presented an analytical model for the transmission spectrum. We then fabricated and characterized a Suzuki-phase lattice photonic crystal slab membrane to validate the model. There is excellent agreement between theoretical predictions and experimental results. We show that the lineshape can be tuned from a highly asymmetric shape to a highly symmetric one, while the corresponding group delay is tuned from negative to positive. The tuning is achieved simply by rotating an external output polarizer. More generally, the polarization angle dependence measurement and surface fitting technique shown in this work provide an accurate method to characterize guided resonance modes. We expect our results to be useful in a range of applications of photonic crystal guided-resonance modes, including filtering and sensing.

Acknowledgments

The authors thank Chenxi Lin and Jing Ma for their help with the fabrication process and fruitful conversation with Yunchu Li. Ningfeng Huang was supported by the Center for Energy Nanoscience, an Energy Frontiers Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Award No. DE-SC0001013. Luis Javier Martínez was supported by the Army Research Office under award no. 56801-PCS. Computing resources were provided by the University of Southern California Center for High Performance Computing and Communication (http://www.usc.edu/hpcc).

References and links

1. A. E. Miroshnichenko, S. Flach, and Y. S. Kivshar, Fano resonances in nanoscale structures, Rev. Mod. Phys. 82, 2257–2298 (2010). [CrossRef]  

2. M. Galli, S. L. Portalupi, M. Belotti, L. C. Andreani, L. O’Faolain, and T. F. Krauss, Light scattering and Fano resonances in high-Q photonic crystal nanocavities, Appl. Phys. Lett. 94, 071101 (2009). [CrossRef]  

3. P. T. Valentim, J. P. Vasco, I. J. Luxmoore, D. Szymanski, H. Vinck-Posada, A. M. Fox, D. M. Whittaker, M. S. Skolnick, and P. S. S. Guimaraes, Asymmetry tuning of Fano resonances in GaAs photonic crystal cavities, Appl. Phys. Lett. 102, 111112 (2013). [CrossRef]  

4. B. Luk’yanchuk, N. I. Zheludev, S. A. Maier, N. J. Halas, P. Nordlander, H. Giessen, and C. T. Chong, The Fano resonance in plasmonic nanostructures and metamaterials, Nat Mater 9, 707–715 (2010). [CrossRef]  

5. U. Fano, Effects of configuration interaction on intensities and phase shifts, Phys. Rev. 124, 1866–1878 (1961). [CrossRef]  

6. S. Fan and J. D. Joannopoulos, Analysis of guided resonances in photonic crystal slabs, Phys. Rev. B 65, 235112 (2002). [CrossRef]  

7. S. Fan, W. Suh, and J. D. Joannopoulos, Temporal coupled-mode theory for the Fano resonance in optical resonators, J. Opt. Soc. Am. A 20, 569–572 (2003). [CrossRef]  

8. J. Song, R. P. Zaccaria, M. B. Yu, and X. W. Sun, Tunable Fano resonance in photonic crystal slabs, Opt. Express 14, 8812–8826 (2006). [CrossRef]   [PubMed]  

9. L. Babić and M. J. A. de Dood, Interpretation of Fano lineshape reversal in the reflectivity spectra of photonic crystal slabs, Opt. Express 18, 26569–26582 (2010). [CrossRef]  

10. W. Suh and S. Fan, Mechanically switchable photonic crystal filter with either all-pass transmission or flat-top reflection characteristics, Opt. Lett. 28, 1763–1765 (2003). [CrossRef]   [PubMed]  

11. Y. Nazirizadeh, U. Bog, S. Sekula, T. Mappes, U. Lemmer, and M. Gerken, Low-cost label-free biosensors using photonic crystals embedded between crossed polarizers, Opt. Express 18, 19120–19128 (2010). [CrossRef]   [PubMed]  

12. M. E. Beheiry, V. Liu, S. Fan, and O. Levi, Sensitivity enhancement in photonic crystal slab biosensors, Opt. Express 18, 22702–22714 (2010). [CrossRef]   [PubMed]  

13. A. R. Alija, L. J. Martínez, P. A. Postigo, J. Sánchez-Dehesa, M. Galli, A. Politi, M. Patrini, L. C. Andreani, C. Seassal, and P. Viktorovitch, Theoretical and experimental study of the Suzuki-phase photonic crystal lattice by angle-resolved photoluminescence spectroscopy, Opt. Express 15, 704–713 (2007). [CrossRef]   [PubMed]  

14. L. J. Martínez, A. R. Alija, P. A. Postigo, J. F. Galisteo-López, M. Galli, L. C. Andreani, C. Seassal, and P. Viktorovitch, Effect of implementation of a Bragg reflector in the photonic band structure of the Suzuki-phase photonic crystal lattice, Opt. Express 16, 8509–8518 (2008). [CrossRef]   [PubMed]  

15. A. Bristow, V. Astratov, R. Shimada, I. Culshaw, M. S. Skolnick, D. Whittaker, A. Tahraoui, and T. Krauss, Polarization conversion in the reflectivity properties of photonic crystal waveguides, IEEE J. Quantum Electron. 38, 880–884 (2002). [CrossRef]  

16. C. Lin, L. J. Martínez, and M. L. Povinelli, Fabrication of transferrable, fully-suspended silicon photonic crystal membranes exhibiting vivid structural color and high-Q guided resonance, J. Vac. Tech. B , in press (2013).

17. A. R. Alija, L. J. Martínez, A. García-Martín, M. L. Dotor, D. Golmayo, and P. A. Postigo, Tuning of spontaneous emission of two-dimensional photonic crystal microcavities by accurate control of slab thickness, Appl. Phys. Lett. 86, 141101 (2005). [CrossRef]  

18. I. Avrutsky, R. Gibson, J. Sears, G. Khitrova, H. M. Gibbs, and J. Hendrickson, Linear systems approach to describing and classifying Fano resonances, Phys. Rev. B 87, 125118 (2013). [CrossRef]  

19. J. Ma, L. J. Martínez, and M. L. Povinelli, Optical trapping via guided resonance modes in a Slot-Suzuki-phase photonic crystal lattice, Opt. Express 20, 6816–6824 (2012). [CrossRef]   [PubMed]  

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 (6)

Fig. 1
Fig. 1 Schematic of system for tuning transmission lineshape.
Fig. 2
Fig. 2 Fabricated structure (a) and simulated transmission spectra of x- and y-polarized incident light (b). The scale bar indicates 1μm.
Fig. 3
Fig. 3 Measured (a) and fitted (b) transmission spectra as a function of the relative angle between polarizers.
Fig. 4
Fig. 4 Measured (a) and fitted (b) spectra for different relative angles between the two polarizers.
Fig. 5
Fig. 5 Pole (red cross) and zeros (blue circles) of transmission coefficients for different θp.
Fig. 6
Fig. 6 Complex amplitude (first column) and phase(second column) and modeled transmission amplitude, phase shift and group delay through second polarizer when θp = 0° ((a),(b),(c)), 60° ((d),(e),(f)), 90° ((g),(h),(i)).

Equations (12)

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

t d 1 = t d 0 ( cos θ i x ^ + sin θ i y ^ )
t g 1 = t g 0 cos θ i x ^
t g 0 = f 1 1 + j ε
t d 2 = t d 0 [ cos θ i cos ( θ i + θ p ) + sin θ i sin ( θ i + θ p ) ] p ^ = t d 0 cos θ p p ^
t g 2 = t g 0 cos θ i cos ( θ i + θ p ) p ^
t = t d 2 + t g 2 = j p ^ cos θ p t d 0 ε + q 1 + j ε
T = | t | 2 = | t d 2 + t g 2 | 2 = cos 2 θ p | t d 0 | 2 | ε + q | 2 1 + ε 2
q = j r d 0 t d 0 cos θ i cos ( θ i + θ p ) cos θ p j cos θ p cos θ i cos ( θ i + θ p ) cos θ p
t = p ^ T ε ε z ε ε p
T = cos θ p t d 0
ε z = q
ε p = j
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.