Issue 
EPJ Applied Metamaterials
Volume 2, 2015
Advanced Metamaterials in Microwaves, Optics and Mechanics



Article Number  10  
Number of page(s)  7  
DOI  https://doi.org/10.1051/epjam/2015017  
Published online  26 January 2016 
https://doi.org/10.1051/epjam/2015017
Research Article
Subwavelength perfect acoustic absorption in membranetype metamaterials: a geometric perspective
Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
^{*} email: sheng@ust.hk
Received:
2
September
2015
Accepted:
3
December
2015
Published online: 26 January 2016
Perfect absorption of low frequency sound with subwavelength absorbers has always been a challenge, owing to the difficulty in achieving impedance matching and the inherent weak absorption coefficients of materials at low frequencies. Recently it was shown that when a membranetype resonator’s modes are hybridized through the addition of a thin airsealed cell with a back reflecting surface, perfect absorption of low frequency acoustic wave can be achieved at a particular tunable frequency. Here we use a geometric perspective, based on the fact that the membrane is very thin and therefore the displacements on both sides of the membrane must be the same, to gain a unified framework for deriving absorption upper bounds as well as for understanding the hybrid resonance and the coherent perfect absorption on the same footing. The latter is another scheme for perfect absorption based on the phase coherence of two counterpropagating waves incident upon the membranetype resonator. Experiments were carried out to verify some relations predicted by the general framework based on this geometric perspective. Excellent agreement between theory and experiment is seen.
Key words: Metamaterials / Decorated membrane resonator / Lowfrequency sounds absorption / Coherent perfect absorption / Hybrid resonance metasurface / Geometric perspective
© M. Yang et al., Published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Even in the 21st century, the absorption of low frequency noise with subwavelength absorbers still presents a challenge. This is due to the fact that linear dynamics of dissipative systems dictates the fractional dissipated power to be linearly proportional to the elastic deformation energy [1], while the latter is usually negligible in the low frequency regime. To enhance the dissipation, it is usually necessary to increase the energy density through resonance. However, radiation coupling to resonance is an alternative that can be effective in reducing dissipation. Therefore, the difficulty of subwavelength perfect absorption may be stated as the realization of resonances in open systems with high energy density, but with the high density regions being nonradiative. In previous studies, localized resonances in membranetype metamaterials [2–7], i.e., decorated membrane resonator (DMR) with tunable weights, have shown efficient sound absorption in the low frequency regime [8]. However, as we shall see below, owing to the geometric constraint, the optimum absorption for a DMR is limited to 50% of the incident energy if the wave is incident from only one side of the sample [9]. More recently, a subwavelength perfect absorber has been realized by hybridizing DMR’s two resonances through the addition of a thin gas layer sealed in a deepsubwavelength cell (with a depth that is ~1/133 of the airborne sound wavelength) [9, 10]. The hybrid resonance is known to have maximum energy density 250 times that of the incoming sound, but radiates out only a negligible amount (<0.6% of the incident energy).
In this work, we present a geometric perspective to establish a unified framework for understanding DMR’s perfect absorption behaviors. Owing to its subwavelength lateral size, it will be shown that only DMR’s surfaceaveraged normal displacement can be coupled to the radiating modes, with all the rest of the deformation energy couples only to the evanescent modes in air [2, 6, 10]. This makes the DMR’s antiresonance, at which the surfaceaveraged normal displacement component vanishes and only the local deformation dominates [2, 11], the perfect candidate for subwavelength absorption since the coupling to radiating modes is minimal in the vicinity of antiresonance. However, as the antiresonance itself is not a resonant state, it cannot couple to the incident wave and hence the absorption at antiresonance is usually negligible. The solution for this dilemma turns out to be simply to place a reflective wall closely behind the DMR, so as to hybridize DMR’s eigenmodes into a new type of hybrid resonances whose frequencies are near DMR’s antiresonance frequencies. While the new hybrid resonance resembles the normal displacement profile near DMR’s antiresonance, its maximum normal displacement can be more than an order of magnitude larger. Since most of the large deformation energy cannot be radiated, perfect absorption is the consequence [10]. With the unified framework based on our geometric perspective, we can also easily see the emergence of coherent perfect absorption (CPA) [12–17], which constitutes another scheme for perfect absorption that utilizes two phasecoherent, counterpropagating waves to be incident on the DMR.
In what follows, Section 2 presents an analysis of the special character of the DMR that arises from its subwavelength lateral dimension, and show how that leads naturally to the antiresonance. The geometric perspective and the derived constraint for DMR’s absorption are presented in Section 3. The CPA is shown to be a natural outcome of the geometric framework, verified by numerical simulations in Section 4. The multiplescatteringbased perfect absorption is analyzed in Section 5, which leads to the emergence of DMR’s hybrid resonance, presented in Section 6. Section 7 describes the experimental realization of the hybrid resonance. Comparison with theoretical predictions is shown to give excellent agreement.
2 Decompositiion of DMR’s motions and antiresonance
An important feature of the DMR is that only its surfaceaveraged normal displacement, 〈W〉, is radiatively coupled to the acoustic waves in air. Here W denotes the displacement normal to the membrane and 〈 〉 represents surface average. The remaining deformation component δW = W − 〈W〉 can only couple to evanescent waves in air and is therefore nonradiative. This is clear from the Fourier wavevectors k_{}’s that delineate the lateral spatial pattern of W. Since acoustic wavelength λ is much larger than the lateral dimension of the DMR, it follows that for the k_{} components of δW, we must have k_{} > 2π/λ. Owing to the displacement continuity and the wave equation, we have for waves in air, where k_{⊥} denotes the wavevector component normal to the membrane. It follows that the δW component couples only to the evanescent waves as its associated k_{⊥} must be imaginary. In contrast, because the k_{} components for 〈W〉 have a distribution that peaks at k_{} = 0, it can couple to the radiation modes. It will be seen in what follows that the deformation component δW is mainly responsible for subwavelength total absorption, whereas the role of 〈W〉 in hybrid resonance is to impedance match with the incident wave so that no reflection is generated.
DMR’s responses under external pressures p are described by its Green function G. For the surface averaged component we have 〈W〉=〈G〉〈p〉. Here the relevant Green function has an explicit expression based on DMR’s eigenstates W_{n} given by [6, 11],(1)where is a parameter related to the displacementweighted mass density for the DMR’s nth eigenmode at angular eigenfrequency ω_{n}, ρ is the local density, and Ω denotes DMR’s volume. We denote the viscosity coefficient of the membrane as η. The dissipation coefficient β_{n} in equation (1) is defined by(2)
It is sometimes convenient to use surface impedance to characterize DMR’s reflection and transmission. Here the overdot represents time derivative, thus Z = i/(ω〈G〉) for harmonic motions with angular frequency ω, i.e., Z and G are inversely related.
Equation (1) has a typical Lorentzian form with the real part changing sign below and above the resonant frequency ω_{n}, with a decreasing magnitude way from ω_{n}. For an incident wave frequency ω that is inbetween the two eigenmodes, both modes will be excited but in opposite phases. At a particular inbetween frequency the net surfaceaverage component can vanish, i.e., 〈W〉 ~ 0 for a weakly dissipative system (small β_{n}). At this frequency the resonator is mostly decoupled from external waves, and it is denoted the “antiresonance” frequency. In contrast, the deformation components δW of the two resonant modes do not cancel. Hence in the vicinity of the antiresonance we have the condition . It follows that if a method can be found to create a resonance near the antiresonance frequency, then the DMR can be an efficient absorber since on the one hand it can very effectively couple to the incident wave, but on the other hand it can have a deformation component of normal displacement that is much larger than the surfaceaveraged component, as necessitated by the resonance frequency being close to the antiresonance frequency. Below we show this to be indeed possible.
3 Geometric constraint on DMR’s absorption
In this section we focus on the implications that follow from the fact that the membrane is very thin and hence the air displacements on two sides of the membrane must be the same.
Consider two counterpropagating incoming waves incident from two sides, acting on the DMR with complex pressure amplitudes and . After scattering, they are converted into two outgoing waves with complex pressure amplitudes and . Here the subscript “− (+)” denotes the left (right)hand side region and the superscript “i(o)” stands for incoming (outgoing) waves. Notice that the sound’s energy flux in air, , is opposite in directions for incident and outgoing waves. Here Z_{0} denotes the air impedance. Therefore, the air’s surfaceaveraged normal displacement velocities on two sides of DMR are given by(3a) (3b)
Since the thickness of the membrane is negligible, we must have , i.e., there is no relative motion between the two sides of the membrane. That immediately implies that the mean complex pressure amplitude (MCPA), , is conserved before and after the scattering:(4)
Conservation of MCPA implies that a portion of the incoming energy has to be preserved from dissipation. This can be made clear by direct analogy to the two (equalmass) particles scattering in classical mechanics. If we regard the p’s as momenta of the two particles before and after the collision, then equation (4) expresses the fact that the center of mass momentum is always conserved before and after the collision. For wave incident from oneside, like the collision between a moving particle and a stationary one, there is half of the incoming energy associated with the MCPA (the center of mass energy), which has to be conserved and hence cannot be dissipated. Therefore the absorption coefficient A cannot exceed A = 50%.
4 Coherent perfect absorption (CPA)
From the above perspective, a natural way to achieve perfect absorption in DMR is to make , i.e., the CPA scenario (see Figure 1c). For perfect absorption, , so that the net pressure applied on DMR is given by . Since the surfaceaveraged normal velocity must be equal to that of air, , thus . In other words, the CPA requirement boils down to . Since Z_{0} is a real number, hence in accordance with equation (1), to satisfy the CPA requirement a DMR must be in resonance so that its impedance Z = i/(ω〈G〉) is real.
Figure 1.
CPA realized by DMR. (a) The reflection (black), transmission (blue), and absorption (red) coefficients for the DMR under oneside incident. (b) The scattering (black) and absorption (red) coefficient under the CPA scenario, i.e., when two beams of sounds come from both sides with the same amplitude but opposite phase as shown in (c). The solid lines are from theoretical calculation based on the eigenmode in (d) and the open circles from numerical simulations. (c) The schematic illustration of the CPA scenario. (d) The simulated crosssectional profile of the DMR’s 1st resonance that satisfies the CPA requirement. 
To demonstrate the CPA, we simulate a DMR mounted on an aluminum frame with an area of 60.8 cm^{2} (Figure 1c). The DMR consists of a circular rubber membrane that is 19 mm in radius and 0.2 mm in thickness. The edge of the membrane is fixed on the aluminum plate, and a 371 mg circular steel platelet with 5 mm radius is attached to the center of the membrane. Numerical simulations were carried out by using COMSOL Multiphysics (with material parameters shown in Ref. [11]). The first eigenmode W_{1} of this DMR is found to be located at ω_{1} = 2π × 212.3 Hz as shown in Figure 1d. With η = 0.042, the impedance of this DMR at ω_{1} is given by equation (1) as Z = 2.004 Z_{0}. This value satisfies the CPA requirement very well.
Since the CPA requirement coincides with the 50% maximum absorption condition for onesided incident wave [16, 18], we first introduce sound wave from one side with normalized intensity, and record the associated reflection R (black) and transmission T (blue) as shown in Figure 1a. An almost 50% absorption peak for A = 1 − R^{2} − T^{2} is confirmed at the first resonance frequency (denoted by blue dashdotted line). We then introduce another sound incoming from the opposite direction with the same intensity but opposite phase (Figure 1c). As shown in Figure 1b (black), for , the two scattered waves satisfy , so that . At the resonance frequency, the relevant absorption coefficient (red) is seen to reach almost perfect absorption (~0.999915). The solid curves in Figures 1a and 1b are predictions from the impedance for each frequency based on the eigenmode shown in Figure 1d, while the open circles are results of the direct numerical simulation for the scattering process. Excellent agreement is seen.
Although CPA is an efficient scheme for perfect absorption, the required counterpropagating control wave is inconvenient in practice. In what follows, we will show a different strategy for realizing perfect absorption without the control wave.
5 Total absorption through multiple scatterings
Although the MCPA energy in the incoming waves has to be conserved in a single scattering event, its complete dissipation is still potentially possible through multiple scatterings. The outgoing MCPA energy after each scattering can serve as the incoming energy for the subsequent scatterings. In the following, we use the geometric perspective to analyze the multiplescattering process, with the aim of deriving the perfect absorption condition.
In the air layer between the membrane and the reflecting surface, the outgoing wave after mth scattering from the membrane, , is reflected by the wall, and becomes the incoming wave , for the (m + 1)th scattering. These multiply scattered waves can be superposed to form two counterpropagating waves: and . Because the reflecting wall presents a velocity node, these two wave components in front of the reflecting wall must be equal in magnitude and differ in phase by 2δ = 2k_{0} s at the position of membrane so as to cancel exactly, i.e.,
Here k_{0} denotes the sound’s wavevector in air. As the MCPA conservation is valid for each scattering, it also holds for their superpositions. Therefore, , where is the superposition for outgoing waves on the incident side after each scattering. By assuming that the incoming energy is completely absorbed by the DMR, there will be no reflection on the incident side, and . Thus the conservation of MCPA, , may be represented by a phasordiagram, shown in Figure 2b, which is an isosceles triangle with and being the two equal sides and being the base. It is clear from the diagram that , which is the backside side pressure on the DMR, is perpendicular to and hence differs from it by a phase of π/2. Therefore . Since the surfaceaveraged velocity is given by , the total pressure applied to the DMR, , is given by . It follows that the impedance condition of the DMR for achieving total absorption is given by(5)
Figure 2.
Schematic illustration of the multiplescattered waves that can lead to perfect absorption on a DMR, with a reflecting wall placed a small distance s behind it. (a) When the incoming wave energy has been totally absorbed, there should be no reflections, i.e., all the scatterings on the incident side cancel each other so that . Incoming wave generates only two net counterpropagating wave components and , in the air layer between the DMR (at z = 0) and the reflecting wall. Here the dotdash lines stand for the location at which the periodic boundary condition is applied. (b) Due to the MCPA conservation, the three complex pressure amplitudes must form an isosceles triangle whose two equal sides are slanted toward each other with an angle given by 2δ = 2k_{0}s. 
For δ = π/2, we have Z = Z_{0}, i.e., impedance matching to air. This occurs at cell’s “drum” resonances when s = λ/4. However, for an air layer thinner than a quarterwavelength, δ → 0, and the required imaginary part of Z approaches −i∞. Since the impedance of the whole structure is the addition of the two serial impedances of the DMR and the sealed cell behind it, perfect absorption requires the net impedance of the composite structure Z_{c} to satisfy the condition Z_{c} = Z_{0}, i.e., impedance matched to air. In this context the imaginary part of Z in equation (5) is noted to just cancels the large (positive) imaginary part of the thin air layer, so that Im(Z_{c}) = 0. But can the condition as stated by equation (5) be satisfied by the DMR, especially the large negative imaginary part of the impedance? It turns out that the DMR’s imaginary part of impedance diverges in the vicinity of antiresonance frequency, as shown below. Hence impedance matching of our composite structure can indeed occur, but only at a frequency close to the antiresonance condition. Also, associated with Im(Z_{c}) = 0 is that the Green function of the composite structure, 〈G_{c}〉 = i/(ωZ_{c}), has a peak in its imaginary part. This is clear sinceso that Im〈G_{c}〉 reaches a peak when Im(Z_{c}) = 0. A peak in the imaginary part of the Green function signifies a resonance mode.
In what follows, we show that a new hybrid resonance of the composite structure precisely fits the required profile deduced from our analysis.
6 Hybrid resonance and total absorption
In order to see the emergence of the hybrid mode, for simplicity we consider the normal displacements of only two relevant eigenmodes, W_{1} and W_{2} (see Figure 3b). The surfaceaveraged Green function of the DMR in equation (1) can be written as(6)where β, the averaged coefficient of β_{1} and β_{2}, is taken to be a small quantity so that the dimensionless β/ω ≪ 1. In anticipation that the hybrid resonance occurs near the antiresonance, we expand equation (6) in the vicinity of DMR’s antiresonance frequency , at which Re(〈G〉) = 0. The expansion can be written as 〈G〉 ≃ 2Ξ(iβ − Δω), where
Figure 3.
(a) Schematic illustration of the unit cell's component and geometry of the hybrid resonator. (b) Schematic crosssectional illustration of the two lowest frequency eigenmodes of the DMR, plotted in reference to the same phase of incident wave, with W being the normal displacement of membrane, normalized to its surface average component 〈W〉. (c) Displacement profile of the DMR’s hybrid mode at 152 Hz, which is clearly a superposition of the two in (b), but with an almost an order of magnitude larger amplitude. However, such large amplitude oscillations are only coupled to the nonradiative evanescent waves. Experimentally measured results (symmetrically plotted from data along the radial direction on one side) are shown in red circles while the simulated results are in black curves. Excellent agreement is seen. 
From the inverse relationship between the Green function and impedance, the associated Z is given by(7)
From equation (7) it is clear that if β = 0, then the imaginary part of the impedance diverges negatively when Δω approaches zero from the positive side. When β is finite but small, then the imaginary part of the impedance can be as large in magnitude as , at a frequency below such that Δω = β. Therefore, perfect absorption can indeed occur for a DMR backed by a reflecting surface. It is noted that the perfect absorption actually occurs under the condition of a small dissipative coefficient β.
Equating equations (5) and (7) yields two equations for the condition of total absorption:(8a) (8b)
Although the dissipation coefficient β is not easily adjustable, this total absorption condition is still robust since we have two other parameters, s and Δω, which can be easily tuned experimentally to satisfy equation (8). Notice that the solution of equation (8) occurs close to the antiresonance (Δω → 0) (Figure 4a). As the hybrid mode frequency is between the two DMR modes, its displacement profile is necessarily the linear superposition of the two DMR eigenmodes. It therefore resembles the profile of membrane motions at the vicinity of the antiresonance frequency. In particular, the magnitudes of the two DMR modes, 〈W_{1}〉 and 〈W_{2}〉, must be comparable here – namely, 〈W_{1}〉/〈W_{2}〉 ≃ 1. However, impedance matching implies 〈W_{1}〉 − 〈W_{2}〉 ≃ W_{s}, the air displacement of incoming sounds. These two conditions can be satisfied only when 〈W_{1,2}〉 ≫ W_{s}, so that 1 − 〈W_{2}〉/〈W_{1}〉 = W_{s}/〈W_{1}〉 is a small number. In other words, the hybrid resonance satisfies the condition , owing to its proximity to the antiresonance frequency.
Figure 4.
Relationships between different parameters at the hybrid resonant frequency with perfect absorption. (a) Dimensionless dissipation parameter (black curve) and dimensionless thickness (red curve) of the sealed air layer, plotted as functions of the dimensionless relative frequency location of the hybrid resonance. (b) Deformation component of the normal displacement at the hybrid resonance frequency, plotted as a function of . Impedance matching means that the surfaceaveraged displacement has to be equal to that of the incident sound – that is, 〈W〉 = 〈W_{1}〉 − 〈W_{2}〉 = W_{s}. Proximity to the antiresonance also means that 〈W_{1}〉 ≃ 〈W_{2}〉. These two conditions imply a large W compared with its surface average, and account for its large variance component of the deformation . In particular, the small dissipation coefficient and the subwavelength scale dovetail with the large W in a selfconsistent manner so as to achieve total absorption. The dotdashed blue lines in the two figures indicate the parameter values relevant to the experimental results shown in Figure 5. 
7 Experimental realization
Experiments reported in reference [10] have corroborated the existence of the hybrid resonance with perfect absorption capability. By placing an aluminumreflecting wall behind a DMR with a distance s (Figure 3a), an absorption peak with magnitude greater than 99% was observed at 152 Hz (Figure 5c). In this case the DMR comprises a 90 mmwide circular elastic membrane decorated by a 800 mg platelet with a radius of 10 mm. The value of s is 21.7 mm, which is inferred from the experimental data of 17 mm of SF6 gas, as the effects of the two thin gas layers are identical if their thicknesses are scaled linearly with their adiabatic indices (SF6 has an adiabatic index of 1.098, compared to 1.4 for air). Note here that the air layer is in the deepsubwavelength regime, s ≃ λ/104, where λ = 2.26 m is the wavelength of 152 Hz sound in air.
Figure 5.
(a), (b) The real and imaginary parts of the DMR’s normalized surface impedance Z/Z_{0}, where Z_{0} denotes the impedance of air, plotted as a function of frequency. They are seen to satisfy the total absorption requirements given by equation (5) at the hybrid resonance frequency, indicated by the intersection of the black curves with the blue dotdashed line. Total absorption is confirmed experimentally, shown in (c). A sharp absorption peak, reaching 0.994, is seen. The solid and dashed arrows indicate the first eigenmode and the antiresonance frequency of the DMR, respectively. Experimental results are shown by the red circles, and the black curves denote the theoretical predictions. 
Numerical simulations by COMSOL Multiphysics (with material parameters shown in Ref. [11]) are based on the two relevant eigenmodes for the DMR as shown in Figure 3b. The values of the relevant parameters are Ξ = 1.268 × 10^{−9} m^{2} s^{3}/kg and Hz. In Figure 4a we plot the dimensionless (black curve) and (red curve) that satisfy condition (8) as a function of . Here γ = 1.4 is the adiabatic index of air and p_{0} is the atmospheric pressure. It is clear that for the hybrid resonant frequency close to antiresonance, we have , and the associated dissipation coefficient and air layer thickness are both much smaller than 1. In particular, for the value , which corresponds with the observed total absorption at 2π × 152 Hz (indicated by blue dotdashed line), one can evaluate the dissipation coefficient (β = 5.14 Hz), and predict the thickness of air layer to be (s = 21.9 mm). The latter is noted to be in excellent agreement with the experimental value of s = 21.7 mm.
The real and imaginary parts of surface impedance of the DMR, given by equation (7) as a function of frequency, are shown in Figures 5a and 5b, respectively. At the total absorption frequency, we have Z = 1.001 – 12.961i, which is noted to match the value required by equation (5), Z_{0}[1 – icot(k_{0}s)] = 1 – 16.342i, reasonably well. Furthermore, from the impedance for the composite system comprising the DMR and the air layer, one can calculate the reflection as R = (Z_{0} − Z_{c})/(Z_{0} + Z_{c}). The absorption coefficient can then be evaluated as 1 − R^{2}, which is plotted as the black curve in Figure 5c. Excellent agreement with the experiment (red circles) is seen.
In Figure 4b we have plotted the calculated magnitude of the deformation component as a function of. In agreement with our expectation, the magnitude of the deformation component is seen to diverge at , at which point both the required s and β approach zero. To confirm the large deformation component in this hybrid resonance, we have used laser vibrometer to measure the vibrational profile of the DMR at the total absorption frequency. The experimental result, normalized to 1 Pa (corresponding to 94 dB) of incident sound amplitude, is shown in Figure 3c as red circles. This hybrid resonance profile is clearly a linear superposition of the membrane’s two eigenmodes (Figure 3b), and resembles the motion profile at DMR’s antiresonance [2]. In the same figure are shown the numerical simulation results (black curve), which are in excellent agreement with the measured profile. From the measured displacement, we can evaluate 〈W〉 ≃ 2.8 µm, which is close to the incident sound amplitude of W_{s} = 2.5 µm – a clear evidence of impedance matching. But 〈W〉 is noted to be one order of magnitude smaller than the maximum displacement of the platelet (over 40 µm). We can also deduce from the measured data that µm, which is consistent with the relation shown in Figure 4b. This significantly amplified displacement component near the antiresonance frequency explains the total absorption functionality for the weakly lossy membrane in subwavelength.
8 Concluding remarks
In this work we use a framework that is based on the geometric perspective to derive upper bounds on absorption by thin membranetype resonators. It is seen that the CPA and hybrid resonance constitute two scenarios in which perfect absorption can be achieved. In particular, the hybrid resonance utilizes multiple scatterings generated by a back reflecting surface to create a new resonance whose frequency is in the vicinity of the antiresonance. The hybrid resonance’s very large deformation component of the membrane, which is nonradiative, is seen to couple with impedance matching to yield perfect absorption. All theoretical predictions are shown to be in excellent agreement with experimental observations.
Acknowledgments
This work is supported by Hong Kong RGC Grants AoE/P02/12.
References
 L. Landau, E. Lifshitz, Theory of elasticity, Pergamon Press, 1970, Chap. 5, p. 153. [Google Scholar]
 Z. Yang, J. Mei, M. Yang, N. Chan, P. Sheng, Phys. Rev. Lett. 101 (2008) 204301. [CrossRef] [PubMed] [Google Scholar]
 S.H. Lee, C.M. Park, Y.M. Seo, C.K. Kim, Phys. Rev. B 81 (2010) 241102. [CrossRef] [Google Scholar]
 C.M. Park, J.J. Park, S.H. Lee, Y.M. Seo, C.K. Kim, S.H. Lee, Phys. Rev. Lett. 107 (2011) 194301. [CrossRef] [Google Scholar]
 J.J. Park, K. Lee, O.B. Wright, M.K. Jung, S.H. Lee, Phys. Rev. Lett. 110 (2013) 244302. [CrossRef] [Google Scholar]
 M. Yang, G. Ma, Z. Yang, P. Sheng, Phys. Rev. Lett. 110 (2013) 134301. [CrossRef] [Google Scholar]
 G. Ma, M. Yang, Z. Yang, P. Sheng, Appl. Phys. Lett. 103 (2013) 011903. [CrossRef] [Google Scholar]
 J. Mei, G. Ma, M. Yang, Z. Yang, W. Wen, P. Sheng, Nat. Commun. 3 (2012) 756. [CrossRef] [PubMed] [Google Scholar]
 M. Yang, Y. Li, C. Meng, C. Fu, J. Mei, Z. Yang, P. Sheng, C. R. Mecanique 343 (2015) 635–644. [CrossRef] [Google Scholar]
 G. Ma, M. Yang, S. Xiao, Z. Yang, P. Sheng, Nat. Mater. 13 (2014) 873. [CrossRef] [Google Scholar]
 M. Yang, G. Ma, Y. Wu, Z. Yang, P. Sheng, Phys. Rev. B 89 (2014) 064309. [CrossRef] [Google Scholar]
 Y. Chong, L. Ge, H. Cao, A.D. Stone, Phys. Rev. Lett. 105 (2010) 053901. [CrossRef] [PubMed] [Google Scholar]
 W. Wan, Y. Chong, L. Ge, H. Noh, A.D. Stone, H. Cao, Science 331 (2011) 889. [CrossRef] [PubMed] [Google Scholar]
 S. Li, J. Luo, S. Anwar, S. Li, W. Lu, Z.H. Hang, Y. Lai, B. Hou, M. Shen, C. Wang, Sci. Rep. 4 (2014) 7369. [CrossRef] [Google Scholar]
 J.R. Piper, S. Fan, ACS Photonics 1 (2014) 347. [CrossRef] [Google Scholar]
 P. Wei, C. Croënne, S.T. Chu, J. Li, Appl. Phys. Lett. 104 (2014) 121902. [CrossRef] [Google Scholar]
 X. Cai, Q. Guo, G. Hu, J. Yang, Appl. Phys. Lett. 105 (2014) 121901. [CrossRef] [Google Scholar]
 M. Yang, C. Meng, C. Fu, Y. Li, Z. Yang, P. Sheng, Appl. Phys. Lett. 107 (2015) 104104. [CrossRef] [Google Scholar]
Cite this article as: Yang M, Ma G, Yang Z & Sheng P: Subwavelength perfect acoustic absorption in membranetype metamaterials: a geometric perspective. EPJ Appl. Metamat. 2015, 2, 10.
All Figures
Figure 1.
CPA realized by DMR. (a) The reflection (black), transmission (blue), and absorption (red) coefficients for the DMR under oneside incident. (b) The scattering (black) and absorption (red) coefficient under the CPA scenario, i.e., when two beams of sounds come from both sides with the same amplitude but opposite phase as shown in (c). The solid lines are from theoretical calculation based on the eigenmode in (d) and the open circles from numerical simulations. (c) The schematic illustration of the CPA scenario. (d) The simulated crosssectional profile of the DMR’s 1st resonance that satisfies the CPA requirement. 

In the text 
Figure 2.
Schematic illustration of the multiplescattered waves that can lead to perfect absorption on a DMR, with a reflecting wall placed a small distance s behind it. (a) When the incoming wave energy has been totally absorbed, there should be no reflections, i.e., all the scatterings on the incident side cancel each other so that . Incoming wave generates only two net counterpropagating wave components and , in the air layer between the DMR (at z = 0) and the reflecting wall. Here the dotdash lines stand for the location at which the periodic boundary condition is applied. (b) Due to the MCPA conservation, the three complex pressure amplitudes must form an isosceles triangle whose two equal sides are slanted toward each other with an angle given by 2δ = 2k_{0}s. 

In the text 
Figure 3.
(a) Schematic illustration of the unit cell's component and geometry of the hybrid resonator. (b) Schematic crosssectional illustration of the two lowest frequency eigenmodes of the DMR, plotted in reference to the same phase of incident wave, with W being the normal displacement of membrane, normalized to its surface average component 〈W〉. (c) Displacement profile of the DMR’s hybrid mode at 152 Hz, which is clearly a superposition of the two in (b), but with an almost an order of magnitude larger amplitude. However, such large amplitude oscillations are only coupled to the nonradiative evanescent waves. Experimentally measured results (symmetrically plotted from data along the radial direction on one side) are shown in red circles while the simulated results are in black curves. Excellent agreement is seen. 

In the text 
Figure 4.
Relationships between different parameters at the hybrid resonant frequency with perfect absorption. (a) Dimensionless dissipation parameter (black curve) and dimensionless thickness (red curve) of the sealed air layer, plotted as functions of the dimensionless relative frequency location of the hybrid resonance. (b) Deformation component of the normal displacement at the hybrid resonance frequency, plotted as a function of . Impedance matching means that the surfaceaveraged displacement has to be equal to that of the incident sound – that is, 〈W〉 = 〈W_{1}〉 − 〈W_{2}〉 = W_{s}. Proximity to the antiresonance also means that 〈W_{1}〉 ≃ 〈W_{2}〉. These two conditions imply a large W compared with its surface average, and account for its large variance component of the deformation . In particular, the small dissipation coefficient and the subwavelength scale dovetail with the large W in a selfconsistent manner so as to achieve total absorption. The dotdashed blue lines in the two figures indicate the parameter values relevant to the experimental results shown in Figure 5. 

In the text 
Figure 5.
(a), (b) The real and imaginary parts of the DMR’s normalized surface impedance Z/Z_{0}, where Z_{0} denotes the impedance of air, plotted as a function of frequency. They are seen to satisfy the total absorption requirements given by equation (5) at the hybrid resonance frequency, indicated by the intersection of the black curves with the blue dotdashed line. Total absorption is confirmed experimentally, shown in (c). A sharp absorption peak, reaching 0.994, is seen. The solid and dashed arrows indicate the first eigenmode and the antiresonance frequency of the DMR, respectively. Experimental results are shown by the red circles, and the black curves denote the theoretical predictions. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.