Search

The Reynolds equation from lubrication theory and the Stokes equations for low Reynolds number flows are distinct models for an incompressible fluid with negligible inertia. Here we investigate the sensitivity of the Reynolds equation to large gradients in the surface geometry. We present an analytic solution to the Reynolds equation in a piecewise-linear domain alongside a more general finite difference solution. For the Stokes equations, we use a finite difference solution for the biharmonic stream-velocity formulation. We compare the fluid velocity, pressure, and resistance for various step bearing geometries in the lubrication and Stokes limits. We find that the solutions to the Reynolds equation do not capture flow separation resulting from large cross-film pressure gradients. Flow separation and corner flow recirculation in step bearings are explored further; we consider the effect of smoothing large gradients in the surface geometry in order to recover limits under which the lubrication and Stokes approximations converge., Comment: 20 pages, 17 figures

Cell size varies between different cell types, and between different growth and osmotic conditions. However, the nuclear-to-cell volume ratio (N/C ratio) remains nearly constant. In this paper, we build on existing deterministic models of N/C ratio homeostasis and develop a simplified gene translation model to study the effect of stochasticity on the N/C ratio homeostasis. We solve the corresponding chemical master equation and obtain the mean and variance of the N/C ratio. We also use a Taylor expansion approximation to study the effects of the system size on the fluctuations of the N/C ratio. We then combine the translation model with a cell division model to study the effects of extrinsic noises from cell division on the N/C ratio. Our model demonstrates that the N/C ratio homeostasis is maintained when the stochasticity in cell growth is taken into account, that the N/C ratio is largely determined by the gene fraction of nuclear proteins, and that the fluctuations in the N/C ratio diminish as the system size increases., Comment: We revised the draft according to reviewers' suggestions. Results remain unchanged

In this paper, we model the configurations of a system of hard rods by viewing each rod in a cell formed by its neighbors. By minimizing the free energy in the model and performing molecular dynamics, where, in both cases, the shape of the cell is a free parameter, we obtain the equilibrium orientational order parameter, free energy and pressure of the system. Our model enables the calculation of anisotropic stresses exerted on the walls of the cell due to shape change of the rod in photoisomerization. These results are a key step towards understanding molecular shape change effects in photomechanical systems under illumination., Comment: 15 pages, 8 figures

We use a theoretical model to explore how fluid dynamics, in particular, the pressure gradient and wall shear stress in a channel, affect the deposition of particles flowing in a microfluidic network. Experiments on transport of colloidal particles in pressure-driven systems of packed beads have shown that at lower pressure drop, particles deposit locally at the inlet, while at higher pressure drop, they deposit uniformly along the direction of flow. We develop a mathematical model and use agent-based simulations to capture these essential qualitative features observed in experiments. We explore the deposition profile over a two-dimensional phase diagram defined in terms of the pressure and shear stress threshold, and show that two distinct phases exist. We explain this apparent phase transition by drawing an analogy to simple one-dimensional models of aggregation in which the phase transition is calculated analytically., Comment: 10 pages, 8 figures including Supplemental Material

The transport of particles in cells is influenced by the properties of intracellular networks they traverse while searching for localized target regions or reaction partners. Moreover, given the rapid turnover in many intracellular structures, it is crucial to understand how temporal changes in the network structure affect diffusive transport. In this work, we use network theory to characterize complex intracellular biological environments across scales. We develop an efficient computational method to compute the mean first passage times for simulating a particle diffusing along two-dimensional planar networks extracted from fluorescence microscopy imaging. We first benchmark this methodology in the context of synthetic networks, and subsequently apply it to live-cell data from endoplasmic reticulum tubular networks., Comment: 18 pages, 11 figures

The geometry of the blood vessel wall plays a regulatory role on the motion of red blood cells (RBCs). The overall topography of the vessel wall depends on many features, among which the endothelial lining of the endothelial surface layer (ESL) is an important one. The endothelial lining of vessel walls presents a large surface area for exchanging materials between blood and tissues. The ESL plays a critical role in regulating vascular permeability, hindering leukocyte adhesion as well as inhibiting coagulation during inflammation. Changes in the ESL structure are believed to cause vascular hyperpermeability and entrap immune cells during sepsis, which could significantly alter the vessel wall geometry and disturb interactions between RBCs and the vessel wall, including the wall-induced migration of RBCs and the thickening of a cell-free layer. To investigate the influence of the vessel wall geometry particularly changed by the ESL under various pathological conditions, such as sepsis, on the motion of RBCs, we developed two models to represent the ESL using the immersed boundary method in two dimensions. In particular, we used simulations to study how the lift force and drag force on a RBC near the vessel wall vary with different wall thickness, spatial variation, and permeability associated with changes in the vessel wall geometry. We find that the spatial variation of the wall has a significant effect on the wall-induced migration of the RBC for a high permeability, and that the wall-induced migration is significantly inhibited as the vessel diameter is increased.

We propose a coarse-grained theoretical model to capture the aging of microfluidic devices under different conditions including constant applied flow rate and constant applied pressure gradient. Microfluidic devices that sort cells by their deformability hold significant promise for medical applications. However, clogging in these microfluidic systems causes their properties to change over time and potentially limits their reliability. We compare the results of the coarse-grained model to those of stochastic simulations and to existing theoretical studies. Lastly, we apply the model to experimental data on the clogging of sickle red blood cells and discuss its wider applicability., Comment: 27 pages, 8 figures

We study the dynamics of membrane vesicle motor transport into dendritic spines, which are bulbous intracellular compartments in neurons that play a key role in transmitting signals between neurons. We consider the stochastic analog of the vesicle transport model in [Park and Fai, The Dynamics of Vesicles Driven Into Closed Constrictions by Molecular Motors. Bull. Math. Biol. 82, 141 (2020)]. The stochastic version, which may be considered as an agent-based model, relies mostly on the action of individual myosin motors to produce vesicle motion. To aid in our analysis, we coarse-grain this agent-based model using a master equation combined with a partial differential equation describing the probability of local motor positions. We confirm through convergence studies that the coarse-graining captures the essential features of bistability in velocity (observed in experiments) and waiting-time distributions to switch between steady-state velocities. Interestingly, these results allow us to reformulate the translocation problem in terms of the mean first passage time for a run-and-tumble particle moving on a finite domain with absorbing boundaries at the two ends. We conclude by presenting numerical and analytical calculations of vesicle translocation., Comment: 33 pages, 10 figures

Within this work we derive and analyse an expression for the free energy of a single-species system in the thermodynamic limit in terms of a generalised cavity volume, that is exact in general, and in principle applicable to systems across their entire range of density, as well as to particles within a general coordinate space. This provides a universal equation of state, and can thus relate the cavity volume to classical results, such as Mayer's cluster expansions. Through this we are able to provide some insight into the connections between cavity volume and free energy density, as well as their consequences. We use examples which permit explicit computations to further probe these results, reclaiming the exact results for a classical Tonks gas and providing a novel derivation of Onsager's free energy for a single species, isotropic system. Given the complexity of the problem we also provide a local lattice ansatz, exact in one dimension, with which we may approximate the cavity volume for hard sphere systems to provide an accurate equation of state in the cases of hard disks and spheres in both dilute regimes as well as beyond the freezing transition.

We study packings of hard spheres on lattices. The partition function, and therefore the pressure, may be written solely in terms of the accessible free volume, i.e. the volume of space that a sphere can explore without touching another sphere. We compute these free volumes using a leaky cell model, in which the accessible space accounts for the possibility that spheres may escape from the local cage of lattice neighbors. We describe how elementary geometry may be used to calculate the free volume exactly for this leaky cell model in two- and three-dimensional lattice packings and compare the results to the well-known Carnahan-Starling and Percus-Yevick liquid models. We provide formulas for the free volumes of various lattices and use the common tangent construction to identify several phase transitions between them in the leaky cell regime, indicating the possibility of coexistence in crystalline materials., Comment: 18 pages, 10 figures

Organelle size control is a fundamental question in biology that demonstrates the fascinating ability of cells to maintain homeostasis within their highly variable environments. Theoretical models describing cellular dynamics have the potential to help elucidate the principles underlying size control. Here, we perform a detailed study of the active disassembly model proposed in [Fai et al, Length regulation of multiple flagella that self-assemble from a shared pool of components, eLife, 8, (2019): e42599]. We construct a hybrid system which is shown to be well-behaved throughout the domain. We rule out the possibility of oscillations arising in the model and prove global asymptotic stability in the case of two flagella by the construction of a suitable Lyapunov function. Finally, we generalize the model to the case of arbitrary flagellar number in order to study olfactory sensory neurons, which have up to twenty cilia per cell. We show that our theoretical results may be extended to this case and explore the implications of this universal mechanism of size control., Comment: 25 pages, 6 figures

We study the dynamics of a model of membrane vesicle transport into dendritic spines, which are bulbous intracellular compartments in neurons driven by molecular motors. We reduce the lubrication model proposed in [Fai et al., Active elastohydrodynamics of vesicles in narrow, blind constrictions. Phys. Rev. Fluids, 2 (2017), 113601] to a fast-slow system, yielding an analytically and numerically tractable equation equivalent to the original model in the overdamped limit. The model's key parameters are the ratio of motors that prefer to push toward the head of the dendritic spine to the ratio of motors that prefer to push in the opposite direction. We perform a numerical bifurcation analysis in these parameters and find that steady-state vesicle velocities appear and disappear through several saddle-node bifurcations. This process allows us to identify the region of parameter space in which multiple stable velocities exist. We show by direct calculations that there can only be unidirectional motion for sufficiently close vesicle-to-spine diameter ratios. Our analysis predicts the critical vesicle-to-spine diameter ratio, at which there is a transition from unidirectional to bidirectional motion, consistent with experimental observations of vesicle trajectories in the literature., Comment: 34 pages, 9 figures

The growth, form, and division of prebiotic vesicles, membraneous bags of fluid of varying components and shapes is hypothesized to have served as the substrate for the origin of life. The dynamics of these out-of-equilibrium structures is controlled by physicochemical processes that include the intercalation of amphiphiles into the membrane, fluid flow across the membrane, and elastic deformations of the membrane. To understand prebiotic vesicular forms and their dynamics, we construct a minimal model that couples membrane growth, deformation, and fluid permeation, ultimately couched in terms of two dimensionless parameters that characterize the relative rate of membrane growth and the membrane permeability. Numerical simulations show that our model captures the morphological diversity seen in extant precursor mimics of cellular life, and might provide simple guidelines for the synthesis of these complex shapes from simple ingredients., Comment: Main text (6 pages, 3 figures) and Supplementary information (14 pages, 5 figures)

Fluid-resistance limited transport of vesicles through narrow constrictions is a recurring theme in many biological and engineering applications. Inspired by the motor-driven movement of soft membrane-bound vesicles into closed neuronal dendritic spines, here we study this problem using a combination of passive three-dimensional simulations and a simplified semi-analytical theory for active transport of vesicles that are forced through such constrictions by molecular motors. We show that the motion of these objects is characterized by two dimensionless quantities related to the geometry and the strength of forcing relative to the vesicle elasticity. We use numerical simulations to characterize the transit time for a vesicle forced by fluid pressure through a constriction in a channel, and find that relative to an open channel, transport into a blind end leads to the formation of an effective lubrication layer that strongly impedes motion. When the fluid pressure forcing is complemented by forces due to molecular motors that are responsible for vesicle trafficking into dendritic spines, we find that the competition between motor forcing and fluid drag results in multistable dynamics reminiscent of the real system. Our study highlights the role of non-local hydrodynamic effects in determining the kinetics of vesicular transport in constricted geometries.

Many biological examples of fluid-structure interaction, including the transit of red blood cells through the narrow slits in the spleen and the intracellular trafficking of vesicles into dendritic spines, involve the near-contact of elastic structures separated by thin layers of fluid. Motivated by such problems, we introduce an immersed boundary method that uses elements of lubrication theory to resolve thin fluid layers between immersed boundaries. We demonstrate 2nd-order accurate convergence for simple two-dimensional flows with known exact solutions to showcase the increased accuracy of this method compared to the standard immersed boundary method. Motivated by the phenomenon of wall-induced migration, we apply the lubricated immersed boundary method to simulate an elastic vesicle near a wall in shear flow. We also simulate the dynamics of a vesicle traveling through a narrow channel and observe the ability of the lubricated method to capture the vesicle motion on relatively coarse fluid grids., Comment: 35 pages, 11 figures

We consider a generalization of the classical Art Gallery Problem, where instead of a light source, the guards, called $k$-transmitters, model a wireless device with a signal that can pass through at most $k$ walls. We show it is NP-hard to compute a minimum cover of point 2-transmitters, point $k$-transmitters, and edge 2-transmitters in a simple polygon. The point 2-transmitter result extends to orthogonal polygons. In addition, we give necessity and sufficiency results for the number of edge 2-transmitters in general, monotone, orthogonal monotone, and orthogonal polygons., Comment: 15 pages, 11 figures