ORIGINAL RESEARCH article
Porous Media Microstructure Determines the Diffusion of Active Matter: Experiments and Simulations
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, CA, United States
Active swimmers are known to accumulate along external boundaries owing to their persistent self-motion, resulting in a significant reduction in their effective mobility through heterogeneous and tortuous materials. The dynamic interplay between the slowdown experienced by the active constituents near boundaries and their long-time diffusivity is critical for understanding and predicting active transport in porous media. In this work, we study the impact of boundary layer accumulation on the effective diffusivity of active matter by analyzing the motion of active Brownian particles in an array of fixed obstacles. We combine Janus particle experiments, Brownian dynamics simulations, and a theoretical analysis based on the Smoluchowski equation. We find that the shape, curvature, and microstructure of the obstacles play a critical role in governing the effective diffusivity of active particles. Indeed, even at dilute packing fractions of obstacles, ϕ = 12%, we observed a 25% reduction in the effective diffusivity of active particles, which is much larger than the hindrance experienced by passive Brownian particles. Our combined experimental and computational results demonstrate a strong coupling between the active force and the porous media microstructure. This work provides a framework to predict and control the transport of active matter in heterogeneous materials.
The transport of living, colloidal-sized species through crowded environments plays a crucial function in many natural and synthetic processes. For example, the transport of bacteria through soil plays a beneficial role in bioremediation [1, 2], and novel drug delivery mechanisms seek to utilize the proliferation of S. typhimurium to access tumor tissues that have been conventionally out of reach [3, 4]. In contrast, the transport of pathogens into wounds sites and mucosa can lead to life-threatening infections without proper treatment [5–7]. The effective transport properties of bacteria in crowded environments depend on the interplay between the swimmer motility and the boundaries that make up the porous material . Many forms of microscopic life enhance their transport via directed self-propulsion, including E. coli bacteria, spermatozoa cells, and C. reinhardtii algae [9–11]. Recent advancements in micro/nanoscale synthesis have also led to the creation of artificial swimmers that are excellent tools in the study of autonomous self-propulsion [12–19]. Understanding the motion of these living and synthetic “active matter” constituents embedded within heterogeneous materials is a challenging problem because of the complex interactions between the swimmer and the material boundaries. For porous materials composed of polymer networks, particle transport may be affected by steric hindrance, nonspecific interactions (hydrophobicity, electrostatics), and specific interactions (ligand-receptor binding) [20, 21]. This behavior is not unique to polymer networks; introducing even a dilute concentration of immobile obstacles with purely excluded-volume interactions provides a substantial slowdown to diffusive flux [22–24].
In addition to the interactions experienced by non-motile Brownian particles, active particles accumulate at physical boundaries due to their persistent self-motion, characterized by a boundary layer near the surface. This accumulation occurs even in the absence of attractive interactions; the active particles propel freely until hitting a surface and continue to propel themselves toward the surface until they reorient and escape into the bulk fluid. This behavior has been observed experimentally and in simulations of rods and spheres, both with and without hydrodynamic interactions [25–40]. A mechanistic understanding of how the local accumulation near physical boundaries affects the macroscopic mobility of active matter through heterogeneous and tortuous materials is lacking. In this work, we study the impact of boundary accumulation on the effective diffusivity of active matter by analyzing the motion of active Brownian particles (ABPs) in a system of rigid 2D obstacles. The presence of boundaries in active systems reduces the effective long-time self diffusivity by an amount that depends on the average swimming speed (U0) and the average reorientation time (τR). This is in sharp contrast to passive Brownian particles, which do not accumulate along boundaries, and whose effective diffusivity depends primarily on the packing fraction of the obstacles [23, 41].
Many studies on active matter transport focus on the effect of alignment along surfaces due to steric or hydrodynamic torques aligning the swimmer parallel or perpendicular to boundaries [27, 36, 38, 40, 42–46]. However, a connection between transport and surface accumulation of active matter [28, 32–34, 37, 47, 48] without any imposed torques has not been fully explored. For active systems, the precise shape and curvature of the boundary can have a strong effect on motility induced accumulation [33, 34, 49]. This increased accumulation corresponds to more time spent “trapped” in the boundary layer, which inhibits the transport of active matter in tight pores. Therefore, we hypothesize that active swimmers experience a strong reduction in the relative diffusivity in porous media due to the synergistic effects of active boundary accumulation and boundary shape.
While many existing theories predict the diffusive transport of passive Brownian particles through porous media [20, 22, 50–54], the unique accumulation of active matter along boundaries—especially at regions of large curvature—leads to unexpected diffusive slowdowns that are not captured in traditional theories. In this work, we combine Janus particle experiments, Brownian dynamics simulations, and theory to demonstrate that the transport of active matter in heterogeneous materials is a strong function of the obstacle shape, curvature, and microstructure. Experimentally, we rely upon optical tracking of active particle trajectories, which has been a powerful tool to study both living and artificial swimmers in porous environments [49, 55–62], and allows for direct comparisons to particle based simulations [63–71]. In addition to advancing our basic understanding of active matter transport, our work provides a mechanism to control the transport of active matter in heterogeneous materials.
To obtain a mechanistic understanding of boundary layer accumulation and slowdown of active matter in heterogeneous materials, we combined Janus particle experiments, Brownian dynamics (BD) simulations, and analytical theory. In our experiments, we immobilized 6.5 μm lipid bilayer-coated silica particles in a random distribution at the bottom of an imaging chamber at ϕ ≈ 12% area fraction. We added a dilute concentration of 4 μm silica Janus particles, coated on one side with a thin layer of platinum and the other side with a lipid bilayer containing fluorescently-labeled lipids (see Materials and Methods). The silica beads sedimented to the bottom of the imaging chamber, so our experiments are conducted in 2D. The lipid bilayers on the obstacles and Janus half-coating contain different fluorescent dyes, which enabled us to track both types of particles simultaneously in different fluorescence channels. Upon adding 2% hydrogen peroxide in Milli-Q water, the Janus particles self-propelled [72–74] with speed U0 = 0.84 ± 0.01 μm/s and reorientation time τR = 14 ± 2 s. The self-propulsive speed and the reorientation time were determined via the mean instantaneous velocity and a fit to the known mean squared displacement (see in Materials and Methods). We conducted time lapse imaging and tracked the positions of the obstacles and the Janus particles using a tracking algorithm . In Figure 1, we show the motion of a single Janus particle moving through a random array of obstacles, punctuated by an obstructed motion of over 1 min in a local grouping of obstacles creating a concave boundary. Eventually, the Janus particle reoriented, propelled away from the concave boundary, and resumed an active random walk (Figure 1H).
FIGURE 1. Active particles experience a significant slowdown near obstacle boundaries and a reduction in effective diffusivity in a random array of fixed obstacles. (A–C): Schematic demonstrating the obstructed motion of an active particle near a grouping of three obstacles. (D–G): Experimental images of a 4 μm Janus particle (green semicircle) moving through 6.5 μm obstacle particles (magenta). (H) Displacement of the active Janus particle (green semicircle) tracked for 100 s. (I) Brownian Dynamics (BD) simulation snapshot of active Brownian particles moving through a recreated copy of the same porous media as in the experiments. More individual particle trajectories are available in the Supplementary Information (Supplementary Figure S1).
To corroborate our experiments, we developed BD simulations in which the motion of ABPs are evolved following the overdamped Langevin equation (see Materials and Methods). We compare the experimental mean squared displacement (MSD) with simulated active Brownian particle MSD to determine if this simple model quantitatively captures the transport behavior observed in the Janus particle experiments in Figure 2. Comparing the two MSDs also allowed us to determine that the entropic effect of immobile obstacles was the cause of the diffusivity reduction, and not some unaccounted for mechanism (e.g. hydrodynamic forces or interparticle attractions). To simulate a dilute system, the active Brownian particles interact with obstacle particles via a purely repulsive potential, but do not interact with each other (“ideal gas” particles). We chose the active particle swimming speed and reorientation time to match our Janus particles, and we placed obstacles of the same size in the same positions as the experiments. By using the experimental obstacle particle positions as inputs into our simulations, we recapitulated our precise experimental system in the simulations (Figure 1I). Consistent with our experimental observations, we also observed a similar accumulation of particles in local groupings of obstacles that form a concave boundary.
FIGURE 2. The effective diffusivity of active particles decreased by 25% in a random array of fixed obstacles at ϕ ≈12% area fraction. The mean squared displacement (MSD) of 4 μm active Janus colloids in experiments (black symbols) agrees with that in our BD simulation of the equivalent system using the active Brownian particle (ABP) model (cyan symbols), for both the control case without obstacles (circles) and with a random distribution of 6.5 μm obstacles (cross marks). Labelled diffusivities correspond to the fit to the slope at times past 50 s in the experiments. Standard error of the mean for the MSDs is available in the Supplementary Figure S2. Additional MSD simulation data for different ABP activities is available in Supplementary Figure S3.
To quantify the effect of active particle accumulation and slowdown near boundaries, we computed the mean squared displacement (MSD) of the active particles in our experiments and BD simulations, MSD(t) = ⟨|r(t) − r (0)|2⟩, where r(t) is the position of the active particle at time t. We obtained the slope of the MSD at large times to find the long-time self diffusivity of the active particles in the experiments and simulations,
In our experiments, we noticed that several 6.5 μm obstacle particles formed local groupings with narrow constrictions (Figure 3A), even at semi-dilute packing fractions (ϕ ≈ 12%). Indeed, we quantified the crowding by finding the number of other obstacle neighbors located within a surface-to-surface distance of 4 μm (Figure 3B). We observed that the Janus particles spent significantly more time in the concave region of these emergent shapes as opposed to the convex side (Figure 3C). We therefore hypothesized that the presence of these structures with curved geometries plays an important role in reducing the effective transport properties of active particles beyond the 1 − ϕ correction observed in passive particles. Furthermore, we anticipated that the role of obstacle geometry on the effective diffusivity is much larger for active particles undergoing persistent self-propulsion compared to passive particles translating due to thermal Brownian motion.
FIGURE 3. Active Janus particles spent a significant amount of time near obstacle clusters that formed curved boundaries with narrow constriction sites. (A) Image of an active Janus particle interacting with a grouping of obstacles that forms a curved boundary with concavity. (B) Histogram quantifying obstacle particle grouping. Number of impassible neighbors is found by counting all the obstacle neighbors within distance 4 μm from a reference obstacle’s surface. Groups of obstacles within this distance are impassible by a Janus swimmer with diameter 4 μm. Error bars are the standard deviation from six independent experiments. (C) Schematic demonstrating the concave shape formed by the packing of obstacles.
To test our hypothesis and to develop a micromechanical understanding of the role of obstacle geometry on active particle diffusion, we analyzed the distribution of active particles near obstacle boundaries, P (x, y, θ, t), which satisfies the Smoluchowski equation
where U0 is the self-propulsive speed of the active particles, q = [ cos(θ), sin(θ)] is the unit orientation vector indicating the direction of self-propulsion, and DT and DR = 1/τR are the translational and rotational diffusivity, respectively. Eq. 1 is subject to the no-flux boundary condition along the obstacle surface,
We obtained steady-state density and polarization fields by taking orientational integrals over the full probability distribution,
We numerically solved the full Smoluchowski equation for a point-sized active particle around fixed obstacles with different shapes. For a circular obstacle, we observed a small accumulation of active particles near the surface, as shown in Figure 4A. In contrast, for a curved obstacle, we observed a significant increase in the number density n (x, y) and polar order vector m (x, y) near the obstacle surface, especially at regions of large concavity (Figures 4D–F). Our results in Figure 4 are presented for a mild activity of U0τR/δ = 1, where
FIGURE 4. Active constituents accumulate along boundaries due to their persistent self-motion, with a significant increase in density and polar order along curved boundaries with large concavity. (A–C): Density and polarization fields around a circular obstacle (radius over run length
Motivated by our micromechanical understanding of active particles near curved obstacles (Figure 4), we hypothesized that the effective diffusivity of active particles in an array of obstacles should depend on the specific shape of the obstacles. To these ends, we conducted BD simulations of active Brownian particles moving through a square lattice of obstacles, carefully varying the curvature while preserving the packing fraction to keep the excluded volume constant within a unit cell. As a model obstacle shape with smoothly-varying curvature, we utilized the “Cassini Oval” (Figures 5B–D), which is described by the equation
where a and b are two shape parameters (a < b), κ* is the maximum curvature in the shape, A is the shape area, and E(x) is the complete elliptic integral of the second kind. The Cassini Oval is a modification of the traditional ellipse with the product of the distance to two foci (located at x = ±a) kept constant at b2. The shape extends laterally and shrinks vertically as it is deformed at constant area, which would generate anisotropies and slowdowns in the effective diffusivity for even passive Brownian particles. Since we wish to isolate the effects of curvature, and not the artifacts from lateral elongation of the shape, we performed BD simulations on both passive and active Brownian particles to quantify the effects of curvature and shape elongation. Passive Brownian particles with purely excluded-volume interactions do not accumulate along boundaries, so any change in their diffusivity is due to shape elongation within the unit cell. We have conducted BD simulations at different activities and obstacle packing fractions, and found that the effective diffusivity is well-approximated by the expression D0 (1 − ϕ) for dilute packing fractions, where D0 is the bulk diffusivity in 2D in the absence of any obstacles. This is a proficient analytical expression for all activities at dilute obstacle densities (See Supplementary Figures S4, S5). However, at larger packing fractions of obstacles greater than ϕ ≈ 5%, we observe deviations in this expression as a function of varying activity parameters. For example, as (U0τR)κ* > 1, the scaled diffusivity decreases compared to the passive case, due to the reduced swim diffusivity in the boundary layer.
FIGURE 5. Curved obstacles with a large concavity generate a significant slowdown in diffusive transport of active particles compared to purely convex shapes at the same packing fraction. (A) Scaled diffusivity deviation as a function of shape concavity for the Cassini Oval shape described by Eq. 3, where bκ* is the non-dimensional curvature at the top of the inclusion. The diffusivity along the direction facing the concavity, Dy, decreases by 24% for the active particles for highly concave obstacles compared to 12.5% for the passive Brownian particles, indicating that the coupling between the active force and curvature plays a key role. Total packing fraction is fixed at ϕ ≈5.5%, and error bars are the standard deviation of three independent trials. (B–D): Schematics of area preserving Cassini Ovals for three different parameter sets. The two foci are separated a distance of 2a. Beyond a critical distance,
To isolate the diffusivity reduction due to obstacle shape, we computed a scaled diffusivity deviation given by
where DE is the effective diffusivity measured from the MSD,
As shown in our data for Dy in Figure 5A, the active-curvature coupling contribution to the effective diffusivity is equally as large as the diffusivity reduction due to obstacle shape elongation. Due to their persistent self-propulsion, the active particles experience a large accumulation of density and polar order fields near boundaries with large concavity, consistent with our Smoluchowski analysis in Figure 4. The magnitude of this boundary layer accumulation is a strong function of curvature and activity, and we have observed that the effect becomes more important at high activity and semi-dilute obstacle packing (see Supplementary Figures S4, S5). Our results validate our hypothesis that the shape and curvature of the obstacles play a critical role in governing the effective mobility of active particles in porous materials. The physical mechanism behind this phenomenon is that active matter accumulates along boundaries, especially along curved and concave surfaces, where active particles are trapped. Therefore, the specific shape and arrangement of obstacles within a porous material modulate the effective diffusivity of the active particles in a manner that is more significant compared to passive Brownian particles.
In this work, we discovered that the obstacle packing fraction alone is insufficient to provide an accurate prediction of effective active particle diffusivity. The specific shape and distribution of physical obstacles plays a critical role in determining active transport. The microscopic details of the external boundary strongly influences the macroscopic observables like the long-time self diffusivity. Both in our experiments and simulations, a random packing of obstacles led to concave structures with narrow constriction sites that gave rise to a significant accumulation of active particles. We showed that the local slowdown of active particles within the boundary layers has a direct effect on their overall mobility across the porous material. We focused on obstacle curvature and microstructure as metrics for predicting the effective diffusivity, which is complementary to other porosity metrics, like chord-length distributions  or tortuosity measurements .
Our scientific basis for focusing on the obstacle shape and curvature was inspired by the strong connection between active matter accumulation and its force generation along boundaries [32, 74, 78, 79]. For example, Burkholder and Brady derived a macrotransport based model to connect fluctuations in surface accumulation and the enhancement in diffusivity of passive spherical tracers . Furthermore, the surface accumulation is highly dependent on surface curvature [28, 31, 33, 34]; the accumulation on certain parts of an asymmetric shape can lead to a pressure imbalance and net translation of anisotropic colloidal tracers [11, 81]. These studies showed that active matter can impart forces on its environment; however, in our work, we focused on how the micromechanical details of the environment, like curved boundaries, can alter the dynamical properties of the active particles. We validated that a strong coupling between surface curvature and active matter accumulation decreases the diffusivity of active particles by a much larger relative degree than the slowdown expected for equivalent passive Brownian particles.
The surface accumulation of active particles around a single obstacle with small curvature ((U0τR), δ ≪ 1/κ) is shown by Yan and Brady :
where n∞ is the bulk concentration of active particles, κ is the curvature in units of inverse obstacle length, ℓ = U0τR is the run length, and
We focused on the effects of immobile obstacles with purely excluded volume interactions. However, soft porous materials present a rich opportunity for future study. Boundary fluctuations of soft surfaces are important for particle transport in mucus, hydrogels, and other polymeric networks [20, 84, 85]. Active particles have been shown to induce large deformations on soft membranes [86, 87], changing the curvature and transport drastically. Hydrogel networks and sediments in 3D provide an additional degree of freedom for a swimmer to avoid obstacles. In 2D, close packing of disks precludes transport; however, a close packing of spheres in 3D allows for bicontinous percolation, resulting in a reduced but nonzero diffusivity.
Our work opens up opportunities for future experimental work to control active matter diffusion via the design of obstacle shape and arrangement. Convex, nonspherical inclusions can be used to control transport and create anisotropic spreading of bacteria along a predefined axis. Novel sorting mechanisms have already been developed using asymmetric blockers , and the method could be extended to sort mixtures of swimmers with different types of motility . Precise consideration of the activity-curvature coupling on transport serves as a promising route to increase the efficacy of these sorting methods, and may lead to more accurate predictions of bacterial transport coefficients.
Materials and Methods
Lipid-coated silica beads were created by coating silica micro-beads with a supported lipid bilayer (SLB). Small unilamellar vesicles (SUVs) were prepared by rehydrating a lipid sheet composed of a mixture of phospholipids with pure deionized water to a concentration of 0.2 mg/ml. For the Janus particles, we used 1,2-dioleoyl-sn-glycero-3-phos-phocholine (DOPC) with 5% of 1,2-dioleoyl-sn-glycero-3-phospho-l-serine (DOPS) and 0.3% of Atto 488-1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE-Atto 488) fluorescent dye. For the obstacle particles, we used DOPC with 5% of 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP) and 0.3% of Atto 647-DOPE (DOPE-Atto 647). After rehydrating the lipids for 30 min, the solution was vigorously vortexed, sonicated at low-power (20% power) using a tip sonicator (Branson SFX250 Sonifier). The resulting SUV solution was buffered with a MOPS buffer (50 mM MOPS, 100 mM sodium chloride, pH 7.5). DOPC (catalog number: 850375), DOPS (catalog number: 840035), and DOTAP (catalog number: 890890) were purchased from Avanti polar lipids. DOPE-Atto 488 and DOPE-Atto 647 were purchased from ATTO-TEC GmbH. Silica microspheres (4.0μm; catalog code: SS05002 and 6.5μm; catalog code: SS06N) were purchased from Bangs Laboratories.
Silica microspheres with diameters 4 and 6.5 μm were cleaned using a 3:2 mixture of sulfuric acid:hydrogen peroxide (Piranha) for 30 min in a bath sonicator, and were spun down at 1000g and washed 3 times before being resuspended in pure water. We fabricated Janus particles from the cleaned 4 μm particles by depositing a monolayer on a glass slide, and coating half of the particle surface with a 2 nm-thick layer of chromium and 8 nm-thick layer of platinum using an E-beam evaporator at a deposition rate of 0.1
For the 4 μm Janus particles, the SLBs coated only half of the particle surface exposed to clean silica. The side with the platinum did not get coated with an SLB. When these SLB half-coated Janus particles were deposited in a 2% solution of hydrogen peroxide, the particles self-propelled pointing away from their platinum half-coating via self-diffusiophoresis. Since silica is more dense than water, the Janus particles moved in 2D along the bottom of the imaging chamber. Within the time frame of our experimental measurements, we did not observe any significant degradation of the SLB from hydrogen peroxide.
For the 6.5 μm obstacle particles, we obtained a uniform SLB across the entire surface of the silica bead. We added a positively-charged DOTAP lipid to the SLB to facilitate a strong electrostatic attraction between the obstacle particles and the borosilicate coverslip substrate. We found that most obstacle particles remained fixed along the bottom substrate with no observable Brownian motion. We did not observe any adhesion or fusion of the SLBs between the obstacles and the Janus particles.
The SLB-coated obstacle particles were added into the imaging chamber at a desired density, followed by the SLB half-coated Janus particles. We added 2% hydrogen peroxide into the chamber and conducted time lapse imaging. All imaging was carried out on an inverted Nikon Ti2-Eclipse microscope (Nikon Instruments) using an oil-immersion objective (Apo ×60, numerical aperture (NA) 1.4, oil). Lumencor SpectraX Multi-Line LED Light Source was used for excitation (Lumencor, Inc.). Fluorescent light was spectrally filtered with emission filters (515/30 and 680/42, Semrock, IDEX Health and Science) and imaged on a Photometrics Prime 95 CMOS Camera (Teledyne Photometrics). Experimental results presented in this work are an average over six independent replicates of systems with obstacles and three independent replicates of systems without obstacles as a control.
Particle Tracking in Experiments
To determine the effect of the porous media on transport, we measured the Janus particle trajectories with and without the presence of the fixed obstacles. We used a modified MATLAB script based on IDL code by Crocker and Grier  to track the individual Janus particles by identifying each particle center and tracking its trajectory over time using an image stack with one frame taken every second. We removed any macroscopic drifts by enforcing that the mean displacement over all particles was zero at any time. We filtered out any Janus particles that were immobile due to defects of the particle (defined as moving less than 30 μm over 100 s or if it moved in only one direction via macroscopic drifts). In all experimental results, we tracked the particles for times t > 7τR. We obtained the Janus particle mean swim speed, U0, and reorientation time, τR, using the control experiments in the absence of obstacles. We obtained the mean swim speed by averaging all tracked particles’ velocity over time,
The mean swim speed was determined to be U0 = 0.84 ± 0.01 μm/s with the reported error as the standard error of the mean. We obtained the reorientation time by measuring the bulk diffusivity for the control experiment in the absence of obstacles,
where DT ≈ 0.1 μm2/s is the thermal Brownian diffusivity for 4 μm diameter spheres using the Stokes-Einstein-Sutherland relation [89–91]. We note that the thermal Brownian diffusivity is negligible compared to the self-propulsive component. Using the experimental measurement of the bulk diffusivity, D0 = 5.1 μm2/s, we obtained the reorientation time and standard deviation τR = 14 ± 2 s using Eq. 7. As a separate measurement, we computed τR using the Janus particle orientation autocorrelation in 2D,
We obtained the Janus particle orientations directly from particle tracking, and we obtained a reorientation time τR ≈ 10 s using Eq. 8. This measurement is similar to the value we obtained using Eq. 7, especially considering the difficulty in finding the orientation using the velocity vector. We concluded that the particles are behaving as active Brownian particles in 2D. In principle, the Janus particles are located along a 2D plane but can reorient in 3D. However, the platinum coating makes the catalytic half-surface more heavy and causes the Janus particles to tilt downwards, making the particle move effectively in 2D. We note that the direction of self-propulsion points away from the platinum half-surface, so the configuration of the platinum half-surface pointing vertically down appears to be an unstable state.
Brownian Dynamics Simulations
The parameters (U0, τR) obtained from the control experiment are then used as inputs into our Brownian dynamics simulations in Figure 2. In Figure 5, the parameters were chosen as described in the caption. We implemented our simulations using HOOMD-blue, a molecular dynamics (MD) simulation package in Python . We focused on the dilute limit of a single active Brownian particle (ABP) in 2D interacting with fixed hard-sphere obstacles. Hydrodynamic interactions are ignored in these simulations. The ABP model describes a swimmer with constant propulsion force but white noise torques [11, 93–96]. Hard-sphere like interactions between the obstacles and the ABPs were implemented using the Weeks-Chandler-Andersen (WCA)  potential (Eq. 10). For the nonspherical Cassini Oval, the structure was formed using overlapping rigid surface particles offset so that the surfaces of the particles formed the boundaries of the Cassini Oval.
The ABPs were initialized and integrated according to the overdamped Langevin equations of motion:
where Fwca is the force on the particle from all potentials and constraints, ζ is the drag coefficient, q = [ cos θ, sin θ] is the unit vector describing particle orientation in 2D, DR = 1/τR is the rotational diffusion coefficient, and (ηi, ξi) are random variables obeying the zero mean and variance consistent with the fluctuation-dissipation theorem. We used a timestep of Δt = 0.001 s, and set the thermal diffusivity to match our experiments at DT = 0.1 μm2/s. The drag coefficient and the energy scale of the potential were chosen such that force induced velocity at contact is
The particle diameters σABP and σobs were preset to 4 and 6.5 μm for Figure 2, and for Figure 5 they were chosen to be small compared to the radius of the disk shown in Figure 5B, σ/R = 1/20. We used periodic boundary conditions to simulate a continuous domain. Simulations are visualized using OVITO .
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
KM and ST conceived of the study; KM and ST designed research; KM performed finite element calculations; KM and YX performed simulations and analyzed data; ST performed experiments; ST supervised the study; and KM, YX, and ST wrote the paper.
This material is based upon work supported by the Air Force Office of Scientific Research under award number FA9550-21-1-0287. KM is supported by the National Science Foundation Graduate Research Fellowship under Grant No. 1650114. Use was made of computational facilities purchased with funds from the National Science Foundation (OAC-1925717) and administered by the Center for Scientific Computing (CSC). The CSC is supported by the California NanoSystems Institute and the Materials Research Science and Engineering Center (MRSEC; NSF DMR 1720256) at UC Santa Barbara.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The authors would like to thank Parth Shah for help in Janus particle synthesis, Dr. Joseph Barakat for a careful reading of this manuscript, and Prof. Ahmad K. Omar for enlightening conversations and insights that helped to conceive of this project.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2022.869175/full#supplementary-material
2. Adadevoh JST, Triolo S, Ramsburg CA, Ford RM. Chemotaxis Increases the Residence Time of Bacteria in Granular Media Containing Distributed Contaminant Sources. Environ Sci Technol (2016) 50:181–7. doi:10.1021/acs.est.5b03956
7. Kao CY, Lin WH, Tseng CC, Wu AB, Wang MC, Wu JJ. The Complex Interplay Among Bacterial Motility and Virulence Factors in Different Escherichia coli Infections. Eur J Clin Microbiol Infect Dis (2014) 33:2157–62. doi:10.1007/s10096-014-2171-2
9. Chaban B, Hughes HV, Beeby M. The Flagellum in Bacterial Pathogens: For Motility and a Whole Lot More. Semin Cel Developmental BiologyBiomineralisation Motorisation Pathog (2015) 46:91–103. doi:10.1016/j.semcdb.2015.10.032
13. Howse JR, Jones RAL, Ryan AJ, Gough T, Vafabakhsh R, Golestanian R. Self-Motile Colloidal Particles: From Directed Propulsion to Random Walk. Phys Rev Lett (2007) 99:048102. doi:10.1103/PhysRevLett.99.048102
55. Perez LJ, Bhattacharjee T, Datta SS, Parashar R, Sund NL. Impact of Confined Geometries on Hopping and Trapping of Motile Bacteria in Porous media. Phys Rev E (2021) 103:012611. doi:10.1103/PhysRevE.103.012611
58. Truong VK, Mainwaring DE, Murugaraj P, Nguyen DHK, Ivanova EP. Impact of Confining 3-D Polymer Networks on Dynamics of Bacterial Ingress and Self-Organisation. J Mater Chem B (2015) 3:8704–10. doi:10.1039/c5tb01880c
64. Kurzthaler C, Mandal S, Bhattacharjee T, Löwen H, Datta SS, Stone HA. A Geometric Criterion for the Optimal Spreading of Active Polymers in Porous media. Nat Commun (2021) 12:7088. doi:10.1038/s41467-021-26942-0
71. Yuan C, Chen A, Zhang B, Zhao N. Activity–crowding Coupling Effect on the Diffusion Dynamics of a Self-Propelled Particle in Polymer Solutions. Phys Chem Chem Phys (2019) 21:24112–25. doi:10.1039/C9CP04498A
73. Howse JR, Jones RAL, Ryan AJ, Gough T, Vafabakhsh R, Golestanian R. Self-motile Colloidal Particles: From Directed Propulsion to Random Walk. Phys Rev Lett (2007) 99:048102. doi:10.1103/PhysRevLett.99.048102
83. Thornton KL, Butler JK, Davis SJ, Baxter BK, Wilson LG. Haloarchaea Swim Slowly for Optimal Chemotactic Efficiency in Low Nutrient Environments. Nat Commun (2020) 11:4453. doi:10.1038/s41467-020-18253-7
84. Bodrenko IV, Salis S, Acosta-Gutierrez S, Ceccarelli M. Diffusion of Large Particles through Small Pores: From Entropic to Enthalpic Transport. J Chem Phys (2019) 150:211102. doi:10.1063/1.5098868
87. Vutukuri HR, Hoore M, Abaurrea-Velasco C, van Buren L, Dutto A, Auth T, et al. Active Particles Induce Large Shape Deformations in Giant Lipid Vesicles. Nature (2020) 586:52–6. doi:10.1038/s41586-020-2730-x
92. Anderson JA, Glaser J, Glotzer SC. HOOMD-Blue: A Python Package for High-Performance Molecular Dynamics and Hard Particle Monte Carlo Simulations. Comput Mater Sci (2020) 173:109363. doi:10.1016/j.commatsci.2019.109363
Keywords: active matter, colloidal transport, porous media, brownian motion, active suspensions, janus particles
Citation: Modica KJ, Xi Y and Takatori SC (2022) Porous Media Microstructure Determines the Diffusion of Active Matter: Experiments and Simulations. Front. Phys. 10:869175. doi: 10.3389/fphy.2022.869175
Received: 04 February 2022; Accepted: 04 March 2022;
Published: 08 April 2022.
Edited by:Sujit Datta, Princeton University, United States
Reviewed by:Tapomoy Bhattacharjee, National Centre for Biological Sciences, India
Christina Kurzthaler, Princeton University, United States
Copyright © 2022 Modica, Xi and Takatori. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sho C. Takatori, email@example.com