Publications

The macroscopic properties of polymeric fluids are inherited from the material properties of the fibers embedded in the solvent. The behavior of such passive fibers in flow has been of interest in a wide range of systems, including cellular mechanics, nutrient acquisition by diatom chains in the ocean, and industrial applications such as paper manufacturing. The rotational dynamics and shape evolution of fibers in shear depends upon the slenderness of the fiber and the non-dimensional “elasto-viscous” number that measures the ratio of the fluid’s viscous forces to the fiber’s elastic forces. For a small elasto-viscous number, the nearly-rigid fiber rotates in the shear, but when the elasto-viscous number reaches a threshold, buckling occurs. For even larger elasto-viscous numbers, there is a transition to a “snaking behavior” where the fiber remains aligned with the shear axis, but its ends curl in, in opposite directions. These experimentally-observed behaviors have recently been characterized computationally using slender-body theory and immersed boundary computations. However, classical experiments with nylon fibers and recent experiments with actin filaments have demonstrated that for even larger elasto-viscous numbers, multiple buckling sites and coiling can occur. Using a regularized Stokeslet framework coupled with a kernel independent fast multipole method, we present simulations that capture these complex fiber dynamics.

The Immersed Boundary (IB) method has been widely used to solve fluid-structure interaction problems, including those where the structure interacts with polymeric fluids. In this paper, we examine the convergence of one such scheme for a well known two-dimensional benchmark flow for the Oldroyd-B constitutive model, and we show that the traditional IB-based scheme fails to adequately capture the polymeric stress near to embedded boundaries. We analyze the reason for such failure, and we argue that this feature is not specific to the case study chosen, but a general feature of such methods due to lack of convergence in velocity gradients near interfaces. In order to remedy this problem, we build a different scheme for the Oldroyd-B system using the Immersed Boundary Smooth Extension (IBSE) scheme, which provides convergent viscous stresses near boundaries. We show that this modified scheme produces convergent polymeric stresses through the whole domain, including on embedded boundaries, and produces solutions in good agreement with known benchmarks.

We examine the interactions between actively rotating proteins moving in a membrane. Experimental evidence suggests that such rotor proteins, like the ATP synthases of the inner mitochondrial membrane, can arrange themselves into lattices. We show that crystallization is possible through a combination of hydrodynamic and repulsive interactions between the rotor proteins. In particular, hydrodynamic interactions induce rotational motion of the rotor protein assembly that, in the presence of repulsion, drives the system into a hexagonal lattice. The entire crystal rotates with an angular velocity which increases with motor density and decreases with lattice diameter - larger and sparser arrays rotate at a slower pace. The rotational interactions allow ensembles of proteins to sample configurations and reach an ordered steady state, which are inaccessible to the quenched nonrotational system. Rotational interactions thus act as a sort of temperature that removes disorder, except that actual thermal diffusion leads to expansion and loss of order. In contrast, the rotational interactions are bounded in space. Hence, once an ordered state is reached, it is maintained at all times.

In this work, we provide a solution to the problem of computing collision stress in particle-tracking simulations. First, a formulation for the collision stress between particles is derived as an extension of the virial stress formula to general-shaped particles with uniform or non-uniform mass density. Second, we describe a collision-resolution algorithm based on geometric constraint minimization which eliminates the stiff pairwise potentials in traditional methods. The method is validated with a comparison to the equation of state of Brownian spherocylinders. Then we demonstrate the application of this method in several emerging problems of soft active matter.

Many subcellular structures contain large numbers of cytoskeletal filaments. Such assemblies underlie much of cell division, motility, signaling, metabolism, and growth. Thus, understanding cell biology requires understanding the properties of networks of cytoskeletal filaments. While there are well established disciplines in biology dedicated to studying isolated proteins — their structure (Structural Biology) and behaviors (Biochemistry) — it is much less clear how to investigate, or even just describe, the structure and behaviors of collections of cytoskeletal filaments. One approach is to use methodologies from Mechanics and Soft Condensed Matter Physics, which have been phenomenally successful in the domains where they have been traditionally applied. From this perspective, collections of cytoskeletal filaments are viewed as materials, albeit very complex, ‘active’ materials, composed of molecules which use chemical energy to perform mechanical work. A major challenge is to relate these material level properties to the behaviors of the molecular constituents. Here we discuss this materials perspective and review recent work bridging molecular and network scale properties of the cytoskeleton, focusing on the organization of microtubules by dynein as an illustrative example.

An important class of fluid-structure problems involve the dynamics of ordered arrays of immersed, flexible fibers. While specialized numerical methods have been developed to study fluid-fiber systems, they become infeasible when there are many, rather than a few, fibers present, nor do these methods lend themselves to analytical calculation. Here, we introduce a coarse-grained continuum model, based on local-slender body theory, for elastic fibers immersed in a viscous Newtonian fluid. It takes the form of an anisotropic Brinkman equation whose skeletal drag is coupled to elastic forces. This model has two significant benefits: (1) the density effects of the fibers in a suspension become analytically manifest, and (2) it allows for the rapid simulation of dense suspensions of fibers in regimes inaccessible to standard methods. As a first validation, without fitting parameters, we achieve very reasonable agreement with 3D Immersed Boundary simulations of a bed of anchored fibers bent by a shear flow. Secondly, we characterize the effect of density on the relaxation time of fiber beds under oscillatory shear, and find close agreement to results from full numerical simulations. We then study buckling instabilities in beds of fibers, using our model both numerically and analytically to understand the role of fiber density and the structure of buckling transitions. We next apply our model to study the flow-induced bending of inclined fibers in a channel, as has been recently studied as a flow rectifier, examining the nature of the internal flows within the bed, and the emergence of inhomogeneous permeability. Finally, we extend the method to study a simple model of metachronal waves on beds of actuated fibers, as a model for ciliary beds. Our simulations reproduce qualitatively the pumping action of coordinated waves of compression through the bed.

The dynamics and deformations of immersed flexible fibers are at the heart of important industrial and biological processes, induce peculiar mechanical and transport properties in the fluids that contain them, and are the basis for novel methods of flow control. Here we focus on the low–Reynolds number regime where advances in studying these fiber–fluid systems have been especially rapid. On the experimental side, this is due to new methods of fiber synthesis, microfluidic flow control, and microscope-based tracking measurement techniques. Likewise, there have been continuous improvements in the specialized mathematical modeling and numerical methods needed to capture the interactions of slender flexible fibers with flows, boundaries, and each other.

In suspensions of microorganisms, pattern formation can arise from the interplay of chemotaxis and the fluid flows collectively generated by the organisms themselves. Here we investigate the resulting pattern formation in square and elongated domains in the context of two distinct models of locomotion in which the chemoattractant dynamics is fully coupled to the fluid flows and swimmer motion. Analyses for both models reveal an aggregative instability due to chemotaxis, independent of swimmer shape and type, and a hydrodynamic instability for “pusher” swimmers. We discuss the similarities and differences between the models. Simulations reveal a critical length scale of the swimmer aggregates and this feature can be utilized to stabilize swimmer concentration patterns into quasi-one-dimensional bands by varying the domain size. These concentration bands transition to traveling pulses under an external chemoattractant gradient, as observed in experiments with chemotactic bacteria.

A key unmet challenge in interpreting omics experiments is inferring biological meaning in the context of public functional genomics data. We developed a computational framework, Your Evidence Tailored Integration (YETI; http://yeti.princeton.edu/ ), which creates specialized functional interaction maps from large public datasets relevant to an individual omics experiment. Using this tailored integration, we predicted and experimentally confirmed an unexpected divergence in viral replication after seasonal or pandemic human influenza virus infection.

The 3D spatiotemporal organization of the human genome inside the cell nucleus remains a major open question in cellular biology. In the time between two cell divisions, chromatin - the functional form of DNA in cells - fills the nucleus in its uncondensed polymeric form. Recent in vivo imaging experiments reveal that the chromatin moves coherently, having displacements with long-ranged correlations on the scale of microns and lasting for seconds. To elucidate the mechanism(s) behind these motions, we develop a novel coarse-grained active-polymer model where chromatin is represented as a confined flexible chain acted upon by molecular motors, which perform work by exerting dipolar forces on the system. Numerical simulations of this model account for steric and hydrodynamic interactions as well as internal chain mechanics. These demonstrate that coherent motions emerge in systems involving extensile dipoles and are accompanied by large-scale chain reconfigurations and nematic ordering. Comparisons with experiments show good qualitative agreement and support the hypothesis that self-organizing long-ranged hydrodynamic couplings between chromatin-associated active motor proteins are responsible for the observed coherent dynamics.