Ultrasonic nodal chains in topological granular metamaterials

TL;DR

This paper introduces a 3D granular metamaterial that exhibits ultrasonic nodal chains, revealing new classical topological phases with robust surface states suitable for ultrasonic testing and material analysis.

Contribution

The study demonstrates a novel 3D mechanical metamaterial with multiple intersecting nodal rings, advancing topological phases in classical systems and overcoming challenges faced in electronic topological materials.

Findings

Unveiled ultrasonic nodal chains in granular metamaterials

Surface states are unaffected by crystal termination

Potential applications in ultrasonic non-destructive testing

Abstract

Three-dimensional (3D) Weyl and Dirac semimetals garner considerable attention in condensed matter physics due to the exploration of entirely new topological phases and related unconventional surface states. Nodal line and ring semimetals on the other hand can facilitate 3D band crossings characterized by nontrivial links such as coupled chains and knots that are protected by the underlying crystal symmetry. Experimental complexities, detrimental effects of the spin-orbit interaction, along with the merging of the underlying surface states into the bulk pose great challenges for growing advancements, but fortunately enable other systems, such as bosonic lattices, as ideal platforms to overcome these obstacles. Here we demonstrate a 3D mechanical metamaterial made of granular beads, which is predicted to provide multiple intersecting nodal rings in the ultrasonic regime. By unveiling these yet unseen classical topological phases, we discuss the resilience of the associated novel surface states that appear entirely unaffected to the type of crystal termination, making them a superb platform in ultrasonic devices for non-destructive testing and material characterization.

Equations

$$m_b\frac{{\partial }^2{\text{u}}^{\alpha }}{\partial t^2}=\sum _{\beta }{\text{F}}^{\beta \alpha },$$

$$I_b\frac{{\partial }^2{\text{w}}^{\alpha }}{\partial t^2}=\sum _{\beta }{\text{M}}^{\beta \alpha }+\frac{1}{2}\sum _{\beta }({\text{R}}^{\beta }-{\text{R}}^{\alpha })\times {\text{F}}^{\beta \alpha },$$

Full text

Topological phases of matter in insulators and superconductors have recently been extended to semimetals, thus broadening the family of exotic topological states Qi and Zhang (2011); Armitage et al. (2018). This active frontier in condensed matter physics recently explored topologically protected degeneracies in Dirac and Weyl semimetals that are identified by topologically robust band-touching manifolds. These unconventional fermionic semimetals are characterized by nontrivial band touching in the form of zero-dimensional discrete points, one-dimensional nodal rings and lines, or two-dimensional nodal surfaces. Peculiar topologically protected k-space geometrical manifolds appear in various interlaced nodal chains and links whose nontrivial linking carry a toroidal Berry phase of $\pi$, provided the loops enclose the aforementioned nodal shapes Burkov et al. (2011); Bzdusek et al. (2016).

The nontrivial Berry flux is responsible for the formation of topologically protected surface states that emerge from the intersecting points of nodal shapes. E.g., bulk Weyl points give rise to surface Fermi arcs, whereas nodal rings host flat drumhead surface states that are relevant for topologically robust transport. Various materials have been predicted to sustain Dirac or Weyl fermions and contemporary experiments report on the observation of anisotropic or negative magnetoresistance in three dimensional semimetals Huang et al. (2015); Arnold et al. (2016); Schoop et al. (2016a); Wang et al. (2016); Hong et al. (2018).

As opposed to nodal points and lines, nodal rings can generate exotic intersecting formations in the form of chains, knots and Hopf links that opens new horizons for unprecedented topological properties, unusual surface states and novel physics in general. A vast amount of theory has been devoted to this frontier in Heusler compounds, alkaline earth metals and electrides Schoop et al. (2016b); Chen et al. (2017); Yan et al. (2017); Bi et al. (2017); Chang et al. (2017); Gong et al. (2018); Hirayama et al. (2018), whereupon experimental verifications of nodal rings have already been observed in non-centrosymmetric superconducting compounds and Zirconium based structures Bian et al. (2016); Neupane et al. (2016); Hu et al. (2016); Lou et al. (2018).

Bosonic settings such as photonic, phononic and sonic crystals are widely studied with the aim to tailor classical wave properties. Accordingly, given the complexity of electronic systems to unveil Weyl and Dirac semimetal physics as well as those peculiar intersecting nodal rings, man-made bosonic system have become increasingly popular to explore in the waves based context, the existence of Dirac-like plasmons Weick et al. (2013), optical Weyl points and Fermi arcs Lu et al. (2015); Noh et al. (2017), hourglass nodal lines Xia et al. (2019), metallic mesh nodal chains Yan et al. (2018) and acoustic Weyl points Xiao et al. (2015); Li et al. (2018) including topologically protected negative refraction He et al. (2018).

In this work, we present an excerpt of the catalogue of nodal-chain semimetals for classical mechanical waves in man-made granular metamaterials. By solving the equations of motion for these artificial discrete media, we are able to engineer exotic mechanical intersecting nodal rings whose surface states at the truncated lattice interfaces display a remarkable durability. The elastic behavior of granular media drastically differs from conventional elastic waves in solids. At frequencies below the first spheroidal resonance of individual beads, the structure is modeled by rigid masses of finite size connected by stiffnesses originating from the contact laws between the beads, thus forming a discrete crystalline lattice. Specifically here, we consider a face centered cubic (FCC) arrangement as illustrated in Fig. 1(a). In the general case in the linear approximation, the contact between two grains is described by one shear stiffness $K_{S}$, one torsion stiffness $G_{T}$ and one bending stiffness $G_{B}$ accounting for the sliding, twisting and rolling resistances at the level of the contact, respectively, in addition to the usual normal stiffness $K_{N}$ as sketched in Fig. 1(a). As a consequence, the rotational degrees of freedom of each individual particle play an important role in the dynamics of granular media. The linear equations of motion for one bead $\alpha$, with its infinitesimal displacement $\textbf{u}^{\alpha}$ and angular rotation $\textbf{w}^{\beta}$ around its equilibrium position in a monodisperse granular assembly, read for translation

[TABLE]

and for rotation

[TABLE]

where the summation of forces $\textbf{F}^{\beta\alpha}$ and torques $\textbf{M}^{\beta\alpha}$ take place over all the beads $\beta$ at the position $\textbf{R}^{\beta}$ in contact with the bead $\alpha$ of mass $m_{b}$, moment of inertia $I_{b}$, and at the position $\textbf{R}^{\alpha}$ Merkel et al. (2010); Merkel and Luding (2017); sup . The FCC structure is formed by spheres of radius $r_{b}=5$ mm made of stainless steel with a Young modulus $E=200$ GPa, Poisson’s ratio $\nu=0.3$ and density $\rho_{b}=7.7\cdot 10^{3}$ kg.m*-3*. The lattice constant of the FCC structure is $a=(2\sqrt{2})r_{b}$. We assume that the contacts between the particles are formed by a solid bridge with a radius $r_{s}=r_{b}/20$ and a length $h_{s}=r_{b}/50$. The different stiffnesses are $K_{N}=\pi Er_{s}^{2}/h_{s}=3.9\cdot 10^{8}$ N/m, $K_{S}=K_{N}/(2+2\nu)=1.5\cdot 10^{8}$ N/m, $G_{T}=K_{S}r_{s}^{2}/2=4.7$ Nm/rad and $G_{B}=K_{N}r_{s}^{2}/4=6.1$ Nm/rad Brendel et al. (2011). The Bloch bands structure shown in Fig. 1(b) depict the propagation of a longitudinal (L), rotational (R), two Transverse-Rotational (TR) and two Rotational-Transverse (RT) modes Merkel et al. (2011). The influence of the rotational degrees of freedom has been experimentally confirmed in granular assembly with both millimeter-sized spheres Merkel et al. (2011) and micron-sized spheres Hiraiwa et al. (2016) showing the wide range of frequencies and size scales where this model remains valid. The RT and R modes are optical-type modes and have no equivalent in continuum elastic solids.

Thanks to the existence of these two type of modes, the bulk bands display accidental degeneracies as marked by squares, triangles, diamonds and hexagons in Fig. 1(b) and Fig. 1(d) in the highlighted zone around 150 kHz. These degeneracies are accidental because of their dependence on the values of the stiffnesses. Owing to the symmetries of the lattice that prohibit a gap opening, the band crossings are not discrete in the reciprocal space but extend to nodal rings in the Brillouin zone, as shown in Fig. 1(d). The formation of the intersecting nodal rings stems from three types of symmetries, namely the mirror reflection symmetry, the time-reversal along with the space inversion symmetries and the nonsymmorphic group with glide plane or screw axes symmetry Yu et al. (2017). On account of the $O_{h}$ point group symmetry of the FCC lattice, the (001) plane, comprising either the $k_{x}$, $k_{y}$ or $k_{z}=0$ planes, is a mirror reflection plane and thus hosts here two overlaying nodal rings with frequencies ranging from 148 to 152 kHz. In the frequency range spanning from 147 to 157 kHz, two nodal rings lie in and are protected by the (110) mirror plane, here shown as the $k_{y}=k_{z}$ plane. The high-symmetry direction $\Gamma$–X joins the (110) and (001) surfaces, therefore the aforementioned inner nodal chains are interconnected in the form of outer nodal chains at the point marked with an orange square as shown in Fig. 1(d) Chang et al. (2017). Similarly, the high symmetry direction $\Gamma$–Z also connects these two symmetry planes, thus, these two nodal chains merge again at the point marked with a triangle. The (111) plane contains the space inversion symmetry along with the time reversal symmetry that is always preserved in our case, consequently, two concentric nodal rings are protected across this surface, however, slightly off-plane Hirayama et al. (2017). Since the bands crossings that are marked with orange diamonds and hexagons occur across the L–U and L–K directions that extend onto the (110) plane, the four nodal rings intersect in an interlaced combination of outer nodal chains and a Hopf link.

We now investigate the existence of topologically protected surface states on the (001) plane by means of numerical simulations in a finite FCC granular metamaterial slab, which remains periodic only in the $x$ and $y$ directions but has stress-free boundary conditions at the upper and lower crystal terminations. The numerical simulations are carried out using a Discrete Particle Method with the code MercuryDPM Thornton et al. (2012); Weinhart et al. (2012); mer and the dispersion bands are found by computing the two-dimensional Fourier transform in time and space along one unique direction Mouraille et al. (2006); Merkel and Luding (2017); sup . In order to unravel these geometrically protected surface excitations in the entire 3D space, we compute them along the usual in-plane components of the wavevector but evaluate the band diagram for various out-of-plane wavevector components $k_{z}$. Technically speaking, we implemented this by an emitting bead-array whose phase difference accounts for different wave propagation momenta. The band diagram of the granular metamaterial slab, which is presented in Fig. 2 is characterized by red circles, also seen in Fig. 1(d), indicating the dispersion along the rim of the nodal ring at $k_{y}=0$. The gray shaded background of the diagrams illustrate the projected bulk bands sustained by the slab containing the stress-free surfaces along the $z$ direction. With seemingly little dispersion with varying out-of-plane momentum $k_{z}$, Fig. 2 illustrates the presence of a surface state within the incomplete band gaps that are confined to the nodal ring. Specifically, as rendered by the white dashed area, the bulk band gap spans from 138 to 146 kHz in between which, an entirely flat surface state resides (141 kHz). Further, this surface excitation that transpires from the degenerate bulk states (red circles) with a near zero group velocity, emerges from the non trivial Berry phase of the nodal rings. In addition, as expected, at lower frequencies trivial surface states that are associated to the stress-free boundary condition are also excited as seen in Fig. 2.

In order to shed light on the topological properties of the surface waves discussed above, we compare three different crystal terminations. In addition to the free boundary condition, we consider the cases where the top layer is in contact with a plane rigid wall and where the top layer is in contact with a layer of fixed (immovable) particles placed in the continuity of the crystalline structure. These three crystal terminations produce an interface between the granular metamaterial and a medium into which sound waves cannot leak. A single chain of excited particles parallel to the $y$ axis acts as the source, hence, all kind of states can be easily excited parallel to the crystal termination for all possible values of $k_{z}$. In the numerical experiments, we probe the mechanical response to the excitation in the nearest vicinity of the termination. As one can see in Fig. 3 within the dashed blue ellipses, the trivial surface states at lower frequencies are highly sensitive to the specific boundary condition employed. Beyond significant shifts and enforced dispersion, the trivial surface state ceases to exist in the presence of a fixed crystal termination. On the contrary, the surface state emerging from the nodal rings remains intact and within the bulk band gap in every scenario as highlighted by the black dotted ellipses, underlining its topological origin and resilience against drastic interface perturbations.

The above mentioned robustness against interface perturbation of the topological surface states emanating the nodal rings is visualized by computing the mechanical motions as seen in Fig. 4. At the ante-penultimate layer we place a source array with a predefined spatial phase profile, as rendered in Fig. 4(a), to momentum match ($k_{x}=1.9\pi/a$, where $a$ is the lattice constant) the two distinct surface states computed in Fig. 3 at their corresponding frequencies. In Fig. 4(b) and Fig. 4(d) we compute the spatial mechanical field of the surface states confined along the stress free crystal termination. The topological non-trivial (trivial) state at 141 kHz (72 kHz) displays a strong field confinement, it appears however that the low frequency excitation displays a smaller penetration depth. Finally, we introduce a fixed termination above the top layer. Apart from an unremarkable frequency shift, the mode shape of the confined topological non-trivial surface state at 145 kHz persists virtually unaffected when comparing Fig. 4(b) and Fig. 4(c). With a fixed termination, as computed in Fig. 3, the trivial state ceases to exist, hence at 72 kHz the field remains solely localized around the excitation source as illustrated in Fig. 4(e).

In conclusion, three-dimensional man-made FCC granular metamaterials have been found to host complex structures of nodal rings interlaced in the form of inner and outer chains. The symmetry protected nodal features give rise to remarkably robust surface states that appear entirely unaffected by the type of crystal termination. Hence, our results provide a convincing platform to implement topological non-trivial nodal rings and chains physics for mechanical waves, potentially useful for robust ultrasonic sensing and transduction.

I Acknowledgements

J. C. acknowledges the support from the European Research Council (ERC) through the Starting Grant No. 714577 PHONOMETA and from the MINECO through a Ramón y Cajal grant (Grant No. RYC-2015-17156).