Source author record

John F. Brady

John F. Brady appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Activity-induced propulsion of a vesicle

Modern biomedical applications such as targeted drug delivery require a delivery system capable of enhanced transport beyond that of passive Brownian diffusion. In this work an osmotic mechanism for the propulsion of a vesicle immersed in a viscous fluid is proposed. By maintaining a steady-state solute gradient inside the vesicle, a seepage flow of the solvent (e.g., water) across the semipermeable membrane is generated which in turn propels the vesicle. We develop a theoretical model for this vesicle-solute system in which the seepage flow is described by a Darcy flow. Using the reciprocal theorem for Stokes flow it is shown that the seepage velocity at the exterior surface of the vesicle generates a thrust force which is balanced by the hydrodynamic drag such that there is no net force on the vesicle. We characterize the motility of the vesicle in relation to the concentration distribution of the solute confined inside the vesicle. Any osmotic solute is able to propel the vesicle so long as a concentration gradient is present. In the present work, we propose active Brownian particles (ABPs) as a solute. To maintain a symmetry-breaking concentration gradient, we consider ABPs with spatially varying swim speed and ABPs with constant properties but under the influence of an orienting field. In particular, it is shown that at high activity the vesicle velocity is $\boldsymbol{U}\sim [K_\perp /(η_e\ell_m) ]\int Π_0^\mathrm{swim} \boldsymbol{n} dΩ$, where $Π_0^\mathrm{swim}$ is the swim pressure just outside the thin accumulation boundary layer on the interior vesicle surface, $\boldsymbol{n}$ is the unit normal vector of the vesicle boundary, $K_\perp$ is the membrane permeability, $η_e$ is the viscosity of the solvent, and $\ell_m$ is the membrane thickness.

preprint2022arXiv

Confined active matter in external fields

We analyze a dilute suspension of active particles confined between walls and subjected to fields that can modulate particle speed as well as orientation. Generally, the particle distribution is different in the bulk compared to near the walls. In the bulk, particles tend to accumulate in the regions of low speed, but in the presence of an orienting field, particles rotate to align with the field and accumulate downstream in the field direction. At the walls, particles tend to accumulate pointing into the walls and thereby exert pressure on walls. But the presence of strong orienting fields can cause the particles to reorient away from the walls, and hence shows a possible mechanism for preventing contamination of surfaces. The pressure at the walls depends on the wall separation and the field strengths. This work demonstrates how multiple fields with different functionalities can be used to control active matter under confinement.

preprint2022arXiv

Forced microrheology of active colloids

Particle-tracking microrheology of dilute active (self-propelled) colloidal suspensions is studied by considering the external force required to maintain the steady motion of an immersed constant-velocity colloidal probe. If the probe speed is zero, the suspension microstructure is isotropic but exhibits a boundary accumulation of active bath particles at contact due to their self-propulsion. As the probe moves through the suspension, the microstructure is distorted from the nonequilibrium isotropic state, which allows us to define a microviscosity for the suspension using the Stokes drag law. For a slow probe, we show that active suspensions exhibit a swim-thinning behavior in which their microviscosity is gradually lowered from that of passive suspensions as the swim speed increases. When the probe speed is fast, the suspension activity is obscured by the rapid advection of the probe and the measured microviscosity is indistinguishable from that of passive suspensions. Generally for finite activity, the suspension exhibits a velocity-thinning behavior -- though with a zero-velocity plateau lower than passive suspensions -- as a function of the probe speed. These behaviors originate from the interplay between the suspension activity and the hard-sphere excluded-volume interaction between the probe and a bath particle.

preprint2022arXiv

Trapped-particle microrheology of active suspensions

In microrheology, the local rheological properties such as viscoelasticity of a complex fluid are inferred from the free or forced motion of embedded colloidal probe particles. Theoretical machinery developed for forced-probe microrheology of colloidal suspensions focused on either constant-force (CF) or constant-velocity (CV) probes while in experiments neither the force nor the kinematics of the probe is fixed. More importantly, the constraint of CF or CV introduces a difficulty in the meaningful quantification of the fluctuations of the probe due to a thermodynamic uncertainty relation. It is known that for a Brownian particle trapped in a harmonic potential well, the product of the standard deviations of the trap force and the particle position is $d k_BT$ in $d$ dimensions with $k_BT$ being the thermal energy. As a result, if the force (position) is not allowed to fluctuate, the position (force) fluctuation becomes infinite. To allow the measurement of fluctuations, in this work we consider a microrheology model in which the embedded probe is dragged along by a moving harmonic potential so that both its position and the trap force are allowed to fluctuate. Starting from the full Smoluchowski equation governing the dynamics of $N$ hard active Brownian particles, we derive a pair Smoluchowski equation describing the dynamics of the probe as it interacts with one bath particle by neglecting hydrodynamic interactions among particles in the dilute limit. From this, we determine the mean and the variance (i.e., fluctuation) of the probe position in terms of the pair probability distribution. We then characterize the behavior of the system in the limits of both weak and strong trap. By taking appropriate limits, we show that our generalized model can be reduced to the well-studied CF or CV microrheology models.

preprint2020arXiv

The "isothermal" compressibility of active matter

We demonstrate that the mechanically-defined "isothermal" compressibility behaves as a thermodynamic-like response function for suspensions of active Brownian particles. The compressibility computed from the active pressure - a combination of the collision and unique swim pressures - is capable of predicting the critical point for motility induced phase separation, as expected from the mechanical stability criterion. We relate this mechanical definition to the static structure factor via an active form of the thermodynamic compressibility equation and find the two to be equivalent, as would be the case for equilibrium systems. This equivalence indicates that compressibility behaves like a thermodynamic response function, even when activity is large. Finally, we discuss the importance of the phase interface when defining an active chemical potential. Previous definitions of the active chemical potential are shown to be accurate above the critical point but breakdown in the coexistence region. Inclusion of the swim pressure in the mechanical compressibility definition suggests that the interface is essential for determining phase behavior.

preprint2019arXiv

Microscopic Origins of the Swim Pressure and the Anomalous Surface Tension of Active Matter

The unique pressure exerted by active particles -- the "swim" pressure -- has proven to be a useful quantity in explaining many of the seemingly confounding behaviors of active particles. However, its use has also resulted in some puzzling findings including an \textit{extremely negative} surface tension between phase separated active particles. Here, we demonstrate that this contradiction stems from the fact that the swim pressure \textit{is not a true pressure}. At a boundary or interface, the reduction in particle swimming generates a net active force density -- an entirely \textit{self-generated body force}. The pressure at the boundary, which was previously identified as the swim pressure, is in fact an elevated (relative to the bulk) value of the \textit{traditional particle pressure} that is generated by this interfacial force density. Recognizing this unique mechanism for stress generation allows us to define a much more physically plausible surface tension. We clarify the utility of the swim pressure as an "equivalent pressure" (analogous to those defined from electrostatic and gravitational body forces) and the conditions in which this concept can be appropriately applied.

preprint2016arXiv

Amorphous and ordered states of concentrated hard spheres under oscillatory shear

Hard sphere colloidal particles are a basic model system for general research into phase behavior, ordering and out-equilibrium glass transitions. Experimentally it has been shown that oscillatory shearing of a monodisperse hard sphere glass, produces two different crystal orientations; a Face Centered Cubic (FCC) crystal with the close packed direction parallel to shear at high strains and an FCC crystal with the close packed direction perpendicular to shear at low strains. Here, using Brownian dynamics simulations of hard sphere particles, we have examined high volume fraction shear-induced crystals under oscillatory shear as well as the same volume fraction glass counterparts. We find that, while the displacements under shear of the glass are isotropic, the sheared FCC crystal structures oriented parallel to shear, are anisotropic due to the cooperative motion of velocity-vorticity layers of particles sliding over each other. These sliding layers generally result in lower stresses and less overall particle displacements. Additionally, from the two crystal types, the perpendicular crystal exhibits less stresses and displacements at smaller strains, however at larger strains, the sliding layers of the parallel crystal are found to be more efficient in minimizing stresses and displacements, while the perpendicular crystal becomes unstable. The findings of this work suggest that the process of shear-induced ordering for a colloidal glass is facilitated by large out of cage displacements, which allow the system to explore the energy landscape and find the minima in energy, stresses and displacements by configuring particles into a crystal oriented parallel to shear.

preprint2015arXiv

Classical Liquids in Fractal Dimension

We introduce fractal liquids by generalizing classical liquids of integer dimensions $d = 1, 2, 3$ to a fractal dimension $d_f$. The particles composing the liquid are fractal objects and their configuration space is also fractal, with the same non-integer dimension. Realizations of our generic model system include microphase separated binary liquids in porous media, and highly branched liquid droplets confined to a fractal polymer backbone in a gel. Here we study the thermodynamics and pair correlations of fractal liquids by computer simulation and semi-analytical statistical mechanics. Our results are based on a model where fractal hard spheres move on a near-critical percolating lattice cluster. The predictions of the fractal Percus-Yevick liquid integral equation compare well with our simulation results.

preprint2015arXiv

Spectral Ewald Acceleration of Stokesian Dynamics for polydisperse suspensions

In this work we develop the Spectral Ewald Accelerated Stokesian Dynamics (SEASD), a novel computational method for dynamic simulations of polydisperse colloidal suspensions with full hydrodynamic interactions. SEASD is based on the framework of Stokesian Dynamics (SD) with extension to compressible solvents, and uses the Spectral Ewald (SE) method [Lindbo & Tornberg, J. Comput. Phys. 229 (2010) 8994] for the wave-space mobility computation. To meet the performance requirement of dynamic simulations, we use Graphic Processing Units (GPU) to evaluate the suspension mobility, and achieve an order of magnitude speedup compared to a CPU implementation. For further speedup, we develop a novel far-field block-diagonal preconditioner to reduce the far-field evaluations in the iterative solver, and SEASD-nf, a polydisperse extension of the mean-field Brownian approximation of Banchio & Brady [J. Chem. Phys. 118 (2003) 10323]. We extensively discuss implementation and parameter selection strategies in SEASD, and demonstrate the spectral accuracy in the mobility evaluation and the overall $\mathcal{O}(N\log N)$ computation scaling. We present three computational examples to further validate SEASD and SEASD-nf in monodisperse and bidisperse suspensions: the short-time transport properties, the equilibrium osmotic pressure and viscoelastic moduli, and the steady shear Brownian rheology. Our validation results show that the agreement between SEASD and SEASD-nf is satisfactory over a wide range of parameters, and also provide significant insight into the dynamics of polydisperse colloidal suspensions.

preprint2014arXiv

Short-time diffusion in concentrated bidisperse hard-sphere suspensions

Diffusion in bidisperse Brownian hard-sphere suspensions is studied by Stokesian Dynamics (SD) computer simulations and a semi-analytical theoretical scheme for colloidal short-time dynamics, based on Beenakker and Mazur's method [Physica 120A, 388 (1983) & 126A, 349 (1984)]. Two species of hard spheres are suspended in an overdamped viscous solvent that mediates the salient hydrodynamic interactions among all particles. In a comprehensive parameter scan that covers various packing fractions and suspension compositions, we employ numerically accurate SD simulations to compute the initial diffusive relaxation of density modulations at the Brownian time scale, quantified by the partial hydrodynamic functions. A revised version of Beenakker and Mazur's $δγ$-scheme for monodisperse suspensions is found to exhibit surprisingly good accuracy, when simple rescaling laws are invoked in its application to mixtures. The so-modified $δγ$ scheme predicts hydrodynamic functions in very good agreement with our SD simulation results, for all densities from the very dilute limit up to packing fractions as high as $40\%$.

preprint2014arXiv

Short-time transport properties of bidisperse suspensions and porous media: a Stokesian Dynamics study

We present a comprehensive computational study of the short-time transport properties of bidisperse neutral colloidal suspensions and the corresponding porous media. Our study covers bidisperse particle size ratios up to $4$, and total volume fractions up to and beyond the monodisperse hard-sphere close packing limit. The many-body hydrodynamic interactions are computed using conventional Stokesian Dynamics (SD) via a Monte-Carlo approach. We address suspension properties including the short-time translational and rotational self-diffusivities, the instantaneous sedimentation velocity, the wavenumber-dependent partial hydrodynamic functions, and the high-frequency shear and bulk viscosities; and porous media properties including the permeability and the translational and rotational hindered diffusivities. We carefully compare the SD computations with existing theoretical and numerical results. For suspensions, we also explore the range of validity of various approximation schemes, notably the Pairwise Additive (PA) approximations with the Percus-Yevick structural input. We critically assess the strengths and weaknesses of the SD algorithm for various transport properties. For very dense systems, we discuss in detail the interplay between the hydrodynamic interactions and the structures due to the presence of a second species of a different size.

preprint2014arXiv

Towards a 'Thermodynamics' of Active Matter

Self-propulsion allows living systems to display unusual collective behavior. Unlike passive systems in thermal equilibrium, active matter systems are not constrained by conventional thermodynamic laws. A question arises however as to what extent, if any, can concepts from classical thermodynamics be applied to nonequilibrium systems like active matter. Here we use the new swim pressure perspective to develop a simple theory for predicting phase separation in active matter. Using purely mechanical arguments we generate a phase diagram with a spinodal and critical point, and define a nonequilibrium chemical potential to interpret the "binodal." We provide a generalization of thermodynamic concepts like the free energy and temperature for nonequilibrium active systems. Our theory agrees with existing simulation data both qualitatively and quantitatively and may provide a framework for understanding and predicting the behavior of nonequilibrium active systems.