## Abstract

We predict that a liquid crystal/silver nanoparticles mixture can be designed so that its effective ordinary and extraordinary permittivities have, in a frequency range, real parts of different signs. We exploit this result to design an optical switch obtained by sandwiching a few hundred nanometers thick slab of the proposed mixture between two silica layers. By resorting to full-wave simulations, we show that, by varying the direction of an externally applied electric field, the transmissivity of the structure can be switched between 0.02 and 0.4 at a wavelength close to the frequency range where the medium is indefinite. The device functionality physically stems from the fact that the orientation of the hyperbola characterizing extraordinary waves within the indefinite medium follows the applied electric field direction and therefore, if the hyperbola asymptote is nearly normal to the slab, full switch between evanescent and homogeneous propagating waves can be achieved within the medium.

©2011 Optical Society of America

## 1. Introduction

Within the rapidly growing field of metamaterials, indefinite media [1], or extremely anisotropic metamaterials characterized by principal permittivities having different signs, are attracting a considerable research interest mainly for the novel potentialities they provide for achieving efficient optical steering and manipulation. The hyperbolic dispersion relation characterizing extraordinary plane waves through indefinite (or hyperbolic) media is the main physical ingredient leading to unusual optical effects. Examples are negative refraction [2–4], which is due to the fact that the hyperbola curvature is opposite to that of the circumference of standard media, and hyperlensing [5–8], arising from to the medium ability of converting vacuum evanescent waves into homogeneous propagating waves [9] (as a consequence of the lack of evanescent waves along the hyperbola major axis). Extreme anisotropy also fundamentally affects the excitation of waves, impinging from vacuum, into an indefinite medium to the point that, by suitably changing the mutual geometric orientation between the hyperbola and the metamaterial interface, novel interesting effects as cancellation of reflection and transmission can be predicted [10]. Consequently, optical devices exploiting indefinite media have been proposed as, for instance, polarization beam splitters [11] and angular filters [12] which are based on the fact that transverse magnetic and transverse electric waves are characterized by different dispersive features (hyperbolic and standard behavior, respectively).

The condition of extreme anisotropy, where the principal permittivities have different signs, can be attained by resorting to suitable metal-dielectric nano-composites with a geometrical structure exhibiting at least one privileged direction. Exploiting the effective medium approach, it is possible to design the composite to show indefinite permittivity in a spectral range whose wavelengths are much greater than the nano-constituents size. Two relevant examples of such structured media are nanowire composites [13–15] (typically silver nanowires embedded in an allumina membrane) and layered media [16–19] (obtained by alternating metal and dielectric layers) whose privileged directions are, evidently, the nanowire and the stacking direction, respectively. Recently, it has been proposed that indefinite permittivity can even be observed in anisotropic natural materials as, for example, graphite [20] and strong anisotropic uniaxial crystals [21] which are shown to be indefinite in the ultraviolet (due to the hybrid electronic transitions) and the infrared region (due to the strong anisotropy of polar lattice vibrations), respectively.

Dynamic light manipulation is a key ingredient of photonics and therefore investigating tunable metamaterials is an essential task. Basically, it has been proposed that if the metamaterial constituents are immersed into a liquid crystal [22, 23] or the liquid crystal is used as one of the metamaterial constituents [24–26], the overall refractive index can be controlled through an external electric field, even allowing the full switch of the refractive index sign. Remarkable is the situation where the liquid crystal [27, 28] hosts dispersed structured nanoparticles since the negative-zero-positive transition of the refractive index sign can be more simply achieved by suitably designing the nanoparticle properties (filling fraction, radius, etc.). Due to its efficiency and manufacturing simplicity, the liquid crystal/nanoparticles mixture has also been investigated by means of refined approaches accounting for the nonlocal and spatial inhomogeneities of the refractive index [29]. Although the present research interest in tunable metamaterials is rapidly growing, it should be stressed that, to the best of our knowledge, an externally driven metamaterial with indefinite permittivity has not been considered yet, not even in the wide and versatile research field of liquid crystal. An example of such an indefinite and tunable metamaterial is the main achievement of the present paper.

In this paper we consider a liquid crystal/silver nanoparticles mixture and we show
that, by suitably choosing the volume filling fraction of the particles, the
anisotropic effective medium theory predicts the existence of a spectral range where
the permittivity tensor is indefinite. Therefore, to the best of our knowledge, we
theoretically discuss the first example of indefinite optical metamaterial which 1)
is obtained by dispersing metal nanoparticles within an anisotropic matrix and 2) is
tunable since, due to the presence of the liquid crystal, an external electric field
is able to change its optical response. These achievements allow us to consider an
active optical device obtained by sandwiching a 300 nm thick slab of the considered
tunable indefinite medium between two 70 nm thick silica layers and operating in the
presence of an externally applied electric field. Exploiting the effective medium
theory, we show that, at optical wavelengths where the inner medium is indefinite,
the slab is completely opaque if the electric field is orthogonal to its interfaces
whereas it becomes gradually more transparent if the angle is increased. Such a
behavior is physically due to the fact that the orientation of the hyperbola
characterizing the dispersion relation of extraordinary waves through indefinite
medium follows the direction of the applied electric field and therefore, when the
hyperbola asymptote is close to be normal to the device interface, a switch between
evanescent and propagating waves within the medium is achieved thus producing the
change of the slab transmissivity. In order to prove that the considered medium does
actually exhibit indefinite permittivity and that the considered device operates the
predicted optical switch, we have performed full-wave (FDTD) simulations and
evaluated the device transmissivity as a function of the wavelength and the applied
electric field direction. We show that a spectral region where the device operates
the switch exists and it is very close to that predicted by the effective medium
theory and, most importantly, that, at a specific wavelength, the actual device
transmissivity can be switched between 0.02 and 0.4 through the externally applied
electric field. It should be stressed that such a variation of the transmissivity
could not be achieved through conventional liquid crystal electro-optical switches
of the same thickness *L* = 300 nm since such length is
too small, at optical frequencies, to produce the required modal phase
difference.

## 2. Liquid Crystal/silver nanoparticles mixture as a tunable indefinite medium

Consider a slab of thickness *L* with its interfaces orthogonal to the
*z*-axis, filled with a nematic liquid crystal (LC) and
sandwiched between to silica layers (see. Fig.1). If the surface region of the LC is weakly anchored to the silica
walls (condition that can be achieved through suitable surface [30] and/or bulk [31] preparation methods) and if the
initial homeotropic director of the LC molecules are along the *z*
axis, an externally applied electric field **E**
* _{S}*, lying in the

*xz*plane and forming an angle

*ϕ*with the

*z*-axis, is able to align the LC homeotropic director along the unit vector

**n̂**= sin

*ϕ*

**ê**

*+ cos*

_{x}*ϕ*

**ê**

*. As a consequence the LC optically behaves as a uniaxially anisotropic crystal (with optical axis along the direction*

_{z}**n̂**) and its permittivity tensor is ${\varepsilon}^{(LC)}\hspace{0.17em}=\hspace{0.17em}R\hspace{0.17em}\text{diag}\hspace{0.17em}\left({\varepsilon}_{o}^{(LC)},\hspace{0.17em}{\varepsilon}_{o}^{(LC)},\hspace{0.17em}{\varepsilon}_{e}^{(LC)}\right)\hspace{0.17em}{R}^{-1}$ where ${\varepsilon}_{o}^{(LC)}$ and ${\varepsilon}_{e}^{(LC)}$ are the ordinary and extraordinary permittivities, respectively and

*y*-axis rotation matrix. The slab hosts silver nanoparticles of radius

*r*= 10 nm which are dispersed within the LC with the volume filling fraction

*f*. If both the the radiation wavelength

*λ*and the slab thickness

*L*are much greater than the nanoparticles radius (

*λ*≫

*r*and

*L*≫

*r*), the overall optical response of the slab can be described by means of the effective medium approach suitably extended to encompass anisotropic constituents [32]. According to this approach, since silver nanoparticles are optically isotropic with scalar permittivity

*ɛ*

^{(Ag)}, the overall effective permittivity tensor is ${\varepsilon}^{(eff)}\hspace{0.17em}=\hspace{0.17em}R\hspace{0.17em}\text{diag}\hspace{0.17em}\left({\varepsilon}_{o}^{(eff)},\hspace{0.17em}{\varepsilon}_{o}^{(eff)},\hspace{0.17em}{\varepsilon}_{e}^{(eff)}\right)\hspace{0.17em}{R}^{-1}$ so that the homogeneous effective medium is a uniaxial crystal with optical axis along the direction

**n̂**and its ordinary (

*j*=

*o*) and extraordinary (

*j*=

*e*) effective permittivities are

*N*are

_{j}*N*= 1/3), the depolarization factors accounting for the effect of the anisotropy. Note that, in analogy to the isotropic situation, the nanoparticles radius does not appear in Eq.(2) so that the mixture design have to be performed by tuning the volume filling fraction

_{j}*f*.

In order to discuss a feasible example of mixture, we have chosen SG-1 as nematic LC
[33] due to its
high birefringence and we have reported its ordinary and extraordinary
permittivities in Fig.2(a) as functions of
the vacuum wavelength *λ* in the spectral range between
410 and 430 nm. For the silver permittivity we have used the Drude model ${\varepsilon}^{(Ag)}\hspace{0.17em}=\hspace{0.17em}{\varepsilon}_{\infty}\hspace{0.17em}-\hspace{0.17em}{\omega}_{p}^{2}/({\omega}^{2}\hspace{0.17em}+\hspace{0.17em}i\mathrm{\Gamma}\omega )$ with *ɛ*
_{∞}
= 4.56, *ω _{p}* = 1.38
· 10

^{16}s

^{−1}and Γ = 0.1 · 10

^{15}s

^{−1}[34] and we have plotted the real and imaginary parts of

*ɛ*

^{(Ag)}in Fig.2(b) for

*λ*between 410 and 430 nm. Note that we have used a scattering rate Γ which is an order of magnitude higher than that of bulk silver in order to accounting for the surface scattering (size effect) occurring within the considered nanometric sized particles. After setting

*f*= 0.05 for the nanoparticles volume filling fraction, we have evaluated from Eq.(2) the principal permittivities of the effective medium and we have plotted their real and imaginary parts in Fig.2(c). Note that for wavelengths between

*λ*

_{1}= 418 nm and

*λ*

_{2}= 427.5 nm, $\text{Re}\hspace{0.17em}\left({\varepsilon}_{e}^{(eff)}\right)\hspace{0.17em}>\hspace{0.17em}0$ and $\text{Re}\hspace{0.17em}\left({\varepsilon}_{o}^{(eff)}\right)\hspace{0.17em}<\hspace{0.17em}0$ so that the effective medium has indefinite permittivity. Combining this result with the fact that the direction of the optic axis can be varied through the externally applied electric field, we conclude that the proposed LC/silver nanoparticles mixture behaves, within the limits of validity of the effective medium theory, as a tunable hyperbolic medium.

The above described design of the anisotropic effective medium has been performed for
the specific value *f* = 0.05 of the nanoparticles volume
filling fraction. In order to complete the analysis and to elucidate the dependence
of the whole phenomenology on the parameter *f* we have evaluated the
wavelengths *λ*
_{1}(*f*) and
*λ*
_{2}(*f*) characterizing
the spectral region *λ*
_{1}(*f*)
< *λ* <
*λ*
_{2}(*f*) where the medium
is indefinite for 0.02 < *f* < 0.1 and
plotted them in Fig.3. Note that, a tolerance
Δ*f* = 0.005 for the volume filling
fraction produces the uncertainties
Δ*λ*
_{1} = 2.3 nm and
Δ*λ*
_{1} = 1.4 nm,
thus proving that the mixture design is not critical.

## 3. Optical switch

The intriguing optical properties of hyperbolic media combined with the tunability,
allow the above described LC/nanoparticles mixture to be used for conceiving optical
devices. As an example, consider the setup of Fig.1 where the LC slab with dispersed nanoparticles has thickness
*L* = 300 nm and it is sandwiched between two
*d* = 70 nm thick silica layers (of permittivity
*ɛ*
^{(S)} =
2.088). If a monochromatic plane wave is made to impinge orthogonally, from vacuum,
onto the device interface (as sketched in Fig.1), it is evident that we are considering an active device since its
optical transmissivity can be altered by varying the direction of the externally
applied electric field **E**
* _{S}*. Using the effective medium approach discussed in Sec.2 to model
LC/nanoparticles mixture, the device transmissivity is easily evaluated and it is
reported in Fig.4(a), as a function of both
the vacuum wavelength

*λ*and the angle

*ϕ*defining the direction of the applied electric field. For clarity purposes, in Fig.4(a) we have reported two transparent gray planes corresponding to the two wavelengths 418 and 427.5 nm between which the effective medium predicts that the mixture is hyperbolic. The most striking feature of Fig.4(a) is the existence of three relevant spectral regions, I)

*λ*< 418 nm, II) 418 nm <

*λ*< 427.5 nm and III)

*λ*> 427.5 nm, where the transmissivity is characterized by qualitatively different behaviors. Specifically within regions I and III, the slab is almost transparent and completely opaque, respectively whereas, within region II the transmissivity strongly depends on

*ϕ*to the point that

*T*≃ 0 for

*ϕ*≃ 0 and, for increasing

*ϕ*, it gradually increases eventually reaching its maximum (dependent on

*λ*) for

*ϕ*=

*π*/2 (see Fig.4(b) where the level plot of the transmissivity of region II is plotted). Therefore, within region II, the considered structure can be regarded as a very thin optical switch which, being driven by the external electric field, is able to completely stop or allow the transmission of the light flow. Physically, such an overall behavior is readily grasped by noting that within region I and III (see Fig.2(c)) the effective medium behaves as a uniaxial crystal and an anisotropic conductor, respectively and therefore, within the same two regions, the plane wave impinging from vacuum excites homogeneous and evanescent waves, respectively. On the other hand, within region II the effective medium is indefinite so that, depending on the value of

*ϕ*, the normally impinging plane wave incoming from vacuum is able to excite either homogeneous or evanescent waves. In order to clarify this crucial point, we start from the dispersion relation of extraordinary plane waves

*k′*and

_{x}*k′*are the components of the wave vector along the principal axis of the uniaxial crystal and

_{z}*k*

_{0}= 2

*π/λ*. Since the plane wave from vacuum normally impinges onto the device interface (see Fig.1), conservation of momentum requires that the wave vector of waves excited within the anisotropic slab is

**K**=

*K*

**ê**

*, i.e. it is along the*

_{z}*z*-axis. Therefore the wavevector components along the principal axis are

*k′*= –

_{x}*K*sin

*ϕ*and

*k′*=

_{z}*K*cos

*ϕ*which, inserted into Eq.(4) yield

*ϕ*on the excited waves is particularly evident. Even though we are not neglecting absorption (see the imaginary parts of the permittivity in Fig.2(c)), it is convenient here to regard the principal permittivities as real quantities with ${\varepsilon}_{o}^{(eff)}\hspace{0.17em}<\hspace{0.17em}0$ and ${\varepsilon}_{e}^{(eff)}\hspace{0.17em}>\hspace{0.17em}0$. Therefore, from Eq.(5), we readily obtain that

*K*(0) is imaginary (evanescent waves) whereas

*K*(

*π*/2) is real (homogeneous waves), the switch between the two kinds of waves occurring at the angle $\phi \hspace{0.17em}=\hspace{0.17em}\text{arctan}\hspace{0.17em}\sqrt{-{\varepsilon}_{e}^{(eff)}/{\varepsilon}_{0}^{(eff)}}$. This explains the dependence of

*T*on

*ϕ*within region II of Fig.4 where, as discussed above, a realistic situation is reported in which absorption is not neglected (the only effects of absorption being the smoothing of the dependence of

*T*on

*ϕ*and the lowering of values of

*T*).

The just described mechanism is pictorially sketched in Fig.5 where we have reported the transmissivity
*T* (*ϕ*) for
*λ* = 420 nm extracted from Fig.4(a). In the insets *A*) and
*B*) of Fig.5 we have plotted
the hyperbolas of Eq.(4) in the two
relevant situations (A) *ϕ* = 20°
and (B) *ϕ* = 80° (regarding the
permittivities as real quantities). The direction of the applied electric field is
able to physically rotate the medium hyperbola and, as a consequence, in the two
situations *ϕ* = 20° and
*ϕ* = 80° the hyperbola does
not intersect and intersects the *z*-axis, respectively, thus
correspondingly allowing the normally impinging vacuum wave to excite evanescent and
homogeneous waves. The switch between the two kinds of waves occurs when the
hyperbola asymptote is nearly normal to the slab interface and this occurs for
*ϕ* close to $\text{arctan}\hspace{0.17em}\sqrt{-{\varepsilon}_{e}^{(eff)}/{\varepsilon}_{0}^{(eff)}}$.

## 4. Full wave simulations

In order to check the homogenization of the LC/nanoparticles mixture and the device
switching functionality, we have performed 3D full-wave simulations by means of a
finite-element solver, using the numerical values of the parameters (constituent
permittivities, nanoparticles radius and volume filling fraction, layers
thicknesses, etc.) considered so far. In Fig.6(a) we have plotted the resulting transmissivity *T* as
a function of the vacuum wavelength *λ* and the angle
*ϕ*. It is worth noting that, even though the surface
of Fig.6(a) does not strictly coincide with
that of Fig.4(a), it is however evident that
the actual structure is almost transparent and practically opaque for wavelengths
*λ* smaller than 430 nm and greater than 445 nm,
respectively. In addition, in the spectral range between 430 and 445 nm the
transmissivity strongly depends on the angle *ϕ* and, at
a given *λ*, *T*
(*ϕ*) is a monotonically increasing function of
*ϕ* (see Fig.6(b)). As a consequence, even though the spectral ranges are slightly
different from those predicted by the effective medium approach, the sample behavior
changes, for increasing radiation wavelengths, from that of a dielectric uniaxial
crystal to that of an anisotropic conductor with an intermediate spectral region
where the medium is hyperbolic. In Fig.6(b)
we have plotted the transmissivity of the structure as a function of
*ϕ* for three relevant wavelengths and the switching
functionality is particularly evident. Moreover, since the transmissivity spans the
range between 0.02 and 0.4 for *λ* = 425 nm,
we conclude that the performance of the actual switching device is better than that
predicted by the effective medium theory (whose transmissivity range is between 0
and 0.15).

The discrepancies between the predictions of the homogenized effective medium theory
and those of the full-wave simulations are a consequence of the fact that, for the
chosen slab thickness *L* = 300 nm, the edge effect due
to the nanoparticles close to the slab interface can not be strictly neglected. The
derivation of the Maxwell Garnett mixing rule (whose anisotropic extension has been
used in this paper) is based on the evaluation of the polarizability of an isolated
nanoparticle experiencing an external field which is uniform within its volume. For
particle belonging to the mixture bulk, these requirements are fulfilled if the
nanoparticle size is much smaller than the wavelength and if the nanoparticles
volume filling fraction is very small. However, for nanoparticles close to the
sample interface the above polarizability evaluation is wrong since it does not take
into account the additional field existing within the nanoparticle volume which is
due to surface polarization charges (i.e. to the discontinuity of the
electromagnetic material properties at the interface) [35]. Since the corrections to the
permittivity scale as (*r/z*)^{3} (*z* being
the distance from the interface) the nanoparticles responsible for the edge effect
are those whose distance from the interface is smaller than 40 or 50 nm (four or
five times the radius *r* = 10 nm), so that we conclude
that roughly a fourth of the considered slab volume can not be described through the
Maxwell Garnett approach.

## 5. Conclusions

We have proposed a novel way for synthesizing indefinite media by dispersing metal nanoparticle within a highly birefringent liquid crystal. In addition to its simplicity and feasibility, the proposed mixture is also tunable since an externally applied electric field is able to change the overall effective permittivity tensor. Tunability of an indefinite medium is a very important property since, as we have shown, the ability of altering the orientation of the extraordinary wave hyperbola allows one to literally choose whether the electromagnetic field within the medium (as due to a give external excitation) has to be composed of propagating homogeneous or evanescent waves. We have exploited such a property characterizing tunable indefinite media to design a nanometric thick optical switch and investigated its functionality both through the effective medium approach and through full wave simulations. The effective medium approach has the sufficient simplicity to highlight the basic physical mechanisms supporting the device functionality. On the other hand, full wave simulations yield results which are slightly different from those predicted by the effective medium approach and the discrepancy are related to the very small device thickness we decided to consider. However, from the point of view of the device functionality, the results of the full wave simulations are even better than those of the effective medium approach.

## References and links

**1. **D. R. Smith and D. Schurig, “Electromagnetic Wave Propagation in
Media with Indefinite Permittivity and Permeability
Tensors,” Phys. Rev. Lett. **90**, 077405 (2003). [CrossRef] [PubMed]

**2. **D. R. Smith, P. Kolinko, and D. Schurig, “Negative refraction in indefinite
media,” J. Opt. Soc. Am. B **21**, 1032–1043
(2004). [CrossRef]

**3. **D. R. Smith, D. Schurig, J. J. Mock, P. Kolinko, and P. Rye, “Partial focusing of radiation by a
slab of indefinite media,” Appl. Phys.
Lett. **84**, 2244–2246
(2004). [CrossRef]

**4. **A. Fang, T. Koschny, and C. M. Soukoulis, “Optical anisotropic metamaterials:
Negative refraction and focusing,” Phys.
Rev. B **79**, 245127 (2009). [CrossRef]

**5. **Z. Jacob, L. V. Alekseyev, and E. Narimanov, “Optical Hyperlens: Far-field
imaging beyond the diffraction limit,”
Opt. Express **14**, 8247–8256
(2006). [CrossRef] [PubMed]

**6. **H. Lee, Z. Liu, Y. Xiong, C. Sun, and X. Zhang, “Development of optical hyperlens
for imaging below the diffraction limit,”
Opt. Express **15**, 15886–15891
(2007). [CrossRef] [PubMed]

**7. **J. Yao, K. T. Tsai, Y. Wang, Z. Liu, G. Bartal, Y. L. Wang, and X. Zhang, “Imaging visible light using
anisotropic metamaterial slab lens,”
Opt. Express **17**, 22380–22385
(2009). [CrossRef]

**8. **B. D. F. Casse, W. T. Lu, Y. J. Huang, E. Gultepe, L. Menon, and S. Sridhara, “Super-resolution imaging using a
three-dimensional metamaterials nanolens,”
Appl. Phys. Lett. **96**, 023114 (2010). [CrossRef]

**9. **Z. Liu, Z. Liang, X. Jiang, X. Hu, X. Li, and J. Zi, “Hyper-interface, the bridge between
radiative wave and evanescent wave,”
Appl. Phys. Lett. **96**, 113507 (2010). [CrossRef]

**10. **J. Yang, X. Hu, X. Li, Z. Liu, X. Jiang, and J. Zi, “Cancellation of reflection and
transmission at metamaterial surfaces,”
Opt. Lett. **35**, 16–18
(2010). [CrossRef] [PubMed]

**11. **J. Zhao, Yan Chen, and Y. Feng, “Polarization beam splitting through
an anisotropic metamaterial slab realized by a layered metal-dielectric
structure,” Appl. Phys. Lett. **92**, 071114 (2008). [CrossRef]

**12. **L. V. Alekseyev, E. E. Narimanov, T. Tumkur, H. Li, Yu. A. Barnakov, and M. A. Noginov, “Uniaxial epsilon-near-zero
metamaterial for angular filtering and polarization
control,” Appl. Phys. Lett. **97**, 131107 (2010). [CrossRef]

**13. **J. Elser, R. Wangberg, and V. A. Podolskiy, “Nanowire metamaterials with extreme
optical anisotropy,” Appl. Phys. Lett. **89**, 261102 (2006). [CrossRef]

**14. **Y. Liu, G. Bartal, and X. Zhang, “All-angle negative refraction and
imaging in a bulk medium made of metallic nanowires in the visible
region,” Opt. Express **16**, 15439–15448
(2008). [CrossRef] [PubMed]

**15. **M. A. Noginov, Yu. A. Barnakov, G. Zhu, T. Tumkur, H. Li, and E. E. Narimanov, “Bulk photonic metamaterial with
hyperbolic dispersion,” Appl. Phys.
Lett. **94**, 151105 (2009). [CrossRef]

**16. **J. Schilling, “Uniaxial metallo-dielectric
metamaterials with scalar positive
permeability,” Phys. Rev. E **74**, 046618 (2006). [CrossRef]

**17. **A. J. Hoffman, L. Alekseyev, S. S. Howard, K. J. Franz, D. Wasserman, V. A. Podolsky, E. E. Narimanov, D. L. Sivco, and C. Gmachl, “Negative refraction in
semiconductor metamaterials,” Nat.
Materials **6**, 946–950
(2007). [CrossRef]

**18. **Y. J. Huang, W. T. Lu, and S. Sridhar, “Nanowire waveguide made from
extremely anisotropic metamaterials,”
Phys. Rev. A **77**, 063836 (2008). [CrossRef]

**19. **D. Korobkin, B. Neuner, C. Fietz, N. Jegenyes, G. Ferro, and G. Shvets, “Measurements of the negative
refractive index of sub-diffraction waves propagating in an indefinite
permittivity medium,” Opt. Express **18**, 22734–22746
(2010). [CrossRef] [PubMed]

**20. **J. Sun, J. Zhou, B. Li, and Feiyu Kang, “Indefinite permittivity and
negative refraction in natural material:
Graphite,” Appl. Phys. Lett. **98**, 101901 (2011). [CrossRef]

**21. **R. Wang, J. Sun, and J. Zhoua, “Indefinite permittivity in uniaxial
single crystal at infrared frequency,”
Appl. Phys. Lett. **97**, 031912 (2010). [CrossRef]

**22. **M. V. Gorkunov and M. A. Osipov, “Tunability of wire-grid
metamaterial immersed into nematic liquid
crystal,” J. Appl. Phys. **103**, 036101 (2008). [CrossRef]

**23. **R. Pratibha, K. Park, I. I. Smalyukh, and W. Park, “Tunable optical metamaterial based
on liquid crystal-gold nanosphere composite,”
Opt. Express **17**, 19459–19469
(2009). [CrossRef] [PubMed]

**24. **X. Wang, D. Kwon, D. H. Werner, and I. C. Khoo, “Tunable optical negative-index
metamaterials employing anisotropic liquid
crystals,” Appl. Phys. Lett. **91**, 143122 (2007). [CrossRef]

**25. **D. H. Werner, D. Kwon, I. C. Khoo, A. V. Kildishev, and V. M. Shalaev, “Liquid crystal clad near-infrared
metamaterials with tunable negative-zero-positive refractive
indices,” Opt. Express **15**, 3342–3347
(2007). [CrossRef] [PubMed]

**26. **A. Minovich, D. N. Neshev, D. A. Powell, I. V. Shadrivov, and Y. S. Kivshar, “Tunable fishnet metamaterials
infiltrated by liquid crystals,” Appl.
Phys. Lett. **96**,
193103
(2010).

**27. **I. C. Khoo, D. H. Werner, X. Liang, and A. Diaz, “Nanosphere dispersed liquid
crystals for tunable negative-zero-positive index of refraction in the
optical and terahertz regimes,” Opt.
Lett. **31**, 2592–2594
(2006). [CrossRef] [PubMed]

**28. **I. C. Khoo, A. Diaz, J. Liou, M. V. Stinger, J. Huang, and Y. Ma, “Liquid Crystals Tunable Optical
Metamaterials,” IEEE J. Sel. Top. Quant.
Electron. **16**, 410–417
(2010). [CrossRef]

**29. **G. Pawlik, M. Jarema, W. Walasik, A. C. Mitus, and I. C. Khoo, “Field-induced inhomogeneous index
distribution of a nano-dispersed nematic liquid crystal metamaterial near
the Freedericksz transition: Monte Carlo
studies,” J. Opt. Soc. Am. B **27**, 567–576
(2010). [CrossRef]

**30. **J. S. Patel and H. Yokoyama, “Continuous anchoring transition in
liquid crystals,” Nature **362**525–527 (1993). [CrossRef]

**31. **G. P. Bryan-Brown, E. L. Wood, and I. C. Sage, “Weak surface anchoring of liquid
crystals,” Nature **399**338–340 (1999). [CrossRef]

**32. **A. Sihvola, “Mixing Rules with Complex
Dielectric Coefficients,” Subsurf. Sens.
Technol. and Appl. **1**, 393–415
(2000). [CrossRef]

**33. **S. Gauza, J. Li, S. T. Wu, A. Spadlo, R. Dabrowski, Y. N. Tzeng, and K. L. Cheng, “High birefringence and high
resistivity isothiocyanate-based nematic liquid crystal
mixtures,” Liq. Cryst. **32**, 10771085 (2005). [CrossRef]

**34. **E. D. Palik, *Handbook of Optical Constants*
(Academic, 1985).

**35. **R. Stepanyan, “Effective dielectric properties of
composite materials in the surface layer,”
Phys. Rev. B **67**, 073403 (2003). [CrossRef]