Straightening of light in a one dimensional dilute photonic crystal.

Light transport in a dilute photonic crystal is considered. The analytical expression for the transmission coefficient is derived. Straightening of light under certain conditions in a one-dimensional photonic crystal is predicted. Such behavior is caused by the formation of a localized state in transversal motion. The main contribution to the central diffracted wave transmission coefficient is due to states, that either close to the conductance band's bottom or deeply localized in the forbidden gap. Both these states suppress mobility in the transverse direction and force light to be straightened. Straightening of light in the optical region along with small reflection make these systems very promising for use in solar cells.

Introduction

Recent developments in material science have made possible the fabrication of photonic crystals that allow the observation of many peculiar effects[1], including perfect reflector for all polarizations over a wide selectable spectrum[2,3], optical Hall effect[4], unidirectional scattering[5] with broken time-reversal symmetry and the propagation of optical beams without spatial spreading. The latter, known as supercollimation, has attracted a lot of attention in various fields of physics. Supercollimation effect has been observed in mesoscopic (~100 μm)[6-11] and macroscopic centimeter-scale photonic crystals[12,13]. The latter results indicate that supercollimation effect is very robust and insensitive to possible irregularities or short-scale disorder in the photonic crystal structure. A standard mechanism that can account for the supercollimation effect of light in two-dimensional photonic crystals is the flat character of dispersion curve in the transverse direction due to the interference of the different transverse components of the wavevectors (see, for example, ref.[14]).

However, as we will show below, in diluted photonic crystals (DPC) exists another physical mechanism, based on photon’s transversal motion restriction, that can lead similar to supercollimation effect[14]. The basic idea of this model can be easily understood by considering the following geometry shown in Fig. 1. First, suppose that a photon falls on the plates normally on 0x direction. In case if the photon’s wavenumber k lies in the photonic band gap region, the transmission coefficient must be suppressed (see, for example, ref.[15]). Next, imagine that the same photon falls down onto the system obliquely, as shown in Fig. 1 and assume that the transversal to plates component k still remains in the photonic band gap. Clearly, for small and moderate scattered angles transversal motion is suppressed due to the appearance of localized or low energy states. Hence, as a direct consequence of the restriction in the transverse direction, photon’s propagation parallel to plates will be enhanced.

Figure 1

Geometry of the problem.

Note that in most papers on 1d photonic crystal photon main motion is normal to the plates, see, for example[16], and for recent reviews[17]. In our manuscript we are considering a different geometry (see Fig. 1) where photon mainly propagates parallel to the plates. The novelty of our consideration is in different geometry and in corresponding theoretical approach. Our approach allows to obtain closed analytical expression for transmission coefficient.

In this paper we aim to present a complete and quantitative theoretical description of the supercollimation[14] and straightening of light[13] effects within simplified DPC model, taking into account the above-mentioned restriction in the transverse direction. We will see, that the simplified DPC models, developed in refs[18,19], can help to provide new insights into the properties of the mentioned effects.

The main simplification is related to the ignoring of backscattering in a DPC. Within this approach, we have investigated the transport of light through a one-dimensional metallic photonic crystal with transverse to incident direction inhomogeneity. Independence of transmission coefficient on the incident light wavelength was found[18]. Beside that, we have predicted and experimentally observed[19] a capsize, a drastic change of polarization to the perpendicular direction in DPC. The present work takes one step further in the study of the supercollimation effect. We will remove the limitations on the normal incident light, discussed in refs[18,19] and consider general oblique incidence case. The latter, as we will see below, can lead to an interesting phenomenon of straightening of light into the normal direction on exit of the photonic crystal. In this sense, the DPC can be used in solar cells to increase their efficiency. The point is that at oblique incidence a large amount of light energy is lost due to reflection. Preliminary straightening would reduce this loss. For other applications of 1d photonic crystals see also[20,21].

Note, that the main difference of DPC model from other photonic crystal systems where supercollimation effects where observed[6-11] is that in the former reflection is negligible due to the fact that the fraction of one component is very small.

Initial Relations

Consider a system with inhomogeneous dielectric permittivity (see Fig. 1) and suppose that a plane wave enters the system from the space at an arbitrary incidence angle θ. Following procedures described in previous work ref.[18] and[22-24], scalar Helmholtz equation is reduced to the following time-dependent Schrödinger equation for a particle with mass k0where

Note, that the parabolic approximation, that is justified, if the characteristic scale L of along z is much longer than a wavelength, that is (see below). The solution of Eq. (1) can be represented through the eigenfunctions of Hamiltonian Eq. (2). Therefore for the solution of Helmholtz equation one haswhere

It follows from Eq. (3) that the local transmission amplitude of a central diffracted wave can be defined as in ref.[18]where L is the system size in the z direction.

The central diffracted wave transmission coefficient that is measured in the experiment can be estimated by the following expressionwhere S is the area of the system. Substituting Eq. (5) into Eq. (6), one has

In order to find the coefficients c let us consider the Eq. (3) for

To proceed, we assume that the wave intruding the system has an amplitude 1 (the region ). From the continuity at , one has . Here we ignore the reflected waves from the dilute system. Within this approach, multiplying both sides of Eq. (8) by and integrating over the surface, one has

Substituting Eq. (9) into Eq. (7), we arrive at the final result for the transmission coefficientwhere is a two dimensional vector on the xy plane. The equation above is the generalization of the previous result, obtained in ref.[18], in case of oblique incidence of light at an arbitrary angle θ. In the succeeding subsections we examine the limitations of the found equation for different models.

Transmission Coefficient

Note that when (E is the bottom value of first energy band), the transmission coefficent (Eq. (10)) is equal to zero. When k0 is inside the allowed band, transmission coefficient can be represented by quasi-momentum in the following waywhere is integrated over the first Brilloin zone. For sake of simplicity and demonstration of the results, we will carry out further consideration in one-dimensional case.

Kronig-Penney model

Assuming that plates are positioned periodically along the x axis (see Fig. 1), we can represent the cross section of the potential as a multitude of square potential wells. Metal layers of PC imitate the potential wells with depth and width b. The vacuum layer is charachterized by potential energy and . So the problem described by Eq. (2) is brought to the Kronig-Penney model[25]. The transmission coefficient T for one-dimensional configuration can be written aswhere L is the system size in the x direction. Here sum is running over the bands with energy . When k0 lies in the energy gap the integration over covers whole Brilloin zone −π/a, π/a. Using Bloch theorem , where is a periodical function, one obtains

Changing the variables one findswhere wave functions in different layers of PC can be found from Eq. (3) (see Supplement). By substituting into Eq. (14), for transmission coefficient we will havewhere . The sum in Eq. (15) include states from different zones with both negative and positive energies with the same quasimomentum . The number of terms in the sum depends on parameters a, k0, b and .

Straightening of Light

Using formula from Supplement we numerically calculate the transmission coefficient of central diffracted wave for the following parameters , , , . The results are shown in Fig. 2. We compare with the vacuum case . In this case the sum Eq. (15) contains only one term with . Clear straightening effect is seen from Fig. 2.

Figure 2

Transmission coefficient dependence on incident angle. Blue line is the theoretical plot with PC. Red line is the background theoretical plot without PC. In the inset the energy band scheme of transverse motion is shown.

Indeed, independent of incident angle, light leaves the system mainly on the 0z direction because T is close to unity. It is clear, as mentioned above, that different energies give contribution to T. However, as indicate our numerical calculations, the main contribution to T comes from relatively small positive and relatively large negative states.

The state, that close to the bottom of the conduction band, characterized by the small group velocity in transverse direction caused by flat dispersion curve and can lead to supercollimation effect, discussed in refs[14]. As for the states with negative energies, they form a very narrow band gap with almost localized states and with limited transport properties in transverse direction, see inset in Fig. 2.

From theoretical point of view, it is clear that under certain conditions when the incident photon wavenumber is trapped in the gap of energy spectrum of transversal motion, in DPCs light can propagate without spreading out. These states also lead to the suppression of mobility in transverse direction and force photons move only parallel to plates. Note, that depending on the incident angle, the first or second type of states contribute differently to straightening effect, see Appendix. However, the resultant T, taking into account the contribution of all states (negative and positive) becomes almost incident angle independent and close to unity value.

Conclusion and Discussion

We have considered the transport of light through a 1d dilute photonic crystal model at oblique incidence. Our theoretical study, based on Maxwell’s equations with a spatially dependent inhomogeneous dielectric permittivity. For certain parameters’ values of DPC the emerging light propagates in the normal direction despite the oblique incidence. This straightening effect is intimately connected with limited transport properties of photons in transverse direction.

When considering the straightening effect we assume continuity of scalar wave field at the interface plane between two media (xy in Fig. 1). Coming back to the em wave case this means that our consideration is correct for s-polarized waves, electric field vector of which is directed on 0y in geometry of Fig. 1. As mentioned above the straightening effect in visible range could be utilized in solar cell elements to make their absorption efficiency higher. It could seem that the polarization dependence of the effect will decrease the application efficiency because the natural light consists of both s- and p-polarizations. However it follows from Fresnel formulae that the p-polarized wave is essentially reflected only for very large incident angles . In contrary a s-polarized wave is reflected even for moderate incident angles. Therefore using DPC to straighten s-polarized light has sense.

As was mentioned above, in the discussed model of dilute photonic crystal, we have ignored the back-scattered waves and reflection. In order to justify formally this approach, we estimate below the value of reflection amplitude r, using the standard definitionwhere the effective dielectric constant defined by . Taking and dielectric fraction , one has and , i.e. neglecting back-scattered waves and reflection is a reasonable approximation while calculating T. As it follows from the above estimates a one-dimensional DPC on the one hand can straighten the oblique incident light in the visible range. On the other hand DPC provides a negligible reflection and diffraction less light propagation, forecasting the potential of DPCs in increasing the absorption efficiency in solar cells. Note that a 1d photonic structure like Fig. 1 could be made by litography on the surface of solar absorbing element at several wavelengths depth.

Supplement. Straightening of light in a one dimensional dilute photonic crystal