Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5547121

Formats

Article sections

Authors

Related links

Sci Rep. 2017; 7: 7445.

Published online 2017 August 7. doi: 10.1038/s41598-017-07914-1

PMCID: PMC5547121

Xiaodong Huang, Email: ua.ude.niws@gnauhx.

Received 2017 May 2; Accepted 2017 July 4.

Copyright © The Author(s) 2017

All-angle negative refraction (AANR) of phononic crystals and its frequency range are dependent on mechanical properties of constituent materials and their spatial distribution. So far, it is impossible to achieve the maximum operation frequency range of AANR theoretically. In this paper, we will present a numerical approach for designing a two-dimensional phononic crystal with broadband AANR without negative index. Through analyzing the mechanism of AANR, a topology optimization problem aiming at broadband AANR is established and solved by bi-directional evolutionary structural optimization method. The optimal steel/air phononic crystal exhibits a record AANR range over 20% and its refractive properties and focusing effects are further investigated. The results demonstrate the multifunctionality of a flat phononic slab including superlensing effect near upper AANR frequencies and self-collimation at lower AANR frequencies.

Since Veselago^{1} first predicted the left-hand material (LHM) with simultaneously negative permittivity and permeability, artificial materials with negative refraction have attracted considerable attention due to their interesting and unusual behaviors. By virtue of double negative effective index, super resolution beyond the traditional diffraction limit was first achieved for electromagnetic waves by Pendry with a flat optical lens, called a superlens^{2}. On the other hand, negative refraction of electromagnetic waves can also result from Bragg scattering effects of photonic crystals, where the effective index is still positive^{3, 4}. Similar to the perfect lens consisting of LHMs, negative refraction in photonic crystals can also lead to a flat superlens^{5, 6}.

Sparked by these pioneering works on electromagnetic waves, increasing efforts have been devoted to exploring similar phenomena in phononic crystals (PnCs) and acoustic metamaterials theoretically and experimentally, such as achieving negative refraction^{7–19} and reflection^{20} of acoustic waves, and very recently realizing acoustic topological insulators^{21–23}. These explorations enable us to manipulate the propagation of acoustic waves better and are promising for engineering devices with unconventional functions, for example, focusing of sound waves with super resolution. Compared with acoustic metamaterials whose lattice constant is unusually two orders shorter than the relevant wavelength, PnCs are more preferable to achieve subwavelength focusing at the ultrasound regime as their lattice constant is comparable to the wavelength. There are two different approaches to realize negative refraction in PnCs. One is exploiting the natural opposite direction of wave vector and energy flow at high-order band, typically at the second band, which has been extensively discussed in the literature^{10, 11, 13, 15, 16}. This mechanism comes along with a backward-wave effect and the negative effective index, and is similar to the negative refraction by metamaterials. The other is achieving negative refraction at the first phononic band by intensive scattering of constituents, resulting in positive effective properties and PnCs that behave like uniform right-handed materials (RHMs)^{8, 9}. Such designs have advantages of single mode, high transmission and possible all-angle negative refraction (AANR), which is of particular interest in this paper. To our knowledge, AANR at the first band was only reported in a 2D square lattice of mercury/water system with a relative AANR frequency range about 6% around frequency 0.24*2*πc*
_{water}/*a*
^{8}.

A critical question is how to attain the maximum operation frequency range of AANR so as to realize its full potential. AANR in PnCs is essentially a band structure controlled behavior, thus influenced by mechanical properties of the constituent materials, and their spatial distributions (or called structure) within the primitive unit cell^{24–27}. Due to the infinite possibility of material distribution, it is impossible to obtain the best structure for AANR from a physical point of view, yet it is anticipated to answer this question by a multidisciplinary approach. In this paper, we establish a numerical approach to achieve AANR behavior in PnCs by integrating the underlying physical mechanism into a topology optimization problem. To seek the best configuration of PnCs with broadband AANR, we employ the bi-directional evolutionary structural optimization (BESO) method which has been proven to be very effective and robust in the design of photonic/phononic band gap crystals, as well as photonic structures with AANR^{28–32}.

The AANR behavior originates from the relation of equifrequency contours (EFCs) for the PnC and air which defines the lower and upper frequency limits. At the first phononic band, the EFCs for the PnC alter from concave to convex as frequencies get higher in the vicinity of *M* point while the EFCs for air are concentric circles centered at *Γ* point with the radius Ω/*c*
_{air}. For these convex EFCs of the PnC, the group velocity pointing along their outward normal will consequently point inwards to *M* as demonstrated in Fig. 1a. As a result, both the incidence and refraction will stand on the same side of the interface normal, i.e. negative refraction is achieved. Thus the alternation point gives the lower limit of AANR Ω_{l}. Notice that the group velocity at first band is never opposite to the wave vector, which differs from the LHMs and negative refraction realized in high-order bands. Therefore the realized negative refraction does not come along with either negative index nor backward wave effects. The ultimate case for AANR is that the EFCs for the PnC and air have the same diameter, so as to guarantee negative refraction for incident beams entering PnC from any direction. This condition defines the upper frequency limit, Ω_{u} of AANR, and AANR can only be achieved when Ω_{u} > Ω_{l}
^{9}. Only PnCs with specially designed inclusions can satisfy these strict conditions, which is a typical topology optimization problem.

Our optimization starts from the popular design consisting of steel cylinders in a square lattice. The filling fraction is 50%. Note that there is no AANR behavior found in the initial configuration. BESO gradually updates the topology and enlarges the AANR frequency range as shown in Supplementary Fig. ^{S2a}. The final configuration along with its band structure and EFCs are presented in Fig. 1b–d. The results demonstrate that a wide AANR frequency range around Ω=0.193~0.237 is obtained, and its relative width, 2(Ω_{u}−Ω_{l})/(Ω_{u}+Ω_{l}), achieves 20.35%, representing a record value in the literature. The optimized shape of steel has changed significantly from the original circle to a concave curve, and can be approximately characterized by two parameters *R*=0.34*a* and *t*=0.03*a*, where *a* is the square lattice constant. The assigned 2*t* air gap between unit cells ensures that only longitudinal waves are supported. The radius of the EFC curvature at the point where EFC and the *ΓM* line intersect is plotted together with the first phononic band and air dispersion line in Fig. 1c. Frequencies Ω are normalized by 2*πc*
_{air}/*a*. The singular point where the radius of EFC curvature goes to positive and negative infinity on two sides gives the lower limit Ω_{l} (indicated as upward-pointing triangles), while Ω_{u} (indicated as downward-pointing triangles) is given by the intersection of the air line and the first band. Initially, AANR behavior does not exist as Ω_{l}>Ω_{u} (see blue triangles). With the evolution of the PnC structure, the first phononic band gets depressed, the wave vector corresponding to the Ω_{u} gets closer to *M* point. At the same time, the wave vector of Ω_{l} approaches *Γ* point. Consequently, the AANR frequency range is enlarged.

We further investigated the influences of parameters *R* and *t* on AANR for the given topology shown in Fig. 1b. The numerical parametic studies revealed that the assigned air gap 2*t* can substantially tune AANR behavior for a fixed *R*. Figure 2 shows the AANR frequency range as a function of the air gap with different *R*. All the AANR frequency ranges decline monotonically with the increase of the air gap (see the straight fitted lines), and finally disappear. In fact, it is even impossible to obtain AANR by means of topology optimization when the air gap is too wide. In comparison, *R* merely shows a slight influence on AANR when the air gap is fixed. For instance, when *R* increases from 0.1*a* to 0.45*a* along with the fixed air gap 2*t*=0.625*a*, the AANR behavior persists with its frequency range varying from 9.55% to 20.35%. We envision that the simple cross structure with any thickness would have a positive AANR frequency range as long as the air gap is carefully chosen. This finding deepens our understanding of the relation between the phononic structure and its AANR behavior and provides useful instructions for practical applications.

Previous research suggested that AANR does not exist in the initial design at a filling fraction lower than 60%, which means that a high filling fraction has to be employed to realize AANR without tailoring the structure^{9}. However, our optimization results indicate that the modification of the geometry can effectively bring about a wide AANR frequency range even at a filling fraction as low as 20% (see Fig. 3a), showing that the AANR is highly dependent on the geometry of PnCs rather than the filling fraction of inclusion. The air gap in all the cases is set to be a constant value 2*t*=0.625*a*. As the filling fraction decreases from 80% to 20%, both the upper and lower limit first drops and then goes back to a higher level. Among all the cases, the widest AANR is achieved at a filling fraction of 50%. The corresponding optimized unit cells are presented in Fig. 3b–g. It is interestingly observed that all the optimized geometries of inclusion at different filling fractions have concave boundaries. With the decreasing of filling fraction, the optimized topology of inclusions finally evolves to a cross-shaped structure. These concave profiles may contribute to the existence of AANR behavior and need further investigations.

In the following section, we will validate the AANR phenomenon by investigating the refractive properties of the optimized PnC and show that the optimized PnC can be constructed as a multifunctional flat lens which can collimate acoustic waves and focus an acoustic source into an image with a resolution beyond the diffraction limit.

Due to the high anisotropy of EFCs for the optimized PnC, the refracted angle depends on the angle of incidence and frequency, as well as the interface normal. Provided that the interface normal is oriented along *ΓM* direction, the dependence of the refracted angle, *β*, on the incident angle, *α*, at different frequencies is presented in Fig. 4. Negative refracted angles indicate that AANR is achieved in all cases. As the incident angle enlarges from 0° to 90°, the refracted angle changes slightly from 0° to 5.5° for frequency Ω=0.200 just above the lower limit while it increases greatly from 0° to 64.7° for Ω=0.230. This phenomenon can be well understood by examining the EFCs of the PnC (see Fig. 1d). The EFC at Ω_{l} is close to a square and evolves to circular when frequency approaches Ω_{u}. The steeper slope of the curve for Ω=0.230 verifies the fact that EFCs at a higher frequency tends to be more isotropic. The inset of Fig. 4 is the simulation result of a plane wave across a 10-layer PnC slab with an incident angle of 45° at frequency Ω=0.220, which clearly demonstrates the negative refraction.

The refracted angle as a function of incident angle at different normalized frequencies. The inset is simulation result for a plane wave through a 10-layer PnC slab with incident angle 45° at frequency Ω=0.220.

To check the focusing effect of the optimized PnC, an 8-layer PnC slab consisting of 100unit cells is constructed with interface normal along *ΓM* direction. The point source is placed at a distance *d*=*a*, close to the left surface of the slab. Absorbing boundary conditions are adopted around the whole simulation domain to avoid reflections. The normalized intensity fields at frequency Ω=0.200, 0.230 and 0.237 are presented in Fig. 5. For the case of Ω=0.200, there is no clear image on the right side of the slab. However, a clear guiding channel is visible inside the phononic slab, which demonstrates a self-collimation effect. The reason is that the refracted angles for all incidences are less than 6°. Thus it may need a wider lens to form an evident image as the refracted waves must intersect in the phononic slab (see Supplementary Fig. ^{S3})^{3}. Stronger self-collimation effect can be observed as frequency approaches the lower limit Ω=0.193 due to the extremely flat portion of EFC (see Supplementary Fig. ^{S4}).

Normalized intensity field of a point source and its image across an 8-layer phononic slab at normalized frequency Ω=0.200 (**a**), Ω=0.230 (**b**) and Ω=0.237 (**c**). The dark (light) **...**

On the other hand, a clear image with high intensity at frequency Ω=0.230 has been formed at an approximate distance of 0.6*a* away from the right slab surface. It locates on the axis of the point source without any deviation. The focusing resolution, which is defined as the full width at half maximum of the intensity peak along the *y* direction, is around 0.46*λ*, where *λ*=*a*/Ω is the sound wavelength in air. It reveals that the super resolution beyond the Rayleigh diffraction limit 0.5*λ* has been achieved by the optimized phononic structure. Similarly, a good focusing effect with super resolution 0.37*λ* has been observed at the upper bound of AANR frequency range Ω_{u}=0.237 as well. Such good lensing phenomena benefit from the so-called canalization mechanism^{33}. When the EFC of PnC is larger than that of air, a great proportion of evanescent modes of the point source can be canalized to the Bloch modes inside the PnC and transported across the slab without attenuation. All the propagating Bloch modes convert back to the corresponding propagating modes in PnC and air on the exit surface. The substantial canalized evanescent waves attribute to enhance the resolution of focusing. The canalization mechanism is different from the amplification of evanescent waves in the second band and LHMs that are achieved by backwards wave effect introduced phase compensations^{13}. As the near circular EFC at higher frequencies of AANR has more uniformly distributed refractive index regardless of the incident direction, the incidences emitted from the point source will be refracted into the same position within the PnC and form a perfect image on the exit side. Meanwhile, the single beam behaviour at frequencies lower than 0.5×2*πc*
_{air}/*a*
_{s} (where ${a}_{s}=\sqrt{2}a$) ensures the high transmission, which is also an important factor for the super resolution.

It should be pointed out that the thickness and surface termination of the slab can strongly tune the focusing frequency range, which is outside the scope of this letter. As for this 8-layer PnC slab, it can focus a clear spot over frequency range Ω=[0.215,0.237] (see Supplementary Fig. ^{S5}). The intensity distributions along the direction perpendicular (*x*) and parallel (*y*) to the slab surface through the peak of focusing point are shown in Fig. 6. Remind that the sound wavelength in air *λ*=*a*/Ω>4*a* as Ω<0.25 in all the cases. It is evident that subwavelength imaging has been achieved for these frequencies. In fact the calculated lateral resolution for these frequencies varies from 0.37*λ* to 0.53*λ*. The focal distances are around *a* from the right lens surface, remaining in the near-field domain. These results demonstrate that the optimized PnC enables us to realize self-collimation and subwavelength focusing of acoustic waves simultaneously with a simple flat lens, and has the advantage of working as a superlens over a broad frequency range.

In summary, we have introduced the topology optimization method for constructing a 2D solid/air PnC displaying AANR phenomenon without negative index. The resulting PnC at a filling fraction 50% exhibits a record AANR frequency range of 20.35% near frequency Ω=0.215. All the optimized PnCs have concave profiles at material interface which may help to open an effective AANR frequency range. By virtue of the AANR behavior, subwavelength focusing with the super resolution was achieved by a phononic slab near the AANR upper limit while self-collimation took place near the AANR lower limit. The wide frequency range of AANR and subwavelength focusing effect makes the proposed PnC a desirable candidate for the design of acoustic superlens.

The 2D PnC considered in this paper is composed of steel inclusion in air background (*ρ*
_{steel}=7800*kg*/*m*
^{3}, *ρ*
_{air}=1.21*kg*/*m*
^{3}, *c*
_{steel}=6100*m*/*s*, *c*
_{air}=334.5*m*/*s*)^{9}, although our approach is general for any type of PnCs. Here, finite element analysis (FEA) is used to calculate the band structure and equifrequency contours since any irregular shape of the inclusion is possible during the optimization process. In the calculation, we ignore the shear modulus of the solid and consider it as metafluid as the solid is nearly rigid in comparison to air^{34}. Such assumption was also validated in the simulation of acoustic focusing effect as presented in Fig. ^{S6} in the Supplementary Materials. Numerical simulations are carried out by COMSOL Multiphysics 5.1 throughout this paper.

The objective is to find the material distribution of PnCs with the maximal frequency range of AANR, which can be mathematically expressed by

$$\begin{array}{ll}\hfill \mathrm{Find}:& \phantom{\rule{1em}{0ex}}\mathbf{X}=[{x}_{1},\phantom{\rule{.25em}{0ex}}{x}_{2},\mathrm{...,}\phantom{\rule{.25em}{0ex}}{x}_{e},\mathrm{...,}\phantom{\rule{.25em}{0ex}}{x}_{N}]\phantom{\rule{1em}{0ex}}{x}_{e}\in \mathrm{[0,1]}\hfill \\ \mathrm{Maximize}:\hfill & \phantom{\rule{1em}{0ex}}f(\mathbf{X})={\mathrm{\Omega}}_{u}(\mathbf{X})-{\mathrm{\Omega}}_{l}(\mathbf{X})\hfill \end{array}$$

1

where *x*
_{e} (*e*=1, 2, ..., *N*) is the design variable of element *e*, where *x*
_{e}=0 represents element *e* is composed of air and *x*
_{e}=1 denotes of steel. *N* is the total number of elements within the primitive unit cell. Thus, the discrete values of design variables can solely represent the material distribution of the unit cell. It should be noted that numerically extracting the lower and upper frequency limits of AANR provides the key for topology optimization.

To establish an optimization algorithm, sensitivity analysis is necessary for evaluating the impact of each element on the frequency range of AANR, which will be used as the optimization criterion. A higher positive value of an elemental sensitivity, i.e. the derivative of the objective function, means that increasing the design variable of the element (physically switching material from air to steel) will increase the AANR frequency range. Following this understanding, BESO will increase design variables (from 0 to 1) for elements with high sensitivities and simultaneously decrease design variables (from 1 to 0) for elements with low sensitivities. As a result, the updated design variables for all elements within the unit cell form a new geometry. The whole optimization procedure including the extraction of the AANR frequency limits, sensitivity analysis, and BESO update scheme is conducted iteratively until the material distribution within the unit cell achieves its optimum. Meanwhile, the AANR frequency range gradually increases from negative, positive to the maximum value. We refer the interested readers to the Supplementary information for the detailed optimization formulation and procedure.

The authors would like to acknowledge the support of Australian Research Council for providing funding (FT130101094).

Author Contributions

Y.F.L. and F.M. conducted the numerical simulations and optimization. X.H. and S.Z. supervised the study, M.H.L. contributed to discussion of the results. The paper was written by Y.F.L. and reviewed and approved by all authors.

The authors declare that they have no competing interests.

**Electronic supplementary material**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-07914-1

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

1. Viktor, G. V. The electrodynamics of substances with simultaneously negative values of and µ. *Soviet Physics Uspekhi***10**, 509 (1968).

2. Pendry JB. Negative refraction makes a perfect lens. Phys. Rev. Lett. 2000;85:3966–3969. doi: 10.1103/PhysRevLett.85.3966. [PubMed] [Cross Ref]

3. Luo, C., Johnson, S. G., Joannopoulos, J. D. & Pendry, J. B. All-angle negative refraction without negative effective index. *Phys*. *Rev*. *B***65**, 2011041–2011044, doi:10.1103/PhysRevB.65.201104.

4. Cubukcu E, Aydin K, Ozbay E, Foteinopoulou S, Soukoulis CM. Negative refraction by photonic crystals. Nature. 2003;423:604–605. doi: 10.1038/423604b. [PubMed] [Cross Ref]

5. Parimi PV, Lu WT, Vodo P, Sridhar S. Imaging by flat lens using negative refraction. Nature. 2003;426:404. doi: 10.1038/426404a. [PubMed] [Cross Ref]

6. Li J, et al. All-angle negative refraction imaging effect with complex two-dimensional hexagonal photonic crystals. J. Appl. Phys. 2007;102:073538. doi: 10.1063/1.2794860. [Cross Ref]

7. Yang S, et al. Focusing of sound in a 3d phononic crystal. Phys. Rev. Lett. 2004;93:024301–1. doi: 10.1103/PhysRevLett.93.024301. [PubMed] [Cross Ref]

8. Zhang X, Liu Z. Negative refraction of acoustic waves in two-dimensional phononic crystals. Appl. Phys. Lett. 2004;85:341–343. doi: 10.1063/1.1772854. [Cross Ref]

9. Feng L, et al. Negative refraction of acoustic waves in two-dimensional sonic crystals. Phys. Rev. B. 2005;72:033108. doi: 10.1103/PhysRevB.72.033108. [Cross Ref]

10. Ke, M. *et al*. Negative-refraction imaging with two-dimensional phononic crystals. *Phys*. *Rev*. *B***72**, 064306, doi:10.1103/PhysRevB.72.064306 (2005).

11. Qiu, C., Zhang, X. & Liu, Z. Far-field imaging of acoustic waves by a two-dimensional sonic crystal. *Phys*. *Rev*. *B***71**, 054302, doi:10.1103/PhysRevB.71.054302 (2005).

12. Fang N, et al. Ultrasonic metamaterials with negative modulus. Nature Materials. 2006;5:452–456. doi: 10.1038/nmat1644. [PubMed] [Cross Ref]

13. Feng, L. *et al*. Acoustic backward-wave negative refractions in the second band of a sonic crystal. *Phys*. *Rev*. *Lett*. **96**, 014301, doi:10.1103/PhysRevLett.96.014301 (2006). [PubMed]

14. Ding, Y., Liu, Z., Qiu, C. & Shi, J. Metamaterial with simultaneously negative bulk modulus and mass density. *Phys*. *Rev*. *Lett*. **99**, 093904, doi:10.1103/PhysRevLett.99.093904 (2007). [PubMed]

15. Sukhovich, A., Jing, L. & Page, J. H. Negative refraction and focusing of ultrasound in two-dimensional phononic crystals. *Phys*. *Rev*. *B***77**, 014301, doi:10.1103/PhysRevB.77.014301 (2008).

16. Sukhovich, A. *et al*. Experimental and theoretical evidence for subwavelength imaging in phononic crystals. *Phys*. *Rev*. *Lett*. **102**, 154301, doi:10.1103/PhysRevLett.102.154301 (2009). [PubMed]

17. Zhang, S., Yin, L. & Fang, N. Focusing ultrasound with an acoustic metamaterial network. *Phys*. *Rev*. *Lett*. **102**, 194301, doi:10.1103/PhysRevLett.102.194301 (2009). [PubMed]

18. Robillard, J. F. *et al*. Resolution limit of a phononic crystal superlens. *Phys*. *Rev*. *B***83**, 224301, doi:10.1103/PhysRevB.83.224301 (2011).

19. Song GY, Huang B, Dong HY, Cheng Q, Cui TJ. Broadband focusing acoustic lens based on fractal metamaterials. Scientific Reports. 2016;6:35929. doi: 10.1038/srep35929. [PMC free article] [PubMed] [Cross Ref]

20. Zhao D, Ye Y, Xu S, Zhu X, Yi L. Broadband and wide-angle negative reflection at a phononic crystal boundary. Appl. Phys. Lett. 2014;104:043503. doi: 10.1063/1.4863691. [Cross Ref]

21. He C, et al. Acoustic topological insulator and robust one-way sound transport. Nat Phys. 2016;12:1124–1129. doi: 10.1038/nphys3867. [Cross Ref]

22. Peng Y-G, et al. Experimental demonstration of anomalous floquet topological insulator for sound. Nature Communications. 2016;7:13368. doi: 10.1038/ncomms13368. [PMC free article] [PubMed] [Cross Ref]

23. Lu J, et al. Observation of topological valley transport of sound in sonic crystals. Nat Phys. 2017;13:369–374. doi: 10.1038/nphys3999. [Cross Ref]

24. Kushwaha MS, Halevi P, Dobrzynski L, Djafari-Rouhani B. Acoustic band structure of periodic elastic composites. Phys. Rev. Lett. 1993;71:2022–2025. doi: 10.1103/PhysRevLett.71.2022. [PubMed] [Cross Ref]

25. Kushwaha MS, Halevi P, MartÃnez G, Dobrzynski L, Djafari-Rouhani B. Theory of acoustic band structure of periodic elastic composites. Phys. Rev. B. 1994;49:2313–2322. doi: 10.1103/PhysRevB.49.2313. [PubMed] [Cross Ref]

26. Kushwaha MS. Classical band structure of periodic elastic composites. Int. J. Mod. Phys. B. 1996;10:977–1094. doi: 10.1142/S0217979296000398. [Cross Ref]

27. Sigalas M, Economou EN. Band structure of elastic waves in two dimensional systems. Solid State Commun. 1993;86:141–143. doi: 10.1016/0038-1098(93)90888-T. [Cross Ref]

28. Huang, X. & Xie, M. *Evolutionary topology optimization of continuum structures: methods and applications* (John Wiley & Sons, Chichester, 2010).

29. Li, Y. F., Huang, X., Meng, F. & Zhou, S. Evolutionary topological design for phononic band gap crystals. *Struct*. *Multidisc*. *Optim*. *54*, 595–617, doi:10.1007/s00158-016-1424-3 (2016).

30. Li YF, Huang X, Zhou S. Topological design of cellular phononic band gap crystals. Materials. 2016;9:186. doi: 10.3390/ma9030186. [PMC free article] [PubMed] [Cross Ref]

31. Meng F, Huang X, Jia B. Bi-directional evolutionary optimization for photonic band gap structures. J. Comput. Phys. 2015;302:393–404. doi: 10.1016/j.jcp.2015.09.010. [Cross Ref]

32. Meng F, Li S, Lin H, Jia B, Huang X. Topology optimization of photonic structures for all-angle negative refraction. Finite Elem. Anal. Des. 2016;117–118:46–56. doi: 10.1016/j.finel.2016.04.005. [Cross Ref]

33. He Z, Cai F, Ding Y, Liu Z. Subwavelength imaging of acoustic waves by a canalization mechanism in a two-dimensional phononic crystal. Appl. Phys. Lett. 2008;93:233503. doi: 10.1063/1.3043684. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |