Maxwell’s equations define and shape our lives, describing the physics from the electromagnetic radiation from the Sun to the interactions of atoms . They are also essential for understanding and controlling light, thereby leading to the advanced photonic devices that are very critical to our lives and global economy: optical fibres, resonators, switches, diffraction gratings, displays and many others , , .
Maxwell’s equations are intrinsically linear, and in most cases we may assume that all optical interactions are also linear. We can say that when red light passes through an optical medium, for example a pair of glasses, it will still be red when it exits the medium. The linearity of optical processes allows us to describe the physics of light separately from the physics of the medium. The medium modifies light, but light does not change the medium.
The possibility of producing nonlinear optical effects with intense light was raised in the 1920s–1940s (see , ). However, nonlinear optics was established only with the invention of the laser , , , . When the intensity of light is sufficiently strong, light can change the properties of the medium it is travelling through, which in turn alters the optical field. In this regime, it is no longer possible to treat the light and medium separately. These mutual light-matter interactions lead to new optical effects  and new nonlinear photonic devices , , , .
One of the most important nonlinear processes is second harmonic generation , , , where high-intensity pump light enters a nonlinear optical material and a weak optical beam is generated at a frequency twice that of the original pump beam. Other important nonlinear processes are third harmonic generation, sum- and difference-frequency generation, and intensity-dependent index of refraction (Kerr effect) to name a few .
High-power lasers required to produce nonlinear effects are, however, expensive and energy-consuming, and the high intensity of light produced by them can cause damage to live cells and tissue, raising significant health and safety issues , . This has motivated the research community to seek alternative solutions for the generation of optical nonlinearities with low-power light.
There have also been mostly sporadic attempts to achieve nonlinear optical phenomena with low-power, incoherent light or at least to reproduce functionality of some nonlinear photonic devices without lasers , , , , , , , , , , . However, even though the number of such works is gradually increasing and they attract keen interest of researchers in the field of nonlinear optics, they have not had a significant impact yet. This is because the fulfilment of the conditions of phase matching and high optical intensity with incoherent light is even more challenging than with laser light.
This does not imply though that phase matching cannot be satisfied with incoherent light. For example, pioneering experiments demonstrated sum-frequency generation via mixing of a ruby laser signal with an incoherent signal produced by a mercury lamp . In that paper, exactly the same technique as in contemporary “all-laser” experiments (special orientation of a KDP crystal ) was employed to satisfy phase matching. However, experimental schemes using non-laser light were mostly abandoned in the late 1960s because lasers operating at different wavelengths became available.
The requirement for high-intensity light is also difficult to satisfy with incoherent light. Consider a 1 W pulsed laser and a 1 W light-emitting diode (or a gas discharge lamp). One can generate optical nonlinearities with the laser because light emitted by it is coherent and can therefore be focused so that the medium experiences very high instantaneous intensities. Even though the light-emitting diode (gas discharge lamp) has the same time-averaged power, the maximum intensity of light focused using the same optical components is fundamentally limited  by laws of thermodynamics and the concept of optical étendue defining how “spread out” the light is in area and angle .
In this paper, we introduce a conceptually different approach – photonic devices capable of generating new optical frequencies with low-power light and nonlinear processes of non-optical origin. Here we suggest to use low-power, incoherent light to harness giant nonlinear acoustic properties of gas bubbles and liquid droplets, and exploit them in the optical domain. Effectively, this will reproduce the behaviour that is usually achieved via a nonlinear optical interaction enabled by high-power laser light.
Sound is described by travelling pressure waves that result from the motion of particles in the medium through which the sound wave is propagating. Although there are fundamental differences between optics and acoustics, it is possible to treat sound propagation using methods analogous to those used for light. In particular, concepts such as wavelength, speed of propagation and intensity all have direct analogues between optics and acoustics. But while light fields typically do not produce large changes in the optical medium, with sound, atoms are being moved relatively large distances even at low intensities, creating non-uniform pressure fields and hence non-uniform speeds of sound in the medium. This in turn means that acoustic nonlinearities are far easier to observe than optical nonlinearities at comparable power levels , .
As a practical example, we first consider the optical frequency comb (Figure 1A). An optical frequency comb is a spectrum consisting of a series of discrete, equally spaced elements that have a well-defined phase relationship between each other , , , . Frequency combs have found important practical applications in precision measurements, microwave generation, telecommunications, astronomy, spectroscopy and imaging , , , . Combs have typically been produced by mode-locked lasers or exploiting nonlinear optical effects in optical fibres , , ,  and nonlinear photonic microresonators , , ,  (see Figure 1A and Section 2).
Acoustic nonlinearities provide an alternate low-power route to generating an optical frequency comb (Figure 1B and C), which we discuss more completely in Section 3.4. This effectively reproduces the result of a classical nonlinear optical interaction (Figure 1A). High-power light is not required for the conversion of acoustic nonlinearities into optical signals. This implies that our approach would, in general, allow the use of low-cost, low-power incoherent optical sources such as light-emitting diodes and discharge lamps for tasks such as optical frequency comb metrology.
It is noteworthy that the terminology used to describe the appearance of new optical frequencies due to the interaction of light with nonlinear acoustic waves is not yet fully stabilised. In the following, we will continue using the term “generation” as this is often employed in nonlinear optics as well as directly relatable to optical frequency combs. This term is also often used to discuss outcomes of Raman and stimulated Brillouin scattering interactions . We note that some researchers use other terms such as “sideband generation” and “higher-order sideband comb” , . Moreover, new physical phenomena were predicted by considering the interaction of light with shock waves , . (We will discuss shock waves below, e.g. in Section 6.) Because suitable terminology did not exist, the authors provided an elegant explanation of those effects by using the terms “shock-like modulation” of light and “transfer of light frequency”.
This remainder paper is organised as follows. In Section 2, we overview the recent advances in the research direction of low-power nonlinear optics. Section 3 introduces the concept of reproduction of nonlinear optical effects with nonlinear sound and low-power light. The discussion in Section 3 is supported by theory and numerical simulations demonstrating the possibility of generating optical frequency combs from nonlinear ultrasound. In Section 4, we discuss the giant acoustic nonlinearities of liquid droplets and gas bubbles in liquids. We also overview experimental techniques enabling the conversion of giant acoustic nonlinearities into the optical domain. Section 5 focuses on the possibility of employing the plasmonic properties of optical nanoantennas and liquid-metal nanoparticles to enhance the interaction between nonlinear sound and light. Finally, in Section 6, we explore novel opportunities enabled by the interaction of light with extreme nonlinear behaviour of acoustically driven gas bubbles.
2 Low-power nonlinear photonics
To achieve nonlinear optical processes with low-power laser light, one can (i) employ nonlinear optical materials with an intrinsically high nonlinear response and/or (ii) create resonant structures to locally increase the light-matter interactions that lead to stronger nonlinear effects. Strong nonlinear optical effects can be achieved in non-resonant photonic structures such as dielectric slot waveguides.
Despite the range of strategies available, in practice, the choice of nonlinear optical architectures is often limited by peculiarities of the available fabrication techniques for the desired optical material, toxicity (important for biomedical applications), cost and other factors. Here, one may employ materials with high optical nonlinearities such as chalcogenide glasses , , silicon ,  and AlGaAs , . Alternatively, one can exploit optical nonlinearities of noble metals such as gold , , . However, factors such as immature fabrication technology for chalcogenide glass, high linear and nonlinear losses for semiconductors and high losses in metal structures due to the excitation of plasmon resonances  motivate the search for alternative materials such as silica-based glasses  and indium tin oxide .
A large and growing body of research investigates high-quality factor (high-Q) nonlinear photonic micro-resonators , , , , , , , , , , , , , , , , , integrated optical waveguides , , , , optical nanoantennas , , , , , , , , , , , , ,  and metasurfaces , , , . A judicious design of such structures enables large local enhancements of the optical field intensity, which in turn leads to strong nonlinear optical effects.
However, there also are significant drawbacks with the use of high-Q micro-resonators. First, high intensity of light trapped inside a high-Q resonator can give rise to considerable multiphoton absorption losses. Indeed such losses are a significant obstacle for silicon-based photonic devices operating at telecommunication and mid-infrared wavelengths , . Secondly, although low-power nonlinear effects are possible mostly with ultra-high-Q microresonators, such devices are sensitive to even minor fabrication imperfections, which compromises their scalability and on-chip integrability. High-Q microresonators excited at their resonance frequency are also sensitive to external perturbations, which significantly complicates their operation in real-life conditions.
Consequently, alternative strategies to achieve optical nonlinearities with low-power light have been proposed. For example, an integrated optical frequency comb device consuming <100 mW of electric power was demonstrated in , where a decrease in power consumption was achieved with an integrated microresonator-laser on-chip structure. It was also suggested that strong nonlinear optical effects suitable for optical comb generation may be achieved with low-power laser light in micro-ring resonators operating in a period-doubling bifurcation regime . In , it was suggested that a single quantum dot integrated with a micropillar resonator  can be used to generate optical frequency combs with ultra-low-power laser light of just of few nW. However, these new ideas still employ regular nonlinear optical processes and materials and require a source of coherent light. Consequently, they also inherit some of the fundamental limitations of modern nonlinear photonic devices such as high optical losses, sensitivity to fabrication imperfection, toxicity of nonlinear optical materials and so on.
3 Photonics with nonlinear sound
It is instructive to start our discussion with a brief overview of the physics that underpins the nonlinear processes in optics and acoustics. We will use similar mathematical approaches to analyse optical and acoustic waves in their respective nonlinear media. This will open up an opportunity to compare the strength of the nonlinear effects.
3.1 Origin of nonlinear optical effects
By considering Maxwell’s equations and taking into account the polarisation term, written for simplicity in scalar form as P=ε0(χ0+χ(1)E+χ(2)E2+χ(3)E3+…), where the χ(j) are the nonlinear susceptibilities, it can be shown that the polarisation of the material responds with powers of the electric field E. If we assume that E=E0exp(iωt), where ω is the frequency of the incident light and t is time, then the polarisation is P(t)=ε0χ0+ε0(χ(1)E0eiωt +χ(2)E02e2iωt +χ(3)E03e3iωt +…). By grouping all the terms in powers of 2 and higher into the nonlinear polarisation PNL, we can arrive to the nonlinear wave equation ∇2E−(n2/c2)∂2E/ ∂t2=−T, where n is the optical refractive index, c is the speed of light and T=−μ0∂2PNL/∂t2 is the source term.
So the system responds as an anharmonic oscillator and the polarisation term contains not just the fundamental frequency ω, but also multiples of this frequency (and in general a dc component). For example, in the original second harmonic generation experiment , due to the anharmonic response of the medium the high-intensity red light produced by a ruby laser (~694 nm wavelength) generates ultraviolet light (~347 nm), which is twice the frequency of the incident light.
Furthermore, if we consider a bichromatic optical field (with the angular frequencies ω1>ω2) incident on a χ(2) material, we can demonstrate that the source term T of the nonlinear wave equation has four components, namely 2ω1, 2ω2, (ω1+ω2) and (ω1−ω2), which give rise to sum- and difference-frequency generation. More generally, one can demonstrate the ability to generate extra terms so that the process can be seen as a natural consequence of the nonlinear interactions. These complex processes may be used, for example, to generate optical frequency combs , , .
3.2 Origin of nonlinear acoustic effects
The origin of acoustic nonlinearities is qualitatively similar, but different in terms of the magnitude of the effect. For a sound wave travelling through air or water, the vibrations of the particles of the medium are best described by longitudinal waves . (Optical waves are transverse, but this distinction does not affect the generality of our discussion because in solids sound can be transmitted as both longitudinal waves and transverse waves, also giving rise to strong nonlinear effects .)
By analogy with the nonlinear wave equation for optical waves, we can write a nonlinear wave equation for acoustic waves by considering the system of Euler equations . We can arrive to the nonlinear wave equation ∂2v/∂t2−c02∂2v/∂x2=−Ta, where v is the particle velocity, c0 is the ambient speed of sound and the source term Ta=L2(v2)+L3(v3)+L4(v4)+…. Here, L2(v2) contains the terms in powers of 2, L3(v3) the terms in powers of 3 and so on. For example, in the explicit form L2(v2)=v2+ 0.5(γ−1)c0−2(∫[∂v/∂t]dx)2 , where γ=cp/cv, cv and cp being the specific heat capacities at constant volume and pressure, respectively. In the following we will demonstrate that γ defines the strength of the nonlinear response.
We find solutions for a sinusoidal plane wave v=v0sin(ωat) originating from the point x=0 and propagating along the x-coordinate direction. Here, ωa is the acoustic angular driving frequency and t is time. By expressing v=v(1)+v(2)+v(3)+…, in the second approximation we can write ∂2v(2)/∂t2−c02∂2v(2)/∂x2=L2v02sin2ωa t′+L3v03sin3ωat′+…, where t′=t−x/c0. Thus, as in the case of optical waves, the nonlinear acoustic medium acts as an anharmonic oscillator and its response contains not just the fundamental frequency ωa, but also multiples of this frequency (and also a dc acoustic component).
Similar to the case of optical waves, considering only the quadratic acoustic nonlinearity term and a bi-harmonic acoustic wave with the angular frequencies ω1(a)>ω2(a), we can demonstrate that the source term Ta of the nonlinear acoustic wave equation has four components, namely 2ω1(a), 2ω2(a), (ω1(a)+ω2(a)) and (ω1(a)−ω2(a)), which give rise to sum- and difference-acoustic frequency generation. More generally, one can demonstrate the ability to generate extra terms so that the process can be seen as a natural consequence of the nonlinear acoustic interactions.
Let us now assume that, without loss of generality, the plane acoustic wave travels in the positive x-direction. Using the equations derived above, one can show  that the propagation speed for any particular point of this wave is given by the local value of (c0+βv), where β=(γ+1)/2 is the nonlinear acoustic parameter. Hence, the points with v>0 (e.g. wave crests) travel faster than c0 and correspond to areas of compression. (Recall that longitudinal acoustic waves can be described as alternating areas of compression and rarefaction.) Conversely, the points with v<c0 (e.g. wave troughs) travel slower and correspond to areas of rarefaction. This is opposed to the case of an idealised, linear wave propagation with β=0, where all points of the wave travel at c0 (compare the solid and dashed curves in Figure 2A).
Thus, an initially sinusoidal sound wave undergoes deformations as it propagates, which implies that its initial monochromatic spectrum (Figure 2B) gets enriched with new higher harmonic frequencies ,  (Figure 2C). This spectral enrichment is stronger in acoustic media with larger values of the nonlinear acoustic coefficient β, the physical meaning of which is as follows. The rate of deformation of the propagating acoustic wave depends on the distance between its crests and troughs, and this depends on both the initial particle velocity v0 and the value of β for the particular fluid. For example, water has β=3.5 and air has β≈0.7 , . As a result of this dependence, a fluid with a higher β will exhibit more rapid wave deformation than a fluid with a lower β, provided the same initial velocity v0 is applied to both.
Significantly, an acoustic medium with β=0 is just an idealisation. Hence, even in the low-power acoustic regime characterised by small values of v0<<c0, there is no ideal linear propagation and, as with the case of nonlinear optics, the acoustic medium responds anharmonically.
Despite isolated attempts to directly compare the strength of nonlinear optical and nonlinear acoustic effects , this task is, in general, challenging due to the different nature of light and sound. Indeed, one could compare the peak amplitude of the electric field E0 with the peak particle velocity v0. However, nonlinear optical effects start to manifest themselves only when optical electric field amplitudes reach high values of 108…1010 V/m . This is in stark contrast with nonlinear acoustics where systematic observation of practicable nonlinear acoustic effects is always possible, even with a very small v0.
Alternatively, one could also compare both types of nonlinearity by comparing the mathematical role of the nonlinear acoustic parameter β and the nonlinear optical susceptibilities χ in their respective wave equations. While successful in some particular cases (e.g. when optical nonlinearities of dielectrics are compared to acoustic nonlinearities of bulk fluids ), this approach fails in other important situations. The strength of nonlinear acoustic processes also depends on the compressibility of the medium, which is a measure of the relative volume change of the medium as a response to a pressure , the density of the medium, temperature and other factors. Fluid media may also form interfaces with solids and other fluids with different density and compressibility. Consequently, the nonlinear response additionally becomes a function of surface and interfacial tension .
For example, this is the case of gas bubbles in water that are characterised by a giant value β that is several orders of magnitude larger than the values of β for water and air . Consequently, whereas in bulk water one can readily observe the generation of five or so higher frequency harmonics of the incident sound (Figure 2C), in water containing gas bubbles one may generate up to 15–20 high-frequency acoustic harmonics by applying the same sound pressure as in the case with bulk water (Figure 2D). These giant nonlinear properties of gas bubbles have no analogues in optics and they arise from an intrinsic property of bubbles to be easily compressible, which will be discussed in more detail in Section 4.
It is noteworthy that, by analogy with optics, in acoustics one can achieve ultra-fast nonlinear response. For example, acoustic pulses as short as a few picoseconds can be generated with pulsed laser light (for a review see, e.g. ). By analogy with ultra-fast nonlinear optics, ultra-fast acoustic allows observing fundamental nonlinear phenomena such as acoustic solitons .
Initially, ultra-fast nonlinear acoustic response was achieved in solids and applications such as nanoscale sonar operating at unprecedentedly high acoustic frequencies of up to a few THz were conceived . Solids, especially in the form of single crystals and at low temperatures, are good conductors of ultra-high-frequency sound (hypersound) . Here, for example, absorption of energy from the pump beam leads to heating of the surface region. Then, coherent strain is generated by thermal expansion of the lattice of the material. The generated strain travels into the material. The wave is also partially reflected at the material interface, which allows detecting it with a delayed, weak probe pulse. At room temperature, the propagation of the high-amplitude acoustic wave in the material gives rise to a nonlinear effect, which leads to the formation of a shock wave (see Section 6). The shock wave can also be detected with a weak probe pulse. At low temperatures, the wave develops into a train of acoustic solitons accompanied by a high-frequency tail.
However, the excitation of hypersonic acoustic waves in fluids is more challenging or even impossible at normal ambient conditions. For example, in air at normal atmospheric pressure and room temperature, acoustic frequencies of >1 GHz correspond to an acoustic wavelength of <3.4×10−7 m, which means that this length becomes comparable with the mean free path of molecules in air under these conditions. Since acoustic waves can mostly propagate in the medium when their wavelengths are considerably longer than the mean free path in gases (or longer than the interatomic distance in liquids and solids), hypersound waves do not propagate in air and gases at normal atmospheric pressure .
However, hypersound can propagate in liquids, even though its attenuation is relatively high and the propagation range is smaller as compared with the case of ultrasound. Consequently, picosecond ultrasonic effects have also been observed in liquids such as water , biological cells  and liquid metals .
Although, to the best of our knowledge, there are no further works explicitly targeting ultra-fast nonlinear acoustic response of gas bubbles, there have been demonstrations of enhanced acoustic cavitation of gas bubbles in the picosecond acoustic regime . Because cavitation is a nonlinear process (see Section 6), it is conceivable that the giant acoustic nonlinearity of gas bubbles can also be accessed in the ultra-fast regime.
3.3 Sensing of nonlinear acoustic waves with light
When a sound wave propagates in a fluid medium, the periodic regions of compression and rarefaction of the particles of the medium lead to changes in the medium density, which also results in changes in the optical refractive index n of the medium . Light perceives these changes as a diffraction grating that moves with the speed of sound .
In general, this model properly describes the coupling between light and nonlinear acoustic waves in fluid media . This model is also useful for the analysis of the interaction of light with acoustic nonlinearities of gas bubbles and liquid droplets.
However, a proper description of many solid-state opto-mechanical systems often requires taking into account electrostriction and piezoelectricity , , , . A piezoelectric displacement in the medium is linearly proportional to the electric field. In contrast, electrostriction leads to contractions and expansions of the dielectric medium in proportion to the square of the electric field. For example, in stimulated Brillouin scattering , the interaction of the electric fields of two optical waves generates a density change in the medium. As a result, the optical waves are scattered by the moving refractive index grating and this leads to the nonlinear coupling between them , , , , , .
Although electrostriction also occurs in fluids , the role of this effect in our further discussions of the interaction between light and acoustic waves in fluids is mostly secondary. For example, this is because we do not consider the scenarios where electrostriction is induced by laser light pulses, which would be the case of laser-induced generation of acoustic waves and cavitation .
The change in the refractive index due to sound opens up opportunities to modulate the optical signal with an acoustic wave. If the acoustic wave exhibits nonlinearities, these will also be imprinted onto the modulated optical signal, thereby effectively enabling the conversion of the acoustic nonlinearity into new optical signals (Figure 1B). For example, this process may reproduce the effect of the nonlinear optical generation of new optical frequencies (Figure 1A).
The physics of the light-sound interaction is in particular important for the development of hydrophones. A hydrophone is a microphone designed to be used underwater for detecting underwater sound. In conventional hydrophones, the acoustic energy mechanically moves a transducer that in turn generates a voltage at the acoustic frequency. Optical hydrophones , , ,  use frequency modulation of a light beam that Doppler shifts in response to acoustic pressure variations, which also gives rise to the effect of Brillouin light scattering (BLS) from sound . They have an increased sensitivity, greater discrimination against noise and improved directionality as compared with conventional hydrophones . Although the invention of laser played an important role in the development of optical hydrophones , in practice, optical hydrophones may also operate with incoherent sources of light , provided the Doppler shift is large enough compared to the linewidth of the source of incoherent light. For example, this condition can be satisfied with gas discharge lamps .
Optical hydrophones are also more compact than mechanical ones. The dimensions of state-of-the-art miniaturised optical ultrasonic hydrophones are set by the wavelength of ultrasound in water (~150 μm at 10 MHz) , , , . The ability to use incoherent light additionally contributes to the development to compact, reliable and inexpensive optical hydrophones.
Optical hydrophones also provide us with a pathway towards using acoustic nonlinearities in optics. However, in this case an optical hydrophone needs to have an effective mechanism for light-sound interaction. We showed that an enhanced strong light-sound interaction can be achieved with a single silver nanorod ,  operating as an antenna for light , , . (See Section 5 for a more detailed discussion of optical antennas.) The nanorod is immersed into water and insonated by MHz-frequency range ultrasound. In contrast to well-researched GHz-frequency range resonant structural deformations of metal nanostructures such as rods, discs, crosses and cubes , , , , in our scenario we investigate the quasi-static regime. This regime of transduction is interesting as some of the dimensions of the nanorod are small compared to the wavelength of light (e.g. the cross-section of the nanorod may be 30 nm×30 nm and it may be 200–400 nm long), and all of them are very small compared to the wavelength of MHz-range sound (~150 μm in water). This means that a plasmonic nanoantenna operating as a hydrophone explores a qualitatively different regime from more conventional hydrophone technologies and nano-opto-mechanical systems , , , .
A special design of the nanorod is essential to enable high sensitivity of light to sound, that is, to changes in the optical refractive index due to sound. The usual operating regime of the nanorod as an optical dipole is unsuitable for sensing of MHz-range ultrasound. This is because in the dipole mode light is localised near the edges of the nanorod (Figure 3A) and hence the interaction of light with sound is small. However, the sensitivity of light to sound is dramatically increased when light bounces back and forth along the nanorod length (Figure 3A). In this case, the nanorod operates as a nanoscale Fabry-Perot resonator, thereby effectively increasing the light-sound interaction.
These increased sensing properties are confirmed by rigorous simulations of BLS . In that particular simulation run, we assume that all acoustic interactions in water are linear. Figure 3B shows the simulated BLS spectrum consisting of the central Rayleigh peak and two Brillouin peaks shifted by ±10 MHz, which is the frequency of the incident sound. Three scenarios are considered: (i) the incident light interacts with sound in bulk water and the nanorod is absent, (ii) the nanorod is present but it is tuned on its fundamental optical mode and (iii) the nanorod operates in the higher-order mode regime. One can see that the intensity of the Brillouin peaks is negligibly small in the first two cases, but the peak intensity, and therefore the sensitivity of light to sound, is much higher when the nanorod is tuned on its higher-order mode. These important observations will be used in the following section to discuss coherent generation of multiple optical frequencies from nonlinear ultrasound.
We note that the intensity of the Brillouin peaks in Figure 3B can be increased by using a regular array of identical plasmonic nanorods fabricated, for example, on top of an optical fibre. Alternatively, one may use a colloidal solution of nanorods. In both cases, a considerable (orders of magnitude) enhancement can be achieved because of the collective plasmonic effects and antenna-like properties of periodically arranged plasmonic nanoparticles (Section 5). In Section 4.2, we also demonstrate that the sensitivity of modern commercial laser Doppler vibrometers suffices to resolve weak peaks appearing in the optical spectrum due to the interaction of light with sound.
3.4 Optical frequency comb generated from nonlinear ultrasound
Recall that an optical frequency comb is a spectrum consisting of a series of discrete, equally spaced peaks that have a well-defined phase relationship between each other, a phenomenon called phase locking , , . In integrated photonic devices, optical combs are often produced by cascading nonlinear-optical wave mixing processes in a high-Q resonator , .
In our approach, when optical frequency combs are produced from acoustic nonlinearities, we first obtain equally spaced acoustic frequency peaks with the same phase, and then convert these peaks into an optical spectrum that inherits the equal frequency spacing and phase locking .
Two physical effects describe the light-sound interaction that underpins the generation of the optical comb: (i) the Bragg reflection of light from the moving diffraction gratings produced by sound waves and (ii) the Doppler shift – change in the frequency of light as the gratings move with respect to the source of light. In a quantum mechanical picture, the interaction of light with sound is an inelastic process where the photon energy is not conserved. As a result, the interaction of light with sound creates Brillouin peaks shifted from the central (Rayleigh) peak by the frequencies of all nonlinear sound waves present in the system (Figure 4A). The central peak is due to the elastic scattering of light (where photon energy is conserved).
As shown in , all the peaks in Figure 4A are phase-matched, which is an essential feature of an optical frequency comb . Consequently, similar to mode-locked optical frequency combs, in the time domain the signals at different frequencies should add constructively at one point resulting in a train of optical pulses spaced by Tp=1/δ , where δ is the frequency spacing between the peaks in Figure 4A. We simulate this physical picture by Fourier-transforming the spectrum in Figure 4A and obtaining time-domain optical intensity signals, which produces well-defined single optical pulses (Figure 4B).
Spectral tuneability of optical frequency combs is required for many real-life applications , , . Spectral composition of optical frequency combs generated from sound can also be tuned by changing the acoustic frequency. For example, in the particular scenario in Figure 4A, the spacing between the comb peaks equals 1 MHz. However, it can be tuned by changing the frequency of the incident sound in a broad spectral range. The lower limit of this range is defined mostly by the linewidth of the optical source. The upper limit is of the order of several GHz, which are the frequencies of hypersound .
In general, the ability to selectively use the broad kHz-to-GHz acoustic frequency spectrum in the optical domain may not only be exploited in frequency combs but would also be useful for sensing, material characterisation and observation of atmospheric phenomena , , , , . Audible acoustic frequencies are also of high interest because of potential applications in optical analysis of acoustic signatures from systems such as cells or molecular machines as well as in localised photoacoustic microscopy .
4 Giant acoustic nonlinearities of bubbles and droplets
4.1 Mechanism of giant acoustic nonlinearities
Liquids have little free space between molecules, which implies that they are not easily compressible. Gases, on the other hand, are easily compressible. Together with the medium’s density ρ, the compressibility of the fluid medium βs defines the speed of sound in the medium as c0=(ρβs)−1/2. We can assume that for water ρ≈1000 kg/m3 and βs≈4.6×10−10 Pa−1. For air, however, ρ≈1.2 kg/m3 and βs≈7×10−6 Pa−1, which implies that sound travels faster in water (1474 m/s) than in air (345 m/s).
This distinction also leads to stronger acoustic nonlinearities in systems composed of two fluids with different compressibility and density, and here gas bubbles in water are the most striking example. Let us consider a single gas bubble in otherwise bulk water. When an acoustic wave propagating in bulk water reaches a gas bubble, due to the low density of the gas inside bubble and its high compressibility as compared with water, the amplitude of the acoustic wave is considerably amplified inside the bubble, resulting in dramatic changes in the bubble volume . (Here, we may assume that the bubble maintains its spherical shape .) The changes in the bubble’s volume result in large acoustic wavefront deformations. As defined above, the deformation of the acoustic wave is quantified by the nonlinear acoustic parameter β. Therefore, whereas the two components of the bubble – water and air – have β=3.5 and 0.7, respectively, when taken separately , an air bubble in water is characterised by β≈5000 . This phenomenon is called giant acoustic nonlinearities . For example, whereas in water one normally observes the generation of about five higher-order acoustic harmonics (Figure 2C), in water with gas bubbles it is possible to observe 15–20 harmonics with the same applied acoustic power and other equal conditions (Figure 2D).
The acoustic response of a single bubble in bulk liquid is described by the nonlinear Rayleigh-Plesset equation . In its simplest form, this equation models a bubble oscillating in an inviscid and incompressible liquid. More complex forms of this equation may take into account the liquid viscosity and surface tension .
For example, Figure 5A shows the solution of the Rayleigh-Plesset equation for a 100-nm-radius air bubble in water excited by an ultrasound pulse with a Gaussian-enveloped sinusoid with the frequency fa=50 MHz. This frequency is detuned from the natural resonance frequency of the bubble f0≈30 MHz, which in the case of water can conveniently be defined as f0R0≈3 m/s , where R0 is the radius of the bubble at rest.
One can see that the radius of the bubble (blue solid curve) changes in response to the driving pressure (red dashed curve) with a phase lag due to the inertia of the surrounding water , . The maximum (minimum) value of the radius reaches ~130 nm (~80 nm). The bubble continues pulsating with a smaller amplitude when the driving pressure signal is turned off. This is because damping is not taken into account in this model. This resonant tail disappears when acoustic losses are taken into account .
Significantly, the behaviour of the pulsating bubble modulates the intensity of light scattered by the bubble. This is shown in Figure 5B, which plots the calculated intensity of light (blue solid curve) transmitted through a hole in a silver film. The hole is filled with water and it contains a 100-nm-radius air bubble. One can see the shape of the acoustic signal (red dashed curve) imprints onto the light intensity transmitted through the hole .
We now turn our attention to the fact that, in general, the physics of liquid droplets and gas bubbles is similar. Indeed, the surface tension of water results in the wall tension required for the formation of the bubble. The natural tendency to minimise the wall tension forces the bubble to maintain its spherical shape. The same surface tension effect determines the shape of liquid droplets. Because the droplet can be easily deformed, it tends to maintain a spherical shape due to the cohesive forces of the surface layer. As in the case of the bubble, the spherical shape minimises the necessary wall tension.
Thus, liquid droplets oscillate in response to sound similar to gas bubbles and they also exhibit strong nonlinear properties analogous to those of oscillating gas bubbles . However, in the case of droplets, we exploit capillary oscillations that arise because of a competition between inertia of the liquid and surface tension , . Unlike acoustic waves that can propagate in all media, capillary oscillations are unique to liquids and they may be driven, for example, electrically or mechanically . Significantly, in contrast to gas bubbles that oscillate and maintain their spherical shape, acoustically driven liquid droplets often assume complex, non-spherical 3D shapes (Figure 6C) corresponding to specific capillary oscillation modes .
Moreover, in the recent work  it has been demonstrated that liquid droplets capped with an elastic film can also assume elliptical and even nearly perfect square shapes (Figure 6D). The aspect ratio of elliptical droplets and the size of square droplets can be tuned by adjusting the tension. In  it was suggested that liquid droplets with unusual shapes would find applications in liquid lenses and liquid patterning systems. Yet, excitation of such droplets with sound may lead to previously unknown oscillation regimes and interesting fundamental phenomena.
It is noteworthy that in droplets we access nonlinear properties because of the softness of the liquid. When a solid-state medium is deformed, the resulting restoring force is due to the stiffness (Young’s modulus) of the material. However, in liquid surfaces, the restoring force relies on surface tension . Moreover, the speed of the capillary wave is three orders of magnitude lower than the speed of the acoustic waves in the same liquid , . Because of all these distinctions, the frequency of capillary oscillation of a liquid droplet (Figure 6A) is about three orders of magnitude lower than that of the correspondent acoustic mode (Figure 6B) of the same droplet , .
4.2 Interaction of light with droplets and gas bubbles
The optical refractive index of air and water is n=1 and n=1.33, respectively. This refractive index contrast allows detecting single spherical bubbles by light extinction and Mie scattering . A water droplet surrounded by air also has the same refractive index contrast and therefore can be detected using similar techniques. Other optically transparent liquids may be used to form droplets, such as octane (n=1.4), and may form droplets in the form of an emulsion when surrounded by an immiscible liquid.
Thus, when a gas bubble (liquid droplet) interacts simultaneously with light and sound, large changes in the volume of the bubble (or the shape of the droplet) lead to large changes in the behaviour of light . More specifically, the strength of light scattering from sound becomes a function of the complex nonlinear process of compression and rarefaction of the gas inside the bubble (or capillary oscillations of the liquid droplet). This effectively converts nonlinearity into optical signals containing the acoustic frequency components.
Although the conversion of acoustic nonlinearity of bubbles and droplets into the optical domain has not been in the focus of previous works, there have been several highly relevant demonstrations of optical signal generation from capillary waves in liquid droplets , . The most relevant results demonstrated in those works will be overviewed below.
Before we do that, we note that conceptually similar functionality can also be achieved with well-known acousto-optic modulators that use the acousto-optic effect to diffract and shift the frequency of light using sound waves. For example, in  a nonlinear acousto-optic modulator operating at GHz acoustic (vibrational) frequencies was demonstrated. However, that nonlinearity originates from resonant vibrations of metal nanostructures. This is a different kind of nonlinearity and, although it has important practical applications in its own right , it is more difficult to access it with low power due to considerable stiffness of solid-state structures. This is in stark contrast with bubbles and droplets whose softness allows accessing their nonlinear properties with low power.
In  optical tweezers were used to manipulate an octane droplet immersed into water. An optical waveguide placed in close proximity of the droplet was used to deliver the pump light. Radiation pressure from the light circulating in the droplet gives rise to capillary oscillations. Shape changes due to these oscillations result in a Doppler shift of the pump light, thereby giving rise to Stokes emission lines when the pump light is slowly scanned through one of the droplet’s optical resonances.
In the resulting frequency spectrum (Figure 7A), the strongest spectral line corresponds to the fundamental mode of the droplet oscillation. (In the notation used in this work, this mode corresponds to the l=2 mode, see the caption to Figure 6.) The oscillations of the fundamental mode give rise to the second, third and fourth high-frequency harmonics by virtue of nonlinear interactions. In agreement with the theoretical prediction for light scattering from oscillating gas bubbles in water (Figure 7B), the corresponding signal in the time domain is represented by a train of well-defined pulses.
A similar result is observed when the direction of the pump light scan is reversed (Figure 7C and D). However, in this case the third higher-order capillary oscillation mode of the droplet is excited (l=4 in Figure 6C). This is because higher-optical-quality-factor modes tend to excite oscillations at a lower rate . The second higher-order mode (l=3 in Figure 6C) is absent because even modes have their node near the equator line  and therefore exhibit smaller opto-capillary coupling.
In , the authors proposed to use this effect to build a ripplon laser. The results presented in their paper also open up novel opportunities to generate optical frequency combs from capillary oscillations. It is noteworthy that capillary oscillations may also be excited mechanically (e.g. by ultrasound propagating in a liquid that surrounds the droplet) or electrically (e.g. by using an electrowetting technique ).
More generally, the ability to use different excitation mechanisms offers additional degrees of freedom in the design and control of novel multiphysics devices. For instance, considerable attention has recently been paid to exploiting novel regimes of strong coupling between different entities such as photons, phonons and magnons , . It has been suggested that multi-resonant photon-phonon-magnon interaction can be achieved with magnetofluidic droplets which combine the softness of fluids with the ability of solids to support resonances of electromagnetic, acoustic and spin waves .
Capillary oscillations of liquid droplets are also often detected with laser Doppler vibrometry (LDV). Commercial LDVs often employ a Mach-Zehnder interferometer (Figure 8A) and usually allow characterising the motion of a fluid-fluid interface at frequencies of up to 100 MHz and displacements of as little as a few tens of picometres . Oscillations of the droplet can be excited by either a thickness-mode piezoelectric transducer (PZT, see Figure 8B) or a surface acoustic wave device .
In , the behaviour of capillary waves in a 2 μl sessile hemispherical water droplet excited by high-frequency acoustic waves oscillating at 500 kHz and 20 MHz was investigated. The droplet was placed on a solid substrate that vibrates to generate capillary oscillation modes (Figure 8B).
In the measured displacement spectral density (Figure 8C), the resonance peaks at ≥200 Hz are in good agreement with the Lamb model of elastic resonance of a spherical capillary surface . In the 100 Hz–2 kHz frequency range, the spectrum exhibits a slope predicted by the wave turbulence theory . Significantly, one can also see higher-order resonances from the excitation frequency of 500 kHz upward in a harmonic cascade f, 2f, 3f, ..., 7f, and perhaps beyond; however, these higher-order resonances are not detected due to the upper limit of the LDV’s measurement range at that particular resolution.
5 Plasmon enhancement of light-sound interactions
5.1 Plasmonic nanoantennas and nanoparticles
By analogy with a conventional roof-top antenna, an optical nanoantenna emits, receives and, more broadly, controls light with nanoscale elements that are much smaller than the wavelength of the incident light , , . In particular, nanoscale dimensions and subwavelength operating regime allow using optical nanoantennas to enhance light-matter interaction which, along with other many interesting phenomena useful for sensing, spectroscopy and imaging, leads to strong nonlinear optical effects , , , , , , , , , , , , , . Optical nanoantennas may be made of metals (most often of gold and silver) or dielectric and semiconductor materials. The former are called plasmonic nanoantennas and the latter are known as dielectric nanoantennas .
In plasmonic nanoantennas, enhancement of nonlinear optical phenomena is achieved due to strong light-matter interactions enabled by localised surface plasmon resonances. Plasmon resonances are collective oscillations of the electron charge around metal nanoparticles in resonance with the frequency of light . Localised plasmon resonances lead to trapping of light in close vicinity to nanoparticles, thereby resulting in stronger local optical electric fields and hence stronger nonlinear optical processes (also see Section 2).
The physics of the enhanced light-matter interaction in dielectric nanoantennas is different . Such antennas do not have metal parts and therefore they do not suffer from high light absorption losses and Joule heat production typical of plasmonic nanoantennas. Indeed, plasmons are oscillations of the electron density in metals, and therefore energy is lost and heat is produced when oscillating electrons collide with the metal lattice. In general, the stronger the optical electric field in the metal due to the plasmon resonance, the stronger the losses and heat.
However, light localisation and enhancement in dielectric nanoantennas may be several orders of magnitude smaller than in their plasmonic counterparts. Nevertheless, the achievable enhancement levels can enable strong nonlinear optical effects, especially because light is localised inside the body of the nanoantenna .
The interaction of light with sound in liquids may be increased when gold or silver nanoparticles operating as optical nanoantennas are present in the liquid . Returning to the example of the optical frequency comb generation from ultrasound propagating in bulk water in the nonlinear regime. As shown in Figure 4A, the presence of a single silver nanorod particle in water gives rise to an ~100-fold increase in the intensity of the peaks in the generated BLS spectrum. An ensemble of such particles acting as an array of plasmonic nanoantennas will further increase the magnitude of the peaks. For example, this may be achieved due to high light intensity “hot spots” in the gaps between the nanoparticles .
As demonstrated above, the interaction of light with acoustically driven gas bubbles in water leads to generation of new optical frequencies. The nanoantenna-assisted increase in the sensitivity of light to sound will also work in this scenario. Indeed, microscopic gas bubbles in liquids can be integrated with metal nanoparticles that are located in the shell of the bubble, thereby leading to the effect of a “metal bubble” , . Such gas bubbles have already been employed as a bimodal biomedical imaging agent that provides contrast for both ultrasonic and optical imaging. In particular, to be useful for optical imaging, the optical resonance of “metal bubbles” falls within the visible-to-near-infrared spectral range .
5.2 UV plasmonic liquid-metal nanoparticles
Thus far, we have discussed the opportunities to use plasmonic resonances to increase the interaction of light with giant acoustic nonlinearities of gas bubbles in liquids. However, it would also be beneficial to combine the plasmonic properties with the capillary oscillations of liquid droplets. Unfortunately, however, water and other common liquids cannot support surface plasmon resonances.
To circumvent this problem, one may employ room-temperature liquid metals such as gallium and its alloys , , . Such metals are excellent candidates for the role of plasmonics nanodroplets (also called liquid-metal nanoparticles).
Nanoparticles made of non-noble metals such as gallium have recently attracted significant interest due to promising plasmonic applications in the UV-visible spectral range , , , , , , , , , , . Unlike gold and silver, gallium has a Drude-like dielectric permittivity function extending from the UV range through the visible and, mostly in the liquid state, into the infrared spectral region , . Due to their high environmental stability and excellent mechanical properties, gallium and its alloys are highly relevant for many emerging applications in electronics, micro-mechanics and chemistry , , .
Significant effort has recently been made to synthetise colloidal liquid-metal nanoparticles , . Consequently, we investigated eutectic gallium-indium (EGaIn, 75% Ga 25% In by weight) nanoparticles suspended in ethanol (Figure 9A and B) and demonstrated their strong plasmonic resonances in the UV spectral region (Figure 9C) . EGaIn has negligible vapour pressure and low viscosity , , which is important for excitation of high-intensity capillary oscillation modes. Consequently, it allows strongly modulating the UV plasmon spectra of liquid-metal nanoparticles with ultrasound (Figure 10A).
Figure 10B shows the calculated light scattering spectra of the liquid-metal nanoparticles tuned on one of the attainable capillary oscillation modes. One can see a dramatic change in the UV plasmonic spectrum of the nanoparticles that oscillates as a function of the driving acoustic pressure. These changes enable the conversion of nonlinear acoustically driven oscillations into optical signals in the UV range.
6 Regimes of extreme nonlinear acoustic events
As demonstrated in Section 4.2, radial oscillations of gas bubbles and effective acoustic response of water containing multiple gas bubbles are highly nonlinear. More broadly, a bubble in a liquid can be considered as a nonlinear dynamical system exhibiting complex resonances, multistability, bifurcations to chaos and other complex processes . For example, by driving a 10-μm radius-at-rest gas bubble with a considerable sound pressure amplitude of 300 kPa and a high frequency of 600 kHz, one may achieve chaotic radial oscillations represented as a complicated set with a fractal structure known as a strange attractor . The oscillations of a gas bubble can also switch to a new oscillation regime with twice the period of the original oscillation system, a process called period doubling . Moreover, with an even stronger acoustic driving pressure, the response of the bubble may develop hysteresis.
Consequently, the conversion of these complex nonlinear phenomena into the optical domain may result in new intriguing fundamental physics that would complement the purely optical, complex nonlinear phenomena such as optical bistability , ,  and low-power optical frequency comb generation exploiting optical period-doubling bifurcation .
Another intriguing possibility associated with the giant nonlinearity of gas bubbles is to exploit extreme processes that accompany the effect of bubble cavitation . Cavitation is the formation of vapour bubbles in a liquid due to forces acting on the liquid. Very often the cavitation occurs when a liquid is subjected to rapid changes of pressure that cause the formation of bubbles where the pressure is relatively low. When subjected to higher pressure, the bubbles implode and can generate an intense shock wave. Experimental results also demonstrate that collapsing bubbles reach temperatures of ~5000 K, pressures of ~1000 atm, and heating and cooling rates of above 1010 K/s . These events can create extreme physical and chemical conditions in otherwise cold liquids, such as ice breaking .
Cavitation also occurs, and can be controlled, in the presence of an acoustic pressure wave. When microscopic gas bubbles are present in the liquid in which the acoustic wave propagates, the bubbles are forced to oscillate due to the applied acoustic pressure. With sufficiently high acoustic pressures, the bubbles will first grow in size and then rapidly collapse , .
A gas bubble collapsing near an interface exhibits especially interesting properties , , . In general, a bubble collapsing near a solid interface develops a jet directed towards the interface (blue curves in Figure 11A). This jet is so strong that it can lead to cavitation damage of the solid surface , . In the case of a free surface and fluid-fluid interfaces , the jet is directed mostly away from the surface but the free surface or the fluid-fluid interface undergoes considerable deformations due to the evolution of the bubble (Figure 11B). The collapse of a bubble near a free surface or a fluid-fluid surface may also result in the formation of new bubbles . The results in Figure 11 were produced by using computational code that simulates asymmetric bubble cavitation . The fluid is assumed to be incompressible, inviscid and irrotational, and surface tension is neglected. The Navier-Stokes equations are solved using the boundary integral method .
The evolution of the bubble near an interface is highly nonlinear , . This implies that the displacement of both the surface of the bubble and free surface is a nonlinear function of the driving acoustic pressure. Light can readily sense these displacements. Furthermore, an even higher sensitivity may be achieved in a system where a gas microbubble collapses near a liquid-metal surface. Here, the liquid-metal is expected to behave as a liquid that responds the oscillations and collapse of the bubble similar to the scenario in Figure 11B. However, at the same time the liquid-metal interface may support surface plasmon modes , which will allow for efficient sensing of small, potentially nanometre-range, fluctuations of the metal surface.
Finally, in this section we highlight the idea to control light and effectively recreate some nonlinear optical effects in photonic crystals interacting with a shock wave , . Photonic crystals are periodic dielectric structures that have a band gap that forbids propagation of light in a certain frequency range. This property opens up new opportunities to manipulate light and achieve effects that are impossible with conventional optics .
Like an ordinary wave, a shock wave carries energy and can propagate through a medium. However, it is characterised by an abrupt, nearly discontinuous change in pressure, temperature, and density of the medium. When a shock wave interacts with a photonic crystal, it creates a shock-like modulation of the dielectric medium, which in turn dramatically affects light trapped in the photonic crystal. As a result, new physical effects arise when light interacts with the moving shock-like modulation of the dielectric medium.
One of the effects predicted in shocked photonic crystals is the generation of new optical frequencies. This effect reproduces certain nonlinear optical interactions without the need for high-power laser light and phase matching. It was theoretically demonstrated that new frequencies can arise from a single-frequency input light of any intensity . It was also shown that this process is tuneable because it depends on the frequency range of the photonic band gap, which in turn is defined by the geometry of the photonic crystal structure .
Another predicted phenomenon is the reversed Doppler effect . The Doppler effect predicts that light emitted by a source onto an object moving towards the source will be reflected with a higher frequency. This effect breaks down for the light reflected from a shocked photonic crystal leading to anomalous Doppler shifts both in sign and magnitude. Normally, the reversed Doppler effect is observed in plasmas that propagate at near-relativistic speeds or it can be achieved in metamaterials (see  and references therein). However, in the shocked photonic crystal, this effect is possible at non-relativistic conditions with linear optical materials.
Although photonic crystals are made of solid dielectric materials , they can be integrated with microfluidic systems , which provides a connection with liquid-state platforms discussed in this review paper. Moreover, the concept of photonic crystals also extends to the acoustics domain where such structures are known as phononic crystals that have a band gap that forbids the propagation of sound in a certain acoustic frequency range . This provides another possibility of exploiting shock waves that arise in liquids due to the acoustic nonlinearity.
7 Conclusions and outlook
Modern nonlinear photonics is mostly impossible without lasers. In this review paper, we have shown how a judiciously engineered interaction between sound, capillarity and low-power incoherent light provides access to exciting and yet largely unexploited nonlinear phenomena of non-optical origin. Understanding, controlling and harnessing these new nonlinear phenomena promises to remove the restriction of high-power laser light operation. This also enables compact and readily deployable photonic devices for telecommunications, sensing and biomedical applications.
For example, although lasers have important applications in biomedicine (e.g. in imaging, skincare, and kidney stone removal ), classical nonlinear optical processes have had relatively little impact in the field of biomedicine  because high laser powers required to achieve optical nonlinearites cause photodamage to living cells and tissues , . Hence, effectively nonlinear, low-power photonics devices would usher a new class of biologically friendly photonic sensors and imaging systems.
We have also demonstrated that acoustic nonlinearity can be accessed with gas bubbles and liquid droplets that oscillate when subjected to acoustic waves. Significantly, gas bubbles have already found important applications in biomedicine and were approved by the Food and Drug Administration and similar organisations as ultrasound contrast agents operating inside the human body , . This dramatic progress will provide a considerable advantage for novel biomedical optical devices and imaging modalities exploiting the giant nonlinear acoustic properties of bubbles.
Of course, the use of fluids as a constituent material for photonic devices presents a number of technological challenges, such as the integration of liquid-state elements into traditional solid-state microphotonic systems , . Nevertheless, liquid-state devices can offer a number of potentially transformative advantages for microphotonic systems , , . Moreover, there have been successful examples of liquid-state and gas-state devices that are compatible with conventional solid-state photonics such as liquid-state optical lenses , liquid-core optical fibres used in spectroscopy  and gas-filled optical fibres used to generate optical frequency combs , . In our Ref. , we further capitalise on the progress in the liquid-core optical fibre development and suggest to employ liquid-core fibre filled with water containing gas micro-bubbles.
Further technological challenges may include the need to align solid metal nanoparticles inside a liquid to optimise the excitation of plasmon resonance modes. This may be required when plasmonic effects are exploited to increase the strength of the light-sound interaction. In solid-state plasmonic devices, the position and also the orientation of nanoparticles with respect to the polarisation of the incident light are fixed. In a liquid-state devices, however, external disturbances and temperature changes may result in movement and changes in the nanoparticle orientation. To circumvent these effects, one may employ magneto-plasmonic nanoparticles made of ferromagnetic metals , , . By applying an external static magnetic field to the resulting magneto-fluid, one may control not only the positioning of the nanoparticles but also their collective optical response, thereby resulting in intriguing magneto-optical effects ,  and important applications in biomedical imaging (see  for a relevant discussion).
Liquid-state devices also have limitations originating from the acoustic and thermal properties of liquids. First, in real liquids, the oscillations of gas bubbles are damped due to acoustic losses, action of viscous stresses at the gas-liquid interface, thermal conduction across the interface, and other phenomena . Losses increase as the acoustic frequency is increased (see Section 3.2). Similar processes also affect oscillations of liquid droplets, which however is not an obstacle for observation of intriguing optical phenomena even in droplets made of liquids that are 50 times more viscous than water .
Secondly, acoustic and thermal losses become especially important when strong nonlinear acoustic effects are induced . Therefore, these effects have to be taken into account to optimise the interaction of light with nonlinear acoustic waves. For example, losses can be neglected when ultrasound propagates in bulk liquid . However, the presence of viscous and thermal boundary layers  in close proximity to solid surfaces and nanostructures used to enhance the light-sound interaction no longer warrants this approximation .
In close proximity to the solid surface, the tangential particle velocity approaches zero, which gives rise to a viscous boundary layer with the thickness dv=[η/(πfρ)]−1/2, where f is the frequency of ultrasound, ρ is the density, and η is the viscosity of the liquid. Moreover, heat exchange between the liquid and the solid surface gives rise to a thermal boundary layer with the thickness dh=[κ/(πfρcp)]−1/2, where κ is the thermal conductivity of the liquid and cp is the corresponding specific heat capacity at constant pressure. One can calculate that at f=50 MHz, dv≈80 nm and dh≈30 nm, which is comparable to the dimensions of the most common plasmonic metal nanoparticles and nanoapertures.
Nevertheless, it was established that the boundary layers mostly affect the interaction of light and sound in nanoapertures  because essentially they are acoustically deep-subwavelength structures , . In the case of nanoparticles, the impact of the boundary layers is small and it does not affect the results presented in Figures 3B and 4. This is because the dimensions of the nanorod antenna used in those calculations are very small as compared with the acoustic wavelength (see supplementary information in ). However, this situation may change when an array or a colloidal solution of nanorods is used. In particular, the impact of losses will depend on the scattering strength of acoustic waves by the nanoparticles .
Finally, we note that the proposed application of acoustic nonlinearities in photonics further contributes to the ongoing effort to complement and strengthen optical nonlinearities with nonlinearities of opto-mechanical systems , ,  and spin waves . Nonlinear opto-mechanical systems open up very interesting opportunities such as second- and higher-order sideband generation , frequency comb generation , and chaos , which warrants their detailed discussion in a separate paper. Further developments in this research direction may lead to interesting fundamental physics and effects such as strong coupling between different physical entities – photons, phonons and magnons , .
We thank Antony Orth, Brant Gibson, Philipp Reineck, Leslie Yeo, Lillian Lee, Amgad Rezk (RMIT University), Mark Hutchinson (University of Adelaide), Michael Dickey (North Carolina State University, USA), Bradley Boyd (University of Canterbury) and Damien Hicks (Swinburne University) for valuable discussions. This work was supported by the Australian Research Council (ARC) through its Centre of Excellence for Nanoscale BioPhotonics (CE140100003), LIEF program (LE160100051) and Future Fellowship (FT160100357, FT180100343). This research was undertaken with the assistance of resources and services from the National Computational Infrastructure (NCI), which is supported by the Australian Government.
Pedrotti FL, Pedrotti LM, Pedrotti LS. Introduction to optics. Cambridge, UK, Cambridge University Press, 2018.
Saleh BEA, Teich MC. Fundamentals of photonics. Hoboken, NJ, USA, Wiley, 2007.
Garmire E. Nonlinear optics in daily life. Opt Express 2013;21:30533–44.
Boyd RW. Nonlinear optics. Cambridge, USA, Academic Press, 2008.
Del’Haye P, Schliesser A, Arcizet O, Wilken T, Holzwarth R, Kippenberg TJ. Optical frequency comb generation from a monolithic microresonator. Nature 2007;450:1214–7.
Wu J, Xu X, Nguyen TG, et al. RF photonics: an optical microcombs’ perspective. IEEE J Sel Top Quantum Electron 2018;24:6101020.
Magidson V, Khodjakov A. Circumventing photodamage in live-cell microscopy. Methods Cell Biol 2013;114:545–60.
Wäldchen S, Lehmann J, Klein T, van de Linde S, Sauer M. Light-induced cell damage in live-cell super-resolution microscopy. Sci Rep 2015;5:15348.
Levenson MD, Pershan PS, Mazur E, Shen YR, eds. Resonances: A volume in honor of the 70th birthday of Nicolaas Bloembergen. Singapore, World Scientific, 1990.
Hakuta K, Suzuki M, Katsuragawa M, Li JZ. Self-induced phase matching in parametric anti-stokes stimulated Raman scattering. Phys Rev Lett 1997;79:209–12.
Katsuragawa M, Liang JQ, Bakuta K. Nonlinear optics with incoherent light. Conference Digest: 2000 International Quantum Electronics Conference, Nice, France, 2000.
Picozzi A. Towards a nonequilibrium thermodynamic description of incoherent nonlinear optics. Opt Express 2007;15:9063–83.
Dylov DV. Coherent-incoherent phenomena in nonlinear optics and imaging. PhD thesis. Princeton University, USA, 2010.
Maksymov IS, Miroshnichenko AE. Active control over nanofocusing with nanorod plasmonic antennas. Opt Express 2011;19:5888–94.
Arahira S, Murai H. Wavelength conversion of incoherent light by sum-frequency generation. Opt Express 2014;22:12944–61.
Kiethe J, Heuer A, Jechow A. Second-order coherence properties of amplified spontaneous emission from a high-power tapered superluminescent diode. Laser Phys Lett 2017;14:086201.
Kurzke H, Kiethe J, Heuer A, Jechow A. Frequency doubling of incoherent light from a superluminescent diode in a periodically poled lithium niobate waveguide crystal. Laser Phys Lett 2017;14:055402.
Hamilton MF, Blackstock DT. Nonlinear acoustics. Cambridge, MA, USA, Academic Press, 1998.
Rudenko OV. Giant nonlinearities in structurally inhomogeneous media and the fundamentals of nonlinear acoustic diagnostic technique. Phys-Usp 2006;49:69–87.
Diels J-C, Rudolph W. Ultrashort laser phenomena: fundamentals, techniques, and applications on a femtosecond time scale. Cambridge, MA, USA, Academic Press, 1996.
Siegman AE. Lasers. Herndon, VA, USA, University Science Books, 1986.
Torres-Company V, Weiner AM. Optical frequency comb technology for ultra-broadband radio-frequency photonics. Laser Photonics Rev 2014;8:368–93.
Temprana E, Myslivets E, Kuo BP, et al. Overcoming Kerr-induced capacity limit in optical fiber transmission. Science 2015;348:1445–8.
Thorpe MJ, Balslev-Clausen D, Kirchner MS, Ye J. Cavity-enhanced optical frequency comb spectroscopy: application to human breath analysis. Opt Express 2008;16:2387–97.
Cherenkov AV, Kondratiev NM, Lobanov VE, Shitikov AE, Skryabin DV, Gorodetsky ML. Raman-Kerr frequency combs in microresonators with normal dispersion. Opt Express 2017;25:31148–58.
Xiong H, Si L-G, Lü X-Y, Yang X, Wu Y. Nanosecond-pulse-controlled higher-order sideband comb in a GaAs optomechanical disk resonator in the non-perturbative regime. Ann Phys 2014;349:43–54.
Reed EJ, Soljaĉić M, Joannopoulos JD. Reversed Doppler effect in photonic crystals. Phys Rev Lett 2003;91:133901.
Ta’eed V, Baker NJ, Fu L, et al. Ultrafast all-optical chalcogenide glass photonic circuits. Opt Express 2007;15:9205–21.
Yeom D-I, Mägi EC, Lamont MRE, Roelens MAF, Fu L, Eggleton BJ. Low-threshold supercontinuum generation in highly nonlinear chalcogenide nanowires. Opt Lett 2008;33:660–2.
Salem R, Foster MA, Turner AC, Geraghty DF, Lipson M, Gaeta AL. Signal regeneration using low-power four-wave mixing on silicon chip. Nat Photon 2008;2:35–8.
Aitchison J, Hutchings D, Kang J, Stegeman G, Villeneuve A. The nonlinear optical properties of AlGaAs at the half band gap. IEEE J Quantum Electron 1997;33:341–8.
Van V, Ibrahim T, Absil P, Johnson F, Grover R, Ho P-T. Optical signal processing using nonlinear semiconductor microring resonators. IEEE J Sel Top Quantum Electron 2002;8:705–13.
Bache M, Lavrinenko AV. Ultrafast nonlinear dynamics of thin gold films due to an intrinsic delayed nonlinearity. J Opt 2017;19:94004.
Klimov V. Nanoplasmonics. Singapore, Pan Stanford, 2014.
Ferrera M, Razzari L, Duchesne D, et al. Low-power continuous-wave nonlinear optics in doped silica glass integrated waveguide structures. Nat Photon 2008;2:737–40.
Alam MZ, De Leon I, Boyd RW. Large optical nonlinearity of indium tin oxide in its epsilon-near-zero region. Science 2016;352:795–7.
Kippenberg TJ, Spillane SM, Vahala KJ. Kerr-nonlinearity optical parametric oscillation in an ultrahigh-Q toroid microcavity. Phys Rev Lett 2004;93:083904.
Agha IH, Okawachi Y, Gaeta AL. Theoretical and experimental investigation of broadband cascaded four-wave mixing in high-Q microspheres. Opt Express 2009;17:16209–15.
Liang W, Savchenkov AA, Matsko AB, Ilchenko VS, Seidel D, Maleki L. Generation of near-infrared frequency combs from a MgF2 whispering gallery mode resonator. Opt Lett 2011;36:2290–2.
Herr T, Hartinger K, Riemensberger J, et al. Universal formation dynamics and noise of Kerr-frequency combs in microresonators. Nat Photon 2012;6:480–7.
Papp SB, Diddams SA. Spectral and temporal characterization of a fused-quartz-microresonator optical frequency comb. Phys Rev A 2011;84:053833.
Li J, Lee H, Chen T, Vahala KJ. Low-pump-power, low-phase-noise, and microwave to millimeter-wave repetition rate operation in microcombs. Phys Rev Lett 2012;109:233901.
Razzari L, Duchesne D, Ferrera M, et al. CMOS-compatible integrated optical hyper-parametric oscillator. Nat Photon 2010;4:41–5.
Jung H, Xiong C, Fong KY, Zhang X, Tang HX. Optical frequency comb generation from aluminum nitride microring resonator. Opt Lett 2013;38:2810–3.
Levy JS, Gondarenko A, Foster MA, Turner-Foster AC, Gaeta AL, Lipson M. CMOS-compatible multiple-wavelength oscillator for on-chip optical interconnects. Nat Photon 2010;4:37–40.
Miller SA, Okawachi Y, Ramelow S, et al. Tunable frequency combs based on dual microring resonators. Opt Express 2015;23:21527–40.
Ferdous F, Miao H, Leaird DE, et al. Spectral line-by-line pulse shaping of on-chip microresonator frequency combs. Nat Photon 2011;5:770–6.
Liu Y, Xuan Y, Xue X, et al. Investigation of mode coupling in normal-dispersion silicon nitride microresonators for Kerr frequency comb generation. Optica 2014;1:137–44.
Xue X, Xuan Y, Liu Y, et al. Mode-locked dark pulse Kerr combs in normal-dispersion microresonators. Nat Photon 2015;9: 594–600.
Griffith AG, Lau RKW, Cardenas J, et al. Silicon-chip mid-infrared frequency comb generation. Nat Commun 2015;6:6299.
Del’Haye P, Coillet A, Loh W, Beha K, Papp SB, Diddams SA. Phase steps and resonator detuning measurements in microresonator frequency combs. Nat Commun 2015;6:5668.
Almeida VR, Xu Q, Barrios CA, Lipson M. Guiding and confining Light in void nanostructure. Opt Lett 2004;29:1209–11.
Baehr-Jones T, Penkov B, Huang J, et al. Nonlinear polymer-clad silicon slot waveguide modulator with a half wave voltage of 0.25 V. Appl Phys Lett 2008;92:163303.
Merklein M, Kabakova IV, Büttner TFS, et al. Enhancing and inhibiting stimulated Brillouin scattering in photonic integrated circuits. Nat Commun 2015;6:6396.
Van Laer R, Kuyken B, Van Thourhout D, Baets R. Interaction between light and highly confined hypersound in a silicon photonic nanowire. Nat Photon 2015;9:199–203.
Palomba S, Novotny L. Nonlinear excitation of surface plasmon polaritons by four-wave mixing. Phys Rev Lett 2008;101:056802.
Grady NK, Knight MW, Bardhan R, Halas NJ. Optically-driven collapse of a plasmonic nanogap self-monitored by optical frequency mixing. Nano Lett 2010;10:1522–8.
Ko KD, Kumar A, Fung KH, et al. Nonlinear optical response from arrays of Au bowtie nanoantennas. Nano Lett 2011;11:61–85.
Schumacher T, Kratzer K, Molnar D, Hentschel M, Giessen H, Lippitz M. Nanoantenna-enhanced ultrafast nonlinear spectroscopy of a single gold nanoparticle. Nat Commun 2011;2:333.
Harutyunyan H, Volpe G, Quidant R, Novotny L. Enhancing the nonlinear optical response using multifrequency gold-nanowire antennas. Phys Rev Lett 2012;108:217403.
Slablab A, Xuan LL, Zielinski M, et al. Second-harmonic generation from coupled plasmon modes in a single dimer of gold nanospheres. Opt Express 2012;20:220–7.
Abb M, Wang Y, Albella P, de Groot CH, Aizpurua J, Muskens OL. Interference, coupling, and nonlinear control of high-order modes in single asymmetric nanoantennas. ACS Nano 2012;6:6462–70.
Navarro-Cia M, Maier SA. Broad-band near-infrared plasmonic nanoantennas for higher harmonic generation. ACS Nano 2012;6:3537–44.
Hentschel M, Utikal T, Giessen H, Lippitz M. Quantitative modeling of the third harmonic emission spectrum of plasmonic nanoantennas. Nano Lett 2012;12:3778–82.
Maksymov IS, Miroshnichenko AE, Kivshar YS. Cascaded four-wave mixing in tapered plasmonic nanoantenna. Opt Lett 2013;38:79–81.
Drachev VP, Kildishev AV, Borneman JD, et al. Engineered nonlinear materials using gold nanoantenna array. Sci Rep 2018;8:780.
Chen P-Y, Argyropoulos C, D’Aguanno G, Alú A. Enhanced second-harmonic generation by metasurface nanomixerand nanocavity. ACS Photon 2015;2:1000–6.
Smirnova D, Kivshar YS. Multipolar nonlinear nanophotonics. Optica 2016;11:1241–55.
Wang C, Li Z, Kim M-H, et al. Metasurface-assisted phase-matching-free second harmonic generation in lithium niobate waveguides. Nat Commun 2017;8:2098.
Chen S, Rahman M, Fai K, et al. Third harmonic generation enhanced by multipolar interference in complementary silicon metasurfaces. ACS Photon 2018;5:1671–5.
Lau RKW, Lamont MRE, Okawachi Y, Gaeta AL. Effects of multiphoton absorption on parametric comb generation in silicon microresonators. Opt Lett 2015;40:2778–81.
Liu M, Wang L, Sun Q, et al. Influences of multiphoton absorption and free-carrier effects on frequency-comb generation in normal dispersion silicon microresonators. Photon Res 2018;6:238–43.
Stern B, Ji X, Okawachi Y, Gaeta AL, Lipson M. Battery-operated integrated frequency comb generator. Nature 2018;562:401–5.
Rogov AS, Narimanov EE. Nonlinear optics at low powers: alternative mechanism of on-chip optical frequency comb generation. Phys Rev A 2016;94:063832.
Li J, Shen S, Qu Y, Zhang D, Wu Y. Generating orthogonally polarized dual frequency combs with slow megahertz repetition rates by a low-nanowatt-level pump. Phys Rev A 2018;98:023848.
Kuznetsov VP. Equations of nonlinear acoustics. Akusticheskij Zhurnal 1970;16:548–53.
Perrin B. Some aspects of the comparison between optics and nonlinear acoustics (Quelques aspects de la comparison entre optique et acoustique non lineaires). J Phys Colloques 1979;11:C8–216.
de Gennes PG, Brochard-Wyart F, Quéré D. Capillarity and wetting phenomena: drops, bubbles, pearls, waves. Berlin, Germany, Springer, 2004.
van Capel PJS, Péronne E, Dijkhuis JI. Nonlinear ultrafast acoustics at the nano scale. Ultrasonics 2015;56:36–51.
Fabelinskii IL. Molecular scattering of light. Berlin, Germany, Springer, 1968.
Störkel U, Vodopyanov KL, Grill W. GHz ultrasound wave packets in water generated by an Er laser. J Phys D: Appl Phys 1998;31:2258–63.
Rossignol C, Chigarev N, Ducousso M, et al. In vitro picosecond ultrasonics in a single cell. Appl Phys Lett 2008;93:123901.
Wright OB, Perrin B, Matsuda O, Gusev VE. Optical excitation and detection of picosecond acoustic pulses in liquid mercury. Phys Rev B 2008;78:024303.
Zijlstra A, Fernandez-Rivas D, Han Gardeniers JGE, Versluis M, Lohse D. Enhancing acoustic cavitation using artificial crevice bubbles. Ultrasonics 2015;56:512–23.
Rakich PT, Reinke C, Camacho R, Davids P, Wang Z. Giant enhancement of stimulated Brillouin scattering in the subwavelength limit. Phys Rev X 2012;2:011008.
Baker C, Hease W, Nguyen D-T, et al. Photoelastic coupling in gallium arsenide optomechanical disk resonators. Opt Express 2014;22:14072–86.
Babchin AJ, Gur Y, Lin IJ. Repulsive interface forces in overlapping electric double layers in electrolyte solutions. Adv Colloid Interface Sci 1978;9:105–41.
Gojani AB, Bejtullahu R, Obayashi S. On two optomechanical effects of laser-induced electrostriction in dielectric liquids. Jpn J Appl Phys 2014;53:092703.
Guan B-O, Jin L, Cheng L, Liang Y. Acoustic and ultrasonic detection with radio-frequency encoded fiber laser sensors. IEEE J Sel Top Quantum Electron 2017;23:5601712.
Wilkens V. Characterization of an optical multilayer hydrophone with constant frequency response in the range from 1 to 75 MHz. J Acoust Soc Am 2003;113:1431–8.
Li G, Guo Z, Chen S-L. Miniature all-optical probe for large synthetic aperture photoacoustic-ultrasound imaging. Opt Express 2017;25:25023–35.
Xia W, West SJ, Finlay MC, et al. Three-dimensional ultrasonic needle tip tracking with a fiber-optic ultrasound receiver. J Vis Exp 2018;138:e57207.
Kessler BV. Method and apparatus for optically detecting acoustic disturbances. US Patent 3,474,253. 1969.
Morris P, Hurrell A, Shaw A, Zhang E, Beard P. A Fabry-Perot fiber-optic ultrasonic hydrophone for the simultaneous measurement of temperature and acoustic pressure. J Acoust Soc Am 2009;125:3611–22.
Chao C-Y, Ashkenazi S, Huang S-W, O’Donnell M, Jay Guo L. High-frequency ultrasound sensors using polymer microring resonators. IEEE Trans Ultrason Ferroelectr Freq Control 2007;54:957–65.
Li H, Dong B, Zhang Z, Zhang HF, Sun C. A transparent broadband ultrasonic detector based on an optical micro-ring resonator for photoacoustic microscopy. Sci Rep 2014;4:4496.
Yakovlev VV, Dickson W, Murphy A, et al. Ultrasensitive non-resonant detection of ultrasound with plasmonic metamaterials. Adv Mater 2013;25:2351–6.
Maksymov IS, Greentree AD. Synthesis of discrete phase-coherent optical spectra from nonlinear ultrasound. Opt Express 2017;25:7496–506.
Agio M, Alù A. Optical antennas. Cambridge, England, Cambridge University Press, 2013.
Ruijgrok PV, Zijlstra P, Tchebotareva AL, Orrit M. Damping of acoustic vibrations of single gold nanoparticles optically trapped in water. Nano Lett 2012;12:1063–9.
O’Brien K, Lanzillotti-Kimura ND, Rho J, Suchowski H, Yin X, Zhang X. Ultrafast acousto-plasmonic control and sensing in complex nanostructures. Nat Commun 2014;5:4042.
Chang WS, Wen F, Chakraborty D, et al. Tuning the acoustic frequency of a gold nanodisk through its adhesion layer. Nat Commun 2015;6:7022.
Wu J, Xiang D, Gordon R. Characterizing gold nanorods in aqueous solution by acoustic vibrations probed with four-wave mixing. Opt Express 2016;24:12458–65.
Loh W, Yegnanarayanan S, Ram RJ, Juodawlkis PW. A nonlinear optoelectronic filter for electronic signal processing. Sci Rep 2014;4:3613.
de Chatellus HG, Cortés LR, Azaña J. Optical real-time Fourier transformation with kilohertz resolutions. Optica 2016;3:1–8.
Anand S, Eryürek M, Karadag Y, et al. Observation of whispering gallery modes in elastic light scattering from microdroplets optically trapped in a microfluidic channel. J Opt Soc Am B 2016;33:1349–54.
Giorgini A, Avino S, Malara P, De Natale P, Gagliardi G. Opto-mechanical oscillator in a nanoliter droplet. Opt Lett 2018;43:3473–6.
Bar-David D, Maayani S, Martin LL, Carmon T. Cavity optofluidics: a μdroplet’s whispering-gallery mode makes a μvortex. Opt Express 2018;26:19115–22.
Maksymov IS, Greentree AD. Acoustically tunable optical transmission through a subwavelength hole with a bubble. Phys Rev A 2017;95:033811.
Maksymov IS, Greentree AD. Dynamically reconfigurable plasmon resonances enabled by capillary oscillations of liquid-metal nanodroplets. Phys Rev A 2017;96:043829.
Schulman RD, Dalnoki-Veress K. Droplets capped with an elastic film can be round, elliptical, or nearly square. Phys Rev Lett 2018;121:248004.
Mitsui T. Photoacoustic excitation of ripplon on liquid droplet adhered to tip of optical fiber. Jpn J Appl Phys 2004; 43:8345–6.
Kaminski S, Martin LL, Maayani S, Carmon T. Ripplon laser through stimulated emission mediated by water waves. Nat Photon 2016;10:758–61.
Laude V, Belkhir A, Alabaid AF, et al. Extraordinary nonlinear transmission modulation in a doubly resonant acousto-optical structure. Optica 2017;4:1245–50.
Carmon T, Vahala KJ. Modal spectroscopy of optoexcited vibrations of a micron-scale on-chip resonator at greater than 1 GHz frequency. Phys Rev Lett 2007;98:123901.
Li J, Zhu S-Y, Agarwal GS. Magnon-photon-phonon entanglement in cavity magnomechanics. arXiv: 1807.07158v3.
Maksymov IS. Perspective: strong microwave photon-magnon coupling in multiresonant dielectric antennas. J Appl Phys 2018;124:150901.
Friend J, Yeo L. Using laser Doppler vibrometry to measure capillary surface waves on fluid-fluid interfaces. Biomicrofluidics 2010;4:026501.
Blamey J, Yeo LY, Friend JR. Microscale capillary wave turbulence excited by high frequency vibration. Langmuir 2013;29:3835–45.
Argo IV TF, Wilson PS, Palan V. Measurement of the resonance frequency of single bubbles using a laser Doppler vibrometer. J Acoust Soc Am 2006;123:EL121–5.
Yoshida K, Yoshikawa T, Koyma D, Nakamura K, Watanabe Y. Experimental measurement of microbubble oscillation by using laser Doppler vibrometer. AIP Conf Proc 2012;1433:304.
Kuo KA, Hunt HEM. Vibrations of bubbles and balloons. Acoust Australia 2012;40:183–7.
Xu L, Rahmani M, Smirnova D, et al. Highly-efficient longitudinal second-harmonic generation from doubly-resonant AlGaAs nanoantennas. Photonics 2018;5:29.
Cavalieri F, Micheli L, Kaliappan S, et al. Antimicrobial and biosensing ultrasound-responsive lysozyme-shelled microbubbles. ACS Appl Mater Interfaces 2013;5:464–71.
Lin CA, Chuang WK, Huang ZY, et al. Rapid transformation of protein-caged nanomaterials into microbubbles as bimodal imaging agents. ACS Nano 2012;6:5111–21.
Lin Y, Liu Y, Genzer J, Dickey MD. Shape-transformable liquid metal nanoparticles in aqueous solution. Chem Sci 2017;8:3832–7.
Daeneke T, Khoshmanesh K, Mahmood N, et al. Liquid metals: fundamentals and applications in chemistry. Chem Soc Rev 2018;47:4073–111.
Wu PC, Kim TH, Brown AS, Losurdo M, Bruno G, Everitt HO. Real-time plasmon resonance tuning of liquid Ga nanoparticles by in situ spectroscopic ellipsometry. Appl Phys Lett 2007;90:103119.
Sanz JM, Ortiz D, Alcaraz de la Osa R, et al. UV plasmonic behavior of various metal nanoparticles in the near- and far-field regimes: geometry and substrate effects. J Phys Chem C 2013;117:19606–15.
McMahon JM, Schatz GC, Gray SK. Plasmonics in the ultraviolet with the poor metals Al, Ga, In, Sn, Ti, Pb, and Bi. Phys Chem Chem Phys 2013;15:5415–23.
Naik GV, Shalaev VM, Boltasseva A. Alternative plasmonic materials: beyond gold and silver. Adv Mater 2013;25: 3264–94.
Yang Y, Akozbek N, Kim T-H, et al. Ultraviolet-visible plasmonic properties of gallium nanoparticles investigated by variable-angle spectroscopic and Mueller matrix ellipsometry. ACS Photon 2014;1:582–9.
Knight MW, Coenen T, Yang Y, et al. Gallium plasmonics: deep subwavelength spectroscopic imaging of single and interacting gallium nanoparticles. ACS Nano 2015;9:2049–60.
Zhang T, Wang Y, Appusamy K, et al. Gallium platinum alloys – a new material system for UV plasmonics. Opt Mater Express 2017;7:2880–7.
Catalán-Gómez S, Redondo-Cubero A, Palomares FJ, Nucciarelli F, Pau JL. Tunable plasmonic resonance of gallium nanoparticles by thermal oxidation at low temperatures. Nanotechnology 2017;28:405705.
Nucciarelli F, Bravo I, Catalan-Gomez S, Vázquez L, Lorenzo E, Pau JL. High ultraviolet absorption in colloidal gallium nanoparticles prepared from thermal evaporation. Nanomaterials 2017;7:172.
Gutiérrez Y, Alcaraz de la Osa R, Ortiz D, Saiz JM, González F, Moreno F. Plasmonics in the ultraviolet with aluminum, gallium, magnesium and rhodium. Appl Sci 2018;8:64.
Morales D, Stoute NA, Yu Z, Aspnes DE, Dickey MD. Liquid gallium and the eutectic gallium indium (EGaIn) alloy: dielectric functions from 1.24 to 3.1 eV by electrochemical reduction of surface oxides. Appl Phys Lett 2016;109:091905.
Morley NB, Burris J, Cadwallader LC, Nornberg MD. GaInSn usage in the research laboratory. Rev Sci Instrum 2008;79:056107.
Dickey MD, Chiechi RC, Larsen RJ, Weiss EA, Weitz DA, Whitesides GM. Eutectic gallium-indium (EGaIn): a liquid metal alloy for the formation of stable structures in microchannels at room temperature. Adv Funct Mater 2008;18:1097–104.
Reineck P, Lin Y, Gibson BC, Dickey MD, Greentree AD, Maksymov IS. UV plasmonic properties of colloidal liquid-metal eutectic gallium-indium alloy nanoparticles. arXiv:1810.06163.
Gibbs HM. Optical bistability: controlling light with light. Cambridge, MA, USA, Academic Press, 1985.
Ikeda K, Akimoto O. Instability leading to periodic and chaotic self-pulsations in a bistable optical cavity. Phys Rev Lett 1982;48:617–20.
Wang QX, Blake JR. Non-spherical bubble dynamics in a compressible liquid. Part 1. Travelling acoustic wave. J Fluid Mech 2010;659:191–224.
Wang QX, Blake JR. Non-spherical bubble dynamics in a compressible liquid. Part 2. Acoustic standing wave. J Fluid Mech 2011;679:559–81.
Fong SW, Klaseboer E, Turangan CK, Khoo BC, Hung KC. Numerical analysis of a gas bubble near bio-materials in an ultrasound field. Ultrasound Med Biol 2006;32: 925–42.
Boyd B, Becker S. Numerical modelling of an acoustically-driven bubble collapse near a solid boundary. Fluid Dyn Res 2018;50:065506.
Supponen O, Obreschkow D, Tinguely M, Kobel P. Scaling laws for jets of single cavitation bubbles. J Fluid Mech 2016;802:263–93.
Blake JR, Leppinen DM, Wang Q. Cavitation and bubble dynamics: the Kelvin impulse and its applications. Interface Focus 2015;5:20150017.
Wang QX, Yeo KS, Khoo BC, Lam KY. Nonlinear interaction between gas bubble and free surface. Comput Fluids 1996;25:607–28.
Joannopoulos JD, Johnson SG, Winn JN, Meade RD. Photonic crystals: molding the flow of light. Princeton, NJ, USA, Princeton University Press, 2008.
Smith CLC, Bog U, Tomljenovic-Hanic S, et al. Reconfigurable microfluidic photonic crystal slab cavities. Opt Express 2008;16:15887–96.
Khelif A, Adibi A. Phononic crystals: fundamentals and applications. New York, USA, Springer, 2016.
Webb CE, Jones JDC. Handbook of laser technology and applications. Volume III: applications. Bristol, UK, Institute of Physics Publishing, 2004.
Masters BR, So PTC. Handbook of biomedical nonlinear optical microscopy. Oxford, UK, Oxford University Press, 2008.
Blomley MJK, Cooke JC, Unger EC, Monaghan MJ, Cosgrove DO. Microbubble contrast agents: a new era in ultrasound. Br Med J 2001;322:1222–5.
Qin S, Caskey CF, Ferrara KW. Ultrasound contrast microbubbles in imaging and therapy: physical principles and engineering. Phys Med Biol 2009;54:R27–57.
Jung EE, Erickson D. Continuous operation of a hybrid solid-liquid state reconfigurable photonic system without resupply of liquids. Lab Chip 2012;12:2575–9.
Phillips CL, Jankowski E, Krishnatreya BJ, et al. Digital colloids: reconfigurable clusters as high information density elements. Soft Matter 2014;10:7468–79.
Altkorn R, Koev I, Van Duyne RP, Litorja M. Low-loss liquid-core optical fiber for low-refractive-index liquids: fabrication, characterization, and application in Raman spectroscopy. Appl Opt 1997;36:8992–8.
Couny F, Benabid F, Roberts PJ, Light PS, Raymer MG. Generation and photonic guidance of multi-octave optical-frequency combs. Science 2007;318:1118–21.
Benoît A, Beaudou B, Alharbi M, et al. Over-five octaves wide Raman combs in high-power picosecond-laser pumped H2-filled inhibited coupling Kagome fiber. Opt Express 2015;23:14002–9.
Stepanov VI, Raikher YL. Dynamic birefringence in magnetic fluids with allowance for mechanical and magnetic degrees of freedom of the particles. J Magn Magn Mater 2002;252:180–2.
Shulyma SI, Tanygin BM, Kovalenko VF, Petrychuk MV. Magneto-optical extinction trend inversion in ferrofluids. J Magn Magn Mater 2016;416:141–9.
Morse PM, Uno Ingard K. Theoretical acoustics. New York, USA, McGraw-Hill, 1968.
Ward GP, Lovelock RK, Murray ARJ, Hibbins AP, Sambles JR, Smith JD. Boundary-layer effects on acoustic transmission through narrow slit cavities. Phys Rev Lett 2015;115: 044302.
Homentcovschi D, Miles RN, Tan L. Influence of viscosity on the diffraction of sound by a periodic array of screens. J Acoust Soc Am 2005;117:2761–71.
Baudoin M, Thomas J-L, Coulouvrat F, Chaneac C. Scattering of ultrasonic shock waves in suspensions of silica nanoparticles. J Acoust Soc Am 2011;129:1209–20.
Xiong H, Si L-G, Zheng A-S, Yang X, Wu Y. Higher-order sidebands in optomechanically induced transparency. Phys Rev A 2012;86:013815.
Cao C, Mi S-C, Wang T-J, Zhang R, Wang C. Optical high-order sideband comb generation in a photonic molecule optomechanical system. IEEE J Quantum Electron 2016;52:7000205.
Vitko VV, Nikitin AA, Ustinov AB, Kalinikos BA. Microwave bistability in active ring resonators with dual spin-wave and optical nonlinearities. IEEE Magn Lett 2018;9:3506304.
Sun Y, Sukhorukov AA. Chaotic oscillations of coupled nanobeam cavities with tailored optomechanical potentials. Opt Lett 2014;39:3543–6.
Childress L, Schmidt MP, Kashkanova AD, et al. Cavity optomechanics in a levitated helium drop. Phys Rev A 2017;98:063842.