# Long-range nonspreading propagation of sound beam through periodic layered structure

Yurii Zubov, Bahram Djafari-Rouhani, Yuqi Jin, Mathew Sofield, Ezekiel

Walker, Arup Neogi, Arkadii Krokhin

► To cite this version:

Yurii Zubov, Bahram Djafari-Rouhani, Yuqi Jin, Mathew Sofield, Ezekiel Walker, et al. Long-range nonspreading propagation of sound beam through periodic layered structure. Communications Physics, 2020, 3 (1), 155, 8 p. 10.1038/s42005-020-00422-1. hal-03321527

HAL Id: hal-03321527

https://hal.science/hal-03321527

Submitted on 17 Aug 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s42005-020-00422-1

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Long-range nonspreading propagation of sound beam through periodic layered structure

Yurii Zubov<fnref n="1" />, Bahram Djafari-Rouhani², Yuqi Jin ID <fnref n="1" />, Mathew Sofield<fnref n="1" />, Ezekiel Walker ID <fnref n="1" />,³, Arup Neogi ID <fnref n="1" />✉ & Arkadii Krokhin<fnref n="1" />✉

Linear spreading of a wave packet or a Gaussian beam is a fundamental effect known in evolution of quantum state and propagation of optical/acoustic beams. The rate of spreading is determined by the diffraction coefficient *D* which is proportional to the curvature of the isofrequency surface. Here, we analyzed dispersion of sound in a solid-fluid layered structure and found a flex point on the isofrequency curve where *D* vanishes for given direction of propagation and frequency. Nonspreading propagation is experimentally observed in a water steel lattice of 75 periods (~1 meter long) and occurs in the regime of anomalous dispersion and strong acoustic anisotropy when the effective mass along periodicity is close to zero. Under these conditions the incoming beam experiences negative refraction of phase velocity leading to backward wave propagation. The observed effect is explained using a complete set of dynamical equations and our effective medium theory.

<footnote n="1">¹Department of Physics, University of North Texas, 1155 Union Circle #311427, Denton, TX 76203, USA. ²Institut d'Electronique, de Microélectronique et Nanotechnologie, UMR CNRS 8520 Université Lille 1, 59653 Villeneuve d'Ascq, France. ³Echonovus Inc., 1800 S. Loop 288 STE 396, Denton, TX 76205, USA. ✉email: Arup.Neogi@unt.edu; Arkadii.Krokhin@unt.edu</footnote>
The dependence of the index of refraction on wavelength, $n=n(\lambda)$ known as dispersion, is responsible for a wide variety of effects in optics and acoustics. In particular, dispersion is the main reason for broadening of a wave packet propagating in a dispersive medium. Only in a hypothetical dispersionless medium does a signal propagate without changing its shape. However, even in this case a collimated beam of initial width $\sigma_0$ conically spreads. Spreading is due to the presence of the Fourier components with wave vectors making angle $\Delta\theta \sim \lambda/(\sigma_0 n)$ with direction of propagation.

It was demonstrated recently that initially spatiotemporally modulated optical pulses may propagate in a linear medium without dispersion and diffraction¹. In these 3D light bullets, the effects of dispersion and diffraction are not manifested and evolution of the packets is reduced to translational displacement without form change. Initial localization of the wave packet in 3D is achieved due to transversal modulation by a Bessel function and due to longitudinal and temporal modulations by an Airy function. Each function represents an exact solution of the corresponding wave equation which exhibits shape-invariant propagation²⁻⁵.

In acoustics, the area of dispersion/diffraction-free propagation of signals is less developed. Generation and diffraction-free propagation of a bottle-shaped localized acoustic pulses in a homogeneous medium was demonstrated by Zhang et al.⁶. A helical structure metamaterial possessing a large effective index of refraction for audible sound was proposed by Zhu et al.⁷. Despite high values of $n_{eff} \sim 10$, the index of refraction was practically independent of frequency, i.e. a slow sound wave did not suffer from dispersion spreading. Dispersionless propagation of mechanical gravity waves with nonlinear spectrum was observed on the surface of a water tank⁸. In all these experiments, as well as in the experiments with light pulses, the nonspreading propagation becomes possible due to a specially prepared initial wave packet.

In practical applications of acoustics in industry, biology, or medicine the most common signal generated by a transducer is a Gaussian beam. Divergence of a beam reduces sound energy concentration, resolution of acoustic imagining, and depth of penetration of the sound wave into a scanning medium. These and other drawbacks of acoustic devices may be improved by a robust and inexpensive method of collimation of acoustic beams.

In this paper, we report long-distance and nonspreading propagation of an acoustic Gaussian beam through a periodic stack of solid-fluid layers with hyperbolic dispersion. Strong suppression of the dispersion broadening of a Gaussian beam is predicted and observed for a set of angles of incidence and frequencies of sound when the dispersion coefficient vanishes exactly. This effect is due to some peculiarities in the dispersion law of the mixed acoustic mode propagating in strongly anisotropic layered environment.

## Results

**Evolution of a Gaussian beam.** Propagation of a Gaussian beam in a dispersive medium is a standard problem of wave dynamics. Let a monochromatic source generate a pressure beam at $x=0$ of width $\sigma_0$ along axis $y$, $p(x=0, y, t) = p_0 \exp(ik_{0y}y - y^2/4\sigma_0^2 - i\omega_0 t)$. The beam propagates in an anisotropic elastic medium with nonlinear dispersion $\omega = \omega(k_x, k_y)$. The direction of the z-axis is chosen in such a way that the wave vector has two nonzero components, $k_0 = (k_{0x}, k_{0y}, 0)$. To emphasize the effect of nonspreading propagation we assume that the beam is homogeneous along $z$, i.e. its width along $z$ is much greater than along $y$. Under this condition, z-dependence of pressure in the wave can be ignored.

The distribution of pressure in the region $x > 0$ is obtained by integration over $k_y$ of the spatial Fourier components of $p(x=0, y, t)$ with phase-accumulation factor $\exp[i(k_x x + k_y y)]$. Asymptotical behavior of this integral is obtained by quadratic expansion of the diffraction relation $k_x = k_x(k_y, \omega_0)$ near $k_{0y}$ (see, e.g. ref. ⁹)

$$p(x, y, t) = \frac{p_0 e^{i(k_0 \cdot x - \omega_0 t)}}{\sqrt{1 - i \frac{Dx}{2\sigma_0^2}}} \exp \left[ - \frac{(y + xV)^2}{4\sigma^2(x)} \left( 1 + i \frac{Dx}{2\sigma_0^2} \right) \right]. \quad (1)$$

Here $\sigma(x) = \sqrt{\sigma_0^2 + (Dx/4\sigma_0)^2}$ is the coordinate-dependent width of the beam, $V(k_{0y}) = \partial k_x/\partial k_y|_{k_y=k_{0y}}$ is the inclination parameter, and $D(k_{0y}) = \partial^2 k_x/\partial k_y^2|_{k_y=k_{0y}}$ is the diffraction coefficient of the medium⁹⁻¹¹. At short distances, $x < \sigma_0^2/|D(k_{0y})|$, the diffraction can be neglected. At longer distances, $x \gg \sigma_0^2/|D(k_{0y})|$ the width of the beam grows linearly with $x$ and the beam approaches a conical shape with the opening angle $\Delta\theta = \arctan(|D|)$. Thus, the diffraction coefficient $D$ defines the rate of beam spreading in spatial domain. The diffraction is normal if $D < 0$, otherwise it is anomalous.

The inclination parameter $V(k_{0y})$ gives the direction of the transverse shift off the center of the Gaussian beam. The center of the beam follows the straight line $y = -Vx$, which is parallel to the vector of group velocity $V_g = (\partial\omega/\partial k_x, \partial\omega/\partial k_y)$. In an anisotropic medium the group $V_g$ and phase $V_{ph} = (\omega/k)\hat{k}$ velocities are not parallel.

Linear diffraction spreading does not occur if $D(k_{0y}) = 0$. This condition can be realized if for certain frequency and direction of propagation the corresponding isofrequency surface has a region with zero curvature. In particular, in a medium with hyperbolic/elliptic dispersion where $\omega(\mathbf{k}) = \sqrt{k_y^2 c_y^2 \mp k_x^2 c_x^2}$ the diffraction coefficient is always negative, $D(k_y) = -(c_y^2/c_x^4)(\omega^2/k_x^3) < 0$. It vanishes if $c_x = \sqrt{B/\rho_x} = \infty$. This can be realized if the dynamic mass $\rho_x=0$, i.e. if the elastic medium exhibits very strong anisotropy, $\rho_y/\rho_x \to \infty$. In what follows we demonstrate that the dynamic mass density $\rho_x$ may be very small in a periodic layered water-steel lattice that gives rise to long-range nonspreading propagation of an acoustic beam.

Nonspreading propagation of an optical beam due to zero curvature of isofrequency surface was realized in a 3D photonic crystal¹². Later, anomalous diffraction and diffractionless propagation of optical beams in discrete optical waveguides were studied in refs. ⁹,¹⁰,¹³. Anomalies in refraction of sound related to curvature have been observed in ref. ¹⁴. A general consideration of anomalous diffraction in metamaterials is proposed in ref. ¹¹.

Diffractionless propagation of acoustic beams was predicted¹⁵,¹⁶ and observed in the experiment on acoustic transmission through perforated layered elastic structures¹⁷. It was shown that for these structures, the isofrequency curves have hyperbolic topology. While the elastic structures studied in refs. ¹⁵⁻¹⁷ were layered, they possess 3D periodicity due to perforation of solid plates. Our goal is to demonstrate that the diffraction coefficient $D(k_{0y})$ vanishes in a layered solid-fluid structure with one-dimensional periodicity.

**Effective dynamical parameters of a periodic layered system.** Let us consider propagation of sound waves in a periodic binary system with layers of materials a and b and lattice spacing $d = a + b$. In the experiment, material a is water and b is steel. The geometrical and physical parameters of the structure are given in Fig. 1. In the long-wavelength limit $kd \ll 1$, a layered structure
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Band structure of periodic water steel lattice. Black lines show the exact dispersion, $\omega = \omega(k_x)|_{k_x=\text{const}}$ at $k_y = 89 \text{ m}^{-1}$ (see Supplementary Note 2). Red lines represent the same dependence obtained in the effective medium approximation given by Eqs. (3) and (2). Solid lines are for passing bands and dashed lines are for band gaps. Inset is the symmetric unit cell with $a = 9.11 \text{ mm}$, $b = 0.89 \text{ mm}$, and $d = 10 \text{ mm}$. Elastic parameters for water (a) and steel (b) are the following: $\lambda_a = 2 \times 10^9 \text{ Pa}$, $\rho_a = 10^3 \text{ kg m}^{-3}$, $\lambda_b = 1.15 \times 10^{11} \text{ Pa}$, $\mu_b = 7.7 \times 10^{10} \text{ Pa}$, $\rho_b = 7.8 \times 10^3 \text{ kg m}^{-3}$. Note that only half of the Brillouin zone is shown.</figcaption></figure>

behaves like a homogeneous anisotropic fluid with effective mass densities $\rho_x$ and $\rho_y = \rho_z$ and effective elastic modulus $B$. Frequency-dependent effective parameters (EPs) of a layered structure can be calculated using the method of homogenization¹⁸ for a binary structure with high acoustic contrast between the constituents. Generalization of the results¹⁸ to a structure with arbitrary contrast is given in Supplementary Note 1:

$$ 
\begin{align} 
\rho_x &= \frac{\rho_a \frac{k_b}{k_a} \sin \frac{k_a a}{2} \cot \frac{k_b b}{2} + \rho_b \cos \frac{k_b a}{2}}{1 + \frac{k_b}{k_a} \sin \frac{k_a a}{2} \cot \frac{k_b b}{2} - 2 \frac{\rho_b}{\rho_a} \sin^2 \frac{k_a a}{4}}, \\ 
\rho_y = \rho_z &= \rho_a \frac{k_a d}{2} \frac{\cos \frac{k_a a}{2} \cot \frac{k_b b}{2} - \frac{\rho_a k_b}{\rho_b k_a} \sin \frac{k_b a}{2}}{\sin \frac{k_a a}{2} \cot \frac{k_b b}{2} + \frac{\rho_a}{\rho_b} \left[ \frac{k_a}{k_b} - 2 \frac{k_b}{k_a} \sin^2 \frac{k_a a}{4} \right]}, \tag{2} \\ 
B &= \lambda_a \frac{k_a d}{2} \frac{\cos \frac{k_a a}{2} \cot \frac{k_b b}{2} - \frac{\rho_a k_b}{\rho_b k_a} \sin \frac{k_a a}{2}}{\sin \frac{k_a a}{2} \cot \frac{k_b b}{2} + \frac{\lambda_a k_a}{\lambda_b k_b} - 2 \frac{k_b \rho_a}{k_a \rho_b} \sin^2 \frac{k_a a}{4}}. 
\end{align} 
$$

Here $k_{a(b)} = \sqrt{(\omega/c_{a(b)})^2 - k_y^2}$ is the x-component of the wave vector in the layer a(b) and $c_{a(b)} = \sqrt{\lambda_{a(b)}/\rho_{a(b)}}$ is the velocity of longitudinal sound. To simplify calculations the transverse mode in the solid layer was ignored in derivation of the EPs Eq. (2). However, this approximation turns out to be sufficient to semi-quantitatively analyze the effect of nonspreading propagation. Note that the EPs depend on $\omega$ and $k$ accounting for temporal and spatial dispersions. In the quasi-static limit, $\omega \to 0$ formulas (2) are reduced to Rytov's results¹⁹, $\rho_x = (\alpha\rho_a + b\rho_b)/d$, $\rho_y^{-1} = (\alpha\rho_a^{-1} + b\rho_b^{-1})/d$, and $B^{-1} = (\lambda_a^{-1} + b\lambda_b^{-1})/d$. Frequency-dependent components $\rho_x$ and $\rho_y$ of the effective density tensor and the effective compliance are shown in Supplementary Figs. 1 and 2.

**Analysis of isofrequency curves.** The dispersion relation

$$ \omega^2/B = k_x^2/\rho_x + k_y^2/\rho_y \quad (3) $$

with the EPs (2) reproduces well the long-wavelength region of the band structure obtained for the water-steel layers. In Fig. 1 the first three transmission bands are plotted for sound Bloch

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Structure of the isofrequency curves. Three isofrequency curves representing elliptic-like (green) and hyperbolic-like (black and blue) dispersion. Bloch waves have elliptic-like dispersion in the first passing band, and hyperbolic-like dispersion in the second band. Transition occurs through the band gap. Dash-dotted blue curve is the result of the effective medium theory (EMT) obtained from Eqs. (2) and (3). The frequencies corresponding to hyperbolic-like dispersion lie slightly above and slightly below the saddle point at $f_0 = 104.8 \text{ kHz}$. The saddle point is the point of crossing of two parabolas marked by red lines in the inset. The arrows show the directions of the phase and group velocities for the Bloch waves with the same projection $k_y$, which is conserved at the water-steel boundaries. The crossing of the vertical line $k_y = 89 (k_y = 71)$ with solid (dash-dotted) blue curve occurs at $k_x = -51 \text{ m}^{-1}$ ($k_x = -71 \text{ m}^{-1}$) that corresponds to the angle $\gamma = 60^\circ$ ($\gamma = 45^\circ$).</figcaption></figure>

wave with Bloch vector $k_x$ varying within the first half of the Brillouin zone and conserved $k_y = 89 \text{ m}^{-1}$ propagating in water-steel periodic lattice. Exact dispersion curves account for transverse and longitudinal modes are shown in black and those obtained in the effective medium approximation in red. The second band has anomalous dispersion where $\partial\omega/\partial k_x < 0$. It starts at frequency $f = 106.8 \text{ kHz}$. When the angle $\gamma = 0$ the band edge is red-shifted to the frequency $f_0 = 104.8 \text{ kHz}$ (see Supplementary Fig. 3).

Near the frequency $f_0$ the dispersion relation, $f=f(k_x, k_y)$ is represented by a saddle surface shown in the inset to Fig. 2. Horizontal planes cutting the saddle surface above and below the level of $f_0$ give two sets of isofrequency curves topologically similar to hyperbolas. These two sets are separated by the isofrequency contour $f(k_x, k_y) = f_0$ as shown in Supplementary Fig. 4. The separatrix is not a straightline since the curves are not canonical hyperbolas. Two orthogonal parabolas passing through the saddle point are highlighted in red in the inset to Fig. 2. Projections of these parabolas on the horizontal plane form a local set of coordinates for the hyperbolic-like curves. Each coordinate line is the axis of symmetry of the hyperbolic-like curves. In Fig. 2 black and blue curves represent the two sets of hyperbolic-like isofrequency curves. The green line represents an isofrequency curve for $f = 24 \text{ kHz}$ with elliptic-like (normal) dispersion corresponding to the first transmitting band in Fig. 1.

For propagation along the lattice axis, $k_y = 0$, the band gap occupies the interval $f_1 < f < f_2$, where the gap edges $f_1 = 46.7 \text{ kHz}$ and $f_2 = 77.5 \text{ kHz}$ are reached at the edge of the Brillouin zone, $k_x = \pi/d = 314 \text{ m}^{-1}$. The isofrequency curves in the first passing band ($f < f_1$) are closed since equation $f=f(k_x, k_y)$ has real solutions for any direction of propagation. For the frequencies slightly above $f_1$ a narrow interval of the directions of propagation almost parallel to the axis, $k_y \ll k_x$ becomes inaccessible. This interval increases with $f$. Thus, the topological transition from closed to opened isofrequency curves occurs at $f=f_1$. Evolution of
isofrequency curves with frequency is shown in Supplementary Figs. 4–6.

The two blue curves are the plots for $f=105.1$ kHz, which belongs to the second transmission band with anomalous dispersion. The solid blue curve is obtained from the exact dispersion relation derived for a solid–liquid layered structure taking into account longitudinal and transverse modes in metal plates (see Supplementary Note 2). The dash-dotted blue line is the result that follows from the effective medium theory, where the transverse mode is neglected. Because of this approximation and also since the parameter $k_x d \approx 0.5$ is not sufficiently small to provide the quantitative validity of the effective medium theory, the solid and dash-dotted blue lines deviate from each other.

Topologically, both blue curves are similar to hyperbolas. However, they are not given by a canonical second-order hyperbolic equation. The dash-dotted blue curve is obtained from Eqs. (2) and (3). Since the EPs in Eq. (2) depend on $k_y$, Eq. (3) becomes a transcendental equation. Unlike canonical hyperbolas, both blue curves contain regions with positive, negative curvatures, and a flex point. Blue solid vertical line in Fig. 2 marks the position of the exact flex point at $k_y = 0.89 \text{ cm}^{-1}$. It corresponds to the direction of the Bloch vector $\mathbf{k}$ that makes an angle $\gamma = \arctan(k_y/k_x) \approx 60^\circ$ with the negative direction of axis $x$. Negative $x$-projections of vector $\mathbf{k}$ and phase velocity $V_{\text{ph,ref}}^{(x)}$ are due to anomalous dispersion. The product $V_{\text{ph,ref}}^{(x)} V_{\text{g,ref}}^{(x)} < 0$ where the group velocity is always positive. Thus, the incoming wave suffers negative phase velocity refraction at the water–lattice interface and backward propagation inside the lattice. Similar effects were predicted for perforated layered structures$^{16}$.

Refraction of a sound beam of $f=105.1$ kHz at water–lattice interface is demonstrated by Supplementary Movie 1. It is seen that the wave fronts in water propagate towards the lattice, i.e. along positive direction of axis $x$. However, entering the lattice under the angle of incidence $\theta = 10^\circ$ the $x$-projection of their velocity is reversed. It is a signature of negative phase velocity refraction. It is interesting to note that negative group velocity refraction cannot be observed at normal incidence. Unlike this, negative phase velocity refraction and backward wave propagation exist even for normal incidence. This Supplementary Movie 2 demonstrates backward wave propagation in steel–water structure for normal incidence of acoustic beam with frequency $f=104.5$ kHz lying in the second transmission band with anomalous dispersion. However, the dispersion is normal in the third band ($f=155.4$ kHz) that results in positive phase velocity refraction (see Supplementary Movie 3).

The flex point on the dash-dotted blue curve in Fig. 2 lies at $k_y \approx 71 \text{ m}^{-1}$ (see Fig. 2). At this point the dynamic density $\rho_x \approx 32 \text{ kg/m}^3$ is much less than the densities of the lattice constituents. It does not vanish exactly (unlike the case of elliptic/hyperbolic dispersion) since the EPs in Eq. (3) are $k_y$-dependent. We may conclude that the nonspreading propagation of acoustic beam occurs under the condition of very strong acoustic anisotropy, $\rho_y/\rho_x \approx 40$ (see Supplementary Note 1 and Supplementary Fig. 1). One more visible signature is that at the flex point, group and phase velocities are almost orthogonal. Also the dispersion curve of the second transmission band in Fig. 1 is quite flat, that leads to the inequality $|V_g| \ll |V_{\text{ph}}|$, which is related to strong scattering at the interfaces$^{20}$.

The isofrequency curves in Fig. 2 obtained from the exact solution of equations of elasticity exhibit discontinuous cuts. Sound waves with frequency and direction corresponding to the cut do not propagate though a given layered structure. The effect of zero acoustic transparency through a solid plate surrounded by

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Nonspreading transmission of acoustic Gaussian beam through 100 periods of the water steel layered structure. 2D map of pressure produced by sound beam. In the incoming beam $V_g = V_{\text{ph}}$. In the refracted beam these vectors are almost orthogonal, in agreement with their orientation in the isofrequency contour for $f=105.1$ kHz. The angle of incidence is $\theta = 10^\circ$ and the angle of refraction for phase velocity is $\pi-\gamma$, where $\gamma = 60^\circ$. The angle of refraction for group velocity is $\gamma_g = 35^\circ$.</figcaption></figure>

a fluid was recently predicted by Quotane et al. $^{21}$. This effect is due to a narrow Fano resonance caused by interference between a discrete trapped mode of a solid plate and a continuum of acoustic modes of a fluid. In the experiment, the direction of propagation of this Fano resonance is realized at must be ruled out. Of course, the blue dash-dotted line does not show Fano resonance since the EPs are obtained neglecting the transverse mode.

Using Snell’s law $\sin \theta/\sin \gamma = c_a/V_{\text{ph,ref}}$ the angle of incidence is calculated to be $\theta \approx 10^\circ$. This angle is counted from the positive direction of axis $x$ since the incident wave propagates in water where $V_{\text{g,inc}} \cdot V_{\text{ph,inc}} = c_a^2 > 0$. Note, that the angle of incidence calculated using the effective medium theory is $9^\circ$ that is close to the exact result (see Supplementary Note 2).

**Numerical modeling and experimental results.** Numerically modeled transmission of a sound beam of frequency $f=105.1$ kHz and angle of incidence $\theta = 10^\circ$ through the steel–water structure is shown in Fig. 3. The width of the beam in the $x-y$ plane practically does not change along its pass through a 100-period long lattice since for the selected parameters the diffraction coefficient $D = \partial^2 k_x / \partial k_y^2$ vanishes. However, in the 3D map of sound intensity shown in Fig. 4, the width of the beam in the $y-z$ plane grows linearly because the corresponding diffraction coefficient $\partial^2 k_x / \partial k_z^2 \neq 0$ (see Supplementary Note 3). The spreading of the beam eventually results in interaction between the beam and the top–bottom boundaries that manifested in unwanted multiple reflections from the acrylic holders. However, the net effect of these reflections is minimal since the impedance mismatch between the acrylic and water produces a low reflection coefficient of ~8%. Also only a small part of the acoustic beam expands over the edges of $17.8 \times 17.8 \text{ cm}$ plates, i.e. very small part of the signal that reaches the receiving transducer suffers from multiple reflections.

The directions of phase and group velocities in the incident and refracted beams are indicated by arrows. The direction of the group velocity in the refracted beam can be calculated from Snell’s law using the definition of the group index of refraction$^{16}$

$$n_g = - \frac{k_y c_a / \omega}{\partial k_x / \partial k_y} \sqrt{1 + \left( \frac{\partial k_x}{\partial k_y} \right)^2} \quad (4)$$
Since in Fig. 2 the derivative $\partial k_x/\partial k_y$ is negative at the flex point ($k_y = 0.89 \text{ cm}^{-1}$ with $V_g^{(x)} > 0$), the beam exhibits positive refraction for group velocity, $n_g > 0$. The angle of refraction calculated from Eq. (4) is $\gamma_g = 35^\circ$.

This angle defines the geometry of the fabricated lattice. Due to strong refraction of the acoustic beam at water-lattice interface, the steel plates in the structure are vertically shifted as shown in

<figure><img src="image_7.png" /><figcaption>**Fig. 4 3D modeling of transmission of acoustic Gaussian beam through water steel layered structure.** 3D distribution of sound intensity in the layered structure. The width of the beam remains almost unchanged along axis y. However, it grows linearly along axis z in the far field. Black dots mark the lines of equal intensity in the transversal cross-section of the beam.</figcaption></figure>

Fig. 5a, b. The input Gaussian beam of width $\sigma_0 \approx 2$ cm propagates a distance $75/\cos 35^\circ \approx 88$ cm through 75-period long lattice. A comparison of the beam spreading in free water versus the lattice is given in the plot inset of Fig. 5b. At 71 cm from the emitter, the beam spreads from $\approx 2.7$ to $\approx 10$ cm through free water. However, the same input spreads to only $\approx 3.5$ cm when propagating through the lattice. Thus, the effect of diffraction broadening is strongly suppressed in the lattice supporting the contention that the diffraction coefficient D is close to zero.

The condition of diffractionless propagation, $D=0$, establishes a dependence among three variables, $k_x, k_y$, and $f$. This dependence can be reduced to the relation between the angle of incidence $\theta$ and frequency $f$ (see Supplementary Note 3). The dependence $f=f(\theta)$ is plotted in Supplementary Fig. 7. The frequency $f=105.1$ kHz is selected as a point on the curve $f=f(\theta)$ corresponding to relatively small angle of incidence ($\theta=10^\circ$), fitting the second transmission band close to the $\Gamma$ point, and lying away from Fano resonance.

The intensity of the signal (at the center of the beam) transmitted through the lattice is reduced by 930 times as compared to the signal intensity at the emitting transducer. The ratio of the intensities averaged over the beam profile is about 500. The reduction of the transmitted signal is due to strong reflection at the first and last metal plates, to the viscous losses and scattering inside the lattice, and to losses outside the lattice. The latter factor is related to the strongly inhomogeneous acoustic field generated by the transducer displaced 8 cm away from the first plate. Sound intensity measured near the first plate is $\approx 5.5$ times lower than at the transducer. The averaged signal

<figure><img src="image_8.png" /><figcaption>**Fig. 5 Experimental measurement of nonspreading transmission of sound.** a 75-period long structure of steel plates immersed in water. b Scheme of the experiment and color scale distribution of pressure in the input and output beams. Inset shows normalized transversal profiles of intensity in the input (black) and output (red) beams. Dashed blue curve is for sound beam propagating in free water measured at the distance of 71 cm from the emitter. The normalizing factors calculated for average sound intensity for red and blue curves are 930 and 6.5, respectively. Fluctuations on the red line are due to sharp decay of the signal away from its axis where amplitude of the signal drops to the level of noise. c Vertical displacement of the collimated beam with slight variation of frequency of sound and the corresponding images of the input and output beams.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Variation of the phase for forward and backward propagation. The phase of propagating wave was measured (solid lines) and numerically calculated (dashed lines) at the positions marked by stars inside the lattice. Frequency 52.7 kHz (112.7 kHz) belongs to the first (second) passing band with normal (anomalous) dispersion. The experiment was performed with a 10-period long structure as schematically shown in the figure. The phases are calculated via fast Fourier transform of the average of 512 independent time windows. The error in phase detection does not exceed $2\pi f \Delta t$, where $\Delta t = 10^{-8}$ s is the time interval used for discretization of the input signal. For both frequencies used in this experiment, the calculated error is $-10^{-3}$, too small to be seen on the scale of 1 radian.</figcaption></figure>

intensity also slightly decays (≈2 times) on the way (4 cm) from the last plate to the receiving hydrophone. Thus, we estimate that, passing through the lattice, the signal intensity is reduced by 500/5.5 × 2 = 45 times. Numerical modeling shown in Fig. 3 gives the decay of the signal intensity by 36 times. This includes the losses due to impedances mismatch and to viscosity. We attribute 45/36 ≈ 1.3 times higher losses in the experiment mostly to scattering at period fluctuations and absorption of sound by small amount of air bubbles which are always present in water. The transmission can be optimized by replacing stainless steel plates by plates of softer material (see Supplementary Note 4 and Supplementary Fig. 8).

The regime of backward propagation of sound inside the lattice was demonstrated by measuring the phases of the wave at six successive unit cells. The results are presented in Fig. 6 and the details are given in Supplementary Note 5. Experimental and numerical results are shown by squares and dots connected by solid and dashed line, respectively. Red and blue symbols show the phases of two passing waves with frequencies lying in the first (red) and second (blue) passing bands. In the experiment the phases of the Bloch wave that exhibit nonspreading propagation at $f = 112.7$ kHz (second passing band) gradually decrease that serves as a signature of backward propagation. Unlike this, the phases of the wave with normal dispersion with $f = 52.7$ kHz (first passing band) gradually increase. Numerically calculated phase behaves similarly. Since in the experiment the nonspreading propagation was observed at the frequency ≈7 kHz higher than the theory predicts, only qualitative agreement is expected. Note that for backward propagation the experimental and numerical phases change much slower than for forward propagation. This occurs because in both cases the frequencies are close to the upper

edge of the bandgap. Exactly at the band edge spatial variation of the phase is compensated by its temporal variation since a propagating Bloch wave becomes a standing wave.

**Discussion.** The divergenceless propagation results in a collimated beam as shown in Fig. 5b. The emerging beam from the lattice does not exhibit strong spreading over the 6 cm from the end of the lattice to the end of the scanned area. From Fig. 5b inset, it can also be seen that propagation through the lattice maintains a strong intensity gradient at the beam edges as compared with the weaker intensity gradient observed at comparable distance without the lattice. Collimation of acoustic beam and energy funneling were previously observed due to diffraction of incoming beam at corrugated surfaces$^{22-24}$. This method partially reduces transverse diffraction of acoustic beam. Resonant interaction of the incoming sound with the surface acoustic wave excited at a corrugated surface also reduces angular spreading of the beam, while numerical simulations$^{22}$ show presence of lateral diffraction maxima. Lateral diffraction disappears if the corresponding isofrequency surface has a flat region$^{11,14,25}$. Usually, the curvature never vanishes completely, leading to some spreading. Experimental realization of a flat band in acoustics was demonstrated in a 2D phononic crystal$^{25}$ of steel rods in water, and more recently in perforated Plexiglas plates in air$^{17}$. Simultaneous collimated propagation of acoustic and electromagnetic waves in phonic crystal of silicon rods in air was numerically predicted$^{26}$. Dissipative losses were not taken into account and high level of collimation is achieved due to existence of a single flex point on isofrequency line.

In these periodic structures where collimation of a sound beam was predicted or observed, sound usually propagates through narrow apertures that strongly enhance viscous losses. Therefore, the length of the collimated beam did not exceed 10 cm. The scheme proposed here is free from such drawback. Much less energy losses allow much longer collimated propagation. Also the angular width of the beam is much less since it is required that curvature vanishes exactly only at a single point. From practical point of view the 1D periodic lattice proposed here is much easier to fabricate.

The exit position of the outgoing beam exhibits strong dependence on frequency. Experimentally, nonspreading propagation occurs at 112.17 kHz. Even small variation of the frequency leads to noticeable vertical displacement of the output beam, as shown by the colored beam profiles in Fig. 5c. This displacement is due to dispersion of the refraction coefficient Eq. (4). The graph in Fig. 5c shows calculated vertical displacement caused by small variation of the refraction coefficient (4), and the angle of refraction, with frequency. Calculated and measured displacements are in reasonable agreement.

While the experimental results confirm most of our theoretical predictions for nonspreading propagation, there is a disagreement over the frequency. The effect is observed at frequency 7 kHz higher than the theory predicts. We attribute this difference to partial vibration coupling and subsequent minor mode shifting of the lattice modes due to the two acrylic plates which hold the metal plates on the top and bottom (see Fig. 5a). These holders are not sufficiently massive to rule out sound transmission, therefore the stiffness of the system as a whole turns out to be higher than that of the steel-water lattice.

We have demonstrated long-range nonspreading propagation of sound through a layered solid-liquid structure with one-dimensional periodicity. Diffraction, usually leading to linear divergence of the beam, is avoided due to vanishing of the diffraction coefficient. Since this coefficient is proportional to the curvature of the isofrequency curve, the possibility of
diffractionless propagation depends on the presence or absence of a flex point on the curve and the direction of propagation is defined by the position of the flex point. The isofrequency curves are obtained from the exact solution of the wave equation in the layered structure and from the developed effective medium theory, which accounts for temporal and spatial dispersions of the EPs. Both methods give close results. We have shown that the effect of nonspreading transmission is accompanied by two other anomalies: negative phase velocity refraction and strong acoustic anisotropy. The reported results may find applications in many areas where robust transmission of narrow acoustic beams is required.

## Methods

**Experiment.** Our sample is assembled from 75 pieces of 17.78 cm square shape grade 409 stainless steel plates with thickness $b = 0.891$ mm. The plates are 10 mm apart, center-to-center with 0.1 mm tolerance. The plates are held in place by specially fabricated acrylic top and bottom plates.

The acoustic experiment was conducted in a large deionized water tank (150 × 90 × 60 cm), in a bistatic setup. Signal generation was achieved with a Licroy WaveStation2012 Arbitrary Waveform Generator set to 20 V peak-to-peak, sweeping from 80 to 120 kHz in 80 s. A 25 mm diameter Ultran plane wave NCG100-D25 immersion transducer served as the emitter. Detection was done using a Müller-Platte 1 mm needle hydrophone connected to a Tektronix MDO 3024b in spectrum analyzer mode. The emitting transducer was rotated in the x-y plane by a Thorlabs CR1-Z6 angular translation stage for variable incident angles. The needle hydrophone was translated in the x-y plane by a two-dimensional 1500 mm by 1200 mm Newmark eBelt linear translation stage controlled by a Newmark NSG-G2-X2 translation stage controller. Calibration of the angle and subsequent alignment was performed for each experiment with angular tolerance of 0.5. An automated raster scan was completed by a pre-prepared MATLAB code with data acquired at each point for 90s. The incident angle with respect to the lattice plane was 10°. For the incident scanned area, the beam profile was measured 60 mm from the incident lattice plate in a 9 cm × 4 cm area in 3 and 10 mm steps, respectively (Fig. 5). The output signal was measured starting 20 mm from the last plate with a scanned area of 9 cm × 4 cm in 3 and 10 mm intervals.

**Theory.** The EPs Eq. (2) are obtained from two dynamical conditions that define the response of elastic medium on external perturbation.

First, the average acceleration of the unit cell in a Bloch wave equals to that calculated for the equal volume of the effective medium under action of propagating plane wave of same amplitude

$$ \begin{aligned} \bar{a}_i &= -\frac{1}{d} \int_{-d/2}^{d/2} \frac{\partial p_{uc}(x, y)}{\partial x_i} \frac{dx}{dx_i} \frac{\rho(x)}{\rho(x_i)} \\ &= -\frac{1}{d\rho_i} \int_{-d/2}^{d/2} \frac{\partial p_h(x, y)}{\partial x_i} dx, \quad i = x, y. \end{aligned} \quad (5) $$

Here $p_{uc}(x, y)$ is the distribution of pressure within the unit cell, $p_h(x, y)$ is the distribution of pressure in a homogeneous effective medium, and $\rho(x)$ is the coordinate-dependent density within the unit cell.

Second, the average deformations are also equal

$$ \begin{aligned} \frac{\Delta V}{V} &= -\frac{1}{d} \int_{-d/2}^{d/2} \frac{p_{uc}(x, y)}{\lambda(x)} dx \\ &= -\frac{1}{Bd} \int_{-d/2}^{d/2} p_h(x, y) dx. \end{aligned} \quad (6) $$

Here $\lambda(x)$ is the coordinate-dependent bulk modulus within the unit cell. Equations (5) and (6) define the EPs not only for a binary structure but for any distribution of density $\rho(x)$ and elastic modulus $\lambda(x)$ within the unit cell.

Derivation of the EPs $\rho_{x}$, $\rho_{y}$ and B is given in Supplementary Note 1.

The dispersion relation for fluid-solid lattice is obtained from the condition of matching the solution of the wave equations for longitudinal mode in layer a (fluid) with the solution for longitudinal and transverse modes in layer b (solid). Together with the Bloch theorem for periodic medium, this leads to a set of 6 × 6 linear homogeneous equations. Numerical simulations are conducted by using an acoustic model in commercial software (COMSOL MULTIPHYSICS).

## Data availability

All data used in this research are available from the corresponding author upon request.

Received: 9 April 2020; Accepted: 10 August 2020;

Published online: 02 September 2020

## References

1. Chong, A., Renninger, W. H., Christodoulides, D. N. & Wise, F. W. Airy-Bessel wave packets as versatile linear light bullets. *Nat. Photon.* **4**, 103–106 (2010).

2. Ignatovich, V. K. Nonspreading wave packets in quantum mechanics. *Found. Phys.* **8**, 565 (1978).

3. Berry, M. V. & Balazs, N. L. Nonspreading wave packets. *Am. J. Phys.* **47**, 264–267 (1979).

4. Durnin, J., Miceli, J. J. & Eberly, J. H. Diffraction-free beams. *Phys. Rev. Lett.* **58**, 1499–1501 (1987).

5. Gori, F., Guattari, G. & Padovani, C. Bessel-Gauss beams. *Opt. Commun.* **64**, 491–495 (1987).

6. Zhang, P. et al. Generation of acoustic self-bending and bottle beams by phase engineering. *Nat. Commun.* **5**, 4316 (2014).

7. Zhu, X. et al. Implementation of dispersion-free slow acoustic wave propagation and phase engineering with helical-structured metamaterials. *Nat. Commun.* **7**, 11731 (2016).

8. Fu, S., Tsur, Y., Zhou, J., Shemer, L. & Arie, A. Propagation dynamics of nonspreading cosine-Gauss water-wave pulses. *Phys. Rev. Lett.* **115**, 254501 (2015).

9. Eisenberg, H. S., Silberberg, Y., Morandotti, R. & Aitchison, J. S. Diffraction management. *Phys. Rev. Lett.* **85**, 1863 (2000).

10. Pertsch, T., Zentgraf, T., Peschel, U., Bräuer, A. & Lederer, F. Anomalous refraction and diffraction in discrete optical systems. *Phys. Rev. Lett.* **88**, 093901 (2002).

11. Paul, T., Rockstuhl, C., Menzel, C. & Lederer, F. Anomalous refraction, diffraction, and imaging in metamaterials. *Phys. Rev. B* **79**, 115430 (2009).

12. Kosaka, H. et al. Self-collimating phenomena in photonic crystals. *Appl. Phys. Lett.* **74**, 1212 (1999).

13. Christodoulides, D. N., Lederer, F. & Silberberg, Y. Discretizing light behaviour in linear and nonlinear waveguide lattices. *Nature* **424**, 817 (2003).

14. Bucay, J. et al. Positive, negative, zero refraction, and beam splitting in a solid/air phononic crystal: theoretical and experimental study. *Phys. Rev. B* **79**, 214305 (2009).

15. Christensen, J. & García de Abajo, F. J. Anisotropic metamaterials for full control of acoustic waves. *Phys. Rev. Lett.* **108**, 124301 (2012).

16. Christensen, J. & García de Abajo, F. J. Negative refraction and backward waves in layered acoustic metamaterials. *Phys. Rev. B* **86**, 024301 (2012).

17. García-Chocano, V. M., Christensen, J. & Sánchez-Dehesa, J. Negative refraction and energy funneling by hyperbolic materials: an experimental demonstration in acoustics. *Phys. Rev. Lett.* **112**, 144301 (2014).

18. Yu, Zubov, Djafari-Rouhani, B. & Krokhin, A. Dynamical effective parameters of elastic superlattice with strong acoustic contrast between the constituents. *Low Temp. Phys.* **44**, 1280 (2018).

19. Rytov, S. M. The acoustical properties of a finely layered media. *Soviet Phys. Acoustics* **2**, 67 (1956).

20. Godin, Yu. A. & Vainberg, B. Slowdown of group velocity in periodic waveguides. Preprint at 1912.06964 [cond-mat.mtrl-sci] (2019).

21. Quotane, I., Houssaine, E., Boudouti, E. H. & Djafari-Rouhani, B. Trapped-mode-induced Fano resonance and acoustical transparency in a one-dimensional solid-fluid phononic crystal. *Phys. Rev. B* **97**, 024304 (2018).

22. Christensen, J., Fernandez-Dominguez, A. I., de Leon-Perez, F., Martin-Moreno, L. & Garcia-Vidal, F. J. Collimation of sound assisted by acoustic surface waves. *Nat. Phys.* **3**, 851 (2007).

23. Zhou, Y. et al. Acoustic surface evanescent wave and its dominant contribution to extraordinary acoustic transmission and collimation of sound. *Phys. Rev. Lett.* **104**, 164301 (2010).

24. Lin, Z. et al. A collimated focused ultrasound beam of high acoustic transmission and minimum diffraction achieved by using a lens with subwavelength structures. *Appl. Phys. Lett.* **107**, 113505 (2015).

25. Espinosa, V., Sánchez-Morcillo, V. J., Staliunas, K., Pérez-Arjona, I. & Redondo, J. Subdiffractive propagation of ultrasound in sonic crystals. *Phys. Rev. B* **76**, 140302 (R) (2007).

26. Tan, Y. et al. Simultaneous beam guides of electromagnetic and acoustic waves in defect-free phoxonic crystals using self-collimation effect. *Appl. Phys. Express* **12**, 06201 (2019).

## Acknowledgements

The authors are thankful to Hyeonu Heo, Tracy Lynch, Mark Lanier for technical support. This work is supported by an EFRI grant from the National Science Foundation (Grant No. 1741677). A.A.K. is thankful to the University of Lille for hospitality where part of this work has been accomplished. M.S. acknowledges the support from the NSF REM summer program. The seed funding from "The Advanced Materials and Manufacturing Processes Institute (AMMPI)" at the University of North Texas is also gratefully acknowledged.
**Authors contributions**

Y.Z. developed the effective medium theory, analyzed the isofrequency contours, and performed numerical simulations. B.D.-R. proposed the model for the effective medium and analyzed the experimental results. Y.J., M.S., and E.W. designed and performed the experiment, A.N. and A.K. managed the project and wrote the manuscript with the input of all co-authors.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** is available for this paper at https://doi.org/10.1038/s42005-020-00422-1.

**Correspondence** and requests for materials should be addressed to A.N. or A.K.

**Reprints and permission information** is available at http://www.nature.com/reprints

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

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020