
Bound states in the continuum in circular clusters of scatterers
Marc Martí-Sabaté,1Bahram Djafari-Rouhani,2and Dani Torrent1, ∗
1GROC, UJI, Institut de Noves Tecnologies de la Imatge (INIT),
Universitat Jaume I, 12071, Castelló de la Plana, Spain
2IEMN, University of Lille, Cité scientifique, 59650 Villeneuve d’Ascq, France
(Dated: October 19, 2022)
In this work, we study the localization of flexural waves in highly symmetric clusters of scatterers. It is
shown that when the scatterers are placed regularly in the perimeter of a circumference the quality factor of the
resonances strongly increases with the number of scatterers in the cluster. It is also found that in the continuous
limit, that is to say, when the number of scatterers tends to infinite, the quality factor is infinite so that the modes
belong to the class of the so called bound states in the continuum or BICs, and an analytical expression for the
resonant frequency is found. These modes have different multipolar symmetries, and we show that for high
multipolar orders the modes tend to localize at the border of the circumference, forming therefore a whishpering
gallery mode with an extraordinarily high quality factor. Numerical experiments are performed to check the
robustness of these modes under different types of disorder and also to study their excitation from the far field.
Although we have focused our study to flexural waves, the methodology presented in this work can be applied
to other classical waves, like electromagnetic or acoustic waves, being therefore a promissing approach for the
design of high quality resonators based on finite clusters of scatterers.
I. INTRODUCTION
Bound states in the continuum (BICs) are eigenmodes of a
system whose energy lies in the radiation part of the spectrum
while remaining localized in a finite part of the system and
with an infinite lifetime. These states were first mathemati-
cally proposed in 1929 by von Neumann and Wigner in the
framework of quantum mechanics [1], although the concept
has been extended to classical waves, like acoustics [2–6], mi-
crowaves [7, 8] or optics [9–11].
Despite the fact that the practical realization of BICs is a
challenging problem, structures based on them present sharp
resonances with extremely high quality factors, which have
as well the advantage, unlike ideal BICs, that can be ex-
cited with external radiative fields. Also named quasi-BICs
(or QBICs), these modes have been widely used in sensing
applications[12–14].
Among the wide variety of geometries and structures used
to find BICs[15], those based on finite structures are specially
interesting for practical applications, since periodic or waveg-
uide BICs will always present finite-size effects which will de-
crease their efficiency. For instance, circular clusters of scat-
terers studied in some recent works[14, 16] are extraordinarily
convenient from the practical point of view. In this work, we
will generalize the study of these circular clusters of scatter-
ers to provide a general schema for the realization of QBICs
based on this geometry.
The manuscript is organized as follows: After this intro-
duction, in section II we study the formation of bound states
in the continuum in open systems by attaching a cluster of
mass-spring resonators to a thin elastic plate. We will find that
when the scatterers in the cluster are arranged in the corners of
a regular polygon the quality factor of the resonances quickly
increases with the number of scatterers in the cluster. In sec-
tion III we perform several numerical experiments to check
∗dtorrent@uji.es
the robustness of these modes, and in section IV their excita-
tion from the far field will be considered. Finally, section V
sumarizes the work.
II. EIGENMODES OF A POLYGONAL CLUSTER OF
SCATTERERS
The propagation of flexural waves in thin elastic plates
where a cluster of Npoint-like resonators has been attached
at positions Rαis described by means of the inhomogeneous
Kirchhoff[17] equation
(∇4−k4
0)ψ(r) =
N
∑
α=1
tαδ(r−Rα)ψ(r)(1)
where ψ(r)is the spatial part of the vertical displacement of
the plate, which is assumed to be harmonic and of the form
W(r,t) = ψ(r)e−iωt.(2)
Also, the free space wavenumber k0is given by
k4
0=ρh
Dω2,(3)
with ρ,hand Dbeing the plate’s mass density, height and
rigidity, respectively. The response of each resonator is given
by the tαcoefficient, which is a resonant quantity whose prop-
erties depend on the geometry of the scatterer attached to the
plate[18]. However, for the purposes of the present work, it
will be assumed that it can take any real value in the range
tα∈(−∞,∞).
A self-consistent multiple scattering solution can be found
for the above equation as
ψ(r) = ψ0(r) +
N
∑
α=1
BαG(r−Rα)(4)
arXiv:2210.09833v1 [physics.class-ph] 17 Oct 2022