# Rayleigh-Ritz Approximation of the Acoustic Vibrations of Clamped Superquadrics-Application to Free Core-Shell Objects

Sajana S, María del Carmen Marco de Lucas, Lucien Saviot

► To cite this version:

Sajana S, María del Carmen Marco de Lucas, Lucien Saviot. Rayleigh-Ritz Approximation of the Acoustic Vibrations of Clamped Superquadrics-Application to Free Core-Shell Objects. Nanomaterials, 2025, 15 (24), pp.1865. 10.3390/nano15241865. hal-05412736

HAL Id: hal-05412736

https://ube.hal.science/hal-05412736v1

Submitted on 12 Dec 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
Article

# Rayleigh-Ritz Approximation of the Acoustic Vibrations of Clamped Superquadrics—Application to Free Core-Shell Objects

Sajana S ID, María del Carmen Marco de Lucas ID and Lucien Saviot * ID

Université Bourgogne Europe, CNRS, Laboratoire Interdisciplinaire Carnot de Bourgogne ICB UMR 6303, 21000 Dijon, France

* Correspondence: lucien.saviot@cnrs.fr

## Abstract

A numerical approach based on the Rayleigh-Ritz method and using a modification of the so-called *xyz* algorithm is introduced to calculate the acoustic vibrations of clamped objects whose shape is delimited by superquadrics. It is then used to improve the convergence for the free vibrations of core-shell objects. The issue in this case is first illustrated in the simpler one-dimensional case of the thickness breathing vibration of an infinite "core-shell" plate. Functions suitable for solving the clamped vibrations of the core are added to the original *xyz* basis of functions to improve the convergence for core-shell superquadrics. The new basis obeys the same symmetry rules as the original one, which allows calculating vibrations for individual irreducible representations when the objects are made of cubic, tetragonal, or orthorhombic materials whose principal axes are aligned with those of the superquadrics. This method is validated for an isotropic spherical core-shell system for which analytic solutions exist.

**Keywords:** Rayleigh-Ritz method; acoustic vibration; core-shell; superquadrics; rigid boundary condition

## 1. Introduction

The acoustic eigenfrequencies of an object depend on its size, shape, the material it is made of, and its environment. Measuring and modeling these frequencies provides valuable information about all these parameters. Resonant ultrasound spectroscopy (RUS) implements this approach and is used mainly in geology and material science to determine the elastic constants of a material [1]. Fast numerical methods are required in this context to calculate the eigenfrequencies quickly enough to allow fitting the measured frequencies. The Rayleigh-Ritz (RR) method fulfills this requirement [2]. This method can be applied to objects with very different dimensions, down to nanomaterials. Different experimental methods allow measuring the eigenfrequencies of nano-objects in the time or frequency domains (for example, transient absorption and low-frequency Raman scattering) [3,4]. Recent experimental results in this domain have concerned more complex nano-objects in terms of shape and composition. In particular, gold-silver core-shell nanocrystals with spherical shapes have been synthesized and investigated with inelastic light scattering in order to assert the epitaxy (or lack thereof) between the gold core and the silver shell [5]. Gold-silver core-shell nanorods have also been synthesized [6] and their vibrations have been investigated with inelastic light scattering [7]. Their geometry is more complex

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

### check for updates

Academic Editor: Francisco Torrens

Received: 12 November 2025

Revised: 28 November 2025

Accepted: 8 December 2025

Published: 11 December 2025

**Citation:** S, S.; Marco de Lucas, M.d.C.; Saviot, L. Rayleigh-Ritz Approximation of the Acoustic Vibrations of Clamped Superquadrics—Application to Free Core-Shell Objects. *Nanomaterials* **2025**, *15*, 1865. https://doi.org/10.3390/nano15241865

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
because the silver coating is not homogeneous. The gold core has a circular cross-section, and the silver shell cross-section tends to a square with increasing thicknesses, and all these parameters must be taken into account to model their optical properties [6]. Superquadrics have been successfully used in this context to model the shapes of the core and the shell [8]. For these two core-shell systems, the gold core and the silver shell have a cubic lattice structure; therefore, anisotropic elasticity must also be taken into account.

Accurate modeling of such nano-objects is required to assign the different eigenfrequencies observed experimentally. The continuum elasticity approximation is valid in this context as long as the objects' sizes are large enough, i.e., the wavelength of the elastic waves is large compared to the interatomic distance. The finite element method (FEM) is a very versatile method to compute eigenfrequencies for arbitrary shapes and elastic properties. However, in the previous examples, the 3D objects are invariant under many symmetry operations, which can be leveraged to help assign the experimental features as usually performed in spectroscopy. For example, transient absorption is essentially sensitive to eigenmodes for which a volume variation occurs during oscillation, i.e., the totally symmetric ones. Raman scattering also obeys selection rules preventing many eigenmodes from contributing to the spectra. Further processing of FEM calculations is generally required to recover the symmetry of the eigenmodes. RR can take advantage of symmetries to calculate only the eigenmodes of interest, i.e., those corresponding to specific irreducible representations. This also speeds up the calculations significantly and lowers the required computing resources.

The present work first extends to superquadrics a recently introduced method to calculate the vibrations of clamped spheres (rigid boundary condition, i.e., vanishing displacement at the surface) [9] by providing suitable basis functions meeting these requirements. This boundary condition is seldom encountered in practice, but it is a very useful limiting case, which helps figure out how the vibrations of objects embedded in a harder and/or denser matrix may differ compared to the free vibrations of the same object [10]. The vibrations of core-shell objects are then considered. They can be calculated using the usual RR method without modifications [5]. However, the convergence of the calculations is not always as good as for homogeneous objects, questioning the validity of this simple approach. This comes from the fact that the RR method is implemented with functions of class $C^\infty$, while the derivative of the displacement is generally different on both sides of an interface in core-shell systems. The present work proposes a method to improve the convergence using additional basis functions similar to those introduced before. This method is first demonstrated in the simpler 1D case of the thickness breathing vibrations of infinite plates. A calculation for a core-shell sphere is then presented and compared to the analytical solution to show the improved convergence.

## 2. Methods

RR calculations are based on the xyz algorithm introduced by Visscher et al. [2]. The main points are summarized below, and readers are referred to this paper for a comprehensive presentation. The three Cartesian components of the displacements are decomposed on the $\phi_\lambda = x^p y^q z^r$ basis of functions with p, q, and r integers, $\lambda = (p, q, r)$, and $p + q + r \le N$. The RR approach turns the dynamic problem into the generalized eigenvalue problem $\omega^2 Ea = \Gamma a$ with $\omega$ the pulsation, $a$ the amplitudes of the $\phi_\lambda$ components, and $E$ and $\Gamma$ two square matrices defined as follows:

$$E_{\lambda i \lambda' i'} = \delta_{ii'} \int_V \phi_\lambda \rho \phi_{\lambda'} dV \quad (1)$$

$$\Gamma_{\lambda i \lambda' i'} = C_{ij i' j'} \int_V \phi_{\lambda, j} \phi_{\lambda', j'} dV \quad (2)$$
where C$_{ijkl}$ is the stiffness tensor, ρ the mass density, and V the volume of the object. The Γ matrix is block-diagonal when the shape and elastic properties are symmetric with respect to the x = 0, y = 0, or z = 0 planes [11]. In those cases, the linear system can be turned into a few smaller linear systems requiring less computer time and resources. The xyz algorithm is well-suited for shapes for which analytic expressions for the volume integrals exist (cubes, spheres, cylinders, ...). It was extended in a previous work [12] to superquadrics, which provides a large variety of shapes from octahedra to spheres and cubes, which are good approximations for the rounded shapes discussed above [6], but also to wires with cross-sections delimited by superellipses.

For the calculations, we considered materials with a strong acoustic impedance mismatch to highlight problematic cases. We chose an inorganic core (gold) and an organic shell (polystyrene). Similar nano-objects are commonly investigated experimentally and theoretically, as in ligand-covered nanoparticles [13] or nanoplatelets [14] and gold nanoplates coated with an organic layer [15]. For gold, we used isotropic parameters with mass density ρ = 19.283 g/cm³ and sound speeds v$_{L}$ = 3330 m/s and v$_{T}$ = 1250 m/s [12]. For polystyrene, we used ρ = 1.05 g/cm³, v$_{L}$ = 2350 m/s, and v$_{T}$ = 1210 m/s [16]. The acoustic impedance for a plane wave at normal incidence is ρv$_{L}$; therefore, the mismatch comes mainly from the very different mass densities of both materials in that case. This large difference is expected to provide a strong impedance mismatch in non-planar geometries too.

## 3. Results

### 3.1. Superquadrics with Rigid Boundary Conditions

We first modify the RR approach presented before to calculate the vibrations of clamped superquadrics whose surface is defined by the implicit Equation (3), where L$_{i}$ is the half-length in the i direction and n$_{i}$ > 0 are the shape parameters (real numbers).

$$ \left| \frac{x}{L_x} \right|^{n_x} + \left| \frac{y}{L_y} \right|^{n_y} + \left| \frac{z}{L_z} \right|^{n_z} = 1 \qquad (3) $$

The vibrations of such objects with the rigid boundary condition can be obtained using a basis of functions that vanish at the surface of the superquadrics. Such a basis was previously introduced for spheres in Nestoklon et al. [9] and is extended to superquadrics in Equation (4). For parallelepipeds (n$_{x,y,z}$ → ∞), Equation (5) can be used instead. The functions are of class C² for n ≥ 2 and positive inside the superquadrics. They also vanish at the surface. Their integrals in Equations (1) and (2) can be computed in the first octant (x, y, z > 0) [12] because they are just linear combinations of the ones in the original xyz method. Values in the other octants are obtained easily by considering the parity of p, q, and r. The symmetry considerations used for the free vibrations also remain valid allowing us to compute independently the frequencies for each irreducible representation.

$$ u_{pqr}(x,y,z) = x^p y^q z^r \times \left( 1 - \left| \frac{x}{L_x} \right|^{n_x} - \left| \frac{y}{L_y} \right|^{n_y} - \left| \frac{z}{L_z} \right|^{n_z} \right) \qquad (4) $$

$$ u_{pqr}(x,y,z) = x^p y^q z^r \times \left(1 - \frac{x^2}{L_x^2}\right) \left(1 - \frac{y^2}{L_y^2}\right) \left(1 - \frac{z^2}{L_z^2}\right) \qquad (5) $$

As an example, Figure 1 shows the frequencies computed for superquadrics made of gold with octahedral symmetry (L$_{x}$ = L$_{y}$ = L$_{z}$ and n$_{x}$ = n$_{y}$ = n$_{z}$ = n) and p + q + r ≤ N = 20. The shape varies from a sphere when the shape parameter n = 2 to a cube as n → ∞. We took into account the symmetry of the problem (point group O$_{h}$) to calculate the vibrations corresponding to each irreducible representation. Only the Raman active
ones are shown: A₁g, Eg, and T₂g. As with the free vibrations [12], these frequencies show
only small variations when the shape varies from a sphere to a cube after normalization by
multiplying with the cubic root of the volume.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Shape variation of the A₁g (**left**), Eg (**center**) and T₂g (**right**) frequencies of nanoparticles made of isotropic gold with the rigid boundary condition. The frequencies are normalized after multiplication by the cubic root of the volume and are expressed in m/s. The superspheroid shape varies from a sphere ($n = 2$) to a cube ($n \to \infty$). Normalized frequencies calculated from the analytical expressions for the sphere are plotted with round symbols. Red disks are for the spheroidal modes with $l = 0$ (A₁g) of $l = 2$ (Eg + T₂g). Other spheroidal or torsional modes are plotted with blue circles. FEM eigenfrequencies for the cube are plotted with empty blue squares.</figcaption></figure>

The convergence of the algorithm was checked for an isotropic sphere for which
analytical solutions exist (markers in Figure 1) [17]. The difference between the frequencies
obtained with both methods is less than 1% up to frequency × V¹/³ ≃ 6900 m/s. When
considering all the possible irreducible representations, more than the first 750 lowest fre-
quencies are reproduced to better than 0.1%. The frequencies for the cube using Equation (5)
were also checked with FEM calculations with the FreeFem++ software (version 4.15) [18]
using a mesh characteristic length 30 times smaller than the edge of the cube. A very
good agreement was obtained with differences of <0.1% up to ≃4740 m/s (300 lowest
frequencies considering all the irreducible representations). Overall the convergence of
the RR calculations for the rigid boundary condition is as good as with the free one for
spheres and cubes. We checked the convergence of the calculations using Equation (4)
as 1/n → 0 toward the values calculated for the cube with Equation (5). We performed a
linear extrapolation of the value obtained for 1/n = 1/50 and 1/n = 1/100 and compared it with
the results discussed above. The difference was <1% below 5000 m/s demonstrating a very
good convergence of this method too, even for rather large values of n.

The Raman active vibrations for the sphere are the spheroidal modes with $l = 0$ (A₁g in O_h) and $l = 2$ (Eg + T₂g in O_h), which are marked with disks in Figure 1. For the free vibrations, the breathing modes can also be identified by calculating the variation of the volume during oscillation. However, this method cannot be applied when the surface is clamped because the volume does not change. We also note that the method introduced here with Equation (4) is applicable when the shape parameter $n > 1$. We only considered $n \ge 2$ in this work because the second derivative of some basis functions diverges at the origin for $1 < n < 2$.

3.2. *Thickness Breathing Vibration of Core–Shell Slabs*

To demonstrate the issues faced with the xyz method for core–shell objects, we start
with a simpler 1D system for which analytic solutions exist and displacements can be
represented more easily. We consider the thickness breathing vibration of an infinite “core-shell” plate with a gold core of thickness $h$ covered on both faces by layers of polystyrene with thicknesses $a$ (total thickness $h + 2a$ between $z = -h/2 - a$ and $z = h/2 + a$). We only consider vibrations at the center of the Brillouin zone $\vec{u}(x, y, z) = u(z)\vec{e}_z$ (vanishing wavevector) with odd longitudinal displacements $u(-z) = -u(z)$. For the analytical solution, $u(z) = A \sin(q_{Au}z)$ inside the gold core, $u(z) = B \cos(q_{ps}z) + C \sin(q_{ps}z)$ for the upper polystyrene layer ($z > h/2$) and $u(z) = -B \cos(q_{ps}z) + C \sin(q_{ps}z)$ for the bottom layer ($z < -h/2$) with $\omega = q_{Au}v_{L,Au} = q_{ps}v_{L,ps}$. Eigenfrequencies calculated by searching the roots of the secular equation obtained from the usual continuity equations at the interfaces and the free boundary condition at the surfaces are given in Table 1 (middle). The displacements of the lowest frequency eigenmodes normalized at the surface are shown in Figure 2 (colored areas).

<table><caption>**Table 1.** Eigenfrequencies calculated using the analytic method and the xyz and xyz + r algorithms for the first overtones of the thickness breathing eigenmodes of an infinite core-shell plate with a gold core (thickness $h$) sandwiched between two polystyrene shell slabs (thicknesses $a = h/2$). The frequencies are multiplied by the total thickness $h + 2a$ and expressed in m/s. $N_b$: the number of functions in the basis, $N_{max}$: largest z exponent used in the basis. The double arrows in the first column indicate the direction for improved accuracy.</caption><thead><tr><th>Method</th><th>N<sub>b</sub></th><th>N<sub>max</sub></th><th colspan="6">Frequency × Total Thickness (m/s)</th></tr></thead><tbody><tr><td rowspan="6">xyz</td><td>6</td><td>11</td><td>2608.98</td><td>3490.04</td><td>9059.21</td><td>9973.30</td><td>16,066.04</td><td>34,164.95</td></tr><tr><td>7</td><td>13</td><td>2539.65</td><td>3472.34</td><td>8263.32</td><td>9969.51</td><td>15,393.52</td><td>17,773.62</td></tr><tr><td>8</td><td>15</td><td>2533.52</td><td>3470.00</td><td>8259.65</td><td>9966.11</td><td>14,698.40</td><td>16,991.05</td></tr><tr><td>↓</td><td>9</td><td>17</td><td>2471.79</td><td>3456.12</td><td>7927.41</td><td>9962.83</td><td>13,944.58</td><td>16,599.73</td></tr><tr><td>10</td><td>19</td><td>2453.82</td><td>3453.18</td><td>7821.38</td><td>9962.68</td><td>13,356.83</td><td>16,568.04</td></tr><tr><td>11</td><td>21</td><td>2445.02</td><td>3451.14</td><td>7787.45</td><td>9961.65</td><td>13,323.51</td><td>16,560.45</td></tr><tr><td>↓</td><td>12</td><td>23</td><td>2412.39</td><td>3445.44</td><td>7645.69</td><td>9961.65</td><td>12,953.16</td><td>16,556.57</td></tr><tr><td>analytic</td><td></td><td></td><td>2247.53</td><td>3423.31</td><td>7039.59</td><td>9958.12</td><td>11,800.09</td><td>16,191.16</td></tr><tr><td rowspan="6">xyz + r</td><td>22</td><td>23</td><td>2247.53</td><td>3423.31</td><td>7039.59</td><td>9958.12</td><td>11,800.09</td><td>16,191.26</td></tr><tr><td>20</td><td>21</td><td>2247.53</td><td>3423.31</td><td>7039.59</td><td>9958.12</td><td>11,800.10</td><td>16,191.58</td></tr><tr><td>18</td><td>19</td><td>2247.53</td><td>3423.31</td><td>7039.60</td><td>9958.12</td><td>11,800.10</td><td>16,196.30</td></tr><tr><td>16</td><td>17</td><td>2247.53</td><td>3423.31</td><td>7039.61</td><td>9958.14</td><td>11,802.25</td><td>16,273.71</td></tr><tr><td>↑</td><td>14</td><td>15</td><td>2247.54</td><td>3423.31</td><td>7039.82</td><td>9958.20</td><td>11,838.04</td><td>16,316.39</td></tr><tr><td>12</td><td>13</td><td>2247.64</td><td>3423.31</td><td>7041.55</td><td>9958.96</td><td>11,890.71</td><td>16,557.73</td></tr><tr><td>↑</td><td>10</td><td>11</td><td>2248.41</td><td>3423.31</td><td>7042.05</td><td>9979.64</td><td>12,185.18</td><td>16,644.00</td></tr><tr><td>8</td><td>9</td><td>2254.85</td><td>3423.44</td><td>7188.40</td><td>9991.92</td><td>16,590.86</td><td>20,576.40</td></tr><tr><td>6</td><td>7</td><td>2313.28</td><td>3430.40</td><td>9743.33</td><td>11,191.98</td><td>17,406.92</td><td>43,458.30</td></tr></tbody></table>

We implemented the xyz algorithm to solve the same problem by also considering only odd displacements with functions $u_i(z) = z^{2i+1}$ and replacing the volume integrals in Equations (1) and (2) by integrals from $z = -h/2 - a$ to $z = h/2 + a$. The results are also shown in Table 1 and Figure 2 (black curves). The convergence is slow, as shown by the variations of the frequencies as the number of functions in the basis $N_b$ increases. The frequency of the fundamental mode deviates by more than 150 m/s (~7%) from the analytic solution, which is large when compared with the accuracy generally for the free vibrations. The deviations for the highest overtones are very large. This slow convergence is confirmed by the corresponding displacements, which all deviate noticeably from the analytical solution. The disagreement is stronger when the change of slope at the interface is more abrupt.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Displacements for the six eigenmodes of the core–shell slab given in Table 1. The normalized frequencies (in m/s) are indicated in the top left part of each plot. The displacements are anti-symmetrical ($u_z(-z) = -u_z(z)$) and represented for $z \ge 0$ only. The analytical solutions are shown with filled curves (yellow for gold and grey for polystyrene). The xyz displacements are plotted with a black line and the xyz + r ones in red.</figcaption></figure>

To overcome this issue, we propose to supplement the original basis with $u'_i$ functions, which would be suitable for the core with the rigid boundary condition and which are extended in the shell with a vanishing value as presented in Equation (6). A few $u_i$ and $u'_i$ functions are plotted in Figure 3 after normalization at their maximum. The $u'_i$ are continuous at the core–shell interface because they vanish at $z = \pm h/2$ and are differentiable inside each domain. They are linear combinations of the original $u_i$ functions because $u'_i(z) = u_i(z) - u_{i+1}(z)/(h/2)^2$, but for $|z| \le h/2$ only. Therefore, the integrals in Equations (1) and (2) still have analytical expressions. Finally, the parity of the $u_i$ and $u'_i$ functions are the same, and the previous symmetry considerations remain unchanged. The key point with this choice is that the left and right derivatives at the interface are different even for small values of $i$, which is expected to help fit the analytical displacement.

$$u'_{i}(z) = \begin{cases} z^{2i+1} \times \left(1 - \left(\frac{z}{h/2}\right)^2\right), & |z| \le \frac{h}{2} \\ 0, & |z| \ge \frac{h}{2} \end{cases} \qquad (6)$$
<figure><img src="image_5.png" /><figcaption>**Figure 3.** $u_i$ and $u_i'$ functions ($i = 0, \dots, 3$) used in the $xyz + r$ algorithm normalized at their maximum. The background color is yellow for gold and grey for polystyrene.</figcaption></figure>

The frequencies and displacements obtained with the new basis (called $xyz + r$ from now on) are also shown in Table 1 (bottom) and Figure 2 (red curves). The convergence is significantly improved when compared with the $xyz$ basis for an identical maximum exponent value ($N_{max} = 23$) or an identical number of functions in the basis ($N_b = 12$). For the fundamental mode, the frequency deviates from the analytical one by less than $3 \times 10^{-5}$ m/s compared to more than 160 m/s with $xyz$. In particular, we note that the first two frequencies obtained for the $xyz + r$ method with the smallest number of functions considered in this work ($N_b = 6$ and $N_{max} = 7$) are closer to the analytical values than the frequencies calculated with $xyz$ up to $N_b = 12$ and $N_{max} = 23$. The $xyz + r$ displacements are also very close to the analytical ones (indistinguishable in Figure 2). The maximum difference in the normalized displacements is $<10^{-5}$ for the two lowest frequency modes, $<4 \times 10^{-5}$ for the third, $<2 \times 10^{-4}$ for the fourth and fifth and $<3 \times 10^{-3}$ for the sixth. It is $>5 \times 10^{-2}$ for all the modes calculated with $xyz$. Adding the rigid boundary functions to improve the convergence of the RR method is therefore very promising and will be tested in the following with 3D objects.

We note that another approach has been proposed to improve the convergence for two other 1D problems with stepped structures (longitudinal vibrations of bars and flexural vibrations of rectangular plates) by Elishakoff et al. [19]. The proposed method introduces additional rigorous contributions to the linear system involving values of the functions at the interface. Its extension to 3D objects is challenging because it would require the calculation of surface integrals. These are unlikely to have analytic expressions for superquadrics for which no such expression exists for the surface area. This is the reason why this approach has not been considered in this work.

### 3.3. 3D Core–Shell

We now apply the method presented above to improve the convergence of the $xyz$ algorithm for 3D core–shell objects. To do so, the basis of functions is supplemented with the $u_{pqr}$ functions presented in Equation (4) inside the core, which are extended with a vanishing value inside the shell. This method was tested for a gold–polystyrene core–shell sphere using the same material parameters as before. The outer diameter is 1.2 times the inner diameter. The results are compared to analytic calculations and presented in Figure 4 for all the irreducible representations. Modes up to the second $S_8$ mode at frequency $\times$ diameter $\simeq 5802.16$ m/s) are shown corresponding to the 593th mode when all the degeneracies are considered.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Relative error $((v_{\text{Rayleigh-Ritz}} - v_{\text{analytic}}) / v_{\text{analytic}})$ for the xyz (black disks) and $xyz + r$ (red triangles) frequencies for a core-shell sphere with an isotropic gold core and a polystyrene shell as a function of the normalized frequency (frequency × diameter). The ratio of the shell to core diameters is 1.2. Errors below $10^{-10}$ are shown as $10^{-10}$. The positions of the Raman active spheroidal modes with $l = 0$ and 2 are indicated with arrows at the top. The position of the lowest frequency mode (torsional) for each angular momentum $l$ is also shown at the bottom.</figcaption></figure>

The differences between computed and analytic values are always smaller with $xyz + r$, showing that the convergence is indeed improved by the additional functions. The frequency difference is less than 1% in most cases. It is noticeably larger for modes with large $l$ and in particular for $l \ge 9$. Indeed, the number of angular oscillations increases with $l$, requiring larger values of $N$ to reproduce the corresponding displacement. The third $S_2$ mode at 3485.73 m/s is not affected by this issue. Still, its frequency is poorly reproduced with $xyz$ (10% error). This is problematic for a Raman-active vibration with a frequency below that of the first breathing mode ($S_0$) where Raman peaks are most frequently observed. This error is reduced by more than 9 orders of magnitude with $xyz + r$. Similar comments apply for many other modes. Overall, $xyz$ underperforms in most cases for this particular core-shell system. The errors are significantly larger than the ones reported above for the sphere with the rigid boundary condition. The addition of rigid functions increases significantly the convergence, resulting in errors less than $10^{-2}$ in most cases and $\le 10^{-4}$ up to the first breathing mode ($S_0$).

We used the largest possible value of $N$ for the $xyz$ and $xyz + r$ calculations in Figure 4. Increasing $N$ further results in errors in the generalized eigenvalue problem solver. Overcoming this issue would require moving to a floating-point format larger than double-precision, which goes beyond the scope of this work. We used $N = 20$ for $xyz$ calculations and $N = 16$ for $xyz + r$, resulting in a slightly larger number of functions in the $xyz + r$ basis. To confirm that the improved convergence was not just a consequence of the larger basis, we also compared with $xyz + r$ calculations performed with $N = 14$. In that case, the bases were smaller for $xyz + r$ for all the irreducible representations. The convergence was still better except for modes with large values of $l$ such as $S_{10}$ and $T_{11}$ because a lower value of $N$ makes it harder to reproduce the angular oscillations in these cases.
## 4. Conclusions

To summarize, we first presented in this work an extension of the RR method to calculate the vibrations of clamped objects applicable to a large variety of shapes, from 3D superquadrics to cylinders with superellipse cross-sections. The convergence of this method is comparable to that of the original method for the free vibrations. The new functions introduced for the clamped vibrations were then shown to help improve the convergence for core–shell systems. This was first illustrated with the 1D case of the thickness breathing vibrations of “core-shell” plates. Its application to a core–shell sphere demonstrated a significantly improved convergence. The calculations presented here are applicable to a large variety of actual systems. Three-dimensional core–shell superquadrics with different parameters ($L_i, n_i$) for the core and the shell are suitable to model existing experimental inelastic light scattering results obtained for gold–silver nanorods [7]. It is also applicable to core–shell nanowires with different superellipse cross-sections for the core and the shell [12], which can model nanorods with lengths much larger than their diameters. The advantages of the RR method in terms of low computer resource requirements and symmetry are preserved, making this approach very useful to fit shape or elasticity parameters to match experimental results.

**Author Contributions:** Conceptualization, L.S.; methodology, L.S.; software, L.S.; validation, L.S., S.S., and M.d.C.M.d.L.; formal analysis, L.S.; investigation, L.S. and S.S.; resources, L.S. and M.d.C.M.d.L.; data curation, L.S.; writing—original draft preparation, L.S.; writing—review and editing, L.S., S.S., and M.d.C.M.d.L.; visualization, L.S.; supervision, L.S. and M.d.C.M.d.L.; project administration, L.S. and M.d.C.M.d.L.; funding acquisition, L.S. and M.d.C.M.d.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** We acknowledge support by ANR DeXCTer (ANR-23-CE29-0003) and the EIPHI Graduate School (ANR-17-EURE-0002) operated by the French National Research Agency (ANR).

**Data Availability Statement:** The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

**Conflicts of Interest:** The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

## References

1. Maynard, J. Resonant Ultrasound Spectroscopy. *Phys. Today* **1996**, *49*, 26–31. [CrossRef]

2. Visscher, W.M.; Migliori, A.; Bell, T.M.; Reinert, R.A. On the normal modes of free vibration of inhomogeneous and anisotropic elastic objects. *J. Acoust. Soc. Am.* **1991**, *90*, 2154–2162. [CrossRef]

3. Jansen, M.; Tisdale, W.A.; Wood, V. Nanocrystal phononics. *Nat. Mater.* **2023**, *22*, 161–169. [CrossRef] [PubMed]

4. Wright, C.; Hartland, G.V. The Science of Nanostructure Acoustic Vibrations. *Annu. Rev. Phys. Chem.* **2025**, *76*, 57–81. [CrossRef] [PubMed]

5. Lee, S.; Portalès, H.; Walls, M.; Beaunier, P.; Goubet, N.; Tremblay, N.; Margueritat, J.; Saviot, L.; Courty, A. Versatile and robust synthesis process for the fine control of the chemical composition and core-crystallinity of spherical core-shell Au@Ag nanoparticles. *Nanotechnology* **2021**, *32*, 95604. [CrossRef] [PubMed]

6. Tebbe, M.; Kuttner, C.; Mayer, M.; Maennel, M.; Pazos-Perez, N.; König, T.A.; Fery, A. Silver-Overgrowth-Induced Changes in Intrinsic Optical Properties of Gold Nanorods: From Noninvasive Monitoring of Growth Kinetics to Tailoring Internal Mirror Charges. *J. Phys. Chem. C* **2015**, *119*, 9513–9523. [CrossRef] [PubMed]

7. Vernier, C. Étude de Nanocristaux Métalliques Anisotropes par Spectroscopies Optique et Vibrationnelle. Ph.D. Thesis, Sorbonne Université, Paris, France, 2024. Available online: https://theses.fr/2024SORUS300 (accessed on 10 October 2025).

8. Vernier, C.; Portalès, H. Impact of tip curvature and edge rounding on the plasmonic properties of gold nanorods and their silver-coated counterparts. *J. Chem. Phys.* **2024**, *161*, 124711. [CrossRef] [PubMed]

9. Nestoklon, M.; Saviot, L.; Goupalov, S.V. Application of perturbation theory to finding vibrational frequencies of a spheroid. *Phys. Rev. B* **2025**, *in press*. [CrossRef]
10. Murray, D.B.; Saviot, L. Phonons in an inhomogeneous continuum: Vibrations of an embedded nanoparticle. *Phys. Rev. B* **2004**, *69*, 94305. [CrossRef]

11. Mochizuki, E. Application of group theory to free oscillations of an anisotropic rectangular parallelepiped. *J. Phys. Earth* **1987**, *35*, 159–170. [CrossRef]

12. Saviot, L. Free vibrations of anisotropic nano-objects with rounded or sharp corners. *Nanomaterials* **2021**, *11*, 1838. [CrossRef] [PubMed]

13. Lee, E.M.Y.; Mork, A.J.; Willard, A.P.; Tisdale, W.A. Including surface ligand effects in continuum elastic models of nanocrystal vibrations. *J. Chem. Phys.* **2017**, *147*, 44711. [CrossRef] [PubMed]

14. Kurus, N.N.; Milekhin, A.G.; Sklyar, R.I.; Saidzhonov, B.M.; Vasiliev, R.B.; Adichtchev, S.V.; Surovtsev, N.V.; Latyshev, A.V.; Zahn, D.R.T. Phonons in Core–Shell CdSe/CdS Nanoplatelets Studied by Vibrational Spectroscopies. *J. Phys. Chem. C* **2022**, *126*, 7107–7116. [CrossRef]

15. Devkota, T.; Yu, K.; Hartland, G.V. Mass loading effects in the acoustic vibrations of gold nanoplates. *Nanoscale* **2019**, *11*, 16208–16213. [CrossRef] [PubMed]

16. Cheng, W.; Wang, J.J.; Jonas, U.; Steffen, W.; Fytas, G.; Penciu, R.S.; Economou, E.N. The spectrum of vibration modes in soft opals. *J. Chem. Phys.* **2005**, *123*, 121104. [CrossRef] [PubMed]

17. Saviot, L.; Champagnon, B.; Duval, E.; Kudriavtsev, I.A.; Ekimov, A.I. Size dependence of acoustic and optical vibrational modes of CdSe nanocrystals in glasses. *J. Non-Cryst. Solids* **1996**, *197*, 238–246. [CrossRef]

18. Hecht, F. New development in FreeFem++. *J. Numer. Math.* **2012**, *20*, 251–265. [CrossRef]

19. Elishakoff, I.; Amato, M.; Ankitha, A.P.; Marzani, A. Rigorous implementation of the Galerkin method for stepped structures needs generalized functions. *J. Sound Vib.* **2021**, *490*, 115708. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.