Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Complex k band diagrams of 3D metamaterial/photonic crystals.

Open Access Open Access

Abstract

A finite element method (FEM) for solving a complex valued k(ω) vs. ω dispersion curve of a 3D metamaterial/photonic crystal system is presented. This 3D method is a generalization of a previously reported 2D eigenvalue method [Opt. Express 15, 9681 (2007)]. This method is particularly convenient for analyzing periodic systems containing dispersive (e.g., plasmonic) materials, for computing isofrequency surfaces in the k-space, and for calculating the decay length of the evanescent waves. Two specific examples are considered: a photonic crystal comprised of dielectric spheres and a plasmonic fishnet structure. Hybridization and avoided crossings between Mie resonances and propagating modes are numerically demonstrated. Negative index propagation of four electromagnetic modes distinguished by their symmetry is predicted for the plasmonic fishnets. By calculating the isofrequency contours, we also demonstrate that the fishnet structure is a hyperbolic medium.

©2011 Optical Society of America

1. Introduction

The numerical simulation of electromagnetic fields inside both metamaterial and photonic crystals is an important tool for analyzing these periodic structures. In particular, the eigenmodes of crystals, defined as freely propagating waves not coupled to external currents, are often of the most interest. The conventional method [2,3] of numerically solving for crystal eigenmodes is to define the geometry of the unit cell of the crystal of interest and the differential equation that the fields must obey in this geometry and then impose Bloch periodic boundary conditions. The partial differential equation problem is then discretized using one of the many standard methods (finite element, finite integral, finite difference, etc.) thereby turning it into an algebraic eigenvalue problem with a finite number of degrees of freedom and the frequency ω as the eigenvalue. This finite sized eigenvalue problem is then solved numerically. An important detail of this method is that the Bloch wavenumber k is chosen beforehand, and the frequency is then computed as a function of the wavenumber, yielding the dispersion curves ω = ω(k). This is the most commonly used method for calculating dispersion curves of the electromagnetic waves propagating in photonic crystals or in closely related metamaterial crystals.

There are however, many instances where it is more convenient to specify the frequency ω and solve for the wavenumber as a function of frequency: k = k(ω). At least five such instances can be identified. First, metamaterials often contain dispersive materials such as metals, where the dielectric function strongly depends on the frequency of the wave. In this case, the eigenfrequency problem becomes a nonlinear eigenvalue problem and must be solved iteratively [4]. In contrast, when solving for the wavenumber as a function of frequency, the resulting eigenvalue problem only needs to be solved once. Second, to describe the microscopic electromagnetic response of a metamaterial inclusion many researchers use tabulated data for the permittivity, such as the data found in Ref. [5]. These tabulated material parameters are functions of a real valued ω and should not be used for ω(k) eigenvalue simulations where ω is often complex valued. Third, it is often useful to solve for the wavenumber as the eigenvalue because of the information contained in the complex wavenumber including the decay lengths of the electromagnetic modes (either due to finite dissipation or because of the evanescent nature of the mode) and the figure of merit [6] of negative index modes. Fourth, in the majority of experiments electromagnetic fields inside metamaterial/photonic crystals are excited by external sources producing time-harmonic fields with real valued frequencies. A complex wavenumber eigenvalue simulation provides the correct field distribution in the photonic crystal relevant to such an experiment. Finally, this approach provides a very natural way of calculating the so-called isofrequency surfaces corresponding to ω(k) = const, where ω is real. Isofrequency diagrams are fundamentally important for predicting wave refraction at the PC interfaces [7] and for calculating the density of states [8].

There are several previously published methods on calculating k(ω) dispersion curves including variations of the plane-wave expansion method [9, 10] and diagonalizing the crystal transfer matrix [11, 12]. A method of solving for complex wavenumber dispersion curves using the FEM has been proposed [13, 14] but only for 2D crystals. The benefits of a complex wavenumber 2D FEM are becoming better appreciated and use of this method is becoming more common [1519]. It is therefore timely to generalize this method of Complex Wavenumber Eigenvalue Simulations (CWES) from two to three dimension, which is the object of this paper. This 3D complex wavenumber eigenvalue simulation was recently used as part of a metamaterial homogenization procedure [20].

The basic theory behind solving for complex wavenumber eigenvalues using the FEM discretization is explained in Sec. 2. The underlying field equations for the electric and magnetic fields and the necessary boundary conditions are discussed. In Sec. 3 we apply this method to a 3D photonic crystal consisting of non-dispersive dielectric spheres. The photonic band structure for electromagnetic waves propagating both parallel and obliquely to the principal axes is calculated, and different types of modes (transverse and longitudinal) are identified. This simple model [21] can be solved with the commercial FEM software package COMSOL Multiphysics. In Sec. 4 we calculate the dispersion curves for a negative index fishnet metamaterial [22] and demonstrate the existence of four distinct negative index modes. We also calculate the two-dimensional isofrequency contours for the least damped transverse mode and demonstrate from the shape of the isofrequency contours that this transverse mode is hyperbolic.

2. The finite element eigenvalue problem

2.1. The field equation

In this section we present the FEM formulation for solving for the magnetic field. Electromagnetic wave propagation is described by the Maxwell equations which can be rearranged into a wave equation for either the electric field E or the magnetic field H. The wave equation for the magnetic field is

×(1ɛ×H)μω2c2H=0.
Here ɛ(x) and μ(x) are the microscopic permittivity and permeability of the metamaterial/photonic crystal of interest. Due to the periodic nature of the crystal both are are assumed to be scalar functions periodic in the crystal lattice. According to Bloch’s theorem [2, 3] the magnetic field can be represented as the product of a periodic function and an exponential factor
H(x)=u(x)exp[i(ωtkx)],
where ω is the frequency of the wave and k is the wavevector of the Bloch-Floquet wave. u(x) is a vector function which is periodic in the crystal lattice. By inserting Eq. (2) into Eq. (1) we obtain an equivalent field equation for u
k2ɛukɛ(ku)ik×(1ɛ×u)i×(1ɛk×u)+×(1ɛ×u)μω2c2u=0,
which can be interpreted as an eigenvalue problem and solved for the Bloch wavenumber k as the eigenvalue. The spatial profile of the eigenmode u(x) is also recovered providing the magnetic field profile according to Eq. (2) and the electric field profile according to
E(x)=1iɛω/c×H=1iɛω/c(ik×u+×u)exp[i(ωtkx)]

2.2. The finite element model

There are several commercial FEM software programs (COMSOL Multiphysics by COMSOL, HFSS by Ansys, Vector Fields Opera by Cobham Technical Services, etc.) that are available for modeling metamaterial/photonic crystals. These commercial software packages provide a convenient graphical user interface for defining a crystal’s geometry, meshing the computational domain, and visualizing the electromagnetic fields. This allows for models to be quickly developed and tested. Of the many commercial FEM codes available, the authors of this paper are only aware of one (COMSOL Multiphysics) that allows the user to specify the field equation to be solved. The simulation examples and results presented here were obtained using COMSOL. The FEM model used for the photonic crystal described in Sec. 3 is available for download [21]. In what follows, only the most essential features of the FEM approach are reviewed; more detailed treatments can be found elsewhere [2325].

The FEM is based on setting the integral of a so-called weak expression over the domain of interest to zero. Doing so ensures the field equation is satisfied and also creates boundary conditions. The weak expression corresponding to Eq. (3) is

FH(v,u)=k2ɛvu1ɛ(kv)(ku)i1ɛv[k×(×u)]i(×v)1ɛ(k×u)+(×v)1ɛ(×u)μω2c2vu,
where v(x) is a test function. When the integral of the weak expression over the unit cell Ω of the crystal is set to zero, integrating by parts gives us two separate integrals (Eq. (6)). The first integral enforces the field equation. The second integral is over the boundary of the domain and represents a natural boundary condition [23, 24],
0=Ωd3xFH(v,u)=Ωd3xv[1ɛk×(k×u)i1ɛk×(×u)i×(1ɛk×u)+×(1ɛ×u)μω2c2u]+ΩdAv[n^×1ɛ(ik×u+×u)],
where is the vector normal to the boundary. On an external boundary, the natural boundary condition enforced by the integral in Eq. (6) over the boundary Ω forces the expression × (−ik × u + ∇ × u)/ɛ to be equal to zero. Recalling Eq. (4) we note that this simply enforces the boundary condition × E = 0. This is known as the perfect electric conductor or PEC boundary condition. This natural boundary condition is the default if no other boundary condition is enforced. On an internal boundary within the unit cell the surface integrals over each side of the boundary must be equal to each other. The effect is that the tangential components of the electric field must be continuous across the internal boundary or × E + = × E where E + and E are the electric fields on opposite sides of the internal boundary.

The periodicity of u is enforced by imposing periodic boundary conditions on the exterior boundaries of the unit cell. In COMSOL, these periodic boundary conditions override the natural boundary condition. However, if a PEC boundary condition is desired inside the unit cell (e.g., on the surface of a metal inclusion) this can be accomplished by removing the subdomain representing the metal inclusion. Now only the exterior side of the metal boundary remains and the natural boundary condition forces the tangential electric fields to zero on this boundary.

If a perfect magnetic conductor or PMC boundary condition ( × H) is desired while solving for the magnetic field, this can be enforced with constraints [23] on the tangential magnetic field on the boundary.

In order to turn Eq. (6) into an eigenvalue problem, the three degrees of freedom that comprise the Bloch wavevector k must be reduced to one by restricting two degrees of freedom. This is accomplished by setting k = k 0 + λ n where λ will be the eigenvalue solved for, k 0 is an offset vector and n is a unit vector ( n · n = 1). Note that since there are many possible ways to restrict these two degrees of freedom, the resulting complex k(ω) vs. ω dispersion curves are in general not unique. The FEM turns the weak form and accompanying boundary conditions into an algebraic problem, in this case a quadratic eigenvalue problem [26]:

Au+λBu+λ2Cu=0,
where A, B and C are N × N matrices and u⃗ is an N × 1 solution vector. N is the number of degrees of freedom of the discretized system. Terms in the weak form (Eq. (5)) that are zero, first and second order in λ contribute to the A, B and C matrices respectively. This quadratic eigenvalue problem can be recast in the form of a generalized eigenvalue equation
(AB01)(uλu)=λ(0C10)(uλu).
When using COMSOL to solve the FEM problem, this linearization is performed automatically during the solution stage.

2.3. Electric field formulation

The previous discussion focused on solving for the magnetic field H or rather the periodic function u equal to the magnetic field with the exponential Bloch factor removed. This is especially convenient when an inclusion requires a PEC boundary condition since that is the natural boundary condition when solving for H. However, solving for the electric field is very similar to solving for the magnetic field. The wave equation for the electric field for a free wave is

×(1μ×E)ɛω2c2E=0.
As before, we replace the electric field with a periodic vector field times an exponential factor
E(x)=u(x)exp[i(ωtkx)],
producing the new field equation
k2μukμ(ku)ik×(1μ×u)i×(1μk×u)+×(1μ×u)ɛω2c2u=0.
The corresponding weak form for this field equation is
FE(v, u)=k2μvu1μ(kv)(ku)i1μv[k×(×u)]i(×v)1μ(k×u)+(×v)1μ(×u)ɛω2c2vu,
which is equivalent to Eq. (5) if ɛ and μ are interchanged. Integrating this weak form over the crystal unit cell by parts and setting its value to zero again gives produces two integrals, one enforcing the field equation and a surface integral enforcing the boundary condition × H = 0. Thus the PMC boundary condition is the natural boundary condition when solving for the electric field.

3. Example: dielectric photonic crystal

For a demonstration of the CWES method of calculating complex k dispersion curves we use a simple photonic crystal as an example. The unit cell, pictured in Fig. 1, is a cube with a dielectric sphere at the center surrounded by vacuum. The sphere has a radius of 0.3a, where a is the lattice constant of the cubic array, and a permittivity of ɛ = 5 – i0.01. Arrays of spheres as well as arrays of point dipoles are both systems that have been studied extensively by metamaterial researchers [2730].

 figure: Fig. 1

Fig. 1 Complex k dispersion curves and field profiles for eigenmodes of the photonic crystal pictured in the inset assuming k 0 = 0 and n = . (a) Real part of kx(ω) for a transversely polarized mode and a diagram of the crystal unit cell. (b) Imaginary part of kx(ω) for a transversely polarized mode and a field profile for the polarized transverse mode. There are two transverse modes, ŷ and electrically polarized, which are degenerate. (c) Real part of kx(ω) for two longitudinally polarized modes and a field profile for the magnetic longitudinal mode. (d) Imaginary part of kx(ω) for two longitudinally polarized modes and a field profile for the electric longitudinal mode. The longitudinal mode with the passband near ω = 4.5c/a is magnetically polarized in the direction and the longitudinal mode with the passband near ω = 5c/a is electrically polarized in the direction. The longitudinal modes correspond to Mie’s dipole resonances. For all dispersion curves the dotted lines are the result of a conventional ω(k) eigenvalue simulation. For all field profiles the frequency is ω = 2c/a with arrows representing Dy and Dz and color representing Dx.

Download Full Size | PPT Slide | PDF

As mentioned in Sec. 2.2, it is necessary to restrict two of the three degrees of freedom of the Bloch wavevector k. There are many possible ways to do this. As the first example, we calculate the dispersion curves corresponding to the propagation along a principal axis. To simulate this we set k 0 = 0 and n = . The results of this eigenvalue simulation for the frequency range 1c/aω ≤ 5.5c/a are plotted in Fig. 1 as ω vs. kx = · k = λ.

For clarity, we have only plotted the three least evanescent (i.e. possessing the smallest values of Im(kz)) eigenmodes. These three k(ω) dispersion curves are plotted as solid lines. For comparison, we have also ploted a conventional ω(k) dispersion curve as dotted line. The three eigenmodes in Fig. 1 are described as either transverse or longitudinal according to their polarization. The symmetry of the dispersion curves is such that for every solution k(ω) there is also the solution –k(ω) indicating that this is a reciprocal crystal. The dispersion curves for the transverse modes plotted in Fig. 1 in fact represent two polarization-degenerate modes because of the symmetry of the crystal. The longitudinal mode with the passband near ω = 4.5c/a is magnetically polarized in the direction making it a magnetic bulk plasmon. The longitudinal mode with the passband near ω = 5c/a is electrically polarized in the direction and is an electric bulk plasmon. The field profiles of both longitudinal modes indicate that the passbands correspond to Mie’s resonances of the dielectric sphere [31].

The transverse mode dispersion curve has a band in the approximate frequency range 4.6c/a < ω < 4.8c/a with a large value of Im(kx), indicating it is an evanescent band, but a Re(kx) that is equal to neither 0 nor π/a as is typical of ω(k) dispersion curves. As described in Refs. [13, 26] for a quadratic eigenvalue problem with hermitian matrices (corresponding to a lossless crystal) the eigenvalues must always be real or come in complex-conjugate pairs. The dielectric photonic crystal under consideration has very low loss, so this lossless condition approximately holds true for the dispersion curves in Fig. 1. The transverse band in the 4.6c/a < ω < 4.8c/a frequency band is one half of a complex conjugate pair, the other half is a transverse doubly degenerate mode not shown here. At the frequency of ω ≈ 4.8c/a the two modes that make up this complex conjugate pair both enter a passband and split, the plotted mode going to the Γ point and the unplotted mode going to the band edge (this unplotted mode corresponds to the dotted lines in Fig. 1 from the ω(k) simulation). Note that in this passband there are two pairs of propagating doubly polarization degenerate modes or four propagating modes in total.

The transverse eigenmodes plotted in Fig. 1 can be excited by a plane wave normally incident onto the vacuum/photonic crystal interface if the interface is parallel to the y-z plane. The longitudinally polarized modes could not be excited by a normally incident wave without the aid of a coupling device at the interface. If the incident beam of light is not normal to the interface, if for example the incident beam has a wavenumber laying in the x-y plane but at a 30° angle from normal then a different set of eigenmodes will be excited at the interface. To simulate these excited eigenmodes we set k 0 = ω/csin(π/6)ŷ and n = and solve the resulting eigenvalue problem. The resulting photonic band structure is plotted in Fig. 2.

 figure: Fig. 2

Fig. 2 Complex wavenumber dispersion curves and field profiles for eigenmodes of the photonic crystal pictured in Fig. 1(a) assuming k 0 = ω/csin(π/6)ŷ and n = . Modes excited by p or s polarized incident light are plotted with solid or dashed lines respectively. (a) Real part of kx(ω) for two transverse hybrid modes and a expanded view of the avoided crossing in Re(kx) space. (b) Imaginary part of kx(ω) for two transverse hybrid modes, an expanded view of the avoided crossing in Im(kx) space (plotting the same modes as the expanded view in Re(kx) space), and a field profile for the Ez polarized transverse hybrid mode. (c) Real part of kx(ω) for two longitudinal hybrid modes and a field profile for the magnetic longitudinal hybrid mode. (d) Imaginary part of kx(ω) for two longitudinal hybrid modes and a field profile for the electric longitudinal hybrid mode. The magnetic longitudinal hybrid mode is excited by s polarized incident light and the electric longitudinal hybrid mode is excited by p polarized incident light. For all field profiles the frequency is ω = 2c/a with arrows representing Dy and Dz and color representing Dx.

Download Full Size | PPT Slide | PDF

The eigenmodes in Fig. 2 are roughly split into predominantly-transverse (transverse hybrid) and predominantly-longitudinal (longitudinal hybrid) modes. The hybridization between the transverse and longitudinal modes is caused by the finite symmetry-breaking ky. Both the transverse and longitudinal hybrid modes in Fig. 2 can be characterized by the polarization of the incident light that couples to them. For a p polarized incident beam (electric field in the xy plane) the p polarized eigenmode is excited (plotted in Fig. 2 with solid lines) and for an s polarized incident beam (electric field in the direction) the s polarized eigenmode is excited (plotted in Fig. 2 with dashed lines).

At the frequency ω ≈ 4c/a the transverse hybrid modes and the longitudinal hybrids modes appear to cross in a propagating band. An expanded view of this region in Fig. 2 plotting both transverse and longitudinal hybrid modes shows that the apparent crossing actually occurs in a band gap. Viewed in complex kx space it is clear that this is actually an avoided crossing and in the band gap the transverse and longitudinal hybrid modes form complex conjugate pairs [13, 26].

We see that even for a simple photonic crystal the CWES method of calculating the dispersion curve produces rich and complex results. In particular, it is not possible to solve for evanescent eigenmodes using the conventional ω(k) method for calculating dispersion curves.

4. Example: negative index plasmonic fishnet

The second example highlighting the versatility of the CWES method is a plasmonic negative index metamaterial (NIM) shown in Fig. 3. Because this metamaterial contains dispersive (plasmonic) components, it is even more convenient to use this method because the dielectric permittivities of metals are tabulated for real-valued frequencies. Recent experiments [22] demonstrated that the so-called fishnet metamaterial supports a negative index eigenmode for near infrared wavelengths of about λ 0 ≈ 1.7μ m. The dimensions and composition of the unit cell taken from Ref. [22] are shown in Fig. 3. The fishnet metamaterial is made of alternating layers of Ag and MgF2 with thicknesses of 30nm and 50nm respectively. This layered structure was milled with a focused ion beam into crisscrossing strips with widths of 265nm and 565nm. The crystal lattice constants are ax = ay = 860nm and az = 80nm. For the dielectric response of the Ag we used the Drude model ɛAg=1ωp2/(ω(ωiγ)) with the same parameters cited in Ref. [22]: ωp = 9eV and γ = 0.054eV. It should be noted that in Ref. [22] the scattering frequency γ was increased by a factor of three from that of standard bulk Ag to account for additional scattering losses. For the dielectric response of MgF2, no values were provided in the original paper so we used an oscillator model of the dielectric permittivity from Ref. [32].

 figure: Fig. 3

Fig. 3 The fishnet metamaterial from Ref. [22]. The lattice constants of the unit cell are ax = ay = 860nm and az = 80nm. The fishnet is made up of alternating layers of Ag with a thickness of 30nm and MgF2 with a thickness of 50nm. The widths of the the crisscrossing fishnet strips are b 1 = 265nm and b 2 = 565nm.

Download Full Size | PPT Slide | PDF

The negative index mode reported in Ref. [22] propagates in the direction. Therefore, we have applied CWES to the specific case of k 0 = 0 and n = . The resulting dispersion curves are plotted in Fig. 4 for the 100THz < ω/2π < 300THz frequency range. For clarity we have only plotted the four eigenmodes that have the smallest values of Im(kz) and therefore the longest propagation lengths. The negative index mode found in Ref. [22] is labeled Ex to indicate that it is transversely polarized with the electric field pointing in the direction. We can see that it is a negative index mode because (according to the convention where fields are Bloch periodic with the exponential factor exp[i(ωtk · x)]) a negative index mode should have a wavenumber whose real and imaginary parts have the same sign. In Fig. 4, that describes the Ex mode for frequencies below 200THz. In addition, we observe three other negative index modes, though they generally have a shorter propagation lengths than the Ex mode. Of these three modes, the two labeled Ey are transversely polarized with the electric field pointing in the ŷ direction and the mode labeled Ez is longitudinally polarized with the electric field pointing in the direction.

 figure: Fig. 4

Fig. 4 TOP: Real (a) and imaginary (b) parts of kz(ω) for several eigenmodes of the fishnet metamaterial shown in Fig. 3. In addition to the transverse mode electrically polarized in the direction labeled Ex which was identified in Ref. [22] we see two other transverse modes electrically polarized in the ŷ direction labeled Ey and a longitudinal mode electrically polarized in the direction labeled Ez. MIDDLE: Field profiles for the transverse Ex mode (c) and the longitudinal Ez mode (d) on a cross-section laying on the x-y plane in the middle of the MgF2 layer. Arrows represent in-plane electric field and color represents the Ez field. BOTTOM: (e) Field profile of the transverse Ex mode on a cross-section laying on the x-z plane halfway between two thin Ag-MgF 2 strips. Arrows represent in-plane electric field and color represents the Hy field. For all field profiles the frequency is 175THz and each region is labeled Ag or MgF 2 according the the material of the region. Unlabeled regions are vacuum.

Download Full Size | PPT Slide | PDF

The figure of merit (FOM) is a number commonly used to quantify the quality of a negative index material. It is often defined [6] as the ratio between the real and imaginary parts of the index of refraction. For eigenmodes of the fishnet crystal with real ω and a complex k pointing in the direction this is equivalent to the ratio between the real and imaginary parts of kz or FOM≡ Re(kz)/Im(kz). Without knowledge of the complex wavenumber of an eigenmode of a metamaterial, the figure of merit must be calculated indirectly. For example, in Ref. [22] the FOM is estimated numerically from simulated transmission through a fishnet sample. With knowledge of the complex wavenumbers of the crystal eigenmodes, it is possible to calculate the FOM directly. The FOM for each of the four eigenmodes identified in Fig. 4 is plotted in Fig. 5(a), along with the propagation lengths plotted in Fig. 5(b). A positive FOM corresponds to the real and imaginary parts of kz having the same sign, therefore indicating a negative index eigenmode. We see again that all four modes are in fact negative index modes in a fairly broad 150THz < ω/2π < 200THz frequency range. Somewhat unexpectedly, the mode experimentally observed in Ref. [22] (labelled Ex) has the lowest FOM despite having the longest propagation length. The theoretically predicted FOM≈ 7.5 of the Ex mode is a factor of 2 higher than the experimentally measured FOM [22], which could be attributed to fabrication imperfections.

 figure: Fig. 5

Fig. 5 The FOM (a) and propagation length (b) for the four fishnet eigenmodes identified in Fig. 4. Note that though the mode labeled Ex has the largest propagation length it also has the smallest FOM. The negative value of FOM for the Ex mode for frequencies above 220THz indicates that it is no longer a negative index mode at these higher frequencies.

Download Full Size | PPT Slide | PDF

The existence of multiple negative index waves in the fishnet structure has important experimental implications. For some frequencies (e.g., at 150THz) the propagation length of the longitudinal mode (Ez) is as long as that of the primary transverse mode (Ex). While the original experiments [22] only excited the Ex mode by using the light normally incident onto the x–y vacuum/fishnet interface, one can envision exciting both Ex and Ez modes using obliquely incident light. For example, if p-polarized light is obliquely incident with the incident wavevector laying in the x-z plane (electric field laying in the x-z plane), then both Ex and Ez modes will be launched into the fishnet with different phase velocities. Their refraction by the fishnet prism [22] would give rise to two distinct beams. The existence of the additional strongly dispersive longitudinal modes (bulk plasmons) also suggests that the fishnet is a metamaterial with significant spatial dispersion that can have a strong effect on the surface waves at the vacuum/fishnet interface [33, 34]. Strong spatial dispersion is not unexpected because the fishnet’s transverse period is about λ/2.

Recent theoretical [35] and experimental [36] studies demonstrated that negative index propagation in fishnet structures is not limited to the optical frequency range and can, in fact, be observed with microwaves. The main physical distinction (excluding their appropriately scaled sizes) between microwave and optical structures is that the metals in the former can be accurately described as PECs. On the other hand, in the optical range metals are described as being ”plasmonic”, which means that optical field penetration into the metal is significant. Our simulation can quantify how important the plasmonic properties of the metal are for the fishnet structure from Ref. [22]. This is done by introducing the plasmonic parameter Tp [37], which is the number that characterizes the plasmonic nature of a metamaterial and it is defined as the ratio of the kinetic energy of the plasmonic electrons to the magnetic energy in the unit cell of a crystal. Strong plasmonic effects and the importance of electrostatic resonances imply Tp ≫ 1. We find that the plasmonic parameter for the mode labelled Ex in Fig. 4 at a frequency of 175THz is Tp = 0.24. Being less than one indicates this mode is not predominately plasmonic in nature because the electromagnetic fields do not significantly penetrate into the Ag. This is consistent with a recent demonstration of a negative index mode in a PEC fishnet structure [35, 36].

Finally, we describe another application of the CWES method: calculation of isofrequency ω(k) = const contours in metamaterial crystals. From the early days of metamaterials reseach, isofrequency diagrams have been used as a simple visual tool for studying refraction at the vacuum/metamaterial interfaces, especially in the context of negative index propagation and negative refraction [38]. Traditionally, isofrequency diagrams are drawn using a conventional ω(k) eigenvalue simulation [2]. This provides no information about propagation lengths which are important in lossy metamaterials such as the fishnet. The conventional approach is also highly laborious for plasmonic fishnets because of the strong frequency dependence of the dielectric permittivities of metals. Other semi-analytic techniques used for analyzing wave propagation in highly symmetric PEC-based fishnets [35] do not apply to plasmonic fishnets of interest.

Figure 6 shows an isofrequency diagram calculated using the CWES approach. A single isofrequency contour was obtain by fixing the real frequency ω, setting k 0 = ky ŷ and n = , and then scanning ky from –π/ay to π/ay. The eigenvalue in this simulation is a complex valued kz. This procedure was repeated for several values of ω from the 150THz to 180THz. The resulting eigenmodes can be be excited by a plane wave incident on an interface between vacuum and the fishnet crystal parallel with the -ŷ plane if the wavevector of the incident wave is confined to the ŷ- plane. Because the ŷ component of the incident wavevector is real valued, the ŷ component of the wave excited in the fishnet crystal must also be real-valued. ky = sinθω/c where θ is the incidence angle with respect to the normal z. The eigenmode excited in the fishnet crystal decays in the direction as indicated by the imaginary part of kz plotted in Fig. 6.

 figure: Fig. 6

Fig. 6 (a) Isofrequency contours for vacuum plotted with respect to kz and ky (b) Isofrequency contours for the Ex polarized eigenmode of the fishnet crystal plotted with respect to Re(kz) and ky. (c) Isofrequency contours for the Ex polarized eigenmode of the fishnet crystal plotted with respect to Im(kz) and ky. The black arrows indicate the direction of the phase velocity and the red arrows indicate the direction of the group velocity at a frequency of 170THz and ky = ω/csin(π/6).

Download Full Size | PPT Slide | PDF

The isofrequency contours in Fig. 6 are hyperbolic in appearance. We can study refraction at the interface between vacuum and the fishnet crystal by comparing the isofrequency contours of the fishnet to those of vacuum, which are also shown in Fig. 6. The vacuum isofrequency contours are circular. As can be seen in Fig. 6 for frequency of 170THz and with an incident angle of 30° the component of the incident wavevector tangential to the interface (ky) must be matched to the eigenmodes of the fishnet. For ky = ω/csin(π/6) and a frequency of 170THz the isofrequency diagram shows two eigenmodes with the correct value of ky and equal and opposite values of kz. We find the correct mode by calculating the group velocity vg = ∂ω/ Re(k) which is by definition normal to the isofrequency contours. In the absence of anomalous dispersion the group velocity indicates the direction of energy flow in the crystal [39]. Choosing the correct solution requires us to choose the solution that has energy flowing in the positive direction (i.e. away from the interface). This selects the solution with a negative Re(kz). This is a negative index mode in the sense that in the direction the phase velocity and group velocity have opposite signs.

However, because the shape of the isofrequency contours is hyperbolic, the phase and group velocities in the ŷ direction have the same sign. Therefore, positive refraction at the x-ŷ vacuum/fishnet interface is expected according to Fig. 6 despite the fact that a negative index eigen-mode is excited. This does not contradict the experiment by Valentine et al. in Ref. [22], where the fishnet structure was cut in the shape of a prism. In that experiment, the wave propagation through the fishnet was entirely in the direction enabling the group and phase velocities to be antiparallel. Negative refraction indeed occurs at a vacuum-fishnet interface tilted with respect to the principal axis of the crystal.

Finally, though we calculated the isofrequency contours in Fig. 6 by varying ky and solving for kz, we could have alternately varied the angle of k by setting k = λcos(θ) +λsin(θ)ŷ, varying the angle θ and solving for λ as the eigenvalue. This would have produced similar but slightly different isofrequency contours. This is an example of how the complex k dispersion curves in this paper are not unique and depend on how one restricts the three degrees of freedom kx, ky and kz to one degree of freedom λ.

5. Conclusion

In conclusion, we have presented a three-dimensional realization of the Complex Wavenumber Eigenvalue Simulation (CWES) approach to calculating photonic band structure of metamaterial/photonic crystals. A detailed implementation of CWES using FEM discretization is described. The CWES approach was applied to two periodic photonic structures: (a) photonic crystal comprised of dielectric spheres, and (b) plasmonic fishnet metamaterial supporting negative refractive index waves. For case (a) we have used the results of CWES to identify both transverse and longitudinal modes and investigated their coupling that gives rise to avoided crossings and bandgap formation. For case (b) we have identified, for the first time, four negative-index modes of the fishnet structure and computed hyperbolic isofrequency surfaces for the least-damped transverse mode. The sphere photonic crystal model is available for [21].

Acknowledgments

We would like to acknowledge the developers of MUMPS [40], ARPACK [41,42], PETSC [43] and SLEPC [44]. Support of the staff of the Texas Advanced Computing Center is gratefully acknowledged. This work was supported by the Air Force Research Laboratory and the Office of Naval Research.

References and links

2. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

3. K. Sakoda, Optical Properties of Photonic Crystals, 2nd ed. (Springer, 2004).

4. G. Shvets and Y. A. Urzhumov, “Engineering the electromagnetic properties of periodic nanostructures using electrostatic resonances,” Phys. Rev. Lett. 93, 243902 (2004). [CrossRef]  

5. E. D. Palik, Handbook of Optical Constants of Solids (Acedemic Press, 1985), vol. 1.

6. S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005). [CrossRef]  

7. 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 (2000). [CrossRef]  

8. R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004). [CrossRef]  

9. T. Suzuki and P. L. Yu, “Tunneling in photonic band strucures,” J. Opt. Soc. Am. B 12, 804 (1995). [CrossRef]  

10. E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).

11. A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995). [CrossRef]  

12. Z. Y. Li and L. L. Lin, “Photonic band structure solved by a plane-wave–based transfer-matrix method,” Phys. Rev. E 67, 046607 (2003). [CrossRef]  

13. M. Davanco, Y. Urzhumov, and G. Shvets, “The complex bloch bands of a 2d plasmonic crystal displaying isotropic negative refraction,” Opt. Express 15, 9681–9691 (2007). [CrossRef]   [PubMed]  

14. C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009). [CrossRef]  

15. X. Zhang, M. Davanco, Y. Urzhumov, and G. Shvets, “A subwavelength near-infrared negative index material,” Appl. Phys. Lett. 94, 131107 (2009). [CrossRef]  

16. X. Zhang, M. Davanco, K. Miller, T. W. J. C. Wu, C. Fietz, D. Korobkin, X. Li, G. Shvets, and S. R. Forrest, “Interferometric characterization of a subwavelength near-infrared negative index metamaterial,” Opt. Express 18, 17788–17795 (2010). [CrossRef]   [PubMed]  

17. C. Fietz and G. Shvets, “Current-driven metamaterial homogenization,” Phys. B 405, 2930–2934 (2010). [CrossRef]  

18. M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27, 1576–1582 (2010). [CrossRef]  

19. C. Fietz and G. Shvets, “Homogenization theory for simple metamaterials modeled as one-dimensional arrays of thin polarizable sheets,” Phys. Rev. B 82, 205128 (2010). [CrossRef]  

20. A. Pors, I. Tsukerman, and S. I. Bozhevolnyi, “Effective constitutive parameters of plasmonic metamaterials: a rigorous approach,” arXiv:1104.2972v1(2011).

21. Both .m and .mph files for the sphere photonic crystal model can be acquired by contacting Chris Fietz at fietz.chris@gmail.com.

22. J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

23. W. B. J. Zimmerman, Process Modelling and Simulation with Finite Element Methods (World Scientific, 2004).

24. J. Jin, The Finite Element Method in Electromagnetics, 2nd ed. (John Wiley & Sons, Inc., 2002).

25. P. P. Silvester and R. L. Ferrari, Finite Elements for Electrical Engineers3rd ed. (Cambridge University Press, 1996).

26. F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review 43, 235–286 (2001). [CrossRef]  

27. C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003). [CrossRef]  

28. R. A. Shore and A. D. Yaghjian, “Traveling waves on two- and three-dimensional periodic arrays of lossless scatterers,” Radio Sci. 42, RS6S21 (2007). [CrossRef]  

29. C. R. Simovski, “Bloch material parameters of magneto-dielectric metamaterials and the concept of Bloch lattices,” Metamaterials 1, 62–80 (2007). [CrossRef]  

30. A. Alú, “First-principle homogenization theory for periodic metamaterial arrays,” arXiv:1012.1351v2(2011).

31. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Sons, 1998). [CrossRef]  

32. J. M. Siqueiros, R. Machorro, and L. E. Regalado, “Determination of the optical constants of MgF2 and ZnS from spectrophotometric measurements and the classical oscillator method,” Appl. Opt. 27, 2549–2553 (1988). [CrossRef]   [PubMed]  

33. P. Halevi, Spatial Dispersion in Solids and Plasmas (Elsevier Science Publishers, 1992).

34. M. A. Shapiro, G. Shvets, J. R. Sirigiri, and R. J. Temkin, “Spatial dispersion in metamaterials with negative dielectric permittivity and its effect on surface waves,” Opt. Lett. 31, 2051 (2006). [CrossRef]   [PubMed]  

35. L. Jelinek, R. Marques, and J. Machac, “Fishnet metamaterials - rules for refraction and limits of homogenization,” Opt. Express 18, 17940–17949 (2010). [CrossRef]   [PubMed]  

36. M. Navarro-Ca, M. Beruete, M. Sorolla, and I. Campillo, “Negative refraction in a prism made of stacked subwavelength hole arrays,” Opt. Express 16, 560–566 (2008). [CrossRef]  

37. Y. A. Urzhumov and G. Shvets, “Optical magnetism and negative refraction in plasmonic metamaterials,” Solid State Commun. 146, 208–220 (2008). [CrossRef]  

38. C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).

39. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, 1960).

40. MUltifrontal Massively Parallel Solver (MUMPS 4.9.2) User’s guide (Lyon, 2009).

41. R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (SIAM, Philadelphia). [PubMed]   [PubMed]  

42. K. J. Maschhoff and D. C. Soerensen, “PARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures,” Lect. Notes Comp. Sci. 1184, 478–486 (1996). [CrossRef]  

43. S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).

44. J. E. Roman, E. Romero, and A. Tomas, “SLEPc users manual,” Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010).

References

  • View by:

  1. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).
  2. K. Sakoda, Optical Properties of Photonic Crystals, 2nd ed. (Springer, 2004).
  3. G. Shvets and Y. A. Urzhumov, “Engineering the electromagnetic properties of periodic nanostructures using electrostatic resonances,” Phys. Rev. Lett. 93, 243902 (2004).
  4. E. D. Palik, Handbook of Optical Constants of Solids (Acedemic Press, 1985), vol. 1.
    [Crossref]
  5. S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).
  6. 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 (2000).
    [Crossref]
  7. R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
    [Crossref]
  8. T. Suzuki and P. L. Yu, “Tunneling in photonic band strucures,” J. Opt. Soc. Am. B 12, 804 (1995).
    [Crossref]
  9. E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).
    [Crossref]
  10. A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995).
  11. Z. Y. Li and L. L. Lin, “Photonic band structure solved by a plane-wave–based transfer-matrix method,” Phys. Rev. E 67, 046607 (2003).
    [Crossref]
  12. M. Davanco, Y. Urzhumov, and G. Shvets, “The complex bloch bands of a 2d plasmonic crystal displaying isotropic negative refraction,” Opt. Express 15, 9681–9691 (2007).
    [Crossref]
  13. C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009).
    [Crossref] [PubMed]
  14. X. Zhang, M. Davanco, Y. Urzhumov, and G. Shvets, “A subwavelength near-infrared negative index material,” Appl. Phys. Lett. 94, 131107 (2009).
    [Crossref]
  15. X. Zhang, M. Davanco, K. Miller, T. W. J. C. Wu, C. Fietz, D. Korobkin, X. Li, G. Shvets, and S. R. Forrest, “Interferometric characterization of a subwavelength near-infrared negative index metamaterial,” Opt. Express 18, 17788–17795 (2010).
    [Crossref]
  16. C. Fietz and G. Shvets, “Current-driven metamaterial homogenization,” Phys. B 405, 2930–2934 (2010).
    [Crossref] [PubMed]
  17. M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27, 1576–1582 (2010).
    [Crossref]
  18. C. Fietz and G. Shvets, “Homogenization theory for simple metamaterials modeled as one-dimensional arrays of thin polarizable sheets,” Phys. Rev. B 82, 205128 (2010).
    [Crossref]
  19. A. Pors, I. Tsukerman, and S. I. Bozhevolnyi, “Effective constitutive parameters of plasmonic metamaterials: a rigorous approach,” arXiv:1104.2972v1(2011).
    [Crossref]
  20. Both .m and .mph files for the sphere photonic crystal model can be acquired by contacting Chris Fietz at fietz.chris@gmail.com.
  21. J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).
  22. W. B. J. Zimmerman, Process Modelling and Simulation with Finite Element Methods (World Scientific, 2004).
  23. J. Jin, The Finite Element Method in Electromagnetics, 2nd ed. (John Wiley & Sons, Inc., 2002).
  24. P. P. Silvester and R. L. Ferrari, Finite Elements for Electrical Engineers3rd ed. (Cambridge University Press, 1996).
  25. F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review 43, 235–286 (2001).
  26. C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
    [Crossref]
  27. R. A. Shore and A. D. Yaghjian, “Traveling waves on two- and three-dimensional periodic arrays of lossless scatterers,” Radio Sci. 42, RS6S21 (2007).
    [Crossref]
  28. C. R. Simovski, “Bloch material parameters of magneto-dielectric metamaterials and the concept of Bloch lattices,” Metamaterials 1, 62–80 (2007).
    [Crossref]
  29. A. Alú, “First-principle homogenization theory for periodic metamaterial arrays,” arXiv:1012.1351v2(2011).
    [Crossref]
  30. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Sons, 1998).
  31. J. M. Siqueiros, R. Machorro, and L. E. Regalado, “Determination of the optical constants of MgF2 and ZnS from spectrophotometric measurements and the classical oscillator method,” Appl. Opt. 27, 2549–2553 (1988).
    [Crossref]
  32. P. Halevi, Spatial Dispersion in Solids and Plasmas (Elsevier Science Publishers, 1992).
    [Crossref] [PubMed]
  33. M. A. Shapiro, G. Shvets, J. R. Sirigiri, and R. J. Temkin, “Spatial dispersion in metamaterials with negative dielectric permittivity and its effect on surface waves,” Opt. Lett. 31, 2051 (2006).
  34. L. Jelinek, R. Marques, and J. Machac, “Fishnet metamaterials - rules for refraction and limits of homogenization,” Opt. Express 18, 17940–17949 (2010).
    [Crossref] [PubMed]
  35. M. Navarro-Ca, M. Beruete, M. Sorolla, and I. Campillo, “Negative refraction in a prism made of stacked subwavelength hole arrays,” Opt. Express 16, 560–566 (2008).
    [Crossref] [PubMed]
  36. Y. A. Urzhumov and G. Shvets, “Optical magnetism and negative refraction in plasmonic metamaterials,” Solid State Commun. 146, 208–220 (2008).
    [Crossref]
  37. C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
    [Crossref]
  38. L. Brillouin, Wave Propagation and Group Velocity (Academic Press, 1960).
  39. MUltifrontal Massively Parallel Solver (MUMPS 4.9.2) User’s guide (Lyon, 2009).
  40. R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (SIAM, Philadelphia).
  41. K. J. Maschhoff and D. C. Soerensen, “PARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures,” Lect. Notes Comp. Sci. 1184, 478–486 (1996).
    [PubMed]
  42. S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
    [Crossref]
  43. J. E. Roman, E. Romero, and A. Tomas, “SLEPc users manual,” Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010).

2011 (1)

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
[Crossref]

2010 (5)

2009 (2)

C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009).
[Crossref] [PubMed]

X. Zhang, M. Davanco, Y. Urzhumov, and G. Shvets, “A subwavelength near-infrared negative index material,” Appl. Phys. Lett. 94, 131107 (2009).
[Crossref]

2008 (3)

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

M. Navarro-Ca, M. Beruete, M. Sorolla, and I. Campillo, “Negative refraction in a prism made of stacked subwavelength hole arrays,” Opt. Express 16, 560–566 (2008).
[Crossref] [PubMed]

Y. A. Urzhumov and G. Shvets, “Optical magnetism and negative refraction in plasmonic metamaterials,” Solid State Commun. 146, 208–220 (2008).
[Crossref]

2007 (3)

R. A. Shore and A. D. Yaghjian, “Traveling waves on two- and three-dimensional periodic arrays of lossless scatterers,” Radio Sci. 42, RS6S21 (2007).
[Crossref]

C. R. Simovski, “Bloch material parameters of magneto-dielectric metamaterials and the concept of Bloch lattices,” Metamaterials 1, 62–80 (2007).
[Crossref]

M. Davanco, Y. Urzhumov, and G. Shvets, “The complex bloch bands of a 2d plasmonic crystal displaying isotropic negative refraction,” Opt. Express 15, 9681–9691 (2007).
[Crossref]

2006 (1)

2005 (1)

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

2004 (2)

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

G. Shvets and Y. A. Urzhumov, “Engineering the electromagnetic properties of periodic nanostructures using electrostatic resonances,” Phys. Rev. Lett. 93, 243902 (2004).

2003 (2)

Z. Y. Li and L. L. Lin, “Photonic band structure solved by a plane-wave–based transfer-matrix method,” Phys. Rev. E 67, 046607 (2003).
[Crossref]

C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
[Crossref]

2001 (1)

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review 43, 235–286 (2001).

2000 (1)

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 (2000).
[Crossref]

1996 (1)

K. J. Maschhoff and D. C. Soerensen, “PARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures,” Lect. Notes Comp. Sci. 1184, 478–486 (1996).
[PubMed]

1995 (2)

T. Suzuki and P. L. Yu, “Tunneling in photonic band strucures,” J. Opt. Soc. Am. B 12, 804 (1995).
[Crossref]

A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995).

1988 (1)

1951 (1)

E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).
[Crossref]

Alú, A.

A. Alú, “First-principle homogenization theory for periodic metamaterial arrays,” arXiv:1012.1351v2(2011).
[Crossref]

Asatryan, A. A.

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

Baker-Jarvis, J.

C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
[Crossref]

Balay, S.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Bartal, G.

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

Beruete, M.

Bohren, C. F.

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Sons, 1998).

Botten, L. C.

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

Bozhevolnyi, S. I.

A. Pors, I. Tsukerman, and S. I. Bozhevolnyi, “Effective constitutive parameters of plasmonic metamaterials: a rigorous approach,” arXiv:1104.2972v1(2011).
[Crossref]

Brillouin, L.

L. Brillouin, Wave Propagation and Group Velocity (Academic Press, 1960).

Brown, J.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Brueck, S. R. J.

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

Buschelman, K.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Campillo, I.

Conforti, M.

Davanco, M.

de Sterke, C. M.

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

Eijkhout, V.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Engström, C.

C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009).
[Crossref] [PubMed]

Fan, W.

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

Ferrari, R. L.

P. P. Silvester and R. L. Ferrari, Finite Elements for Electrical Engineers3rd ed. (Cambridge University Press, 1996).

Fietz, C.

C. Fietz and G. Shvets, “Homogenization theory for simple metamaterials modeled as one-dimensional arrays of thin polarizable sheets,” Phys. Rev. B 82, 205128 (2010).
[Crossref]

X. Zhang, M. Davanco, K. Miller, T. W. J. C. Wu, C. Fietz, D. Korobkin, X. Li, G. Shvets, and S. R. Forrest, “Interferometric characterization of a subwavelength near-infrared negative index metamaterial,” Opt. Express 18, 17788–17795 (2010).
[Crossref]

C. Fietz and G. Shvets, “Current-driven metamaterial homogenization,” Phys. B 405, 2930–2934 (2010).
[Crossref] [PubMed]

Forrest, S. R.

Genov, D. A.

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

Green, A. A.

E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).
[Crossref]

Gropp, W. D.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Guasoni, M.

Hafner, C.

C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009).
[Crossref] [PubMed]

Halevi, P.

P. Halevi, Spatial Dispersion in Solids and Plasmas (Elsevier Science Publishers, 1992).
[Crossref] [PubMed]

Holloway, C. L.

C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
[Crossref]

Huffman, D. R.

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Sons, 1998).

Istrate, E.

E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).
[Crossref]

Jelinek, L.

Jin, J.

J. Jin, The Finite Element Method in Electromagnetics, 2nd ed. (John Wiley & Sons, Inc., 2002).

Joannopoulos, J. D.

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
[Crossref]

J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

Johnson, S. G.

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
[Crossref]

Kabos, P.

C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
[Crossref]

Kaushik, D.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Knepley, M. G.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Korobkin, D.

Kuester, E. F.

C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
[Crossref]

Lehoucq, R. B.

R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (SIAM, Philadelphia).

Li, X.

Li, Z. Y.

Z. Y. Li and L. L. Lin, “Photonic band structure solved by a plane-wave–based transfer-matrix method,” Phys. Rev. E 67, 046607 (2003).
[Crossref]

Lin, L. L.

Z. Y. Li and L. L. Lin, “Photonic band structure solved by a plane-wave–based transfer-matrix method,” Phys. Rev. E 67, 046607 (2003).
[Crossref]

Luo, C.

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
[Crossref]

Machac, J.

Machorro, R.

Malloy, K. J.

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

Marques, R.

Maschhoff, K. J.

K. J. Maschhoff and D. C. Soerensen, “PARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures,” Lect. Notes Comp. Sci. 1184, 478–486 (1996).
[PubMed]

McInnes, L. C.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

McOrist, J.

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

McPhedran, R. C.

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

Meade, R. D.

J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

Meerbergen, K.

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review 43, 235–286 (2001).

Miller, K.

Navarro-Ca, M.

Nicorovici, N. A.

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

Notomi, M.

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 (2000).
[Crossref]

Osgood, R. M.

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

Palik, E. D.

E. D. Palik, Handbook of Optical Constants of Solids (Acedemic Press, 1985), vol. 1.
[Crossref]

Panoiu, N. C.

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

Pendry, J. B.

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
[Crossref]

A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995).

Pors, A.

A. Pors, I. Tsukerman, and S. I. Bozhevolnyi, “Effective constitutive parameters of plasmonic metamaterials: a rigorous approach,” arXiv:1104.2972v1(2011).
[Crossref]

Regalado, L. E.

Roman, J. E.

J. E. Roman, E. Romero, and A. Tomas, “SLEPc users manual,” Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010).

Romero, E.

J. E. Roman, E. Romero, and A. Tomas, “SLEPc users manual,” Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010).

Sakoda, K.

K. Sakoda, Optical Properties of Photonic Crystals, 2nd ed. (Springer, 2004).

Sargent, E. H.

E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).
[Crossref]

Schmidt, K.

C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009).
[Crossref] [PubMed]

Shapiro, M. A.

Shore, R. A.

R. A. Shore and A. D. Yaghjian, “Traveling waves on two- and three-dimensional periodic arrays of lossless scatterers,” Radio Sci. 42, RS6S21 (2007).
[Crossref]

Shvets, G.

C. Fietz and G. Shvets, “Homogenization theory for simple metamaterials modeled as one-dimensional arrays of thin polarizable sheets,” Phys. Rev. B 82, 205128 (2010).
[Crossref]

C. Fietz and G. Shvets, “Current-driven metamaterial homogenization,” Phys. B 405, 2930–2934 (2010).
[Crossref] [PubMed]

X. Zhang, M. Davanco, K. Miller, T. W. J. C. Wu, C. Fietz, D. Korobkin, X. Li, G. Shvets, and S. R. Forrest, “Interferometric characterization of a subwavelength near-infrared negative index metamaterial,” Opt. Express 18, 17788–17795 (2010).
[Crossref]

X. Zhang, M. Davanco, Y. Urzhumov, and G. Shvets, “A subwavelength near-infrared negative index material,” Appl. Phys. Lett. 94, 131107 (2009).
[Crossref]

Y. A. Urzhumov and G. Shvets, “Optical magnetism and negative refraction in plasmonic metamaterials,” Solid State Commun. 146, 208–220 (2008).
[Crossref]

M. Davanco, Y. Urzhumov, and G. Shvets, “The complex bloch bands of a 2d plasmonic crystal displaying isotropic negative refraction,” Opt. Express 15, 9681–9691 (2007).
[Crossref]

M. A. Shapiro, G. Shvets, J. R. Sirigiri, and R. J. Temkin, “Spatial dispersion in metamaterials with negative dielectric permittivity and its effect on surface waves,” Opt. Lett. 31, 2051 (2006).

G. Shvets and Y. A. Urzhumov, “Engineering the electromagnetic properties of periodic nanostructures using electrostatic resonances,” Phys. Rev. Lett. 93, 243902 (2004).

Silvester, P. P.

P. P. Silvester and R. L. Ferrari, Finite Elements for Electrical Engineers3rd ed. (Cambridge University Press, 1996).

Simovski, C. R.

C. R. Simovski, “Bloch material parameters of magneto-dielectric metamaterials and the concept of Bloch lattices,” Metamaterials 1, 62–80 (2007).
[Crossref]

Siqueiros, J. M.

Sirigiri, J. R.

Smith, B. F.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Soerensen, D. C.

K. J. Maschhoff and D. C. Soerensen, “PARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures,” Lect. Notes Comp. Sci. 1184, 478–486 (1996).
[PubMed]

Sorensen, D. C.

R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (SIAM, Philadelphia).

Sorolla, M.

Stewart, W. J.

A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995).

Suzuki, T.

Temkin, R. J.

Tisseur, F.

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review 43, 235–286 (2001).

Tomas, A.

J. E. Roman, E. Romero, and A. Tomas, “SLEPc users manual,” Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010).

Tsukerman, I.

A. Pors, I. Tsukerman, and S. I. Bozhevolnyi, “Effective constitutive parameters of plasmonic metamaterials: a rigorous approach,” arXiv:1104.2972v1(2011).
[Crossref]

Ulin-Avila, E.

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

Urzhumov, Y.

X. Zhang, M. Davanco, Y. Urzhumov, and G. Shvets, “A subwavelength near-infrared negative index material,” Appl. Phys. Lett. 94, 131107 (2009).
[Crossref]

M. Davanco, Y. Urzhumov, and G. Shvets, “The complex bloch bands of a 2d plasmonic crystal displaying isotropic negative refraction,” Opt. Express 15, 9681–9691 (2007).
[Crossref]

Urzhumov, Y. A.

Y. A. Urzhumov and G. Shvets, “Optical magnetism and negative refraction in plasmonic metamaterials,” Solid State Commun. 146, 208–220 (2008).
[Crossref]

G. Shvets and Y. A. Urzhumov, “Engineering the electromagnetic properties of periodic nanostructures using electrostatic resonances,” Phys. Rev. Lett. 93, 243902 (2004).

Valentine, J.

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

Ward, A. J.

A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995).

Winn, J. N.

J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

Wu, T. W. J. C.

Yaghjian, A. D.

R. A. Shore and A. D. Yaghjian, “Traveling waves on two- and three-dimensional periodic arrays of lossless scatterers,” Radio Sci. 42, RS6S21 (2007).
[Crossref]

Yang, C.

R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (SIAM, Philadelphia).

Yu, P. L.

Zentgraf, T.

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

Zhang, H.

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

Zhang, S.

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. Panoiu, and R. M. Osgood, “Near-infrared double negative metamaterials,” Opt. Express 613, 4922 (2005).

Zhang, X.

Zimmerman, W. B. J.

W. B. J. Zimmerman, Process Modelling and Simulation with Finite Element Methods (World Scientific, 2004).

Appl. Opt. (1)

Appl. Phys. Lett. (1)

X. Zhang, M. Davanco, Y. Urzhumov, and G. Shvets, “A subwavelength near-infrared negative index material,” Appl. Phys. Lett. 94, 131107 (2009).
[Crossref]

IEEE Trans. Antennas Propag. (1)

C. L. Holloway, E. F. Kuester, J. Baker-Jarvis, and P. Kabos, “A double negative (DNG) composite medium composed of magnetodielectric spherical particles embedded in a matrix,” IEEE Trans. Antennas Propag. 51, 2596–2603 (2003).
[Crossref]

J. Comput. Theor. Nanosci. (1)

C. Engström, C. Hafner, and K. Schmidt, “Computations of lossy bloch waves in two-dimensional photonic crystals,” J. Comput. Theor. Nanosci. 6, 1–9 (2009).
[Crossref] [PubMed]

J. Opt. Soc. Am. B (2)

J. Phys.: Condens. Matter (1)

A. J. Ward, J. B. Pendry, and W. J. Stewart, “Photonic dispersion surfaces,” J. Phys.: Condens. Matter 7, 2217–2224 (1995).

Lect. Notes Comp. Sci. (1)

K. J. Maschhoff and D. C. Soerensen, “PARPACK: An efficient portable large scale eigenvalue package for distributed memory parallel architectures,” Lect. Notes Comp. Sci. 1184, 478–486 (1996).
[PubMed]

Metamaterials (1)

C. R. Simovski, “Bloch material parameters of magneto-dielectric metamaterials and the concept of Bloch lattices,” Metamaterials 1, 62–80 (2007).
[Crossref]

Opt. Express (5)

Opt. Lett. (1)

Phys. B (1)

C. Fietz and G. Shvets, “Current-driven metamaterial homogenization,” Phys. B 405, 2930–2934 (2010).
[Crossref] [PubMed]

Phys. Rev. B (4)

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 (2000).
[Crossref]

E. Istrate, A. A. Green, and E. H. Sargent, “Behavior of light at photonic crystal interfaces,” Phys. Rev. B 71, 195122 (2005).
[Crossref]

C. Fietz and G. Shvets, “Homogenization theory for simple metamaterials modeled as one-dimensional arrays of thin polarizable sheets,” Phys. Rev. B 82, 205128 (2010).
[Crossref]

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B 65, 201104 (2002).
[Crossref]

Phys. Rev. E (2)

Z. Y. Li and L. L. Lin, “Photonic band structure solved by a plane-wave–based transfer-matrix method,” Phys. Rev. E 67, 046607 (2003).
[Crossref]

R. C. McPhedran, L. C. Botten, J. McOrist, A. A. Asatryan, C. M. de Sterke, and N. A. Nicorovici, “Density of states functions for photonic crystals,” Phys. Rev. E 69, 016609 (2004).
[Crossref]

Phys. Rev. Lett. (1)

G. Shvets and Y. A. Urzhumov, “Engineering the electromagnetic properties of periodic nanostructures using electrostatic resonances,” Phys. Rev. Lett. 93, 243902 (2004).

Radio Sci. (1)

R. A. Shore and A. D. Yaghjian, “Traveling waves on two- and three-dimensional periodic arrays of lossless scatterers,” Radio Sci. 42, RS6S21 (2007).
[Crossref]

Science (1)

J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, and G. Bartal, “Three-dimensional optical metamaterial with a negative refractive index,” Science 455, 376–380 (2008).

SIAM Review (1)

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review 43, 235–286 (2001).

Solid State Commun. (1)

Y. A. Urzhumov and G. Shvets, “Optical magnetism and negative refraction in plasmonic metamaterials,” Solid State Commun. 146, 208–220 (2008).
[Crossref]

Other (16)

A. Alú, “First-principle homogenization theory for periodic metamaterial arrays,” arXiv:1012.1351v2(2011).
[Crossref]

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Sons, 1998).

P. Halevi, Spatial Dispersion in Solids and Plasmas (Elsevier Science Publishers, 1992).
[Crossref] [PubMed]

W. B. J. Zimmerman, Process Modelling and Simulation with Finite Element Methods (World Scientific, 2004).

J. Jin, The Finite Element Method in Electromagnetics, 2nd ed. (John Wiley & Sons, Inc., 2002).

P. P. Silvester and R. L. Ferrari, Finite Elements for Electrical Engineers3rd ed. (Cambridge University Press, 1996).

E. D. Palik, Handbook of Optical Constants of Solids (Acedemic Press, 1985), vol. 1.
[Crossref]

J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University Press, 2008).

K. Sakoda, Optical Properties of Photonic Crystals, 2nd ed. (Springer, 2004).

A. Pors, I. Tsukerman, and S. I. Bozhevolnyi, “Effective constitutive parameters of plasmonic metamaterials: a rigorous approach,” arXiv:1104.2972v1(2011).
[Crossref]

Both .m and .mph files for the sphere photonic crystal model can be acquired by contacting Chris Fietz at fietz.chris@gmail.com.

L. Brillouin, Wave Propagation and Group Velocity (Academic Press, 1960).

MUltifrontal Massively Parallel Solver (MUMPS 4.9.2) User’s guide (Lyon, 2009).

R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (SIAM, Philadelphia).

S. Balay, J. Brown, K. Buschelman, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, B. F. Smith, and H. Zhang, “PETSc users manual,” Tech. Rep. ANL-95/11 - Revision 3.1, Argonne National Laboratory (2010).
[Crossref]

J. E. Roman, E. Romero, and A. Tomas, “SLEPc users manual,” Tech. Rep. DSIC-II/24/02 - Revision 3.1, D. Sistemas Informáticos y Computación, Universidad Politécnica de Valencia (2010).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 Complex k dispersion curves and field profiles for eigenmodes of the photonic crystal pictured in the inset assuming k 0 = 0 and n = . (a) Real part of kx (ω) for a transversely polarized mode and a diagram of the crystal unit cell. (b) Imaginary part of kx (ω) for a transversely polarized mode and a field profile for the polarized transverse mode. There are two transverse modes, ŷ and electrically polarized, which are degenerate. (c) Real part of kx (ω) for two longitudinally polarized modes and a field profile for the magnetic longitudinal mode. (d) Imaginary part of kx (ω) for two longitudinally polarized modes and a field profile for the electric longitudinal mode. The longitudinal mode with the passband near ω = 4.5c/a is magnetically polarized in the direction and the longitudinal mode with the passband near ω = 5c/a is electrically polarized in the direction. The longitudinal modes correspond to Mie’s dipole resonances. For all dispersion curves the dotted lines are the result of a conventional ω(k) eigenvalue simulation. For all field profiles the frequency is ω = 2c/a with arrows representing D y and D z and color representing D x .
Fig. 2
Fig. 2 Complex wavenumber dispersion curves and field profiles for eigenmodes of the photonic crystal pictured in Fig. 1(a) assuming k 0 = ω/csin(π/6)ŷ and n = . Modes excited by p or s polarized incident light are plotted with solid or dashed lines respectively. (a) Real part of kx (ω) for two transverse hybrid modes and a expanded view of the avoided crossing in Re(kx ) space. (b) Imaginary part of kx (ω) for two transverse hybrid modes, an expanded view of the avoided crossing in Im(kx ) space (plotting the same modes as the expanded view in Re(kx ) space), and a field profile for the E z polarized transverse hybrid mode. (c) Real part of kx (ω) for two longitudinal hybrid modes and a field profile for the magnetic longitudinal hybrid mode. (d) Imaginary part of kx (ω) for two longitudinal hybrid modes and a field profile for the electric longitudinal hybrid mode. The magnetic longitudinal hybrid mode is excited by s polarized incident light and the electric longitudinal hybrid mode is excited by p polarized incident light. For all field profiles the frequency is ω = 2c/a with arrows representing D y and D z and color representing D x .
Fig. 3
Fig. 3 The fishnet metamaterial from Ref. [22]. The lattice constants of the unit cell are ax = ay = 860nm and az = 80nm. The fishnet is made up of alternating layers of Ag with a thickness of 30nm and MgF2 with a thickness of 50nm. The widths of the the crisscrossing fishnet strips are b 1 = 265nm and b 2 = 565nm.
Fig. 4
Fig. 4 TOP: Real (a) and imaginary (b) parts of kz (ω) for several eigenmodes of the fishnet metamaterial shown in Fig. 3. In addition to the transverse mode electrically polarized in the direction labeled E x which was identified in Ref. [22] we see two other transverse modes electrically polarized in the ŷ direction labeled E y and a longitudinal mode electrically polarized in the direction labeled E z . MIDDLE: Field profiles for the transverse E x mode (c) and the longitudinal E z mode (d) on a cross-section laying on the x-y plane in the middle of the MgF2 layer. Arrows represent in-plane electric field and color represents the E z field. BOTTOM: (e) Field profile of the transverse E x mode on a cross-section laying on the x-z plane halfway between two thin Ag-MgF 2 strips. Arrows represent in-plane electric field and color represents the H y field. For all field profiles the frequency is 175THz and each region is labeled Ag or MgF 2 according the the material of the region. Unlabeled regions are vacuum.
Fig. 5
Fig. 5 The FOM (a) and propagation length (b) for the four fishnet eigenmodes identified in Fig. 4. Note that though the mode labeled Ex has the largest propagation length it also has the smallest FOM. The negative value of FOM for the Ex mode for frequencies above 220THz indicates that it is no longer a negative index mode at these higher frequencies.
Fig. 6
Fig. 6 (a) Isofrequency contours for vacuum plotted with respect to kz and ky (b) Isofrequency contours for the E x polarized eigenmode of the fishnet crystal plotted with respect to Re(kz ) and ky . (c) Isofrequency contours for the E x polarized eigenmode of the fishnet crystal plotted with respect to Im(kz ) and ky . The black arrows indicate the direction of the phase velocity and the red arrows indicate the direction of the group velocity at a frequency of 170THz and ky = ω/csin(π/6).

Equations (12)

Equations on this page are rendered with MathJax. Learn more.

× ( 1 ɛ × H ) μ ω 2 c 2 H = 0 .
H ( x ) = u ( x ) exp [ i ( ω t k x ) ] ,
k 2 ɛ u k ɛ ( k u ) i k × ( 1 ɛ × u ) i × ( 1 ɛ k × u ) + × ( 1 ɛ × u ) μ ω 2 c 2 u = 0 ,
E ( x ) = 1 i ɛ ω / c × H = 1 i ɛ ω / c ( i k × u + × u ) exp [ i ( ω t k x ) ]
F H ( v , u ) = k 2 ɛ v u 1 ɛ ( k v ) ( k u ) i 1 ɛ v [ k × ( × u ) ] i ( × v ) 1 ɛ ( k × u ) + ( × v ) 1 ɛ ( × u ) μ ω 2 c 2 v u ,
0 = Ω d 3 x F H ( v , u ) = Ω d 3 x v [ 1 ɛ k × ( k × u ) i 1 ɛ k × ( × u ) i × ( 1 ɛ k × u ) + × ( 1 ɛ × u ) μ ω 2 c 2 u ] + Ω d A v [ n ^ × 1 ɛ ( i k × u + × u ) ] ,
A u + λ B u + λ 2 C u = 0 ,
( A B 0 1 ) ( u λ u ) = λ ( 0 C 1 0 ) ( u λ u ) .
× ( 1 μ × E ) ɛ ω 2 c 2 E = 0 .
E ( x ) = u ( x ) exp [ i ( ω t k x ) ] ,
k 2 μ u k μ ( k u ) i k × ( 1 μ × u ) i × ( 1 μ k × u ) + × ( 1 μ × u ) ɛ ω 2 c 2 u = 0.
F E ( v ,  u ) = k 2 μ v u 1 μ ( k v ) ( k u ) i 1 μ v [ k × ( × u ) ] i ( × v ) 1 μ ( k × u ) + ( × v ) 1 μ ( × u ) ɛ ω 2 c 2 v u ,

Metrics

Select as filters


Select Topics Cancel
© Copyright 2022 | Optica Publishing Group. All Rights Reserved