Exploiting Modern Chladni Plates to Analogously Manifest the Point Interaction

Abstract

A truncated basis is employed to analyze the influence of the point interaction on the eigenvalues and eigenfunctions in quantum billiards. The point interaction is numerically confirmed to cause the shift in the eigenvalue that leads the original eigenfunctions to be superposed to form the new eigenfunction. The amplitude and nodal-line patterns of the eigenfunctions are found to vary significantly with the coupling strength. It is further verified that the point-driven Chladni plates can be employed to analogously manifest the noticeable dependence of the nodal-line pattern on the coupling strength. More importantly, the dispersion relation between the frequency and the wave number for the flexural wave can be precisely determined in the process of utilizing the point interaction to model the modern Chladni plates.

1. Introduction

In the 18th century, Ernst Chladni employed a bow to excite a thin plate, and used small particles to visualize the nodal lines of the resonance [1,2], as shown in Figure 1a. Chladni’s work is the first real attempt to explore acoustic properties through experiments [3,4,5,6]. The resonant modes excited by a bow are usually the eigenmodes of the plate. Nowadays, the modern Chladni plate, as shown in Figure 1b, is excited by an electronically controlled point oscillator to generate the resonant modes systematically. Owing to the point interaction between the plate and the driving source, the nodal patterns observed in the modern Chladni plate have been observed to be entirely different from the eigenmodes of the free plate [7]. The formation of modern Chladni patterns has been explored via Helmholtz equation, RLC networks, and Kirchhoff–Love equation [7,8,9,10,11]. However, there is no theoretical model that includes the action of the point interaction between the plate and the driving source.

The point interaction in terms of the Dirac d-function potential has attracted great interest in the development and advancement of quantum mechanics [12,13,14]. The delta-function potential in two or three dimensions has been widely considered in the literature because of its usefulness in nuclear [15], atomic [16], solid-state [17], and particle physics [18]. Šeba [19] firstly explored the coupling interaction between a two-dimensional integrable billiard and a d-function potential. The strong coupling in quantum billiards might cause the transition from integrable to chaotic features [20,21,22,23,24,25]. As the two- or three-dimensional delta-function potential is plagued with the presence of ultraviolet divergences, renormalization analyses and cutoff regularization [26,27,28,29,30] are usually used to deal with infinities. An alternative approach to remove the difficulty of the divergence is to utilize a truncated basis for expressing the point interaction [31]. From the viewpoint of physical perspectives, the truncation of the basis can be linked to the real systems in which the interaction ordinarily originates from a point-like scatterer with a small, but finite size. So far, the explorations regarding to the point interaction almost all focus on the variation of eigenvalues, and rarely on the structure of eigenfunctions. It is scientifically important to verify the influence of the point interaction on the pattern formation of eigenfunctions through experiment.

In this work, we study the influence of point interaction on the resonant modes in modern Chladni plates. The structure of the paper is as follows. Section 2 presents the theoretical model for point interaction. We employ a truncated basis to explore the point interaction in quantum billiards. Numerical analyses reveal that the point interaction gives rise to the shift of the eigenvalue that causes the superposition of the unperturbed eigenfunctions to form the eigenfunction. The variation in the amplitude patterns with the coupling strength is confirmed to be very significant, as with the nodal-line patterns. Section 3 presents our findings in modern Chladni plates, followed by applications in Section 4 and conclusions in Section 5. We systematically verify that the perceptible dependence of the nodal-line pattern on the coupling strength can be analogously manifested from the point-driven Chladni plates. Furthermore, the dispersion relation between the frequency and the wave number for the flexural wave can be precisely determined in using the point-interaction model to analyze the modern Chladni patterns. Since the dispersion relation contains the information of the flexural rigidity and the Young modulus, the proposed model is believed to be not only pedagogically beneficial to understanding the point interaction in quantum systems but also practically useful for measuring the acoustic properties of plates.

2. Theoretical Model for the Point Interaction

The present analysis follows the model developed by Šeba [19]. The two-dimensional Helmholtz equation for the domain Ω with boundary shape δΩ is given by

(1)$\left({\nabla }_{2D}^2+k_n^2\right){\mathsf{\psi }}_n(r)=0,$

where

(2)${\nabla }_{2D}^2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}$

$r=(x,y)$ and $k_n$ and ${\mathsf{\psi }}_n$, with the indices of $n=1,2,3\cdots$, are the eigenvalues and eigenfunctions, respectively. The issue of quantum billiards is associated with the Dirichlet boundary conditions with ${\mathsf{\psi }}_n=0$. Considering a d-function potential at $r_s$ inside the domain, the Helmholtz equation in Equation (1) can be expressed as

(3)$\left({\nabla }_{2D}^2+\alpha \delta (r-r_s)+k^2\right)\mathsf{\Psi }(r;k)=0,$

where the parameter α is a real number and its absolute value is directly related to the coupling strength. Although the singular potential is not an ordinary function, Equation (3) can be physically analyzed using a truncated basis. The truncation of a basis certainly leads the effective potential to be equivalent to a scatterer with small finite size accompanied with Gibbs-like oscillations. Nevertheless, the finiteness of the scatterer size should be applicable to realistic systems from a physical perspective. In terms of a truncated basis $\{ {\mathsf{\psi }}_n(r) \}$ ($n=1,2,3\cdots$), the functions $\delta (r-r_s)$ and $\mathsf{\Psi }(r)$ in Equation (3) can be approximately expressed as

(4)$\delta (r-r_s)={\displaystyle \sum _{n=0}^N {\mathsf{\psi }}_n^{\ast }(r_s) {\mathsf{\psi }}_n(r)},$

(5)$\mathsf{\Psi }(r;k)={\displaystyle \sum _{n=0}^N a_n(k) {\mathsf{\psi }}_n(r)},$

Substituting Equations (4) and (5) into Equation (3) and using Equation (1), the eigenfunction of Equation (3) can be derived as

(6)$\mathsf{\Psi }(r;k)=-\alpha \mathsf{\Psi }(r_s;k) G(r,r_s;k),$

where $G(r,r_s;k)$ is the Green function given by

(7)$G(r,r_s;k)={\displaystyle \sum _{n=0}^N \frac{{\mathsf{\psi }}_n^{\ast }(r_s) {\mathsf{\psi }}_n(r)}{k^2-k_n^2} }.$

The identification $\delta (r-r_s)\mathsf{\Psi }(r;k)=\delta (r-r_s)\mathsf{\Psi }(r_s;k)$ was used in deriving Equation (6). Setting $r=r_s$ in both sides of Equation (6), the transcendental equation for solving the eigenvalues k is given by $1+\alpha \xi (k)=0$, where $\xi (k)=G(r_s,r_s;k)$ is a meromorphic function. It has been discussed [30,31] that a d-function potential in the Helmholtz equation is problematic in the limit of $N\to \infty$ for a complete basis. The renormalization process from the self-adjoint extension theory [30,31] is employed to treat short-range singularities in a proper manner. Based on the renormalization process, the eigenvalues are determined by

(8)${\alpha }_b^{-1}={\displaystyle \sum _{n=0}^{\infty } {\left| {\mathsf{\psi }}_n(r_s) \right|}^2\left(\frac{1}{k^2-k_n^2}+\frac{k_n^2}{k_n^4+1}\right) },$

where ${\alpha }_b$ is called the renormalized coupling constant. It has been fully discussed [30] that ${\alpha }_b$ can be related to the physical coupling constant α defined in Equation (3) with a truncated basis.

We use a square-shape quantum billiard to illustrate the relationship between the coupling parameter α and the spectrum of eigenvalues. For a free square-shape billiard with the region in $0\le x,y\le L$, the eigenfunctions are given by

(9)${\mathsf{\psi }}_{n,m}(r)=\frac{2}{L}\sin \left(\frac{n\pi }{L}x\right)\sin \left(\frac{m\pi }{L}y\right),$

where n = 1,2,3··· and m = 1,2,3···. Using Equation (9), the meromorphic function $\xi (k)$ for a Dirac’s delta potential at $r_s$ can be expressed as

(10)$\xi (k)={\displaystyle \sum _{n=0}^N {\displaystyle \sum _{m=0}^N \frac{{\left|{\mathsf{\psi }}_{n,m}(r_s)\right|}^2}{k^2-\left[\left(n^2+m^2\right){\pi }^2/L^2\right]}} }.$

The index N is chosen by the criterion that the eigenvalue k is less than $\pi N/2L$. Under this criterion, the eigenvalues k are numerically confirmed to approach some asymptote. Figure 2 shows the calculated result for $\xi (k)$ in Equation (10) as a function of k with $r_s=(L/2,L/2)$. For a given α, the eigenvalues are determined by the condition of $\xi (k)=-1/\alpha$. In the case of $\alpha \to \infty$, the spectrum of eigenvalues corresponds to the zeros $k^{(z)}$ of the function $\xi (k)$, while the spectrum of the free system is given by the poles $k^{(p)}$ of $\xi (k)$. The poles $k^{(p)}$ for the square billiard are given by $k_{n,m}^{(p)} L/\pi =\sqrt{(n^2+m^2)}$. For a positive α increasing from zero to infinity, the eigenvalue can be seen to shift from a pole to a zero by decreasing.

To demonstrate the trend of approach for computing eigenvalues with respect to the index N, the calculated results for ${\left[1+\alpha \xi (k)\right]}^{-1}$ as a function of k for different N increasing from 10 to 50 are shown in Figure 3a for a given $\alpha =5$ and Figure 3b for a given $\alpha =10$. The eigenvalues k can be found to be located at the abrupt changes of the spectrum ${\left[1+\alpha \xi (k)\right]}^{-1}$. The convergence of the numerical eigenvalues can be clearly seen as the value of the index N increases from 10 to 50 in the range of $kL/\pi <20$. Furthermore, the redshift of the eigenvalue subject to the increasing of the coupling strength α can be observed by comparing the results shown in Figure 3a,b for $\alpha =5$ and $\alpha =10$, respectively. We also investigated the influence of the spreading function for the interaction by considering a Gaussian-type function $\exp (-r^2/a^2)$. The numerical results revealed that the results obtained with Gaussian-type function are almost the same as those obtained with d function for $a<L/2N$.

Figure 4a shows the shift in the eigenvalue with the coupling parameter α for a specific case of $k_{7,7}^{(p)}$. The dependence of the wave pattern of the eigenfunction $\mathsf{\Psi }(r;k)$ on the parameter α is shown in Figure 4b. The structure of the nodal line corresponding to the eigenfunction is illustrated in Figure 4c. The wave pattern $\left|\mathsf{\Psi }(r;k)\right|$ at $\alpha =0$ can be seen to be the feature of the square array, corresponding to the original eigenfunction $\left|{\mathsf{\psi }}_{7,7}(r)\right|$; its nodal-line pattern displays the characteristic of the square lattice. The increase in the coupling strength α causes the red shift of the eigenvalue, leading the eigenfunction to be the superposition of the original eigenfunctions ${\mathsf{\psi }}_{n,m}(r)$. The weighting coefficient for the state ${\mathsf{\psi }}_{n,m}(r)$ is inversely proportional to $k^2-{(k_{n,m}^{(p)})}^2$. The closer the value of $k_{n,m}^{(p)}$ to the eigenvalue k, the larger the contribution of the state ${\mathsf{\psi }}_{n,m}(r)$. The eigenvalue can be seen to be significantly decreasing for $\alpha <5$ and to gradually saturate for $\alpha >5$. As shown in Figure 4b,c, the variation in the amplitude patterns with the coupling strength is very conspicuous, as with the nodal-line patterns. Figure 4d shows the contribution of original eigenfunctions ${\mathsf{\psi }}_{n,m}(r)$ in the perturbed eigenfunction $\mathsf{\Psi }(r;k)$ shown in Figure 4c. The number of contributed eigenfunctions ${\mathsf{\psi }}_{n,m}(r)$ in the perturbed eigenfunction $\mathsf{\Psi }(r;k)$ can be clearly seen to increase with the increasing coupling strength α. In the following, we will verify that the perceptible dependence of the nodal-line pattern on the coupling strength can be analogously manifested from the point-driven modern Chladni plates.

3. Exploring the Point Interactions in Modern Chladni Plates

The transverse vibration of a plate is governed by the biharmonic equation [32]

(11)$\left({\nabla }_{2D}^4-k^4\right)\mathsf{\psi }(r)=0,$

where ${\nabla }_{2D}^4$ is the biharmonic operator given by

(12)${\nabla }_{2D}^4 \mathsf{\psi }=\frac{{\partial }^4\mathsf{\psi }}{\partial x^4}+2\frac{{\partial }^4\mathsf{\psi }}{\partial x^2\partial y^2}+\frac{{\partial }^4\mathsf{\psi }}{\partial y^4}.$

Equation (11) can be factorized as

(13)$\left({\nabla }_{2D}^2+k^2\right)\left({\nabla }_{2D}^2-k^2\right)\mathsf{\psi }(r)=0,$

in which the first factor describes propagating dispersive waves, whereas the second one describes evanescent waves. As the aspect ratio of the thickness to the lateral dimension for the studied plate is generally less than 0.02, the vibrating eigenmodes can be approximated with the 2D Helmholtz equation [33]. Even so, the free edge boundary conditions make the problem particularly difficult, as noted by Rayleigh [34]. As a first approximation, one can model the plate as a tightly stretched thin elastic membrane by assuming the constant μ (the ratio of lateral contraction to longitudinal elongation) to be zero, i.e., every point of the circumference is free to move along lines perpendicular to the plane of the plate. Under this approximation, the boundary condition can be simplified as the Neumann boundary conditions with $\partial {\mathsf{\psi }}_n/\partial n=0$ on δΩ. Strictly, the Neumann boundary condition is practically applicable to a stretched membrane, but not to a plate vibrating in virtue of rigidity. Nevertheless, it has been demonstrated [35] that a hypothetical free membrane can be used to deduce some of classical Chladni figures successfully. Based on this confirmation, a free membrane with a point scatter is proposed to model the modern Chladni figures.

A square-shaped aluminum plate with a side length of L = 32 cm and a thickness of h = 1.0 mm was used to perform the modern Chladni experiment. The frequency response was measured by a sine function generator with an automatic scanning system with a resolution of 0.1 Hz. As the conductance of the vibrating system displayed a sharp drop at each resonant frequency, we accurately measured the resonant frequency by exceeding the change in conductance of the system [36]. Chladni patterns at resonant frequencies were exhibited with silica sand of 0.3 mm grain size and were stored with a digital camera. Fixing the center of the plate on the mechanical oscillator, several resonant patterns could be precisely observed by scanning the frequency between 200 and 3000 Hz, as shown in Figure 5a for the nodal-line patterns. All the observed nodal-line patterns are obviously different from the characteristics of eigenfunctions of free plate.

From Equations (6) and (7), the eigenfunctions of the square plate with the point interaction at $r_s$ can be directly expressed as the Green function:

(14)$G(r,r_s;k)={\displaystyle \sum _{n=0}^N {\displaystyle \sum _{m=0}^N \frac{{\tilde{\mathsf{\psi }}}_{n,m}(r_s) {\tilde{\mathsf{\psi }}}_{n,m}(r)}{k^2-\left[\left(n^2+m^2\right){\pi }^2/L^2\right]}} },$

where the free eigenfunctions ${\tilde{\mathsf{\psi }}}_{n,m}(r)$ are given by

(15)${\tilde{\mathsf{\psi }}}_{n,m}(r)=\frac{2}{L}\cos \left(\frac{n\pi }{L}x\right)\cos \left(\frac{m\pi }{L}y\right),$

where n = 0,1,2,3··· and m = 0,1,2,3···. The eigenvalues are determined from the transcendental equation $1+\alpha \xi (k)=0$ with $\xi (k)=G(r_s,r_s;k)$. For a given α, the eigenvalue k can be solved to one-to-one correspond to the free eigenvalue of $k_{n,m}^{(p)} =\pi {(n^2+m^2)}^{1/2}/L$. We numerically confirmed that all experimental nodal-line patterns can be excellently reconstructed by scanning the value of a for the best fit and using $r_s=(0.5L, 0.5L)$. Figure 5b shows the numerical nodal-line patterns corresponding to the experimental results shown in Figure 5a. The value of α for the best reconstruction and the eigenvalue k are specified for each case in Figure 5b. The excellent agreement not only confirms the point-interaction model, but also validates the present analysis for explaining the coupling strength of the modern Chladni figure.

The square plate was experimentally driven by a mechanical oscillator located at $r_s=(0.57L, 0.57L)$ instead of the center to further confirm the theoretical model. Using Equations (14) and (15), numerical calculations were performed under the same procedure of finding the value of α and the corresponding k to make the best reconstruction for the experimental results. Figure 6a,b show experimental and numerical nodal-line patterns, respectively. The off-center excitation can be found to introduce the effect of symmetry breaking into the driven system, leading the nodal-line structure to be more intricate. Numerical patterns can be found to be similar to the experimental results for global structures, except for some differences in the crossing and anti-crossing of nodal lines. Nevertheless, the overall theoretical nodal-line patterns are rather in agreement with the counterparts of experimental figures. This agreement indicates that the proposed model can be universally applied to an arbitrary location for the point interaction. Figure 7 depicts the dependence of the coupling strength α on the resonant frequency f by means of the best reconstruction of experimental patterns. The coupling strength α can be seen to be linearly proportional to the resonant frequency. The value of α is significantly greater than 5 for the frequency higher than 1000 Hz. This result indicates that the point interaction in modern Chladni plate generally approaches the strong-coupling regime. Based on the quantum-acoustic analogy, it has been outlined [37] that the interpretation of Chladni figures of irregularly shaped plates is intimately connected with the quantum mechanics of chaotic billiards. In this work, we originally demonstrate that the modern Chladni plate can be employed to analogously manifest the point interaction in quantum billiards.

4. Applications

More importantly, the point-interaction model can be exploited to precisely determine the dispersion relation between the frequency and the wave number for the flexural wave. From the Kirchhoff–Love plate theory, the dispersion relation between the frequency and the wave number is given by [38]

(16)$f(k)=\frac{C}{2\pi } k^2,$

where $C=\sqrt{D/\rho h}$, D is the flexural rigidity given by

(17)$D=\frac{E{h}^3}{12(1-{\nu }^2)},$

where E is the Young modulus, ν is the Poisson ratio, ρ is the mass density, and h is the thickness of the plate. The material properties for aluminum are as follows: E = 70 GPa, ν = 0.33, and ρ = 2700 kg/m³ [38]. Using these material properties and h = 1 mm, the theoretical coefficient $C/2\pi$ can be calculated to be 0.248. We employed the parabolic formula in Equation (16) to make the best fit to the experimental results obtained in Figure 5 and Figure 6. The fitting coefficient $C/2\pi$ is approximately 0.228, as shown in Figure 8. The experimentally determined coefficient $C/2\pi$ can be found to agree with the theoretical value very well. This good agreement validates that the point-interaction model can be utilized not only to explore the formation of modern Chladni patterns, but also to determine the dispersion relation between the resonant frequency and the wave number for the vibration of plates. As the dispersion relation in Equation (16) includes the information of the flexural rigidity D as well as the Young modulus E, the proposed theoretical model is believed to be practically useful for measuring the acoustic properties of plates.

5. Conclusions

In summary, we have exploited a truncated basis to explore the influence of the point interaction on the eigenvalues and eigenfunctions in quantum billiards. We have numerically found that the shift in the eigenvalue leads the eigenfunction of the point interaction to be the superposition of the unperturbed eigenfunctions. It has been confirmed that the variation in the amplitude patterns with the coupling strength is very conspicuous, as with the nodal-line patterns. We have further verified that the noticeable dependence of the nodal-line pattern on the coupling strength can be analogously manifested from the point-driven Chladni plates. Moreover, utilizing the point-interaction model to analyze the modern Chladni patterns can precisely obtain the dispersion relation between the frequency and the wave number for the flexural wave. The proposed theoretical model is believed to be useful for understanding the point interaction in quantum systems and for measuring the acoustic properties of plates [39,40].

Author Contributions

Conceptualization, K.-F.H. and Y.-F.C.; software, Y.-C.T. and Y.-H.H.; validation, Y.-T.Y. and H.-C.L.; formal analysis, Y.-H.L., Y.-T.Y., and Y.-F.C.; resources, K.-F.H. and Y.-F.C.; writing—original draft preparation, Y.-H.L.; writing—review and editing, H.-C.L. and Y.-F.C.; supervision, K.-F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Ministry of Science and Technology of Taiwan (Contract No. 109-2119-M-009-015-MY3).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data underlying the results presented in this paper are not publicly available at this time, but may be obtained from the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. Experimental setups of (a) classical and (b) modern Chladni plates.

Figure 2. Calculated result for [Forumla omitted. See PDF.] in Equation (8) as a function of k with [Forumla omitted. See PDF.] and N = 50.

Figure 3. Calculated results for [Forumla omitted. See PDF.] as a function of k for different N increasing from 10 to 50 for (a) [Forumla omitted. See PDF.] and (b) [Forumla omitted. See PDF.].

Figure 4. (a) Eigenvalue with the coupling parameter α for a specific case of [Forumla omitted. See PDF.]. (b) Dependence of the wave patterns of the eigenfunctions on the parameter α. (c) Nodal-line patterns corresponding to the wave patterns in (b). (d) Contribution of original eigenfunctions [Forumla omitted. See PDF.] in the perturbed eigenfunction [Forumla omitted. See PDF.].

Figure 5. (a) Experimental nodal-line patterns for fixing the plate at [Forumla omitted. See PDF.](on the center). (b) Numerical nodal-line patterns corresponding to the experimental results.

Figure 6. (a) Experimental nodal-line patterns for fixing the plate at [Forumla omitted. See PDF.] (off the center). (b) Numerical nodal-line patterns corresponding to the experimental results.

Figure 7. Dependence of the coupling strength on the resonant frequency by means of the best reconstruction of experimental patterns.

Figure 8. Dispersion relation between the frequency and the wave number determined from the experimental reconstruction shown in Figure 5 and Figure 6.

References

1. Chladni, E.F.F. Entdeckungen über die Theorie des Klanges; Breitkopf und Härtel: Leipzig, Germany, 1787.

2. Chladni, E.F.F. Die Akustik; Breitkopf und Härtel: Leipzig, Germany, 1802.

3. Gough, C. The violin: Chladni patterns, plates, shells and sounds. Eur. Phys. J. Spec. Top.; 2007; 145, pp. 77-101. [DOI: https://dx.doi.org/10.1140/epjst/e2007-00149-0]

4. Dorrestijn, M.; Bietsch, A.; Açıkalın, T.; Raman, A.; Hegner, M.; Meyer, E.; Gerber, C. Chladni figures revisited based on nanomechanics. Phys. Rev. Lett.; 2007; 98, 026102. [DOI: https://dx.doi.org/10.1103/PhysRevLett.98.026102]

5. Wood, C.D.; Cunningham, J.E.; O’Rorke, R.; Wälti, C.; Linfield, E.H.; Davies, A.G.; Evans, S.D. Formation and manipulation of two-dimensional arrays of micron-scale particles in microfluidic systems by surface acoustic waves. Appl. Phys. Lett.; 2009; 94, 054101. [DOI: https://dx.doi.org/10.1063/1.3076127]

6. Jenny, H. Cymatics: A Study of Wave Phenomena and Vibration; Macromedia: Newmarket, NH, USA, 2001.

7. Tuan, P.H.; Wen, C.P.; Chiang, P.Y.; Yu, Y.T.; Liang, H.C.; Huang, K.F.; Chen, Y.F. Exploring the resonant vibration of thin plates: Reconstruction of Chladni patterns and determination of resonant wave numbers. J. Acoust. Soc. Am.; 2015; 137, 2113. [DOI: https://dx.doi.org/10.1121/1.4916704]

8. Tuan, P.H.; Liang, H.C.; Tung, J.C.; Chiang, P.Y.; Huang, K.F.; Chen, Y.F. Manifesting the evolution of eigenstates from quantum billiards to singular billiards in the strongly coupled limit with a truncated basis by using RLC networks. Phys. Rev. E; 2015; 92, 062906. [DOI: https://dx.doi.org/10.1103/PhysRevE.92.062906] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26764773]

9. Tuan, P.H.; Tung, J.C.; Liang, H.C.; Chiang, P.Y.; Huang, K.F.; Chen, Y.F. Resolving the formation of modern Chladni figures. Europhys. Lett.; 2015; 111, 64004. [DOI: https://dx.doi.org/10.1209/0295-5075/111/64004]

10. Tuan, P.H.; Lai, Y.H.; Wen, C.P.; Huang, K.F.; Chen, Y.F. Point-driven modern Chladni figures with symmetry breaking. Sci. Rep.; 2018; 8, pp. 1-13. [DOI: https://dx.doi.org/10.1038/s41598-018-29244-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30022128]

11. Tuan, P.H.; Wen, C.P.; Yu, Y.T.; Liang, H.C.; Huang, K.F.; Chen, Y.F. Exploring the distinction between experimental resonant modes and theoretical eigenmodes: From vibrating plates to laser cavities. Phys. Rev. E; 2014; 89, 022911. [DOI: https://dx.doi.org/10.1103/PhysRevE.89.022911]

12. Albeverio, S.; Gesztesy, F.; Høegh-Krohn, R.; Holden, H. Solvable Models in Quantum Mechanics; 2nd ed. AMS Chelsea Publishing: Providence, RI, USA, 2004.

13. Schmidt, A.G.M.; Cheng, B.K.; da Luz, M.G.E. Green functions for generalized point interactions in one dimension: A scattering approach. Phys. Rev. A; 2002; 66, 062712. [DOI: https://dx.doi.org/10.1103/PhysRevA.66.062712]

14. Arnbak, H.; Christiansen, P.L.; Gaididei, Y.B. Non-relativistic and relativistic scattering by short-range potentials. Philos. Trans. R. Soc. A; 2011; 369, pp. 1228-1244. [DOI: https://dx.doi.org/10.1098/rsta.2010.0330] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21320914]

15. Kruppa, A.T.; Varga, K.; Révai, J. Local realizations of contact interactions in two-and three-body problems. Phys. Rev. C; 2001; 63, 064301. [DOI: https://dx.doi.org/10.1103/PhysRevC.63.064301]

16. Demkov, Y.; Ostrovskii, V.N. Zero-Range Potentials and Their Applications in Atomic Physics; Plenum: New York, NY, USA, 1989.

17. Doniach, S.; Sondheimer, E.H. Green’s Functions for Solid State Physicists; Benjamin: Reading, MA, USA, 1974.

18. Thorn, C. Quark confinement in the infinite-momentum frame. Phys. Rev. D; 1979; 19, 639. [DOI: https://dx.doi.org/10.1103/PhysRevD.19.639]

19. Šeba, P. Wave Chaos in Singular Quantum Billiard. Phys. Rev. Lett.; 1990; 64, 1855. [DOI: https://dx.doi.org/10.1103/PhysRevLett.64.1855]

20. Shigehara, T. Conditions for the appearance of wave chaos in quantum singular systems with a pointlike scatterer. Phys. Rev. E; 1994; 50, 4357. [DOI: https://dx.doi.org/10.1103/PhysRevE.50.4357]

21. Šeba, P.; Exner, P. Point interactions in two and three dimensions as models of small scatterers. Physics Letters A; 1996; 222, pp. 1-4.

22. Bogomolny, E.; Gerland, U.; Schmit, C. Singular statistics. Phys. Rev. E; 2001; 63, 036206. [DOI: https://dx.doi.org/10.1103/PhysRevE.63.036206]

23. Berkolaiko, G.; Keating, J.P.; Winn, B. Intermediate wave function statistics. Phys. Rev. Lett.; 2003; 91, 134013. [DOI: https://dx.doi.org/10.1103/PhysRevLett.91.134103] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14525308]

24. Tudorovskiy, T.; Kuhl, U.; Stöckmann, H.-J. Singular statistics revised. New J. Phys.; 2010; 12, 123021. [DOI: https://dx.doi.org/10.1088/1367-2630/12/12/123021]

25. Rudnick, Z.; Ueberschär, H. Statistics of wave functions for a point scatterer on the torus. Commun. Math. Phys.; 2012; 316, pp. 763-782. [DOI: https://dx.doi.org/10.1007/s00220-012-1556-2]

26. Jackiw, R.M.A.B. Bég Memorial Volume; Ali, A.; Hoodbhoy, P. World Scientific: Singapore, 1991.

27. Phillips, D.R.; Beane, S.R.; Cohen, T.D. Nonperturbative regularization and renormalization: Simple examples from nonrelativistic quantum mechanics. Ann. Phys.; 1998; 263, pp. 255-275. [DOI: https://dx.doi.org/10.1006/aphy.1997.5771]

28. Henderson, R.J.; Rajeev, S.G. Renormalized contact potential in two dimensions. J. Math. Phys.; 1998; 39, pp. 749-759. [DOI: https://dx.doi.org/10.1063/1.532350]

29. Adhikari, S.K.; Frederico, T. Renormalization group in potential scattering. Phys. Rev. Lett.; 1995; 74, 4572. [DOI: https://dx.doi.org/10.1103/PhysRevLett.74.4572] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10058544]

30. Weaver, R.L.; Sornette, D. Range of spectral correlations in pseudointegrable systems: Gaussian-orthogonal-ensemble statistics in a rectangular membrane with a point scatterer. Phys. Rev. E; 1995; 52, 3341. [DOI: https://dx.doi.org/10.1103/PhysRevE.52.3341] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9963809]

31. Shigehara, T.; Cheon, T. Wave chaos in quantum billiards with a small but finite-size scatterer. Phys. Rev. E; 1996; 54, 1321. [DOI: https://dx.doi.org/10.1103/PhysRevE.54.1321]

32. Chakraverty, S. Vibration of Plates; CRC Press: Taylor and Francis Group: Boca Raton, FL, USA, 2009.

33. Ventsel, E.; Krauthammer, T. Thin Plates and Shells; Dekker: New York, NY, USA, 2004; pp. 1-14.

34. Rayleigh, L. Theory of Sound; Dover: New York, NY, USA, 1945; pp. 367-380.

35. Waller, M.D. Vibrations of free square plates: Part I. Normal vibrating modes. Proc. Phys. Soc.; 1939; 51, pp. 831-844. [DOI: https://dx.doi.org/10.1088/0959-5309/51/5/312]

36. Snowdon, J.C. Forced vibration of internally damped rectangular and square plates with simply supported boundaries. J. Acoust. Soc. Am.; 1974; 56, pp. 1177-1184. [DOI: https://dx.doi.org/10.1121/1.1903404]

37. Stöckmann, H.-J. Quantum Chaos. An Introduction; University Press: Cambridge, UK, 1999.

38. Van Vlack, L. Elements of Material Science and Engineering; Addison-Wesley: Reading, MA, USA, 1980.

39. Luo, Y.; Feng, R.; Li, X.D.; Liu, D.H. A simple approach to determine the mode shapes of Chladni plates based on the optical lever method. Eur. J. Phys.; 2019; 40, 065001. [DOI: https://dx.doi.org/10.1088/1361-6404/ab2e2b]

40. Lei, J.; Cheng, F.; Liu, G.; Li, K.; Gou, Z. Dexterous formation of unconventional Chladni patterns using standing bulk acoustic waves. Appl. Phys. Lett.; 2020; 117, 184101. [DOI: https://dx.doi.org/10.1063/5.0032304]