Abstract
The concept of valley pseudospin, labelling quantum states of energy extrema in momentum space, is attracting attention1,2,3,4,5,6,7,8,9,10,11,12,13 because of its potential as a new type of information carrier. Compared with the non-topological bulk valley transport, realized soon after predictions1,2,3,4,5, topological valley transport in domain walls6,7,8,9,10,11,12,13 is extremely challenging owing to the inter-valley scattering inevitably induced by atomic-scale imperfections—but an electronic signature was recently observed in bilayer graphene12,13. Here, we report the experimental observation of topological valley transport of sound in sonic crystals. The macroscopic nature of sonic crystals permits a flexible and accurate design of domain walls. In addition to a direct visualization of the valley-selective edge modes through spatial scanning of the sound field, reflection immunity is observed in sharply curved interfaces. The topologically protected interface transport of sound, strikingly different from that in traditional sound waveguides14,15, may serve as the basis for designing devices with unconventional functions.
Similar content being viewed by others
Main
The macroscopic controllability enables the sonic crystals (SCs) to be tractable classical counterparts for exploring complex quantum physics requiring atomic-scale manipulations. Two-dimensional (2D) hexagonal crystals, such as graphene, bilayer graphene, and transition-metal dichalcogenides, exhibit a pair of degenerate states at the inequivalent K and K′ valleys in momentum space. Inversion symmetry breaking generically gaps the degeneracy and gives rise to a tantalizing valley Hall insulator1,2,3,4,5,6,7,8,9,10,11,12,13, in which valley-chiral electrons counter-propagate at the boundary in the absence of inter-valley scattering. To observe topological valley transport in the acoustic system, an experimental knob is required to trigger the valley Hall phase transition. Instead of breaking the inversion symmetry, as in graphene systems7,8,9,10,11,12,13, here we introduce a mirror-symmetry-breaking mechanism to induce the acoustic analogue of the topological semimetal–insulator transition. To do this, anisotropic scatterers, which have been employed to optimize the performance of the SCs16,17, are used to close and reopen the bandgaps18,19. Interestingly, by simply rotating the anisotropic scatterers, we can control the frequency gap of the acoustic insulators and the shape of the phase domain wall flexibly.
As depicted in Fig. 1a, the SC consists of a triangular-lattice array (with lattice constant 4.3 cm) of rod-like scatterers placed in a 2D air waveguide formed by two parallel plastic plates (Methods). Here the top plate is removed for visualization. The rod has the shape of a regular triangle (with side length 3.0 cm), and its orientation degree of freedom, characterized by the rotation angle α, enables the system to exhibit different symmetries. For α = mπ/3, with m being an integer, the point group at K (and K′) is featured by C3v symmetry, due to the perfect match of the mirrors of individual scatterers to those of the triangular-lattice. Apart from the specific angles, the mirrors are mismatched and the symmetry reduces to C3. Therefore, a two-fold Dirac degeneracy is protected at K (and K′) for any SC with α = mπ/3, whereas the degeneracy would be lifted for any other rod orientation because of the broken mirror symmetry. This is exemplified by the dispersions for the SCs with α = 0° and −10° (Fig. 1b). To reveal the characteristics of the Dirac degeneracy, Fig. 1c presents the field patterns for the two degenerate states at K. Both states, labelled by ψp−0 and ψq+0, display typical vortex profiles centred at the two inequivalent triangular-lattice centres p and q, where the sign ± indicates an anticlockwise or clockwise energy flow. The vortex chirality can be viewed as a pseudospin analogous to the A–B sublattice or the top–bottom layer pseudospin in graphene systems6,7,8. The counterparts at valley K′, denoted by ψp+0 and ψq−0, possess invariant vortex cores but opposite chirality, as required by time-reversal symmetry. Below we focus on the physics at valley K, while that at K′ can be derived directly from time-reversal symmetry.
The acoustic valley Hall (AVH) phase transition has been convincingly observed in our simulations and experiments. Figure 1d provides a continuous evolution of the band-edge frequencies ωp− and ωq+ versus the rotation angle α. Evidently, the frequency order of the first two bands, locked to the opposite vortex pseudospins, is inverted at α = 0°, a signal of the AVH phase transition. This is analogous to the semimetal–insulator transition induced by an electric field in bilayer graphene7,8. We have measured the band-edge frequencies for SCs with different α. The experimental data in Fig. 1d (circles) confirm the gap closure and reopening as the scatterer is rotated through α = 0°. To further identify different AVH insulators, the phase profiles of the K valley states are scanned around the phase singular point p or q. As exemplified by the SCs with α = ±10° (Fig. 1e), the overall slopes of the experimentally measured phase distributions (circles) agree well with those extracted from the eigenfields (lines). Intriguingly, the vortex states labelled by p− and q+ carry quantized angular momenta −1 and +1, respectively, owing to the three-fold rotation symmetry ofthe SC.
The AVH phase transition can be captured by an α-dependent continuum Hamiltonian, δH = vDδkxσx + vDδkyσy + mvD2σz, derived from the k ⋅ p perturbation method and spanned by the degenerate vortex pseudospins ψp−0 and ψq+0 (Methods). Here vD is the Dirac velocity of the conic dispersion at α = 0°, δk is the momentum deviation from the K point, and σi are Pauli matrices of the vortex pseudospins. Consistent with the above band inversion picture (Fig. 1d), the sign of the effective mass m = (ωq+ − ωp−)/2vD2 characterizes two different AVH insulators separated by the Dirac semimetal phase with m = 0 in the phase diagram. Similarly to the case of graphene, the massive Dirac Hamiltonian δH produces a nontrivial Berry curvature Ω(δk) = (1/2)mvD(δk2 + m2vD2)−3/2 (refs 2,8) in the first band, which can be integrated into a topological charge CK = sgn(m)/2 (refs 6,7,8,9,10,11). Therefore, for a SC interface separating distinct AVH insulators, the difference in the topological charge across the interface is quantized (that is, |ΔCK| = 1), which predicts a chiral edge mode propagating along the interface6,7,8,9,10,11. This bulk-boundary correspondence has been confirmed by analytically solving the boundary problem (Methods). For an interface oriented along an arbitrary x′ direction and located at y′ = 0, the edge state has the general form EquationSource math miφ mo= mrow mo( mrow msub mic mn1 msubsup miψ mrow msup mip mo− mn0 mo+ msub mic mn2 msubsup miψ mrow msup miq mo+ mn0 mo) msup mtexte mrow mii miδ msub mik msup mix mo′ msup mix mo′ mo- mo| mim msub miv mtextD msup miy mo′ mo| , with a linear dispersion δω = vDδkx′ centred at the projected K valley onto the x′ direction. Owing to time-reversal symmetry, the K′ valley also hosts one chiral AVH edge mode, but travelling along the opposite direction.
To confirm the above picture, we first simulate the dispersion for the SC interface selected, for example, along the x direction. Two different systems are studied comparatively: one is constructed by the SCs with α = 10° and 50°, and the other is constructed by the SCs with α = −10° and 10° (referred to as domains I and II below). These SCs are considered intentionally since they share the same bandgap between 3.82 and 4.34 kHz (Fig. 1d). As shown clearly in Fig. 2a, for the former system the edge spectrum is completely gapped since the two SCs belong to the same AVH phase; for the latter case, however, a pair of valley-chiral edge states counter-propagates at each interface (green lines) because of the mass inversion. The edge states φI, II± and φII, I± respectively label those gapless modes hosted by the interfaces I–II and II–I displayed in Fig. 2b, travelling along the ±x directions. The dispersions are linear near the crossing points, as predicted by the continuum model (red lines).
The presence of the topological AVH edge states has been solidly validated in experiments. As shown in Fig. 2b, a sandwich structure made of the SCs with α = −10° (phase I) and 10° (phase II) is used to study two different horizontal interfaces simultaneously. The sound signal is launched from a deep-subwavelength-sized tube placed inside the left entrance of the SC interface, and probed by a movable microphone in the same channel. From the pressure distributions scanned separately along the interfaces I–II and II–I, we obtain the corresponding interface spectra through a Fourier transform. As shown in Fig. 2c, where the dark and bright colours indicate the Fourier amplitude of low and high values, the experimental data (bright colour) capture well the numerical dispersions (green lines) for the right-moving states φI, II+ and φII, I+, even in the frequency regimes beyond the bulk gap (since the bulk mode is not well excited). The time-reversal counterparts, φI, II− and φII, I−, have also been checked by sending sound signals from the entrances on the right. Note that a key feature of the topological edge state is the exponential decay of the field amplitude away from the SC interface, for which the decay length is determined only by the bulk parameter |mvD|. To confirm this, we scan the pressure field along a straight line vertically traversing the interface II–I. As shown in Fig. 2d, the experimental result agrees excellently with the model prediction and full-wave simulation, which is indeed insensitive to frequency (see inset).
As aforementioned, the topological AVH edge state stems inherently from the single-valley physics. There are two implications: each chiral edge mode is projected from a specific valley and is also a linear superposition of the two basis states. The inter-valley decoupling can thus lead to many fascinating transport phenomena—for example, angularly selective excitation by external sound and negligible backscattering in sharply curved SC interfaces.
To stimulate the valley-projected interface state by a spatial Gaussian beam, an angular selection rule can be established according to the conservation of the momentum parallel to the sample boundary, k∥ = k0 sinγ. Here k∥ is a projection of K (or K′) on the boundary, k0 is the wavevector in air space, and γ characterizes the incidence direction of the sound beam. The angular selectivity is exemplified in Fig. 3a by a vertical SC interface. As predicted, the AVH edge mode (projected from the K valley) is well excited at the incidence angle γ = sin−1(k∥/k0) ≈ −42° (upper panel), in striking contrast to the deep suppression at γ = 0° (lower panel). To validate this phenomenon, we have experimentally measured the transmissions for a wide range of incidence angles. The angular selectivity can be observed clearly in Fig. 3b (black circles), where the optimized transmission occurs in the predicted incidence angle, in good agreement with the simulation (black line). Similarly, an optimal excitation emerges at γ = 42° if the positions of the two AVH phases are switched (red line and circles), and the edge mode is now contributed from the K′ valley. The broadening of the angular peaks stems from the presence of the exponential decay of the edge states. Remarkably, the angular selection rule is robust for the entire bulk gap, as displayed consistently by the numerical and experimental data in Fig. 3c. The horizontal SC interfaces in Fig. 2b deserve special attention. As shown in Fig. 3d, e, for an incident Gaussian beam wide enough to cover both interfaces simultaneously, only the edge state in the lower channel is excited despite the fact that the upper one is also allowed by the criterion of momentum conservation. This can be understood from the parities of the edge modes (Methods).
Below we demonstrate the negligibly weak backscattering of the AVH edge mode propagating along sharply twisted SC interfaces. Such curved interfaces emerge frequently in electronic systems12 and may cause inter-valley scattering to destroy the edge states10. Figure 4a shows the sound transport in a zigzag bending channel. As exemplified by the field pattern simulated at 4.06 kHz (inset), the sound travels smoothly in the curved path despite suffering two sharp corners (bent by 120°). To give an exact description for the bending corners, the reflection (as revealed by the interference pattern in the inlet of the waveguide) originating from the impedance mismatch between the sample and the free space must be precluded. For this purpose, we have derived the transmission and reflection through a one-dimensional scattering matrix approach (Supplementary Information). Figure 4a shows a negligible reflection in the entire bulk gap. This behaviour has not been observed previously in the sound waveguide designed by a SC14,15. Note that the accuracy of the scattering matrix method depends sensitively on the pressures detected (Supplementary Information). In practical experiments, we have measured the pressure in the output channel and compared it with the result for a sample containing a straight channel of the same length. As shown in Fig. 4b, the transmitted pressures for the two samples agree well in the frequency range of the bulk gap, in contrast to the remarkable difference beyond the gap. This confirms the weak influence of the bending corners on the wave transport of the interface modes. More complex configurations have been further checked numerically (Supplementary Information), even for those simultaneously sustaining the symmetric and antisymmetric edge modes. Again, the unusual phenomenon is closely related to the valley-projected topological origin of the edge mode. As deduced from the theoretical model, the forward-moving modes are always projected from the same (K) valley, since the relative positions of the two AVH insulators remain invariant with the propagation of sound. In other words, the field profiles match well between the edge states of adjacent channels, which contribute to the preference of high transmission. (For a general defect located in the curved or straight waveguide, the inter-valley coupling could become sizable, depending on the property of the defect.) Recently, the reflection immunity has also been observed for electromagnetic waves travelling in a twisted interface separated by different photonic topological insulators20,21. In those systems, the gaps are produced by a magneto-electric coupling (dubbed bianisotropy) equivalent to the electronic spin–orbital interaction, which is absent in our sonic systems.
It is worth noting that the rotating-scatterer mechanism enables easily tunable operation bandwidth and reconfigurable shape of the SC interface. These merits, plus the intriguing valley transport properties, absent in the conventional SC-based waveguides14,15, could be very useful in designing exceptional devices (for example, for sound signal processing). Our finding may also pave the way for exploring controllable topological phases and valley-dependent phenomena in various classical systems, which have been proved to be excellent macroscopic platforms in revealing topological properties22,23,24,25,26,27,28,29,30,31,32,33,34,35 proposed originally in electronic systems—for example, quantum Hall insulators22,23,24,25,26,27,28, topological insulators20,21,29,30,31,32,35, and topological semimetals33,34. Finally, the study provides a special insight into realizing topological insulators that require internal degrees of freedom, which is particularly important for neutral scalar sound that lacks an intrinsic polarization and is uncoupled from external fields.
Methods
Modelled Hamiltonian.
It has been proved14 that, for the SC with α = 0° the perturbation Hamiltonian δk ⋅ p, spanned by the degenerated states ψp−0 and ψq+0, yields conic dispersions centred at the hexagonal Brillouin zone corners due to the protection of the C3v symmetry, where p is a vectorial operator determined by the density distribution of the SC. As the scatterer is rotated, the mirror symmetry is broken, and the deterministic degeneracy is thus removed. This produces a Dirac mass term—namely, δij(ωi2 − ωD2), where ωi = ωp− or ωq+ is the α-dependent band-edge frequency, and ωD is the Dirac frequency of the SC with α = 0°. Using the detailed form of the p matrix—that is, p11 = p22 = 0 and —we obtain a compact form of the perturbation Hamiltonian δH (see text) that satisfies the eigen-problem δHψ = δωψ, where δω is the frequency deviation from ωD. The preciseness of the modelled Hamiltonian has been confirmed by duplicating the numerical dispersion near the Brillouin zone corner.
Derivation of the AVH edge states.
For simplicity, we consider a SC interface oriented along the x direction, formed by two SCs with m < 0 for y > 0 and m > 0 for y < 0. Substituting a trial decay solution into the eigen-problem, an edge state travelling along the +x direction can be derived EquationSource math mrow msup miφ mo+ mo= mrow mo( mrow msubsup miψ mrow msup mip mo− mn0 mo+ msubsup miψ mrow msup miq mo+ mn0 mo) msup mtexte mrow mii miδ mrow msub mik mix mix mo- mo| mim mrow msub miv mtextD miy mo| , together with a gapless dispersion δω = vDδkx centred at the projection of K onto the interface. Interestingly, because of the mirror symmetry between the two basis states (Fig. 1c), the edge state is locally symmetric with respect to a specific horizontal axis, apart from the exponential decay factor EquationSource math msup mtexte mrow mim msub miv mtextD mtexty mo= msup mtexte mrow mo− mo| mim msub miv mtextD mtexty mo| . Note that the property of the edge state is much different if the positions of the two AVH insulators are switched, since the detailed geometrical structure of the SC interface is inherently different from the original one. Now the edge state EquationSource math mrow msup miφ mo+ mo= mrow mo( mrow msubsup miψ mrow msup mip mo+ mn0 mo− msubsup miψ mrow msup miq mo− mn0 mo) msup mtexte mrow mo− mii miδ mrow msub mik mix mix mo- mo| mim mrow msub miv mtextD miy mo| , projected from the K′ valley, is locally antisymmetric in each unit cell. This leads to the deep suppression of the edge mode in the upper channel in Fig. 3d, e. Similarly, the valley-projected edge state can be derived for a SC interface with any arbitrary orientation.
Simulations.
All full-wave simulations are accurately carried out by a commercial finite-element solver (COMSOL Multiphysics), where the triangular polymethyl methacrylate rods used in real experiments are modelled as acoustically rigid, considering the great impedance mismatch with respect to air. Note that, in practical experiments, the rods (of height 1.2 cm) are closely sandwiched between two acoustically rigid parallel plates. The whole structure can be safely modelled as a 2D system, since the planar waveguide supports only the propagating mode uniform in the z-direction for the wavelength under consideration. To capture the weak reflection exactly from the two bent corners (Fig. 4), we must remove the influence of the multi-reflections induced by the impedance mismatch between the sample and free space. Thus, we collect the forward-moving and backward-moving wave information in the input and output channels, according to the Bloch wave connection between two spatially equivalent locations in each channel. From the incoming and outgoing wave information for the corners, a simple one-dimensional scattering model (Supplementary Information) can be established as long as each phase domain is thick enough to prevent energy leakage.
Experimental measurements.
In experiments, several hundreds of triangular polymethyl methacrylate rods fabricated by laser cutting are arranged into desired sample configurations. The sound signal is launched from a narrow tube (of diameter ∼0.8 cm) and scanned by a movable microphone (of diameter ∼0.7 cm, B&K Type 4187), together with an identical microphone fixed for phase reference. The sound signal is analysed by a multi-analyser system (B&K Type 3560B), from which both the phase and amplitude of the pressure field can be obtained. To prepare a Gaussian beam, the sound emitted from the tube is reflected by a carefully designed parabolic concave mirror, where the beam width is controlled by the mirror’s size. In all measurements, absorbers are placed at the ends of the sample to reduce the unwanted reflection generated by the impedance mismatch between the sample and free space.
Geometric parameters of the experimental samples are listed below. The sample (see Fig. 2b) used to measure the dispersion is formed by three parts with identical areas, where the central one (domain II) corresponds to the SC with the rotation angle α = 10°, and the lateral two (domain I) correspond to the SC with α = −10°. Each part consists of 16 × 16 rods—that is, 16 1ayers along the y direction and 16 rods for each layer along the x direction. The sample has a total size ∼0.7 × 1.8 m2. The sample used to confirm the angularly sensitive excitation of the edge modes (Fig. 3b, c) is made of 16 layers and 22 rods for each, and the sample involved in Fig. 3e is completely identical to that mentioned in Fig. 2b. In Fig. 4b, the sample with a zigzag (or straight) path consists of 22 × 40 (or 22 × 46) rods.
Data availability.
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon request.
References
Rycerz, A., Tworzydlo, J. & Beenakker, C. W. J. Valley filter and valley valve in graphene. Nat. Phys. 3, 172–175 (2007).
Xiao, D., Yao, W. & Niu, Q. Valley-contrasting physics in graphene: magnetic moment and topological transport. Phys. Rev. Lett. 99, 236809 (2007).
Gorbachev, R. V. et al. Detecting topological currents in graphene superlattices. Science 346, 448–451 (2014).
Xu, X., Yao, W., Xiao, D. & Heinz, T. F. Spin and pseudospins in layered transition metal dichalcogenides. Nat. Phys. 10, 343–350 (2014).
Mak, K. F., McGill, K. L., Park, J. & McEuen, P. L. The valley Hall effect in MoS2 transistors. Science 344, 1489–1492 (2014).
Semenoff, G. W., Semenoff, V. & Zhou, F. Domain walls in gapped graphene. Phys. Rev. Lett. 101, 087204 (2008).
Martin, I., Blanter, Y. M. & Morpurgo, A. F. Topological confinement in bilayer graphene. Phys. Rev. Lett. 100, 036804 (2008).
Zhang, F., Jung, J., Fiete, G. A., Niu, Q. & MacDonald, A. H. Spontaneous quantum Hall states in chirally stacked few-layer graphene systems. Phys. Rev. Lett. 106, 156801 (2011).
Qiao, Z., Jung, J., Niu, Q. & Macdonald, A. H. Electronic highways in bilayer graphene. Nano Lett. 11, 3453–3459 (2011).
Zhang, F., MacDonald, A. H. & Mele, E. J. Valley Chern numbers and boundary modes in gapped bilayer graphene. Proc. Natl Acad. Sci. USA 110, 10546–10551 (2013).
Vaezi, A., Liang, Y., Ngai, D. H., Yang, L. & Kim, E.-A. Topological edge states at a tilt boundary in gated multilayer graphene. Phys. Rev. X 3, 021018 (2013).
Ju, L. et al. Topological valley transport at bilayer graphene domain walls. Nature 520, 650–655 (2015).
Li, J. et al. Gate-controlled topological conducting channels in bilayer graphene. Nat. Nanotech. 11, 1060–1065 (2016).
Khelif, A., Choujaa, A., Benchabane, S., Djafari-Rouhani, B. & Laude, V. Guiding and bending of acoustic waves in highly confined phononic crystal waveguides. Appl. Phys. Lett. 84, 4400 (2004).
Khelif, A. & Adibi, A. Phononic Crystals: Fundamentals and Applications (Springer, 2016).
Goffaux, C. & Vigneron, J. P. Theoretical study of a tunable phononic band gap system. Phys. Rev. B 64, 075118 (2001).
Moiseyenko, R. P., Liu, J., Declercq, N. F. & Laude, V. Blazed phononic crystal grating. Appl. Phys. Lett. 102, 034108 (2013).
Lu, J. et al. Dirac cones in two-dimensional artificial crystals for classical waves. Phys. Rev. B 89, 134302 (2014).
Lu, J., Qiu, C., Ke, M. & Liu, Z. Valley vortex states in sonic crystals. Phys. Rev. Lett. 116, 093901 (2016).
Ma, T., Khanikaev, A. B., Mousavi, S. H. & Shvets, G. Guiding electromagnetic waves around sharp corners: topologically protected photonic transport in metawaveguides. Phys. Rev. Lett. 114, 127401 (2015).
Cheng, X. et al. Robust reconfigurable electromagnetic pathways within a photonic topological insulator. Nat. Mater. 15, 542–548 (2016).
Haldane, F. D. M. & Raghu, S. Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry. Phys. Rev. Lett. 100, 013904 (2008).
Wang, Z., Chong, Y., Joannopoulos, J. D. & Soljačić, M. Observation of unidirectional backscattering-immune topological electromagnetic states. Nature 461, 772–775 (2009).
Fang, K., Yu, Z. & Fan, S. Realizing effective magnetic field for photons by controlling the phase of dynamic modulation. Nat. Photon. 6, 782–787 (2012).
Rechtsman, M. C. et al. Photonic Floquet topological insulators. Nature 496, 196–200 (2013).
Yang, Z. et al. Topological acoustics. Phys. Rev. Lett. 114, 114301 (2015).
Khanikaev, A. B., Fleury, R., Mousavi, S. H. & Alù, A. Topologically robust sound propagation in an angular-momentum-biased graphene-like resonator lattice. Nat. Commun. 6, 8260 (2015).
Ni, X. et al. Topologically protected one-way edge mode in networks of acoustic resonators with circulating air flow. New J. Phys. 17, 053016 (2015).
Hafezi, M., Demler, E. A., Lukin, M. D. & Taylor, J. M. Robust optical delay lines with topological protection. Nat. Phys. 7, 907–912 (2011).
Khanikaev, A. B. et al. Photonic topological insulators. Nat. Mater. 12, 233–239 (2013).
Chen, W. J. et al. Experimental realization of photonic topological insulator in a uniaxial metacrystal waveguide. Nat. Commun. 5, 6782 (2014).
Süsstrunk, R. & Huber, S. D. Observation of phononic helical edge states in a mechanical topological insulator. Science 349, 47–50 (2015).
Lu, L. et al. Experimental observation of Weyl points. Science 349, 622–624 (2015).
Xiao, M., Chen, W., He, W. & Chan, C. T. Synthetic gauge flux and Weyl points in acoustic systems. Nat. Phys. 11, 920–924 (2015).
Lu, L. et al. Symmetry-protected topological photonic crystal in three dimensions. Nat. Phys. 12, 337–340 (2016).
Acknowledgements
We thank Chuanwei Zhang and Meng Xiao for fruitful discussions. This work is supported by the National Basic Research Program of China (Grant No. 2015CB755500); National Natural Science Foundation of China (Grant Nos. 11674250, 11374233, 11534013, 11574233, 11547310); Postdoctoral Innovation Talent Support Programme (BX201600054). F.Z. is supported by the UT-Dallas research enhancement funds. F.Z. is also grateful to the Kavli Institute for Theoretical Physics for their hospitality during the finalization of this work, which is supported in part by the National Science Foundation under Grant No. PHY11-25915.
Author information
Authors and Affiliations
Contributions
C.Q. and J.L. conceived the original idea. J.L. performed the simulations. J.L., L.Y., X.F. and M.K. carried out the experiments. C.Q. and Z.L. supervised the project. C.Q., J.L., F.Z. and Z.L. analysed the data and wrote the manuscript. All authors contributed to scientific discussions of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary information
Supplementary information (PDF 907 kb)
Rights and permissions
About this article
Cite this article
Lu, J., Qiu, C., Ye, L. et al. Observation of topological valley transport of sound in sonic crystals. Nature Phys 13, 369–374 (2017). https://doi.org/10.1038/nphys3999
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys3999
This article is cited by
-
Elastic three-dimensional phononic topological insulators with Dirac hierarchy
Science China Physics, Mechanics & Astronomy (2024)
-
Topology optimization of acoustic bandgap crystals for topological insulators
Engineering with Computers (2024)
-
Topological pseudospin switcher: Mechanism and acoustic realization
Science China Physics, Mechanics & Astronomy (2024)
-
Breakdown of conventional winding number calculation in one-dimensional lattices with interactions beyond nearest neighbors
Communications Physics (2023)
-
Orbital topological edge states and phase transitions in one-dimensional acoustic resonator chains
Nature Communications (2023)