Abstract
Controlling acoustic fields is crucial in diverse applications such as loudspeaker design, ultrasound imaging and therapy or acoustic particle manipulation. The current approaches use fixed lenses or expensive phased arrays. Here, using a process of analoguetodigital conversion and wavelet decomposition, we develop the notion of quantal metasurfaces. The quanta here are small, premanufactured threedimensional units—which we call metamaterial bricks—each encoding a specific phase delay. These bricks can be assembled into metasurfaces to generate any diffractionlimited acoustic field. We apply this methodology to show experimental examples of acoustic focusing, steering and, after stacking single metasurfaces into layers, the more complex field of an acoustic tractor beam. We demonstrate experimentally singlesided airborne acoustic levitation using metalayers at various bitrates: from a 4bit uniform to 3bit nonuniform quantization in phase. This powerful methodology dramatically simplifies the design of acoustic devices and provides a keystep towards realizing spatial sound modulators.
Introduction
The shaping of light using spatial light modulators (SLMs) is an established technology for advanced threedimensional (3D) displays^{1} and micromanipulation^{2}. In the SLM an incident beam of coherent light is transformed via amplitude and phase manipulation into a wide range of reflected or transmitted optical distributions. Crucially the SLM is computer controlled and it is possible to reconfigure the optical field produced in almost realtime.
The acoustic equivalent of an SLM does not exist. Traditionally, the control of an acoustic field distribution was achieved by using fixed lenses^{3}, which perform a single function, or phased arrays^{4,5}, where the amplitudes and phases of the individual array elements are independently controlled. However, phased arrays are often bulky and expensive, with cost and complexity scaling linearly with the number of channels. Despite these limitations, phased arrays are in widespread use. In High Frequency Focused Ultrasound (HIFU), for example, sparse arrays of transducers are used to treat a variety of tumours^{6} or functional brain disorders^{7}, inducing a localized heating effect, even behind the ribs^{8}. In industrial applications, focusing and steering of ultrasonic waves is required to find small cracks in metallic components, which can be complex in geometry and highly anisotropic^{9}. New applications that require precise control of acoustic waves include parametric loudspeakers^{10}, ultrahaptics^{11}, caustic engineering and acoustic levitation^{12}.
Acoustic metamaterials are an emerging class of engineered materials designed to control, direct and manipulate acoustic waves^{13}. Typically, they are made up of a collection of subwavelength structures (that is, unit cells), and are characterized by their effective mass density and bulk modulus. The possibility of manufacturing metamaterials with negative effective parameters^{14} has led to effects such as negative refraction^{15} and subdiffraction focusing. Particularly interesting for beamshaping applications are twodimensional (2D) planar metasurfaces: closelypacked structures of phase shifters whose thickness is comparable to the wavelength of operation. Very recent examples of acoustic metasurfaces include the use of labyrinthine structures^{16}, helical structures^{17}, spacecoiling^{18,19}, multislits^{20,21} and Helmholtz resonators^{22,23,24,25}. These diverse metasurfaces have always been built using a limited number of unit cells, and as such the optimal number of phase discretization levels required for each application has been ignored. For instance, previous studies^{16,17,18,19,20,21,22,23,24,25,26,27,28,29} suggest a wide range of discretization levels, varying from 2 to 256, selected on an adhoc basis.
Here we develop the notion of quantal metasurfaces to demonstrate a different metamaterial concept, based on the use of a small set of premanufactured 3D unit cells, termed metamaterial bricks, which can be assembled into 2D structures ondemand. The bricks become, in isolation, the building blocks of an assembly, encoding prerequisite phase delays. This operation is a form of analoguetodigital conversion (Fig. 1): the desired acoustic pressure field is sampled at a certain distance from a metasurface and used as input for acoustic holography, leading to a phase distribution that gets quantized in the spatial and phase domains, whose values are then mapped into a series of premanufactured metamaterial bricks. Starting from a limited set of unique bricks, we use a discrete wavelet transform based method to synthesize the metasurface needed in a given application, optimizing the number of bricks needed. We obtain a reconfigurable metasurface, which transforms an incident sound wave into an arbitrary range of diffractionlimited acoustic fields. Furthermore, once a number of metasurfaces are available, each performing a given transformation (for example, steering and focusing), we show that single metasurfaces can be stacked into layers to perform more complex transformations, thus creating the acoustic equivalent of optical components. This simple, yet powerful, concept simplifies the design of acoustic devices and systems, and lays the foundations for realization of spatial sound modulators (SSMs).
Results
Metasurface quantization
A number of studies have explored how to encode a uniform phase distribution ϕ(x, y) in a metasurface to produce given acoustic transformations^{13,14}. Realizing a quantal metasurface, however, requires AD conversion with two parameters: one in the spatial domain, which—for a fixed size of metasurface—depends on the size of its unit cells and on how the phase ϕ_{ij} is assigned to each location (x_{i}, y_{j}), and one in the phase domain, which governs the number of different phases. Practically, the AD step defines the parts list and assembly instructions for a particular metasurface. We chose a spatial resolution of λ_{0}/2, which is a good compromise between the easeofmanufacture and the need to realize diffractionlimited fields without spatial aliasing^{4}; here, λ_{0} denotes the operating wavelength.
Different uniform mappings of the phase domain have been attempted in the past, with acoustic examples raging between 2 and 3bit^{16,18,19} and electromagnetic studies down to 1bit^{26,27,29}. None of these works, however, have discussed how the choice of the quantization level impacts on the fidelity with which the desired field is reproduced (Supplementary Note 1). Here we treat the phase distribution ϕ_{ij} like a 2D image, so that optimization of the AD conversion stage—that is, the process that, given a desired precision in the acoustic field, minimizes the number of different phases to be used and, possibly, the number of elements for each phase (while maintaining the given precision in the reproduced field)—is then analogous to a vector quantization in image compression. Different compression methods are possible (everyone being familiar with classical JPEG protocol), but waveletbased methods^{30} are specifically aimed at determining the lowest number of coefficients necessary for a specified reconstruction quality of localized features^{31,32}, like the abrupt changes of phase in a focusing metasurface (Supplementary Note 1). A discrete wavelet transform (DWT) represents the image over different scales, selecting at each step the key features (low spatial frequencies) and the residual ones (highspatial frequencies). This generates a hierarchical tree of matrices, where the spatial resolution doubles at each step (Supplementary Fig. 1). Once the tree is obtained, the compression procedure works as follows: is computed, up to level K, and the coefficients below a certain threshold value δ are set to zero. The inverse transform is finally computed, containing in general less information and a smaller number of required phases, which are crucially not uniformly distributed either spatially or in the phase domain.
Figure 2a,b shows the phase distribution to form a focus at (0,0,100) mm (see Methods), based on a 16 × 16 grid, where we calculated the DWT using the classic Haar function as parent wavelet (Supplementary Fig. 1). Note here each pixel in Fig. 2a, that is, the original uncompressed image, is λ_{0}/2 in size. As shown in Fig. 2c, the number of unique phases, that is, quanta, needed to realize , decrease with increasing δ: 8 quanta (3bit) are sufficient for a compression rate of 4:1, 6 quanta at 4.6:1 and 4 quanta (2bit) at 8:1. As expected, for all cases, the quanta values are nonuniformly distributed in the interval [0, 2π). To select the thresholding levels to use in the following experiments, we compute the error in approximating the continuous phase distribution with the sum of the squared differences:
Since this quantity is directly related with the precision with which the desired field can be realized (Supplementary Note 2), we target the value ϕ_{E}=0.25 radians, which corresponds to the error in the phase distribution obtained using a uniform 4bit quantization, leading to a 0.1 dB error on the pressure distribution of a focusing metasurface (Supplementary Note 2; Supplementary Fig. 2). As shown in Fig. 2d, both δ=1/8 (11 quanta) and δ=3/16 (8 quanta) fulfil this condition and will be realized experimentally in the following sections. Errors of 1 dB are expected with 8 uniformly distributed quanta (3 dB with 4 uniformly distributed quanta), and they may become as large as 8 dB when 1bit solutions are used (Supplementary Fig. 2).
To properly assign the desired quantized phase values resulting from this analysis, we sculpt next the geometry of our quanta, which we call metamaterial bricks. As we will show in the rest of the presented work, these metamaterial bricks can be used to build various different structures and acoustic devices, emphasizing the power of the simplicity of this concept.
Metamaterial bricks
To form a desired acoustic field with exceptional performance, the metamaterial bricks should possess the ability to transmit sound effectively, locally shift phase with a 2π range, and hold subwavelength spatial resolution to avoid spatial aliasing effects^{33}. By spatially tailoring the geometry of the bricks in a planar metasurface, one can modulate the transmitted waves in nearly arbitrary ways in a specific frequency range. While most of the other studies explore the audible range, our metamaterial bricks are designed for operation in the ultrasonic range at 40 kHz (wavelength λ_{0}≈8.66 mm in air at 25 °C). As shown in Fig. 3a, each brick appears geometrically like a rectangular cuboid with a square baseshape of side λ_{0}/2 and a height of λ_{0}, and consists of an open central channel that delays the incident wave, hence shifting the relative phase of the output. The channel topology is designed to be suitable for microfabrication (see Methods) and was inspired by a pool of alternative designs^{16,17,18,19,20,21,22,23,24,25}, where some of these are based on the pioneering work on coiling space by Liang and Li^{34}. Fully 3D bricks are constructed by extruding the surface of four parallel bars (of variable spacing and length) positioned orthogonally to the wave direction, thereby creating a labyrinth meander. The interbar spacing b_{s} and bar length b_{l} can be effectively tuned, resulting in a phase shift covering a 2π span (Supplementary Fig. 3). Simulations showed that a uniform 4bit quantization of the phase space can reproduce any focused field (with a focal length between 3λ_{0} and 47λ_{0}) with an error of <0.1 dB (Supplementary Fig. 2), so sixteen bricks were designed—corresponding to delays of in phase. Figure 3b shows fullwave simulations (see Methods) of the pressure field distributions when a plane wave travels through each of the 16 selected bricks. While the design of the bricks is not unique (Supplementary Note 3), an important feature of the selected geometry is that the effective acoustic impedance of each brick is matched to that of air, increasing the efficiency of wave transmission and suppressing reflection. In fact, the average transmission magnitude over all units, according to our fullwave simulations, is >98% (Supplementary Table 1). We then manufactured in advance, by rapid prototyping, a set of 16 types of bricks, which could be mounted into a lasercut grid frame structure, where each gridsquare contains 4 bricks in a 2 × 2 assembly (Fig. 3c; Supplementary Video 1). Transmission measurements for three selected bricks (see Methods) confirmed an experimental average transmission of 97±5%, within 10% from the predicted values. Note that functionally, each brick presents a directivity pattern, which depends on its internal topology and on the frequency used. Hence, to reduce unwanted steering effects, the bricks were assembled with a random orientation.
Through a series of numerical simulations and accompanying experiments, we now demonstrate the utility of the concept of quantal metasurfaces in the design of several structures with special functionalities. Specifically, we synthesize metasurfaces, which in isolation apply simple transformations (for example, beam steering, focusing) to a plane incident beam. The lowtransmission loss means that quantal metasurfaces can also be stacked into layers to perform additive transformations; for example, a focusing layer can be combined with a beamsteering layer, to form an offaxis focus, or with an inverting annulus, thus creating a bottleshaped field. The latter configuration is used to successfully levitate a small polystyrene bead at different bitrates (from 4bit to 3bit quantization). While any diffractionlimited field can be theoretically created with just a single layer, stacking brings about additional conveniences for the concept of spatial sound modulators. Due to the additivity of phase delays, each metasurface in the stack can in fact be realized with a lower bitrate, so that a smaller number of brick types is eventually needed. This advantage is particularly evident when the phase quantization is nonuniform.
Steered and focused beams
Using the metamaterial bricks shown in Fig. 3, we now design a blazed transmission grating. Previous studies on reflection/refraction/beam steering with metasurfaces^{16,19,21,25} were typically performed using a sawtooth phase gradient along the output interface. The observed effects were explained in terms of a generalized version of the Snell’s law^{27}, valid only for slow variations of the phase ϕ(x). However, spatial discretization of a metasurface means that the phase variation is a series of steps, and so the behaviour is more completely explained in terms of diffraction theory^{35}.
The space discretized phaseramp is an acoustic diffraction grating that produces beam steering at an angle dependent on the wavelength and on the spatial period, d, over which the range is covered, that is, a ‘line’^{36}. In this configuration, the energy is split between the desired order (typically m=1) and the other orders of diffraction (mainly m=0 and m=−1), which correspond to unwanted directions of propagation. It is worth noting that, when the energy goes mostly into m=−1, one gets what is called negative refraction; in many studies, this is certainly a desired effect. Tuning the gradient of the phase over the line width d, a very efficient energy transfer to the first order can be achieved under certain conditions (that is, blazing), so that in practice only one beam is present, as described by generalized Snell’s law. In optics, blazed gratings are similarly realized with a nematic SLM: the steering angle is determined by the spacing between the lines d, but the gradient of the phase in each line can be optimized to maximize the energy in a single output beam.
Here, we demonstrate the strength of the diffractionbased approach by realizing an experiment that cannot be explained by the generalized Snell’s law. Each line contains a linear phaseramp and a constant phase section of variable length, covering the full 2πrange. We use two λ_{0}/2bricks (respectively of phases ϕ_{1}=5/8π and ϕ_{2}=11/8π) and a channel of variable width (ϕ_{3}=2π). In this way the gradient over the first two bricks stays fixed at , while the line width d varies between λ_{0} and 2λ_{0}. Results in Fig. 4a–c confirm that the diffraction angle decreases as line spacing is increased, but also show that energy gets progressively shifted into secondary lobes as this happens. The farfield diffraction pattern of the grating for varying d is compared in Fig. 4d with the theoretical transmitted angles θ_{t} calculated with the grating equation, mλ_{0}=d sin θ_{t}, and excellent agreement is found. This is also confirmed by experimental results reported in Fig. 4d, which agree perfectly with both numerical and theoretical results.
Another example with special functionality, to which we apply the notion of quantal metasurface, is the steering of a focused beam. To achieve an offaxis beam focusing, we combine a focusing metasurface with a blazed grating in a multilayered structure. Optical studies^{26} have investigated theoretically the concept of multilayered metamaterials, indicating how stacks of metasurfaces can be used to create an arbitrary Green’s function, to either transform a known input into a desired field and/or to perform mathematical operations (that is, differentiation). Here we build on this concept, providing the experimental realization (in acoustics) thanks to the excellent transmission capabilities of the bricks. Here, we first synthesized a metasurface with a uniform 4bit phase quantization for focusing along the axis, at F_{0}=100 nm (see Methods). We then mounted a blazed grating on top of the focusing metasurface; the interlayer separation was arbitrarily fixed to 3λ_{0}/4. The grating metasurface is a result of the previous analysis, and is formed by three λ_{0}/2bricks, respectively of phases ϕ_{1}=5/8π, ϕ_{2}=11/8π and ϕ_{3}=2π, whose performance is simulated in Fig. 5a. The measurements of Fig. 5b show that the 4bit experimental realization performs as predicted by the fullwave simulations. Naturally, we could have used a single metasurface and 16 brick types to perform the desired operation, but stacking makes it possible to achieve the same result using only 8 unique bricks (identified by the waveletbased compression algorithm), as shown in Supplementary Fig. 4. Both fullwave simulations and experiments (Supplementary Fig. 4a,b) suggest in fact only small differences in the focusing field obtained with a 4bit uniform and the 3bit nonuniform phase quantization identified by wavelet compression; in particular, we note a prolonged focal region with a slightly shifted focus with the lowbit quantal metasurface. Also, interestingly, the size of the focal region perpendicular to the axis depends on the lateral dimensions of the metasurface (8λ_{0} in this study): the larger the latter, the tighter the focus. According to the Rayleigh criterion, our source should give a resolution limit of 0.61λ_{0}, whereas the ½width of the spot in experiments was found to be 0.85λ_{0} (Supplementary Fig. 4c). It is also worth noting that our focus is at 50 mm from the axis and therefore outside the boundaries of the metasurface, bringing focused energy where other studies^{16} did not manage.
Bottleshaped beams and acoustic levitation
The final example deals with the use of quantal metasurfaces in realizing a bottle trap, which can also be used as a tractor beam. This field can be created by the superposition of a focus and an additional phase screen in which a central circular region is driven out of phase with the outer regions^{12}. Such beams have been shown to offer stable singlesided acoustic levitation near the focal region. Figure 6 plots the simulated and measured (see Methods) transmitted field maps corresponding to a stack comprising a focusing quantal metasurface on top of a phaseinverting annulus. Here, we compare two levels of quantization for the phase distribution: 4bit and 3bit. In both cases we observed the expected lowpressure ‘quiet’ region near the focus (F_{0}=100 nm. As in Supplementary Fig. 4, the quiet region was shifted further away from the top metasurface when the 3bit phase quantization was applied. These effects were captured both in simulations and measurements, with excellent agreement between the two. For the phase distributions of the synthesized metasurfaces shown here, Supplementary Fig. 4d. A snapshot demonstrating the acoustic levitation of a polystyrene bead with a diameter in the range of 2.1–2.5 mm is shown in Fig. 6d.
The procedure of assembling metamaterial bricks into metasurfaces, and layers to achieve acoustic levitation of a polystyrene sphere is shown in Supplementary Video 1. Herein, in addition to 4 and 3bit quantization (16 and 8 unique bricks, respectively), we also show an acoustic levitation demonstration with 11 unique bricks corresponding to a 2.5:1 compression rate of a focusing metasurface.
Discussion
The results reported in this paper reveal an algorithm for the fabrication of quantal metasurfaces with wellcontrolled arbitrary 2D phase distributions using 3D metamaterial bricks. The concept of our method was inspired from analoguetodigital conversion and image compression. Here, the AD conversion denotes a process in which a continuously variable (analogue) phase distribution is changed, without altering its essential content, into a multilevel (digital) phase. A quantal metasurface becomes in effect a digital phase processing system operating as an SLM. As summarized in Fig. 1, three possible routes can be considered to realize a quantal metasurface. First, we start with a spatial resolution of λ_{0}/2 and 4bit uniform representation of the analogue phase (Fig. 1a). This method is similar to what has been reported by other authors. Second, we treat the phase distribution as an image and use a wavelet transform to perform a lossy compression (Fig. 1b). In this way, the number of bits necessary to represent the analogue phase is reduced to 3, while maintaining the error in reproducing the field below 0.1 dB: this has been explored in Fig. 6 and Supplementary Fig. 4. There is however a lossless alternative (Fig. 1c), where we can exploit the additivity of phase and the hierarchical properties of wavelet transforms. These transforms result in a branched tree to represent the structure of the decomposition; at a given level, each of the branch (that is, metasurface) can be layered to form a stack with a smaller number of unique bricks. An example where this route is applied to a focusing field is reported in Supplementary Fig. 5. Here, the analogue phase is represented with 3.5 bits and a stack of 4 quantal metasurfaces. In this solution, each metasurface is made of ‘blocks’—that is, unit cells of dimension λ_{0}, each made of a 2 × 2 structure of identical bricks. While the spatial resolution of the bricks has to remain lower than λ_{0}/2 to reduce higher order diffraction, the presence of blocks makes mechanical assembly much simpler, which is a key aspect of costeffective devices based, for example, on active metamaterials^{26,27,28}.
Our metamaterial bricks are designed for operating at 40 kHz, nevertheless the presence of the labyrinth meander means that each brick can also work at lower frequencies^{34}. In particular, since the additional phase delay depends on the ratio between the effective length of the meander and the wavelength λ_{0} of the incident wave, each selected brick will have the same transmission performance at frequencies ; here, j=0, 1, 2… is an integer, f_{0} is the design frequency and c_{0} is the speed of sound. Although not demonstrated herein, fullwave simulations confirm this behaviour. The frequency response, and the potential use of amplitude modulation, would require further studies that go beyond the scope of the present work. If realized, full 3D arbitrary acoustic field distributions will open the door to new acoustic devices combining diffraction, scattering and refraction. As a result of advances in metamaterials and the rapidly increasing capabilities for fabricating materials, it is expected that the traditional notion of what constitutes an acoustic device continues to evolve. For instance, in a recent work^{37}, modern rapid prototyping was used to craft monolithic acoustic holograms, enabling complete control over phase, albeit in a static device. Our quantization approach, along with the functional metamaterial bricks, may enable the future development of fully digital spatial sound modulators, which can be controlled in realtime with minimal resources.
Methods
Manufacturing and assembly
The metamaterial bricks were manufactured from thermoplastics using a 3D printer (ProJet HD 3000 Plus), which has a print resolution of 25 μm. Our bricks are then constituted by two materials^{29}: the thermoplastic itself and air. The shoulder bar filets in each brick increased stability during manufacture, as well as contributing to impedance matching. For each position in the grid, we select the phase, among the available ones, which is nearest to the desired analogue value ϕ(x, y). In doing this, we account for the presence of the grid (1 mm wide in Fig. 3c) by making sure that the phase to assign ϕ_{ij} is the one corresponding to an imaginary point at the centre of each brick.
Numerical simulations
Fullwave numerical simulations are carried out by the finite element solver in commercial software COMSOL Multiphysics v5.2. Except for the results regarding each brick, all the other simulations are done in 2D. Perfectly matched layers are imposed on the outer boundaries of simulation domains to prevent reflections.
Transmission measurements
Transmission amplitude was measured for the bricks by mounting 64 identical bricks in an 8 × 8 arrangement. The rest of the grid was filled with absorbing material (to mimic the perfectly matched layers described above). Amplitude was measured with an insertion loss technique using a calibrated B&K microphone (model 4138A015) positioned 30 mm above the source array of a 40 kHz transducer (MA40S4S, Murata Electronics, Japan). Results presented in the text refer to bricks #5, #11 and #13.
Field mapping measurements
Quantitative measurements were obtained using a calibrated B&K microphone (model 4138A015) and a scanning 3D linear stage built in our laboratory. A 2D planar (8 × 8) array of transducers was used to generate a plane wave. Each of the transducers (MA40S4S, Murata Electronics, Japan) has a central frequency of 40 kHz, a beam spread angle of ±40°and sound pressure levels of 120±3 dB (measured on the axis at a z=30 cm). Microphone scans at z=10 mm over the array showed a maximum variation of the phase of 6% across the array, justifying the hypothesis of input plane wave. In our experiments, the transducers are driven at 20 Vpp and operate below the onset of nonlinear effects in air. The measured data are normalized with the maximum amplitude of the corresponding simulated data to aid comparison.
Focusing arrangement
In this study, we consider the centre of coordinates to be at the centre of the top surface of the higher metasurface in a stack. To obtain a focus at distance F_{0} along the axis for an axissymmetrical arrangement, we used the analogue phase distribution:
where r= is the distance of the selected point from the axis, located at (x_{0}, y_{0}). Since the phase ϕ(x, y) should be modulo2π, we selected the central phase ϕ_{0}=π to alter the phase distribution to the centre of Fig. 3b and reduce the phase changes across the metasurface. A similar method, which here simplifies manufacturing, has been used with SLMs to increase their operational speed^{38}. Note that when using a point source, this phase needs to be complemented with an additional phase distribution, but with the transducer source 4 mm below the metasurface, this correction was not necessary.
Data availability
The authors declare that all data supporting the findings of this study are available within the article and its Supplementary Information files. Further information is available from the corresponding author upon reasonable request.
Additional information
How to cite this article: Memoli, G. et al. Metamaterial bricks and quantization of metasurfaces. Nat. Commun. 8, 14608 doi: 10.1038/ncomms14608 (2017).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
 1
Geng, J. Threedimensional display technologies. Adv. Opt. Photon. 5, 456–535 (2013).
 2
Padgett, M. & Di Leonardo, R. Holographic optical tweezers and their relevance to lab on chip devices. Lab Chip 11, 1196–1205 (2011).
 3
Lalonde, R. J., Worthington, A. & Hunt, J. W. Field conjugate acoustic lenses for ultrasound hyperthermia. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 40, 592–602 (1993).
 4
Wooh, S.C. & Shi, Y. Influence of phased array element size on beam steering behavior. Ultrasonics 36, 737–749 (1998).
 5
Turnbull, D. H. D. H. & Foster, F. S. S. Beam steering with pulsed twodimensional transducer arrays. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 38, 320–333 (1991).
 6
Cranston, D. A review of high intensity focused ultrasound in relation to the treatment of renal tumours and other malignancies. Ultrason. Sonochem. 27, 654–658 (2015).
 7
Lipsman, N. et al. MRguided focused ultrasound thalamotomy for essential tremor: A proofofconcept study. Lancet Neurol. 12, 462–468 (2013).
 8
Gélat, P. et al. A comparison of methods for focusing the field of a HIFU array transducer through human ribs. Phys. Med. Biol. 59, 3139–3171 (2014).
 9
Lane, C. J. L., Dunhill, A. K., Drinkwater, B. W. & Wilcox, P. D. The inspection of anisotropic singlecrystal components using a 2D ultrasonic array. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 57, 2742–2752 (2010).
 10
Yang, J., Gan, W. S., Tan, K. S. & Er, M. H. Acoustic beamforming of a parametric speaker comprising ultrasonic transducers. Sensors Actuators A Phys. 125, 91–99 (2005).
 11
Carter, T., Seah, S. A., Long, B., Drinkwater, B. & Subramanian, S. UltraHaptics : multipoint midair haptic feedback for touch surfaces. Proc. UIST 2013, 505–514 (2013).
 12
Marzo, A. et al. Holographic acoustic elements for manipulation of levitated objects. Nat. Commun. 6, 8661 (2015).
 13
Cummer, S. A. et al. Controlling sound with acoustic metamaterials. Nat. Rev. Mater 1, 16001 (2016).
 14
Ma, G. & Sheng, P. Acoustic metamaterials: from local resonances to broad horizons. Sci. Adv. 2, e1501595 (2016).
 15
Pendry, J. B. Negative refraction makes a perfect lens. Phys. Rev. Lett. 85, 3966–3969 (2000).
 16
Xie, Y. et al. Wavefront modulation and subwavelength diffractive acoustics with an acoustic metasurface. Nat. Commun. 5, 5553 (2014).
 17
Zhu, X. et al. Implementation of dispersionfree slow acoustic wave propagation and phase engineering with helicalstructured metamaterials. Nat. Commun. 7, 11731 (2016).
 18
Li, Y., Liang, B., Gu, Z., Zou, X. & Cheng, J. Reflected wavefront manipulation based on ultrathin planar acoustic metasurfaces. Sci. Rep. 3, 2546 (2013).
 19
Li, Y. et al. Experimental realization of full control of reflected waves with subwavelength acoustic metasurfaces. Phys. Rev. Appl. 2, 064002 (2014).
 20
Al Jahdali, R. & Wu, Y. High transmission acoustic focusing by impedancematched acoustic metasurfaces. Appl. Phys. Lett. 108, 031902 (2016).
 21
Mei, J. et al. Controllable transmission and total reflection through an impedancematched acoustic metasurface. New J. Phys. 16, 123007 (2014).
 22
Li, Y. & Assouar, M. B. Threedimensional collimated selfaccelerating beam through acoustic metascreen. Sci. Rep. 5, 17612 (2015).
 23
Jiang, X., Li, Y., Liang, B., Cheng, J. & Zhang, L. Convert acoustic resonances to orbital angular momentum. Phys. Rev. Lett. 117, 034301 (2016).
 24
Li, Y., Jiang, X., Liang, B., Cheng, J. & Zhang, L. Metascreenbased acoustic passive phased array. Phys. Rev. Appl. 4, 024003 (2015).
 25
Li, Y. et al. Theory of metascreenbased acoustic passive phased array. New J. Phys. 18, 043024 (2016).
 26
Silva, A. et al. Performing mathematical operations with metamaterials. Science 343, 160–163 (2014).
 27
Yu, N. et al. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science 334, 333–337 (2011).
 28
Cui, T. J., Qi, M. Q., Wan, X., Zhao, J. & Cheng, Q. Coding metamaterials, digital metamaterials and programming metamaterials. Light Sci. Appl. 3, e218 (2014).
 29
Della Giovampaola, C. & Engheta, N. Digital metamaterials. Nat. Mater. 13, 1115–1121 (2014).
 30
Marcellin, M. W. et al. An overview of quantization in JPEG 2000. Signal Process. Image Commun. 17, 73–84 (2002).
 31
Antonini, M., Barlaud, M., Mathieu, P. & Daubechies, I. Image coding using wavelet transform. IEEE Trans. Image Process 1, 205–220 (1992).
 32
DeVore, R. A., Jawerth, B. & Lucier, B. J. Image compression through wavelet transform coding. IEEE Trans. Inf. Theory 38, 719–746 (1992).
 33
Johnson, D. H. & Dudgeon, D. E. Array Signal Processing:Concepts and Techniques P T R Prentice Hall (1993).
 34
Liang, Z. & Li, J. Extreme acoustic metamaterial by coiling up space. Phys. Rev. Lett. 108, 114301 (2012).
 35
Larouche, S. & Smith, D. R. Reconciliation of generalized refraction with diffraction theory. Opt. Lett. 37, 2391 (2012).
 36
Born, M. & Wolf, E. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light Cambridge University Press (1999).
 37
Melde, K., Mark, A. G., Qiu, T. & Fischer, P. Holograms for acoustics. Nature 537, 518–522 (2016).
 38
Thalhammer, G., Bowman, R. W., Love, G. D., Padgett, M. J. & RitschMarte, M. Speeding up liquid crystal SLMs using overdrive with phase change reduction. Opt. Express 21, 1779 (2013).
Acknowledgements
We would like to thank Tom LlewellynJones from the University of Bristol, who helped with the pressure measurements and Luis F. Veloso from Sussex University, who helped with the supplementary video. We acknowledge funding from the EPSRC project EP/N014197/1, ‘User Interaction with selfsupporting freeform physical objects’.
Author information
Affiliations
Contributions
G.M. and S.S. conceived the study and the quantization analysis. Experiments and data analysis were led by G.M., with contributions from all authors. M.C. led the metamaterial simulations with contributions from M.A. and B.W.D. D.R.S. Optimized the manufacturing technique. G.M. wrote the paper, with contributions from all authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures, Supplementary Notes and Supplementary Table (PDF 1295 kb)
Supplementary Movie 1
This video summarises the content of the paper. (MOV 23761 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Memoli, G., Caleap, M., Asakawa, M. et al. Metamaterial bricks and quantization of metasurfaces. Nat Commun 8, 14608 (2017). https://doi.org/10.1038/ncomms14608
Received:
Accepted:
Published:
Further reading

Three dimensional acoustic tweezers with vortex streaming
Communications Physics (2021)

Emitting longdistance spiral airborne sound using lowprofile planar acoustic antenna
Nature Communications (2021)

Analogue computing with metamaterials
Nature Reviews Materials (2021)

Creation of acoustic vortex knots
Nature Communications (2020)

Design and characterization of an acoustic composite lens with highintensity and directionally controllable focusing
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.