## How to make a giant bubble

Stephen Frazier, Xinyi Jiang, and Justin C. Burton ∗ Department of Physics, Emory University, Atlanta, Georgia 30322, USA

(Dated: November 25, 2019)

Using mixtures of soap, water, and long chain polymers, free-floating soap bubbles can be formed with volumes approaching 100 m . 3 Here we investigate how such thin films are created and maintained over time. We show how the extensional rheology is the most important factor in creating the bubble, and how polydispersity in molecular weight of the solvated polymers leads to better performance at lower concentrations. Additionally, using IR absorption, we measure soap film thickness profiles and film lifetimes. Although the initial thickness mostly depends on the choice of detergent, polymers can dramatically increase film lifetime at high molecular weights and high concentrations, although such high concentrations can inhibit the initial film formation. Thus, the ideal concentration of polymer additives for making giant bubbles requires a robust viscoelastic rheology during extension, and is aided by long film lifetimes during gravitational drainage and evaporation.

∗ justin.c.burton@emory.edu

## I. INTRODUCTION

For the last 2 decades, giant bubble enthusiasts have been creating soap film bubbles with ever increasing volumes. As of 2019, the world record for a free-floating soap bubble stands at 96.27 m 3 [1]. For a spherical bubble, this corresponds to a diameter of 5.7 m and a surface area of 101 m 2 . A simple glance at the multitude of colors from the reflecting bubble film suggests a film thickness of order a few microns. Thus, the extension-to-thickness ratio of this film is nearly 5 × 10 , which is quite staggering considering that a single hole can lead to the film's demise. 5 How are such large films created, and how do they remain stable? These questions pose interesting fluid mechanics problems, but are also important considering the vast role that soap bubble films have and continue to play in physics research [2-6] and education [7-9]. Additionally, the formation, mechanics, and stability of foams with various additives is an important environmental problem. Such foams contribute to long-lasting pollution in rivers and waterways contaminated with industrial run-off [10, 11].

For those interested in making giant bubbles, the Soap Bubble Wiki [12] contains a wealth of empirical information and recipes for optimal bubble solutions. Most solutions use industrial dish detergents as a surfactant for the soap films. For soap films supported solely by surface tension forces, one may estimate the maximum size of a bubble based on the Bond number, Bo ≡ ρgzR /σ b , where R b is the radius of the bubble, z is the thickness of the film, ρ is the liquid density, σ is the surface tension, and g is the acceleration due to gravity [13]. The Bond number measures the ratio of gravitational forces to surface tension forces. Assuming Bo = 1 for films with d ≈ 2 µ m, the maximum size is R b ≈ 1 5 m. . Although Marangoni forces can help support larger films, up to a few meters [14-16], additional forces are necessary to make world-record giant bubbles.

The key ingredient is the addition of long chain polymer molecules. Two of the most common types are guar gum and polyethylene oxide (PEO). When solvated in water, guar forms long polysaccharide chains of polydisperse molecular weight, typically between 0.25-5.0 × 10 6 g/mole [17]. PEO is more well-controlled, as samples of monodisperse molecular weights can be obtained, although most industrial samples are polydisperse. Both additives drastically alter the viscoelastic rheology of the soap solution in different ways, and are partially confined in soap films since their radius of gyration is fractions of a micron. Despite this wealth of information, the physics of making giant bubbles with polymer solutions is poorly understood. Properties such as molecular weight, concentration, and polydispersity are exceedingly important, as well as the solution's overall rheology and choice of surfactant. It is no wonder a complete on-line Wiki exists that is devoted to the fine-tuning of giant bubble recipes.

In this paper we identify some of the underlying physical mechanisms that give rise to giant bubbles. We study both the shear and extensional rheology of common solutions of water, surfactant, and PEO or guar. The most robust solutions for making bubbles have intermediate concentrations and a polydisperse mixture of polymers of various molecular weights, allowing a large volume of liquid to be continuously drawn into a film without breaking. Additionally, we measure the thickness of soap films over time using IR absorption and find that the lifetimes increase at larger concentrations than those often used in bubble solutions, indicating that polymers may also enhance film longevity after evaporation and gravitational drainage have occurred. We suggest that an optimal polymer solution making giant bubbles comes from a combination of the extensional rheology of isolated chains, and cooperative interactions between polymers of differing molecular weights.

## II. EXPERIMENTAL METHODS

The basic solutions used in our experiments consisted of a mixture of detergent and water. Unless otherwise noted, all solutions were a mixture of 4% v/v Dawn Pro Dish Detergent and deionized water. For some experiments, we used a detergent common in many soap film experiments, sodium dodecyl sulfate (SDS), which was obtained from MilliporeSigma. Long-chain polymers were then added to each soap solution. We measured the surface tension of all soap solutions using axisymmetric drop shape analysis [18]. Regardless of polymer concentration, all solutions had a surface tension of 32 ± 2 mN/m. Guar powder was obtained from MilliporeSigma. The powder was directly added to the soap solution, then magnetically stirred at 50 ◦ C for 4 hrs to allow for full dissolution of the particles. Guar is often dissolved into an alcohol slurry first to prevent clumping. However we found that heating and stirring worked just as well. Solutions were ready to use after cooling to room temperature (22 ◦ C). For making giant bubbles, guar is often added in concentrations ranging from 1.0-2.0 g/l [12], although we explored a wider range in our experiment.

PEO was also added to soap solutions. First, we used a common industrial PEO lubricant for making giant bubbles, J-Lube (Jorgensen Labs), which is composed of 25% polydisperse PEO (up to 8 × 10 6 g/mole) and 75% sucrose. We did not characterize the composition of the J-Lube sample beyond information from the supplier. For making giant bubbles, J-Lube is often added in concentrations of 0.1-0.4 g/l [12], 25% of which is actually PEO, so the concentrations are much lower than guar. In addition to commercial PEO, we used monodisperse samples of PEO from MilliporeSigma with molecular weights of 0.1M, 0.6M, 1.0M, 2.0M, 4.0M, and 8.0M g/mole, where M = 10 . 6

FIG. 1. ( A ) Shear rheology of the soap + water + polymer solutions at 25 ◦ . The black squares are the control (no polymer), and the blue and red points represent soap solutions with polymers added in the indicated concentrations. ( B ) Schematic of a soap film being pulled at characteristic velocity v from a rope with a circular cross-section of diameter d . Away from the rope, the flow is mostly extensional.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-69432d74-3538-4e2d-b5c2-366bce074f3b.jpg)

These are viscosity-averaged molecular weights ( M v ), which lies between the number-averaged and weight-averaged molecular weights. The polydispersity index was not available for these samples. All powdered samples of PEO and guar were stored in opaque containers and placed in the refrigerator until use. Some PEO was intentionally degraded by aging for 6 months in room temperature conditions, as noted.

## III. SHEAR RHEOLOGY AND EXTENSIONAL PROPERTIES

For an initial characterization of the bubble solutions, we employed steady-state rheometry to measure the concentration and shear rate dependence of the shear viscosity. We used a TA Instruments AR2000 with a parallel-plate configuration and a 0.5 mm gap. A Peltier plate was used to control the temperature of the sample. Figure 1A shows the results from various concentrations of guar and J-Lube, as well as the pure soap solution. For the J-Lube solutions, the viscosity varies by less than a factor of 2, even for concentrations larger than those used for making giant bubbles ( ≈ 0.1-0.4 g/l). All the guar solutions tested showed some degree of shear-thinning, and the shear viscosity increased dramatically with concentration at low strain rates. Given that both 1.5 g/l and 2.4 g/l are concentrations used for giant bubbles, and their viscosities are approximately 10 times that of the J-Lube solutions, we conclude that the shear viscosity is not the leading factor in determining a solution's ability to create giant bubbles.

Figure 1B shows a representation of a soap film being pulled from a rope of diameter d with a characteristic velocity v . Near the rope where the velocity is zero at the boundary, the fluid is sheared with a characteristic shear rate ˙ γ ≈ v/d . For a 1 mm rope with a film velocity of U = 1 m/s, the shear rate is 1000 s -1 . This is an upper bound and the highest shear rate we tested in our shear rheology. However, away from the rope, the flow is mostly extensional due to slow variation in film thickness in the pulling direction [19]. Thus, we expect the bulk extensional viscosity to be more important for creating giant bubbles. For solutions of long chain polymer molecules, it is well-known that

FIG. 2. ( A ) Image of a soap + water solution at the moment the connecting thread ruptures near the top of the drop. ( B ) A solution with 3.0 g/l of guar added. The thread is highly extended, and will reach ≈ 3 cm in length before rupture. The scale bar applies to both images. ( C ) Thread rupture length vs. polymer concentration for J-Lube and guar. The gray and light red rectangles indicate the range often used to make giant bubble solutions. Error bars represent the standard deviation of multiple drops, although the experiments were extremely repeatable.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-cbd97607-48e2-46b8-bfc6-473a2a762e39.jpg)

the extensional properties change dramatically at very low concentrations where changes in the shear viscosity are nearly immeasurable [20-24].

To characterize the extensional viscosity of each sample, we used high-speed video to analyze the falling of droplets under gravity, and the dynamics of the viscoelastic thread that they pulled in their wake. Although previous studies have used the break-up rate of a liquid bridge to characterize non-Newtonian extensional rheology [25-27], we chose to measure the length of the viscoelastic thread prior to its final rupture since this is more analogous to the continuous pulling of a soap film. Bulk liquid was filled into a glass burette, and the flow rate was adjusted so that one drop fell every 10 s. To image the drops, we used a Vision Research Phantom v7.11 camera at 2000 frames-per-second with a Tonkia macro lens, resulting in a resolution of 70 pixels/mm. The liquid drops were illuminated from behind with a 150 W tungsten-halogen lamp.

Figure 2A-B shows images of two drops immediately after detachment from the joining thread. Without the addition of polymers, the thread resembles that of pure water and other low-viscosity pure liquids [28]. When 3.0 g/l of guar is added, the thread length increases dramatically. At higher concentrations, the length of the thread prior to rupture seems to diverge, and the drop leaves the imaging region before the thread ruptures. While this behavior is monotonic for the guar solutions, the PEO-based J-Lube solutions oscillate between long and short threads before finally diverging in length around 0.7 g/l (Fig. 2B). We attribute this non-monotonic behavior to the formation of 'beads-on-a-string' during extension in the lower-viscosity J-Lube solutions [21, 29]. Guar solutions typically have a 10 fold larger shear viscosity, which suppresses the formation of beads through inertial forces [29]. Rupture of the thread can then occur where beads are connected to the thread rather than at the main drop.

Although the particular placement of the oscillations shown in Fig. 2C are likely specific to our dripping experiment, the data shows that the addition of either J-Lube or guar has a significant effect on the extensional properties of the solutions for concentrations used to make giant bubbles. However, the molecular weight distribution of both polymer sources is unknown. Both contain high molecular weight ( &gt; 1 M ) components, but this alone is not sufficient, and is potentially detrimental to making giant bubbles. Many enthusiasts prefer to 'age' their PEO in the powder form for many weeks or months prior to usage [12]. Solutions of newly-purchased J-Lube and other industrial sources of PEO are often very sensitive to concentration and can become quite 'stringy' and elastic. This is evident in the Fig. 2B where the rupture length begins to diverge.

Aging the PEO in solid-state involves photo-induced or thermal degradation where larger chains break into smaller ones, and the distribution of molecular weight broadens [30, 31]. Degradation can also occur in solution as a result of aging or strong fluid flow [32-34]. The concentration at which chains begin to overlap and interact, denoted by c ∗ , depends on the molecular weight. At lower concentrations, the rheology of very dilute solutions depends strongly on the nonlinear extensional properties of individual, solvated polymer chains. Thus we investigated the extensional properties of soap solutions with monodisperse PEO of different molecular weights.

Table I lists the molecular weight, radius of gyration ( R g ), and overlap concentration ( c ∗ ) for the samples used in our experiments. The radius of gyration for PEO was calculated using the formula R g = 0 02 . M 0 58 . v [nm], in accordance with references [35-37], and the overlap concentration was calculated as c ∗ = M / v (4 πR N / 3 g A 3) [g/m ], where 3 N A is

FIG. 3. ( A ) Images of a falling drop with 50.0 g/l of 0.6M PEO. The formation of a distinct bead induces rupture of the thread. ( B ) Images of a falling drop with 0.05 g/l of 8.0M PEO. The drop exits the viewing area before rupture of the thread. ( C ) Rupture length vs. concentration for 5 different molecular weights of PEO. Error bars represent the standard deviation of multiple drops.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-94d58774-28ae-47e2-9bd6-1d748631e8dc.jpg)

TABLE I. Molecular weights of PEO used in the experiments and their corresponding radii of gyration ( R g ) and overlap concentration ( c ∗ ) assuming they are well-solvated. The quantities R g and c ∗ were computed as described in the text and references [35-37].

| PEO MW   |   R g (nm) |   c ∗ (g/l) |
|----------|------------|-------------|
| 0.1M     |         16 |        9.9  |
| 0.6M     |         45 |        2.63 |
| 1.0M     |         60 |        1.8  |
| 2.0M     |         90 |        1.08 |
| 4.0M     |        135 |        0.65 |
| 8.0M     |        202 |        0.39 |

Avogadro's number. We have assumed that the polydispersity of each sample is small enough so that the viscosity and weight-averaged molecular weights are the same (i.e. M v = M w ).

In analogy to Fig. 2, Fig. 3A-B shows images of falling drops of PEO solutions of different molecular weights. For lower molecular weights, 1.0M or less, the increase in the rupture length coincides more closely with c ∗ , indicating that polymer blob overlap and transient entanglements may contribute to the extensional rheology. High molecular weights, above 2.0M, lead to divergent rupture lengths above a certain concentration less than c ∗ (Tab. I). The nonlinear rheology is apparent; an 8.0M solution produces extremely long threads at concentrations 1000 times less than a 0.6M solution. This shows that for very high molecular weights, the entropic cost of extension for individual chains in the flow leads to viscoelastic properties and the formation of long threads.

One major difference between most of the data shown in Fig. 3 and the data for J-Lube shown in Fig. 2 is the nearly monotonic behavior of the monodisperse PEO, and the oscillations in the polydisperse J-Lube solutions. A polydisperse solution will certainly contain a broader range of relaxation times in its rheological response [23, 38-42], possibly giving rise to a non-monotonic extensional response. The rheology of a solution at very low concentrations

FIG. 4. ( A ) Rupture length vs. concentration for freshly-purchased 2M PEO and for 2M PEO that has been aged for 6 months at room temperature. ( B ) Rupture length vs. total concentration for 4M PEO, 2M PEO, and a 50/50 mixture of each sample. Error bars represent the standard deviation of multiple drops in both graphs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-a9518d24-22c5-4180-ab20-0b56a00d3d65.jpg)

would depend mostly on the longest polymers in the polydisperse mixture. However, one would expect that as a collection of polymers ages or degrades, the individual polymers would not lengthen, and the distribution would only widen toward shorter lengths.

Our results with aged and polydisperse mixtures suggest a more cooperative behavior. Figure 4A shows the rupture length for falling drops in a freshly-purchased sample of 2M PEO, and a sample that has been aged for 6 months at room temperature in an opaque bottle. Surprisingly, the aged sample shows a stronger extensional response at lower concentrations. We can also observe the same effect by simply mixing a 2M and 4M PEO sample together. For the same total weight of polymer, the mixture displays an increase in the length of the thread at lower concentrations than either monodisperse sample. This behavior is only possible if there are interactions between the longer and shorter chain polymers, i.e. they can not be considered as dilute, individual molecules.

Although the exact nature of these interactions is unclear, we hypothesize that some degree of clustering of the longer-chain polymers may give rise to the measured results. The smaller polymers may act as depletants, leading to a weak, attractive force between the longest and largest polymers in the solution. Transient entanglements between long polymers in close proximity could then lead to extensional resistance at low bulk concentrations, with an ultimate elastic resistance limited by the strength of the transient entanglement. To our knowledge, evidence for this behavior has not been reported in the literature, and is left for future studies focused on extensional rheology. Nevertheless, it is evident that the polydispersity in molecular weight makes the extensional properties of bubble solutions less sensitive to changes in concentration, and thus more robust for making giant bubbles.

## IV. LIFETIME AND THICKNESS OF SOAP FILMS

There is approximately 300 ml of liquid in a giant spherical bubble with a surface area of 100 m 2 and an average thickness of 3 µ m. Having the right balance of extensional properties that resist thinning and breaking but allow for sufficient flow rates is one part of the puzzle in making giant bubbles. These bubbles must also survive for a considerable amount of time. However, the thinning of soap films is quite complex [43]. Although vertical soap films will drain and thin under the influence of gravity, evaporation also plays a key role [44]. Low humidity and high temperatures are well-known barriers to long-lived, giant bubbles [12]. Assuming a constant thickness, both the volume of the liquid and mass evaporation rate should be proportional to the surface area of the film, so to first approximation, evaporation should affect both small and large films in the same manner, provided that the diffusion of water vapor away from the surface is similar.

Marangoni forces and film stretching can also significantly affect the lifetime of soap films without polymers [45]. The type of surfactant and its overall concentration will determine the flow profile in the film by altering the the degree of 'slip' at the liquid-air interface [46-48], and positive vertical gradients in the surfactant concentration at the surface will pull upwards on the film [14, 15]. Furthermore, the 2D surface rheology can influence the flow and overall drainage of the film [19, 49-51]. Finally, polymers have been shown to affect the initial thickness of the film at high concentrations, near or above c ∗ , displaying strong deviations from Frankel's law [52-54]. Polymer-surfactant interactions may also be important [53, 55], these systems require more quantitative studies. Nevertheless, nearly all

FIG. 5. ( A ) Schematic of the experimental apparatus used for thickness measurements. The infrared light from the LED is modulated by the optical chopper and detected using a lock-in amplifier. The height of the measurement ( h ) can be adjusted for each experiment. ( B ) Three separate experiments of the thickness vs. time at h = 7.5 cm for a soap + water film with no polymers. The films popped after ≈ 25-30 s. ( B ) Film thickness profiles at t = 5 s for three different bubble solutions. Each data point is a separate experiment. The thickness error for each measurement is ≈ 100 nm, and is mostly due to the variation between experiments.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-2c31a998-6982-4cfe-b888-832ce931e8b4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-916fb923-78f2-4466-9d77-9ee3751147b1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-e11273fc-acb0-4dcb-b883-31a336a9d01c.jpg)

of these experiments use soap films of a few centimeters, well below the size of giant bubbles.

It is not clear if large soap films with polymers last longer by producing thicker films, slowing gravitational drainage, or slowing evaporation. To answer some of these questions, we performed time-dependent thickness measurements of soap films using a custom-built infrared absorption apparatus. Our design is similar to those used in previous measurements of soap films [56]. Figure 5A shows a schematic of the setup. Two cotton strings (1 mm diameter) were anchored in a liquid reservoir containing the soap solution, connected by a third cotton string at the top. The strings were submerged into the solution and then quickly raised and made tight in order to produce a rectangular soap film with dimensions 10 cm × 15 cm. The film was allowed to drain under gravity until it popped.

An infrared LED (Boston Electronics) with a peak wavelength λ = 3 µ m was used as a semi-collimated source. The light was passed through the center line of the film and detected by a photoconductor (Thor Labs) on the other side. The signal was modulated at 600 Hz by an optical chopper and detected by a lock-in amplifier (SRS 830, 1 MΩ input impedance). The impedance of the photoconductor is approximately 1 MΩ, but varies by a small amount based on the intensity of the incoming light. The capacitor shown in the figure is used to block the DC component of the signal. The optical elements were mounted on an adjustable stage so that the height of the beam from the

FIG. 6. Film thickness profiles at 3 different times for a soap + water solution ( A ), a solution with 0.5 g/l of 1M PEO ( B ), and a solution of 0.5 g/l of 8M PEO ( C ). Each data point is a separate experiment. The thickness error for each measurement is ≈ 100 nm, and is mostly due to the variation between experiments.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-f42973ed-0022-4146-9136-5035795233d4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-4b57890f-5a93-49c4-bc5c-befb5e85e3c8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-ea0288d3-1a44-4dc2-a3db-1c036f793d08.jpg)

liquid bath, h , could be varied for each experiment. A large acrylic box was placed around the entire experimental apparatus in order to control for ambient air currents and to keep the humidity stable.

Water has a strong, sharp absorption peak at λ = 3 µ m, so that the exponential extinction length is z 0 /similarequal 0 9 . µ m [56]. Taking into account the reflections at the incident interface and within the film (index n = 1.17), the ratio of the transmitted to incident intensity is

$$\frac { I _ { T } } { I _ { 0 } } = \frac { ( 1 - R ) ^ { 2 } e ^ { - z / z o } } { 1 - R ^ { 2 } e ^ { - 2 z / z o } },$$

where R = ( n -1) 2 / n ( +1) 2 and z is the local film thickness. However, as discussed in reference [56], R can be safely ignored since R = 0 0061 for . n = 1 17. Thus, after verifying the linearity of our optical detector using optical density . filters to reduce the intensity of the incident beam, we used the following relationship to convert the voltage measured to the film thickness:

$$V ( z ) = V _ { 0 } e ^ { - z / z o },$$

where V 0 is the voltage measured in the absence of the film.

Figure 5B shows thickness vs. time data for three separate soap films with no polymers at h = 7.5 cm. The initial rise is the raising of the cotton strings to create the film, and the slow decay is the thinning of the film until it eventually pops at t ≈ 25 s. We found that there was often more variability in the final film lifetime than in the initial thickness and drainage rate. This may be expected since the film can pop due to an instability that forms anywhere in the film. Most films popped when some portion became thinner than 500 nm, although some films with polymers reached smaller thicknesses. Perhaps surprisingly, the overall film thickness did not change much with the addition of polymers. Figure 5C shows initial profiles of the film thickness vs. height for pure soap + water, a solution with J-Lube, and a solution with guar. Data for each height was taken with different, identical soap films since the apparatus had to be raised and lowered between each experiment. As will be shown later, the film thickness depended much more strongly on the choice of surfactant rather than the addition of any polymers.

The overall lifetime of the film strongly depended on the presence of polymers in the solution. Figure 6 shows the thickness profiles at 3 different times for soap + water (A), and 0.5 g/l solutions of 1M and 8M PEO. The qualitative shape of the profiles was basically the same, although there was some variation near the top of the film for the 8M PEO sample (Fig. 6D), possibly due to a two-dimensional 'beads-on-a-string' effect where thick and thin regions exist simultaneously. However, the concentrations required to significantly increase the lifetime were distinctively larger than that needed to produce long, viscoelastic threads in falling drops, as shown in Fig. 3. The slow, linear drainage and long timescale (more than 10 seconds) suggests that the flow profile in the film is more parabolic than uniform, possibly due to an immobile surfactant boundary for large concentrations of soap [43, 46-48]. In this case, the polymers experience a simple Poiseuille-type flow since shear-thinning behaviors require higher concentrations (Fig. 1) where individual polymer chains overlap above c ∗ . If the velocity profile in the film was more like 'plug' flow, then we would expect and exponential dependence on time [57].

This concentration and molecular weight dependence can be seen in Fig. 7. The optimal concentrations for making giant bubbles (J-Lube and guar; gray and light red regions in Fig. 7A) show that film lifetime is not the most

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-68e0d2e2-82b8-45c1-8c9b-afec35454bb2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-0fa9a7f4-a795-4162-a311-9ad31b303d4f.jpg)

FIG. 7. ( A ) Film lifetime vs. concentration for both J-Lube and guar bubble solutions. The gray and light red boxes indicate the concentrations typically used to make giant bubbles. ( B ) Film lifetime vs. concentration for monodisperse PEO polymer solutions. Reported error bars represent the variation among repeat experiments.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00537v2-figure-cac0390e-0534-4813-9bd1-3fb469b8096c.jpg)

400

B

film lifetime (s)

350

300

250

200

150

100

50

0

0.3

0.4

1.5 g/l guar

0.5

0.6

0.7

0.8

0.9

1.0

relative humidity

FIG. 8. ( A ) Average lifetimes of films formed from different bubble solutions. The guar concentration used was 2.4 g/l, and was 0.5 g/l for all J-Lube and PEO solutions. Note that J-Lube is only 25% PEO by weight. ( B ) Film lifetime vs. relative humidity for a 1.5 g/l guar solution. Error bars for all data represent the variation among multiple experiments.

important property since the lifetime for guar is much longer than for J-Lube at the optimal concentrations. However, the correlation between the extensional properties and the maximum lifetime for guar solutions may be what makes it an excellent choice for making large, stable soap bubbles. The lifetime also depends on the molecular weight of the individual polymers. Figure 7B shows that the lifetime of a film typically increases with concentration, then decreases at high concentrations. For 1M PEO, the lifetime can be longer than 80 s, whereas a film with 8M PEO never lives longer than 60 s. We suggest that the increase in the shear viscosity of the solution alongside the increase in extensional properties for lower molecular weight polymers is the primary reason for such long lifetimes, whereas for high molecular weight polymers, the film is more elastic, and may be more susceptible to breakage and failure.

A summary of lifetimes from various bubble solutions is shown in Fig. 8A. Also shown are lifetimes for films made with a commonly used detergent for making soap films, SDS. For small SDS soap films of order ∼ 1-2 cm , lifetimes 2 of the films can reach 30 s or more [48, 58], often with the help of Marangoni stresses. We found that for films of order ∼ 100 cm , SDS resulted in extremely short lifetimes. 2 Even though the volume concentrations are similar, Dawn Pro dish detergent dramatically increases the stability of the film. The main conclusion from these data is that although lifetime itself does not determine a bubble solution's ability to make giant bubbles, it is certainly not a detriment. As soon as a film is pulled, gravity will act to drain the liquid from parts of the film. Additionally, the water will begin to evaporate from the huge surface area, also thinning the film. The latter will increase the concentration of any polymer, and help to increase the lifetime and slow the drainage. As Fig. 7B shows, a polydisperse mixture of molecular weights can help to increase the lifetime over a broader range of concentrations.

Finally, as is known to many bubble enthusiasts, evaporation will decrease the lifetime of bubbles no matter which polymer is used. A recent experimental study showed that the time to rupture after the formation of a centimeterscale soap film strongly depends on humidity [44]. In our experiments the effect is dramatic. This is shown in Fig. 8B for a 2.4 g/l guar solution. The relative humidity (RH) in the experimental enclosure was increased using a cup of warm water exposed to the air. A slow increase in the lifetime was observed for low humidity, and then at ≈ 75% RH, the lifetime began to increase, reaching almost 250 s at 90% RH. Considering this significant increase, we can conclude that much of the film thinning occurs through evaporation. As mentioned previously, the liquid volume and mass evaporation rate both scale with surface area, so for larger soap films, the removal of the evaporative boundary layer through convection or ambient wind is likely to be important. It is no wonder that many bubble enthusiasts prefer cool, humid days for making giant bubbles.

## V. CONCLUSION

In conclusion, we have investigated the properties of surfactant and polymer solutions commonly used to make giant bubbles with surface areas approaching 100 m 2 and thicknesses of a few microns. The most important additives for making these bubbles are hydrophillic, high molecular weight polymers, such as guar gum and PEO. We found that it is the extensional rheology of these dilute, polymer solutions that is the most important factor in creating the films. Although PEO is mildly surface-active and can absorb to the air-water interface [53, 59], we are not aware of any studies on its surface rheology. Given that the range of polymer concentrations used for giant bubbles agree well with the bulk extensional properties, we suspect that surface viscoelastic effects play a minor role in the formation of giant bubbles.

Assuming a pulling velocity of U ≈ 1 m/s, the capillary number associated with pulling a giant bubble is Ca ≡ Uη/σ ∼ 0.05, were η is the bulk viscosity, assuming it is Newtonian. This is much larger than those used in typical lab experiments [43, 47, 53]. If a soap film is pulled too fast, then surfactants do not have time to migrate to the newly-created free surface and stabilize the film [14]. Additional elastic forces from added polymers are thus necessary to prevent film breakage during the formation of a giant bubble. Although the lifetime of the film does increase with polymer concentration, this is a secondary effect that occurs at higher concentrations, and may influence the stability of the film during evaporation and thinning. Other than the concentration, molecular weight, and polydispersity of the polymers, there are other factors that are often tuned to make the best bubble solutions, such as pH buffering, which is usually done with the addition of baking powder. Although, we expect this mostly affects the agglomeration of solvated polymers, and the ensuing shelf life of the solution.

One main finding of this work that requires further investigation is the cooperative influence of polydisperse solutions. Mixtures of two different molecular weights, or aged solid samples of polymers would produce noticeable extensional properties at concentrations much lower than required in monodisperse solutions. This can not be explained by the properties of dilute, polymer chains. In polydisperse polymer melts, cooperative effects have been identified, such as nematic phase separation [60]. Interactions between the surfactant micelles and polydisperse polymers may also lead to a rich set of emergent interactions [61], which may also depend on the ambient flow of the solvent [62, 63]. Dedicated extensional rheology experiments and molecular dynamics simulations may help to shed light on the nature of these cooperative effects in future studies.

## ACKNOWLEDGMENTS

This work was supported by the NSF DMR Grant No. 1455086.

[1] Guinness World Records, https://www.guinnessworldrecords.com/world-records/largest-free-floating-soap-bubble , [Accessed: 29-June-2019].

- [2] R. E. Goldstein, H. K. Moffatt, A. I. Pesci, and R. L. Ricca, Soap-film m´ 'obius strip changes topology with a twist singularity, Proc. Nat. Acad. Sci. 107 , 21979 (2010).
- [3] M. Rivera, P. Vorobieff, and R. E. Ecke, Turbulence in flowing soap films: Velocity, vorticity, and thickness fields, Phys. Rev. Lett. 81 , 1417 (1998).
- [4] L. Ristroph and J. Zhang, Anomalous hydrodynamic drafting of interacting flapping flags, Phys. Rev. Lett. 101 , 194502 (2008).
- [5] L. Salkin, A. Schmit, P. Panizza, and L. Courbin, Generating soap bubbles by blowing on soap films, Phys. Rev. Lett. 116 , 077801 (2016).

- [6] S. Poulain, E. Villermaux, and L. Bourouiba, Ageing and burst of surface bubbles, J. Fluid. Mech. 851 , 636 (2018).
- [7] C. Isenberg, Soap films and bubbles, Phys. Educ. 16 , 218 (1981).
- [8] G. R ´ amme, Reflected laser light from a soap bubble - a demonstration experiment, Phys. Educ. ' 27 , 282 (1992).
- [9] G. R ´ amme, Surface tension from deflating a soap bubble, Phys. Educ. ' 32 , 191 (1997).
- [10] K. Schilling and M. Zessner, Foam in the aquatic environment, Water Res. 45 , 4355 (2011).
- [11] Toxic foam floods the streets of Bangalore, https://www.cnn.com/2017/05/31/asia/india-toxic-foam-lake/index. html , [Accessed: 11-July-2019].
- [12] Soap Bubble Wiki, https://soapbubble.fandom.com , [Accessed:
- 29-June-2019].
- [13] Physics of Continuous Matter , 2nd ed. (CRC Press, 2011).
- [14] P. G. de Gennes, 'young' soap films, Langmuir 17 , 2416 (2001).
- [15] Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves , 1st ed. (Springer, 2004).
- [16] C. Cohen, B. D. Texier, E. Reyssat, J. H. Snoeijer, D. Quéré, and C. Clanet, On the shape of giant soap bubbles, Proc. Natl. Acad. Sci USA 114 , 2515 (2017).
- [17] D. Mudgil, S. Barak, and B. S. Khatkar, Guar gum: processing, properties and food applications-a review, J. Food Sci. Technol. 51 , 409 (2014).
- [18] J. C. Burton, F. M. Huisman, P. Alison, D. Rogerson, and P. Taborek, Experimental and numerical investigation of the equilibrium geometry of liquid lenses, Langmuir 26 , 15316 (2010).
- [19] E. A. van Nierop, B. Scheid, and H. A. Stone, On the thickness of soap films: an alternative to frankel's law, J. Fluid Mech. 602 , 119 (2008).
- [20] B. H. Zimm, Dynamics of polymer molecules in dilute solution: Viscoelasticity, flow birefringence and dielectric loss, J. Chem. Phys. 24 , 269 (1956).
- [21] C. Wagner, Y. Amarouchene, D. Bonn, and J. Eggers, Droplet detachment and satellite bead formation in viscoelastic fluids, Phys. Rev. Lett. 95 , 164504 (2005).
- [22] C. Clasen, J. P. Plog, W.-M. Kulicke, M. Owens, C. Macosko, L. E. Scriven, M. Verani, and G. H. McKinley, How dilute are dilute solutions in extensional flows?, J. Rheol. 50 , 849 (2006).
- [23] L. Palangetic, N. K. Reddy, S. Srinivasan, R. E. Cohen, G. H. McKinley, and C. Clasen, Dispersity and spinnability: Why highly polydisperse polymer solutions are desirable for electrospinning, Polymer 55 , 4920 (2014).
- [24] F. Del Giudice, S. J. Haward, and A. Q. Shen, Relaxation time of dilute polymer solutions: A microfluidic approach, J. Rheol. 61 , 327 (2017).
- [25] G. H. McKinley and T. Sridhar, Filament-stretching rheometry of complex fluids, Ann. Rev. Fluid Mech. 34 , 375 (2002).
- [26] J. Dinic, L. N. Jimenez, and V. Sharma, Pinch-off dynamics and dripping-onto-substrate (dos) rheometry of complex fluids, Lab Chip 17 , 460 (2017).
- [27] J. Dinic and V. Sharma, Macromolecular relaxation, strain, and extensibility determine elastocapillary thinning and extensional viscosity of polymer solutions, Proc. Nat. Acad. Sci. 116 , 8766 (2019).
- [28] J. C. Burton, J. E. Rutledge, and P. Taborek, Fluid pinch-off in superfluid and normal he , Phys. Rev. E 4 75 , 036311 (2007).
- [29] P. P. Bhat, S. Appathurai, M. T. Harris, M. Pasquali, G. H. McKinley, and O. A. Basaran, Formation of beads-on-a-string structures during break-up of viscoelastic filaments, Nature Phys. 6 , 625 (2010).
- [30] P. de Sainte Claire, Degradation of peo in the solid state: A theoretical kinetic model, Macromolecules 42 , 3469 (2009).
- [31] S. Morlat and J.-L. Gardette, Phototransformation of water-soluble polymers. i: photoand thermooxidation of poly(ethylene oxide) in solid state, Polymer 42 , 6071 (2001).
- [32] A. Dupas, I. Hénaut, J.-F. Argillier, and T. Aubry, Mechanical degradation onset of polyethylene oxide used as a hydrosoluble model polymer for enhanced oil recovery, Oil &amp; Gas Sci. Tech. 67 , 931 (2012).
- [33] B. A. Buchholz, J. M. Zahn, M. Kenward, G. W. Slater, and A. E. Barron, Flow-induced chain scission as a physical route to narrowly distributed, high molar mass polymers, Polymer 45 , 1223 (2004).
- [34] A. J. M ´ uller, J. A. Odell, and S. Carrington, Degradation of semidilute polymer solutions in elongational flows, Polymer ' 33 , 2598 (1992).
- [35] R. Holyst, A. Bielejewska, J. Szymański, A. Wilk, A. Patkowski, J. Gapiński, A. Żywociński, T. Kalwarczyk, E. Kalwarczyk, M. Tabaka, N. Ziębacz, and S. A. Wieczorek, Scaling form of viscosity at all length-scales in poly(ethylene glycol) solutions studied by fluorescence correlation spectroscopy and capillary electrophoresis, Phys. Chem. Chem. Phys. 11 , 9025 (2009).
- [36] K. Devan, , and J. C. Selser, Asymptotic behavior and long-range interactions in aqueous solutions of polyethylene oxide), Macromolecules 24 , 5943 (1991).
- [37] Q. Ying and B. Chu, Overlap concentration of macromolecules in solution, Macromolecules 20 , 362 (1987).
- [38] M. M. Cross, Polymer rheology: Influence of molecular weight and polydispersity, J. Appl. Polym. Sci. 13 , 765 (1969).
- [39] W.-M. Kulicke, M. Elasabee, C. D. Eisenbach, and M. Peuscher, Effect of molecular weight and molecular weight distribution on the rheological properties of aqueous poly(ethylene oxide) solution, Polymer Bulletin 9 , 190 (1983).
- [40] P. K. Bhattacharjee, J. P. Oberhauser, G. H. McKinley, L. G. Leal, and T. Sridhar, Extensional rheometry of entangled solutions, Macromolecules 35 , 10131 (2002).
- [41] X. Ye, R. G. Larson, C. Pattamaprom, and T. Sridhar, Extensional properties of monodisperse and bidisperse polystyrene solutions, J. Rheol. 47 , 443 (2003).
- [42] X. Ye and T. Sridhar, Effects of the polydispersity on rheological properties of entangled polystyrene solutions, Macromolecules 38 , 3442 (2005).
- [43] Soap films: studies of their thinning , 1st ed. (Pergamon Press, 1959).

- [44] L. Champougny, J. Miguet, R. Henaff, F. Restagno, F. Foulogne, and E. Rio, Influence of evaporation on soap film rupture, Langmuir 34 , 3221 (2018).
- [45] L. Saulnier, L. Champougny, G. Bastien, F. Restagno, D. Langevin, and E. Rio, A study of generation and rupture of soap films, Soft Matter 10 , 2899 (2014).
- [46] S. Naire, R. J. Braun, and S. A. Snow, An insoluble surfactant model for a vertical draining free film, J. Colloid. Int. Sci. 230 , 91 (2000).
- [47] S. Berg, E. A. Adelizzi, and S. M. Troian, Experimental study of entrainment and drainage flows in microscale soap films, Langmuir 21 , 3867 (2005).
- [48] L. W. Schwartz and R. V. Roy, Modeling draining flow in mobile and immobile soap films, J. Colloid Interface Sci. 218 , 309 (1999).
- [49] D. Langevin, Influence of interfacial rheology on foam and emulsion properties, Adv. Colloid Int. Sci. 88 , 209 (2000).
- [50] J. Seiwert, B. Dollet, and I. Cantat, Theoretical study of the generation of soap films: role of interfacial visco-elasticity, J. Fluid Mech. 739 , 124 (2014).
- [51] A. A. Sonin, A. Bonfillon, and D. Langevin, Thinning of soap films: The role of surface viscoelasticity, J. Colloid. Int. Sci. 162 , 323 (1994).
- [52] R. Bruinsma, J. di Meglio, D. Quéré, and S. Cohen-Addad, Stabilization of aqueous foam by hydrosoluble polymers. 1. sodium dodecyl sulfate-poly(ethylene oxide) system, Langmuir 8 , 324 (1992).
- [53] E. A. Adelizzi and S. M. Troian, Interfacial slip in entrained soap films containing associating hydrosoluble polymer, Langmuir 20 , 7482 (2004).
- [54] R. Bruinsma, J. di Meglio, D. Quéré, and S. Cohen-Addad, Formation of soap films from polymer solutions, Langmuir 8 , 3161 (1992).
- [55] S. Cohen-Addad and J.-M. di Meglio, Stabilization of aqueous foam by hydrosoluble polymers. 2. role of polymer/surfactant interactions, Langmuir 10 , 773 (1994).
- [56] X. L. Wu, R. Levine, M. Rutgers, H. Kellay, and W. I. Goldburg, Infrared technique for measuring thickness of a flowing soap film, Rev. Sci. Instrum. 72 , 2467 (2001).
- [57] G. Debrégeas, P.-G. de Gennes, and F. Brochard-Wyart, The life and death of 'bare' viscous bubbles, Science 279 , 1704 (1998).
- [58] M. S. Bhamla, C. Chai, M. A. Àlvarez Valenzuela, J. Tajuelo, and G. G. Fuller, Interfacial mechanisms for stability of surfactant-laden films, PLoS One 12 , e0175753 (2017).
- [59] B. H. Cao and M. W. Kim, Molecular weight dependence of the surface tension of aqueous poly(ethylene oxide) solutions, Faraday Discuss. 98 , 245 (1994).
- [60] W. Zhang and R. G. Larson, Tension-induced nematic phase separation in bidisperse homopolymer melts, ACS Cent. Sci. 4 , 1545 (2018).
- [61] S. M. Clegg, P. A. Williams, P. Warren, and I. D. Robb, Phase behavior of polymers with concentrated dispersions of surfactants, Langmuir 10 , 3390 (1994).
- [62] E. Helfand and G. H. Fredreickson, Large fluctuations in polymer solutions under shear, Phys. Rev. Lett. 62 , 2468 (1989).
- [63] S. T. Milner, Hydrodynamics of semidilute polymer solutions, Phys. Rev. Lett. 66 , 1477 (1991).
Stephen Frazier Xinyi Jiang Justin C. Burton
How to make a giant bubble
Using mixtures of soap, water, and long chain polymers, free-floating soap
bubbles can be formed with volumes approaching 100 m$^3$. Here we investigate
how such thin films are created and maintained over time. We show how the
extensional rheology is the most important factor in creating the bubble, and
how polydispersity in molecular weight of the solvated polymers leads to better
performance at lower concentrations. Additionally, using IR absorption, we
measure soap film thickness profiles and film lifetimes. Although the initial
thickness mostly depends on the choice of detergent, polymers can dramatically
increase film lifetime at high molecular weights and high concentrations,
although such high concentrations can inhibit the initial film formation. Thus,
the ideal concentration of polymer additives for making giant bubbles requires
a robust viscoelastic rheology during extension, and is aided by long film
lifetimes during gravitational drainage and evaporation.