Negative refraction can make non-diffracting
beams

doi:10.1364/OE.16.014582

« Show journal navigation

Negative refraction can make non-diffracting beams

M. Kaliteevski, S. Brand, R. A. Abram, A. J. Gallant, and J. M. Chamberlain »View Author Affiliations

Optics Express, Vol. 16, Issue 19, pp. 14582-14587 (2008)
http://dx.doi.org/10.1364/OE.16.014582

View Full Text Article

Acrobat PDF (196 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract

– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (15)
Cited By
Figures (4)
Metrics

Abstract

We report the results of simulations relating to the illumination of a structure consisting of a slab constructed from a 2-D hexagonal array of metal rods with a terahertz frequency source. As a consequence of negative refraction an essentially non-divergent beam pattern is observed. Although the results presented relate to the terahertz regime they should also be applicable at other frequencies.

The diffraction of waves is a ubiquitous feature of electromagnetism but at the same time is a process which can severely restrict the performance of electromagnetic devices [1

Kishan Dholakia, “Against the spread of the light,” Nature 451, 413 (2008) [CrossRef] [PubMed]

]. In essence diffraction causes a beam of radiation with wavelength λ and width w to diverge by an angle δ which satisfies the relation sin δ≈λ/w. In an attempt to eliminate, or at least to reduce the divergence of beams, Bessel beams [2

J Durnin, J J Mieceli, and J H Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987) [CrossRef] [PubMed]

], and more recently Airy packets [3

G A Siviloglou, J Broky, A Dogariu, and D N Christodoulides, “Observation of accelerating Airy beams,” Phys.Rev.Lett. 99, 213901 (2007) [CrossRef]

], have been employed, where the transverse profile of the beam intensity is a non-monotonic function of position. However, it is still necessary for the characteristic width of the beam to be a substantial number of wavelengths to provide something approaching non-divergent behaviour over large distances. Here we demonstrate that an essentially non-divergent beam of terahertz radiation with a width of a few wavelengths can be produced by the use of a flat lens made from a metallic photonic crystal in the form of a 2-D hexagonal array of metal rods operating in the negative refraction regime.

In 1968 Veselago [4

V G Veselago, “The electrodynamics of substances with simultaneously negative values of ε and µ,” Sov. Phys. Usp. 10, 509–514 (1968) [CrossRef]

] pointed out that for materials in which both the dielectric constant ε and the magnetic permeability µ have negative values, E×H is antiparallel to the wave vector k. Such materials have been called left-handed materials (LHMs) in contrast to conventional materials, which can be called right-handed materials (RHMs). LHMs are also called negative index materials or negative refraction materials because the negative angle of refraction that occurs when radiation is incident from an RHM can be formally described by Snell’s law using a negative refractive index. Later, it was recognized that negative refraction can occur near an optical resonance in conventional material with positive magnetic permeability[5

J B Pendry, “Negative refraction makes a perfect lens,” Phys.Rev. Lett. 85, 3966–3969 (2000) [CrossRef] [PubMed]

], in quasi-homogeneous metamaterials [6

R A Shelby, D R Smith, and S Shultz, “Experimental verification of a negative index of refraction,” Science 292, 77–79 (2002) [CrossRef]

] and in photonic crystals [7

V. Patanjali, V. Parimi, Wentao T. Lu, Plarenta Vodo, and Srinivas Sridhar, “Photonic crystals: Imaging by flat lens using negative refraction,” Nature 426, 404 (2003) [CrossRef]

–9

M A Kaliteevski, S Brand, J Garvie-Cook, R A Abram, and J M Chamberlain, “Terahertz filter based on refractive properties of metallic photonic crystal,” Opt. Express 16, 7330–7335 (2008) [CrossRef] [PubMed]

]. Various lensing–related applications of LHMs have been discussed, and in particular it has been shown that a flat slab of LHM having a refractive index of the same absolute value as the surrounding material, but opposite sign, acts as a superlens, providing imaging with subwavelength resolution [5

J B Pendry, “Negative refraction makes a perfect lens,” Phys.Rev. Lett. 85, 3966–3969 (2000) [CrossRef] [PubMed]

R A Shelby, D R Smith, and S Shultz, “Experimental verification of a negative index of refraction,” Science 292, 77–79 (2002) [CrossRef]

,10

D O S Melville, R J Blaikie, and C R Wolf, “Submicron imaging with a planar silver lens,” Appl. Phys. Lett. 84, 4403–4405 (2004) [CrossRef]

]. However, in this paper we consider the case illustrated in Fig. 1 where the slab of LHM has a negative refractive index n with a magnitude different from that of the surrounding RHM. The slab is illuminated by a point source in the 2-D plane (a line source in 3-D, parallel to the axis of the metallic rods in the slab).

Fig. 1. Propagation of rays from a 2-D point source through a slab of negatively refracting material.

In Fig. 1, a ray emerging from a source located at distance A from the right side of the slab crosses the optic axis at a distance D from the other side of the slab, where D depends on the angle α between the ray and the optical axis. To be specific, we take the medium surrounding the slab to be vacuum when Snell’s law gives the relationship

\sin α = ∣ n ∣ \sin β

(1)

and it follows that

D + A = \frac{L \tan β}{\tan α} = \frac{L \sqrt{1 - \sin^{2} α}}{\sqrt{{∣ n ∣}^{2} - \sin^{2} α}} .

(2)

It is interesting to note that if |n|<1, the radiation incident on the slab can experience total internal reflection, and only those rays for which the angle of incidence α is smaller than a critical angle α _tot , given by sin α_tot =|n|, can enter the slab. When α→0 rays cross the optic axis near D ₀=L/|n|- A, but when the incident angle corresponds to α_tot , D ₀→∞. Hence, the image of the source will not be a single point (in the 2-D plane), as would be the case if |n|= 1, but will be stretched from the position D ₀ to infinity. Since D ₀ can be made arbitrarily large by reducing |n|, the result suggests that it should be possible to use a slab with a negative effective refractive index of small magnitude to produce an essentially non-divergent beam of radiation from a point source in the 2-D plane. In this paper we illustrate how a non-divergent beam with a width of a few wavelengths can be produced at terahertz frequencies using a slab in the form of a metallic photonic crystal having the required properties.

Figure 2 shows the bandstructure, f(k), for the two-dimensional photonic crystal considered, which is in the form of gold (plasma frequency=8.9 eV) circular rods of diameter 80 µm arranged in a hexagonal lattice with period 200 µm. For radiation in the E polarization, with electric field parallel to the rods, such a structure exhibits an effective plasma frequency [11

S Brand, R A Abram, and M A Kaliteevski, “Complex photonic bandstructure and effective plasma frequency of a two-dimensional array of metal rods,” Phys. Rev. B75, 035102, (2007)

,12

A J Gallant, M A Kaliteevski, D Wood, M C Petty, R A Abram, S Brand, G P Swift, D A Zeze, and J M Chamberlain, “Passband filters for terahertz radiation based on dual metallic photonic structures,” Appl. Phys. Lett. 91, 161115 (2007) [CrossRef]

] of about 1 THz. The dispersion in the upper band is such that negative refraction occurs when radiation is incident on the crystal from vacuum [9

]. In particular near the top of the upper band the equifrequency surfaces are circular and the group velocity and the Poynting vector averaged over a unit cell are antiparallel to the Bloch wavevector. It follows that when radiation is incident on the photonic crystal, the conditions of continuity of transverse wavevector and flow of energy away from the interface and into the crystal require a negative angle of refraction. Further, straightforward trigonometry shows that the negative refraction can be described by a scalar effective refractive index n=-k(f)/k ₀ where k(f) and k ₀ are the magnitudes of the wavevectors at frequency f in the photonic crystal and the vacuum respectively. Also, k(f) decreases as the frequency is increased, and near the top of the band the effective refractive index can be very small in magnitude, making large values of D ₀ possible [13

Dispersion relations have been calculated using a complex photonic bandstructure method described in [11], and finite difference time domain software OMNISIM© has been employed to simulate the beam propagation. Due to minor convergence problems the frequency scales for the two methods are slightly different. In order to provide matching of the two scales, the bandstructure has been renormalized to provide matching of the frequencies at the J point which can be determined for both calculation techniques.

Fig. 2. Band structure [13

] of the hexagonal photonic crystal with lattice constant of 200 µm formed by metallic rods of diameter 80 µm. Horizontal arrows indicate the frequencies for which the modelling of the field pattern is shown in Fig. 3.

Figures 3(a–c) show the 2-D field patterns formed by a flat slab of photonic crystal of thickness L=1200 µm irradiated with line sources (perpendicular to the plane shown) of frequencies 1.621 THz, 1.667 THz and 1.715 THz placed at a distance A=1.8 mm to the right of the structure. For a frequency of 1.621 THz (wavelength λ=185 µm) the beam has a width of about 3λ and it propagates for about 10 mm (50λ) before there is significant divergence. The divergence of the beam is illustrated in Fig. 4a by the black solid line, which shows the intensity profile at a distance of 7.5 mm to the left of the slab (approximately in the middle of Fig. 3), and by the black dashed line which shows the profile at 16 mm (at the left edge of Fig. 3). With increasing frequency, |n| reduces, as shown in Fig. 2, and D ₀ increases, as is implicit in Figs. 3(b) and 3(c).

Fig. 3. The field patterns formed when a line source irradiates the photonic crystal slab at frequencies of (a) 1.621 THz, (b) 1.667 THz and (c) 1.715 THz. The distance between the source and the edge of the slab is A=1.8 mm.

For a frequency of 1.667 THz the divergence can hardly be seen in Fig. 3b, but the spreading of the Gaussian-like profile is clearly seen in Fig. 4b. For a frequency of 1.715 THz, which is close to the top of the band, the beam has a waist of about 5λ and shows no sign of divergence at a propagation distance corresponding to 100λ.

For all three beam patterns shown in Figs. 3 and 4 the beam profile is Gaussian-like in the region close to the slab (at distance 7.5 mm). The wavelengths, effective refractive indexes, angles of total internal reflection, distances D ₀ and waist widths w are shown in Table 1.

Table 1. Parameters relating to the results shown in Figs. 3 and 4.

f (THz)	λ(µm)	n	α _tot (degrees)	D ₀ (µm)	w (µm)
1.621	185	-0.15	8.6	5250	400
1.667	180	-0.095	5.5	9300	550
1.715	175	-0.04	2.3	24500	810

Figure 3 shows that there are two side-lobes in the region where the beams emerge from the left side of the slab in Fig. 3. These are evanescent in nature and there is little energy associated with them, but they are an inherent feature of the intensity pattern. Conventional diffracting behaviour occurs if a screen is placed to block the side-lobes. We should point out that the collimation of the beam is not due to the same self-collimating/self-guiding effect described by Kosaka et al [14

H Kosaka, T Kawashima, A Tomita, M Notomi, T Tamamura, T Sato, and S Kawakami, Appl. Phys. Lett. 74, 1212–1214 (1999) [CrossRef]

] and Chigrin et al. [15

D N Chigrin, S Enoch, C M Sotomayor Torres, and G Tayeb, Opt. Express 11, 1203–1211 (2003) [CrossRef] [PubMed]

] for the dielectric structures they considered. The collimation there is a result of the existence of a straight-line section of the equifrequency contours, meaning that the group velocity is in the same direction for a range of different wavevectors. In the present case, with essentially circular equifrequency contours, the collimation of the beam and the form of the side lobes are a consequence of the detailed interference effects at the particular frequency.

The red and green lines in Fig. 4 show the beam profiles when the source is separated from the slab by the distances A=4.5 mm and A=5.3 mm respectively. It is apparent that the width of the beam is not increased in proportion to A, as would be expected on purely geometrical grounds, and in fact for all three source distances considered, the parameters of the non-diffracting beam differ only slightly (although when the distance is smaller than 1 mm, results not shown here demonstrate that the divergence of the beam becomes significantly more pronounced). However, the data in Table 1 do show a clear increase in the waist of the beam with decreasing angle of total reflection, which is the opposite to what would be expected on geometrical grounds. These results can be explained, at least intuitively, by a Fourier argument where the width of the beam is determined by the spread of its transverse wavevectors within the slab, which in turn is determined principally by sin α_tot . Then, as α_tot is decreased, the wavevector spread is decreased and, through the usual Fourier relationship, there is an increase in the beam width in real space.

Fig. 4. Intensity profile of the beam for frequencies (a) 1.621 THz, (b) 1.667 THz and (c) 1.715 THz at distances from the left-hand edge of the slab of 7.5 mm (solid line), 16 mm (dashed line), and for distances of the source from the right-hand edge of the slab of A= 1.8 mm (black), A=4.5 mm (red), A=5.3 mm (green).

Since there is omnidirectional emission from the line source, the fraction of total power impinging on the slab with an angle of incidence less than α_tot is only α_tot /180, and of this, somewhat less than 10% is transferred into the collimated beam. However, more efficient illumination schemes could be employed to increase the fraction of emitted power that is transmitted.

In conclusion, we have shown that left-handed materials can be employed to produce an essentially non-diffracting beam. Although this has been demonstrated for a system operating in the THz frequency regime, a region in which we have a particular interest, it should also be applicable at other frequencies.

Acknowledgement

The work was funded by EPSRC, UK, via grant no. EP/C534263.

References and links

1.	Kishan Dholakia, “Against the spread of the light,” Nature 451, 413 (2008) [CrossRef] [PubMed]
2.	J Durnin, J J Mieceli, and J H Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987) [CrossRef] [PubMed]
3.	G A Siviloglou, J Broky, A Dogariu, and D N Christodoulides, “Observation of accelerating Airy beams,” Phys.Rev.Lett. 99, 213901 (2007) [CrossRef]
4.	V G Veselago, “The electrodynamics of substances with simultaneously negative values of ε and µ,” Sov. Phys. Usp. 10, 509–514 (1968) [CrossRef]
5.	J B Pendry, “Negative refraction makes a perfect lens,” Phys.Rev. Lett. 85, 3966–3969 (2000) [CrossRef] [PubMed]
6.	R A Shelby, D R Smith, and S Shultz, “Experimental verification of a negative index of refraction,” Science 292, 77–79 (2002) [CrossRef]
7.	V. Patanjali, V. Parimi, Wentao T. Lu, Plarenta Vodo, and Srinivas Sridhar, “Photonic crystals: Imaging by flat lens using negative refraction,” Nature 426, 404 (2003) [CrossRef]
8.	E Cubukcu, K Aydin, E Ozbay, S Foteinopoulou, and C M Soukoulis, “Electromagnetic waves: Negative refraction by photonic crystals,” Nature 423, 604–605 (2003) [CrossRef] [PubMed]
9.	M A Kaliteevski, S Brand, J Garvie-Cook, R A Abram, and J M Chamberlain, “Terahertz filter based on refractive properties of metallic photonic crystal,” Opt. Express 16, 7330–7335 (2008) [CrossRef] [PubMed]
10.	D O S Melville, R J Blaikie, and C R Wolf, “Submicron imaging with a planar silver lens,” Appl. Phys. Lett. 84, 4403–4405 (2004) [CrossRef]
11.	S Brand, R A Abram, and M A Kaliteevski, “Complex photonic bandstructure and effective plasma frequency of a two-dimensional array of metal rods,” Phys. Rev. B75, 035102, (2007)
12.	A J Gallant, M A Kaliteevski, D Wood, M C Petty, R A Abram, S Brand, G P Swift, D A Zeze, and J M Chamberlain, “Passband filters for terahertz radiation based on dual metallic photonic structures,” Appl. Phys. Lett. 91, 161115 (2007) [CrossRef]
13.	Dispersion relations have been calculated using a complex photonic bandstructure method described in [11], and finite difference time domain software OMNISIM© has been employed to simulate the beam propagation. Due to minor convergence problems the frequency scales for the two methods are slightly different. In order to provide matching of the two scales, the bandstructure has been renormalized to provide matching of the frequencies at the J point which can be determined for both calculation techniques.
14.	H Kosaka, T Kawashima, A Tomita, M Notomi, T Tamamura, T Sato, and S Kawakami, Appl. Phys. Lett. 74, 1212–1214 (1999) [CrossRef]
15.	D N Chigrin, S Enoch, C M Sotomayor Torres, and G Tayeb, Opt. Express 11, 1203–1211 (2003) [CrossRef] [PubMed]

OCIS Codes
(260.2110) Physical optics : Electromagnetic optics
(160.5298) Materials : Photonic crystals

ToC Category:
Physical Optics

History
Original Manuscript: July 11, 2008
Revised Manuscript: August 29, 2008
Manuscript Accepted: August 29, 2008
Published: September 2, 2008

Citation
M. Kaliteevski, S. Brand, R. A. Abram, A. J. Gallant, and J. M. Chamberlain, "Negative refraction can make non-diffracting beams," Opt. Express 16, 14582-14587 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-19-14582

Sort: Author | Year | Journal | Reset

References

K. Dholakia, "Against the spread of the light," Nature 451, 413 (2008) [CrossRef] [PubMed]
J .Durnin, J. J. Mieceli, and J. H. Eberly, "Diffraction-free beams," Phys. Rev. Lett. 58, 1499-1501 (1987). [CrossRef] [PubMed]
G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Observation of accelerating Airy beams," Phys.Rev.Lett. 99, 213901 (2007) [CrossRef]
V. G. Veselago, "The electrodynamics of substances with simultaneously negative values of ? and ?," Sov. Phys. Usp. 10, 509-514 (1968) [CrossRef]
J. B. Pendry, "Negative refraction makes a perfect lens," Phys.Rev. Lett. 85, 3966-3969 (2000) [CrossRef] [PubMed]
R. A. Shelby, D. R. Smith, and S. Shultz, "Experimental verification of a negative index of refraction," Science 292, 77-79 (2002) [CrossRef]
V. Patanjali V. Parimi, W. T. Lu, Plarenta Vodo and Srinivas Sridhar, "Photonic crystals: Imaging by flat lens using negative refraction," Nature 426, 404 (2003) [CrossRef]
E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou and C. M. Soukoulis, "Electromagnetic waves: Negative refraction by photonic crystals," Nature 423, 604-605 (2003) [CrossRef] [PubMed]
M. A. Kaliteevski, S. Brand, J. Garvie-Cook, R. A. Abram, and J. M. Chamberlain, "Terahertz filter based on refractive properties of metallic photonic crystal," Opt. Express 16, 7330-7335 (2008) [CrossRef] [PubMed]
D. O. S. Melville, R. J .Blaikie, and C. R. Wolf, "Submicron imaging with a planar silver lens," Appl. Phys. Lett. 84, 4403-4405 (2004) [CrossRef]
S. Brand, R. A. Abram and M. A. Kaliteevski, "Complex photonic bandstructure and effective plasma frequency of a two-dimensional array of metal rods," Phys. Rev. B 75, 035102, (2007)
A. J. Gallant, M. A. Kaliteevski, D. Wood, M. C. Petty, .R A. Abram, S. Brand, G. P. Swift, D. A. Zeze and J. M. Chamberlain, "Passband filters for terahertz radiation based on dual metallic photonic structures," Appl. Phys. Lett. 91, 161115 (2007) [CrossRef]
Dispersion relations have been calculated using a complex photonic bandstructure method described in [11], and finite difference time domain software OMNISIM© has been employed to simulate the beam propagation. Due to minor convergence problems the frequency scales for the two methods are slightly different. In order to provide matching of the two scales, the bandstructure has been renormalized to provide matching of the frequencies at the J point which can be determined for both calculation techniques.
H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. Tamamura, T. Sato, and S. Kawakami, Appl. Phys. Lett. 74, 1212-1214 (1999) [CrossRef]
D. N. Chigrin, S. Enoch, C. M. Sotomayor Torres and G. Tayeb, "Self-guiding in two-dimensional photonic crystals," Opt. Express 11, 1203-1211 (2003). [CrossRef] [PubMed]

Appl. Phys. Lett.

D. O. S. Melville, R. J .Blaikie, and C. R. Wolf, "Submicron imaging with a planar silver lens," Appl. Phys. Lett. 84, 4403-4405 (2004) [CrossRef]
A. J. Gallant, M. A. Kaliteevski, D. Wood, M. C. Petty, .R A. Abram, S. Brand, G. P. Swift, D. A. Zeze and J. M. Chamberlain, "Passband filters for terahertz radiation based on dual metallic photonic structures," Appl. Phys. Lett. 91, 161115 (2007) [CrossRef]
H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. Tamamura, T. Sato, and S. Kawakami, Appl. Phys. Lett. 74, 1212-1214 (1999) [CrossRef]

Nature

K. Dholakia, "Against the spread of the light," Nature 451, 413 (2008) [CrossRef] [PubMed]
V. Patanjali V. Parimi, W. T. Lu, Plarenta Vodo and Srinivas Sridhar, "Photonic crystals: Imaging by flat lens using negative refraction," Nature 426, 404 (2003) [CrossRef]
E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou and C. M. Soukoulis, "Electromagnetic waves: Negative refraction by photonic crystals," Nature 423, 604-605 (2003) [CrossRef] [PubMed]

Opt. Express

Phys. Rev. B

S. Brand, R. A. Abram and M. A. Kaliteevski, "Complex photonic bandstructure and effective plasma frequency of a two-dimensional array of metal rods," Phys. Rev. B 75, 035102, (2007)

Phys. Rev. Lett.

J .Durnin, J. J. Mieceli, and J. H. Eberly, "Diffraction-free beams," Phys. Rev. Lett. 58, 1499-1501 (1987). [CrossRef] [PubMed]

Phys.Rev. Lett.

J. B. Pendry, "Negative refraction makes a perfect lens," Phys.Rev. Lett. 85, 3966-3969 (2000) [CrossRef] [PubMed]

Phys.Rev.Lett.

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Observation of accelerating Airy beams," Phys.Rev.Lett. 99, 213901 (2007) [CrossRef]

Science

R. A. Shelby, D. R. Smith, and S. Shultz, "Experimental verification of a negative index of refraction," Science 292, 77-79 (2002) [CrossRef]

Sov. Phys. Usp.

V. G. Veselago, "The electrodynamics of substances with simultaneously negative values of ? and ?," Sov. Phys. Usp. 10, 509-514 (1968) [CrossRef]

Other

Dispersion relations have been calculated using a complex photonic bandstructure method described in [11], and finite difference time domain software OMNISIM© has been employed to simulate the beam propagation. Due to minor convergence problems the frequency scales for the two methods are slightly different. In order to provide matching of the two scales, the bandstructure has been renormalized to provide matching of the frequencies at the J point which can be determined for both calculation techniques.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper