Rational social distancing policy during epidemics with limited healthcare capacity.

Epidemics of infectious diseases posing a serious risk to human health have occurred throughout history. During recent epidemics there has been much debate about policy, including how and when to impose restrictions on behaviour. Policymakers must balance a complex spectrum of objectives, suggesting a need for quantitative tools. Whether health services might be 'overwhelmed' has emerged as a key consideration. Here we show how costly interventions, such as taxes or subsidies on behaviour, can be used to exactly align individuals' decision making with government preferences even when these are not aligned. In order to achieve this, we develop a nested optimisation algorithm of both the government intervention strategy and the resulting equilibrium behaviour of individuals. We focus on a situation in which the capacity of the healthcare system to treat patients is limited and identify conditions under which the disease dynamics respect the capacity limit. We find an extremely sharp drop in peak infections at a critical maximum infection cost in the government's objective function. This is in marked contrast to the gradual reduction of infections if individuals make decisions without government intervention. We find optimal interventions vary less strongly in time when interventions are costly to the government and that the critical cost of the policy switch depends on how costly interventions are.

Direct numerical simulations of a microswimmer in a viscoelastic fluid.

This study presents the application of the smoothed profile (SP) method to perform direct numerical simulations for the motion of both passive and active "squirming" particles in Newtonian and viscoelastic fluids. We found that fluid elasticity has a significant impact on both the transient behavior and the steady-state velocity of the particles. Specifically, we observe that the swirling flow generated by the squirmer's surface velocity significantly enhances their swimming speed as the Weissenberg number increases, regardless of the swimming type. Furthermore, we find that pushers outperform pullers in Oldroyd-B fluids, suggesting that the speed of a squirmer depends on the swimmer type. To understand the physical origin of the phenomenon of swirling flow enhancing the swimming speed, we investigate the velocity field and polymer conformation around non-swirling and swirling neutral squirmers in viscoelastic fluids. Our investigation reveals that the velocity field around the neutral swirling squirmers exhibits pusher-like extensional flow characteristics, as well as an asymmetric polymer conformation distribution, which gives rise to this increased propulsion. This is confirmed by the investigation of the force on a fixed squirmer, which revealed that the polymer stress, particularly its diagonal components, plays a critical role in enhancing the swimming speed of swirling squirmers in viscoelastic fluids. Additionally, our results demonstrate that the maximum swimming speeds of swirling squirmers occur at an intermediate value of the fluid viscosity ratio for all swimmer types. These findings have important implications for understanding the behavior of particles and micro-organisms in complex fluids.

Learning to swim efficiently in a nonuniform flow field.

Microswimmers can acquire information on the surrounding fluid by sensing mechanical queues. They can then navigate in response to these signals. We analyze this navigation by combining deep reinforcement learning with direct numerical simulations to resolve the hydrodynamics. We study how local and nonlocal information can be used to train a swimmer to achieve particular swimming tasks in a nonuniform flow field, in particular, a zigzag shear flow. The swimming tasks are (1) learning how to swim in the vorticity direction, (2) learning how to swim in the shear-gradient direction, and (3) learning how to swim in the shear-flow direction. We find that access to laboratory frame information on the swimmer's instantaneous orientation is all that is required in order to reach the optimal policy for tasks (1) and (2). However, information on both the translational and rotational velocities seems to be required to accomplish task (3). Inspired by biological microorganisms, we also consider the case where the swimmers sense local information, i.e., surface hydrodynamic forces, together with a signal direction. This might correspond to gravity or, for microorganisms with light sensors, a light source. In this case, we show that the swimmer can reach a comparable level of performance to that of a swimmer with access to laboratory frame variables. We also analyze the role of different swimming modes, i.e., pusher, puller, and neutral.

Rational social distancing in epidemics with uncertain vaccination timing.

During epidemics people may reduce their social and economic activity to lower their risk of infection. Such social distancing strategies will depend on information about the course of the epidemic but also on when they expect the epidemic to end, for instance due to vaccination. Typically it is difficult to make optimal decisions, because the available information is incomplete and uncertain. Here, we show how optimal decision-making depends on information about vaccination timing in a differential game in which individual decision-making gives rise to Nash equilibria, and the arrival of the vaccine is described by a probability distribution. We predict stronger social distancing the earlier the vaccination is expected and also the more sharply peaked its probability distribution. In particular, equilibrium social distancing only meaningfully deviates from the no-vaccination equilibrium course if the vaccine is expected to arrive before the epidemic would have run its course. We demonstrate how the probability distribution of the vaccination time acts as a generalised form of discounting, with the special case of an exponential vaccination time distribution directly corresponding to regular exponential discounting.

Mechanochemical subcellular-element model of crawling cells.

Constructing physical models of living cells and tissues is an extremely challenging task because of the high complexities of both intra- and intercellular processes. In addition, the force that a single cell generates vanishes in total due to the law of action and reaction. The typical mechanics of cell crawling involve periodic changes in the cell shape and in the adhesion characteristics of the cell to the substrate. However, the basic physical mechanisms by which a single cell coordinates these processes cooperatively to achieve autonomous migration are not yet well understood. To obtain a clearer grasp of how the intracellular force is converted to directional motion, we develop a basic mechanochemical model of a crawling cell based on subcellular elements with the focus on the dependence of the protrusion and contraction as well as the adhesion and de-adhesion processes on intracellular biochemical signals. By introducing reaction-diffusion equations that reproduce traveling waves of local chemical concentrations, we clarify that the chemical dependence of the cell-substrate adhesion dynamics determines the crawling direction and distance with one chemical wave. Finally, we also perform multipole analysis of the traction force to compare it with the experimental results. Our present work sheds light on how intracellular chemical reactions are converted to a directional cell migration under the force-free condition. Although the detailed mechanisms of actual cells are far more complicated than our simple model, we believe that this mechanochemical model is a good prototype for more realistic models.

Smoothed profile method for direct numerical simulations of hydrodynamically interacting particles.

A general method is presented for computing the motions of hydrodynamically interacting particles in various kinds of host fluids for arbitrary Reynolds numbers. The method follows the standard procedure for performing direct numerical simulations (DNS) of particulate systems, where the Navier-Stokes equation must be solved consistently with the motion of the rigid particles, which defines the temporal boundary conditions to be satisfied by the Navier-Stokes equation. The smoothed profile (SP) method provides an efficient numerical scheme for coupling the continuum fluid mechanics with the dispersed moving particles, which are allowed to have arbitrary shapes. In this method, the sharp boundaries between solid particles and the host fluid are replaced with a smeared out thin shell (interfacial) region, which can be accurately resolved on a fixed Cartesian grid utilizing a SP function with a finite thickness. The accuracy of the SP method is illustrated by comparison with known exact results. In the present paper, the high degree of versatility of the SP method is demonstrated by considering several types of active and passive particle suspensions.

Dynamics of a chiral swimmer sedimenting on a flat plate.

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the dynamics of a single squirmer with and without gravity to clarify its motion in the vicinity of a flat plate. In the absence of gravity and chirality, the usual dynamics of a squirmer near a wall are recovered. The introduction of chirality modifies the swimming motion of squirmers, adding a component of motion in the third direction. When sedimentation is considered, different dynamics emerge for different gravity strength regimes. In a moderate gravity regime, neutral squirmers and pullers eventually stop moving and reorient in a direction perpendicular to the plate; by contrast, pushers exhibit continuous motion in a tilted direction. In the strong gravity regime, all squirmers sediment and reorient perpendicular to the plate. In this study, chirality is introduced to model realistic microswimmers, and its crucial effects on the nature of squirmer trajectories, which change from straight to circular, are discussed.

Control of cell colony growth by contact inhibition.

Contact inhibition is a cell property that limits the migration and proliferation of cells in crowded environments. Here we investigate the growth dynamics of a cell colony composed of migrating and proliferating cells on a substrate using a minimal model that incorporates the mechanisms of contact inhibition of locomotion and proliferation. We find two distinct regimes. At early times, when contact inhibition is weak, the colony grows exponentially in time, fully characterised by the proliferation rate. At long times, the colony boundary moves at a constant speed, determined only by the migration speed of a single cell and independent of the proliferation rate. Further, the model demonstrates how cell-cell alignment speeds up colony growth. Our model illuminates how simple local mechanical interactions give rise to contact inhibition, and from this, how cell colony growth is self-organised and controlled on a local level.

Direct observation of the attachment behavior of hydrophobic colloidal particles onto a bubble surface.

The attachment of solid particles to the surface of immersed gas bubbles plays a fundamental role in surface science, and hence plays key roles in various engineering fields ranging from industrial separation processes to the fabrication of functional materials. However, detailed investigation from a microscopic view on how a single particle attaches to a bubble surface and how the particle properties affect the attachment behavior has been so far scarcely addressed. Here, we observed the attachment of a single particle to a bubble surface using a high-speed camera and systematically investigated the effects of the wettability and shape of particles. We found that hydrophobic particles abruptly "jumped into" the bubble while sliding down the bubble surface to eventually satisfy their static contact angles, the behavior of which induced a much stronger attachment to the bubble surface. Interestingly, the determinant factor for the attachment efficiency of spherical particles was not the wettability of the spherical particles but the location of the initial collision with the bubble surface. In contrast, the attachment efficiency of anisotropically-shaped particles was found to increase with the hydrophobicity caused by a larger contact area to the bubble surface. Last but not least, a simple formulation is suggested to recover the contact angle based on the jump-in behavior.

Spontaneous spatiotemporal ordering of shape oscillations enhances cell migration.

The migration of cells is relevant for processes such as morphogenesis, wound healing, and invasion of cancer cells. In order to move, single cells deform cyclically. However, it is not understood how these shape oscillations influence collective properties. Here we demonstrate, using numerical simulations, that the interplay of directed motion, shape oscillations, and excluded volume enables cells to locally "synchronize" their motion and thus enhance collective migration. Our model captures elongation and contraction of crawling ameboid cells controlled by an internal clock with a fixed period, mimicking the internal cycle of biological cells. We show that shape oscillations are crucial for local rearrangements that induce ordering of neighboring cells according to their internal clocks even in the absence of signaling and regularization. Our findings reveal a novel, purely physical mechanism through which the internal dynamics of cells influences their collective behavior, which is distinct from well known mechanisms like chemotaxis, cell division, and cell-cell adhesion.

Modeling the mechanosensitivity of fast-crawling cells on cyclically stretched substrates.

The mechanosensitivity of cells, which determines how they are able to respond to mechanical signals, is crucial for the functioning of biological systems. Experimentally, this is investigated by studying the reorientation of cells on cyclically stretched substrates. The reorientation depends on the type of cell and on the stretching protocol, but the mechanisms responsible for the response are still not completely understood. Here, we introduce a computational model for fast crawling cells on cyclically stretched substrates that accounts for the sub-cellular elements responsible for cell shape and motility. This includes the dynamics of the cell membrane, the actin cytoskeleton, and the focal adhesions with the stretching substrate. These processes evolve over characteristic time scales that can vary by orders of magnitude and naturally give rise to the frequency dependent reorientation observed experimentally. Depending on which processes are being probed by the stretching and on the type of coupling with the substrate, our simulations predict either no reorientation, a bi-stability in the parallel and perpendicular directions, or a complete reorientation in either the parallel or perpendicular direction. In particular, we show that an asymmetry in the adhesion dynamics during the loading and unloading phases of the stretching, whether it comes from the response of the cell itself or from the precise stretching protocol, can be used to selectively align the cells. Our results provide further evidence for the importance of focal adhesion dynamics in determining the mechanosensitive response of cells, as well as a way to interpret recent experiments.

Metallonanobelt: A Kinetically Stable Shape-Persistent Molecular Belt Prepared by Reversible Self-Assembly Processes.

A triptycene-based shape-persistent belt-shaped macrocycle, metallonanobelt, was synthesized by the self-assembly of 2,3,6,7-tetraaminotriptycene L and square planar Pd2+. The pentamer was selectively formed by the complexation of L with Pd2+ in the presence of the pillar[6]arene derivative P6 having triethylene glycol pendant as a template, whereas a mixture of a trimer, tetramer, and pentamer was obtained in the absence of the template. The pentamer was successfully isolated based on the solubility difference between the metallonanobelt and the template. It was also revealed that the isolated pentamer was remarkably stable in solutions such as acetonitrile or methanol thanks to the relatively inert planar chelate metal complex, [Pd( o-phenylenediamine)2] unit. Thus, we can handle the metallonanobelt almost as a static organic nanobelt that is synthesized covalently.

Reynolds-number-dependent dynamical transitions on hydrodynamic synchronization modes of externally driven colloids.

The collective dynamics of externally driven N_{p}-colloidal systems (1≤N_{p}≤4) in a confined viscous fluid have been investigated using three-dimensional direct numerical simulations with fully resolved hydrodynamics. The dynamical modes of collective particle motion are studied by changing the particle Reynolds number as determined by the strength of the external driving force and the confining wall distance. For a system with N_{p}=3, we found that at a critical Reynolds number a dynamical mode transition occurs from the doublet-singlet mode to the triplet mode, which has not been reported experimentally. The dynamical mode transition was analyzed in detail from the following two viewpoints: (1) spectrum analysis of the time evolution of a tagged particle velocity and (2) the relative acceleration of the doublet cluster with respect to the singlet particle. For a system with N_{p}=4, we found similar dynamical mode transitions from the doublet-singlet-singlet mode to the triplet-singlet mode and further to the quartet mode.

Field-induced dipolar attraction between like-charged colloids.

The field induced anisotropic interactions between like-charged colloidal particles is studied using direct numerical simulations, where the polarization of the electric double layer is explicitly computed under external AC electric fields. These interactions are found to depend on the magnitude E0 and frequency ω of the applied field, as well as the zeta potential, the Debye length, and the relative orientation of the particles. We also determined the range of E0 and ω over which a dipolar attraction is induced between a pair of like-charged colloids. Finally, we performed simulations for systems of six and twelve colloidal particles to study the stability of pear-chain-like configurations.

Do hydrodynamically assisted binary collisions lead to orientational ordering of microswimmers?

We have investigated the onset of collective motion in systems of model microswimmers, by performing a comprehensive analysis of the binary collision dynamics using three-dimensional direct numerical simulations (DNS) with hydrodynamic interactions. From this data, we have constructed a simplified binary collision model (BCM) which accurately reproduces the collective behavior obtained from the DNS for most cases. Thus, we show that global alignment can mostly arise solely from binary collisions. Although the agreement between both models (DNS and BCM) is not perfect, the parameter range in which notable differences appear is also that for which strong density fluctuations are present in the system (where pseudo-sound mound can be observed (N. Oyama et al., Phys. Rev. E 93, 043114 (2016))).

Collective motion of cells crawling on a substrate: roles of cell shape and contact inhibition.

Contact inhibition plays a crucial role in cell motility, wound healing, and tumour formation. By mimicking the mechanical motion of cells crawling on a substrate, we constructed a minimal model of migrating cells that naturally gives rise to contact inhibition of locomotion (CIL). The model cell consists of two disks, a front disk (a pseudopod) and a back disk (cell body), which are connected by a finite extensible spring. Despite the simplicity of the model, the collective behaviour of the cells is highly non-trivial and depends on both the shape of the cells and whether CIL is enabled. Cells with a small front disk (i.e., a narrow pseudopod) form immobile colonies. In contrast, cells with a large front disk (e.g., a lamellipodium) exhibit coherent migration without any explicit alignment mechanism in the model. This result suggests that crawling cells often exhibit broad fronts because this helps facilitate alignment. After increasing the density, the cells develop density waves that propagate against the direction of cell migration and finally stop at higher densities.

Direct numerical simulation of an arbitrarily shaped particle at a fluidic interface.

A consistent formulation is presented for the direct numerical simulation of an arbitrarily shaped colloidal particle at a deformable fluidic interface. The rigid colloidal particle is decomposed into a collection of solid spherical beads and the three-phase boundaries are replaced with smoothly spreading interfaces. The major merit of the present formulation lies in the ease with which the geometrical decomposition of the colloidal particle is implemented, yet allows the dynamic simulation of intricate three-dimensional colloidal shapes in a binary fluid. The dynamics of a rodlike, a platelike, and a ringlike particle are presently tested. It is found that platelike particles attach more rapidly to a fluidic interface and are subsequently harder to dislodge when subject to an external force. Using the Bond number, i.e., the ratio of the gravitational force to the reference capillary force, a spherical particle with equal affinity for the two fluids breaks away from a fluidic interface at the critical value Bo=0.75. This value is in line with our numerical experiments. It is here shown that a plate and a ring of equivalent masses detach at greater critical Bond numbers approximately equal to Bo=1.3. Results of this study will find applications in the stabilization of emulsions by colloids and in the recovery of colloidal particles by rising bubbles.

Purely hydrodynamic origin for swarming of swimming particles.

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the collective motion of model swimmers in bulk and confinement. Calculating the dynamic structure factor, we clarified that the swarming in bulk systems can be understood as a pseudoacoustic mode. Under confinement between flat parallel walls, this pseudoacoustic mode leads to a traveling wavelike motion. This swarming behavior is due purely to the hydrodynamic interactions between the swimmers and depends strongly on the type and strength of swimming (i.e., pusher or puller).

Mutual information reveals multiple structural relaxation mechanisms in a model glass former.

Among the key challenges to our understanding of solidification in the glass transition is that it is accompanied by little apparent change in structure. Recently, geometric motifs have been identified in glassy liquids, but a causal link between these motifs and solidification remains elusive. One 'smoking gun' for such a link would be identical scaling of structural and dynamic lengthscales on approaching the glass transition, but this is highly controversial. Here we introduce an information theoretic approach to determine correlations in displacement for particle relaxation encoded in the initial configuration of a glass-forming liquid. We uncover two populations of particles, one inclined to relax quickly, the other slowly. Each population is correlated with local density and geometric motifs. Our analysis further reveals a dynamic lengthscale similar to that associated with structural properties, which may resolve the discrepancy between structural and dynamic lengthscales.

Fungal hemicellulose-degrading enzymes cause physical property changes concomitant with solubilization of cell wall polysaccharides.

MAIN CONCLUSION: Physical properties of wheat coleoptile segments decreased after treatment with hemicellulose-degrading enzymes, indicating that hemicellulosic polysaccharides function to control the strength of primary cell walls. Changes in the physical properties of plant cell walls, a viscoelastic structure, are thought to be one of the growth-limiting factors for plants and one of the infection-affecting factors for fungi. To study the significance of hemicellulosic polysaccharides that form cross-bridges between cellulose microfibrils in controlling cell wall strength in monocot plants, the effects of hemicellulose degradation by recombinant Magnaporthe oryzae xylanase and 1,3-1,4-ß-glucanase, and recombinant Aspergillus oryzae xyloglucanase on the physical properties and polysaccharide solubilization were investigated using wheat (Triticum aestivum L.) coleoptiles. Treatments with xylanase or 1,3-1,4-ß-glucanase significantly decreased the viscosity and elasticity of wheat coleoptile segments. In addition, xyloglucanase treatment slightly decreased the viscoelasticity. Furthermore, 1,3-1,4-ß-glucan polymer was solubilized during hydrolysis with xylanase and xyloglucanase, even though neither enzyme had hydrolytic activity towards 1,3-1,4-ß-glucan. These results suggest that xylan and xyloglucan interact with 1,3-1,4-ß-glucan and that the composites and hemicellulosic polysaccharides form inter-molecular bridges. Degradation of these bridges causes decreases in the physical properties, resulting in increased extensibility of the cell walls. These findings provide a testable model in which wheat coleoptile cell walls are loosened by the degradation of hemicellulosic polysaccharides and hemicellulose-degrading enzymes play a significant role in loosening the walls during fungal infection.