## Abstract

As a consequent work of the two-dimensional (2D) dendritic metamaterial which has been systematically studied in our previous work, a model of three-dimensional (3D) sphere-rod shaped structure is presented. Based on Drude model of the dielectric function of silver in the visible region, the parametric curves of electromagnetic response to the incident fields have been retrieved from detailed simulations. It is shown that the simultaneously negative values of permittivity and permeability in the optical range lead to a negative refractive index (NIM) through adjusting structural parameters, only the dimensions of the unit cells satisfy the effective medium theory. We therefore conclude that the proposed model offers a feasible route to fabricating 3D optical NIMs by ‘bottom-up’ approach.

©2010 Optical Society of America

## 1. Introduction

Artificially engineered negative-index metamaterials (NIMs) of which compositions and morphology can be rationally designed and tailored play an increasingly important role in the physics, optics, and engineering communities, owing to the potential applications [1,2] resulting from their extraordinary properties [3–5]. Artificial magnetism and negative refractive index are two specific types of behavior in this man-made materials which have been shifted from gigahertz frequencies all the way to terahertz frequencies over the past decade [6–8], and now approach the visible part of the spectrum [9,10]. The real realization of negative refraction for optical frequencies would generate revolutionary impacts on present-day optical technologies. However, there exist some difficulties in the procedure of fabrication and it is still a challenge expected to be overcome at present.

Up to now, a relatively typical and extensively studied system is the so-called fishnet structure consisting of two perforated metal films separated by a thin layer of dielectric [7,8,11], pushing negative refraction up to the red-light wavelengths. This structure is often fabricated with the ‘top-down’ approach such as electron-beam lithography [12] or focused ion-beam technology [13]. According to the effective medium theory, the interior structure of the medium must be on a sub-wavelength scale. This requires that the structural unit of artificial metamaterials should be constructed on the nanometer scale when operating at visible wavelengths. However, in current time this ‘top-down’ method could hardly meet this small precision, and moreover it often needs expensive equipments and involves complicated processes. Over the past several years, our group has proposed a ‘bottom-up’ approach of double template-assisted electrochemical deposition to fabricate NIMs which is convenient to operate and also can satisfy the nano-scale precision. In our previous work, the two-dimensional (2D) isotropic dendritic NIM has been successfully fabricated by using the ‘bottom-up’ approach and the negative refraction has been experimentally demonstrated at infrared frequencies. We have performed explicit study on 2D dendritic NIM. Negative permeability and subwavelength focusing of quasi-periodic dendritic cell metamaterial were studied in [14]. The dendritic structure with simultaneously negative permeability and permittivity has been verified [15–17]. The infrared dendritic NIM was successfully fabricated by means of double template-assisted electrochemical deposition [18]. NIM composed of random dendritic cells displays a multi-band resonance, negative refractive index and a high intensity of the pass-band at infrared frequencies [19]. Recently we have experimentally demonstrated the trapped rainbow in tapered negative-index heterostructures at visible frequencies [20].

In addition to realizing negative refraction at optical frequencies, designing and fabricating three-dimensional (3D) NIMs is another challenging task. Until now, few 3D models have been reported, the 3D fishnet structure in [21] which can achieve negative refraction in the infrared part of the spectrum has attracted extensive attention of researchers worldwide. But it is not perfect for its absence of isotropy and the complexity in the procedure of fabrication. In this paper, we design a 3D isotropic NIM composed of sphere-rod shaped cells. We have performed numerical simulations with adjusting its geometrical parameters to study its electromagnetic response to the applied electromagnetic fields. It is found that this structure can yield a negative refractive index in the visible region. This model is so highly symmetric about itself that it possesses the property of isotropy, hence it is more propitious for actual applications. Because of the simplicity of the 3D structure and the convenient ‘bottom-up’ approach of fabricating 2D dendrtic NIM mentioned above, we are now attempting to fabricate this sphere-rod shaped NIM by means of chemical synthesization. So the purpose of this paper is to present a new 3D model for NIMs and also afford the theoretic support for after ongoing experiment.

## 2. Numerical simulation

Silver is chosen for the presented structure, because of its significantly low losses in the visible range, as well as easy to be fabricated with chemical synthesization. Numerical simulations are performed with a 3D electromagnetic simulative software based on the finite-difference time-domain method. At lower frequencies such as the gigahertz range, the silver behaves as nearly perfect conductor. Whereas in the terahertz range, its dielectric function should be characterized by Drude model with plasma frequency ω_{p} = 1.37 × 10^{16} s^{−1} and collision frequency ω_{c} = 8.5 × 10^{13} s^{−1} [22]. The electromagnetic radiation is set as the transverse electric and magnetic (TEM) mode as shown in Fig. 1
. Because the unit cells are periodically repeated in the 3D space, we simulate one unit cell by applying the perfect electric and magnetic boundary conditions. After obtaining the reflection (S_{11}) and transmission (S_{21}) coefficients, we first extract refractive index n and impedance z from the two parameters, then the effective permittivity ε and the effective permeability μ can be directly calculated from μ = nz and ε = n/z [23].

The structural unit and its geometric dimensions are depicted in Fig. 1. It is a 3D sphere-rod shaped structure with the rods distributing uniformly around the nucleus. There are four geometrical parameters (L, r, t, a) on the order of nanometer, characterizing the length of the rod originated from the central nucleus, the radius of the nucleus which is assumed to be spherical, the diameter of the cross section at the end of the rod, and the lattice constant, respectively.

## 3. Results and analysis

#### 3.1 Parametric optimization

Recently, parametric study shows that the spectra of reflection and transmission depend strongly on the structural parameters of the periodic structures [24]. Optimization of the shape and size of metamaterial is essential to get an appropriate structure with better electromagnetic performance. We have studied the dependence of the reflection and transmission spectra on the structural parameters. In this 3D model, among the four geometric parameters (L, r, t, a) is the radius of the nucleus r playing a primary role in affecting the electromagnetic response of the structure, this is because that r is directly associated with the overall morphology of the model. In simulation, we start with choosing appropriate parameters L, t, a, and then change r from zero to the value of the nucleus fully covering the rods. Figures 2(a) and 2(b) display the curves of the retrieved permeability of the structure composed of only the rods and only the nucleus, respectively. It can be seen that the real parts of the permeability in these two cases are positive over the entire spectral bands, it obviously cannot give rise to a negative refractive index. Next, we fix a proper value of r and respectively change other three parameters (L, t, a) in order to get a primitive model. At last synthetically considering the effects of the four parameters, and further adjusting them is to find out the expected structure exhibiting a negative index. This procedure is intricate and time-consuming, however, this paper is purely focused on displaying the negative refraction in the optical domain, so we only take a brief depiction on adjusting parameters in this paragraph.

#### 3.2 Response of visible region

Based on the strategy of geometric optimization mentioned above, rationally adjusting the four parameters (L, r, t, a) can give rise to a negative refractive index from the red-light to the blue-light region. Owing to space limitations, the detailed and complicated adjustment of parameters are not displayed. Below, the spectra of electromagnetic response to the red-light and to the blue-light are displayed in Fig. 3 and Fig. 4 , respectively.

#### 3.2.1 Negative refraction for red-light

In this case, the specific dimensions of the structure are as follows: L *=* 65 nm, r *=* 39 nm, t *=* 19 nm, a *=* 140 nm. In Fig. 3 (a), it can be seen that there is a transmission peak of S_{21} and a reflection trough of S_{11} at bout 436.3 THz. In the vicinity of this frequency, the real part of the refractive index is negative with simultaneously negative real parts of permeability and permittivity as shown in Figs. 3(b), 3(c), 3(d). The frequency band of the negative refractive index is within the red part of the spectrum (630-780nm). The ratio of the operation wavelength to the lattice constant is equal to 4.9, satisfying the effective medium theory.

#### 3.2.2 Negative refraction for blue-light

The following numerical results for the structure with L *=* 50 nm, r *=* 39 nm, t *=* 18 nm, a *=* 110 nm. In Fig. 4(a), there is a transmission peak with a high value close to 0.6 and the corresponding reflection trough has a low value of less than 0.05 at about 608.6 THz. As the same of Fig. 3, in the neighboring area of this frequency all the real parts of permeability, permittivity and refractive index simultaneously have a negative value, as illustrated in Figs. 4(b), 4(c), 4(d). The ratio of the operation wavelength over the lattice constant is 4.5, so the system can be regarded as an effective medium. The negative refraction is achieved within the wavelength rang of blue-light (470-500nm).

As Fig. 3 and Fig. 4 have been shown, let us analyze the underlying physical mechanism. The electric response is attributed to the electric dipole oscillation. With regard to magnetic resonance, there are two arc currents exist between the two neighboring rods of either side of the structure, every arc current excites a magnetic fields (Fig. 5 ), which generates a magnetic field opposing the original field when its phase is reversed. This diamagnetic response yields a negative permeability, further combining the negative permittivity at the same frequencies leads to a negative refractive index. In the other hand, this sphere-rod shaped structure can be thought of as a 3D dendritic structure with a nucleus, its detailed theoretical analysis is similar to that of the 2D dendritic structure mentioned above.

Further observing Fig. 3 and Fig. 4, we note here that the spectral width of negative index is broader than that of the simultaneously negative permeability and permittivity in the lower band of the overlapped frequencies. This can be explained that designing an artificial negative-index metamaterial usually requires the real parts of the permeability (μ = μ' + iμ”) and the permittivity (ε = ε' + iε”) are simultaneously negative over a certain frequency band. But strictly speaking, it is not a necessary condition. There exists a more relaxed relationship P = μ'ε” + μ”ε' < 0 can also achieve a negative refractive index [25]. So we calculated the value of P for the two cases as shown in Fig. 6 , it can be found that the wavelength range of the negative refractive index is perfectly corresponding to that of the negative P. The curves of the impedance of the two cases can also interpret this point, as depicted in Fig. 7 . The real part of impedance more than zero implies this area is a pass-band which agrees with the negative refraction pass-band, the other regions of zero value of impedance are opaque bands which correspond with the single negative value of the permeability or the permittivity.

## 4. Conclusions

In summary, we propose a new 3D optical NIM consisting of sphere-rod shaped cells. Under the condition that the ratio of lattice constant to operating wavelength meets the requirement of the effective medium theory, adjusting its geometrical parameters can make the structure give rise to a negative refractive index. By further scaling down and optimizing its geometric parameters, it is found that the response of negative refraction can be shifted from red-light to blue-light frequencies. The spectral band of the calculated negative refractive index is broader than that of simultaneously negative permeability and permittivity, this phenomenon can be explained by calculating the value of P. The model has the advantage of simplicity and isotropy, it is expected to be obtained from the ‘bottom-up’ approach of chemical synthesization, which also offers a feasible route to massive and low-cost fabrication of optical NIMs.

## Acknowledgments

We acknowledge support from the National Natural Science Foundation of China under Grant No. 50632030, 50872113, 50936002.

## References and links

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

**2. **A. Alù and N. Engheta, “Achieving transparency with plasmonic and metamaterial coatings,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **72**(1), 016623 (2005). [CrossRef] [PubMed]

**3. **R. A. Shelby, D. R. Smith, S. C. Nemat-Nasser, and S. Schultz, “Microwave transmission through a two-dimensional, isotropic, left-handed metamaterial,” Appl. Phys. Lett. **78**(4), 489–491 (2001). [CrossRef]

**4. **N. Seddon and T. Bearpark, “Observation of the inverse Doppler effect,” Science **302**(5650), 1537–1540 (2003). [CrossRef] [PubMed]

**5. **J. Lu, T. M. Grzegorczyk, Y. Zhang, J. Pacheco Jr, B. I. Wu, J. A. Kong, and M. Chen, “Cerenkov radiation in materials with negative permittivity and permeability,” Opt. Express **11**(7), 723–734 (2003). [CrossRef] [PubMed]

**6. **R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science **292**(5514), 77–79 (2001). [CrossRef] [PubMed]

**7. **G. Dolling, C. Enkrich, M. Wegener, C. M. Soukoulis, and S. Linden, “Low-loss negative-index metamaterial at telecommunication wavelengths,” Opt. Lett. **31**(12), 1800–1802 (2006). [CrossRef] [PubMed]

**8. **G. Dolling, C. Enkrich, M. Wegener, C. M. Soukoulis, and S. Linden, “Simultaneous negative phase and group velocity of light in a metamaterial,” Science **312**(5775), 892–894 (2006). [CrossRef] [PubMed]

**9. **G. Dolling, M. Wegener, and S. Linden, “Realization of a three-founctional-layer negative-index photonic metamaterial,” Opt. Lett. **32**(5), 551–553 (2007). [CrossRef] [PubMed]

**10. **U. K. Chettiar, A. V. Kildishev, H. K. Yuan, W. Cai, S. Xiao, V. P. Drachev, and V. M. Shalaev, “Dual-band negative index metamaterial: double negative at 813 nm and single negative at 772 nm,” Opt. Lett. **32**(12), 1671–1673 (2007). [CrossRef] [PubMed]

**11. **S. Zhang, W. J. Fan, N. C. Panoiu, K. J. Malloy, R. M. Osgood, and S. R. J. Brueck, “Experimental demonstration of near-infrared negative-index metamaterials,” Phys. Rev. Lett. **95**(13), 137404 (2005). [CrossRef] [PubMed]

**12. **C. Enkrich, M. Wegener, S. Linden, S. Burger, L. Zschiedrich, F. Schmidt, J. F. Zhou, Th. Koschny, and C. M. Soukoulis, “Magnetic metamaterials at telecommunication and visible frequencies,” Phys. Rev. Lett. **95**(20), 203901 (2005). [CrossRef] [PubMed]

**13. **C. Enkrich, F. Perez-Williard, D. Gerthsen, J. Zhou, T. Koschny, C. M. Soukoulis, M. Wegener, and S. Linden, “Focused-ion-beam nanofabrication of near-infrared magnetic metamaterials,” Adv. Mater. **17**(21), 2547 (2005). [CrossRef]

**14. **X. Zhou, Q. H. Fu, J. Zhao, Y. Yang, and X. P. Zhao, “Negative permeability and subwavelength focusing of quasi-periodic dendritic cell metamaterials,” Opt. Express **14**(16), 7188–7197 (2006). [CrossRef] [PubMed]

**15. **Y. Yao and X. P. Zhao, “Multilevel dendritic structure with simultaneously negative permeability and permittivity,” J. Appl. Phys. **101**(12), 124904 (2007). [CrossRef]

**16. **Y. Yao, Q. H. Fu, and X. P. Zhao, “Three-level dendritic structure with simultaneously negative permeability and permittivity under normal incidence of electromagnetic wave,” J. Appl. Phys. **105**(2), 024911 (2009). [CrossRef]

**17. **X. Zhou and X. P. Zhao, “Resonant condition of unitary dendritic structure with overlapping negative permittivity and permeability,” Appl. Phys. Lett. **91**(18), 181908 (2007). [CrossRef]

**18. **H. Liu, X. P. Zhao, Y. Yang, Q. W. Li, and J. Lv, “Fabrication of infrared left-handed metamaterials via double template-assisted electrochemical deposition,” Adv. Mater. **20**(11), 2050–2054 (2008). [CrossRef]

**19. **B. Q. Liu, X. P. Zhao, W. R. Zhu, W. Luo, and X. C. Cheng, “Multiple Pass-Band Optical Left-Handed Metamaterials Based on Random Dendritic Cells,” Adv. Funct. Mater. **18**(21), 3523–3528 (2008). [CrossRef]

**20. **X. P. Zhao, W. Luo, J. X. Huang, Q. H. Fu, K. Song, X. C. Cheng, and C. R. Luo, “Trapped rainbow effect in visible light left-handed heterostructures,” Appl. Phys. Lett. **95**(7), 071111 (2009). [CrossRef]

**21. **J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, G. Bartal, and X. Zhang, “Three-dimensional optical metamaterial with a negative refractive index,” Nature **455**(7211), 376–379 (2008). [CrossRef] [PubMed]

**22. **P. B. Johnson and R. W. Christy, “Optical Constants of the Noble Metals,” Phys. Rev. B **6**(12), 4370–4379 (1972). [CrossRef]

**23. **X. D. Chen, T. M. Grzegorczyk, B.-I. Wu, J. Pacheco Jr, and J. A. Kong, “Robust method to retrieve the constitutive effective parameters of metamaterials,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **70**(1), 016608 (2004). [CrossRef] [PubMed]

**24. **V. D. Lam, J. B. Kim, S. J. Lee, and Y. P. Lee, “Left-handed behavior of combined and fishnet structures,” J. Appl. Phys. **103**(3), 033107 (2008). [CrossRef]

**25. **R. A. Depine and A. A. Lakhtakia, “New condition to identify isotropic dielectric-magnetic materials displaying negative phase velocity,” Microw. Opt. Technol. Lett. **41**(4), 315–316 (2004). [CrossRef]