Trapping, gliding, vaulting: transport of semiflexible polymers in periodic post arrays

Transport and clogging dynamics of flexible rods in pore constrictions

The transport and clogging behavior of flexible particles in confined flows is a complex interplay between elastic and hydrodynamic forces and wall interactions. While the motion of non-spherical particles in unbounded flows is well understood, their behavior in confined spaces remains less explored. This study introduces a coupled computational fluid dynamics-discrete element method (CFD-DEM) approach to investigate the transport and clogging dynamics of flexible rod-shaped particles in confined pore constrictions. The spatio-temporal analysis reveals the influence of the rod's initial conditions and flexibility on its transport dynamics through a pore constriction. The simulation results demonstrate an increase in the lateral drift of the rod upon exiting the pore that can be scaled with channel height confinement. The clogging dynamics are explored based on hydrodynamic and mechanical forces, unveiling conditions for mechanical clogging through sieving. The developed method allows for the deconvolution of the forces that contribute to particle trajectories in confined flow, which is highly relevant in particle separation processes, fibrous-shaped virus filtration, biological flows, and related applications. The method is embedded into the open-source CFDEM framework, facilitating future extensions to explore multiple particle dynamics, intermolecular forces, external influences, and complex geometries.

Non-equilibrium structural and dynamic behaviors of active polymers in complex and crowded environments

Active matter systems, which convert internal chemical energy or energy from the environment into directed motion, are ubiquitous in nature and exhibit a range of emerging non-equilibrium behaviors. However, most of the current works on active matter have been devoted to particles, and the study of active polymers has only recently come into the spotlight due to their prevalence within living organisms. The intricate interplay between activity and conformational degrees of freedom gives rise to novel structural and dynamical behaviors of active polymers. Research in active polymers remarkably broadens diverse concepts of polymer physics, such as molecular architecture, dynamics, scaling and so on, which is of significant importance for the development of new polymer materials with unique performance. Furthermore, active polymers are often found in strongly interacting and crowded systems and in complex environments, so that the understanding of this behavior is essential for future developments of novel polymer-based biomaterials. This review thereby focuses on the study of active polymers in complex and crowded environments, and aims to provide insights into the fundamental physics underlying the adaptive and collective behaviors far from equilibrium, as well as the open challenges that the field is currently facing.

Cross-streamline migration and near-wall depletion of elastic fibers in micro-channel flows

The complex dynamics of elastic fibers in viscous fluids are central to many biological and industrial systems. Fluid-structure interactions underlying these dynamics govern the shape and transport of flexible fibers, and understanding these interactions can help tune flow properties in applications such as microfluidic separation, printing and clogging. In this work, we use slender-body theory to study micromechanical dynamics that arise from the coupling between the elastic backbone of a fiber and the local straining flow that contributes to filament flipping and cross-streamline migration. The resulting transverse drift is unbiased in either direction in simple shear flow. However, a non-uniform shear rate results in bias towards regions of high shear, which we connect to the shape transitions during flips. We discover a depletion layer that forms near the boundaries of pressure-driven channel flow due to the competition between such a cross-streamline drift and steric exclusion from the walls. Finally, we develop scaling laws for the curvature of filaments during flip events, demonstrating the origin of the drift bias in non-uniform flows, and confirm this behavior from our simulations. Put together, these results shed light on the role of a local and dominant coupling between elasticity and viscous resistance in dictating long-term dynamics and transport of elastic fibers in confined flows.

Parallel computing for mobilities in periodic geometries

We examine methods for calculating the effective mobilities of molecules driven through periodic geometries in the context of particle-based simulation. The standard formulation of the mobility, based on the long-time limit of the mean drift velocity, is compared to a formulation based on the mean first-passage time of molecules crossing a single period of the system geometry. The equivalence of the two definitions is derived under weaker assumptions than similar conclusions obtained previously, requiring only that the state of the system at subsequent period crossings satisfy the Markov property. Approximate theoretical analyses of the computational costs of estimating these two mobility formulations via particle simulations suggest that the definition based on first-passage times may be substantially better suited to exploiting parallel computation hardware. This claim is investigated numerically on an example system modeling the passage of nanoparticles through the slit-well device. In this case, the traditional mobility formulation is found to perform best when the Péclet number is small, whereas the mean first-passage time formulation is found to converge much more quickly when the Péclet number is moderate or large. The results suggest that, given relatively modest access to modern GPU hardware, this alternative mobility formulation may be an order of magnitude faster than the standard technique for computing effective mobilities of biomolecules through periodic geometries.

Obstacle-induced giant jammed aggregation of active semiflexible filaments

Filaments driven by bound motor proteins and chains of self-propelled colloidal particles are a typical example of active polymers (APs). Due to deformability, APs exhibit very rich dynamic behaviors and collective assembling structures. Here, we are concerned with a basic question: how APs behave near a single obstacle? We find that, in the presence of a big single obstacle, the assembly of APs becomes a two-state system, i.e. APs either gather nearly completely together into a giant jammed aggregate (GJA) on the surface of the obstacle or distribute freely in space. No partial aggregation is observed. Such a complete aggregation/collection is unexpected since it happens on a smooth convex surface instead of, e.g., a concave wedge. We find that the formation of a GJA experiences a process of nucleation and the curves of the transition between the GJA and the non-aggregate state form hysteresis-like loops. Statistical analysis of massive data on the growing time, chirality and angular velocity of both the GJAs and the corresponding nuclei shows the strong random nature of the phenomenon. Our results provide new insights into the behavior of APs in contact with porous media and also a reference for the design and application of polymeric active materials.

Jammed solids with pins: Thresholds, force networks, and elasticity

The role of fixed degrees of freedom in soft or granular matter systems has broad applicability and theoretical interest. Here we address questions of the geometrical role that a scaffolding of fixed particles plays in tuning the threshold volume fraction and force network in the vicinity of jamming. Our two-dimensional simulated system consists of soft particles and fixed "pins," both of which harmonically repel overlaps. On the one hand, we find that many of the critical scalings associated with jamming in the absence of pins continue to hold in the presence of even dense pin latices. On the other hand, the presence of pins lowers the jamming threshold in a universal way at low pin densities and a geometry-dependent manner at high pin densities, producing packings with lower densities and fewer contacts between particles. The onset of strong lattice dependence coincides with the development of bond-orientational order. Furthermore, the presence of pins dramatically modifies the network of forces, with both unusually weak and unusually strong forces becoming more abundant. The spatial organization of this force network depends on pin geometry and is described in detail. Using persistent homology, we demonstrate that pins modify the topology of the network. Finally, we observe clear signatures of this developing bond-orientational order and broad force distribution in the elastic moduli which characterize the linear response of these packings to strain.

Compression-induced buckling of a semiflexible filament in two and three dimensions

The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. For filaments in the low-temperature limit and under a constant compressive force, Euler buckling theory predicts a sudden transition from a compressed to a bent state in these slender rods. In this paper, we use a mean-field theory to show that if a semiflexible chain is compressed at a finite temperature with a fixed end-to-end distance (permitting fluctuations in the compressive forces), it exhibits a continuous phase transition to a buckled state at a critical level of compression. We determine a quantitatively accurate prediction of the transverse position distribution function of the midpoint of the chain that indicates this transition. We find the mean compressive forces are non-monotonic as the extension of the filament varies, consistent with the observation that strongly buckled filaments are less able to bear an external load. We also find that for the fixed extension (isometric) ensemble, the buckling transition does not coincide with the local minimum of the mean force (in contrast to Euler buckling). We also show the theory is highly sensitive to fluctuations in length in two dimensions, and that the buckling transition can still be accurately recovered by accounting for those fluctuations. These predictions may be useful in understanding the behavior of filamentous biomolecules compressed by fluctuating forces, relevant in a variety of biological contexts.

Single-walled carbon nanotube reptation dynamics in submicron sized pores from randomly packed mono-sized colloids

Studying the Brownian motion of fibers and semi-flexible filaments in porous media is the key to understanding the transport and mechanical properties in a variety of systems. The motion of semi-flexible filaments in gel-like porous media including polymer networks and cell cytoskeleton has been studied theoretically and experimentally, whereas the motion of these materials in packed-colloid porous media, advanced foams, and rock-like systems has not been thoroughly studied. Here we use video microscopy to directly visualize the reptation and transport of intrinsically fluorescent, semiflexible, semiconducting single-walled carbon nanotubes (SWCNTs) in the sub-micron pores of packed colloids as fixed obstacles of packed-colloid porous media. By visualizing the filament motion and Brownian diffusion at different locations in the pore structures, we study how the properties of the environment, like the pore shape and pore structure of the porous media, affect SWCNT mobility. These results show that the porous media structure controls SWCNT reorientation during Brownian diffusion. In packed-colloid pores, SWCNTs diffuse along straight pores and bend across pores; conversely, in gel pores, SWCNTs consistently diffuse into curved pores, displaying a faster parallel motion. In both gel and packed-colloid porous media, SWCNT finite stiffness enhances SWCNT rotational diffusion and prevents jamming, allowing for inter-pore diffusion.

Chromatographic separation of active polymer-like worm mixtures by contour length and activity

The convective transport rate of polymers through confined geometries depends on their size, allowing for size-based separation of polymer mixtures (chromatography). Here, we investigate whether mixtures of active polymers can be separated in a similar manner based on their activity. We use thin, living Tubifex tubifex worms as a model system for active polymers and study the transport of these worms by an imposed flow through a channel filled with a hexagonal pillar array. The transport rate through the channel depends strongly on the degree of activity, an effect that we assign to the different distribution of conformations sampled by the worms depending on their activity. Our results demonstrate a unique way to sort mixtures of active polymers based on their activity and provide a versatile and convenient experimental system to investigate the hydrodynamics of active polymers.

Lévy Walks and Path Chaos in the Dispersal of Elongated Structures Moving across Cellular Vortical Flows

In cellular vortical flows, namely arrays of counterrotating vortices, short but flexible filaments can show simple random walks through their stretch-coil interactions with flow stagnation points. Here, we study the dynamics of semirigid filaments long enough to broadly sample the vortical field. Using simulation, we find a surprising variety of long-time transport behavior-random walks, ballistic transport, and trapping-depending upon the filament's relative length and effective flexibility. Moreover, we find that filaments execute Lévy walks whose diffusion exponents generally decrease with increasing filament length, until transitioning to Brownian walks. Lyapunov exponents likewise increase with length. Even completely rigid filaments, whose dynamics is finite dimensional, show a surprising variety of transport states and chaos. Fast filament dispersal is related to an underlying geometry of "conveyor belts." Evidence for these various transport states is found in experiments using arrays of counterrotating rollers, immersed in a fluid and transporting a flexible ribbon.