## Abstract

When a wave impinges obliquely to the interface of a Photonic Crystal (PhC), the wave can be completely reflected in counter-propagating direction instead of the usually expected specular direction. However the beam is totally specularly reflected with a simple modification of the surface termination. The analysis of the time average Poynting vector evidences that PhC termination modifies the energy flow and determines the reflection properties.

© 2007 Optical Society of America

## 1. Introduction

Propagation properties in Photonic Crystals (PhCs) [1–2], in particular those related to dispersion such as negative refraction [3] self-collimating [4] or superprism effect [5], are usually analyzed using dispersion surfaces, also known as Equi-Frequency Surfaces (EFS).

The electromagnetic field in a PhC can be decomposed using a Bloch waves expansion. EFS are the locus of wavevectors, in the reciprocal space and for a given frequency, associated to the propagating modes of the Bloch expansion. EFS represent the wavevectors of the self-consistent optical field in the unbounded structure and are obtained considering the PhC as an infinitely periodic structure.

The wavevectors excited inside the PhC are determined applying the continuity of the tangential component of the wavevectors across the boundary between the PhC and the external medium [3]. No propagating modes are excited inside the PhC when there are no points on the EFS that satisfy the continuity condition. In such case, the incident wave experiences a photonic gap (a partial gap when it happens in a limited range of incident angles) and the wave is totally externally reflected by Bragg diffraction.

From a general point of view the application of wavevectors tangential continuity
ensures the linear momentum conservation and, following the Jackson [6] classification, it describes a *kinematic*
property i.e. it determines the light propagation direction.

On the other hand the application of continuity of the tangential component of the
electric and magnetic fields determines the *dynamic* properties of
the electromagnetic field, i.e. it determines the intensities of the refracted and
reflected-diffracted beams.

Actually it is not surprising that a dynamic relationship, such as Fresnel relations that are obtained for homogeneous media, cannot generally be applied to PhCs [7], even when their kinematic properties, such as refraction-like behaviour [3], can be applied to PhCs. For instance uncoupled modes that do not match the symmetry between the incident wave and the modes in PhC, can give zero transmission even without a band-gap [8]–[9]. This is a typical kinematic property that cannot be analysed simply using EFS even if other phenomena, such as changes of refraction angles obtained by rotating directors of liquid crystals, usefully can be studied with help of EFS [10]. Essentially, even if inside the first Brillouin zone the EFS has the shape of a homogenous medium (circle in 2D, sphere in 3D), the PhC does not generally behave as an effective refractive medium. In addition to the previously mentioned symmetry mismatch, also the propagation inside the crystal depends strongly on the interface termination that plays a fundamental role for surface localised states at PhC boundaries [11]–[14], that is useful to improve focusing with negative refraction PhC [14–16] or the beaming effects [18]–[19].

In the following sections we show that an incident wave inclined respect to the surface of a PhC can completely be back-reflected in a counter-propagating direction. This can be interpreted as an enhancement of a diffraction order, where the specular reflected wave is completely suppressed.

The surface termination plays a key-role in the distribution of energy of the incident wave between the specularly reflected wave and the back-diffracted wave in a counter-propagating direction. Simply changing the surface termination the wave can be completely back reflected or completely specularly reflected.

The observed effect is somehow related to the strong variation of diffraction efficiency experimentally observed in Ref. [20] using the Littrow geometry, which substantially is the same configuration of the present study. However in Ref. [20] the external wave can couple with internal Bloch waves because in such case authors works far from a bandgap. Then in Ref. [20] this behaviour is interpreted as an interference effect between the light diffracted by the surface corrugation and the light scattered from the Bloch modes of the bulk photonic crystal. The surface termination in Ref. [20] acts like a very special antireflection coating.

In our case we are involved with a partial bandgap and the interpretation of the phenomenon cannot be performed in analogy with [20]. As explained earlier, EFS can simply give information on the propagation but cannot distinguish between two different surface terminations. However also other approaches, even taking into account the whole Bloch expansion including real-propagating and complex-evanescent wavevectors [8] cannot simply explain the surface termination dependence of the reflected power.

We will show that an analysis of the energy flow allows an interpretation of the influence of the surface termination in terms of matching between symmetry of the energy flow inside and outside the PhC.

Finally we underline that the energy inside the PhC close to the interface, is distributed over vortexes reproducing the hexagonal symmetry of the PhC, which is itself a special relevant result that will be analysed in a forthcoming paper.

## 2. EFS design procedure: the bulk waves

We have analysed a 2D PhC with a background dielectric
*ε _{Si}* =11.95 (silicon at telecom infrared
wavelength) and air-holes arranged on a hexagonal lattice. We have limited our
attention to the TM polarization (electric field out of the plane, directed along
the holes axis).

The PhC is cut along the direction ΓK and the normal of the PhC with
external homogenous silicon is along the ΓM direction (Fig. 1). In the next sections we have studied the influence
of the cut position, in particular the cut in two relevant symmetric positions: in
the middle and on the limit of the elementary cell, as defined in Fig. 1(a). The incident plane wave, with a Gaussian
transverse profile, incomes from the silicon background with a direction of
30° respect to the normal of the PhC. The considered incident wavevector
is highly symmetric respect to the hexagonal lattice symmetries, Fig. 1(b). A partial gap appears when there are no
intersections between EFS and the line intersecting the end of the
*k _{in}* vector normal to the PhC along the ΓM
direction. In Fig. 2(a) we studied the dependence of this partial gap by
the normalised hole radius

*r*/

*a*.

The incident wave satisfies the Bragg condition respect to the
*h*⃗_{11} = (1,1) lattice vector as defined in Fig.
1(b):|*k _{in}*| ≡

*n*2

_{Si}*π*/

*λ*= 2

*π*/

*a*so that normalized frequency is

*ω*≡

_{n}*a*/

*λ*= 1/

*n*= 0.289. In fact the wavevector

_{Si}*k*⃗

_{nr}satisfies the Bragg relation with the

*h*⃗

_{11}= (1,1) vector of the reciprocal lattice :

*k*⃗

_{i}=

*k*⃗

_{nr}+

*h*⃗

_{11}, Fig. 1(b).

In correspondence of the normalised frequency
*ω _{n}*=0.289, for TM polarization, a partial gap
opens in the range 0.29<r/a<0.42. In particular for
r/a<0.33 small regions of EFS appear in correspondence of the K points in
the reciprocal space, as a consequence of the partial overlap between the first and
the second band, the green circles in Fig. 2(b). As a general property the exponential decay of the
evanescent wave inside a gap is at maximum in the middle of the gap and increases
when the gap increases in the reciprocal space [21]. Because of the bands overlap, shown in Fig. 2(b), the partial gap is quite small in the reciprocal
space in the range 0.29<r/a<0.33. In this range the
exponential decay is small too and the evanescent wave penetrates significantly in
the PhC for a length equivalent to many lattice constants. This enhances the
interaction of the incident wave with the PhC.

## 3. Negative reflection

To maximize the penetration of the wave in the PhC as previously explained, we have
fixed r/a=0.31. The incident wave incomes from the homogeneous silicon; the incident
wavevector *k _{i}* is defined in Fig. 1(b) and Fig. 2(b). Two different PhC surface terminations will be
considered. First a termination as in Fig. 1(a) was analyzed, i.e. with a truncation of the PhC at
the end of an elementary cell without any modification of the holes along the
surface line. We found out that - as expected - the wave penetrates in the PhC for
many periods and, after an interaction-time, it is completely reflected in the
silicon, according to the partial gap for the selected radius.

In Fig. 3(a) a movie versus time is shown, obtained using a
Finite Difference Time Domain (FDTD) code, of a wave packet centred on the
wavelength *λ* = *a*/0.289 and with a
packet length of 12 *λ*. We have observed that more than
80 % of the incident energy is transferred to the diffracted wave that propagates in
the same direction as the incident wave but they have opposite orientations
*k*⃗_{nr}: in the following, we refer to this wave as a negatively reflected beam. We
have also observed that a very small part of the energy is transferred to the
specularly reflected beam, which actually coincides with the family of all the
diffracted waves centred in reciprocal lattice nodes (0, n), where n=0,-1, 1,-2, 2.

However if the surface is terminated as in position 2 of Fig. 1(a), i.e. cutting in the middle the holes of the surface line, the wave is totally reflected in a specular direction, Fig. 3(b).

## 4. Analysis of the surface termination influence on the negative reflection

Using analysis of section 2 we found modes that can propagate in PhC and external homogeneous silicon satisfying the tangential wavevector continuity (i.e. that satisfy the photon linear momentum continuity). EFS are a bulk property of media and such analysis couples bulk properties inside and outside the PhC. However, EFS analysis cannot determine the distribution of energy across each mode. This requires the additional application of the tangential continuity of the electric and magnetic fields, as for instance in Ref. [8] or [21].

On the other hand, an energy flow picture gives a lot of interesting information useful for the interpretation of the surface termination dependence of reflected beams. Figure 4 shows the time-average Poynting vector calculated for both the studied surface terminations close to the interface. In Fig. 4 the incident wave has the same characteristics of Fig. 3 except for a continuous wave excitation, instead of the finite length pulse considered in Fig. 3.

Inside the PhC, which is in the upper part of Fig. 4(a) and Fig. 4(b), the energy flow is distributed along the vortex, arranged following the hexagonal symmetry. In the middle of the holes the flow is parallel to the interface with external medium.

In Fig. 4 we can see that symmetry of the energy flow is preserved across the interface. When the PhC surface terminates in the middle of holes, the energy flow is parallel to the interface in the PhC, Fig. 4(b), and this flow is conserved in the external homogenous silicon. This enhances the specular reflected beam: in fact in the external medium the sum of the incident with the specular reflected beam originate a net flow parallel to the surface, whereas the components along the normal to the interface are opposite each other and their interference establishes a standing wave in the external medium along the normal direction.

Conversely when surface termination does not cut the holes, the vortex arrangement of the energy flow in the PhC close to the interface is preserved in the external medium. This arrangement of the time average Poynting vector gives a zero flow across the lines perpendicular to the incident wavevector. Such topology property clearly help an enhancement of the negative reflected beam, where incident and reflected beam travel in an opposite direction, Fig. 4(a). Considering the plane waves decomposition of the wavefield inside and outside the PhC the vortex position depends nonlinearly on the amplitude and on the direction of the plane wave components, and will be deeply analyzed in a forthcoming paper. We underlined here that a vortex close to the surface inside and outside the PhC reveals the complexity of the wavefield in the interface region that cannot simply be resumed to the interference of two plane waves travelling in opposite directions [22].

The energy flow topology suggests the way to further enhance the negatively reflected
beam. Looking carefully at Fig. 4(a) we see that the vortexes inside the PhC and outside
are not exactly aligned. The misalignment is slight and displacing the holes line
along the interface by a quantity equal to 0.2 *a* we obtain a
totally reflected beam, as shown in Fig. 5. The displacement of last holes lines is a blazing
effect, that concentrates all the energy across a diffraction order [23]. Electromagnetic waves crossing the photonic crystal
interact strongly with the boundary, but it’s still a difficult problem
to establish a rigorous link between finite and infinite PhCs [24]–[25].

## 5. Conclusions

We have shown that surface termination strongly influences the reflection from a PhC. In particular a counter-propagating reflected beam can be obtained also for an oblique incident angle. A such negatively reflected beam is very sensitive to the PhC termination. We have proved that a study of the energy flow close to the interface is extremely useful to analyze this surface termination influence. In the case of a negatively reflected beam the energy flow is distributed across vortexes arranged with the same hexagonal symmetry of the PhC. When the surface termination breaks the symmetry of the energy flow, the negatively reflected beam is suppressed. Conversely, if the termination is slightly modified in accordance with the symmetry of the vortex flow, the negatively reflected beam is further enhanced.

## References and links

**1. **J. D. Joannopulos, R. D. Mead, and J. N. Winn, *Photonic crystal: Molding the flow of light*,
(Princeton University Press, Princeton1995).

**2. **K. Sakoda, *Optical Properties of Photonic Crystals*,
(Springer Verlag2001).

**3. **M. Notomi. “Theory of light propagation in
strongly modulated photonic crystals: Refractionlike behavior in the
vicinity of the photonic band gap,” Phys.
Rev. B **62**, 10696–10705
(2000). [CrossRef]

**4. **P. T Rakich, M. S. Dahlem, S. Tandon, M. Ibanescu, M. Soljacic, G. S. Petrich, J. D. Joannopoulos, L. A. Kolodziejski, and E. P. Ippen, “Achieving centimetre-scale
supercollimation in a large-area two-dimensional photonic
crystal,” Nat. Mater. **5**, 93–96
(2006). [CrossRef] [PubMed]

**5. **H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. Tamamura, T. Sato, and S. Kawakami, “Superprism phenomena in photonic
crystals,” Phys. Rev. B **58**, R10 096–099
(1998). [CrossRef]

**6. **J. D. Jackson, *Classical Electrodynamics*, 2nd ed.
(Wiley, 1975)

**7. **B. T. Schwartz “Dynamic properties of photonic
crystals and their effective refractive index,”
J. Opt. Soc. Am. B **22**, 2018–2026
(2005) [CrossRef]

**8. **K. Sakoda, “Transmittance and Bragg reflectivity
of two-dimensional photonic lattices,”
Phys. Rev. B **52**, 8992–9002
(1995) [CrossRef]

**9. **Z. Y. Yuan, J. W. Haus, and K. Sakoda, “Eigenmode symmetry for simple cubic
lattices and the transmission spectra,”
Opt. Express **3**, 19–27
(1998). [CrossRef] [PubMed]

**10. **H. Takeda and K. Yoshino, “Tunable refraction effects in 2D
photonic crystals utilizing liquid crystals,”
Phys. Rev. E **67**, 056607 (2003). [CrossRef]

**11. **R. D. Meade, K. D. Brommer, A. M. Rappe, and J. D. Joannopulos, “Electromagnetic Bloch waves at the
surface of a photonic crystal,” Phys.
Rev. B **44**, 10961–10964
(1992). [CrossRef]

**12. **W. M. Robertson, G. Arjavalingam, R. D. Meade, K. D. Brommer, A. M. Rappe, and J. D. Joannopoulos, “Observation of Surface Photons on
Periodic Dielectric Arrays,” Opt. Lett. **18**, 528–530
(1993). [CrossRef] [PubMed]

**13. **F. Ramos-Mendieta and P. Halevi, “Surface electromagnetic waves in
two-dimensional photonic crystals: Effect of the position of the surface
plane,” Phys. Rev. B. **59**, 15112–15120
(1999). [CrossRef]

**14. **Y. A. Vlasov, N. Moll, and S. J. McNab, “Observation of surface states in a
truncated photonic crystal slab,” Opt.
Lett. **29**, 2175–2177
(2004). [CrossRef] [PubMed]

**15. **S. S. Xiao, M. Qiu, Z. C. Ruan, and S. L. He, “Influence of the surface termination
to the point imaging by a photonic crystal slab with negative
refraction,” Appl. Phys. Lett. **85**, 4269–4271
(2004). [CrossRef]

**16. **X. D. Zhang, “Effect of interface and disorder on
the far-field image in a two-dimensional photonic-crystal-based flat
lens,” Phys. Rev. B. **71**, 165116 (2005). [CrossRef]

**17. **A. Martinez and J. Marti, “Negative refraction in
two-dimensional photonic crystals: Role of lattice orientation and interface
termination,” Phys. Rev. B. **71**, 235115 (2005). [CrossRef]

**18. **E. Moreno, F. J. Garcia-Vidal, and L. Martin-Moreno, “Enhanced transmission and beaming of
light via photonic crystal surface modes,”
Phys. Rev. B. **69**, 121402 (2004). [CrossRef]

**19. **P. Kramper, M. Agio, C. M. Soukoulis, A. Birner, F. Muller, R. B. Wehrspohn, U. Gosele, and V. Sandoghdar, “Highly directional emission from
photonic crystal waveguides of subwavelength
width,” Phys. Rev. Lett. **92**, 113903 (2004). [CrossRef] [PubMed]

**20. **G. von Freymanna, W. Koch, D. C. Meisel, M. Wegener, M. Diem, A. Garcia-Martin, S. Pereira, K. Busch, J. Schilling, R. B. Wehrspohn, and U. Gösele, “Diffraction properties of
two-dimensional photonic crystals,” Appl.
Phys. Lett. **83**, 614–616
(2003). [CrossRef]

**21. **E. Istrate and E. H. Sargent, “Photonic crystal heterostructures
and interfaces,” Rev. Mod. Phys. **78**, 455–481
(2006). [CrossRef]

**22. **K. O’Holleran, M. Padgett, and M. R. Tennis, “Topology of optical vortex lines
formed by three, four, and five plane waves,”
Opt. Express **14**, 3039–3041
(2006). [CrossRef] [PubMed]

**23. **E. Popov, B. Bozhkov, and M. Nevière, “Almost perfect blazing by photonic
crystal rod gratings,” Appl. Opt. **40**2417–2422
(2001). [CrossRef]

**24. **D. Felbacq and R, Smaâli, “Density of states for finite
photonic crystals,” Phys. Rev. B **67**, 085105 (2003). [CrossRef]

**25. **P. A. Belov and C. R. Simovski, “Boundary conditions for interfaces
of electromagnetic crystals and the generalized Ewald-Oseen extinction
principle,” Phys. Rev. B **73**, 045102 (2006). [CrossRef]