Exploring the Origin of Maximum Entropy States Relevant to Resonant Modes in Modern Chladni Plates

The resonant modes generated from the modern Chladni experiment are systematically confirmed to intimately correspond to the maximum entropy states obtained from the inhomogeneous Helmholtz equation for the square and equilateral triangle plates. To investigate the origin of maximum entropy states, the inhomogeneous Helmholtz equation is modified to consider the point interaction coming from the driving oscillator. The coupling strength associated with the point interaction is characterized by a dimensionless factor α. The δ potential of the point interaction is numerically modelled by a truncated basis with an upper index N. The asymptotic behavior for the upper index N is thoroughly explored to verify that the coupling strength of α = 1.0 can make the theoretical resonant modes agree excellently with the maximum entropy states as N→∞. It is further authenticated that nearly the same resonant modes can be obtained by using a larger coupling strength α when a smaller upper index N is exploited in the calculation.


Introduction
Ernst Chladni in the 18th century demonstrated a classical experiment related to acoustic phenomenon by means of vibrating a thin plate with a bow and manifesting the resonant nodal-line pattern with small sands [1,2]. Chladni figures have stimulated numerous scientific explorations in physics including quantum chaos [3], self-organization of granular media [4,5], microscale acoustofluidics [6,7], and pattern formation [8,9]. The classical Chladni figures obtained by a bow are found to considerably different from the resonant nodal-line patterns in the modern Chladni plate which is driven by a mechanical oscillator whose vibrating frequency can be precisely controlled to generate the resonant modes. Thanks to the advantages of high precision, reproducibility, and simple performing, the modern Chladni plate is a promising experiment for creating advanced applications such as automated patterning of micro-objects [10,11], non-contaminated positioning of biomolecules [12,13], and sorting different particles [14,15].
The formation of modern Chladni patterns has been successfully elucidated with the inhomogeneous Helmholtz equation and the maximum entropy principle [16][17][18][19][20]. The coincidence of the experimental resonant frequencies and the maximum entropy states can be comprehended from the principle of energy equipartition in statistical mechanics. Nevertheless, it is of great importance to explore the origin of the resonant modes in modern Chladni plates for progressive developments.
In this work, we systematically identify the intimate correspondence between the resonant modes generated from the modern Chladni experiment and the maximum entropy states solved from the inhomogeneous Helmholtz equation for the square as well as equilateral triangle plates. The point interaction arising from the driving oscillator is considered in the inhomogeneous Helmholtz equation to explore the origin of maximum entropy states. A dimensionless factor α is used to characterize the coupling strength associated with the point interaction. Moreover, a truncated basis with an upper index N is employed to model the δ potential of the point interaction. We numerically verify that as N → ∞ , the coupling strength of α = 1.0 can make the theoretical resonant modes agree excellently with the maximum entropy states. We further confirm that when a smaller upper index N is used in the calculation, a larger α can be effectively exploited to obtain almost the same theoretical resonant modes. The developed model with the point interaction is believed to be useful for creating novel applications. The whole structure of this paper includes the experimental results for resonant frequencies and modes of modern Chladni plates, the theorical maximum entropy states, and the point-interaction model for exploring the origin of maximum entropy states relevant to modern Chladni plates. Figure 1 depicts the experimental scheme for the modern Chladni plate driven with an electronically controlled mechanical oscillator [16]. A sinusoidal function generator with an automatic scanning system was designed to measure the spectrum of resonant frequencies with a resolution of 0.1 Hz. Two different plates were used to investigate the resonant Chladni figures. One is a square plate with a side length of 320 mm; the other is an equilateral triangle plate with a side length of 365 mm. Two thin plates were manufactured with aluminum sheets with a thickness of 1 mm. The center of the thin plate was fixed at the mechanical oscillator with a screw supporter. The mechanical oscillator was excited with an amplified sinusoidal voltage source for the frequency ranging from 200 to 5000 Hz. A digital galvanometer was connected in series to the mechanical oscillator to measure the electric current through the vibrating system. Since the resonance causes the electric resistance of the vibrating system to increase abruptly, the resonant frequencies can be precisely determined via probing the change in the electric resistance of the system [16]. We employed silica sand of 0.3 mm grain size to exhibit Chladni figures at resonant frequencies and exploited a digital camera to store the experimental images.

Frequency Spectrum and Resonant Modes
Entropy 2022, 24, x FOR PEER REVIEW 2 of 15 considered in the inhomogeneous Helmholtz equation to explore the origin of maximum entropy states. A dimensionless factor α is used to characterize the coupling strength associated with the point interaction. Moreover, a truncated basis with an upper index N is employed to model the δ potential of the point interaction. We numerically verify that as N → ∞ , the coupling strength of α = 1.0 can make the theoretical resonant modes agree excellently with the maximum entropy states. We further confirm that when a smaller upper index N is used in the calculation, a larger α can be effectively exploited to obtain almost the same theoretical resonant modes. The developed model with the point interaction is believed to be useful for creating novel applications. The whole structure of this paper includes the experimental results for resonant frequencies and modes of modern Chladni plates, the theorical maximum entropy states, and the point-interaction model for exploring the origin of maximum entropy states relevant to modern Chladni plates. Figure 1 depicts the experimental scheme for the modern Chladni plate driven with an electronically controlled mechanical oscillator [16]. A sinusoidal function generator with an automatic scanning system was designed to measure the spectrum of resonant frequencies with a resolution of 0.1 Hz. Two different plates were used to investigate the resonant Chladni figures. One is a square plate with a side length of 320 mm; the other is an equilateral triangle plate with a side length of 365 mm. Two thin plates were manufactured with aluminum sheets with a thickness of 1 mm. The center of the thin plate was fixed at the mechanical oscillator with a screw supporter. The mechanical oscillator was excited with an amplified sinusoidal voltage source for the frequency ranging from 200 to 5000 Hz. A digital galvanometer was connected in series to the mechanical oscillator to measure the electric current through the vibrating system. Since the resonance causes the electric resistance of the vibrating system to increase abruptly, the resonant frequencies can be precisely determined via probing the change in the electric resistance of the system [16]. We employed silica sand of 0.3 mm grain size to exhibit Chladni figures at resonant frequencies and exploited a digital camera to store the experimental images.    Figure 2a shows the measured spectrum for the variation of the resistances as a function of the driving frequency in the vibrating square plate. The resonant frequencies can be precisely identified from the peaks of the frequency spectrum. The nodal-line patterns of the resonant modes are shown completely in Figure 2b. All nodal-line patterns can be found to be obviously distinct from the chessboard structures of eigenmodes. Experimental results for the frequency spectrum and individual resonant nodal-line patterns obtained with the equilateral triangle plate are shown completely in Figure 3a,b. Like the results for the square plates, resonant nodal-line patterns are dissimilar to the nodal-line structures of eigenmodes.

Maximum Entropy States of the Vibrating Plate
The transverse vibration of a plate is governed by the biharmonic equation [21] ( ) ( , ) x y = r , and k is the wave number. Equation (1) can be factorized as

Maximum Entropy States of the Vibrating Plate
The transverse vibration of a plate is governed by the biharmonic equation [21] where r = (x, y), and k is the wave number. Equation (1) can be factorized as in which the first factor describes propagating dispersive waves, whereas the second one describes evanescent waves. Since the aspect ratio h/L is generally less than 0.02, the vibrating eigenmodes can be approximated with the 2D Helmholtz equation [22], where h and L are the thickness and lateral size of the studied plate, respectively. Even so, the free edge boundary conditions make the problem particularly difficult, as noted by Rayleigh [23]. As a first approximation, one can model the plate as a tightly stretched thin elastic membrane by assuming the constant m (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 ∂ψ/∂n = 0 on the periphery. 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 [24] that a hypothetical free membrane can be used to deduce some of the classical Chladni figures successfully. Based on this confirmation, a free membrane with a point scatter is proposed to model the modern Chladni figures. The two-dimensional Helmholtz equation for the domain Ω with boundary shape dΩ is given by where k n and ψ n with the indices of n = 0, 1, 2, 3 · · · are the eigenvalues and eigenfunctions, respectively. Considering the exciter as a point source, the inhomogeneous Helmholtz equation for characterizing the modern Chladni plate is given by where r s is the position of the exciter and A is the amplitude of the driving force. Using the expansion of the eigenfunctions {ψ n (r)} (n = 1, 2, 3 · · · ), the source function δ(r − r s ) and the response function Ψ(r, r s ; k) in Equation (5) can be expressed as Ψ(r, r s ; k) = lim N→∞ N ∑ n=0 a n (r s ; k) ψ n (r), Substituting Equations (6) and (7) into Equation (5) and using Equation (4), the coefficients a n (r s ; k) in the response function Ψ(r, r s ; k) can be derived as The response function is given by Ψ(r, r s ; k) = A G(r, r s ; k), where G(r, r s ; k) is the wellknown Green function Practically, a finite upper index N is chosen for numerical analysis. The use of a finite index N physically corresponds to a truncated basis. The truncation of a basis is relevant to a realistic system with a finite range of frequency response. The influence of the value of the index N is discussed in all the following numerical analyses. The resonant frequencies in modern Chladni experiment have been found to be intriguingly correlated to the maximum entropy states [16,19]. From the Shannon theory and in terms of the coefficients a n (r s ; k), the information entropy of the response function Ψ(r, r s ; k) can be evaluated by [25,26]  where the probability of the eigenmode ψ n (r) is given by For bipartite states, the information entropy is often used to evaluate the degree of spatial entanglement [27][28][29][30]. Here we demonstrate complete results for the correlation between the resonant frequencies and the maximum entropy states for the square and equilateral triangle plates. For a square-shape plate with the region in 0 ≤ x, y ≤ L, the eigenfunctions are given by The eigenvalue corresponding to the eigenfunction ψ n,m (r) is given by Using Equations (8), (10) and (11), the probability of the eigenfunction ψ n,m (r) in the information entropy for a given upper index N is given by In terms of the probabilities p n,m (r s ; k), the information entropy can be expressed as For an equilateral triangle plate with vertices at (0, 0), (L/2, √ 3L/2), and (−L/2, √ 3L/2), the eigenmodes for the driving source with even symmetry can be expressed as [31] ψ n,m (r) = 16 for m ≥ 2n. The eigenvalue for the eigenmode ψ n,m (r) is given by The information entropy for the equilateral triangle plate with an upper index N is given by The information entropy can be exploited to evaluate the effective number N e f f of participated eigenfunctions in the response function, given by N e f f = exp[S(r s ; k)]. It has been verified that the resonant wave numbers can be properly attained from local maxima of the spectrum N eff [16,19].
Based on thorough computation, the entropy S(r s ; k) is confirmed to be nearly independent of the choice of the upper index N on condition that N > kL is satisfied. Since the range of kL for investigation is less than 60, the results obtained with N = 100 are used for convenience. Figure 4a depicts the calculated results for exp[S(r s ; k)] as a function of k with r s = (L/2, L/2). The nodal-line patterns corresponding to the local maxima are shown completely in Figure 4b. The nodal-line patterns can be one-to-one seen to agreeably resemble the experimental results for the resonant modes shown in Figure 2b. Not only the square plate, the same correspondence between the resonant frequencies and the maximum entropy states can be also confirmed for the equilateral triangle plate, as shown in Figure 5a,b for calculated results with r s = (0, L/ √ 3). verified that the resonant wave numbers can be properly attained from local maxima of the spectrum Neff [16,19]. Based on thorough computation, the entropy ( ; ) s S k r is confirmed to be nearly independent of the choice of the upper index N on condition that N kL > is satisfied. Since the range of kL for investigation is less than 60, the results obtained with 100 N = are used for convenience. Figure  . The nodal-line patterns corresponding to the local maxima are shown completely in Figure 4b. The nodal-line patterns can be one-to-one seen to agreeably resemble the experimental results for the resonant modes shown in Figure 2b. Not only the square plate, the same correspondence between the resonant frequencies and the maximum entropy states can be also confirmed for the equilateral triangle plate, as shown in Figure 5a,b for calculated results with   , E is the Young's modulus, ν Poisson ratio, ρ is the mass density, and h is the thickness of plate. By using the m properties of E = 70 GPa, ν = 0.33, ρ = 2700 kg/m 3 [33], and h = 1 mm, the theoretica ficient / 2 C π can be found to be 0.248. From the correspondence between the re frequencies and the maximum entropy states, the coefficient / 2 C π can be found approximately 0.226 and 0.220 for the square and equilateral triangle plates, respec The good agreement confirms that the information of the flexural rigidity D a Young's modulus E can be deduced by using the maximum entropy principle as w the dispersion relation.
The maximum entropy principle has been developed to analyze the collective ior in multimode systems such as maximum emission for lasers [34], self-organizat Moreover, the dispersion relation for the flexural wave of the plate can be directly deduced from the correspondence between the resonant frequencies f and the resonant wave numbers k of maximum entropy states. From the Kirchhoff-Love plate theory [32], the dispersion relation between f and k is given by f (k) = Ck 2 /2π, where C = D/ρh, D is the flexural rigidity given by D = Eh 3 /[12(1 − ν 2 )], E is the Young's modulus, ν is the Poisson ratio, ρ is the mass density, and h is the thickness of plate. By using the material properties of E = 70 GPa, ν = 0.33, ρ = 2700 kg/m 3 [33], and h = 1 mm, the theoretical coefficient C/2π can be found to be 0.248. From the correspondence between the resonant frequencies and the maximum entropy states, the coefficient C/2π can be found to be approximately 0.226 and 0.220 for the square and equilateral triangle plates, respectively. The good agreement confirms that the information of the flexural rigidity D and the Young's modulus E can be deduced by using the maximum entropy principle as well as the dispersion relation.
The maximum entropy principle has been developed to analyze the collective behavior in multimode systems such as maximum emission for lasers [34], self-organization for complex systems [35], wave function localization for disordered systems [36], and phase transitions for open quantum systems [37]. However, the coincidence of the resonant frequencies and the maximum entropy states can be comprehended from the principle of energy equipartition in statistical mechanics. Nevertheless, it is of great importance to disclose the origin of the resonant modes in modern Chladni plates for developing further applications.

Coupling Interaction between Source and Plate
Considering the mechanical oscillator to play as not only a driving source but also a coupling interaction, the characteristic equation in Equation (5) is modified as where the parameter α is a real number. The term α δ(r − r s )Ψ α (r, r s ; k) in Equation (20) is to take account of the coupling effect between the mechanical oscillator and the thin plate. The coupling interaction comes from the reaction force of the thin plate subjected to the driving force. The dimensionless parameter α is directly related to the coupling strength. Intuitionally, the coupling strength is proportional to the relative tightness of the fixed force between the plate and the oscillator. The point interaction [38][39][40] has attracted much attention in the research of the nuclear [41], atomic [42], solid-state [43], and particle physics [44]. Furthermore, Šeba [45] studied the coupling interaction of a delta-function potential in two-dimensional integrable billiards to verify that the strong coupling can cause the transition from integrable to chaotic feature [46][47][48][49][50][51][52][53]. On the other hand, a novel method for analytically solving field distribution in two-dimensional inhomogeneous waveguides is presented in Ref. [54].
Using the expansion forms in Equations (6) and (7), the response function Ψ α (r, r s ; k) in Equation (20) can be derived as Ψ α (r, r s ; k) = A G(r, r s ; k) 1 + αG(r s , r s ; k) where G(r s , r s ; k) is a meromorphic function. Note that the identification of δ(r − r s ) Ψ α (r, r s ; k) = δ(r − r s )Ψ α (r s , r s ; k) has been used in the derivation of Equation (21). The norm |Ψ α (r s , r s ; k)| indicates the response amplitude of the plate at the driving position r s for a given wave number k. Consequently, the local maxima of |Ψ α (r s , r s ; k)| as a function of the driving wave number k can be identified to correspond to the resonant wave numbers for the modern Chladni plate locally excited at the position r s . From Equation (21), the amplitude of the response function |Ψ α (r s , r s ; k)| can be expressed as Equation (22) indicates that if G(r s , r s ; k) = 0, the resonant wave numbers corresponding to the local maxima of |Ψ α (r s , r s ; k)| are determined by the transcendental equation 1 + α G(r s , r s ; k) = 0. Note that the overall structure of the distribution |Ψ α (r s , r s ; k)| is independent of the driving amplitude A. Without the coupling effect for the case α = 0, the resonant wave numbers are determined by the poles of G(r s , r s ; k), which are exactly the eigenvalues k n of the free plates. For a given coupling factor α, the resonant wave numbers are determined from G(r s , r s ; k) = −1/α. The key issue to validate the present model is to verify that the distribution |Ψ α (r s , r s ; k)| can be consistent with the experimental spectrum as well as the maximum entropy states.
Unlike the numerical analyses for information entropy, the calculated results for the distribution |Ψ α (r s , r s ; k)| are found to be not only dependent on the coupling factor α but also on the upper index N. Nevertheless, the asymptotic behavior for the index N can be attained from the computation. The dependence of the distribution |Ψ α (r s , r s ; k)| on the coupling factor α is first discussed. Figure 6 shows the calculated result for |Ψ α (r s , r s ; k)| by using Equation (22) as a function of k for the square plate with r s = (L/2, L/2) and N = 50 for several different coupling factors α = 1.0, 1.5, 2.0 and 2.5. The coupling factor α can be clearly seen to cause the redshifts of the resonant wave numbers. The larger the coupling factor, the greater the redshift. On the other hand, the dependence of the distribution |Ψ α (r s , r s ; k)| on the upper index N is shown in Figure 7 for the square plate with a fixed α = 1.0. The overall distribution |Ψ α (r s , r s ; k)| can be found to be a little redshifted with increasing the upper index N.     The resonant wave numbers are confirmed to approach some asymptote for increasing the upper index N. Numerical analyses reveal that as N → ∞ , the coupling factor of α = 1.0 can make the distribution |Ψ α (r s , r s ; k)| agree the best with the experimental spectrum as well as the maximum entropy states. In practice, when a smaller upper index N is used in the calculation, a larger value for the coupling factor α can be effectively utilized to lead the distribution |Ψ α (r s , r s ; k)| to be nearly the same, as shown in Figure 8. The resonant peaks of the distributions |Ψ α (r s , r s ; k)| for different upper index N can be seen to excellently coincide together. More importantly, all resonant peaks agree very well with the maximum entropy states. In Figure 8, the information entropy is multiplied by a factor of 5 for convenience of displaying the comparison. Figure 9 shows the numerical result for the relationship between the upper index N and the effective coupling factor α for the distributions |Ψ α (r s , r s ; k)| with the best fit to the maximum entropy states. By using the deduced dispersion relation in Section 3, the distribution |Ψ α (r s , r s ; k)| with N = 400 and α = 2.0 is plotted in Figure 10 to compare with the experimental spectrum. For clear comparison, the experimental spectrum is multiplied by a factor of 100. The theoretical resonant spectrum generally agrees very well with the experimental result. A little discrepancy is mainly attributed to the simplified model for the dispersion relation between f and k.
The resonant peaks of the distributions | ( , ; )| s s k α Ψ r r for different upper index N can be seen to excellently coincide together. More importantly, all resonant peaks agree very well with the maximum entropy states. In Figure 8, the information entropy is multiplied by a factor of 5 for convenience of displaying the comparison. Figure 9 shows the numerical result for the relationship between the upper index N and the effective coupling factor α for the distributions | ( , ; )| s s k α Ψ r r with the best fit to the maximum entropy states. By using the deduced dispersion relation in Section 3, the distribution | ( , ; )| s s k α Ψ r r with 400 N = and 2.0 α = is plotted in Figure 10 to compare with the experimental spectrum. For clear comparison, the experimental spectrum is multiplied by a factor of 100. The theoretical resonant spectrum generally agrees very well with the experimental result. A little discrepancy is mainly attributed to the simplified model for the dispersion relation between f and k.  The resonant peaks of the distributions | ( , ; )| s s k α Ψ r r for different upper index N can be seen to excellently coincide together. More importantly, all resonant peaks agree very well with the maximum entropy states. In Figure 8, the information entropy is multiplied by a factor of 5 for convenience of displaying the comparison. Figure 9 shows the numerical result for the relationship between the upper index N and the effective coupling factor α  Figure 10 to compare with the experimental spectrum. For clear comparison, the experimental spectrum is multiplied by a factor of 100. The theoretical resonant spectrum generally agrees very well with the experimental result. A little discrepancy is mainly attributed to the simplified model for the dispersion relation between f and k.   Figure 12. For clear comparison, the experimental spectrum is multiplied by a factor of 100. Once again, the theoretical resonant spectrum is in good agreement with the experimental result. The good agreement further confirms that the point-interaction model can be used to manifest the formation of modern Chladni patterns.  With the best fit to the maximum entropy states for the equilateral triangle plate, the relationship between the upper index N and the effective coupling factor α for the distributions |Ψ α (r s , r s ; k)| is intriguingly found to be almost the same as the case of square plate. Figure 11 depicts the calculated results for the distributions |Ψ α (r s , r s ; k)| with four different best pairs of N and α to compare with the distribution of the information entropy. For convenience of displaying the comparison, the information entropy is multiplied by a factor of 10. Just like the plot in Figure 8, the peaks of the distributions |Ψ α (r s , r s ; k)| can be clearly seen to coincide together and to agree quite well with the maximum entropy states. By using the deduced dispersion relation in Section 3, the comparison between the experimental spectrum and the distribution |Ψ α (r s , r s ; k)| with N = 400 and α = 2.0 is shown in Figure 12. For clear comparison, the experimental spectrum is multiplied by a factor of 100. Once again, the theoretical resonant spectrum is in good agreement with the experimental result. The good agreement further confirms that the point-interaction model can be used to manifest the formation of modern Chladni patterns.

Conclusions
In summary, we have demonstrated that the resonant modes in the modern Chladni experiment can be completely reconstructed with the maximum entropy states solved from the inhomogeneous Helmholtz equation for the square and equilateral triangle plates. We have judiciously considered the point interaction coming from the driving oscillator to explore the relevance between the maximum entropy states and experimental resonant modes. A dimensionless factor α and a truncated basis with an upper index N were used to characterize the coupling strength and the δ potential for the point interaction. We have employed an asymptotic way to verify that the coupling strength of

Conclusions
In summary, we have demonstrated that the resonant modes in the modern Chladni experiment can be completely reconstructed with the maximum entropy states solved from the inhomogeneous Helmholtz equation for the square and equilateral triangle plates. We have judiciously considered the point interaction coming from the driving oscillator to explore the relevance between the maximum entropy states and experimental resonant modes. A dimensionless factor α and a truncated basis with an upper index N were used to characterize the coupling strength and the δ potential for the point interaction. We have employed an asymptotic way to verify that the coupling strength of α = 1.0 can make the theoretical resonant modes agree very well with the maximum entropy states as N → ∞ . We have further validated that a larger coupling strength α can lead to nearly the same resonant spectrum when a smaller upper index N is used in the calculation.