# Coupled acoustoplasmonic resonators: the role of geometrical symmetries

Beatriz Castillo López de Larrinzar, Jorge M García, Chushuang Xiang, Norberto D Lanzillotti-Kimura, Antonio García-Martín

► To cite this version:

Beatriz Castillo López de Larrinzar, Jorge M García, Chushuang Xiang, Norberto D Lanzillotti-Kimura, Antonio García-Martín. Coupled acoustoplasmonic resonators: the role of geometrical symmetries. Nanophotonics, 2025, 14 (1), pp.75-80. 10.1515/nanoph-2024-0482. hal-05355344

HAL Id: hal-05355344

https://hal.science/hal-05355344v1

Submitted on 8 Nov 2025

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>

Research Article

Beatriz Castillo López de Larrinzar, Jorge M. García, Chushuang Xiang, Norberto D. Lanzillotti-Kimura and Antonio García-Martín*

# Coupled acoustoplasmonic resonators: the role of geometrical symmetries

https://doi.org/10.1515/nanoph-2024-0482

Received September 13, 2024; accepted November 28, 2024;

published online January 3, 2025

**Abstract:** Acoustoplasmonic resonators, such as nanobars and crosses, are efficient acousto-optical transducers. The excitation of mechanical modes in these structures strongly depends on the spatial profile of the eigenmodes of the resonator. Using a system of two identical gold elongated bars placed on a silicon dioxide substrate, we examine how breaking mirror symmetries affects the optical and acoustic properties to provide insights in the design of acoustoplasmonic metasurfaces for nonsymmetric acousto-optical transducers. Our findings show that the absence of mirror symmetries affects differently the optical and nanomechanical response. Broken mirror symmetries not only couple nanomechanical modes existing in single bars but introduce new torsional resonant modes.

**Keywords:** optics; nanoacustics; transducers; torsional modes; symmetry breaking

## 1 Introduction

Symmetries are paramount to understand the physical properties of a given system [1]. These symmetries span to all magnitudes, as they impose restrictions, e.g., time reversal invariance makes the scattering matrix symmetric, or the absence of absorption makes it unitary. Others are

"simply" related to the geometrical aspects of the nanostructures. In many cases, the physics is determined by the presence of symmetries; in many others, the interest lies in breaking them. A beautiful example is provided by chiral systems [2]–[4].

In optics, simple geometries, when brought into close contact, might reveal scattering processes that are highly nontrivial due to the break of symmetry, e.g., chiral behavior of nonchiral interacting elements [5], [6].

In nanoacoustics [7], most of the reported results are based on simple structures or metasurfaces, which present a high degree of spatial symmetry. For example, bars, crosses, arrays, micropillars, and microdisks all present mirror symmetry planes [8]–[13]. Even the engineering of more complex devices, such as directional acoustic nanoantennas, presents these symmetry planes [9], [14].

It is expected that breaking these symmetry planes would open a plethora of different phenomena to tailor the fundamental modes in nanoacoustic resonators. For example, two concatenated superlattices presenting inverted symmetries [15] give rise to a topologically protected nanoacoustic mode. However, despite recent efforts [6], how broken symmetries enable nanoacoustic modes in three-dimensional nanostructures remain an open question.

In this work, we make use of numerical simulations to theoretically explore the idea of breaking spatial symmetries to engineer new nanoacoustic modes. Through a simple geometry, consisting of two identical coplanar elongated Au bars, we develop a comprehensive toolbox for the design of acoustoplasmonic metasurfaces for nonsymmetric hypersound transducers. To that end, we analyze separately the optical and mechanical behavior. First, we consider the optical response to establish the coupling characteristics of the system and how symmetries are lifted. Afterward, we study the nanoacoustic resonances to unveil the phonon dynamics in the case of systems with no planemirror symmetries. Then, we make analogies whenever the physical picture is similar in both systems and point out the main discrepancy points.

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

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

<footnote>*Corresponding author: Antonio García-Martín, Instituto de Micro y Nanotecnología IMN-CNM, CSIC, CEI UAM-CSIC, Tres Cantos, Spain, E-mail: a.garcia.martin@csic.es. https://orcid.org/0000-0002-3248-2708

Beatriz Castillo López de Larrinzar and Jorge M. García, Instituto de Micro y Nanotecnología IMN-CNM, CSIC, CEI UAM-CSIC, Tres Cantos, Spain.

https://orcid.org/0000-0002-3346-0025 (B. Castillo López de Larrinzar)

Chushuang Xiang and Norberto D. Lanzillotti-Kimura, Université Paris-Saclay, CNRS, Centre de Nanosciences et de Nanotechnologies, Paris, France. https://orcid.org/0000-0002-6056-5551 (N.D. Lanzillotti-Kimura)</footnote>
## 2 Description of the system and methodology

The system under study consists of two gold bars with dimensions $L = 100 \text{ nm} \times W = 30 \text{ nm} \times H = 30 \text{ nm}$ (length $\times$ width $\times$ height, as the one depicted in Figure 1(a)). To analyze the effects of symmetries, we explore the interactions of two coplanar bars attached to a $\text{SiO}_2$ substrate when they are separated $S = 20 \text{ nm}$ along the long-axis direction, Figure 1(b), or short-axis direction, Figure 1(c). These two structures present clear mirror symmetries with respect to the planes depicted as doted-dashed lines in Figure 1(b) and (c). A simple displacement $D_y$ along the long axis (short axis) for the bars aligned along the short axis (long axis) will break the mirror symmetry, leaving only a rotational symmetry of 180°, also known as C2-symmetry.

The Lumerical® software suite has been utilized to address the scattering and absorption cross sections. The computational solutions are based on the Finite Difference

Time Domain (FDTD) method, with the dispersive optical constants for gold referenced in Ref. [16]. A linearly polarized plane wave – incident along either the long or short axis of the nanostructure – with a uniform amplitude across the entire simulation cell has been employed. The simulation cell dimensions ($6 \ \mu\text{m} \times 6 \ \mu\text{m} \times 10 \ \mu\text{m}$) were chosen to ensure that perfectly absorbing boundary conditions exert a negligible influence on the resultant electromagnetic fields. A highly refined mesh has been implemented within the nanostructures and the proximal near-field region ($2.2 \ \mu\text{m} \times 2.2 \ \mu\text{m} \times 3.7 \ \mu\text{m}$) with spatial resolutions of $1 \text{ nm} \times 1 \text{ nm} \times 1 \text{ nm}$ ($dx-dy-dz$, respectively), which continuously expands to a maximum size of 10 nm as the distance from the near-field region increases toward the simulation boundaries, thus ensuring numerical convergence to the highest possible accuracy.

The Finite Element Method (FEM) in COMSOL® (Structural Mechanics Module) offers an effective framework to simulate the elastic response resulting from the isotropic thermal expansion in a metallic nanostructure [14]. The FEM simulation utilizes material parameters from COMSOL®'s materials library. The nanobars are meshed with a maximum element size of 5 nm. The $\text{SiO}_2$ substrate is a 1,600 nm diameter semi-spherical region, surrounded by nonreflective domains (Perfectly Matched Layers, PML) that truncate the physical domain. The substrate is attached to the nanobars and oscillates upon the excitation of the bars, causing the mutual interactions between the nanobars (see Supplementary Material for details). To adequately represent the acoustic interactions between bars, a 300 nm diameter semi-spherical region of the substrate beneath the gold bars is meshed with a minimum element size of 4 nm and a maximum element size of 50 nm.

The thermoelastic excitation of phonons that would be produced by photon absorption is simulated by raising the temperature of a nanobar by a 5K initial thermal excitation, which is depicted by a dark-orange illumination around the heated bar in Figure 1. The center of mass of each system, whether a single bar or double bars, is positioned at the origin of the coordinates. In studies involving two coupled bars, the geometric arrangement is designed to permit interactions (such as multiple wave scattering), in the optical and acoustical descriptions. We would like to emphasize that the optical and mechanical studies are carried out separately. In order to allow counter-phase oscillations of the coupled bars, thermal excitation is applied only to one of the bars. Phonon detection is modeled by computing the integral of the root mean square (RMS) of the displacement on the two bars.

<figure><img src="image_5.png" /><figcaption>**Figure 1:** The basic system under study in this work consists on a Au bar (a) with dimensions $L \times W \times H = 100 \times 30 \times 30 \text{ nm}^3$. The interaction of two parallel coplanar bars on a silicon oxide, $\text{SiO}_2$, substrate is studied when they are spaced $S = 20 \text{ nm}$ along the long-axis direction (b) or the short-axis direction (c and d). Additionally, a displacement $D_y = 75 \text{ nm}$ along the long-axis direction is shown in (d), which breaks the mirror-symmetry of the system. A dark-orange illumination around the heated bar represents an initial 5K thermal excitation.</figcaption></figure>
# 3 Results

## 3.1 Optical response

The optical properties of metallic bars or nanoantennae have been the subject of extensive research during the past decades. Recently, studies have focused on systems where two or more bars interact, leading to new optical phenomena [6], [17], [18]. Despite the vast literature on plasmonics in interacting bars [6], [17], [19]–[21], let us briefly describe the optical response of the structures depicted in Figure 1.

In Figure 2, we present the scattering cross section as a function of the wavelength of the incident light for the structures depicted in Figure 1 with $D_y = 75$ nm. In Figure 2(a), we present the results when the incident light is polarized along the short axis of the bars, and in Figure 2(b), when the polarization is along the long axis. The spectral location of the plasmon resonance for a single bar (or nanoparticle) is known to depend on the size of the bar aligned with the polarization direction, bigger size leads to longer resonance wavelength. That is nicely shown by the black line, which presents a resonance at ca. 550 nm for the polarization aligned along the short axis and at ca. 750 nm when aligned along the long axis. This could be understood as one particle with an electric dipolar polarization with two resonances one for each principal axis. When two resonant particles are in interacting distance this picture is modified, and resembling molecular coupling, or coupled classical springs, the two resonances give rise to a coupled system whose resonances essentially depend on the separation [19], [22]. In

our case, only one of the resulting hybridized resonances lie in the visible range of the spectrum, for both polarizations and separation distances along the short and long axis. For the two highly symmetric cases in Figure 1(b) and (c), where mirror symmetry is present, the interaction gives rise to red or blue shifts of the resonances with respect to the isolated bar (see red and blue lines in Figure 2), but the equivalent hybridized electric dipoles would be aligned along the principal axis. However, when the mirror symmetry is broken, Figure 1(d), the equivalent hybridized dipoles would be misaligned with respect to the principal axis of the individual bars. That is equivalent to having a rotated dipole, and thus both resonances would be excited for both polarizations, along the short and long axis, but with different strengths. This is exactly the situation displayed in Figure 2(a), where it can be seen that the resonance along the long axis is excited when the incident polarization is along the short one. It must be noted that a multipolar expansion of the induced currents in the bars shows that (see Supplementary Material), the resonance corresponding to the short axis is also excited for polarization along the long axis, but the intensity is so low that it can barely be appreciated in Figure 2(b).

## 3.2 Nanomechanical response

In order to unveil the physics introduced by the absence of a mirror symmetry, it is convenient to first understand the resonances associated to the most elementary systems presenting clear symmetry planes. When considering two

<figure><img src="image_6.png" /><figcaption>**Figure 2:** Optical cross sections. (a) Scattering cross section for polarization along the short axis, depicting the single resonance at ca. 550 nm for all the mirror-symmetric cases (with the corresponding red and blue shifts due to interactions), and the simultaneous excitation of the long-axis resonance for the shifted case due to the broken mirror symmetry. (b) The same for polarization along the long axis. The resonance is ca. 750 nm and the blue and red shifts are more pronounced, but the weaker short axis resonance is within the numerical noise level.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 3:** Integral of the RMS displacement as a function of frequency for the structures in Figure 1(a)–(c). The insets depict the dynamics of the vibrational modes for the low-frequency region. Bars aligned along the short (long) axis show a vibrational mode schematically depicted by black arrows that correspond to a frequency peak marked with a red (green) arrow.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4:** Integral of the RMS displacement as a function of frequency for the structures in Figure 1(d) for $D_y = 75$ nm. The peak around 4.6 GHz (red arrow) has a similar origin as the one marked also with a red arrow in Figure 3. Blue arrows mark torsional vibrational modes schematically depicted in the corresponding insets and are the result of a crosstalking of horizontal and vertical alignment of the bars. Above these two resonances, there is another mechanical mode related to the peak marked with a green arrow in Figure 3.</figcaption></figure>

nanobars, to find the allowed in-phase and counter-phase resonances for the interacting bars [23]–[25], we thermally excite only one of the bars (see Figure 1). We present in Figure 3 the resonances found in the high- and in the low-frequency regions (4–6.5 GHz -left- and 10–11.5 GHz -right-), for the mirror-symmetric structures. In the high-frequency region, the single bar presents only one resonance at 10.8 GHz corresponding to a longitudinal-like mode. However, the resonance at high frequency splits into two when the bars are aligned along the short axis: in-phase (lower frequency) and counter-phase (higher frequency) longitudinal modes. The alignment along the long axis leads to a single, in-phase, resonance at the same region.

Let us now focus on the low-frequency range, where the single bar shows no resonances. In this region, two distinct peaks can be observed for the other symmetric alignments, along the short and long axes, at ca. 4.6 GHz and ca. 5.7 GHz. These peaks result from the interactions between the two bars through the substrate. The mode is essentially the same for both alignments, manifesting as low-frequency modes oscillating along the coupling direction, maintaining the appropriate mirror symmetry.

Once the systems with mirror symmetries have been analyzed, we will consider the geometry depicted in Figure 1(d), where the shift introduced by $D_y$ lifts the mirror symmetry and the only remaining one is a C2 in-plane rotational symmetry. In Figure 4, we show the simulation results for the structure with $D_y = 75$ nm.

The high-frequency range presents a behavior similar to the short-axis alignment, mirror-symmetric geometry, but with a smaller splitting. The mechanical mode profile of these two resonances is equivalent to the short-axis alignment.

The low-frequency range presents a richer behavior. Similarly to the optical case, this geometry exhibits the two resonances present in the two mirror-symmetric ones, with a small frequency shift and are signaled with a red (green) arrow for the resonance associated to the short-axis (long-axis) alignment. The mechanical deformation profiles of these two resonances are equivalent to those schematically presented in Figure 3, maintaining to a very high degree the mirror symmetries, and there is essentially no torsion in each of the bars. In contrast to the optical case, when mirror symmetries are lifted not only cross-coupling between long-axis and short-axis modes is allowed, but extra mechanical resonances at low frequencies are clearly visible at ~5.7 GHz and ~5.8 GHz. These new resonances display a mechanical deformation profile different to the previous ones. In contrast to the previous ones, the deformation is characterized by torsion. The two bars experience a torsion since, in both cases, the displacement of the ends of the bars is in the direction of the short axis and in opposite directions. The two resonances, as depicted as insets in Figure 4, correspond to in-phase and counter-phase torsional deformations of each bar.
It is important to note that the “counterphase” res-
onances at ~11 GHz in Figure 3 and at ~5.8 GHz and in
~10.8 GHz Figure 4 are not allowed for a simultaneous ini-
tial thermal expansion of the two bars. The Q-factor that
these two latter resonances present, being the Q factor
defined as f/Δf, reveals that Q is larger for the fundamen-
tal modes appearing around 11 GHz. The exact values for
the Q-factor at the frequencies of 5.8 GHz and 10.8 GHz are
Q<sub>5.8GHz</sub> = 339 and Q<sub>10.8GHz</sub> = 602, respectively.

**4 Summary and conclusion**

We examined the effects of broken mirror symmetries in
the optical and nanomechanical response of acoustoplas-
monic resonators using a system composed of two identical
elongated metallic bars. Our findings show that the absence
of mirror symmetries affects differently the optical and
nanomechanical response.

The analysis of the optical resonances through the opti-
cal cross sections of mirror-symmetric and C2-symmetric
coupled bars demonstrates that broken mirror symme-
tries induces a crosstalking between the resonances of the
mirror-symmetric geometries. This is manifested by the
excitation of resonances corresponding to the long axis
when the light incidence is cross-polarized along the short
axis.

In the nanomechanical case, the situation is quite dif-
ferent. Interactions between bars produce a splitting of the
nanomechanical resonances when aligned along the short
axis. Additionally, new unconventional modes, preserving
the mirror symmetry imposed by the symmetry, appear in
the low-frequency region. Broken mirror symmetries not
only couple these nanomechanical modes but introduce
new resonant ones, whose deformation displays a distinct
torsional profile.

To conclude, we show that broken mirror symmetries
affect optical and nanomechanical resonances in a different
way. Even in simple systems, such as the one considered
here involving mirror-symmetric but interacting elements,
it is possible to excite highly nontrivial resonances mediated
by interactions. These new nanomechanical modes present
a marked torsional deformation profile. The elastic energy
will eventually be transferred to the substrate, showing a
spatiotemporal profile directly determined by the excited
bar mode.

**Research funding:** We acknowledge financial support from several founding agencies. BCLdL and AG-M acknowledge the Spanish Ministry of Science and Innovation (MCIN),

AEI and FEDER (UE) through projects PID2022-137569NB-C41 and TED2021-131417B-I00, and JGM acknowledges support through project PID2021-124193OB-C22. NDL-K and CX acknowledge support from the European Research Council Consolidator Grant No. 101045089 (T-Recs).

**Author contributions:** AG-M and JMG conceived the idea and the geometry to be employed. BCLdL carried the optical calculations that were discussed with AG-M. AG-M, JMG, and CX performed the acoustic calculations that were additionally discussed with BCLdL and NDL-K. All authors participated in the discussion and interpretation of the results. All authors cowrote the manuscript, have read, and agreed to the published version of the manuscript. All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

**Conflict of interest:** Authors state no conflict of interest.

**Data availability:** All data needed to evaluate the conclusions in the paper are present in the paper.

**References**

[1] K. Brading, E. Castellani, and N. Teh, “Symmetry and symmetry breaking,” in *The Stanford Encyclopedia of Philosophy*, Fall 2023 ed E. N. Zalta and U. Nodelman, Eds., Stanford, Metaphysics Research Lab, Stanford University, 2023.

[2] G. Zheng, et al., “Discrete metal nanoparticles with plasmonic chirality,” *Chem. Soc. Rev.*, vol. 50, no. 6, pp. 3738–3754, 2021.

[3] B. Ni, et al., “Chiral seeded growth of gold nanorods into fourfold twisted nanoparticles with plasmonic optical activity,” *Adv. Mater.*, vol. 35, no. 1, p. 2208299, 2023.

[4] B. Castillo López de Larrinzar, et al., “Photonic and nanomechanical modes in acoustoplasmonic toroidal nanopropellers,” *Nanomaterials*, vol. 14, no. 15, 2024, https://doi.org/10.3390/nano14151276.

[5] Z. Han, et al., “Recent advances in ultrathin chiral metasurfaces by twisted stacking,” *Adv. Mater.*, vol. 35, no. 3, 2022. https://doi.org/10.1002/adma.202206141.

[6] B. Castillo López de Larrinzar, et al., “Towards chiral acousto-plasmonics,” *Nanophotonics*, vol. 12, no. 11, pp. 1957–1964, 2023.

[7] Priya, E. R. Cardozo de Oliveira, and N. D. Lanzillotti-Kimura, “Perspectives on high-frequency nanomechanics, nanoacoustics, and nanophononics,” *Appl. Phys. Lett.*, vol. 122, no. 14, p. 140501, 2023.

[8] N. D. Lanzillotti-Kimura, et al., “Polarization-controlled coherent phonon generation in acoustoplasmonic metasurfaces,” *Phys. Rev. B*, vol. 97, no. 23, p. 235403, 2018.

[9] F. Della Picca, et al., “Tailored hypersound generation in single plasmonic nanoantennas,” *Nano Lett.*, vol. 16, no. 2, pp. 1428–1434, 2016.

[10] K. O’Brien, et al., “Ultrafast acousto-plasmonic control and sensing in complex nanostructures,” *Nat. Commun.*, vol. 5, p. 4042, 2014.

[11] O. Ortiz, et al., “Fiber-integrated microcavities for efficient generation of coherent acoustic phonons,” *Appl. Phys. Lett.*, vol. 117, no. 18, p. 183102, 2020.
[12] L. Ding, et al., “High frequency gaas nano-optomechanical disk resonator,” *Phys. Rev. Lett.*, vol. 105, no. 26, p. 263903, 2010.

[13] D. K. Armani, et al., “Ultra-high-Q toroid microcavity on a chip,” *Nature*, vol. 421, pp. 925–928, 2003.

[14] A. Bragas, et al., “Nanomechanics with plasmonic nanoantennas: ultrafast and local exchange between electromagnetic and mechanical energy,” *J. Opt. Soc. Am. B*, vol. 40, no. 5, 2023. https://doi.org/10.1364/josab.482384.

[15] O. Ortiz, et al., “Topological optical and phononic interface mode by simultaneous band inversion,” *Optica*, vol. 8, no. 5, pp. 598–605, 2021.

[16] P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” *Phys. Rev. B*, vol. 6, no. 12, p. 4370, 1972.

[17] M. Cotrufo, C. I. Osorio, and A. F. Koenderink, “Spin-dependent emission from arrays of planar chiral nanoantennas due to lattice and localized plasmon resonances,” *ACS Nano*, vol. 10, no. 3, pp. 3389–3397, 2016.

[18] A. G. Curto, et al., “Unidirectional emission of a quantum dot coupled to a nanoantenna,” *Science*, vol. 329, no. 5994, pp. 930–933, 2010.

[19] E. Prodan, et al., “A hybridization model for the plasmon response of complex nanostructures,” *Science*, vol. 302, no. 5644, pp. 419–422, 2003.

[20] B. Wang, et al., “Chiral metamaterials: simulations and experiments,” *J. Opt. A: Pure Appl. Opt.*, vol. 11, no. 11, p. 114003, 2009.

[21] B. C. L. de Larrinzar, N. D. Lanzillotti-Kimura, and A. García-Martín, “Interaction effects in chiral acoustoplasmonic nanostructures,” in *Nanophotonics X*, vol. 12991, 129910B, International Society for Optics and Photonics, SPIE, Bellingham, Washington, USA, 2024.

[22] E. Cortés, et al., “Plasmonic hot electron transport drives nano-localized chemistry,” *Nat. Commun.*, vol. 8, p. 14880, 2017.

[23] J. Lermé, J. Margueritat, and A. Crut, “Vibrations of dimers of mechanically coupled nanostructures: analytical and numerical modeling,” *J. Phys. Chem. C*, vol. 125, no. 15, pp. 8339–8348, 2021.

[24] A. Girard, et al., “Acoustic mode hybridization in a single dimer of gold nanoparticles,” *Nano Lett.*, vol. 18, no. 6, pp. 3800–3806, 2018.

[25] A. Crut, “Substrate-supported nano-objects with high vibrational quality factors,” *J. Appl. Phys.*, vol. 131, no. 24, p. 244301, 2022.

**Supplementary Material:** This article contains supplementary material (https://doi.org/10.1515/nanoph-2024-0482).