Subwavelength perfect acoustic absorption in membrane-type metamaterials : a geometric perspective

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 membrane-type resonator’s modes are hybridized through the addition of a thin air-sealed 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 counter-propagating waves incident upon the membrane-type 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.


Introduction
Even in the 21st century, the absorption of low frequency noise with sub-wavelength 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 non-radiative.In previous studies, localized resonances in membrane-type metamaterials [2][3][4][5][6][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 deep-subwavelength 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 surface-averaged 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 anti-resonance, at which the surface-averaged 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 anti-resonance itself is not a resonant state, it cannot couple to the incident wave and hence the absorption at anti-resonance 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 anti-resonance frequencies.While the new hybrid resonance resembles the normal displacement profile near DMR's anti-resonance, 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][13][14][15][16][17], which constitutes another scheme for perfect absorption that utilizes two phase-coherent, counter-propagating 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 anti-resonance.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 multiple-scattering-based 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.

Decompositiion of DMR's motions and anti-resonance
An important feature of the DMR is that only its surfaceaveraged normal displacement, hWi, is radiatively coupled to the acoustic waves in air.Here W denotes the displacement normal to the membrane and h i represents surface average.The remaining deformation component dW = W À hWi 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 k is much larger than the lateral dimension of the DMR, it follows that for the k || components of dW, we must have |k || | > 2p/k.Owing to the displacement continuity and the wave equation, we have jk jj j 2 þ k 2 ?¼ ð2p=kÞ 2 for waves in air, where k ?denotes the wavevector component normal to the membrane.It follows that the dW component couples only to the evanescent waves as its associated k ?must be imaginary.In contrast, because the k || components for hWi 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 dW is mainly responsible for subwavelength total absorption, whereas the role of hWi 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 hWi=hGihpi.Here the relevant Green function has an explicit expression based on DMR's eigenstates W n given by [6,11], where q n R X qjW n j 2 dV is a parameter related to the displacement-weighted mass density for the DMR's nth eigenmode at angular eigen-frequency x n , q is the local density, and X denotes DMR's volume.We denote the viscosity coefficient of the membrane as g.The dissipation coefficient b n in equation ( 1) is defined by It is sometimes convenient to use surface impedance Z hpi=h _ W i to characterize DMR's reflection and transmission.Here the over-dot represents time derivative, thus Z = i/(xhGi) for harmonic motions with angular frequency x, 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 x n , with a decreasing magnitude way from x n .For an incident wave frequency x that is in-between the two eigen-modes, both modes will be excited but in opposite phases.At a particular in-between frequency the net surface-average component can vanish, i.e., hWi ~0 for a weakly dissipative system (small b n ).At this frequency the resonator is mostly decoupled from external waves, and it is denoted the ''anti-resonance'' frequency.In contrast, the deformation components dW of the two resonant modes do not cancel.Hence in the vicinity of the anti-resonance we have the condition ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi h dW ð Þ 2 i q > hW i.It follows that if a method can be found to create a resonance near the anti-resonance 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 surface-averaged component, as necessitated by the resonance frequency being close to the anti-resonance frequency.Below we show this to be indeed possible.

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 counter-propagating incoming waves incident from two sides, acting on the DMR with complex pressure amplitudes p i À and p i þ .After scattering, they are converted into two outgoing waves with complex pressure amplitudes p o À 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 surface-averaged normal displacement velocities on two sides of DMR are given by Since the thickness of the membrane is negligible, we must have h _ , there is no relative motion between the two sides of the membrane.That immediately implies that the mean complex pressure amplitude (MCPA), " p, is conserved before and after the scattering: 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 (equal-mass) 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 one-side, 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%.

Coherent perfect absorption (CPA)
From the above perspective, a natural way to achieve perfect absorption in DMR is to make , the CPA scenario (see Figure 1c).For perfect absorption, p o À ¼ p o þ ¼ 0, so that the net pressure applied on DMR is given by Since the surface-averaged normal velocity h _ W i must be equal to that of air, h _ In other words, the CPA requirement boils down to Z ¼ hp tot i=h _ W i ¼ 2Z 0 .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/(xhGi) is real.
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 x 1 = 2p • 212.3 Hz as shown in Figure 1d.With g = 0.042, the impedance of this DMR at x 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 one-sided 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 dash-dotted 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 " p ¼ 0, the two scattered waves satisfy 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 counter-propagating control wave is inconvenient in practice.In what follows, we will show a different strategy for realizing perfect absorption without the control wave.

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 multiple-scattering 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, p om þ , is reflected by the wall, and becomes the incoming wave p i mþ1 þ , for the (m + 1)th scattering.These multiply scattered waves can be superposed to form two counter-propagating waves: 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 2d = 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 p o À ¼ P 1 m¼0 p om À 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 p o À ¼ 0. Thus the conservation of MCPA, , may be represented by a phasor-diagram, shown in Figure 2b, which is an isosceles triangle with p i þ and p o þ being the two equal sides and p i À being the base.It is clear from the diagram that which is the backside side pressure on the DMR, is perpendicular to p i À and hence differs from it by a phase of p/2.Therefore , is given by hp tot i ¼ p i À ð1 À i cot dÞ.It follows that the impedance condition of the DMR for achieving total absorption is given by For d = p/2, we have Z = Z 0 , i.e., impedance matching to air.This occurs at cell's ''drum'' resonances when s = k/4.However, for an air layer thinner than a quarter-wavelength, d !0, and the required imaginary part of Z approaches Ài1.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 anti-resonance frequency, as shown below.Hence impedance matching of our composite structure can indeed occur, but only at a frequency close to the anti-resonance condition.Also, associated with Im(Z c ) = 0 is that the Green function of the composite structure, hG c i = i/(xZ c ), has a peak in its imaginary part.This is clear since so that ImhG c i 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.

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 eigen-modes, W 1 and W 2 (see Figure 3b).The surface-averaged Green function of the DMR in equation ( 1) can be written as hGi ' where b, the averaged coefficient of b 1 and b 2 , is taken to be a small quantity so that the dimensionless b/x ( 1.In anticipation that the hybrid resonance occurs near the anti-resonance, we expand equation ( 6) in the vicinity of DMR's anti-resonance frequency x, at which Re(hGi) = 0.The expansion can be written as hGi ' 2N(ib À Dx), where From the inverse relationship between the Green function and impedance, the associated Z is given by From equation (7) it is clear that if b = 0, then the imaginary part of the impedance diverges negatively when Dx approaches zero from the positive side.When b is finite but small, then the imaginary part of the impedance can be as large in magnitude as À1=ð4 xNbÞ, at a frequency below x such that Dx = b.Therefore, perfect absorption can indeed occur for a DMR backed by a reflecting surface.It is noted that the perfect in the air layer between the DMR (at z = 0) and the reflecting wall.Here the dot-dash 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 2d = 2k 0 s.Equating equations ( 5) and ( 7) yields two equations for the condition of total absorption: Although the dissipation coefficient b is not easily adjustable, this total absorption condition is still robust since we have two other parameters, s and Dx, which can be easily tuned experimentally to satisfy equation (8).Notice that the solution of equation ( 8) occurs close to the anti-resonance (Dx !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 eigen-modes.It therefore resembles the profile of membrane motions at the vicinity of the anti-resonance frequency.In particular, the magnitudes of the two DMR modes, hW 1 i and hW 2 i, must be comparable here -namely, hW 1 i/hW 2 i ' 1.However, impedance matching implies hW 1 i À hW 2 i ' W s , the air displacement of incoming sounds.These two conditions can be satisfied only when hW 1,2 i ) W s , so that 1 À hW 2 i/hW 1 i = W s /hW 1 i is a small number.In other words, the hybrid resonance satisfies the condition ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi hðdW Þ 2 i q > hW i, owing to its proximity to the antiresonance frequency.

Experimental realization
Experiments reported in reference [10] have corroborated the existence of the hybrid resonance with perfect absorption capability.By placing an aluminum-reflecting 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 mm-wide 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 deep-subwavelength regime, s ' k/104, where k = 2.26 m is the wavelength of 152 Hz sound in air.
Numerical simulations by COMSOL Multiphysics (with material parameters shown in Ref. [11]) are based on the two relevant eigen-modes for the DMR as shown in Figure 3b.The values of the relevant parameters are N = 1.268 • 10 À9 m 2 s 3 /kg and x ¼ 2p Â 162:3 Hz.In Figure 4a we plot the dimensionless " b ¼ b= x (black curve) and " s ¼ s=ðcp 0 Z À1 0 xÀ1 Þ (red curve) that satisfy condition (8) as a function of Á " x ¼ Dx= x.Here c = 1.4 is the adiabatic index of air and p 0 is the atmospheric pressure.It is clear that x.Impedance matching means that the surfaceaveraged displacement has to be equal to that of the incident sound -that is, hWi = hW 1 i À hW 2 i = W s .Proximity to the anti-resonance x also means that hW 1 i ' hW 2 i.These two conditions imply a large W compared with its surface average, and account for its large variance component of the deformation ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi hðdW Þ 2 i q .In particular, the small dissipation coefficient and the subwavelength scale dovetail with the large W in a self-consistent manner so as to achieve total absorption.The dot-dashed blue lines in the two figures indicate the parameter values relevant to the experimental results shown in Figure 5. for the hybrid resonant frequency close to anti-resonance, we have Á " x !0, and the associated dissipation coefficient " b and air layer thickness " s are both much smaller than 1.In particular, for the value Á " x ¼ 0:06, which corresponds with the observed total absorption at 2p • 152 Hz (indicated by blue dot-dashed line), one can evaluate the dissipation coefficient " b ¼ 0:00502 (b = 5.14 Hz), and predict the thickness of air layer to be " s ¼ 0:0772 (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 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi hðdW Þ 2 i q =W s as a function of Á " x.In agreement with our expectation, the magnitude of the deformation component is seen to diverge at xðÁ " x ¼ 0Þ, at which point both the required s and b 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 anti-resonance [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 hWi ' 2.8 lm, which is close to the incident sound amplitude of W s = 2.5 lm -a clear evidence of impedance matching.But hWi is noted to be one order of magnitude smaller than the maximum displacement of the platelet (over 40 lm).We can also deduce from the measured data that ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi hðdW Þ 2 i q ¼ 5:2hW i ' 14:5 lm, which is consistent with the relation shown in Figure 4b.This significantly amplified displacement component near the anti-resonance frequency explains the total absorption functionality for the weakly lossy membrane in subwavelength.

Concluding remarks
In this work we use a framework that is based on the geometric perspective to derive upper bounds on absorption by thin membrane-type 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 anti-resonance.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.

Figure 1 .
Figure 1.CPA realized by DMR.(a) The reflection (black), transmission (blue), and absorption (red) coefficients for the DMR under one-side 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 cross-sectional profile of the DMR's 1st resonance that satisfies the CPA requirement.

Figure 2 .
Figure 2. Schematic illustration of the multiple-scattered 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 thatp o À exp Àik 0 z ð Þ¼ P 1 m¼0 p om À exp Àik 0 z ð Þ¼0.Incoming wave p i À exp Àik 0 z ð Þ generates only two net counter-propagating wave components p o þ exp Àik 0 z ð Þ¼ P 1 m¼0 p om þ exp Àik 0 z ð Þ and p i 2 exp Àik 0 z ð Þ¼ P 1 m¼1 p im þ exp Àik 0 z ð Þ,in the air layer between the DMR (at z = 0) and the reflecting wall.Here the dot-dash 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 2d = 2k 0 s.

Figure 4 .
Figure 4. Relationships between different parameters at the hybrid resonant frequency with perfect absorption.(a) Dimensionless dissipation parameter " b (black curve) and dimensionless thickness " s (red curve) of the sealed air layer, plotted as functions of the dimensionless relative frequency location Á " x of the hybrid resonance.(b) Deformation component ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi hðdW Þ 2 i q of the normal displacement at the hybrid resonance frequency, plotted as a function of Á "x.Impedance matching means that the surfaceaveraged displacement has to be equal to that of the incident sound -that is, hWi = hW 1 i À hW 2 i = W s .Proximity to the anti-resonance x also means that hW 1 i ' hW 2 i.These two conditions imply a large W compared with its surface average, and account for its large

Figure 3 .
Figure 3. (a) Schematic illustration of the unit cell's component and geometry of the hybrid resonator.(b) Schematic cross-sectional illustration of the two lowest frequency eigen-modes 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 hWi.(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 non-radiative 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.

Figure 5 .
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 dot-dashed 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 anti-resonance frequency of the DMR, respectively.Experimental results are shown by the red circles, and the black curves denote the theoretical predictions.