Source author record

Chiu Fan Lee

Chiu Fan Lee 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

25works
14topics
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

25 published item(s)

preprint2022arXiv

An Infinite Set of Integral Formulae for Polar, Nematic, and Higher Order Structures at the Interface of Motility-Induced Phase Separation

Motility-induced phase separation (MIPS) is a purely non-equilibrium phenomenon in which self-propelled particles phase separate without any attractive interactions. One surprising feature of MIPS is the emergence of polar, nematic, and higher order structures at the interfacial region, whose underlying physics remains poorly understood. Starting with a model of MIPS in which all many-body interactions are captured by an effective speed function and an effective pressure function that depend solely on the local particle density, I derive analytically an infinite set of integral formulae (IF) for the ordering structures at the interface. I then demonstrate that half of these IF are in fact exact for a wide class of active Brownian particle systems. Finally, I test the IF by applying them to numerical data from direct particle dynamics simulation and find that all the IF remain valid to a great extent.

preprint2022arXiv

Diversity of phase transitions and phase separations in active fluids

Active matter is not only indispensable to our understanding of diverse biological processes, but also provides a fertile ground for discovering novel physics. Many emergent properties impossible for equilibrium systems have been demonstrated in active systems. These emergent features include motility-induced phase separation, long-ranged ordered (collective motion) phase in two dimensions, and order-disorder phase co-existences (banding and reverse-banding regimes). Here, we unify these diverse phase transitions and phase co-existences into a single formulation based on generic hydrodynamic equations for active fluids. We also reveal a novel co-moving co-existence phase and a putative novel critical point.

preprint2022arXiv

Model of inverse bleb growth explains giant vacuole dynamics during cell mechanoadaptation

Cells can withstand hostile environmental conditions manifest as large mechanical forces such as pressure gradients and/or shear stresses by dynamically changing their shape. Such conditions are realized in the Schlemm's canal of the eye where endothelial cells that cover the inner vessel wall are subjected to the hydrodynamic pressure gradients exerted by the aqueous humor outflow. These cells form fluid-filled dynamic outpouchings of their basal membrane called \textit{giant vacuoles}. The inverse of giant vacuoles are reminiscent of cellular blebs, extracellular cytoplasmic protrusions triggered by local temporary disruption of the contractile actomyosin cortex. Inverse blebbing has been first observed experimentally during sprouting angiogenesis, but its underlying physical mechanisms are poorly understood. Here, we identify giant vacuole formation as inverse blebbing and formulate a biophysical model of this process. Our model elucidates how cell membrane mechanical properties affect the morphology and dynamics of giant vacuoles and predicts coarsening akin to Ostwald ripening between multiple invaginating vacuoles. Our results are in qualitative agreement with observations from the formation of giant vacuoles during perfusion experiments. Our model not only elucidates the biophysical mechanisms driving inverse blebbing and giant vacuole dynamics, but also identifies universal features of the cellular response to pressure loads that are relevant to many experimental contexts.

preprint2021arXiv

Polar Fluctuations Lead to Extensile Nematic Behavior in Confluent Tissues

How can a collection of motile cells, each generating contractile nematic stresses in isolation, become an extensile nematic at the tissue-level? Understanding this seemingly contradictory experimental observation, which occurs irrespective of whether the tissue is in the liquid or solid states, is not only crucial to our understanding of diverse biological processes, but is also of fundamental interest to soft matter and many-body physics. Here, we resolve this cellular to tissue level disconnect in the small fluctuation regime by using analytical theories based on hydrodynamic descriptions of confluent tissues, in both liquid and solid states. Specifically, we show that a collection of microscopic constituents with no inherently nematic extensile forces can exhibit active extensile nematic behavior when subject to polar fluctuating forces. We further support our findings by performing cell level simulations of minimal models of confluent tissues.

preprint2021arXiv

Uncovering novel phase transitions in dense dry polar active fluids using a lattice Boltzmann method

The dynamics of dry active matter have implications for a diverse collection of biological phenomena spanning a range of length and time scales, such as animal flocking, cell tissue dynamics, and swarming of inserts and bacteria. Uniting these systems are a common set of symmetries and conservation laws, defining dry active fluids as a class of physical system. Many interesting behaviours have been observed at high densities, which remain difficult to simulate due to the computational demand. Here, we show how two-dimensional dry active fluids in a dense regime can be studied using a simple modification of the lattice Boltzmann method. We apply our method on a model that exhibits motility-induced phase separation, and an active model with contact inhibition of locomotion, which has relevance to collective cell migration. For the latter, we uncover multiple novel phase transitions: two first-order and one potentially critical. We further support our simulation results with an analytical treatment of the hydrodynamic equations obtained via a Chapman-Enskog coarse-graining procedure.

preprint2020arXiv

A novel nonequilibrium state of matter: a $d=4-ε$ expansion study of Malthusian flocks

We show that "Malthusian flocks" -- i.e., coherently moving collections of self-propelled entities (such as living creatures) which are being "born" and "dying" during their motion -- belong to a new universality class in spatial dimensions $d>2$. We calculate the universal exponents and scaling laws of this new universality class to $O(ε)$ in a $d=4-ε$ expansion, and find these are different from the "canonical" exponents previously conjectured to hold for "immortal" flocks (i.e., those without birth and death) and shown to hold for incompressible flocks with spatial dimensions in the range of $2 < d \leq 4$. We also obtain a universal amplitude ratio relating the damping of transverse and longitudinal velocity and density fluctuations in these systems. Furthermore, we find a universal separatrix in real (${\bf r}$) space between two regions in which the equal time density correlation $\langleδρ({\bf r}, t)δρ(0, t)\rangle$ has opposite signs. Our expansion should be quite accurate in $d=3$, allowing precise quantitative comparisons between our theory, simulations, and experiments.

preprint2020arXiv

Active Lévy Matter: Anomalous Diffusion, Hydrodynamics and Linear Stability

Anomalous diffusion, manifest as a nonlinear temporal evolution of the position mean square displacement, and/or non-Gaussian features of the position statistics, is prevalent in biological transport processes. Likewise, collective behavior is often observed to emerge spontaneously from the mutual interactions between constituent motile units in biological systems. Examples where these phenomena can be observed simultaneously have been identified in recent experiments on bird flocks, fish schools and bacterial swarms. These results pose an intriguing question, which cannot be resolved by existing theories of active matter: How is the collective motion of these systems affected by the anomalous diffusion of the constituent units? Here, we answer this question for a microscopic model of active Lévy matter -- a collection of active particles that perform superdiffusion akin to a Lévy flight and interact by promoting polar alignment of their orientations. We present in details the derivation of the hydrodynamic equations of motion of the model, obtain from these equations the criteria for a disordered or ordered state, and apply linear stability analysis on these states at the onset of collective motion. Our analysis reveals that the disorder-order phase transition in active Lévy matter is critical, in contrast to ordinary active fluids where the phase transition is, instead, first-order. Correspondingly, we estimate the critical exponents of the transition by finite size scaling analysis and use these numerical estimates to relate our findings to known universality classes. These results highlight the novel physics exhibited by active matter integrating both anomalous diffusive single-particle motility and inter-particle interactions.

preprint2020arXiv

Hydrodynamics of Active Lévy Matter

Collective motion is often modeled within the framework of active fluids, where the constituent active particles, when interactions with other particles are switched off, perform normal diffusion at long times. However, in biology, single-particle superdiffusion and fat-tailed displacement statistics are also widespread. The collective properties of interacting systems exhibiting such anomalous diffusive dynamics, which we call active Lévy matter, cannot be captured by current active fluid theories. Here, we formulate a hydrodynamic theory of active Lévy matter by coarse-graining a microscopic model of aligning polar active particles that perform superdiffusion akin to Lévy flights. Applying a linear stability analysis on the hydrodynamic equations at the onset of collective motion, we find that, in contrast to its conventional counterpart, the order-disorder transition can become critical. We then estimate the corresponding critical exponents by finite size scaling analysis of numerical simulations. Our work highlights the novel physics in active matter that integrates both anomalous diffusive motility and inter-particle interactions.

preprint2020arXiv

Moving, reproducing, and dying beyond Flatland: Malthusian flocks in dimensions $d>2$

We show that "Malthusian flocks" -- i.e., coherently moving collections of self-propelled entities (such as living creatures) which are being "born" and "dying" during their motion -- belong to a new universality class in spatial dimensions $d>2$. We calculate the universal exponents and scaling laws of this new universality class to $O(ε)$ in an $ε=4-d$ expansion, and find these are different from the "canonical" exponents previously conjectured to hold for "immortal" flocks (i.e., those without birth and death) and shown to hold for incompressible flocks in $d>2$. Our expansion should be quite accurate in $d=3$, allowing precise quantitative comparisons between our theory, simulations, and experiments.

preprint2016arXiv

Comment on "Anomalous Discontinuity at the Percolation Critical Point of Active Gels"

In their recent work Sheinman et al. [Phys. Rev. Lett. 114, 098104 (2015)] introduce a variation of percolation which they call no-enclaves percolation (NEP). The main claims are 1) the salient physics captured in NEP is closer to what happens experimentally; 2) The Fisher exponent of NEP, is $τ=1.82(1)$; 3) Due to the different Fisher exponent, NEP constitutes a universality class distinct from random percolation (RP). While we fully agree with 1) and found NEP to be a very interesting variation of random percolation, we disagree with 2) and 3). We will demonstrate that $τ$ is exactly $2$, directly derivable from RP, and thus there is no foundation of a new universality class.

preprint2016arXiv

Percolation with trapping mechanism drives active gels to the critically connected state

Cell motility and tissue morphogenesis depend crucially on the dynamic remodelling of actomyosin networks. An actomyosin network consists of an actin polymer network connected by crosslinker proteins and motor protein myosins that generate internal stresses on the network. A recent discovery shows that for a range of experimental parameters, actomyosin networks contract to clusters with a power-law size distribution [Alvarado J. et al. (2013) Nature Physics 9 591]. Here, we argue that actomyosin networks can exhibit robust critical signature without fine-tuning because the dynamics of the system can be mapped onto a modified version of percolation with trapping (PT), which is known to show critical behaviour belonging to the static percolation universality class without the need of fine-tuning of a control parameter. We further employ our PT model to generate experimentally testable predictions.

preprint2016arXiv

Surprising mappings of 2D polar active fluids to 2D soap and 1D sandblasting

Active fluids and growing interfaces are two well-studied but very different non-equilibrium systems. Each exhibits non-equilibrium behavior quite different from that of their equilibrium counterparts. Here we demonstrate a surprising connection between these two: the ordered phase of incompressible polar active fluids in two spatial dimensions without momentum conservation, and growing one-dimensional interfaces (that is, the 1+1-dimensional Kardar-Parisi-Zhang equation), in fact belong to the same universality class. This universality class also includes two equilibrium systems : two-dimensional smectic liquid crystals, and a peculiar kind of constrained two-dimensional ferromagnet. We use these connections to show that two-dimensional incompressible flocks are robust against fluctuations, and exhibit universal long-ranged, anisotropic spatio-temporal correlations of those fluctuations. We also thereby determine the exact values of the anisotropy exponent $ζ$ and the roughness exponents $χ_{_{x,y}}$ that characterize these correlations.

preprint2014arXiv

Critical Phenomenon of the Order-Disorder Transition in Incompressible Flocks

We study incompressible systems of motile particles with alignment interactions. Unlike their compressible counterparts, in which the order-disorder (i.e., moving to static) transition, tuned by either noise or number density, is discontinuous, in incompressible systems this transition can be continuous, and belongs to a new universality class. We calculate the critical exponents to $O(ε)$in an $ε=4-d$ expansion, and derive two exact scaling relations. This is the first analytic treatment of a phase transition in a new universality class in an active system.

preprint2014arXiv

Thermal breakage of a semiflexible polymer: Breakage profile and rate

Understanding fluctuation-induced breakages in polymers has important implications for basic and applied sciences. Here I present for the first time an analytical treatment of the thermal breakage problem of a semi-flexible polymer model that is asymptotically exact in the low temperature and high friction limits. Specifically, I provide analytical expressions for the breakage propensity and rate, and discuss the generalities of the results and their relevance to biopolymers.

preprint2014arXiv

Thermodynamic costs of information processing in sensory adaption

Biological sensory systems react to changes in their surroundings. They are characterized by fast response and slow adaptation to varying environmental cues. Insofar as sensory adaptive systems map environmental changes to changes of their internal degrees of freedom, they can be regarded as computational devices manipulating information. Landauer established that information is ultimately physical, and its manipulation subject to the entropic and energetic bounds of thermodynamics. Thus the fundamental costs of biological sensory adaptation can be elucidated by tracking how the information the system has about its environment is altered. These bounds are particularly relevant for small organisms, which unlike everyday computers operate at very low energies. In this paper, we establish a general framework for the thermodynamics of information processing in sensing. With it, we quantify how during sensory adaptation information about the past is erased, while information about the present is gathered. This process produces entropy larger than the amount of old information erased and has an energetic cost bounded by the amount of new information written to memory. We apply these principles to the E. coli's chemotaxis pathway during binary ligand concentration changes. In this regime, we quantify the amount of information stored by each methyl group and show that receptors consume energy in the range of the information-theoretic minimum. Our work provides a basis for further inquiries into more complex phenomena, such as gradient sensing and frequency response.

preprint2013arXiv

Active particles under confinement: Aggregation at the wall and gradient formation inside a channel

I study the confinement-induced aggregation phenomenon in a minimal model of self-propelled particles inside a channel. Starting from first principles, I derive a set of equations that govern the density profile of such a system at the steady-state, and calculate analytically how the aggregation at the walls varies with the physical parameters of the system. I also investigate how the gradient of the particle density varies if the inside of the channel is partitioned into two regions within which the active particles exhibit distinct levels of fluctuations in their directions of travel.

preprint2013arXiv

Fluctuations and the role of collision duration in reaction-diffusion systems

In a reaction-diffusion system, fluctuations in both diffusion and reaction events, have important effects on the steady-state statistics of the system. Here, we argue through extensive lattice simulations, mean-field type arguments, and the Doi-Peliti formalism that the collision duration statistics -- i.e., the time two particles stay together in a lattice site -- plays a leading role in determining the steady state of the system. We obtain approximate expressions for the average densities of the chemical species and for the critical diffusion coefficient required to sustain the reaction.

preprint2012arXiv

Length distribution of stiff, self-assembled polymers at thermal equilibrium

We investigate the length distribution of self-assembled, long and stiff polymers at thermal equilibrium. Our analysis is based on calculating the partition functions of stiff polymers of variable lengths in the elastic regime. Our conclusion is that the length distribution of this self-assembled system follows closely the exponential distribution, except at the short length limit. We then discuss the implications of our results on the experimentally observed length distributions in amyloid fibrils.

preprint2011arXiv

Singular perturbation analysis of a reduced model for collective motion: A renormalization group approach

In a system of noisy self-propelled particles with interactions that favor directional alignment, collective motion will appear if the density of particles is beyond a critical density. Starting with a reduced model for collective motion, we determine how the critical density depends on the form of the initial perturbation. Specifically, we employ a renormalization-group improved perturbative method to analyze the model equations, and show analytically, up to first order in the perturbation parameter, how the critical density is modified by the strength of the initial angular perturbation in the system.

preprint2010arXiv

Fluctuation-induced collective motion: A single-particle density analysis

In a system of noisy self-propelled particles with interactions that favor directional alignment, collective motion will appear if the density of particles increases beyond a certain threshold. In this paper, we argue that such a threshold may depend also on the profiles of the perturbation in the particle directions. Specifically, we perform mean-field, linear stability, perturbative and numerical analyses on an approximated form of the Fokker-Planck equation describing the system. We find that if an angular perturbation to an initially homogeneous system is large in magnitude and highly localized in space, it will be amplified and thus serves as an indication of the onset of collective motion. Our results also demonstrate that high particle speed promotes collective motion.

preprint2010arXiv

Predicting rare events in chemical reactions: application to skin cell proliferation

In a well-stirred system undergoing chemical reactions, fluctuations in the reaction propensities are approximately captured by the corresponding chemical Langevin equation. Within this context, we discuss in this work how the Kramers escape theory can be used to predict rare events in chemical reactions. As an example, we apply our approach to a recently proposed model on cell proliferation with relevance to skin cancer [P.B. Warren, Phys. Rev. E {\bf 80}, 030903 (2009)]. In particular, we provide an analytical explanation for the form of the exponential exponent observed in the onset rate of uncontrolled cell proliferation.

preprint2009arXiv

Network Automata: Coupling structure and function in real-world networks

We introduce Network Automata, a framework which couples the topological evolution of a network to its structure. It is useful for dealing with networks in which the topology evolves according to some specified microscopic rules and, simultaneously, there is a dynamic process taking place on the network that both depends on its structure but is also capable of modifying it. It is a generic framework for modeling systems in which network structure, dynamics, and function are interrelated. At the practical level, this framework allows for easy implementation of the microscopic rules involved in such systems. To demonstrate the approach, we develop a class of simple biologically inspired models of fungal growth.

preprint2009arXiv

Self-assembly of protein amyloid: a competition between amorphous and ordered aggregation

Protein aggregation in the form of amyloid fibrils has important biological and technological implications. Although the self-assembly process is highly efficient, aggregates not in the fibrillar form would also occur and it is important to include these disordered species when discussing the thermodynamic equilibrium behavior of the system. Here, we initiate such a task by considering a mixture of monomeric proteins and the corresponding aggregates in the disordered form (micelles) and in the fibrillar form (amyloid fibrils). Starting with a model on the respective binding free energies for these species, we calculate their concentrations at thermal equilibrium. We then discuss how the incorporation of the disordered structure furthers our understanding on the various amyloid promoting factors observed empirically, and on the kinetics of fibrilization.

preprint2008arXiv

Density functional methods for polymers: a coil-globule transition case study

We consider a free energy functional on the monomer density function that is suitable for the study of coil-globule transition. We demonstrate, with explicitly stated assumptions, why the entropic contribution is in the form of the Kullback-Leibler distance, and that the energy contribution is given by two-body and three-body terms. We then solve for the free energy analytically on a set of trial density functions, and reproduce de Gennes' classical theory on polymer coil-globule transition. We then discuss how our formalism can be applied to study polymer dynamics from the perspective of dynamical density function theory.

preprint2005arXiv

Renormalization scheme for a multi-qubit-network

We present a renormalization scheme which simplifies the dynamics of an important class of interacting multi-qubit systems. We show that a wide class of M+1 qubit systems can be reduced to an equivalent n+1 qubit system with n equal to, or greater than, 2, for any M. Our renormalization scheme faithfully reproduces the overall dynamics of the original system including the entanglement properties. In addition to its direct application to atom-cavity and nanostructure systems, the formalism offers insight into a variety of situations ranging from decoherence due to a spin-bath with its own internal entanglement, through to energy transfer processes in organic systems such as biological photosynthetic units.