Resonant Phenomena of Faraday Waves: A Classification According to New and Existing Terms
Rona Geffen1, Amira Val Baker1*, Daniela Gentile1, and Paul Oomen1
1*The Works Research Institute, H-1044 Budapest, Hungary
*Correspondence: avbaker@theworks.info
ORCID ID (in order of authors): 0000-0002-8241-9870, 0000-0002-5436-4825, 0000-0002-8340-0397, 0000-0003-3001-4385
Keywords: Aluminum; Faraday Waves; Pattern Morphology; Resonant Phenomena; Cymatics; Vibrational Modes; CymaScope; Taxonomy
Submitted: October 9, 2023
Revised: May 5, 2025
Accepted: June 5, 2025
Published: September 9, 2025
Abstract
Nonlinear standing waves known as Faraday waves are formed in vibrating liquids enclosed by a container. To further our understanding of the relation between frequency and wave propagation within spatial boundaries, we investigated Faraday wave phenomena utilizing the CymaScope instrument, which works by transposing sonic periodicities to water wavelet periodicities. This paper presents the resulting observations and suggests a unified language for describing sound-induced resonant wave phenomena.
Introduction
Sound-related vibratory phenomena are observed in the formation of geometric patterns of either Faraday waves formed in a vertically vibrated liquid-filled cuvette or Chladni figures formed by vibrating grains on a flat medium. The work presented here investigates Faraday waves obtained using the CymaScope instrument that was invented and developed by John Stuart Reid (Web ref. 1). The trade name, “CymaScope,” is derived by simply combining “kyma,” meaning “wave,” with the ancient Greek “skopós,” meaning “watcher.”
Faraday waves, first reported by Michael Faraday in an appendix to a document in 1831, are standing waves resulting from parametric instability (where a change in any parameter significantly changes the outcome) of the capillary gravity waves (where both gravity and surface tension are factors in creating the waves) (Faraday, 1831; Douady and Fauve,1988; Fauve, 1998; Gallaire and Brun, 2017). This instability is termed the “Faraday instability” and relates to standing waves created at a critical amplitude and driving frequency of a vertically vibrated liquid-filled cuvette. When the parameters that create the instability are graphed, the shape of the unstable region in the graph can manifest in the shape of a tongue, known as an instability tongue (Benjamin and Ursell, 1954; Douady, 1990; Bechhoefer and Johnson, 1996; Perlin and Schultz, 2000; Gallaire and Brun, 2017).
Chladni figures, also referred to as nodal patterns, are observed as emanating symmetric patterns following excitation of a flat medium sprinkled with small grains (usually sand). The patterns emerge as the granular substance sprinkled on the plate collects at the nodal points, which are the points of equilibrium or close to zero amplitude. For example, see Figure 1. In the wake of previous experiments, conducted in 1630 by Galileo Galilei (Galilei, 1638) and in 1680 by Robert Hooke (Andrade, 1950), the German musician and physicist Ernest Chladni observed the vibrations of matter through various experiments. Chladni documented this phenomenon of sound frequencies organizing matter in symmetric patterns and demonstrated that these patterns get more intricate as the frequency increases (Chladni, 1787). Mary Desiree Waller continued the research into Chladni figures and in her book, “Chladni Figures, a study in symmetry,” she reported observing the symmetry in Chladni figures and was the first person to deduce a law of symmetry of vibrating two-dimensional systems (Waller, 1938, 1961). Chladni figures are primarily determined by the excitation frequency, the boundary conditions and the excited medium (Escaler and De La Torre, 2018; Tuan et al., 2015), inherently leading to a length-wavelength relation to the pattern’s symmetry, allowing for the successful theoretical determination of the nodal line patterns (e.g., see Val Baker et al. 2024 and references therein).

While in Chladni figures the distinction is between nodes and antinodes, in Faraday wave phenomena the distinction is between the crests and troughs. The crests of the wavelets are convex and therefore diffuse light away from the eye or camera, thus appearing darker, while the concave troughs of the wavelets act as reflectors, sending light back to the eye or camera thus appearing brighter. The antinodes therefore exist as either darker crests or brighter troughs and the nodal points are the transition, or equilibrium points (e.g., see Figure 2).

vibrating liquid and the corresponding antinodes and
nodes observed in Faraday wave phenomena at a
frequency of 92 Hz.
Since their discovery, Faraday waves have been the subject of successive investigations allowing for advanced applications across a myriad of fields (Liu et al., 2022). These include, for example: measuring the surface tension of soft materials (Shao et al., 2018); developing new photonic devices (Tarasov et al., 2016; Huang et al., 2017); metamaterials (Domino et al., 2016; Francois et al., 2017); applications in cell culture patterns (Takagi et al., 2002; Hong et al., 2020); detecting physiological processes of organisms (Maksymov and Pototsky, 2020); monitoring earthworms (Mitra et al., 2008); modifying soil structure to increase crop yields (Blakemore and Hochkirch, 2017; Lacoste et al., 2018; Ruiz and Or, 2018); medical ultrasound and photoacoustic imaging modalities (Emelianov, 2019); Brillouin light scattering spectroscopy (Ballmann et al., 2019); imaging cancer and healthy cell sounds in water (Reid et al., 2019); laser vibrometry (Blamey et al., 2013); liquid-based templates for directed assembly of microscale materials, including biological entities such as individual cells (Chen et al., 2014); and the development of new methods for eradicating viruses and bacteria (Zinin et al., 2006; Ivanova et al., 2013; Boyd et al., 2020). Furthermore, as noted by Sheldrake and Sheldrake (2017), the behavior of the observed resonant phenomena offers a model for analogous behavior observed in both physical (e.g., Schrödinger, 1926; Xu et al., 1983; Moon et al., 2008; Britton et al., 2012) and biological systems, (e.g., Matsuhashi et al., 1998; Gruntman and Novoplansky, 2004; Gholami et al., 2015; Liu et al., 2017). Indeed, the similarity in the pattern morphology with both biological and physical systems is prevalent and intriguing; examples range from the structure of DNA, the trilobite, the starfish, the phyllotaxis of spiral aloe, and snowflakes, to the hexagon storm of Saturn (Web ref. 1; Gallaire and Brun, 2017). The idea that this pattern morphology is somehow inherent in the universe is not new. Examples range from studies on the structure of the universe (Broadhurst et al., 1990; Battaner and Florida, 1998; Eisenstein and Bennett, 2008) to studies on an archaeal virus (Rice et al., 2004), where notably, in both these cases, sound is the driving force. Understanding Faraday wave phenomena and the details of their formation, morphology and dynamics therefore has huge implications for both scientific knowledge and potential applications.
Previous observations of Faraday waves in a vertically vibrated liquid-filled cuvette concluded that pattern formation is primarily determined by the excitation frequency (Abraham, 1976; Simonelli and Gollub, 1989; Kumar and Tuckerman, 1994; Sheldrake and Sheldrake, 2017), the boundary condition of the observed medium (Kumar and Tuckerman, 1994; Topaz et al., 2004; Francois et al., 2014; Sheldrake and Sheldrake, 2017) and the viscosity of the liquid (Escaler and De La Torre, 2018; Tuan et al., 2015; Westra et al., 2003), (Web ref. 1). The resulting pattern morphology, observed at different frequencies, is shown to form a myriad of patterns, e.g., lattices and stripes, hexagonal shapes, and cubes, etc. (Francois et al., 2014; Topaz et al., 2004). Furthermore, when excitation is by more than one frequency, akin to sounds in nature, there is an indication of pattern selection, where either one frequency pattern is more dominant or a new pattern is formed (Batson et al., 2015).
However, although pattern morphology is somewhat deterministic, it is not entirely predictable (Benjamin and Ursell, 1954; Kumar and Tuckerman, 1994; Bechhoefer et al., 1995). And while the seemingly unpredictable patterns observed at higher frequencies can likely be attributed to mode mixing, and the variations observed at the lower frequencies may be understood in terms of chaos theory and/or catastrophe theory, it is evident that further investigation is required (Thom, 1975; Douady and Fauve, 1988; Sheldrake and Sheldrake, 2017).
Improving our understanding of Faraday wave phenomena is therefore of great importance. Specifically, with the aim of identifying the set of initial conditions that will lead to pattern formation of a specific morphology and dynamic, i.e., the inverse problem. To investigate this further and identify the inverse problem, we first need to define a unified classification language and methodology.
In this paper we report on the observations of Faraday wave phenomena following the excitation by one or more frequencies of a water-filled cuvette observed with the CymaScope Pro Mid Frequency (MF+) instrument. This study explores the varied morphology of patterns observed for different frequencies. In each case the pattern morphology is defined, and where required, a new terminology is suggested. In this paper, we report on both new phenomena and those which were previously reported.
Methodology
The observations were conducted on the CymaScope Pro MF+ instrument that makes wave patterns visible by transforming the periodicity of audio input signals to water wavelet periodicity via acoustic excitation of a fused-quartz water-filled cuvette. The CymaScope Pro MF+ comprises a fuse-quartz cuvette direct-coupled to a Voice Coil Motor (VCM) with a vertically driven piston. The fused-quartz cuvette is acoustically excited by the VCM which, in this configuration, receives audio input signals generated by a Max/MSP patch used within Ableton Live (v11). Electronic filtering is applied in the signal path through a DBX 231 graphic equalizer to ensure that resonances inherent in the VCM-cuvette assembly are negated, resulting in a characteristic flat amplitude response curve. The resulting resonant phenomena for all audio samples were observed via direct ocular viewing while simultaneously recorded using a Canon EOS 5D Mk IV camera (see Figures 3a-c).



The experiments presented here used fused quartz cuvettes of the following shape and size: a triangle of 50 mm length; a rectangle of 24.25 mm length and 15 mm width; a circle of 24.25 mm diameter; and a square of 25.2 mm length. Each cuvette was filled with 2.1 ml; 1.1 ml; 2.3 ml; and 3.3 ml of distilled water, respectively (see Figures 4a-d).
Observations were made for each frequency in accordance with the ~ 250 Hz compatibility range of the CymaScope instrument (see Appendix 1 – 4). In each case, measurements of the wave pattern formations were made, and the following environmental factors were monitored: humidity; air pressure; water temperature; and room temperature. To account for any effects due to prior oscillation (Montagnier et al., 2009), new water was used for each measurement and the order in which the different frequencies were observed was varied and random. The amplitude was manually increased using the decibel level controls and was recorded as audio input in decibel relative to full scale (dBFS). As a unit of measurement, dBFS measures amplitude levels in digital systems, where 0 dBFS is assigned to the maximum level. Throughout this paper all values in dB indicate the input amplitude in dBFS. It should be noted that the videos, and corresponding images, are only able to capture a moment in time, and thus, depending on the compatibility between the Faraday waves and the equipment, are not always able to fully represent the wave dynamics.
Before taking measurements, the following apparatus checks were made: the CymaScope instrument level check with the in-built level meter; the camera mounting level via the camera in live view with the diode reflection as a guide; dust removal from the CymaScope instrument including the cuvette and the camera lens as well as ensuring no dust remained once the cuvette had been filled with water.
Results and Discussion
The resulting observations are described in terms of: (1) pattern formation; (2) pattern morphology; and (3) pattern dynamics. In each case, existing or newly defined terminology is shown in bold.
Pattern Formation
Faraday wave phenomena formed in a vertically vibrated liquid-filled cuvette are dependent on the frequency and the boundary conditions. Moreover, for stable pattern formation to occur: (i) The excitation frequency needs to be compatible with the boundary conditions; (ii) the amplitude needs to be of sufficient magnitude; and (iii) the required time for pattern formation to occur needs to be realized.
In this work we investigated frequencies in the range 8 – 256 Hz (see Appendix 1 – 4). The frequencies where pattern formation occurred are referred to as reference frequencies. Note, this does not mean that the same pattern formation will always occur at these reference frequencies. Sometimes the same pattern will instead be observed at a frequency of +/- a few Hz of the reference frequency and sometimes a completely new pattern will be observed. As well, for the frequencies that showed no stable pattern formation, it does not mean that pattern formation will never occur and/or that the same pattern will occur. The results presented here and in previous work (Sheldrake and Sheldrake, 2017) showed that other factors, internal and external, could affect the occurrence of pattern formation. For example, we see that pattern formation occurs at a specific amplitude, which we refer to as the reference amplitude. This reference amplitude was observed to increase with frequency, for every shape of cuvette investigated, e.g., see Figure 5. However, it should be noted that the pattern morphology, observed under the same conditions, for a given reference amplitude and reference frequency, exhibits variations. For example, at a frequency of 128.4 Hz and an amplitude of -19.4 dB, we see pattern formation occurs in one trial (see Figure 6a) and does not occur in another trial, apparently under the same conditions (see Figure 6b), but clearly an unidentified parameter change was responsible for this anomaly. We refer to this as an unidentified anomaly and to identify the cause of such variations further investigations are required where both internal and external factors are monitored in more detail.

We also see that for pattern formation to occur, a specific time needs to be realized. This time is referred to as the Time to Full Expression (TFE) and was originally defined in 2017 by Sheldrake and Sheldrake. Pattern formation can be observed either by the naked eye, or more quantitatively determined from the grey pixel count (see Sheldrake and Sheldrake, 2017, for a detailed explanation of the grey pixel count methodology). Sheldrake and Sheldrake also note that the TFE varies with a dependence on the frequency and amplitude, e.g., see Figure 7, which shows the TFE decreases with both amplitude and frequency.

as a function of frequency for three frequencies: 56; 111; and 180 Hz.
The top plot shows how the TFE decreases with increasing amplitude and frequency, and the bottom plot shows the variation between trials.
(Figure 5, Sheldrake and Sheldrake, 2017).
Pattern Morphology
Faraday wave phenomena exhibit various pattern morphologies, depending on the driving frequency and boundary conditions of the excited medium. This morphology will vary with time due to changes in both external and internal parameters. We can therefore view the pattern morphology as individual moments in time, which are referred to as cymaglyphs. The term cymaglyph was coined by CymaScope inventor John Stuart Reid and is defined as a cymatic representation of a moment of sound. For example, see Figure 8.

Cymaglyphs can show different degrees of expression of wave pattern formation that we identify and outline below:
Fully expressed cymaglyph – shows a wave pattern morphology with visible crests and troughs. The observed patterns can exhibit either a symmetric wave pattern formation, referred to as a coherent cymaglyph (e.g., see Figure 9), or a non-symmetric, irregular wave pattern formation, referred to as an incoherent cymaglyph (e.g., see Figure 10). In some cases, the wave pattern morphology forms a complex and dynamic structure, which shows a winding and twisting pattern (e.g., see Figure 11).



Chaotically expressed cymaglyph – shows a chaotic wave pattern morphology and occurs when a critically high excitation amplitude is applied. Eventually, as the amplitude is further increased, this chaotic behavior will result in water droplets jumping out of the cuvette (e.g., see Figure 12).

Pre-cymaglyph – shows no visible wave movement to the naked eye, even when the water is optimally excited or overexcited (e.g., see Figure 13). However, in such cases the presence of wave formation within the system cannot be ruled out. Further research investigating the system, from different angles, rather than the surface, may reveal the system morphology in greater detail.

Underexpressed cymaglyph – shows either no wave pattern formation or an incomplete pattern morphology, i.e., a lower degree of expression of the pattern than previously observed. For example, see Figures 14 and 15.


Various aspects of the pattern morphology can also be described by means of the radial symmetry. For example, we can see radial symmetry of 180 degrees, defined as 2-fold or bilateral symmetry; 120 degrees, defined as 3-fold or trilateral symmetry; 90 degrees, defined as 4-fold or quadrilateral symmetry; etc. (e.g., see Figure 16). Sheldrake and Sheldrake (2017) refer to the radial symmetry as n-fold symmetry. We adopt this terminology throughout this paper and any subsequent papers on this topic of research. This radial symmetry pattern morphology is also observed in the triangular, square and rectangular cuvettes, e.g., see Figure 17.


Another aspect of the pattern morphology is a ring pattern, which is usually observed in the center of the cuvette (e.g., see Figure 9c and 9d). For the circular cuvettes we also observe concentric rings. These concentric rings are typically observed in pre-cymaglyphs (e.g., see Figure 13) before sufficient amplitude has been reached to form a cymaglyph and are due to reflections of the CymaScope’s light ring, which consists of 32 diodes. For example, see Figure 18, which shows: the reflection of the light ring before excitation; concentric rings after initial excitation, lower than the critical amount; and the 4-fold pattern formation after sufficient excitation (critical) has been reached. However, in some cases we see a cymaglyph with concentric rings, which can also be referred to as 1-fold symmetry (e.g., see Figure 19). For the square and rectangular cuvettes, we also see patterns of parallel and perpendicular lines; e.g., see Figure 20. We refer to this as a grid cymaglyph.



The variety of observed pattern morphologies is dependent on the frequency, amplitude and boundary conditions. However, due to the complex and delicate nature of Faraday waves we cannot confidently attribute pattern morphologies to specific parameters. For example, for specific frequencies measured under the same conditions, including amplitude conditions, a variety of cymaglyphs are produced (i.e., different symmetry fold, pattern morphology, etc.). For example, see Figure 21, which for a circular cuvette at a frequency of 40 Hz, shows a 14-fold symmetry pattern formation in one trial and a 4-fold symmetry pattern in another trial. We refer to this as an unidentified anomaly and, to identify the cause of such variations in the manifestation of cymaglyphs, further investigations are required where both internal and external factors are monitored in more detail.

Pattern Dynamics
When investigating how the pattern morphology changes with time, we observed both stable and unstable pattern morphology. In the cases with stable pattern morphology, a coherent and steady wave formation is observed, resulting in a series of stable cymaglyphs (e.g., see Figure 22). However, unstable pattern morphology shows a dynamic and changing wave formation, resulting in a series of unstable cymaglyphs (e.g., see Figures 23-28). Note, these unstable cymaglyphs can be coherent and incoherent.

In some cases, this unstable pattern morphology could be due to the fact that no standing wave formation is compatible within the boundary conditions at that specific frequency. For example, see Figure 23, which shows an incoherent morphing pattern in a series of unstable cymaglyphs.

In other cases, it could be that the specific frequency is between two stable pattern formation frequencies and thus oscillations between alternate patterns occur. This is referred to by Sheldrake and Sheldrake as a transition between stability points and addressed by Thom as a “competition of resonance” (Sheldrake and Sheldrake, 2017; Thom, 1975). In the latter case, where the patterns appear to oscillate between two or more alternating symmetry folds, the resulting series of unstable cymaglyphs may appear as a combination of two or more symmetry folds, e.g., see Figures 24 and 25. Sheldrake and Sheldrake observed such a phenomenon and referred to it as a composite of two alternating phases of oscillation, whereby crests become troughs and troughs become crests (Sheldrake and Sheldrake, 2017).


However, these oscillations could also be alternating between the same pattern but of different degrees of expression (e.g., see Figure 26) and/or morphing between different patterns (e.g., see Figure 27). Note these oscillations between alternating patterns are not necessarily regular in time and can instead exhibit an irregular dynamic.


We also see coherent dynamics such as rotation where clear circular motions are visible, either clockwise, counterclockwise, or both, and for varying degrees of rotation, i.e., 10 degrees, 45 degrees, 360+ degrees (see Figures 28-30). We conjecture that the observed rotational effects are produced by interference between the sinusoidal driving frequency and inherent resonances in the CymaScope instrument’s fused quartz cuvette. The quartz resonances cause reflections of acoustic energy into the water that are out of phase with the driving frequency, causing a type of beat frequency effect that manifests as hunting rotational effects in the cymaglyph. Although the natural resonances in the characteristic frequency response curve of the CymaScope are identified by a calibrated accelerometer system and tuned out by the graphic equaliser, this form of frequency response correction cannot take account of dynamic resonances inherent within the fused quartz of the cuvette. Interference between the opposing waveforms, i.e., those emanating from the centre outwards and those reflected from the outer boundary and directed inwards, will also contribute to the observed rotational effects. Further investigation of these rotations and the system dynamics in general, could yield insights on similar systems seen in nature, such as whirlpools.



In other cases, we observed the pattern morphology to change in such a way that it appeared as a toroidal dynamic, i.e., a spherical inwards and outwards dynamic, which is visible on the water surface, either observed in parts or in the entire wave formation (e.g., see Figure 31). It would be interesting to investigate this further and verify if such a toroidal movement is indeed occurring. This would require a new methodology that allows for 3D viewing, as well as a suspension of particles that could be used to track the motion. Again, like rotational dynamics, improving our understanding of the observed toroidal dynamics could yield insights on similar systems seen in nature.

So far, we have described the resulting Faraday wave phenomena when a liquid-filled cuvette is excited by a single frequency. However, it is known that distinct phenomena are observed when two or more frequencies are applied mechanically; for example, Topaz et al., (2004), Francois et al. (2014), and Batson et al. (2015) report such distinct behaviors. Furthermore, natural sounds are made up of multiple frequencies, so studying cymaglyphs created by two or more frequencies is of great interest and a first step toward understanding more complex sounds akin to nature. As well, such investigations can reveal how interference leads to pulsating and/or rotational effects, which has numerous applications, from cell manipulation to water containment and filtration methods. To investigate this further, we simultaneously applied dyad frequencies of different ratios. We refer to the resulting cymaglyph as an interference cymaglyph. In each case we observe a morphing of the respective patterns to various degrees where the wave pattern will settle into either: a summation pattern; a new pattern, which in some cases incorporates elements from patterns associated with the individual frequencies; or one of the individual frequency patterns (e.g., see Figures 32-34). However, it should be noted, the distinction between patterns is not always obvious and highly subjective. For example, Figure 32c and 33c could both be exhibiting either a summation pattern or a new pattern. To investigate this further, more data is needed as well as the utilisation of image analysis software.



Summary
In this paper we identified different phenomena in Faraday wave formation following sound excitation as observed with the CymaScope instrument. We suggest a unified terminology of classification to describe these phenomena, solidifying previously coined and newly defined vocabulary. Notably, an observed pattern can be defined with multiple terms describing every aspect of the pattern formation, morphology and dynamics. For example, Figure 22 can be described as a fully expressed coherent cymaglyph of 10-fold symmetry at a reference frequency of 69 Hz with a stable dynamic showing both clockwise and counterclockwise rotation.
Faraday wave phenomena are intrinsically delicate and susceptible to variations, yet we also observe consistency in pattern formation determined by frequency and boundary conditions. However, variation in pattern formation and morphology are observed, which we refer to as an unidentified anomaly. Further investigations with the aim of identifying the inverse problem are therefore needed.
With the help of a unified terminology and methodology, future studies can be improved and yield further insights on Faraday wave phenomena.
References
Abraham R (1976). Vibrations and the Realization of Form. In: Evolution and Consciousness: Human Systems in Transition (ed. E Jantsch, CH Waddington), Reading, MA, USA: Addison-Wesley, pp. 134–148.
Andrade ENDC (1950). Wilkins Lecture: Robert Hooke. Proceedings of the Royal Society of London. Series B, Biological Sciences 137(887): 153–187.
Ballmann CW, Meng Z, Yakovlev VV (2019). Nonlinear Brillouin spectroscopy: What makes it a better tool for biological viscoelastic measurements. Biomed. Opt. Express 10(4): 1750–1759.
Batson W, Zoueshtiagh F, Narayanan R (2015). Two-frequency excitation of single-mode Faraday waves. J. Fluid Mech. 764: 538–571
Battaner E and Florido E (1998). The egg-carton Universe. arXiv: Astrophysics/9802009v1.
Bechhoefer J, Ego V, Manneville S, Johnson B (1995). An experimental study of the onset of parametrically pumped surface waves in viscous fluids. J. Fluid Mech. 288: 325–350.
Bechhoefer J and Johnson B (1996). A simple model for Faraday waves. Am. J. Phys. 64:1482–1487.
Benjamin TB and Ursell F (1954). The stability of the plane free surface of a liquid in vertical periodic motion. Proc. Roy. Soc. A 225(1163): 505–515.
Blakemore R and Hochkirch A (2017). Restore earthworms to rebuild topsoil. Nature 545:30.
Blamey J, Yeo LY, Friend JR (2013). Microscale Capillary Wave Turbulence Excited by High Frequency Vibration. Langmuir 29(11): 3835–3845.
Boyd B, Suslov SA, Becker S, Greentree AD, Maksymov IS (2020). Beamed UV sonoluminescence by aspherical air bubble collapse near liquid-metal microparticles. Sci. Rep 10: 1501.
Britton JW, Sawyer BC, Keith AC, Wang CCJ, Freerick JK, Uys H, Biercuk MJ, Bollinger JJ (2012). Engineered two-dimensional Ising interactions in a trapped-ion quantum simulator with hundreds of spins. Nature 484: 489–492.
Broadhurst TJ, Ellis RS, Koo DC, Szalay AS (1990). Large-scale distribution of galaxies at the Galactic poles. Nature 343: 726-728.
Chen P, Luo Z, Guven S, Tasoglu S, Weng A, Ganesan AV, Demirci U (2014). Microscale assembly directed by liquid-based template. Advanced Materials 26(34): 5936-5934.
Chladni EFF (1787). Entdeckungen über die Theorie des Klanges. Leipzig, Germany: Weibmanns, Erben und Reich.
Domino L, Tarpin M, Patinet S, Eddi A (2016). Faraday wave lattice as an elastic metamaterial. Phys. Rev. E. 93:1–5.
Douady S and Fauve S (1988). Pattern Selection in Faraday Instability. EPL 6(3): 221–226.
Douady S (1990). Experimental study of the Faraday instability. J. Fluid Mech. 221: 383-409.
Eisenstein DJ and Bennett CL (2008). Cosmic sound waves rule. Physics Today 61(4): 44–50.
Emelianov S (2019). Light, sound, nanobubbles: New approach to contrast-enhanced ultrasound and photoacoustic imaging. J. Acoust. Soc. Am. 145(3): 1779.
Escaler X and De La Torre O (2018). Axisymmetric vibrations of a circular Chladni plate in air and fully submerged in water. Journal of Fluids and Structures 82: 432–445.
Faraday M (1831). XVII. On a peculiar class of acoustical figures; and on certain forms assumed by groups of particles upon vibrating elastic surfaces. Philosophical Transactions of the Royal Society of London 121: 299–340.
Fauve S (1998). Pattern forming instabilities. In: Hydrodynamic and nonlinear instabilities (ed. C Godrèche, P Manneville), Cambridge, UK: Cambridge University Press, pp. 387–491.
Francois N, Xia H, Punzmann H, Ramsden S, Shats M (2014). Three-Dimensional Fluid Motion in Faraday Waves: Creation of Vorticity and Generation of Two-Dimensional Turbulence. Phy. Rev. X. 4(2): 021011–021021.
Francois N, Xia H, Punzmann H, Fontana PW, Shats M (2017). Wave-based liquid-interface metamaterials. Nat. Commun 8: 14325.
Galilei G (1638). Dialogues Concerning Two New Sciences, trans. H. Crew and A de Salvio, 1914. London, UK: Macmillan.
Gallaire F and Brun PT (2017). Fluid dynamic instabilities: theory and application to pattern forming in complex media. Phil. Trans. R. Soc. A 375: 20160155.
Gholami A, Steinbock O, Zykov V, Bodenschatz E (2015). Flow-Driven Waves and Phase-Locked Self-Organization in Quasi-One-Dimensional Colonies of Dictyostelium discoideum. Phys. Rev. Lett 114(1): 018103–018105.
Gruntman M and Novoplansky A (2004). Physiologically mediated self/non-self-discrimination in roots. PNAS 101(11): 3863–3867.
Hong S-H, Gorce JB, Punzmann H, Francois N, Xia H (2020). Surface waves control bacterial attachment and formation of biofilms in thin layers. Sci. Adv 6(22): 1654-1685.
Huang S-W, Vinod AK, Yang J, Yu M, Wong CW (2017). Multispectral optical frequency comb based on microresonator Faraday instability. In IEEE Photonics Conference (IPC), Orlando, FL, USA.
Ivanova EP, Hasan J, Webb HK, Gervinskas G, Crawfoed RJ (2013). Bactericidal activity of black silicon. Nat. Commun 4: 2838: 1–7.
Kumar K and Tuckerman LS (1994). Parametric instability of the interface between two fluids. J. Fluid Mech. 279: 49-68.
Lacoste M, Ruiz SA, Or D (2018). Listening to earthworms burrowing and roots growing-acoustic signatures of soil biological activity. Sci. Rep 8:10236.
Liu J, Martinez-Corral R, Prindle A, Lee D-Y, Larkin J, Gabalda-Sagarra M, Garcia-Ojalvo J, Süel GM (2017). Coupling between distant biofilms and emergence of nutrient time-sharing. Science 356(6338): 638-642.
Liu J, Song W, Ma G, Li K (2022). Faraday Instability in Viscous Fluids Covered with Elastic Polymer Films. Polymers 14: 2334.
Maksymov IS and Pototsky A (2020). Excitation of Faraday-like body waves in vibrated living earthworms. Sci. Rep 10: 8564.
Matsuhashi M, Pankrushina AN, Takeuchi S (1998). Production of sound waves by bacterial cells and the response of bacterial cells to sound. J. Gen. Appl. Microbiol 44(1): 49-55.
Mitra O, Callaham MA, Smith ML, Yack JE (2008). Grunting for worms: Seismic vibrations cause Diplocardia earth-worms to emerge from the soil. Biol. Lett 5(1): 16–19.
Montagnier L, Aïssa J, Ferris S, Montagnier J-L, Lavallée C (2009). Electromagnetic signals are produced by aqueous nanostructures derived from bacterial DNA sequences. Interdiscip Sci 2: 81-90.
Moon CR, Mattos LS, Foster BK, Zeltzer G, Ko W, Manoharan HC (2008). Quantum Phase Extraction in Isospectral Electronic Nanostructures. Science 319(5864): 782–787.
Perlin M, Schultz WW (2000). Capillary effects on surface waves. Annu. Rev. Fluid Mech. 32(1), 241–274.
Reid JS, Park BJ, Ji S (2019). Imaging Cancer and Healthy Cell Sounds in Water by Cymascope, Followed by Quantitative Analysis by Planck-Shannon Classifier. WATER Journal 11.
Rice G, Tang L, Stedman K, Roberto F, Spuhler J, Gillitzer E, Johnson JE, Douglas T, Young M (2004). The structure of a thermophilic archaeal virus shows a double-stranded DNA viral capsid type that spans all domains of life. Proc. Natl. Acad. Sci. U.S.A., 101(20): 7716-7720.
Ruiz SA and Or D (2018). Biomechanical limits to soil penetration by earthworms: Direct measurements of hydroskeletal pressures and peristaltic motions. J. R. Soc. Interface 15(144).
Schrödinger E (1926). An undulatory theory of the mechanics of atoms and molecules. Phys. Rev. 28(6): 1049-1070.
Shao X, Saylor JR, Bostwick JB (2018). Extracting the surface tension of soft gels from elastocapillary wave behavior. Soft Matter 14: 7347–7353.
Sheldrake M and Sheldrake R (2017). Determinants of Faraday Wave-Patterns in Water Samples Oscillated Vertically at a Range of Frequencies from 50-200 Hz. WATER Journal 9: 1–27.
Simonelli F and Gollub JP (1989). Surface wave mode interactions: effects of symmetry and degeneracy. J. Fluid Mech. 199: 471-494.
Takagi S, Krinsky V, Pumir A, Frelin C (2002). The use of Faraday instability to produce defined topological organization in cultures of mammalian cells. Int. J. Bifurc. Chaos 12(9): 2009-2019.
Tarasov N, Perego AM, Churkin DV, Staliunas K, Turitsyn SK (2016). Mode-locking via dissipative Faraday instability in a fiber laser. Nat. Commun 7: 12441.
Thom R (1975). Structural Stability and Morphogenesis, trans. DH Fowler with a Foreword by CH Waddington. Reading, MA, USA: W. A. Benjamin, INC.
Topaz CM, Porter J, Silber M (2004). Multifrequency control of Faraday wave patterns. Phys. Rev. E. 70(6): 066206–066216.
Tuan PH, Wen CP, Chiang PY, Yu YT, Liang HC, Huang KF, Chen YF (2015). Exploring the resonant vibration of thin plates: Reconstruction of Chladni patterns and determination of resonant wave numbers. J. Acoust. Soc. Am. 137(4): 113–2123.
Val Baker A, Csanad M, Fellas N, Atassi N, Mgvdliashvili I, Oomen P (2024). Exploration of Resonant Modes for Circular and Polygonal Chladni Plates. Entropy 26 (3):264.
Waller MD (1938). Vibrations of free circular plates. Part 3: A study of Chladni’s original figures. Proceedings of the Physical Society 50(1): 83–86.
Waller MD (1961). Chladni Figures: A Study in Symmetry. London, UK: G. Bell & Sons.
Westra M-T, Binks DJ, Van De Water W (2003). Patterns of Faraday waves. J. Fluid Mech. 496: 1-32.
Xu Y, Vest CM, Murray JD (1983). Holographic interferometry used to demonstrate a theory of pattern formation in animal coats. App. Opt 22: 3479-3483.
Zinin PV, Allen JS III, Levin VM (2006). Mechanical resonances of bacteria cells. Phys. Rev. E. 72(6): 061907.
Web References
Web reference 1: Reid JS (2023). Retrieved from https://cymascope.com/ on [14/17/2023].
Discussion with Reviewers
Reviewer 1: In your manuscript, you describe a “toroidal wave,” but it is not clear to me from looking at the figures and videos that there is a toroidal dynamic or that there is any such thing as a toroidal “‘wave.” This throw-away term, toroidal wave, needs unpacking in much more detail if you want to describe the wave patterns in the CymaScope in terms of toroids.
Please also explain how this might relate to the diagrams in Figures 8a and b and discuss whether these patterns of flow are in fact happening in the CymaScope, giving more evidence than an uninterpreted video.
Authors: The images, and their corresponding descriptions, have now been amended to clarify the observations presented. As well, the suggestion to observe the videos at a slow speed has been included. However, as you say, without further investigations we cannot be sure that what we are observing is a toroidal dynamic. We have therefore amended this to read as, ‘… appeared as a toroidal dynamic …’ and included a sentence that suggests a potential experimental verification. Furthermore, this definition has been moved to the new section, Pattern Dynamics.
Reviewer 1: When describing pattern morphology and a ring pattern, you say that the concentric rings may not necessarily orbit in the same direction or at the same speed, but do not give examples to support this assertion. You also say they may be a result of concentric ripples in the cuvette. Surely, they must be a result of concentric ripples or standing waves.
Authors: This definition has now been amended for clarification and moved to the new section, Pattern Morphology. In this section we describe ring patterns and make the distinction between ring patterns observed in the center of the cuvette; “concentric rings” observed due to reflections of the CymaScope light ring; and concentric rings observed as a 1-fold symmetry pattern.
Reviewer 1: You describe the effects of playing two frequencies at the same time, compared with the patterns produced by the single frequencies, or what I refer to as the “interference cymaglyph.” This would be a good subject for systematic investigation rather than giving random examples. For example, you refer to patterns from video 18 (see Appendix). The first shows a 3/6-fold pattern and the second shows a 4-fold pattern. When combined they give a 3/6-fold pattern. Here, one of the two patterns seems to predominate over the other. By contrast, video 19 combines two patterns, the first of which shows a 2/4-fold pattern and the second a 1/2-fold pattern that, when combined, produce a ring pattern.
Here is an opportunity for you to provide a systematic investigation of the interaction of different frequencies, rather than just a few random examples. It would be good to insert a systematic study of these frequency combinations illustrating in how many cases one pattern predominates over the other, in how many cases both are combined in unstable mixtures, and in which cases the two patterns produce something that is different from both of these, like the rings in video 19.
Authors: This definition has now been amended for clarity and moved to the new section, Pattern Dynamics. It is noted, as well, that these are just preliminary observations and that an in-depth study on this topic should be done in the future. However, we decided to keep the name interference cymaglyph for both harmonic and inharmonic frequencies to distinguish between cymaglyphs that are from single frequencies and those that are from multiple frequencies.
Reviewer 1: In your paper, “Amplitude anomaly” hardly required an entry in a taxonomy. They would simply observe that sometimes it requires a higher amplitude to elicit a pattern than others. It is also not clear why decibel readings have been given here as negatives, e.g. – 23.3dB and -17.8dB.
In your paper, “cymaglyph anomaly” hardly merits a separate taxonomic entry. You could simply mention that sometimes under those same conditions different patterns formed. However, the conditions might not be the same because the atmospheric pressure varies from day to day and there are other variables that affect the expression of the pattern. In Sheldrake and Sheldrake (2017), Figure 11 shows 3 replicates at every frequency between 50 and 199 Hz, observed on three different days. On some days the patterns appeared at slightly different frequencies for unknown reasons. But one uncontrolled variable was atmospheric pressure, which could have had some effect and accounted for these “cymaglyph anomalies.” However, the discussion that follows sheds no light on the reason for these anomalies and tells us nothing more than can already be observed from Sheldrake and Sheldrake (2017), Figure 11. This kind of variation is one of the main features discussed that simply makes the point that several different factors could affect these discrepancies, without illuminating what is going on any further.
Authors: The term amplitude anomaly describes the observed variation in pattern formation for a specific amplitude, despite being under the same conditions. For example, for a given amplitude, pattern formation occurs in one trial and not in another. This has now been moved to the new section, Pattern Formation, and is referred to as an unidentified anomaly. Regarding the decibel readings being negative, this is due to the signal level in the audio system being lower than the 0 dB full scale. A note to clarify this has now been included in the Methodology section. The term cymaglyph anomaly describes the observed variation in cymaglyphs (i.e., different symmetry fold, pattern morphology, etc.) for specific frequencies measured under the same conditions, including amplitude. For example, for a given frequency a multitude of pattern morphologies can be observed, even for the same amplitude. This has now been moved to the new section, Pattern Morphology, and is referred to as an unidentified anomaly.
Reviewer 1: You do not refer in this paper to a key paper on Faraday waves that gives a far more systematic understanding of the patterns formed and the radial and circumferential modes, together with a mathematical analysis in terms of Bessel functions, namely Shao et al., (2021). Surface pattern wave formation in a cylindrical container. Journal of Fluid Mechanics, Vol. 915, A19. Reading this paper, you will be able to make more sense of your observations and produce a description of them that is less confused and confusing.
Authors: After reading the suggested paper, and others, 12 new references have been added: Abraham, 1976; Bechloefer, 1995; 1996; Benjamin and Ursell, 1954; Douady, 1988; Fauve, 1998; Gallaire and Brun, 2017; Kumar and Tuckerman, 1994; Perlin and Schultz, 2000; Simonelli and Gollub, 1989; Val Baker, 2024, and Montagnier, 2009. We’ve also referred to multiple studies and applications in our paper and referenced them accordingly.
Reviewer 2: The rotation “Rotation” is interesting. Is there any physical explanation behind or beyond simple chance?
Authors: The rotation is interesting and something that we hope to investigate further. We don’t think the rotational effects are simply chance. Instead, it is likely the result of the resonant characteristics of the mechanical system as well as interference effects between the driving wave and the reflected wave. These rotational effects, along with a note on the possible causes has now been moved to the new section, Pattern Dynamics.
Reviewer 2: The key question for me is whether these have any relevance to pattern formation in biology or are just interesting patterns. The motivation of this taxonomy paper is to try and unify the terms that are used, which is well worthwhile as it can be confusing if different papers use the different terminologies to refer to the same thing, and the taxonomy discussion appears to be a welcome contribution to the field. I think that the paper might be improved if the Taxonomy section was more clearly identified as the results section as the results of your experiments are clearly described in this section. The taxonomic terms only refer to where such effects are seen, but most of them have already been reported by other authors, so the results section should also note, with references, specific examples where the effect being described has been seen by others. It may then be possible to more clearly tie together the discussion and the results, perhaps by using bold font whenever a specific taxonomic term was used so that it is easier to see the terms being used and easier to go back to the results section and identify the more detailed results that are being discussed.
Authors: The Taxonomy section has now been updated to the Results and Discussion section. The terminology is now clearly shown in bold and in cases where the definition has previously been used, we refer to the relevant reference. As well we have now noted the connection between biology and the patterns observed.
Reviewer 3: What is the rationale for experimenting with triangular, square, and rectangular cuvettes? For instance, is it to reveal forms not visible or possible in circular configurations, to study them as biological analogues, or to go beyond what math or simulation can show us? Do we gain any understanding, or can we go on simply trying things and describing them without providing insights or understanding? I would appreciate a systematic study of a triangular boundary as a function of frequency (like what Sheldrake and Sheldrake, 2017, did for circular cuvettes) or a systematic study of rectangular shapes as the sides are progressively changed from a square into a straighter rectangle.
Authors: The rationale for experimenting with different shaped cuvettes is to investigate the effects of boundary conditions on resonant wave phenomena. The work presented in this paper is qualitative and shows the resulting resonant phenomena for a broad range of frequencies for 4 different shaped cuvettes. We hope that these observations will be the basis to build upon when future work allows us to investigate further, gaining understanding of their behavior and ultimately describing them mathematically.
Reviewer 3: What is the rationale for using two frequencies? Again, a systematic study of frequencies, rather than “being all over the place,” would be valuable.
Authors: The rationale for using two frequencies is to observe the resulting effects on the resonant wave phenomena. How do the frequencies interact with each other? What other effects does it lead to? Again, this paper is qualitative and looks at a broad range of parameters that can affect the observed resonant phenomena. We identified 3 different effects observed when using two frequencies instead of one. We hope that this can be investigated further in future and more insights can be gained. Understanding resonant wave phenomena, especially with water, has huge implications ranging from water containment, manipulation, and understanding how biological forms take hold.
Reviewer 4: What I found most interesting in this paper was the different Taxonomies and especially “Over Excited Water” and “Toroidal Rotation.” From my point of view the method is rather flat, i.e., 2-dimensional and the Taxonomy of Over Excited Water gives a thought that there is something happening in 3-dimensions as well. The Taxonomy of Toroidal Rotation gives the idea of using sound as a means for propulsion. If we can get the correct combination of frequencies, we can get the water or maybe air to shape a rotating torus to have a new type of propulsion for airplanes.
Authors: The cymaglyphs show the surface morphology, so indeed, it would be interesting to observe the pattern morphology and dynamics from all angles and within the container so that the full 3D aspect of these phenomena can be observed. To achieve this, a new methodology that allows for 3D viewing would be required. As well. the suspension of particles within the liquid could be utilized so that the motion could be tracked. This is a great idea for a future study and has now been mentioned in the new section, Pattern Dynamics. As you mention, a future study that can help us better understand toroidal dynamics could not only yield insights on similar systems seen in nature, but could also help us design systems such as propulsion systems.
Reviewer 5: A key question for me is whether these have any relevance to pattern formation in biology or are just interesting patterns; a question that has been raised at the water conference.
Authors: Yes, the pattern morphologies we observe with the CymaScope are relevant to biology and one of the driving factors in studying them. This has been the case since Faraday waves have been studied and has yielded myriad applications across all fields of science, including biology. We have now included this in the introduction, where an extensive list of potential applications is given. Throughout the paper, when relevant, the similarity to biological or physical systems and dynamics is mentioned.
Appendix 1: Resonant wave phenomena in a 50 mm triangular cuvette for frequencies 8 – 256 Hz.













Appendix 2: Resonant wave phenomena in a 24.25 x 15 mm rectangular cuvette for frequencies 16 – 256 Hz.










Appendix 3: Resonant wave phenomena in a 25.2 mm square cuvette for frequencies 16 – 256 Hz.










Appendix 4: Resonant wave phenomena in a 24.25 mm circular cuvette for frequencies 16 – 256 Hz.










