Tunable negative-index photonic crystals using colloidal magnetic fluids*
Geng Tao†, Wang Xin, Wang Yan, Dong Xiang-Mei
Engineering Research Center of Optical Instrument and System, Ministry of Education, Shanghai Key Laboratory of Modern Optical System, University of Shanghai for Science and Technology, Shanghai 200093, China

Corresponding author. E-mail: Tao_Geng@hotmail.com

*Project supported by the National Basic Research Program of China (Grant No. 2015CB352001), the Shanghai Rising-Star Program, China (Grant No. 12QA1402300), the China Scholarship Council (CSC) Program, and the Basic Research Program of Shanghai, China (Grant No. 14ZR1428500).

Abstract

The model of using colloidal magnetic fluid to build tunable negative-index photonic crystal is established. The effective permittivity ɛe and permeability μe of the two-dimensional photonic crystal are investigated in detail. For transverse magnetic polarization, both ɛe and μe exhibit a Lorentz-type anomalous dispersion, leading to a region where ɛe and μe are simultaneously negative. Then, considering a practical case, in which the thickness of photonic crystal is finite, the band structures for odd modes are calculated by the plane wave expansion method and the finite-difference time-domain method. The results suggest that reducing the external magnetic field strength or slab thickness will weaken the periodic modulation strength of the photonic crystal. Simulation results prove that the negative-index can be tuned by varying the external magnetic field strength or the slab thickness. The work presented in this paper gives a guideline for realizing the flat photonic crystal lens with tunable properties at optical frequencies, which may have potential applications in tunable near-field imaging systems.

PACS: 42.70.Qs; 78.20.Bh; 78.20.Ci
Keyword: tunable negative-index; photonic crystals; colloidal magnetic fluids
1. Introduction

The negative-index material (NIM) was first theoretically proposed by Veselago[1] in 1968. The NIM, which is proposed as a perfect lens of unlimited resolution, has attracted a great deal of attention from researchers.[2] Metamaterials and photonic crystals (PCs) have been proposed to construct NIM in the last ten years.[36] However, most NIM constructions usually have a narrow operating bandwidth and the effective index is not tunable. In recent years, some researchers used magnetic materials to construct tunable NIM at microwave frequencies.[79] The negative index can be realized based on the ferromagnetic resonance effect and the resonance frequency can be tuned by varying the external magnetic field strength. Unfortunately, the ferromagnetic resonance effect is suppressed dramatically as the operating frequency increases from microwave frequency to optical frequency. It means that the tunable NIM, which is based on the ferromagnetic resonance effect, cannot be used at optical frequency.

On the other hand, researchers have developed many selfassembled PCs to provide a photonic response in the visible regime.[10, 11] Magnetic fluid, as a kind of easy-made mature nanoscale magnetic material, has proved to be an excellent candidate for fabricating self-assembled PCs, [11, 12] especially the band-tunable PC with a fast and reversible response to external magnetic fields. In Ref.  [13], tunable negative refraction at optical frequencies was theoretically demonstrated in ferrofluids containing Fe3O4 nanoparticles coated by an isotropic Ag shell. The proposed tunable metamaterial can exhibit negative refraction due to the hyperbolic dispersion.[14] Thus, it is hard to use in the imaging system because of the strong anisotropy. In addition, the Ag shell results in the loss significantly.

In this paper, the proposed work focuses on the tunability of the negative refractive index via theoretical simulations based on the experimental data in the literature. The results show that the tunable negative-index PC can be realized easily by using general colloidal magnetic fluid in the optical regime, which may be desirable for tunable near-field imaging systems.

2. Model

Under an external magnetic field H, which is perpendicular to the magnetic fluid film surface, the magnetic nanoparticles are agglomerated due to the attraction among particles. The magnetic columns will appear when H is beyond a critical value, and then the fluid becomes a two-phase system (solid and liquid phases). Structural patterns will evolve from a onedimensional ordered array to a two-dimensional (2D) ordered array as H increases. Finally, a two-dimensional PC slab can be formed[15, 16] as shown in Fig.  1. The self-assembled PC usually has a triangular lattice structure due to the condition of minimum total energy.[17] Strictly speaking, it is a quasi-twodimensional PC because of the finite thickness of the sample in practice.

Fig.  1. Schematic diagram of the PC slab formed in a magnetic fluid under externally perpendicular magnetic field and confined by the substrate and superstratum along the direction of the magnetic field. a and h are the periods and heights of the magnetic columns, respectively.

In this work, a set of experimental data in Refs.  [18] and [19] for typical water-based Fe3O4 magnetic fluid samples are referred to. In the self-assembled process of PC, the permittivity of magnetic columns is close to that in the solid state and does not vary with the increase of H. However, the permittivity of the liquid phase decreases with the increase of H due to more magnetic nanoparticles aggregating into magnetic columns. Nevertheless, it is hard to measure the index of the liquid phase directly. For calculating the band structure, the relationship is needed between the concentrations of nanoparticles in the liquid phase and H. In practice, the relationship can be derived by measuring the field-dependent magnetization M of the liquid phase, because M is dependent on the number of magnetic nanoparticles per unit volume in the liquid phase. The relationship between M and H can be measured directly. The empirical formula can be expressed as[19]

Here, ɛ liq is the permittivity of the liquid phase and the magnetic permeability is 1 for optical frequency. Thus, according to Eq.  (1), the field-dependent index of liquid phase nliq can be calculated by using the experimental data in Ref.  [19].

The area ratio Acol/A is usually used to express the density of magnetic columns in practice and Acol is the area occupied by the magnetic columns with an area of A. For triangular structural patterns, the area ratio is expressed as

where r is the radius of the cross-section of the magnetic column and a is the perpendicular distance between the centers of adjacent magnetic columns, namely, the lattice constant. Then, according to Eq.  (2), the value of r/a can be calculated through by using the experimental data in Ref.  [18]. In general, according to the previous analysis, listed in Table  1 are the calculation results of fielddependent nliq and r/a under different magnetic field strengths by using the experimental data in Refs.  [18] and [19]. In the near-infrared and visible regime, the index of the magnetic column is a constant of 2.2 for Fe3O4.[20]

Table 1. The calculated results of nliq and r/a under different magnetic field strengths in the magnetic fluid film. The unit 1  Oe= 79.5775  A/m.
3. Simulation results and discussion

The PC slab can be regarded as an ideal two-dimensional (2D) PC when the column height h is much larger than the lattice constant a (about 10 times). Figures  2(a) and 2(b) show the band structures for transverse magnetic (TM) polarization (electric field is parallel to the columns) and transverse electric (TE) polarization (magnetic field is parallel to the columns) under external magnetic field H = 250  Oe by the plane wave expansion (PWE) method. Here, the frequency is normalized as λ /a. A band gap between the first and second band is presented in TM mode, but it does not appear in the TE mode.

Fig.  2. Band structures of the 2D PC under external magnetic field 250  Oe for (a) TM polarization and (b) TE polarization.

The refraction in PC is more complicated than that in uniform medium. In PC, an incident wave can produce multiple refracted Bloch-wave orders, if there are multiple branches of dispersion contours at the same frequency. Therefore, the candidate PC should guarantee the absence of any high-order Bragg reflections for any angle of incidence. Furthermore, the periodic modulation effect in the candidate PC should be strong enough in order to make the equal frequency surface (EFS) keep a circular shape.[21] Figures  3(a) and 3(b) show the EFS contours for TM polarization and TE polarization at several relevant frequencies in the second band. It is clear that the EFS shape of the TM mode is closer to the isotropic case than that of TE modes. It indicates that the periodic modulation effect for TM polarization in the second band is stronger. The EFS shape of TE modes cannot keep a circular shape even near the Γ points around f = 0.6551. Therefore, the analysis of TM modes is conducted in detail. Figure  4 shows the first band gap of TM modes as a function of external magnetic field H. With H increasing, the gap width increases, but the mid-frequency of the gap decreases. The feature has also been demonstrated experimentally.[22, 23] The results show that the periodic modulation strength will be weakened as H decreases.

Fig.  3. EFS contours of the 2D PC under external magnetic field 250  Oe for (a) TM polarization and (b) TE polarization at several relevant frequencies in the second band. The gray line shows the first Brillouin zone.

Fig.  4. The first band gap of TM modes as a function of magnetic field. Square symbols denote the middle of the band gaps and the vertical bars mark the spans of the band gaps.

From Fig.  3(a), it is clear that EFS contours shrink as frequency increases. This means that the group velocities are opposite to the phase velocities. A developing approach of effective medium theory beyond the static regime[24, 25] has been employed to investigate the behavior of effective permittivity ɛ e and permeability μ e. Figure  5(b) shows the calculated results of ɛ e and μ e with H = 250  Oe for TM polarization. Both ɛ e and μ e exhibit a Lorentz-type anomalous dispersion. The resonant frequencies appear at values 0.389 and 0.369 of ɛ e and μ e respectively. The regions for ɛ e < 0 and μ e < 0 are located at 0.389 < f < 0.627 and 0.369 < f < 0.701, respectively. It indicates that the magnetic resonance has a larger resonance width than the electric resonance. In this work, since both ɛ e and μ e are a real function of f, the effective refractive index neff can be either real (ɛ e > 0 and μ e > 0 or ɛ e < 0 and μ e < 0) or imaginary (ɛ e > 0 and μ e < 0 or ɛ e < 0 and μ e > 0). When neff is real, the figure of merit (FOM), which is defined as − Re(neff)/Im(neff), [26, 27] is infinite. When neff is imaginary, FOM is zero. The band gap appears at the resonant frequencies as shown in Fig.  5(a). The fields are localized around individual dielectric cylinders with rather weak interactions between neighboring cells at the resonant frequencies. The localized modes were discussed in detail in Ref.  [27]. In the region of the second band, neff is negative due to ɛ e < 0 and μ e < 0 simultaneously as shown in Fig.  5(b).

Fig.  5. (a) Band structure, (b) effective permittivity ɛ e and permeability μ e of 2D PC with external magnetic field 200  Oe for TM polarization. In panel (a), the shadow area indicates the band gap.

Figure  6(a) shows the calculated effective indices neff as a function of magnetic field at f = 0.45 along the Γ M direction. The effective index increases from – 1 to – 0.964 as H increases from 150  Oe to 250  Oe. Figures  6(b) and 6(c) show the normalized average intensities over a period at H = 150  Oe and H = 250  Oe, respectively. A continuous wave point source with f = 0.45 is placed at a distance d0 = 0.65λ from the lower interface of the PC lens. The width of PC lens d is 1.34λ . The images are located at distances di_150  Oe = 0.44λ and di_250  Oe = 0.49λ from the upper interface, respectively. It indicates clearly that the effective refractive index of such a PC lens can be manipulated by the external magnetic field strength. For an ideal superlens with refractive index equal to – 1, d should be equal to d0 + di.[28] The mismatch between d and d0 + di is due to the anisotropy of the EFS. For example, the EFS shape of the PC lens at H = 250  Oe is not circular when 0.39 < f < 0.55 as shown in Fig.  3(a). Similarly, the full widths at half-maximum (FWHMs) of both images at H = 150  Oe (FWHM = 0.83λ ) and H = 250  Oe (FWHM = 0.63λ ) are both larger than λ /2 as shown in Fig.  6(d). The FWHM of the image at H = 250  Oe is smaller than that at H = 150  Oe due to a stronger periodic modulation effect. Fortunately, the resolution problem can be settled by fixing a subwavelength metallic grating before the PC lens. The grating can inhibit the zero-order diffraction and enhance the high-order diffraction to achieve super-resolution.[29] Figure  6(e) shows the transmission of light from a point source with f = 0.45 through the PC lens as a function of H. It can be seen that the transmittivities are more than 0.45 for all the cases.

Fig.  6. (a) Calculated effective indices of 2D PCs versus external magnetic field at f = 0.45 along the Γ M direction. The normalized average intensities of a point source with f = 0.45 located at 0.65λ from the PC lens at (b) H = 150  Oe and (c) H = 250  Oe. In panels  (b) and (c), the solid lines show the position of the PC lens. The z positions of the images are marked by the dashed lines. (d) The normalized average intensities along the x direction at the image plane at H = 150  Oe and H = 250  Oe. (e) The transmission of light from a point source with f = 0.45 through the PC lens as a function of external magnetic field.

When the column height h is of the same order with the lattice constant a, the band structure will be different from the ideal 2D case. The states in PC slabs are not purely TM nor TE polarized because of the absence of translational symmetry in the vertical direction (perpendicular to the plane of periodicity). Instead, the guided bands can be classified as odd or even modes. However, the odd and even modes have strong similarities to TM and TE modes, respectively, in two dimensions.[30]

For PC slab with finite thickness, there is a light cone in the band structure, which consists of all the states existing in the bulk background. When a guided band reaches the edge of the light cone, it extends infinitely far into the background, and cannot be used to permanently confine the light within the slab. In other words, only the bands below the light cone can be guided in the PC slab. The position of the light cone depends on the refractive indices of substrate and superstratum. It should be tried to use a low refractive index material (such as aerogel glass, whose index is nearly 1) as the cladding layer. Here, the index of the cladding layer is assumed to be 1. The band structures, which lie above in the light cone, are retained here, because these modes could be important for some applications such as the photonic crystal LED.[31] The band structures for PC slabs are calculated using both the finitedifference time-domain (FDTD) method and PWE method. Since PWE requires periodicity in all directions, the photonic crystal slab is already two-dimensionally periodic, and it can be imposed by a three-dimensional periodicity through a periodic sequence of slabs separated by a sufficient amount of background region. Because the guided modes are localized within the slab, the addition of periodic slabs at large intervals exerts little effect on their frequencies. However, this supercell approach inevitably results in spurious modes and subsequent spurious bands.

Figure  7 shows the band structure of PC slab with h = a at H = 250  Oe for odd modes as an example. The results of the FDTD method (empty circles) and PWE method (solid line) agree well with each other. It can be seen that spurious modes arise from the fake slabs of the supercell approach in the PWE method (solid line), and all the bands move to a slightly higher frequency compared with the case of the 2D PC (black dash line). Figure  8 shows the first band gap of odd modes as a function of H and h. The gap width decreases as either H or h decreases, on the contrary, the mid-frequency of the gap increases. The calculated results indicate that the decrease of h weakens the periodic modulation strength, thereby reducing the imaging quality of the PC slab lens.

Fig.  7. Band structures of the PC slab with h = a at H = 250  Oe for odd modes using the FDTD (empty circles) and PWE methods (solid line). Gray lines specify the light line. Dashed line shows the band structure of the ideal 2D PC at H = 250  Oe for TM modes using the PWE method for comparison.

Fig.  8. The first band gap of odd modes as a function of H and h. Here, H changes from 150  Oe to 250  Oe and h changes from 2a to 6a. Solid symbols denote the middle of the band gaps and the vertical bars mark the spans of the band gaps.

Figures  9(a)– 9(c) show the normalized average intensities over a period at H = 150  Oe with h = 3a, H = 250  Oe with h = 3a, and H = 250  Oe with h = 4a, respectively. The calculation conditions are the same as those in Figs.  6(b) and 6(c), except that the normalized frequency of the point source is 0.44. The images are located at di1 = 0.144λ , di2 = 0.281λ , and di3 = 0.349λ from the upper interface in Figs.  9(b)– 9(c), respectively. The results show that the effective refractive index of the PC lens can also be manipulated by h, besides H.

Fig.  9. Normalized average intensities of a point source with f = 0.44 located at 0.65λ from the PC slab lens at (a) H = 150  Oe with h = 3a, (b) H = 250  Oe with h = 3a, and (c) H = 250  Oe with h = 4a. In panels (a), (b), and (c), the solid lines show the position of the PC slab lens. The z positions of the images are marked by the dashed lines.

4. Conclusions

In this work, the model of a self-assembled PC slab lens is established by using colloidal magnetic fluids. The effective refractive index of the PC slab lens can be tuned easily by varying the external magnetic field strength. The results show that around frequency 0.45, the effective refractive index is nearly – 1 and the FWHM of the image is close to the diffraction limit. It may have great potential for tunable nearfield imaging systems. Considering the fact that the periodic modulation strength of the PC slab lens will decrease as the column height decreases, which leads to the degradation of image quality, a thick container filled with colloidal magnetic fluids should be used to form the self-assembled PC slab lens.

Reference
1 Veselago V G 1968 Sov. Phys. Usp. 10 509 DOI:10.1070/PU1968v010n04ABEH003699 [Cited within:1]
2 Jin L, Zhu Q Y, Fu Y Q and Yu W Q 2013 Chin. Phys. B 22 104101 DOI:10.1088/1674-1056/22/10/104101 [Cited within:1]
3 Xu T, Agrawal A, Abashin M, Chau K J and Lezec H J 2013 Nature 497 470 DOI:10.1038/nature12158 [Cited within:1]
4 Li Y N, Gu P F, Zhang J L, Li M Y and Liu X 2006 Acta Phys. Sin. 55 4918 (in Chinese) 118.145.16.217/magsci/article/article?id=17386984 [Cited within:1]
5 Zhang S, Fan W, Panoiu N C, Malloy K J, Osgood R M and Brueck S R J 2005 Phys. Rev. Lett. 95 137404 DOI:10.1103/PhysRevLett.95.137404 [Cited within:1]
6 Fang Y T and Shen T G 2005 Chin. Phys. 14 949 [Cited within:1]
7 Zhao H, Zhou J, Kang L and Zhao Q 2009 Opt. Express 17 13373 DOI:10.1364/OE.17.013373 [Cited within:1]
8 Zhao H, Li B, Zhou J, Kang L, Zhao Q and Li W 2011 Opt. Express 19 15679 DOI:10.1364/OE.19.015679 [Cited within:1]
9 Liu S, Chen W, Du J, Lin Z, Chui S T and Chan C T 2008 Phys. Rev. Lett. 101 157407 DOI:10.1103/PhysRevLett.101.157407 [Cited within:1]
10 Liu Z Q, Feng T H, Dai Q F, Wu L J and Lan S 2009 Chin. Phys. B 18 179 DOI:10.1088/1674-1056/18/1/028 [Cited within:1]
11 Ge J and Yin Y 2008 Adv. Mater. 20 3485 DOI:10.1002/adma.200800657 [Cited within:2]
12 Kim H, Ge J, Kim J, Choi S, Lee H, Park W, Yin Y and Kwon S 2009 Nat. Photon. 3 534 DOI:10.1038/nphoton.2009.141 [Cited within:1]
13 Gao Y, Huang J P, Liu Y M, Gao L, Yu K W and Zhang X 2010 Phys. Rev. Lett. 104 034501 DOI:10.1103/PhysRevLett.104.034501 [Cited within:1]
14 Fang A, Koschny T and Soukoulis C M 2009 Phys. Rev. B 79 245127 DOI:10.1103/PhysRevB.79.245127 [Cited within:1]
15 Yang S Y, Horng H E, Hong C Y, Yang H C, Chou M C, Pan C T and Chao Y H 2003 J. Appl. Phys. 93 345 [Cited within:1]
16 Horng H E, Yang S Y, Lee S L, Hong C Y and Yang H C 2001 Appl. Phys. Lett. 79 350 DOI:10.1063/1.1385796 [Cited within:1]
17 Wen W, Zhang L and Sheng P 2000 Phys. Rev. Lett. 85 5464 DOI:10.1103/PhysRevLett.85.5464 [Cited within:1]
18 Horng H E, Hong C Y, Yang S Y and Yang H C 2003 Appl. Phys. Lett. 82 2434 DOI:10.1063/1.1568147 [Cited within:3]
19 Yang S Y, Chieh J J, Horng H E, Hong C Y and Yang H C 2004 Appl. Phys. Lett. 84 5204 DOI:10.1063/1.1765201 [Cited within:4]
20 Buchenau U and Müller I 1972 Solid State Commun. 11 1291 DOI:10.1016/0038-1098(72)90845-9 [Cited within:1]
21 Notomi M 2000 Phys. Rev. B 62 10696 DOI:10.1103/PhysRevB.62.10696 [Cited within:1]
22 Ge J, Hu Y and Yin Y 2000 Angewand te Chemie 119 7572 [Cited within:1]
23 Yang S Y, Horng H E, Shiao Y T, Hong C Y and Yang H C 2006 J. Magn. Magn. Mater. 307 43 DOI:10.1016/j.jmmm.2006.03.038 [Cited within:1]
24 Wu Y, Li J, Zhang Z Q and Chan C T 2006 Phys. Rev. B 74 085111 DOI:10.1103/PhysRevB.74.085111 [Cited within:1]
25 Jin J, Liu S, Lin Z and Chui S T 2009 Phys. Rev. B 80 115101 DOI:10.1103/PhysRevB.80.115101 [Cited within:1]
26 Zhang S, Fan W, Panoiu N C, Malloy K J, Osgood R M and Brueck S R J 2006 Opt. Express 14 6778 DOI:10.1364/OE.14.006778 [Cited within:1]
27 Chern R L and Chen Y T 2009 Phys. Rev. B 80 075118 DOI:10.1103/PhysRevB.80.075118 [Cited within:2]
28 Pendry J B 2000 Phys. Rev. Lett. 85 3966 DOI:10.1103/PhysRevLett.85.3966 [Cited within:1]
29 Shi P, Huang K and Li Y 2012 Opt. Lett. 37 359 DOI:10.1364/OL.37.000359 [Cited within:1]
30 Johnson S G, Fan S, Villeneuve P R, Joannopoulos J D and Kolodziejski L A 1999 Phys. Rev. B 60 5751 DOI:10.1103/PhysRevB.60.5751 [Cited within:1]
31 Kim D H, Cho C O, Roh Y G, Jeon H, Park Y S, Cho J, Im J S, Sone C, Park Y, Choi W J and Park Q H 2005 Appl. Phys. Lett. 87 203508 DOI:10.1063/1.2132073 [Cited within:1]