The dimer under study is a dielectric structure formed of two identical sub-units. Dimer interactions with electromagnetic waves are widely studied in connection with electromagnetic properties of complex systems. The dimer in the form of two closely spaced grapes is also a subject of high public interest because its gap region sparks curiously in a household microwave oven. A recent paper presented the first scientific interpretation of this long-standing puzzle. It attributed the sparking phenomenon to an electromagnetic hotspot in the gap created by electromagnetic resonances in the dimer. This study has opened up a fertile ground for further research as well as motivating the current study on an independent mechanism. Our simulation and experimental results consistently point to an electrical origin for the sparks. The triggering mechanism is a two to three orders-of-magnitude buildup of localized electric field in the narrow gap, which results from the mutual enhancement of polarization charges on opposite sides of the gap. Consequently, sparks are observed at the frequency of 27 MHz, at which no electromagnetic resonance in the dimer is possible. Results also indicate a broad frequency range of the electrical mechanism, which persists even when the dimer is in strong electromagnetic resonance with the first few higher-order modes. These quantitative characterizations of basic dimer properties, in particular the broad frequency range of the polarization-charge enhancement effect, may be helpful for the understanding of collective behavior of multi-particle systems under electromagnetic radiations.

Dimer interactions with electromagnetic waves offer a research tool to obtain microscopic insights into macroscopic matter behavior. They also generate a popularly known household phenomenon in the form of fiery sparks between two closely spaced grapes in a microwave oven. However, the cause has not been scientifically examined until recently. In a detailed study on aqueous dimers,1 Khattak, Bianucci, and Slepkov concluded that the merging of resonant modes in two closely spaced watery spheres results in a hot electromagnetic spot in their gap, which causes the sparks. This work and the proposed applications have been widely reported2–5 because of their novelty and public interest, in particular, the claimed demonstration of optical writing with a resolution better than λ0/80.

In Ref. 1, the authors put a long-standing household puzzle on a rich academic basis. Their work has also motivated the current study on an independent cause of a much different nature. Both theory and experiment indicate that with the wave polarized along the dimer axis (same model as in Ref. 1), mutual enhancement of polarization charges on both sides of the narrow gap can result in an electric field hundreds of times greater than that of the wave, hence triggering the sparks through air arcing. This is evidenced by (1) the predicted and observed sparks in an essentially magnetic-field free environment at 27 MHz, (2) the negligibly small magnetic field in the dimer gap even when the dimer body is in strong electromagnetic resonance with a 2.45-GHz microwave, and (3) a video display of the attraction by electric force between the two spheres of the dimer, rather than the repulsion expected from the radiation pressure of an electromagnetic hotspot.

It is well recognized that bound molecular charges in a dielectric medium will be slightly displaced by an electric field to form polarization charges. The simplest example is a uniform dielectric sphere immersed in a static, uniform, external electric field (Eext), for which the exact analytical solutions are well-documented (Ref. 6, Sec. 4.4). The main feature is the induced polarization charges (σpol) on the spherical surface [see Fig. 1(a)] given by

(1)

where ε and ε0 are the permittivity of the sphere and free space, respectively. θ is the angle between Eext and the observation direction.

FIG. 1.

Electric field patterns of a dielectric sphere immersed in a static, uniform, external electric field, Eext. (a) Analytically calculated electric field lines for a sphere with ε/ε0 = 4. Polarization charges are induced on the surface. (b) Simulated electric field E(x) for a water sphere with a much larger ε/ε0 of 80. In (b), the color code gives the normalized E(x)/Eext in logarithmic scale: 10 × log10[E(x)/Eext].

FIG. 1.

Electric field patterns of a dielectric sphere immersed in a static, uniform, external electric field, Eext. (a) Analytically calculated electric field lines for a sphere with ε/ε0 = 4. Polarization charges are induced on the surface. (b) Simulated electric field E(x) for a water sphere with a much larger ε/ε0 of 80. In (b), the color code gives the normalized E(x)/Eext in logarithmic scale: 10 × log10[E(x)/Eext].

Close modal

Surface polarization charges partially cancel the Eext in the interior of the sphere, while still keeping the total interior field (Ein) uniform. The larger the dielectric constant, the greater the cancelation. For example, Ein/Eext = 0.5 for ε/ε0 = 4 and Ein/Eext = 0.037 for ε/ε0 = 80. Two complementary methods are used for graphic illustrations. Figure 1(a) displays a field-line plot for ε/ε0 = 4 to show field directions and discontinuities (at σpol), while Fig. 1(b) displays a simulated color plot for ε/ε0 = 80 to show large field variations (see  Appendix A, Part 1, for the software used in all simulations). Here and in all subsequent simulations, air is assumed to be free space.

Consider next a dielectric sphere (of radius R) exposed to an electromagnetic wave with free-space wavelength λfree. The wavelength in the dielectric (λd) is λfree divided by the index of reflection [= Re(ε/ε0)1/2 (Ref. 6, p. 314)]. If λd ≫ R, the region in and around the sphere falls in the near zone. The induced fields are, thus, dominated by a quasi-static electrical component, which oscillates harmonically in time with a spatial profile well approximated by that of the static equations (Ref. 6, Sec. 10.1), i.e., the fields are essentially the same as in Fig. 1.

Because water has a very large ε/ε0 at and below microwave frequencies, the induced ±σpol cancels most of Eext in the sphere and approximately doubles it immediately out of the right and left ends of the sphere [Fig. 1(b)]. A weakly loaded, 1-kW microwave oven has a peak electric field (Eext) of 300–400 V/cm,7 while the air breakdown strength is about 3 × 104 V/cm at 1 atm.8 So a mere doubling of Eext in the microwave oven is far from enough to trigger a spark (as pointed out in Ref. 1).

The situation is markedly different in the gap region between two water spheres aligned side-by-side along Eext [Fig. 2(a)]. First, the induced polarization charges are no longer exactly on the surface. They penetrate slightly into the spheres to form a volume charge density (ρpol). However, the penetration is so shallow that ρpol still effectively resides on the surface as in Eq. (1). Second and more importantly, ρpol on opposite sides of the gap is of opposite signs (denoted below by ±ρpol). Facing each other, ±ρpol is mutually enhanced, i.e., more molecules are polarized on the +ρpol sphere by the action of −ρpol on the opposite sphere and vice versa. A narrower gap leads to greater enhancement, consequently a greater gap electric field. Such an effect is well recognized. The relevant question here is whether it can drive up the gap electric field to the air breakdown strength, hence triggering the sparks in a household microwave oven.

FIG. 2.

Simulation model and results at 27 MHz under Eext (instantaneous) = 400 V/cm. (a) Model of a dimer composed of two R = 7 mm water spheres separated by a variable gap width d. (b) The simulation model and x-y plane E(x) of a capacitor loaded with a dimer (with d = 0.5 mm) aligned along Eext. The x-y plane cuts through the centers of the two spheres. (c) Close-up view of the dimer in (b). In (b) and (c), the color code gives the normalized E(x)/Eext in logarithmic scale 10 × log10[E(x)/Eext]. (d) Eaxis profiles along the x-axis for several values of d under the model in (b). (e) Egap as a function of d for d =0.005–3 mm under the model in (b). (f) For comparison, the dimer is now aligned perpendicular to Eext. Other parameters are unchanged. (g) Close-up view of the dimer in (f) showing no enhancement of the gap electric field.

FIG. 2.

Simulation model and results at 27 MHz under Eext (instantaneous) = 400 V/cm. (a) Model of a dimer composed of two R = 7 mm water spheres separated by a variable gap width d. (b) The simulation model and x-y plane E(x) of a capacitor loaded with a dimer (with d = 0.5 mm) aligned along Eext. The x-y plane cuts through the centers of the two spheres. (c) Close-up view of the dimer in (b). In (b) and (c), the color code gives the normalized E(x)/Eext in logarithmic scale 10 × log10[E(x)/Eext]. (d) Eaxis profiles along the x-axis for several values of d under the model in (b). (e) Egap as a function of d for d =0.005–3 mm under the model in (b). (f) For comparison, the dimer is now aligned perpendicular to Eext. Other parameters are unchanged. (g) Close-up view of the dimer in (f) showing no enhancement of the gap electric field.

Close modal

We first address this question at 27 MHz, which belongs to an ISM (industrial, scientific, and medical) band. The dimer under study is formed of two identical water spheres with a variable gap width d [Fig. 2(a)]. Here and subsequently, the radius (R) of each sphere is fixed at 7 mm, so λd (≈ 124 cm)≫R at 27 MHz and the fields are quasi-static with a negligibly small magnetic component in the near zone.

In the 27 MHz simulation, a voltage is applied to a capacitor formed of two circular plates, 10 cm in radius and 4 cm apart, with their axis coinciding with the x-axis [Fig. 2(b)]. The on-axis electric field, Eext, is nearly uniform (in the ±x-direction) between the plates. Two water spheres, each with R (radius) = 7 mm and separated by d (gap width) = 0.5 mm, are aligned along Eext [Fig. 2(b)]. The complex permittivity of water at 27 MHz is assumed to be ε = (78.4 + 0.1i)ε0.9 

Figure 2(b) shows the simulated x-y plane electric field strength, E(x), of a capacitor loaded with a dimer (d =0.5 mm) aligned along Eext, with an instantaneous field strength of 400 V/cm, approximately the maximum field in an 1-kW microwave oven.7 The close-up view of the dimer [Fig. 2(c)] features a strong ρpol-enhanced electric field in the gap region.

Simulations have also been performed for the dimer in the same instantaneous Eext for several other gap widths. Figure 2(d) displays the on-axis electric field Eaxis along the x-axis, each for a fixed d. In all profiles, Eaxis peaks sharply in the gap region where it is much greater than Eext and more so for a smaller d.

We further define Egap = Eaxis (x =0) to quantify the gap electric field strength. Figure 2(e) plots Egap as a function of d for d =0.005–3 mm. With air modeled as free space, air breakdown cannot be simulated. However, it can be seen that Egap reaches the air breakdown strength (∼3 × 104 V/cm) at d 0.13 mm. With d <0.13 mm, sparking is expected, and the required energy is localized within the gap region [Figs. 2(c) and 2(d)].

As a point of interest, when the two spheres are in contact (d =0), there is still a gap region around the contact point. Furthermore, bound to their respective molecules, polarization charges are not free to move across the contact point. Hence, the transition from separation to point contact is a smooth one, and the electric field near the contact point is accurately given by the d→0 data in Fig. 2(e).

For the purpose of comparison, in Fig. 2(f), we rotate the dimer in Fig. 2(b) by 90° to make the dimer axis aligned perpendicular to Eext. Figure 2(g) gives the close-up view of the dimer in Fig. 2(f). It shows no enhancement of the gap electric field. This is expected because |ρpol| is now at a minimum near the gap [θ = ±π/2 in Eq. (1)].

This is borne out in experiments at 27 MHz (see  Appendix A, Part 4, for the experimental setup). Figure 3(a) shows a capacitor [same dimensions as in Fig. 2(b)] loaded with two R =7 mm hydrogel spheres aligned along Eext. The hydrogel spheres are fixed in a sample holder (see  Appendix A, Part 2, for the sample preparation). The capacitor in Fig. 3 is plotted vertically, so that its Eext is along the x-axis as in Fig. 2. In actual experiments, the capacitor is horizontal, while the sample holder can be rotated by 90°. An infrared camera is used to obtain thermal images during the experiment (see  Appendix A, Part 3).

FIG. 3.

Experimental configuration and in situ thermal images of the dimer in a 27 MHz capacitor. (a) Two R= 7 mm hydrogel spheres are aligned parallel to Eext of the capacitor. (b) Thermal images of the dimer in (a) for three gap widths. (c) Two R= 7 mm hydrogel spheres are aligned perpendicular to Eext. (d) Thermal images of the dimer in (c) for the three gap widths. In (b) and (d), Eext (peak) = 400 V/cm, and the exposure time is 35 s.

FIG. 3.

Experimental configuration and in situ thermal images of the dimer in a 27 MHz capacitor. (a) Two R= 7 mm hydrogel spheres are aligned parallel to Eext of the capacitor. (b) Thermal images of the dimer in (a) for three gap widths. (c) Two R= 7 mm hydrogel spheres are aligned perpendicular to Eext. (d) Thermal images of the dimer in (c) for the three gap widths. In (b) and (d), Eext (peak) = 400 V/cm, and the exposure time is 35 s.

Close modal

The 27 MHz arcing typically started as soon as Eext = 400 V/cm was applied (the field used in theory), followed shortly by spark ignitions (see supplementary videos). However, the gap surface temperature did not rise as rapidly as the optically bright sparks would suggest. At 27 MHz, the conventional AC circuit theory applies to our circuit dimensions. The impedance between each capacitor plate and the dimer [see Fig. 3(a)] is capacitive (∼1/ω), hence higher at lower frequencies. A relatively high circuit impedance at 27 MHz may have resulted in the slow heating observed.

Figure 3(b) displays the thermal images of the dimer surface for three gap widths, taken in situ by an infrared camera after a 35-sec exposure to Eext (peak) = 400 V/cm. The initial dimer temperature is 23.7 °C. For d =0.7 mm, the gap is too far apart to trigger an arc. For d =0.2 mm, the gap starts to arc/spark with Tmax (spatial maximum) = 41.3 °C on the gap surfaces. For d =0, the gap arcs/sparks strongly with Tmax = 72.1 °C.

Assuming an air breakdown strength of 3 × 104 V/cm,8 the simulated threshold (i.e., longest) d for arcing is 0.13 mm [Fig. 2(e)], while the measured value is 0.2 mm [Fig. 3(b)], indicating a lower actual air breakdown strength. There is a ±10% uncertainty in the capacitor voltage measurement. On the other hand, the AC air breakdown strength is significantly lower than the DC value (by as much as 20% at 30 kHz as shown in Ref. 10, which also reports a weaker dependence on humidity). Considering these factors, the measured threshold d for arcing (∼0.2 mm) is in reasonable agreement with the prediction (0.13 mm).

To further confirm that the gap electric field originates from the ± ρpol buildup, we rotate the dimer by 90° with respect to the Eext [Fig. 3(c)]. As expected, no arcing is observed [Fig. 3(d)] since |ρpol| is now at a minimum near the gap [θ = ±π/2 in Eq. (1)].

To conclude, the electrical origin for gap sparks has, thus, been established in an essentially magnetic-field free environment at 27 MHz.

It should be pointed out that the effect under study is not specific to aqueous dimers. E-field intensification in dielectric or conducting dimers has been commonly studied to explain rapid or non-uniform processing of materials by microwaves.11–15 The treated systems are rather complex, and multiple effects are involved in addition to electrical polarization. The current model of an isolated dimer allows a detailed quantitative study on various effects such as the polarization-charge enhancement.

An important issue to address is the frequency range for the validity of the electrical mechanism, especially when the dimer is in strong electromagnetic resonances with the incident wave.

As an electromagnetic wave enters a dielectric object, multiple reflections off the inside walls can constructively superpose into a resonant standing wave at certain discrete frequencies. In 1908, Mie developed the analytical theory for a plane wave incident onto a dielectric sphere.16 A dielectric sphere of radius R is known to exhibit Mie resonances in an electromagnetic wave if λd ≤ 2 R. For example, at 2.45 GHz, water has a dielectric constant of ε/ε0= 77.5 + 10i.17 Thus, λfree = 12.25 cm is shortened to λd 1.39 cm in water. By documented formula,18 at 2.45 GHz, the lowest-order water sphere resonance occurs at R 6.9 mm and the second resonance at R ≈ 9.8 mm, each with a rather broad resonant width due to the large loss tangent [see, for example, Fig. 5(b) of Ref. 1].

Instead of varying R to bring out the electromagnetic resonances, we may equivalently vary the incident wave frequency (f) at the fixed sphere radius of R =7 mm. In this case, the first resonance occurs at ∼2.4 GHz and the second resonance at ∼3.4 GHz. Assume that a plane wave with electric field linearly polarized in the x-direction is propagating in the z-direction. We denote the electric field amplitude of the incident plane wave by Eext. The magnetic field Hext is in phase with the electric field with an amplitude given by Hext = Eext/377 in SI unit (377 being the impedance of free space in Ohm).

Let the wave be incident on a dimer of two R =7 mm water spheres on the x-y plane with its axis aligned along Eext (the most favorable orientation for polarization charge formation across the gap). As in the 27 MHz simulation, we assume Eext = 400 V/cm for all frequencies simulations. However, the dimer has a frequency-dependent complex permittivity given by the Debye function,19 

(2)

for exp(−ωt) time dependence, where τ is the temperature-dependent relaxation time. Value of ε(ω = ), ε(ω = 0), and τ can be obtained by a fitting procedure with measured values.20 For water at 25 °C, ε(ω = ) = 5.2ε0, ε(ω = 0) = 78.36ε0, and τ = 8.27 × 10−12 (Ref. 20, Table III).

For water, we have used ε/ε0 = 78.4 + 0.1i for 27 MHz (Refs. 9 and ε/ε0 = 77.5 + 10i for 2.45 GHz.17 Question arises as to whether Eq. (2) can yield similar ε/ε0 at these two frequencies. Simple substitutions give ε/ε0 = 78.36 + 0.10i for 27 MHz and ε/ε0 = 77.19 + 9.17i for 2.45 GHz, both being very close to what we used in single-frequency cases.

When the quasi-static condition (λd ≫ R) is not satisfied at high frequencies, fields at different points will not oscillate in phase. However, the peak value (amplitude) is of primary interest. Hence, we denote the on-axis electric and magnetic field amplitudes by Eaxis and Haxis, respectively. Figure 4(a) plots Eaxis (left figure) and Haxis (right figure), normalized to the incident Eaxis and Hext, respectively, as functions of x for f =10 MHz, 1 GHz, and 1.5 GHz for a dimer with d =0.2 mm. The x-coordinate range covers the whole dimer, so the gap appears as a narrow region in the figure. It is seen that Eaxis/Eext peaks sharply in the gap region, where we also find Eaxis/Eext≫ Haxis/Hext.

FIG. 4.

Normalized on-axis electric field amplitude Eaxis/Eext (left column) and magnetic field amplitude Haxis/Hext (right column) as functions of x. (a) f =10 MHz, 1 GHz, and 1.5 GHz. (b) f =2.45 GHz. (c) f =6 GHz. (D) f =10 GHz. The model is a plane wave incident on a dimer of two R =7 mm water spheres separated by a gap width of 0.2 mm on the x-y plane with its axis aligned along Eext. For all data, Eext = 400 V/cm (Hext = 1.06 A/cm) and ε (dimer) = (77.5 + 10i)ε0.

FIG. 4.

Normalized on-axis electric field amplitude Eaxis/Eext (left column) and magnetic field amplitude Haxis/Hext (right column) as functions of x. (a) f =10 MHz, 1 GHz, and 1.5 GHz. (b) f =2.45 GHz. (c) f =6 GHz. (D) f =10 GHz. The model is a plane wave incident on a dimer of two R =7 mm water spheres separated by a gap width of 0.2 mm on the x-y plane with its axis aligned along Eext. For all data, Eext = 400 V/cm (Hext = 1.06 A/cm) and ε (dimer) = (77.5 + 10i)ε0.

Close modal

The frequency range considered in Fig. 4(a) (f <1.5 GHz) has an in-water wavelength (λd > 23 mm) satisfying λd ≫ R (= 7 mm). Thus, the dimer falls in the quasi-static near zone and is dominated by the electric field. The quasi-static regime is characterized by a spatial electric field profile essentially independent of f. For example, 27 MHz is deep in the quasi-static regime. It is indistinguishable from the 10 MHz curve in Fig. 4(a) (also true in Fig. 5 below).

FIG. 5.

Simulated gap electric field Egap/Eext in log scale (a) and gap magnetic field Hgap/Hext in linear scale (b) as a function of f from 10 MHz to 10 GHz. The model is a plane wave with Eext = 400 V/cm incident on a dimer of two R =7 mm spheres on the x-y plane with its axis aligned along Eext. For ε (of spheres) = (77.5 + 10i)ε0 and four gap distances d.

FIG. 5.

Simulated gap electric field Egap/Eext in log scale (a) and gap magnetic field Hgap/Hext in linear scale (b) as a function of f from 10 MHz to 10 GHz. The model is a plane wave with Eext = 400 V/cm incident on a dimer of two R =7 mm spheres on the x-y plane with its axis aligned along Eext. For ε (of spheres) = (77.5 + 10i)ε0 and four gap distances d.

Close modal

Figures 4(b)–4(d) plot Eaxis/Eext (left column) and Haxis/Hext (right column) for f =2.45, 6, and 10 GHz for the same gap width of 0.2 mm. As f rises above ∼1.5 GHz, electromagnetic resonances start to appear one after another (band edge first). The x-axis profiles of these resonant modes in the two spheres can be seen in Figs. 4(b)–4(d). In the presence of resonances, the gap electric field is a mixture of polarization-charge field and fringe electric fields of the resonances. However, we may readily identify the polarization-charge part of the field. The steep drop of the electric field on both sides of the gap [Figs. 4(a)–4(d)] indicates a large x-gradient of Eaxis; hence, a concentration of polarization charges on the gap surfaces. In other words, the sharp peak of Eaxis in the gap region is predominantly due to surface polarizations charges.

Figure 4 also shows a rather broad frequency range of the polarization-charge effect. Eaxis/Eext in the gap decreases from 40 at f <2.45 GHz, to 30 at f =6 GHz, then to 12 at f =10 GHz, an indication that the polarization-charge effect rapidly fades away beyond 10 GHz.

We may look at the problem from a different angle to find the influence of sphere resonances on the polarization-charge field. For simplicity, we quantify the gap electric and magnetic fields by a single value, Egap and Hgap, respectively, where Egap = Eaxis (x =0) and Hgap = Haxis (x =0). Figures 5(a) and 5(b) plot, respectively, Egap/Eext (in log scale) and Hgap/Hext (in linear scale) as a function of f from 10 MHz to 10 GHz for four gap distances d. It can be seen that Egap/Eext is the superposition of a smooth polarization-charge field and a much smaller, fluctuating fringe electric field of the resonances. The smallness and fluctuating behavior of the fringe fields are evidenced by Hgap/Hext in Fig. 5(b).

In Fig. 5(a), Egap/Eext at low frequencies (<1.5 GHz) is independent of f, as is consistent with the quasi-static data in Fig. 4(a). As resonances set in at higher frequencies, the fringe fields influence the polarization-charge field positively or negatively. If suffices to focus on the region (∼3 GHz), where the maximum Egap/Eext occurs. The fringe fields are seen to add a d-dependent, 0.45–3.4% positive contribution to the quasi-static polarization charge field. Thus, even the strong electromagnetic resonances contribute to the gap electric field in a minor way.

For a detailed simulation of the electrical origin in the presence of electromagnetic resonances in the dimer body, we again consider a dimer oriented parallel to Eext. To relate to the subsequent experiment, a 2.45-GHz microwave linearly polarized in the x-direction is inputted via a rectangular horn antenna to propagate in the z-direction onto the sample [Fig. 6(a)]. Two R =7 mm water spheres with a variable gap width are placed side by side along the incident Eext at a distance of 13 cm in front of the antenna opening. At this position, the incident power varies by ∼1% from the center to the far edge of the dimer. So, the Eext on the dimer is well represented by its value at the dimer center.

FIG. 6.

Simulation model and results at 2.45 GHz for a dimer with its axis aligned along Eext. (a) The simulation model and instantaneous electric field pattern on the x-z plane. (b) The x-y plane electric field amplitude pattern near the dimer for three gap widths. (c) Corresponding x-y plane magnetic field amplitude pattern near the dimer for three gap widths. (d) The x-z plane electric field amplitude pattern near the dimer for three gap widths. (e) Corresponding x-z plane magnetic field amplitude pattern near the dimer for the three gap widths. For all figures, the dimer with ε = (77.5 + 10i)ε0 is exposed to Eext (peak) = 400 V/cm. The x-y plane in (b) and (c) and x-z plane in (d) and (e) cut through the centers of the two spheres. The color code gives the electric and magnetic fields (normalized to Eext and Hext, respectively) in logarithmic scale.

FIG. 6.

Simulation model and results at 2.45 GHz for a dimer with its axis aligned along Eext. (a) The simulation model and instantaneous electric field pattern on the x-z plane. (b) The x-y plane electric field amplitude pattern near the dimer for three gap widths. (c) Corresponding x-y plane magnetic field amplitude pattern near the dimer for three gap widths. (d) The x-z plane electric field amplitude pattern near the dimer for three gap widths. (e) Corresponding x-z plane magnetic field amplitude pattern near the dimer for the three gap widths. For all figures, the dimer with ε = (77.5 + 10i)ε0 is exposed to Eext (peak) = 400 V/cm. The x-y plane in (b) and (c) and x-z plane in (d) and (e) cut through the centers of the two spheres. The color code gives the electric and magnetic fields (normalized to Eext and Hext, respectively) in logarithmic scale.

Close modal

The simulated instantaneous electric field pattern is shown on the x-z plane in Fig. 6(a). The x-y plane electric field amplitude patterns around the dimer are shown in Fig. 6(b) for d =2, 1, and 0.2 mm. As expected, an electromagnetic resonant mode now appears in each sphere, and there are two important features as discussed below.

First, as the gap width d narrows, there is no observable change of the resonant field profiles (pattern and strength) within each sphere [Figs. 6(b)–6(e)], and the gap magnetic field remains negligibly small [Figs. 6(c) and 6(e)]. Thus, there is no evidence of mode coupling into an electromagnetic hotspot in the gap. This is because the resonant mode under study is formed of a standing wave in the body of the sphere and fringe fields on the outside. The body field of one sphere can hardly be affected by the fringe fields of a neighboring sphere.

Second, in contrast to the gap magnetic field, the gap electric field rises sharply with a decreasing d to values far greater than the body field [Figs. 6(b) and 6(d)], which again shows the polarization charges effect. Thus, in addition to generating electromagnetic resonances in the two spheres, the incident traveling wave has also induced a large polarization-charge electric field in their gap, which can cause an AC air breakdown similar to that in a capacitor.

As an additional note, these results are obtained with the dimer axis aligned along Eext. We have examined other relative orientations between the dimer and Eext. The gap electric field always maximizes when the dimer is aligned along Eext.

In a multimode cavity, the magnitude and polarization of Eext vary spatially and are highly load-dependent (see Fig. 2 of Ref. 7), which explains the stochastic nature of the sparks taking place in a microwave oven.1 The resulting difficulty in data reproducibility can be remedied by traveling-wave irradiation [see Fig. 6(a)], which is ∼99% uniform on the sample with a known magnitude and polarization.

For a traveling-wave experiment at 2.45 GHz (see  Appendix A, Part 5, for the experimental setup) in the same configuration as in Fig. 6(a), two R =7 mm hydrogel spheres are aligned along Eext and placed 13 cm in front of the antenna opening, where the magnetron power is only sufficient to produce an Eext (peak) of 100 V/cm (instead of 400 V/cm as in the simulation). Extrapolating from Fig. 5(a), we find that arcing occurs only when the two spheres are almost in contact.

Figure 7(a) displays thermal image of the dimer surface for d =0.7, 0.1 mm, and 0 after a 10-s exposure to Eext (peak) ≈ 100 V/cm, which is incident along the z-direction and linearly polarized in the x-direction. Arcing barely starts at d =0.1 mm and becomes stronger at d =0. The maximum temperature measured in situ by an IR camera on the side surfaces of the gap is 52.7 °C for d =0.1 mm and 83.2 °C for d =0 vs an initial temperature of 22 °C. Because of insufficient magnetron power, no sparks have been observed for both d =0.1 mm and d =0. However, arcing is witnessed by the thermal hotspots on gap surfaces.

FIG. 7.

In situ thermal images of dimers exposed to a 2.45 GHz Eext. (a) Thermal image of the surface of the hydrogel dimer aligned along Eext for three gap widths. (b) Thermal image of the surface of the hydrogel dimer aligned perpendicular to Eext for three gap widths. In all figures, the R = 7 mm dimer is exposed to Eext (peak) ≈ 100 V/cm, which is incident along the z-direction and linearly polarized in the x-direction. The exposure time is 10 s. Two tips of the sample holder are seen below the images.

FIG. 7.

In situ thermal images of dimers exposed to a 2.45 GHz Eext. (a) Thermal image of the surface of the hydrogel dimer aligned along Eext for three gap widths. (b) Thermal image of the surface of the hydrogel dimer aligned perpendicular to Eext for three gap widths. In all figures, the R = 7 mm dimer is exposed to Eext (peak) ≈ 100 V/cm, which is incident along the z-direction and linearly polarized in the x-direction. The exposure time is 10 s. Two tips of the sample holder are seen below the images.

Close modal

To check the wave polarization effect, we rotate the dimer in Fig. 7(a) by 90° on the x-y plane, so that its axis is aligned perpendicular to the Eext. As in the 27 MHz experiment in Fig. 3(d), no temperature rise has been observed after a 10-s exposure [Fig. 7(b)].

As a final note, the infrared camera gives the thermal image of the surface temperature. While a hotspot on gap surfaces can provide evidence of arcing, the surface temperature profile (especially after arcing) by no means represents the actual electromagnetic field distribution [e.g., Fig. 7(b) above or Figs. 3(b), 3(e), and 3(h) in Ref. 1]. For a direct display of field profiles, numerical simulations should be resorted to [e.g., Figs. 6(b)–6(e) above or Figs. 3(c), 3(f), and 3(i) in Ref. 1].

The observation of sparks in aqueous dimers in a dominantly electrical environment at 27 MHz proves the existence of an electrical origin for the phenomenon. Specifically, the strong electric field required for air arcing comes from a 2–3 orders-of-magnitude buildup of ±ρpol on opposite sides of the dimer gap. This is supported by the consistency of predicted and observed threshold conditions as well as by the absence of sparks upon a 90° rotation of the dimer relative to Eext. Further simulations over a broad frequency range as well as the observed gap arcing and sphere attraction at 2.45 GHz indicate the viability of the electrical origin even in the presence of the first few electromagnetic resonances in the dimer. These results, to our knowledge, shed new light into the basic properties of isolated dimers. It may hopefully be a useful reference for current investigation of the electromagnetic behavior of complex systems.

See the supplementary material for optical movies showing the sparks for the d =0.2 mm and d =0 at 27 MHz with the legends “Real-time sparks between two R = 7 mm hydrogel spheres with 0.2 mm gap in a 27 MHz capacitor” and “Real-time sparks between two R = 7 mm hydrogel spheres with 0 gap in a 27 MHz capacitor,” respectively. An optical movie showing the attractive motion and the eventual sparks is presented in a supplementary movie with the legend “Real-time attraction between two R =10 mm hydrogel spheres on a flat surface, d (initial) =1.5 mm and f =2.45 GHz.” Two larger-radius (R =10 mm) spheres on a flat horizontal surface were placed 2 cm from the horn opening for a greater attractive force and clearer attractive motion. The sphere attraction at 2.45 GHz serves as further evidence for the electrical nature of the gap field. The R =10 mm spheres are in strong electromagnetic resonance in the second mode. If an electromagnetic hotspot were present in the gap, its radiation pressure would have generated a repulsive force.

Another point of interest concerns the role of water vapor under 2.45 GHz radiation. The heating rate is proportional to f, ε″ (imaginary part of ε) and the square of the local electric field (Ref. 6, Chap. 7). Heating of the hydrogel spheres at 2.45 GHz, thus, evaporates far more water vapor than at 27 MHz. Either the Eext or air arcing can subsequently create a water vapor plasma, which then plays a role commensurate with its density. An experimental resolution of the water vapor effect is beyond the capability of our facility.

The authors are grateful to ANSYS, Inc. for generously providing the HFSS software used throughout this research, to Wave Power Technology for providing technical assistance on their 2.45 GHz magnetron, and to H. H. Teng for helpful discussions on the HFSS software. This work was funded by the Ministry of Science and Technology, Taiwan, under Grant No. MOST109–2112-M-002–018.

The authors have no conflicts to disclose.

K.R.C. conceived the idea and wrote the manuscript with input from all authors. L.C.L. and M.S.L., who contributed equally to this work, performed the simulation and experimental research. L.R.B. designed and fabricated the 27 MHz experimental system and supervised its use. Y.F.T. refined the experimental systems and participated in the measurements.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1. Simulation software

All simulations were performed using the Ansys HFSS (high-frequency structure simulator) software, which solves the full set of Maxwell equations by the finite element method (http://anlage.umd.edu/HFSSv10UserGuide.pdf). In our simulations, PMLs (perfectly matched layers) are assumed on a closed boundary, which encloses the entire simulated structure. PMLs are fictitious materials that fully absorb the electromagnetic fields acting upon them. This is consistent with the rapidly decaying fringe fields out of the 27 MHz capacitor and the absorbing walls of the anechoic chamber in the 2.45 GHz experiment.

2. Experimental sample

The only samples used in the experiment are the commonly available hydrogel spheres, which are polymer networks extensively swollen with water. The average weight of each is 0.019 g (dry) and 1.437 g (swollen with water to 7 mm radius). So, the water content is ∼99% by weight. Dry polymer networks were placed in tap water at room temperature for about 2 days before the experiment. The hydrogel spheres were taken out for the experiment as soon as they had expanded to the desired radius. Each sphere was used only once. Figure 8(a) shows the relative sizes of the dry polymer network and the swollen hydrogel sphere.

FIG. 8.

(a) The hydrogel sample (left) used for all experiments is a dry polymer network (right) extensively swollen with water. (b) Two hydrogel spheres (7 mm in radius) placed in a sample holder for the 27 MHz experiment. (c) Two hydrogel spheres (7 mm in radius) placed in a sample holder for the 2.45 GHz experiment.

FIG. 8.

(a) The hydrogel sample (left) used for all experiments is a dry polymer network (right) extensively swollen with water. (b) Two hydrogel spheres (7 mm in radius) placed in a sample holder for the 27 MHz experiment. (c) Two hydrogel spheres (7 mm in radius) placed in a sample holder for the 2.45 GHz experiment.

Close modal

Figures 8(b) and 8(c) show two hydrogel spheres placed in the sample holder for the 27 MHz and 2.45 GHz experiment, respectively. The sample holder for the 27 MHz experiments [Fig. 8(b)] is 4 cm in height, 4 cm in width, and 2 cm in thickness. It can be inserted between the capacitor plates (4 cm in separation) in either parallel or perpendicular orientation to Eext. The sample holder for the 2.45 GHz experiments is shown in Fig. 8(c). The radiation (with Eext parallel or perpendicular to the dimer axis) comes vertically down from the horn antenna.

The sample holder is made of a rigid polyurethane foam (Eccostock SH4). According to the manufacturer's data sheet, it is a heat insulating material that can withstand a temperature up to 135 °C. It has a density of ∼0.08 g/cc, dielectric constant of ∼1.08 (at 1 MHz), loss tangent of ∼0.0015 (at 1 MHz), and dielectric strength of ∼1.6 × 104 V/cm. For the dielectric constant and loss tangent, the manufacturer gives only 1 MHz data, but General Plastics tested their polyurethane foam at 1, 5, and 10 GHz (https://docplayer.net/50164587-Polyurethane-foam-dielectric-materials-for-use-in-radomes-and-other-applications.html), which fall within ±5% of the 1 MHz data for Eccostock SH4.

3. The infrared camera

We used a Flir Tau2 336 infrared camera to obtain thermal images of the sample during the experiment. The infrared camera has a 35 mm lens, a resolution of 336 × 256 pixels, a sensitivity of <60 mK, and a 30 Hz maximum frame rate (https://www.flir.com/products/tau-2/?model=46336035H). It is always placed in a weak field area and is also enclosed in a shielded aluminum box (except for an aperture for an observation port) for extra protection [see Figs. 9(c) and 10(b)]. We have never observed any effect by RF or microwaves on our IR camera.

FIG. 9.

The 27 MHz experimental system. (a) An overall view of the system. (b) Inside view of the capacitor box. (c) At the center is a close-up view of the capacitor plates, and on the left is an infrared camera attached to the outside wall. The infrared camera is interchangeable with a conventional camera.

FIG. 9.

The 27 MHz experimental system. (a) An overall view of the system. (b) Inside view of the capacitor box. (c) At the center is a close-up view of the capacitor plates, and on the left is an infrared camera attached to the outside wall. The infrared camera is interchangeable with a conventional camera.

Close modal
FIG. 10.

The 2.45 GHz experimental system. (a) The 2.45 GHz magnetron and diagnostic circuit. (b) The 60 × 60 × 60 cm3 box with dark walls is the anechoic chamber. The rectangular horn antenna enters the chamber through the upper side. The lens of a height-adjustable infrared camera is seen on the left side. The sample (not shown) will be placed directly below the horn antenna. The infrared camera is interchangeable with a conventional camera.

FIG. 10.

The 2.45 GHz experimental system. (a) The 2.45 GHz magnetron and diagnostic circuit. (b) The 60 × 60 × 60 cm3 box with dark walls is the anechoic chamber. The rectangular horn antenna enters the chamber through the upper side. The lens of a height-adjustable infrared camera is seen on the left side. The sample (not shown) will be placed directly below the horn antenna. The infrared camera is interchangeable with a conventional camera.

Close modal
4. The 27 MHz experimental system

An overview of the 27-MHz applicator system is shown in Fig. 9(a). The applicator shielding box with tuner controls is on the right. A 1.5 kW CW power tube amplifier is on the left. A power meter measures forward and reflected powers with a trip relay for system protection. A 5-W transmitter source, a variable attenuator, and a 300-W solid state amplifier are used to drive the power tube amplifier. An rf SPDT switch connects the power source to the applicator or a load. All power transmission lines are LMR-400 type low loss coaxial cable with 100% shield coverage.

Figure 9(b) shows the inside of the applicator box with capacitor plates. The rf matching configuration is a pi C-L-C matching circuit that can transform 50 Ω impedance power up to at least 105 Ω. The input variable capacitor is 25 to 800 pF rated for 1 kV peak, and the variable rf power inductor is 0–15 μH modified to withstand 15 kV peak. The unloaded Q of this circuit is ∼150 but varies with the capacitor size and spacing. When there is no rf absorbing material between the capacitor plates, nearly 100% of the rf loss is in the inductor. The blower-cooled inductor can handle 300 W of power absorption. When an rf absorbing load is placed between the plates, higher input power can be applied. A capacitor divider probe in the bottom plate is calibrated by applying a few watts of rf power and measuring the total capacitor voltage at a particular plate spacing with a high frequency 300+ MHz scope probe. The estimated measurement error is ±10%. All components inside the shielding box [Fig. 9(b)] are rigidly fixed. The box can be rotated by 90°, so that the capacitor becomes vertical.

5. The 2.45 GHz experimental system

The 2.45 GHz experimental system (Fig. 10) is the modified version of a 24 GHz system.21 The microwave source is a commercial 3-kW, 2.45-GHz magnetron. The 2.45 GHz wave passes through a protection device (a circulator) and a diagnostic circuit for power and frequency measurements [Fig. 10(a)] before entering into a 60 × 60 × 60 cm3 anechoic chamber [Fig. 10(b)] via a rectangular horn antenna [same dimensions as in Fig. 6(a)].

The chamber walls are lined with pyramid-shaped radiation C-RAM SFC-HC-3 absorbers (>99% absorption) to prevent reflections [Fig. 10(b)]. So, the sample is irradiated by a nearly un-interfered traveling wave. A height-adjustable infrared camera is mounted outside the chamber wall to take thermal images of the sample during the experiment.

Compared with the multimode standing waves in a cavity, the traveling wave allows polarization and magnitude control of the fields on the sample, hence much higher data reproducibility. However, this is at the expense of field strength. For the same input power, the field strength of a traveling wave can be much weaker than that of a standing wave in an unloaded resonant structure.

1.
H. K.
Khattak
,
P.
Bianucci
, and
A. D.
Slepkov
, “
Linking plasma formation in grapes to microwave resonances of aqueous dimers
,”
Proc. Natl. Acad. Sci. U.S.A.
116
,
4000
3005
(
2019
).
2.
R.
Bansal
, “
The wrath of grapes?
,”
IEEE Microwave Mag.
20
,
11
12
(
2019
).
3.
Research Highlight
. “
Why microwaving a grape sparks a fiery glow
,”
Nature
566
,
429
(
2019
).
4.
A.
Fox
, https://www.sciencemag.org/news/2019/02/why-sparks-fly-when-you-microwave-grapes for “
Why sparks fly when you microwave grapes
,” Science Online (
2019
).
5.
A.
Slepkov
, “
Fruit photonics and the shape of water
,”
Phys. Today
73
(
6
),
62
63
(
2020
).
6.
J. D.
Jackson
,
Classical Electrodynamics
, 3rd ed. (
John Wiley
,
New York
,
1998
).
7.
J.
Monteiro
,
L. C.
Costa
,
M. A.
Valente
,
T.
Santos
, and
J.
Sousa
, “
Simulating the electromagnetic field in microwave ovens
,” in
SBMO/IEEE MTT-S International Microwave and Optoelectronics Conference (IMOC 2011), Natal.
(IEEE, Piscataway, NJ,
2011
), pp.
493
497
.
8.
D.
Halliday
,
J.
Walker
, and
R.
Resnick
,
Fundamentals of Physics
, 10 ed. (
Wiley
,
Hoboken, NJ
,
2014
), Chap. 25.
9.
D. R.
Lide
,
CRC Handbook of Chemistry and Physics
, 85th ed. (
CRC
,
Boca Raton, FL
,
2004
).
10.
D.
Rodriguez
,
R. S.
Gorur
, and
P. M.
Hansen
, “
Effect of humidity on the breakdown characteristics of air in non-uniform fields at 30 kHz
,”
IEEE Trans. Dielectr. Electr. Insul.
17
,
45
52
(
2010
).
11.
A.
Birnboim
,
J. P.
Calame
, and
Y.
Carmel
, “
Microfocusing and polarization effects in spherical neck ceramic microstructures during microwave processing
,”
J. Appl. Phys.
85
,
478
482
(
1999
).
12.
S.
Horikoshi
,
A.
Osawa
,
S.
Sakamoto
, and
N.
Serpone
, “
Control of microwave-generated hot spots. Part V. Mechanisms of hot-spot generation and aggregation of catalyst in a microwave-assisted reaction in toluene catalyzed by PD-loaded AC particulates
,”
Appl. Catal. A
460–461
,
52
60
(
2013
).
13.
X.
Qiao
and
X.
Xie
, “
The effect of electric field intensification at interparticle contacts in microwave sintering
,”
Sci. Rep.
6
,
32163
(
2016
).
14.
N.
Haneishi
,
S.
Tsubaki
,
E.
Abe
,
M. M.
Maitani
,
E.
Suzuki
,
S.
Fujii
,
J.
Fukushima
,
H.
Takizawa
, and
Y.
Wada
, “
Enhancement of fixed-bed flow reactions under microwave irradiation by local heating at the vicinal contact points of catalyst particles
,”
Sci. Rep.
9
,
222
(
2019
).
15.
M.
Liu
,
L.
Zhang
,
W.
Zhou
,
J.
Gao
,
P.
Wang
, and
G.
Zhao
, “
The mechanism of microwave-induced discharge between submillimeter active coke
,”
Plasma Sources Sci. Technol.
29
,
075015
(
2020
).
16.
G.
Mie
, “
Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen
,”
Ann. Phys.
330
,
377
445
(
1908
).
17.
C.
Gabriel
,
S.
Gabriel
,
E. H.
Grant
,
B. S. J.
Halstead
, and
D. M. P.
Mingos
, “
Dielectric parameters relevant to microwave dielectric heating
,”
Chem. Soc. Rev.
27
,
213
223
(
1998
).
18.
R. F.
Harrington
,
Time-Harmonic Electromagnetic Fields
(
IEEE Press
,
New Jersey
,
2001
), Chap. 6.
19.
P. J. W.
Debye
,
Polar Molecules
(
Chemical Catalog Co
.,
New York
,
1929
).
20.
U.
Kaatze
, “
Complex permittivity of water as a function of frequency and temperature
,”
J. Chem. Eng. Data
34
,
371
374
(
1989
).
21.
W. Y.
Chiang
,
M. H.
Wu
,
K. L.
Wu
,
M. H.
Lin
,
H. H.
Teng
,
Y. F.
Tsai
,
C. C.
Ko
,
E. C.
Yang
,
J. A.
Jiang
,
L. R.
Barnett
, and
K. R.
Chu
, “
A microwave applicator for uniform irradiation by circularly polarized waves in an anechoic chamber
,”
Rev. Sci. Instrum.
85
,
084703
(
2014
).

Supplementary Material