lopscience

¡opscience.iop.org

Home Search Collections Journals About Contact us My IOPscience Acoustic coherent perfect absorbers

This content has been downloaded from IOPscience. Please scroll down to see the full text. View the table of contents for this issue, or go to the journal homepage for more

Download details:

IP Address: 131.111.164.128

This content was downloaded on 09/08/2015 at 11:49

Please note that terms and conditions apply.

New Journal of Physics

The open access journal at the forefront of physics

Deutsche PhysikalischeGeseUschaft DPG IOP Institute Of PhySjCS

Acoustic coherent perfect absorbers

J Z Song, P Bai, Z H Hang and Yun Lai

Department of Physics, Soochow University, 1 Shizi Street, Suzhou 215006, People's Republic of China

E-mail: laiyun@suda.edu.cn

Received 15 December 2013, revised 10 February 2014 Accepted for publication 13 February 2014 Published 21 March 2014

New Journal of Physics 16 (2014) 033026

doi:10.1088/1367-2630/16/3/033026

Abstract

In this paper, we explore the possibility of achieving acoustic coherent perfect absorbers. Through numerical simulations in two dimensions, we demonstrate that the energy of coherent acoustic waves can be totally absorbed by a fluid absorber with specific complex mass density or bulk modulus. The robustness of such absorbing systems is investigated under small perturbations of the absorber parameters. We find that when the resonance order is the lowest and the size of the absorber is comparable to the wavelength in the background, the phenomenon of perfect absorption is most stable. When the wavelength inside both the background and the absorber is much larger than the size of the absorber, perfect absorption is possible when the mass density of the absorber approaches the negative value of the background mass density. Finally, we show that by using suitable dispersive acoustic metamaterials, broadband acoustic perfect absorption may be achieved.

Keywords: metamaterial, acoustic metamaterial, coherent perfect absorber, acoustic absorption

Recently, acoustic metamaterials [1-36], which can exhibit almost any value of mass density and modulus, have provided a new approach to manipulating acoustic waves. Similar to their electromagnetic counterparts, acoustic metamaterials can exhibit a rich variety of new material types, such as negative mass density [1, 2], negative modulus [3], double negativity [4-13], controlled anisotropy [14, 15], etc. The origins of such special properties are attributed to the

Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal

citation and DOI.

New Journal of Physics 16 (2014) 033026 1367-2630/14/033026+13$33.00

©2014 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft

local resonance [1-10], or coiled space [11-13], or anisotropic design in metamaterials [14, 15]. Based on these unique parameters, unprecedented wave phenomena such as low-frequency sound blocking [1, 2], acoustic negative refraction and lensing [7-22], and acoustic cloaking [23-34] have been proposed. Recently, the concept of acoustic metamaterials has also been further applied to enhance the absorption of sound energy [35, 36]. In [35], a type of 'dark' acoustic metamaterial composed of a resonant membrane structure has been designed, which can absorb sound waves of wavelength much larger than its thickness. In [36], broadband large absorption has been achieved by designing porous lamella-crystals. It has been shown that acoustic metamaterials have great potential in the applications of sound absorption.

Recently, in the field of optics, a new concept named optical coherent perfect absorbers (CPAs) [37] which can achieve perfect absorption (100%) of light have attracted a lot of research interest [37-46]. Optical CPAs may be regarded as the reverse process of laser generation. Besides the function of achieving total absorption, they also demonstrate a unique way of controlling light absorption by interference effects. Optical CPA imposes strict conditions on the material and loss parameters of the absorber. For most optical materials, the permeability is unity. Therefore, both the real parts and the imaginary parts of the permittivity have to be adjusted to specific values to realize the phenomenon of optical CPA.

In this paper, we investigate the possibility and conditions of achieving acoustic CPAs, which can totally absorb sound energy instead of light. Unlike optical materials, both of the two parameters of the acoustic materials, i.e. the mass density and modulus, can exhibit many different values. In particular, acoustic metamaterials have provided an approach to engineering the effective mass density and modulus in an almost arbitrary manner. Here, through numerical analysis and simulations in two dimensions, we show that acoustic CPA can be achieved by designing the mass density or bulk modulus of the absorber to be specific complex values. CPAs can be classified by angular momentum m and resonance order n. We also investigated the robustness of such absorbing systems under small perturbations of the absorber parameters. Among different resonance orders, we find that the CPA corresponding to the lowest resonance order is the most stable. When the size of the absorber is much larger or much smaller than the wavelength in the background, the CPA effect turns out to be quite sensitive to the ratio change of the real parts of the absorber parameters. If the wavelengths both outside and inside the absorber are much larger than the size of the absorber, coherent perfect absorption is possible when the effective mass density of the absorber approaches the negative value of the background mass density. Finally, we investigate the dependence of the CPA parameters as a function of frequency. By using suitable dispersive acoustic metamaterials that match the CPA parameter dispersion, broadband acoustic perfect absorption may be achieved.

The acoustic wave equation in fluids is written as:

' 1 N 2

+ ^p = 0, (1)

where p is the pressure field, p is the mass density, k is the bulk modulus. The time average of the divergence of energy flux I = p~v, where V is velocity, gives

£ V • I dt = 1 ^Re (-imp* |Vf) + Re ^—^ |p|2 j j. If the medium is absorptive, then

I V • I dt < 0. Therefore, it is easy to see that Im (p) > 0 and Im (k) < 0 indicate energy

absorption. Therefore, the absorption of sound energy shall be attributed to the imaginary parts of mass density and bulk modulus in theory.

We consider the case of cylindrical absorbers in two dimensions. A cylindrical absorber with radius R, mass density p and bulk modulus k1 is embedded in a background fluid of mass density p0 and bulk modulus k0. In two dimensions, the pressure wave can be expanded into Hankel and Bessel functions as:

^inside

m=-œ (2)

= 2 CmJm (V ) e

2 KHfOkor) + ^(kor)]

in which am, /} and cm are the coefficients of Hankel and Bessel functions for angular momentum quantum number m. The wave vectors in the cylinders and background are k1 = m/Cj = rn^p^K-t and k0 = m/c0 = m^p0/K0, respectively. By substituting (2) into the

boundary conditions, i.e. pmside\r=R = Pouside lr=R and ^

, we obtain:

CmJm (kR) = [amHm2)(koR) + ¡iH^ (kR)

r n' (

Cmk-Po J'm (k-R) = p- [amkH '£> (k R) + Pnkfi (k R ) ]

and the scattering coefficient sm is obtained as

s = l = cpJm (kR)Hf(k0R) - cpJ!m (kR)Hm2)(kR) (4)

Sm «m C0P0 J'm (kR)H™ (kR) - cpJm (kR)H(kR) '

When the energy of the incident acoustic waves with angular momentum m is perfectly absorbed, the scattered waves should vanish, i.e. sm = 0. This is the condition for acoustic coherent perfect absorbers. However, sm = 0 is difficult to solve analytically. In this paper, we mainly rely on numerical analysis and simulations to solve the problem.

Comparing with most of the previous optical CPAs in which the loss parameter is the imaginary part of the permittivity, here, the loss parameter can be the imaginary parts of either the mass density or bulk modulus or both. Here we consider the following simple cases: (1) Im (p) > 0 and Im (k ) = 0; (2) Im (p) = 0 and Im (k ) < 0. For acoustic metamaterials, these cases correspond to absorptive meta-atoms with dipolar or monopolar resonances, respectively.

We first consider case 1 of Im (p) > 0 and Im (k) = 0. The background medium is set to be water (k0 = 2.24 x 109 Pa, p0 = 998 kg m-3) and the cylinder radius R is 0.037 m. The bulk

modulus of the cylinder k1 is set to be 3 x 109 Pa. The frequency is chosen to be 2000 Hz and the wavelength in the background is X0 = 0.75 m » R. Then, we can numerically obtain the complex value of p1 according to sm = 0 for different angular momentum m. In figure 1, we show the field maps of four cases of coherent perfect absorption corresponding to m = 0, 1, 2, 3, respectively. For m = 1, 2 , 3, a spiral wave front is produced due to the lack of reflection from

Figure 1. (a)-(d) The field maps of total coherent perfect absorption for cases of m = 0, 1, 2, 3. The background is water and the incident wave frequency is 2000 Hz. The radius and the bulk modulus of the absorber are 0.037 m and 3 X 109 Pa, respectively. The complex mass densities of the absorber are obtained as (a) 77 830.2 + 2728.28i-kgm-3, (b) -1133.23+170.56i kg m-3, (c) -1030.92+ 1.91i kg m-3, (d) -1009.38+ 0.008i kg m-3.

the absorber placed at the center. The corresponding p are found to be 77 830.2 + 2 728.28i-

kgm-3, -1 133.23+ 170.56i kg m-3, -1 030.92+1.91i kg m-3, -1 009.38+ 0.008i kg m-3, respectively. In this case, one interesting observation is that for Re (p) > 0 for m = 0 and Re (p) < 0 for m = 1, 2, 3. The negative real parts of p for m # 0 indicate that absorption is actually induced by surface-plasmon-like surface acoustic waves [20]. It can be realized by using acoustic metamaterials with negative mass density [1, 2].

Next, we consider case 2 of Im (p) = 0 and Im (k) < 0. The background medium is still water but the cylinder radius R is increased to be 0.371 m. The mass density of the cylinder p is set to be 1300 kg m-3. The frequency is still 2000 Hz. Thus, the wavelength in the background X0 = 0.75m ~ 2R. Then, we obtain the complex value of k1 according to sm = 0. In figure 2, we show the field maps of four cases of coherent perfect absorption for m = 1. It can be seen that they correspond to different resonance orders n = 1, 2, 3, 4, which may be defined as the number of amplitude peaks from the center to the surface of the cylinder. The corresponding k1 are

Figure 2. (a)-(d) The field maps of total coherent perfect absorption for cases of resonance order n = 1, 2, 3, 4. Angular momentum quantum number m = 1. The background is water and the incident wave frequency is 2000 Hz. The radius and the

mass density of the absorber are R = 0.371m and 1300 kg m 3, respectively. The complex bulk moduli of the absorber are obtained as (a) 2.1-1.2i (10 Pa), (b) 0.82

-0.24i (109 Pa), (c) 0.37-0.042i (109Pa), (d) 0.2-0.012i (109Pa).

found to be 2.1-1.2i (109Pa), 0.82-0.24i (109Pa), 0.37-0.042i (109Pa), 0.2-0.012i (109Pa), respectively. We note that the imaginary parts are all negative, indicating absorption. Similarly, for angular momentum m = 2, we also observe coherent perfect absorption corresponding to different resonance orders, which are shown in figure 3.

Figures 1, 2 and 3 show that acoustic CPA can be achieved under different angular momentum m and resonance order n by using absorbers with specific complex p or complex k1. Such complex p1 or complex k1 may be realized by using composite acoustic materials or acoustic metamaterials. However, in practice, there are always perturbations and the absorber parameters may deviate from the required precise value. Therefore, an important question is: how robust is the coherent perfect absorption system? To answer this question, we calculate the variation of |sm| within a small perturbation range of the specific complex px or complex k1 of CPA (sm = 0) and show it in figure 4. Generally, it is seen that |sm| is more sensitive to perturbations in the real parts than in the imaginary parts. In figures 4(e), (f), the sensitivity of

Figure 3. (a)-(d) The field maps of total coherent perfect absorption for cases of resonance order n = 1, 2, 3, 4. Angular momentum quantum number m = 2. The background is water and the incident wave frequency is 2000 Hz. The radius and the mass density of the absorber are R = 0.371m and 1300 kg m-3, respectively. The complex bulk moduli of the absorber are obtained as (a) 1.4-0.67i (109 Pa), (b) 0.58 -0.094i (109 Pa), (c) 0.28-0.02i (109Pa), (d) 0.16-0.006i (109Pa).

\sm\ on Re (Kj) is observed to increase when the resonance order n increases. Therefore, although there are infinite resonance orders that can achieve CPA, the lowest one, i.e. n = 1, is the most stable one.

In the following, we investigate how the complex pj or complex Kj changes as a function of frequency or the other absorber parameter. The frequency dependence is significantly important because if we can engineer the frequency dispersion of the absorber material to exhibit the CPA dispersion, then we can realize broadband CPA in principle. From (4), we can see that sm is dependent on k0R instead of m alone. In figures 5(a)-(c), we demonstrate the dependence of complex k1 as a function of kR for m = 0, 1, 2, respectively. The background medium is water. The mass density of the absorber is set as 1300 kg m-3. The resonance order n is chosen to be 1. It is clearly seen that when k0R increases, bulk modulus Kj also increases in both the real and imaginary parts (negatively). This is reasonable because the wave vector inside the absorber

Figure 4. (a)-(d) and (e)-(h) show, respectively, the absolute value of the scattering coefficient, i.e. |sm|, as a function of the real and imaginary parts of complex parameters near coherent perfect absorption for the cases shown in figures 1 and 2. Generally, |sm| is found to be more sensitive to ratio changes in the real parts than in the imaginary parts. For different n, the case of n =1 is the most stable.

Figure 5. The background medium is water. The mass density of the absorber is fixed as 1300 kg m-3. (a), (b), (c) show the dependence of complex bulk modulus as a function of kR for n =1 and m = 0, 1,2, respectively. (d) shows the dependence of complex bulk modulus as a function of the mass density for n =1, m = 1, and kR = n.

kR = (k1lk0)kR = <JK0pjlKp0 • kR should not change much so as to maintain resonance order n =1. In figures 5(d), we show the dependence of complex k1 as a function of p1 when kR = n and m = 1. It is seen that k increases when p1 increases. kR = (k1lk0)kR = ^K0pjlK1p0 • kR is

also changed less to maintain n =1.

In figure 6(a)-(c), we demonstrate the dependence of complex px as a function of kR for m = 0, 1, 2, respectively. The background medium is water. The bulk modulus of the absorber is set as 3 x 109 Pa. The resonance order n is chosen to be 1. For m = 0, it is clearly seen that when kR increases, mass density px decreases, in both the real and imaginary parts. This is reasonable

because the wave vector inside the absorber kR = (k1lk0)k0R = ^K0p1/Kjp0 • kR is almost unchanged so as to maintain resonance order n = 1. However, for m >0, when kR is large, p1 also decreases when kR increases; but when kR is small, Re (px) turns negative and approaches —p0, and Im p) decreases to zero when kR goes to zero. Such behavior is totally different from the cases in figures 5(a)-(c) and 6(a). The physical origin of this change in kR dependence lies in the different behaviors of the resonant mode. It is known that when Re p) is negative and m >0, there exists surface-plasmon-like surface acoustic waves [20]. Such a resonant mode can induce a strong field on the surface of the absorber. Therefore, even when the wavelengths inside and outside the cylinder are much larger than the size of the cylinder, CPA can still be

Figure 6. The background medium is water. The bulk modulus of the absorber is fixed as 3x109Pa. (a), (b), (c) show the dependence of the complex mass density as a function of kR for n = 1 and m = 0, 1, 2, respectively. (d) shows the dependence of the complex mass density as a function of the bulk modulus for n =1, m = 1, and kR = n.

achieved. In this case, for m > 0, we can expand the Bessel functions in (4) as Jm (x)

(m - 1)!

H2)(x)

[n (k qR)

and obtain

2m + i22m (m - 1) !m!]

n(kqR) - i22m (m - 1) !m!

Obviously when k0 i 0, p1 i - p0. Figure 6(d) shows the dependence of complex p1 as a function of k1 when kR = n and m = 1. It is seen that for large k1, Re (p1) < 0. This is consistent with the above analysis because large k1 implies small kR, and when kR i 0, p1 i - p0. Such negative mass density can be achieved by using acoustic metamaterials [1, 2]. For small k1, we find Re (p1) > 0 and p1 decreases when k1 goes to zero, so that

kR = (k1/k0)kR = ^K0p1/K1p0 • kR is less changed and n = 1 is maintained.

Figures 5 and 6 show that in order to have broadband acoustic CPA, the material of the absorber must be dispersive. Such materials may be constructed using acoustic metamaterials. It is well known that locally resonant metamaterials (see, e.g. [1]) are inherently dispersive. Here, we use the Lorentz model to describe the dispersion of the complex mass density and bulk

modulus as p1(f) = p1inf -

and K1-1(f) = K~n -

f2 - fo2 + ifo

f - fo2 + fr

and try to fit

Figure 7. (a), (b) The dispersion fitting for the complex mass density and bulk modulus in the cases of figures 6(b) and 5(b) using the following Lorentz models:

& 1 1 D

p1(f) = Pl.inf - f2 _ f 2 + ifr' Kl~ (f) = K1_inf - f 2 _ f2 + ifr' reSpectively, where f f0 + lfl0 f f0 + lfl0 f0 is the resonant frequency, r0 is the damping constant, f is the plasma frequency.

Radius of the cylinder is set as 0.2 m. In (a), a good fit region is observed between

4000 Hz and 6000 Hz with p (f) = 205 _ -:-. In (b), a good fit

1 f2 _ 1.3e7 + l5214.38f

region is observed between 3500 Hz and 4500 Hz with K1_1(f) = 0.37_2e6_

f2 _ 1.3e7 + l2627.22f '

the dispersions of CPA parameters shown in figures 6(b) and 5(b), respectively, where f is the resonant frequency, r0 is the damping constant and f is the plasma frequency. In figure 7, we

can see that a good match is achieved between 4000 Hz and 6000 Hz by fitting the complex mass density in figure 7(a), and achieved between 3500 Hz and 4500 Hz by fitting the complex bulk modulus in figure 7(b). Acoustic metamaterials with such Lorentz responses should be capable of extending the phenomenon of acoustic CPA from a single frequency to a broadband regime.

Finally, we discuss the robustness dependence of acoustic CPAs for cases under various frequency regimes. In figure 8, we show the sensitivity of |im| on the complex p or complex k1 for cases of m =1, n =1, but with different frequencies of k0 R = 10 X 2n, k0 R = 2n, k0R = 2n/10, respectively. The real k1 and real p has been chosen as 3 x 109Pa and 1300

kg m-3, respectively. It is seen that |im| becomes much more sensitive on the real parts of the complex p1 or k1 in the low and high frequencies than in the middle frequency. In other words, the CPA effect is most stable when the background wavelength is comparable to the size of the absorber. Interestingly, the dependence of |im| on the imaginary parts of the complex p1 or k1 is much more robust. Even when the imaginary parts change 20% from the CPA value, very large absorption (|5m| < 0.1) is still observed.

We note that the condition of CPA can actually be regarded as the issue of impedance matching of the boundary fields of the dissipative resonant mode inside the absorber and the incident waves in the background. Here we have discussed the characteristics of some typical resonant solutions that demonstrate acoustic CPA, but not all the solutions. For instance, under

Figure 8. (a-c) and (d-f) show the absolute value of the scattering coefficient, i.e. |sm , as a function of the real and imaginary parts of complex mass density (Re (pj) and Im (p-)) or complex bulk moduli (Re (k1) and Im (k1)) for the coherent perfect absorption cases all corresponding to m =1, n =1, but at different frequencies of k0 R = 10 X 2n, k0R = 2n, k0R = 2n/10, respectively. |sm| turns more sensitive on the real parts of the complex p or k1 in both low and high frequency regimes, but its dependence on the imaginary parts of the complex p or k1 is changed much less with frequency.

some extreme conditions such as the long wavelength limit, other interesting solutions of CPA may exist, as has been shown in [35].

In summary, we have shown that acoustic CPAs can be achieved by using absorbers with specific mass density and bulk modulus. CPAs may be classified by angular momentum m and resonance order n. Under a small disturbance of the absorber parameters, this effect is most robust when the resonance order is the lowest and the size of the absorber is comparable to the wavelength in the background. When both the wavelengths inside and outside the absorber are much larger than the size of the absorber, perfect absorption requires that the mass density of the absorber be the negative value of the background mass density. The dependence of the CPA parameters as a function of frequency is also investigated. By using suitable dispersive acoustic metamaterials that match the CPA parameter dispersion, broadband acoustic perfect absorption may be achieved.

Acknowledgments

We thank H Y Chen and J Luo for helpful discussions. This work is supported by the State Key Program for Basic Research of China (no. 2012CB921501), the National Natural Science Foundation of China (no. 11374224, 11104196, 11304215), the Natural Science Foundation of Jiangsu Province (no. BK2011277, BK20130281), the Program for New Century Excellent Talents in University (NCET), and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

References

[1] Liu Z Y etal 2000 Science 289 1734

[2] Yang Z Y etal 2008 Phys. Rev. Lett. 101 204301

[3] Fang N et al 2006 Nat. Mater. 5 452

[4] Li J and Chan C T 2004 Phys. Rev. E 70 055602

[5] Wu Y, Lai Y and Zhang Z Q 2007 Phys. Rev. B 76 205313

[6] Ding Y Q etal 2007 Phys. Rev. Lett. 99 093904

[7] Wu Y, Lai Y and Zhang Z-Q 2011 Phys. Rev. Lett. 107 105506

[8] Zhou X M and Hu G K 2009 Phys. Rev. B 79 195109

[9] Lee S H etal 2010 Phys. Rev. Lett. 104 054301

[10] Lai Y etal 2011 Nat. Mater. 10 620

[11] Liang Z X and Li J 2012 Phys. Rev. Lett. 108 114301

[12] Xie Y B et al 2013 Phys. Rev. Lett. 110 175501

[13] Liang Z X et al 2013 Sci. Rep. 3 1614

[14] Torrent D and Sanchez-Dehesa J 2010 Phys. Rev. Lett. 105 174301

[15] Christensen J and de Abajo F J G 2012 Phys. Rev. Lett. 108 124301

[16] Zhang X D and Liu Z Y 2004 Appl. Phys. Lett. 85 341

[17] Yang S X etal 2004 Phys. Rev. Lett. 93 024301

[18] Feng L et al 2006 Phys. Rev. Lett. 96 014301

[19] Lu M H etal 2007 Nat. Mater. 6 744

[20] Ambati M et al 2007 Phys. Rev. B 75 195447

[21] Zhang S, Yin L L and Fang N 2009 Phys. Rev. Lett. 102 194301

[22] Li J et al 2009 Nat. Mater. 8 931

[23] Cummer S A and Schurig D 2007 New J. Phys. 9 45

[24] Chen H Y and Chan C T 2007 Appl. Phys. Lett. 91 183518

[25] Cummer S A etal 2008 Phys. Rev. Lett. 100 024301

[26] Zhang S etal 2011 Phys. Rev. Lett. 106 024301

[27] Pendry J B and Li J 2008 New J. Phys. 10 115032

[28] Popa B-I, Zigoneanu L and Cummer S A 2011 Phys. Rev. Lett. 106 253901

[29] Garcia-Chocano V M et al 2011 Appl. Phys. Lett. 99 074102

[30] Cheng Y et al 2008 Appl. Phys. Lett. 92 151913

[31] Farhat M et al 2008 Phys. Rev. Lett. 101 134501

[32] Farhat M et al 2009 Phys. Rev. Lett. 103 024301

[33] Zhu X F etal 2011 Phys. Rev. Lett. 106 014301

[34] Stenger N, Wilhelm M and Wegener M 2012 Phys. Rev. Lett. 108 014301

[35] Mei J et al 2012 Nat. Commun. 3 756

[36] Christensen J et al 2013 arXiv:1311.1933

[37] Chong Y D et al 2010 Phys. Rev. Lett. 105 053901

[38] Wan W J et al 2011 Science 331 889-92

[39] Longhi S and Valle G D 2012 Phys. Rev. A. 85 053838

[40] PuMB et al 2012 Opt. Express. 20 2246-54

[41] Noh H et al 2012 Phys. Rev. Lett. 108 186805

[42] Noh H, Popoff S M and Cao H 2013 Opt. Express. 21 17435

[43] Zhang J F, Macdonald K F and Zheludev N I 2012 Light 1 e18

[44] Kang M et al 2013 Opt. Lett. 38 3086

[45] Ramakrishnan G et al 2013 Opt. Express. 21 16784

[46] Pirruccio G et al 2013 ACS Nano 7 4810