Abstract
Active fluid droplets surrounded by oil can spontaneously develop circulatory flows. However, the dynamics of the surrounding oil and their influence on the active fluid remain poorly understood. To investigate interactions between the active fluid and the passive oil across their interface, kinesin-driven microtubule-based active fluid droplets were immersed in oil and compressed into a cylinder-like shape. The droplet geometry supported intradroplet circulatory flows, but the circulation was suppressed when the thickness of the oil layer surrounding the droplet decreased. Experiments with tracers and network structure analyses and continuum models based on the dynamics of self-elongating rods demonstrated that the flow transition resulted from flow coupling across the interface between active fluid and oil, with a millimeter–scale coupling length. In addition, two novel millifluidic devices were developed that could trigger and suppress intradroplet circulatory flows in real time: one by changing the thickness of the surrounding oil layer and the other by locally deforming the droplet. This work highlights the role of interfacial dynamics in the active fluid droplet system and shows that circulatory flows within droplets can be affected by millimeter–scale flow coupling across the interface between the active fluid and the oil.
Similar content being viewed by others
Introduction
Active fluids flow without external energy input owing to force generated by active entities that consume local fuel to generate kinetic energy1,2,3,4,5,6,7,8,9,10. Active fluids can self-organize into circulatory flows that are sensitive to confinement shape11,12,13,14,15,16,17,18,19,20,21,22,23,24,25. However, little is known about the role of boundary conditions in the self-organization of confined active fluids, especially fluid boundaries such as water–oil interfaces. Fluid boundaries are known to induce the coupling of the fluid dynamics on both sides of the boundary owing to hydrodynamic coupling26,27. For example, when active fluid is confined in a droplet immersed in liquid crystal, the liquid crystal develops oscillating rings surrounding the droplet28,29,30,31, which indicates that the active fluid alters the passive fluid configuration on the opposite side of the boundary. However, the principles underlying such coupling and the associated fluid mechanics have not yet been elucidated. Understanding such coupling dynamics is essential to unravel the dynamic role of fluid boundaries at interfaces between passive and active fluids and their impact on active fluid flows. Here, we investigated the hydrodynamic coupling between active and passive fluids in a water-in-oil active fluid droplet system with an approach that combines experiments and modeling. In our experiments, we confined an active fluid in a water-in-oil droplet that was compressed into a cylinder-like shape32 and characterized the flow coupling between active fluid and oil (passive fluid) near water–oil interface. We focused on how active fluids drive the oil through interface and on how the oil configuration can, in return, influence the self-organization of active fluid. We also determine the characteristic length scale of this active–passive fluid coupling. To gain deeper insight into our experimental results, we developed a continuum complex fluid simulation based on established active fluid models26 and explored methods of directing active fluid flows with novel millifluidic devices that can manipulate interfacial dynamics and droplet shapes in real time.
Results and discussion
Kinesin-driven, microtubule-based active fluid
To investigate coupling between active and passive fluids across a fluid boundary, we selected microtubule–kinesin complex active fluid because it is tunable and reproducible and has established characterizations and models to serve as references33,34,35,36,37. Microtubule–kinesin active fluid is driven by extensile microtubule bundles that are assembled by depletion and extended by kinesin motor clusters (Fig. 1a)38. The bundles repeatedly extend, brake, and anneal, resulting in an active gel network whose structure is constantly rearranged32. The rearranging active gel stirs the surrounding liquid (96% water), which causes flows. In this paper, we study the fluid dynamics of this active fluid when confined in a fluid cylindrical boundary.
Confining active fluid in an immobilized active droplet
To create a fluid cylindrical boundary, we confined the active fluid in water-in-oil droplets compressed by plates at the top and bottom. The combination of compression and interfacial tension deformed the droplets into cylinder-like shapes39. However, the compressed droplets were motile26,32,40, so we curved the top plate into a dome-like shape to fix the droplet at the dome center and thus allow the long-term observation of intradroplet flows (Fig. 1b and Supplementary Discussion S1). The flows within the compressed droplets self-organized into droplet-wise circulatory flows (Supplementary Video S1)13 whose circulation depended on the thickness of the oil layer surrounding the droplet (Fig. 1c). The intradroplet circulatory flows were more pronounced with a thicker layer of oil but suppressed with a thinner layer of oil (Supplementary Video S2).
To quantify the coherence of the circulatory flows, we added fluorescent tracer particles to the active fluid. We defined the circulation order parameter (COP) as the average fraction of the azimuthal component in the flow velocity: COP ≡ \(\langle v_{{i,\theta }} /\left| {\user2{v}_{i} } \right|\rangle_{i}\), where vi,θ is the azimuthal component of velocity vi of the ith tracer, and \(~\left\langle {\quad} \right\rangle _{i}\) indicates averaging of the tracer particles11. A COP of 1 indicates perfect counterclockwise circulation, a COP of − 1 indicates perfect clockwise circulation, and a COP of 0 indicates chaotic flow. The COP analyses revealed that, in droplets with a thick oil layer, circulatory flow developed over the first hour and remained at a steady state for several hours. By contrast, in droplets with a thin oil layer the COP fluctuated (Fig. 1d). These results demonstrate that self-organization of active fluid confined in fluid boundaries is sensitive to parameters outside the boundaries.
Intradroplet circulatory flows depended on the geometries of droplets and oil
To investigate how the intradroplet flows were influenced by the shapes of the droplet and oil, we first immersed compressed active droplets with the same height (\(h =\) 1 mm) but different radii (\(r =\) 0.5–5 mm) in oil baths (\(R =\) 5 mm) and measured the COP (Fig. 2a). Our measurements revealed that circulatory flow developed when the droplet radius was smaller than a critical threshold (\(r_{c} \approx\) 2.4 mm); enlarging the droplet above this limit suppressed the circulatory flow (Fig. 2a). Moreover, we found that this critical radius was reduced (\(r_{c} \approx\) 1.4 mm) when the oil bath was smaller (\(R =\) 3.5 mm; Fig. 2b). These results highlight the role of oil configuration in the formation of intradroplet circulatory flows.
Circulatory flows within droplets depended on the oil layer thickness
To characterize the influence of the oil configuration on intradroplet flows, we fixed the droplet shape (\(r \approx\) 2.4 mm, \(h =\) 2 mm), varied the size of the oil bath (R) or the thickness of the oil layer surrounding the droplet (Δ ≈ 0–11 mm) and measured the COP as a function of oil layer thickness (Fig. 3a). Our measurements revealed that the intradroplet flows were sensitive to the oil configuration when the oil layer was thinner than a critical thickness (Δc ≈ 4 mm); oil layers thicker than this limit did not influence intradroplet flows. We hypothesized that the critical thickness of the oil layer originates from active fluid–oil coupling near the interface, which should not be sensitive to droplet shape40. To test this hypothesis, we repeated the measurements with droplet shapes with different r and h values (Fig. 3b–d). The results showed that across our explored droplet geometries, the critical thickness was similar. When the oil layer was thinner than this thickness, the COP varied, whereas when the oil layer was thicker than the critical thickness, the circulatory flows remained intact. The existence of this universal critical thickness implied that the intradroplet active fluid was coupled to the oil within the critical thickness from the interface, which is consistent with our hypothesis.
To investigate how this coupling affected the net flow rates of intradroplet circulation, we analyzed the flow profiles of averaged azimuthal velocities (Fig. 3e–h) and found that the intradroplet flows that had a higher COP flowed more coherently in the azimuthal direction so had a higher net azimuthal flow rate. Conversely, the flows that had a lower COP flowed more chaotically and so had a lower net flow rate. A slower flow rate does not indicate a slower flow speed; rather, it indicates that the flows are more chaotic and thus the positive (counterclockwise) and negative (clockwise) azimuthal flows cancel each other out. The flow speed of intradroplet active fluid was independent of oil layer thickness. Our analysis revealed that thickening the oil layer irregularly supported, suppressed, or even ceased the net azimuthal flows, which agreed with our COP data (Fig. 3a–d). This irregularity revealed a nonlinear influence of the active fluid–oil coupling on circulatory flows. Despite the complex nature of the active fluid–oil coupling, our data demonstrated that the oil configuration could direct the self-organization of intradroplet active fluid through the water–oil interface.
The oil layer thickness influenced microtubule network structure in the droplet
To explore whether the microtubule network structure was affected by the oil configuration, we analyzed confocal microscopy images of the active fluid droplet with the snake algorithm to extract the network structure and reveal the bundle orientational distributions (Fig. 4b–d)41. Our analyses revealed that when the droplet (r ≈ 2.4 mm, ℎ = 1 mm) developed circulatory flows, the microtubule bundles near the water–oil interface tended to align at an angle ~ 15° from parallel with the interface (blue solid curve in Fig. 4d). The alignment decreased with increasing distance from the interface (red solid curve in Fig. 4d), and at the droplet center the microtubule bundles were oriented randomly (green solid curve in Fig. 4d). This indicates that circulatory flows within droplets were accompanied by a thin nematic layer of microtubule bundles near the water–oil interface11. Furthermore, this formation of a nematic layer was suppressed when the thickness of the oil layer was decreased and the intradroplet circulatory flow was suppressed (dashed curves in Fig. 4d). These results suggest that fluid dynamics in the oil penetrate the water–oil interface and influence the self-arrangement of the microtubule network inside the droplet. This is consistent with our observations that circulatory flows are dependent on the thickness of the oil layer (Fig. 3).
The oil layer thickness influenced active stress distribution in the droplet
To gain insight into the impact of oil layer thickness from the perspective of fluid dynamics, we analyzed the time-averaged flow fields, vorticity maps, director fields, nematic order parameter distributions, and active stress distributions in active fluid near the water–oil interface (Fig. 4e–j). We found that, for a droplet immersed in a thicker oil layer (Δ = 2.4 mm), the directors were mostly aligned with the fluid flow and the nematic order parameter decreased with distance from the interface (Fig. 4e, f). This variation in alignment order led to a gradient in active stress near water–oil interface (Fig. 4g). This active stress gradient generated forces that directed the coherent flow (Fig. 4e). Conversely, for the droplet immersed in a thinner oil layer (Δ = 1.1 mm), the directors oriented more chaotically with nearly zero nematic order parameter near the water–oil interface (Fig. 4i). As such, the active stress was uniform (Fig. 4j) and the net flow velocity was nearly zero (Fig. 4h). These results show flow coupling across the interface between active fluid and oil that influences the stress distribution in the active fluid and thus can direct active fluid flows.
Active fluid in the droplet induced chaotic flows in the oil
Our data showed that the thickness of oil layer surrounding the droplet influenced the self-organization of the intradroplet active fluid flows, which suggests an interaction between flows in the droplet and flows in the oil. To reveal such an interaction, we monitored the flows in the oil (Δ ≈ 2.4 mm) as well as in the droplet (\(r \approx\) 2.4 mm, h = 1 mm). To distinguish the flows in both regimes, we doped the oil with 1-µm tracers and doped the droplet with 3-µm tracers (Fig. 5a) and monitored the tracers for 1 h. Time-averaged velocity fields and vorticity maps (Fig. 5b) and flow profiles of azimuthal velocities (Fig. 5c) showed no observable net flows in oil in either circulating or noncirculating droplets. However, the absence of net flow does not necessarily imply that the oil is quiescent, as zero net flows can result from chaotic flows whose velocities canceled out over a time average11. Flow speed profiles revealed that oil near the interface developed 2- to 4-µm/s chaotic flows that decayed with distance from the interface with a decay length of ~ 0.5 mm (Fig. 5d). This decay length suggests that the dynamics of the intradroplet active fluid were coupled to the oil near the interface with a millimeter–scale coupling length. This implies that disturbances to the oil (such as stirring) within this coupling range might influence the intradroplet active fluid flows even if the disturbance does not physically contact the droplet. Conversely, a disturbance outside this coupling range might not affect the intradroplet circulatory flows. This suggestion is consistent with our observation that changing the thickness of oil layers within a critical thickness affects the formation of intradroplet circulatory flows and vice versa (Fig. 3).
These results characterizing flows in the oil surrounding active fluid droplets are consistent with the model prediction by Young et al. that the surrounding oil will remain quiescent when an active droplet develops circulating flows and will be driven to flow when the active fluid flows are noncirculating (such as extensile or quadruple flows)27. Our data showed that when the droplet was in a circulating state (Fig. 5d, blue curve), flows in the oil only developed near the water–oil interface, but when the droplet was in a chaotic state (red curve) the flows in the oil extended to the container surface (\(\rho =\) 1.1 mm). However, Young et al. predicted that intradroplet circulation would be accompanied by a counter-rotation of active fluid near the water–oil interface, which was observed in bacteria-based active droplets14, but our data did not show this counter-rotation (Figs. 3e–h, 5b, c). A possible explanation for this discrepancy is system dimensionality; Young et al.’s model and bacteria-based active droplet system were both two dimensional, whereas our active droplet is a three-dimensional system in which confinement by the ceiling and floor might have induced additional friction that inhibited the development of counter-rotation near the interface.
A continuum simulation qualitatively agreed with experimental outcomes
Our experimental results show that the coupling of flows within and outside droplets influenced the self-organization of flows within the droplet. To gain deeper insight into this flow coupling, we modeled the active fluid droplet system with an existing active fluid model developed by Gao et al. in 201726. We selected Gao et al.’s model because, while active fluids have been modeled using swimmer-based simulations42,43 and continuum equations of mean fields of active particles44,45,46,47,48,49,50, Gao et al.’s framework not only includes multiphase fluids (oil and water) along with associated interfaces, but it also succeeds in describing the self-propelling and self-rotating characteristics of active droplets26. Moreover, Gao et al.’s model shows that active fluid encapsulated in a water-in-oil droplet can induce flows in the surrounding oil, which was observed in our experiments (Fig. 5)26. Therefore we adopted Gao et al.’s model to test its capability to describe our experimental outcomes.
Model description
The model considered two main forces: (1) interfacial tension force from the droplet surface, \(\user2{T} \equiv \left( {\gamma K/\epsilon } \right)\nabla c\), where \(\gamma\) is the water–oil interfacial tension; \(\epsilon\) is the interface thickness; \(K \equiv c\left( {c - 1} \right)\left( {c - 1/2} \right) - \epsilon^{2} \nabla ^{2} c\), the chemical potential that characterizes the phase variation within the interface region; and c is a phase function with \(c = 1\) representing water and \(c = 0\) representing oil51 and (2) active stress exerted by extensile microtubule-based bundles in active fluid that was proportional to the orientational order of bundles, \(\user2{\sigma }_{a} \equiv \alpha c\user2{D}\)52, where \(\alpha\) is an activity coefficient; \(\user2{D} \equiv \mathop \smallint \limits_{\user2{p}}^{~} \user2{pp}\Psi d\user2{p}\), the local nematic order of bundles53; \(\user2{p}\) represents the bundle orientation; and \(\Psi\) represents the probability distribution of bundles that satisfies the Smoluchowski equation44,54:
where \({x}\) represents the center of mass of the bundle and \(\nabla _{\user2{p}} \equiv \partial \left( {\user2{I} - \user2{pp}} \right)/\partial \user2{p}\) is the surface derivative on the unit sphere. These two forces were exerted on incompressible fluids (\(\nabla \cdot \user2{u} = 0\)) to create flows (u) that satisfied the Navier–Stokes equation:
where \(\user2{F} \equiv \user2{T} + \nabla \cdot \user2{\sigma }_{a}\), the net body forces from interfacial tension and extensile bundles, \(\rho\) is the fluid density, p is the fluid pressure, and \(\mu\) is the dynamic viscosity of fluids. The fluids were confined in no-slip boundaries whose geometries were identical to experimental containers and consisted of a circular side wall, a flat floor, and a curved ceiling (Fig. 1b). The boundary was filled with oil within which a compressed active droplet was immersed.
To solve the equations so as to determine the evolution of fluid flows, u, we initialized the flow field as quiescent fluids \({{u=0}}\) under uniform pressure \({{p=0}}\) with uniformly suspended isotropic bundles whose translational and orientational distributions were perturbed with 15 random modes (details of the random modes are provided in Supplementary Discussion S2)44,55. Then we evolved the fluid flows for 3 h with the finite element method based on the computational fluid dynamics software COMSOL Multiphysics™56,57. We made the assumption that the geometry of the interface remained invariant over time and water–oil interfacial fluctuation was negligible, based on our experimental results (Supplementary Video S1) and previous studies32 showing that the water–oil interfacial tension is strong enough that the geometry of the droplet interface remains nearly unchanged over time. Thus, the phase function is independent of time:
where \({{H}}\) is the Heaviside step function. Details of the model, including the chosen values of parameters and explicit forms of each equation in three-dimensional components, are provided in Supplementary Discussion S2.
Comparison between model prediction and experimental measurements
To test the model’s ability to describe our experimental system, we arranged two simulation systems with identical droplets (\(r =\) 2.4 mm, \(h =\) 2 mm) immersed in oil layers of different thicknesses (Δ = 1.1 and 2.6 mm). The simulation predicted that the droplet immersed in the thicker oil layer would develop a steady intradroplet circulatory flow, whereas the droplet immersed in the thinner oil layer would have chaotic flows (Fig. 6a, b inset). We then systematically varied the oil layer thickness (Δ = 0–9.6 mm) in the model while maintaining the droplet geometry and analyzed the time-averaged COP within each droplet (Fig. 6b). The COP was sensitive to the oil layer thickness when the layer was thinner than ~ 2.2 mm, which suggests that the flows within and outside of the droplets were coupled. To reveal such coupling, we analyzed the flow profiles of azimuthal velocities across the water–oil interface (Fig. 6c), which showed that circulatory flows within the droplet induced a thin layer of circulatory flow in the oil with a layer thickness of 0.3 to 2 mm, whereas chaotic flows in the droplet did not induce net flows in the oil. However, analysis of the flow speed profiles revealed that oil near the interface developed flows with a thickness of ~ 1 mm (magenta curve in Fig. 6c inset), suggesting that flows within the active fluid droplet induce flows in the oil near the water–oil interface regardless of flowing state of the intradroplet active fluid (circulating or noncirculating).
Characterization of active fluid–oil flow coupling with a cross-correlation function
To gain deeper insight into coupling between the active fluid and the oil, we analyzed how the flows of active fluid near the interface (100 µm from the interface), \(\user2{v}_{w}\), were correlated to oil flows, \(\user2{v}_{o}\), by calculating the normalized same-time spatial cross-correlation function between \(\user2{v}_{w}\) and \(\user2{v}_{o}\):
where \( {\Delta \user2{X}}\) represents the separation between a pair of active fluid and oil elements and \(\left\langle {\quad} \right\rangle _{{\user2{x},t}}\) indicates averaging over time in the active fluid region within 100 µm of the interface (\(\left| {r - x} \right| \le 100\) µm). To minimize the influence of the top and bottom boundaries12,58,59,60 in our correlation analysis, we only considered the flows at the midplane (\(z = h/2\)). To reveal the coupling range between flows of active fluid and oil, we averaged the correlation function over the orientation:
Our analysis revealed that the correlation function decayed nearly exponentially with increasing distance between the active fluid element and the oil element (\(C\sim e^{{ - \Delta X/L}}\), where L is the correlation length), which suggests that the active fluid–oil interaction is short ranged (Fig. 7a). To quantify the interaction range, we extracted the correlation length, L, and then analyzed the correlation length as a function of oil layer thickness (Fig. 7b). Our analysis revealed that the correlation length linearly increased with the oil layer thickness (\(L \approx \Delta\)), eventually reaching saturation (\(L \approx 1.4\) mm). The linear increase indicated that the active fluid flows were coupled to the oil flows throughout the oil region. Hence, the geometry of the oil (such as oil layer thickness) affected the active fluid flows, which is consistent with our observation that the COP in the active fluid changed rapidly with the oil layer thickness when the thickness was small (\(\Delta < 2.2\) mm in Fig. 6b). The saturation indicated that the range of active fluid–oil interaction has an upper limit (~ 1.4 mm) above which the motion of oil elements did not affect the active fluid. This is consistent with our observation that the COP in the active fluid was independent of the oil layer thickness when the oil layer was sufficiently thick (\(\Delta > 2.2\) mm in Fig. 6b). Moreover, the scale of the analyzed correlation length (~ 1.4 mm; Fig. 7b) was consistent with the observed coupling length (~ 1 mm; Fig. 6c inset). This consistency, along with simulated shear stress analysis (Supplementary Discussion S3), supports the assertion that active fluid and oil interact across the water–oil interface with a millimeter–scale interaction range. Further analyses of the role of interfacial properties, such as viscosity contrast and interfacial tension, on the simulated results are provided in Supplementary Discussion S4. Overall, the simulation outcomes were qualitatively consistent with our experimental characterizations of the coupling of intradroplet and extradroplet flows (Figs. 3, 5). This consistency demonstrates that the model can describe active fluid–oil coupling and how the coupling influences the flows inside the active fluid droplet.
Limitations of the model
The simulation failed to match two of our experimental outcomes: First, the simulation showed that the intradroplet fluid flows became chaotic when oil layers were thicker than ~ 5.6 mm (Fig. 6b), whereas in our experimental data, the intradroplet fluid flows were insensitive to the oil arrangement when the oil layers were thicker than ~ 4 mm (Fig. 3a–d). Second, the simulation predicted the induction of circulatory flows in oil driven by intradroplet circulation (Fig. 6c), whereas in the experiments, oil did not develop net flows regardless of how the fluid flowed within the droplets (Fig. 5c). It is possible that these discrepancies could be mitigated by allowing the interface to deform following the rules of spontaneous phase separation of oil and water (the Cahn–Hilliard model)61.
Intradroplet circulatory flows were triggered and suppressed in real time with novel millifluidic devices
Our experimental results show that the formation of intradroplet circulatory flows depends on droplet geometry (Fig. 2) and oil layer thickness (Fig. 3). This suggests that millifluidic devices could control the formation and suppression of circulatory flows in real time by changing the droplet shape or oil layer thickness. We demonstrated this by designing and testing two such devices.
Device that changes oil layer thickness with a movable wall
To manually tune the oil layer thickness, we developed a cylindrical container that compressed an active droplet (ℎ = 2 mm, r ≈ 2.4 mm) and had one movable wall (pink blade in Fig. 8a, b) that could be moved toward the droplet to reduce the thickness of the oil layer adjacent to one part of the droplet (minimum oil layer thickness \(\Delta _{m} \approx\) 1.2 mm; Fig. 8b left). When the wall was near the droplet, the active fluid flowed chaotically (\(\left| {{\text{COP}}} \right|\) ≲ 0.2; Fig. 8c). When we increased the oil thickness by moving the wall away from the droplet (\(\Delta _{m} \approx\) 2.6 mm; Fig. 8b middle), the active fluid developed circulatory flows in ~ 30 min (COP = 0.4–0.6; Fig. 8c). The circulatory flows lasted for ~ 1 h and then transitioned to chaotic flows after we moved the wall back toward the droplet (\(\Delta _{m} \approx\) 1.2 mm; Fig. 8b right, Supplementary Video S3). These results show that it is possible to develop and inhibit intradroplet circulatory flows locally in real time without physically contacting the droplet.
Device that deforms droplet with a movable ceiling
To manually shape the droplet, we compressed the droplet in a cylindrical container with a movable ceiling (Fig. 8d). We first compressed the droplet to a short cylinder-like shape (\(r \approx\) 2.0 mm, \(h =\) 2 mm), which supported the formation of circulatory flows (COP = 0.4–0.6, Fig. 8e). The circulatory flow persisted for ~ 40 min before it was manually suppressed (\(\left| {{\text{COP}}} \right|\) ≲ 0.2) by lifting the ceiling (\(h =\) 3 mm), which shaped the droplet into a taller cylinder-like shape (\(r \approx\) 1.7 mm, \(h =\) 3 mm). Conversely, intradroplet circulatory flow could be manually triggered by deforming the droplet from a taller (\(r \approx\) 1.7 mm, \(h =\) 3 mm) to a shorter (\(r \approx\) 2.0 mm, \(h =\) 2 mm) cylinder-like shape (Fig. 8f, Supplementary Video S4). These results demonstrate that manually shaping the droplet can turn intradroplet circulatory flows on and off locally. These findings pave the way for designing fluidic devices that can shape deformable boundaries to direct the self-organization of confined active fluids in real time.
Limitations of the study
This study focused on microtubule-based active fluid and thus the results may not be generalizable to other active fluids. Another limitation of this study is that both the models and experiments neglected interfacial fluctuation, because interfacial tension in this system was strong enough to inhibit interfacial fluctuation. Low interfacial tension can distort interfaces62,63 which could alter flow coupling across the interface. Future research could investigate flow coupling in systems with different interfacial tensions by incorporating vesicles62 or varying interfacial surfactant concentrations64.
Conclusions
This work demonstrates that the self-organization of a water-in-oil droplet of microtubule-based active fluid is influenced by flow coupling across the water–oil interface with a millimeter–scale coupling length. Our experimental data and simulation results indicate that active fluid within the droplet can induce flows in the oil within this coupling range and that disturbances outside of the droplet, such as reducing the thickness of the oil layer to be less than the coupling range, can influence the microtubule network structure and active stress distribution inside the droplet and thus impact the intradroplet flows. While previous studies have shown that the formation of circulatory flows depends on the confinement boundary geometry11,13,14,15,18, this work is the first to highlight the role of boundary conditions on the formation of circulatory flows—specifically the role of fluid boundaries and hydrodynamic coupling across active–passive boundaries in confined active fluid systems.
We also developed two millifluidic devices that can trigger and suppress intradroplet circulatory flows in real time: one suppresses intradroplet circulatory flows without contacting the droplet by manipulating the oil layer thickness and thus disturbing the active fluid–oil coupling and the other compresses droplets to desired height and radius combinations that support or suppress intradroplet circulatory flows. These novel devices provide the biology community with in vitro model systems to probe how the deformation of cell membranes or the disturbance of fluids around cells influences intracellular activities. These approaches could also be used in the development of treatment modalities for cells that are sensitive to biomedical approaches and can only be treated with physical methods, such as deforming cells and imposing shear flows65,66. In mechanical engineering, these systems pave the way for designing machines driven by active fluid with adjustable power output67.
Methods
Fabricate a millifluidic device to confine active water-in-oil droplets
To confine the active fluid in a cylinder-like water-in-oil droplet, we designed a millifluidic device that compressed the droplet between a pair of plates separated by a height (ℎ) of 1 to 2 mm. Because the compressed droplet was self-propelling26,32,40, we immobilized the droplet by curving the upper plate surface into a half-oblate spheroidal dome with a shallow cylindrical well (height 0.2 mm; radius 1 mm) at the dome center (Supplementary Fig. S1a). The curved surface and shallow well immobilized the droplet without significantly impacting the intradroplet active fluid behaviors (Supplementary Discussion S1). Finally, to load the oil and active fluid into the millifluidic device, we drilled a 3.6-mm-long loading channel 2 mm wide and 1.7 mm high, and joined the channel to a 2.4-mm-long neck that had the same width (2 mm) but shorter height (1.4 mm). To fabricate the device, we sketched the 3D design in SolidWorks, programmed the corresponding tool paths in Esprit, and used the tool paths to end-mill a 610 × 38 × 6.4 mm3 acrylic rectangular bar with computer numerical control (McMaster 1227T222). The milled chip was then cleaned with sequential 10-min sonication in detergent (Sigma-Aldrich Z805939), ethanol, and 100 mM potassium hydroxide solution and glued to a fluorophilically-treated glass slide (VWR GWBJ17) with epoxy (Bob Smith Industries BSI-201) to complete the device fabrication68.
Prepare the compressed water-in-oil droplet
We prepared microtubule–kinesin active fluid according to our previous protocols32,69 and pipetted the active fluid through the loading channel into the chamber of the fabricated millifluidic device that was filled with oil (hydrofluoroether, 3 M Novec 7500; Supplementary Fig. S1a). To prevent the microtubule and kinesin proteins from contacting the oil (and thus denaturing), we doped the oil with 1.8% surfactant (perfluoropolyether–polyethylene glycol–perfluoropolyether, RAN Biotechnologies 008-FluoroSurfactant)68 to stabilize the protein near the water–oil interface. (Previous studies of systems with microtubule-based active fluid interfacing with oil showed that microtubules could be centrifuged to the water–oil interface and would then form a two-dimensional active nematic layer18,30,32,62. In our system, we did not centrifuge the samples, and though a small portion of microtubules were spontaneously deposited onto the interface and formed 2D active nematics [Supplementary Discussion S5], the majority of microtubules remained in the bulk and induced active fluid flows [Supplementary Video S1].) The active fluid loaded into the channel then formed a water-in-oil droplet that was compressed between the ceiling and floor of the chamber. The compression deformed the droplet into a cylinder-like shape39 whose height (h) depended on the ceiling–floor separation of the chamber and whose radius (r) depended on the pipetted fluid volume. After injecting the active droplet, we sealed the channel with epoxy. However, after the channel was sealed, air bubbles sometimes formed in the sample and affect the experimental outcomes. To keep the bubbles away from the droplet, we tilted the sample to direct the air bubbles to exit the chamber into the loading channel through the neck. The neck’s smaller opening prevented the bubbles from reentering the chamber (Supplementary Fig. S1a close-up).
In our experiments, we first varied the droplet radius in millifluidic devices with chambers of radii R = 3.5 and 5 mm (Fig. 2) to examine the role of the droplet radius in the formation of the intradroplet circulatory flows. Each of the devices had a half-spheroidal dome with a vertical semi-axis of Ra = 0.25 mm and a horizontal semi-axis matching the chamber radius, Rb = R (Supplementary Fig. S1a). We then immersed the droplets in oil layers of various thicknesses (Fig. 3). The oil was contained in the chamber whose radius was the sum of the droplet radius and oil layer thickness R = r + Δ (Fig. 1b). To minimize the influence of the ceiling shape on the experimental outcomes, we chose a fixed half-spheroidal ceiling (Ra = 0.5 mm, Rb = 5 mm) to cover the chamber. When the chamber was smaller than the ceiling (2R < 2Rb), the ceiling was trimmed to fit into the chamber; conversely, when the chamber was larger than the ceiling (2R > 2Rb), the ceiling was extended horizontally to match the chamber size (Supplementary Fig. S1b&c).
Image and analyze flows
To observe the flows of the active fluid and oil, we doped the active fluid with 0.0004% Alexa 488-labeled 3-µm tracer particles (Polysciences 18861-1), which could be imaged with a green fluorescent protein filter cube (excitation: 440–466 nm, emission: 525–550 nm, Semrock 96372). To characterize the flow behaviors, we imaged the tracers for 3 h and tracked their trajectories \(\user2{x}_{i} \left( t \right)\) with the Lagrangian algorithm70. The trajectories revealed the evolution of the particles’ velocities, \(\user2{v}_{i} \left( t \right) \equiv d\user2{x}_{i} \left( t \right)/dt\), and allowed us to calculate the COP, COP(t) ≡ \(\langle v_{{i,\theta }} \left( t \right)/\left| {\user2{v}_{i} \left( t \right)} \right|\rangle_{i}\), where vi,θ is the azimuthal component of velocity vi of the ith tracer and \(\left\langle {\quad } \right\rangle _{i}\) indicates averaging of the tracer particles (Fig. 1d). To quantify the coherence of the intradroplet circulatory flows, we measured the time-averaged COP: \(\langle {\text{COP}}\left( t \right)\rangle_{t}\) (Figs. 2, 3a–d). To reveal flow rates of intradroplet circulatory flows, we measured the flow profiles of the azimuthal velocities across the droplet, \(\langle v_{\theta } \left( d \right)\rangle \equiv \langle v_{{i,\theta }} \left( {d,t} \right)\rangle _{{t,i}}\), where d is the distance from the water–oil interface (Fig. 3e–h). To reveal the structure of the flows in the droplet, we measured the time-averaged, normalized velocity fields, \(\user2{V}\left( \user2{x} \right) \equiv \left\langle \frac{{\user2{v}\left( {\user2{x},t} \right)}}{{\langle\left| {\user2{v}\left( {\user2{x},t} \right)} \right|\rangle_{\user2{x}} }}\right\rangle_{t}\), and vorticity distributions, \(\Omega \left( \user2{x} \right) \equiv \left\langle {\frac{{\omega \left( {\user2{x},t} \right)}}{{3\sigma \left( t \right)}}} \right\rangle _{t}\), where \(\user2{v}\left( {\user2{x},t} \right)\) is the velocity field calculated from the sequential tracer images with the particle image velocimetry algorithm71; \(\omega \left( {\user2{x},t} \right) \equiv \left[ {\nabla \times \user2{v}\left( {\user2{x},t} \right)} \right]_{z}\), the corresponding vorticity distribution; and \(~\sigma \left( t \right) \equiv {\text{std}}\left[ {\omega \left( {\user2{x},t} \right)} \right]\), the standard deviation of the vorticity (Fig. 1c)11.
To characterize flows in oil, we doped the oil with 0.002% Alexa 488-labeled 1-µm tracer particles (Polysciences 18860-1) and conducted the same observations and analyses for active fluid flows (Fig. 5b, c). To reveal flow activities near water–oil interfaces, we measured the flow speed profiles \(\langle\left| {\user2{v}\left( \rho \right)} \right|\rangle \equiv \langle\left| {\user2{v}_{i} \left( {\rho ,t} \right)} \right|\rangle_{{t,i}}\), where \(\rho\) is the radial coordinate relative to the water–oil interface (Fig. 5c, d).
Image and analyze microtubule network structure
To characterize the influence of the oil layer thickness on the intradroplet microtubule network structure, we imaged the microtubules at the droplet midplane with confocal microscopy (Leica SP5 point scanning confocal microscope). The microtubules were labeled with Alexa 647 (according to our previous protocol69), excited with a 633-nm helium–neon laser, and observed in a 633–647 nm window. To reveal the microtubule arrangement in a circulating active fluid, we used our data to select the droplet shape (r ≈ 2.4 mm, ℎ = 1 mm) and oil layer thickness (\(\Delta \approx\) 2.4 mm) that supported the formation of intradroplet circulatory flows (Fig. 3b). To observe the time-averaged network structure, we imaged the network every 2 s for 1 h and then analyzed the images with the snake algorithm to extract the network structure, which consisted of unit-length segments (Fig. 4b, c)41. We stacked the segment orientations from each image and then counted these orientations to reveal the orientational distribution of the microtubule bundles (Fig. 4d). The microtubule bundle orientations were measured near the water–oil interface, ~ 400 µm from the interface, and at the droplet center (Fig. 4a) and the measurements were repeated on another sample with the same droplet shape but different oil layer thickness (Δ ≈ 1.1 mm) where intradroplet circulation was suppressed.
Analyze microtubule network dynamics
To further characterize the influence of oil layer thickness on the dynamics of the microtubule network, we first measured time-averaged velocity fields and vorticity maps by analyzing the sequential confocal images of the microtubule network near water–oil interfaces with the particle image velocimetry algorithm (Fig. 4e, h)71. To reveal the corresponding bundle alignments in these flows, we measured the time-averaged director fields \(\user2{\hat{n}}\) and nematic order parameter maps (Fig. 4f, i) by first analyzing the bundle orientational tensor, \(\user2{D} \equiv \langle \user2{pp}\rangle_{t}\), where p represents the extracted bundle orientation from the confocal images (Fig. 4c) and \(\left\langle {\quad} \right\rangle _{t}\) represents averaging over time. Then, we determined the nematic order tensor, \(\user2{Q} \equiv \user2{D} - \user2{I}/2\), calculated the maximum eigenvalue, \(\lambda _{{\text{m}}}\), and determined the nematic order parameter as NOP \(= 2\lambda _{{\text{m}}}\) and the director \(\user2{\hat{n}}\) as the corresponding eigenvector53. To show the dynamics that resulted from these bundle configurations, we measured the time-averaged active stress maps (Fig. 4g, j) by calculating the magnitude of the bundle orientational tensor, \(\sqrt {\user2{D}:\user2{D}}\). \(\sqrt {\user2{D}:\user2{D}}\) represents active stress because, according to Gao et al.’s model26, active stress tensor is defined as \(\user2{\sigma }_{a} \equiv \alpha \user2{D}\) where \(\alpha\) is an activity coefficient, which is a constant in our active fluid, and active stress magnitude is determined as \(\sqrt {\user2{\sigma }_{a} :\user2{\sigma }_{a} } = \alpha \sqrt {\user2{D}:\user2{D}}\), which is proportional to \(\sqrt {\user2{D}:\user2{D}}\).
Design fluidic device with movable wall
To change the thickness of oil surrounding a droplet in real time, we designed a fluidic device with a movable wall, inspired by the mechanical iris that is used to adjust the aperture size of optical devices72. However, a conventional mechanical iris consists of at least six blades, and fabricating a six-blade mechanical iris at the micron scale was challenging because at this scale the blades bent spontaneously in our fabrication process and could not be assembled. Therefore, we simplified the design to contain only one blade that was thick enough (1 mm) to remain flat (pink components in Supplementary Fig. S8). The single blade was pinned to a cylindrical container where the pin was constrained in a radial groove in the container and would move toward the chamber center when the blade was rotated counterclockwise. This blade rotation caused the midpoint of the blade to move 1.4 mm closer to the chamber center, thus shrinking the size of the chamber (Supplementary Fig. S8, right). Conversely, when the blade was rotated clockwise, the chamber expanded (Supplementary Fig. S8, left). We enclosed the container by gluing fluorophilic glass to the bottom of the chamber and a curved ceiling to the top of the container. Once the glue was cured, oil and active fluid were pipetted into the chamber via the loading channel and the channel was then sealed with epoxy. While the channel was sealed, the chamber remained semi-open because the movable blade required a gap between the ceiling and floor. This gap would allow oil to evaporate and create air bubbles that could influence experimental outcomes. To minimize the influence of oil evaporation, we overflowed the gap and blade with oil so that the blade was below the oil surface. This arrangement allowed us to rotate the blade without exposing the chamber to air.
Design fluidic device with movable ceiling
To control the shape of the droplet in real time, we designed a millifluidic device whose ceiling could be tuned manually. To fabricate the device, we separated the ceiling from the rest of the device and attached a handle to manually move the ceiling vertically. We also attached a hanger to hold the ceiling on top of the device chamber (Supplementary Fig. S9a). To control the vertical position of the ceiling, we designed three platforms with different heights that were placed outside the edge of the device chamber where we could hang the ceiling to adjust its height from ℎ = 1 to 3 mm (Supplementary Fig. S9b). To minimize the influence of oil evaporation, we immersed the system (including the device chamber and ceiling) in an oil bath enclosed in a Petri dish.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
References
Palacci, J., Sacanna, S., Steinberg, A. P., Pine, D. J. & Chaikin, P. M. Living crystals of light-activated colloidal surfers. Science 339, 936–940 (2013).
Duclos, G. et al. Spontaneous shear flow in confined cellular nematics. Nat. Phys. 14, 728–732 (2018).
Herminghaus, S. et al. Interfacial mechanisms in active emulsions. Soft Matter 10, 7008–7022 (2014).
Schaller, V., Weber, C., Semmrich, C., Frey, E. & Bausch, A. R. Polar patterns of driven filaments. Nature 467, 73–77 (2010).
Farhadi, L., Fermino Do Rosario, C., Debold, E. P., Baskaran, A. & Ross, J. L. Active self-organization of actin-microtubule composite self-propelled rods. Fronti. Phys. 6, 75 (2018).
Yan, J. et al. Reconfiguring active particles by electrostatic imbalance. Nat Mater 15, 1095–1099 (2016).
Driscoll, M. et al. Unstable fronts and motile structures formed by microrollers. Nat. Phys. 13, 375 (2017).
Morin, A., Desreumaux, N., Caussin, J.-B. & Bartolo, D. Distortion and destruction of colloidal flocks in disordered environments. Nat Phys 13, 63–67 (2017).
Soni, V. et al. The odd free surface flows of a colloidal chiral fluid. Nat. Phys. 15, 1188 (2019).
Blair, D. L. & Kudrolli, A. Collision statistics of driven granular materials. Phys. Rev. E 67, 041301 (2003).
Wu, K.-T. et al. Transition from turbulent to coherent flows in confined three-dimensional active fluids. Science 355, eaal1979 (2017).
Fan, Y., Wu, K.-T., Aghvami, A., Fraden, S. & Breuer, K. Effects of confinement on the dynamics and correlation scales in active fluids. arXiv:2103.06334 [cond-mat.soft] (2021).
Suzuki, K., Miyazaki, M., Takagi, J., Itabashi, T. & Ishiwata, S. i. Spatial confinement of active microtubule networks induces large-scale rotational cytoplasmic flow. Proc. Natl. Acad. Sci. 114, 2922–2927 (2017).
Wioland, H., Woodhouse, F. G., Dunkel, J., Kessler, J. O. & Goldstein, R. E. Confinement stabilizes a bacterial suspension into a spiral vortex. Phys. Rev. Lett. 110, 268102 (2013).
Wioland, H., Lushi, E. & Goldstein, R. E. Directed collective motion of bacteria under channel confinement. N. J. Phys. 18, 075002 (2016).
Doxzen, K. et al. Guidance of collective cell migration by substrate geometry. Integr. Biol. 5, 1026–1035 (2013).
Norton, M. M. et al. Insensitivity of active nematic liquid crystal dynamics to topological constraints. Phys. Rev. E 97, 012702 (2018).
Opathalage, A. et al. Self-organized dynamics and the transition to turbulence of confined active nematics. Proc. Natl. Acad. Sci. 116, 4788–4797 (2019).
Bricard, A. et al. Emergent vortices in populations of colloidal rollers. Nat. Commun. 6, 7470 (2015).
Morin, A. & Bartolo, D. Flowing active liquids in a pipe: Hysteretic response of polar flocks to external fields. Phys. Rev. X 8, 021037 (2018).
Delphine, G., David, M., Julien, T. & Bartolo, D. Freezing a flock: Motility-induced phase separation in polar active liquids. Phys. Rev. X 9, 031043 (2019).
Theillard, M. & Saintillan, D. Computational mean-field modeling of confined active fluids. J. Comput. Phys. 397, 108841 (2019).
Narayan, V., Ramaswamy, S. & Menon, N. Long-lived giant number fluctuations in a swarming granular Nematic. Science 317, 105–108 (2007).
Buhl, J. et al. From disorder to order in marching locusts. Science 312, 1402–1406 (2006).
Chen, S., Gao, P. & Gao, T. Dynamics and structure of an apolar active suspension in an annulus. J. Fluid Mech. 835, 393–405 (2018).
Gao, T. & Li, Z. Self-driven droplet powered by active nematics. Phys. Rev. Lett. 119, 108002 (2017).
Young, Y.-N., Shelley, M. J. & Stein, D. B. The many behaviors of deformable active droplets. Math. Biosci. Eng. 18, 2849–2881 (2021).
Hardoüin, J., Guillamat, P., Sagués, F. & Ignés-Mullol, J. Dynamics of ring disclinations driven by active nematic shells. Front. Phys. 7, 169 (2019).
Guillamat, P., Hardoüin, J., Prat, B. M., Ignés-Mullol, J. & Sagués, F. Control of active turbulence through addressable soft interfaces. J. Phys.: Condens. Matter 29, 504003 (2017).
Guillamat, P., Ignés-Mullol, J. & Sagués, F. Taming active turbulence with patterned soft interfaces. Nat. Commun. 8, 564 (2017).
Guillamat, P., Ignés-Mullol, J., Shankar, S., Marchetti, M. C. & Sagués, F. Probing the shear viscosity of an active nematic film. Phys. Rev. E 94, 060602 (2016).
Sanchez, T., Chen, D. T. N., DeCamp, S. J., Heymann, M. & Dogic, Z. Spontaneous motion in hierarchically assembled active matter. Nature 491, 431–434 (2012).
Henkin, G., DeCamp, S. J., Chen, D. T. N., Sanchez, T. & Dogic, Z. Tunable dynamics of microtubule-based active isotropic gels. Phil. Trans. R. Soc. A. 372, 20140142 (2014).
Chandrakar, P. et al. Microtubule-based active fluids with improved lifetime, temporal stability and miscibility with passive soft materials. arXiv:1811.05026 [cond-mat.soft] (2018).
Lemma, L. M., Decamp, S. J., You, Z., Giomi, L. & Dogic, Z. Statistical properties of autonomous flows in 2D active nematics. Soft Matter 15, 3264 (2019).
Bate, T. E., Jarvis, E. J., Varney, M. E. & Wu, K.-T. Collective dynamics of microtubule-based 3D active fluids from single microtubules. Soft Matter 15, 5006–5016 (2019).
Putzig, E., Redner, G. S., Baskaran, A. & Baskaran, A. Instabilities, defects, and defect ordering in an overdamped active nematic. Soft Matter 12, 3854–3859 (2016).
Hilitski, F. et al. Measuring cohesion between macromolecular filaments one pair at a time: depletion-induced microtubule bundling. Phys. Rev. Lett. 114, 138102 (2015).
Alabuzhev, A. A. & Lyubimov, D. V. Effect of the contact-line dynamics on the oscillations of a compressed droplet. J. Appl. Mech. Tech. Phys. 53, 9–19 (2012).
Tjhung, E., Marenduzzo, D. & Cates, M. E. Spontaneous symmetry breaking in active droplets provides a generic route to motility. Proc. Natl. Acad. Sci. 109, 12381–12386 (2012).
Xu, T. et al. SOAX: A software for quantification of 3D biopolymer networks. Sci. Rep. 5, 9081 (2015).
Vicsek, T., Czirók, A., Ben-Jacob, E., Cohen, I. & Shochet, O. Novel type of phase transition in a system of self-driven particles. Phys. Rev. Lett. 75, 1226–1229 (1995).
Hernandez-Ortiz, J. P., Stoltz, C. G. & Graham, M. D. Transport and collective dynamics in suspensions of confined swimming particles. Phys. Rev. Lett. 95, 204501 (2005).
Saintillan, D. & Shelley, M. J. Instabilities, pattern formation, and mixing in active suspensions. Phys. Fluids 20, 123304 (2008).
Dunkel, J. et al. Fluid dynamics of bacterial turbulence. Phys. Rev. Lett. 110, 228102 (2013).
Bratanov, V., Jenko, F. & Frey, E. New class of turbulence in active fluids. Proc. Natl. Acad. Sci. 112, 15048–15053 (2015).
Toner, J., Tu, Y. & Ramaswamy, S. Hydrodynamics and phases of flocks. Ann. Phys. 318, 170–244 (2005).
Thampi, S. P., Golestanian, R. & Yeomans, J. M. Vorticity, defects and correlations in active turbulence. Philos. Trans. R. Soc. A: Math. Phys. and Eng. Sci. 372, 20130366 (2014).
Giomi, L. Geometry and topology of turbulence in active nematics. Phys. Rev. X 5, 031003 (2015).
Urzay, J., Doostmohammadi, A. & Yeomans, J. M. Multi-scale statistics of turbulence motorized by active matter. J. Fluid Mech. 822, 762–773 (2017).
Anderson, D. M., McFadden, G. B. & Wheeler, A. A. Diffuse-interface methods in fluid mechanics. Annu. Rev. Fluid Mech. 30, 139–165 (1998).
Aditi Simha, R. & Ramaswamy, S. Hydrodynamic fluctuations and instabilities in ordered suspensions of self-propelled particles. Phys. Rev. Lett. 89, 058101 (2002).
Chaikin, P. M. & Lubensky, T. C. Principles of Condensed Matter Physics. (Cambridge University Press, 2000).
Gao, T., Blackwell, R., Glaser, M. A., Betterton, M. D. & Shelley, M. J. Multiscale polar theory of microtubule and motor-protein assemblies. Phys. Rev. Lett. 114, 048101 (2015).
Saintillan, D. & Shelley, M. J. Instabilities and pattern formation in active particle suspensions: kinetic theory and continuum simulations. Phys. Rev. Lett. 100, 178103 (2008).
Mittal, S. & Kumar, B. A stabilized finite element method for global analysis of convective instabilities in nonparallel flows. Phys. Fluids 19, 088105 (2007).
Bendaraa, A., Charafi, M. M. & Hasnaoui, A. Numerical study of natural convection in a differentially heated square cavity filled with nanofluid in the presence of fins attached to walls in different locations. Phys. Fluids 31, 052003 (2019).
Chandrakar, P. et al. Confinement controls the bend instability of three-dimensional active liquid crystals. Phys. Rev. Lett. 125, 257801 (2020).
Varghese, M., Baskaran, A., Hagan, M. F. & Baskaran, A. Confinement-induced self-pumping in 3D active fluids. Phys. Revi. Lett. 125, 268003 (2020).
Chandragiri, S., Doostmohammadi, A., Yeomans, J. M. & Thampi, S. P. Flow states and transitions of an active nematic in a three-dimensional channel. Phys. Rev. Lett. 125, 148002 (2020).
Cahn, J. W. & Hilliard, J. E. Free energy of a nonuniform system. I. Interfacial free energy. J. Chem. Phys. 28, 258–267 (1958).
Keber, F. C. et al. Topology and dynamics of active nematic vesicles. Science 345, 1135–1139 (2014).
Soni, H., Luo, W., Pelcovits, R. A. & Powers, T. R. Stability of the interface of an isotropic active fluid. Soft Matter 15, 6318–6330 (2019).
Li, X. & Pozrikidis, C. The effect of surfactants on drop deformation and on the rheology of dilute emulsions in stokes flow. J. Fluid Mech. 341, 165–194 (1997).
Jaros, D. et al. Shear treatment of starter culture medium improves separation behavior of streptococcus thermophilus cells. Eng. Life Sci. 18, 62–69 (2018).
Diz-Muñoz, A., Fletcher, D. A. & Weiner, O. D. Use the force: membrane tension as an organizer of cell shape and motility. Trends Cell Biol. 23, 47–53 (2013).
Needleman, D. & Dogic, Z. Active matter at the interface between materials science and cell biology. Nat. Rev. Mater. 2, 17048 (2017).
DeCamp, S. J., Redner, G. S., Baskaran, A., Hagan, M. F. & Dogic, Z. Orientational order of motile defects in active nematics. Nat. Mater. 14, 1110–1115 (2015).
Bate, T. E., Jarvis, E. J., Varney, M. E. & Wu, K.-T. Controlling flow speeds of microtubule-based 3D active fluids using temperature. J. Vis. Exp. (153), e60484 (2019).
Ouellette, N. T., Xu, H. & Bodenschatz, E. A quantitative study of three-dimensional Lagrangian particle tracking algorithms. Exp. Fluids 40, 301–313 (2005).
Thielicke, W. & Stamhuis, E. J. PIVlab–Towards user-friendly, affordable and accurate digital particle image velocimetry in MATLAB. J. Open Res. Softw. 2, e30 (2014).
Wang, Z., Ding, H., Lu, G. & Bi, X. Use of a mechanical iris-based fiber optic probe for spatially offset Raman spectroscopy. Opt. Lett. 39, 3790–3793 (2014).
Acknowledgements
We thank Dr. Zvonimir Dogic for the gift of K401-BCCP-H6 plasmids which were used for expressing the kinesin motor proteins for driving the microtubule-based active fluid. We thank Drs. Yuan-Nan Young and David B. Stein for the insightful discussion on interpreting our experimental data. We thank Victoria M. Bicchieri for her assistance on collecting confocal data with the Leica SP5 point scanning confocal microscope in the Life Sciences and Bioengineering Center at Worcester Polytechnic Institute. We thank Ellie Lin from Lin Life Science for her assistance on editing the manuscript to enhance its flow and readability. K.-T.W. acknowledges support from the National Science Foundation (NSF-CBET-2045621). This research was performed with computational resources supported by the Academic & Research Computing Group at Worcester Polytechnic Institute. We acknowledge the Royal Society of Chemistry for adapting the figure from Bate et al. on Soft Matter [36]. We acknowledge the Brandeis Materials Research Science and Engineering Center (NSF-MRSEC-DMR-2011486) for use of the Biological Materials Facility. C.-C.C. acknowledges support for the numerical studies from the Headquarters of University Advancement at the National Cheng Kung University, sponsored by the Ministry of Education, Taiwan.
Author information
Authors and Affiliations
Contributions
Y.-C.C., B.J., C.-C.C., and K.-T.W. performed the research and designed the experiments; Y.-C.C. fabricated the millifluidic devices, prepared the samples, tracked the flows, and reconstructed the microtubule network structures; Y.-C.C., B.J., and C.-C.C. established the continuum simulation platform on modeling active droplet systems; Y.-C.C. and K.-T.W. organized and analyzed the data; Y.-C.C., C.-C.C., and K.-T.W. wrote the manuscript; and K.-T.W. supervised the research. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Supplementary Video 1.
Supplementary Video 2.
Supplementary Video 3.
Supplementary Video 4.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chen, YC., Jolicoeur, B., Chueh, CC. et al. Flow coupling between active and passive fluids across water–oil interfaces. Sci Rep 11, 13965 (2021). https://doi.org/10.1038/s41598-021-93310-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-93310-9
This article is cited by
-
Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity
Nature Communications (2022)
-
A review of tracer testing techniques in porous media specially attributed to the oil and gas industry
Journal of Petroleum Exploration and Production Technology (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.