SciPost Submission Page
Nonexistence of motility induced phase separation transition in one dimension
by Indranil Mukherjee, Adarsh Raghu, P. K. Mohanty
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Pradeep Kumar Mohanty 
Submission information  

Preprint Link:  scipost_202209_00024v1 (pdf) 
Date submitted:  20220912 15:45 
Submitted by:  Mohanty, Pradeep Kumar 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approach:  Theoretical 
Abstract
We introduce and study a model of hardcore particles obeying runandtumble dynamics on a onedimensional lattice, where particles run in either +ve or ve xdirection with an effective speed v and tumble (change their direction of motion) with a constant rate ω. We show that the coarsegrained dynamics of the system can be mapped to a beadsinurn model called misanthrope process where particles are identified as urns and vacancies as beads that hop to a neighbouring urn situated in the direction opposite to the current. The hop rate is same as the magnitude of the particle current;we calculate it analytically for two particle system and show that it does not satisfy the criteria required for a phase separation transition. Nonexistence of phase separation in this model, where tumbling dynamics is rather restricted, necessarily imply that motility induced phase separation transition can not occur in other models in one dimension with unconditional tumbling.
Current status:
Reports on this Submission
Report #4 by Anonymous (Referee 4) on 20221028 (Invited Report)
 Cite as: Anonymous, Report on arXiv:scipost_202209_00024v1, delivered 20221028, doi: 10.21468/SciPost.Report.5999
Report
This is an interesting paper that addresses the question of whether motilityinduced phase separation (MIPS) can occur in a simple model of runandtumble particles in one dimension. By mapping the dynamics to an urn model, coarsegraining and then solving with a matrixproduct ansatz, the authors conclude that MIPS is not possible in this system. This corroborates certain existing findings, for example the work of Ref [14] which indicates that the domain size grows continuously as the tumble rate decreases, i.e., without any nonanalyticities of the type that might be associated with a phase transition. There may however be some regimes, such as that considered in Ref [27], where some kind of phase separation may occur, which are briefly discussed at the end of the manuscript.
I find this to be a timely contribution that address a question of current interest in the active matter / nonequilibrium statistical physics community. A particular strength is that it provides some analytical results for the manybody problem which, although not exact, go some way to improving our understanding of these driven systems.
Unfortunately, I found aspects of the presentation rather too terse for the manuscript to be easily followable, and I say this as someone who has worked extensively on 1d systems driven far from equilibrium. In particular the following sections of the manuscript would benefit from a considerable expansion to make the presentation digestible:
1. It took me some time to realise that the third update in (1) allows only the lefthand particle to change direction ("tumble"), and then only when the right neighbour is moving the same way. I think it would help first of all to write out the $\omega_+$ and $\omega_$ rules separately, and furthermore to state clearly in the text which updates are allowed and not allowed.
2. Related to 1, I think some discussion of the rationale for this set of rules is required. I understand that this allows analytical progress later in the manuscript, but did not follow at what point this becomes important. Furthermore, one may worry that the rule breaks CP symmetry for those cases where $r_+=r_$ for each of the transition rates $r\in\{p,q,\omega\}$ and such symmetry might be expected. I understand that setting certain tumbling rates to zero moves in the direction of making the system more likely to phase separate, so I doubt this affects the final conclusion, but one does feel nervous about removing a symmetry from a model that might initially be present.
3. I did not follow the setup of the coarsegrained model at all. The exact urn model has particles hopping in both directions, whereas the coarsegrained model seems to have them moving only to the right (or perhaps this is a misleading impression given by Fig 1c)? I did not see how the hop rates $u(m)$ are defined, whether these depend only on the occupancy of one site (as suggested by Fig 1c) or two sites (as suggested by the text). I really could not understand why some of the domains of size $m$ in Fig 1b overlap, or what determines the start and end of a domain. This entire mapping needs to be laid out completely and methodically.
4. In Fig 2, where hop rates are obtained from numerical simulations, it was unclear what was being simulated here. Is it the full dynamics of the original model, or one of the derivatives described in the main text.
5. I would like to see a bit more discussion of the comparison with Ref 27. Although the latter relies on a specific limit, I understood that the current analysis applied to any parameter combination so one ought be able to take the same scaling in the present work to compare more directly. Moreover, the $1/L$ separation of the run and tumbling rates I think corresponds to the scaling limit considered in Ref 15, which makes me wonder if something different happens in this scaling limit and whether this can be seen from the present analysis.
I have recorded these as "minor revisions" as I think all can be adequately dealt with by the authors simply expanding the text with information that should be easily accessible to them (as opposed to conducting further research). However, I think the first four are essential for a proper understanding of the manuscript.
Report #3 by Anonymous (Referee 3) on 20221020 (Invited Report)
 Cite as: Anonymous, Report on arXiv:scipost_202209_00024v1, delivered 20221020, doi: 10.21468/SciPost.Report.5935
Strengths
1. The manuscript addresses an important physical question : whether or not certain runandtumble models can show phase separation in one dimension.
2. The authors introduce a lattice model of restricted tumbling which is amenable to well established tools of analytical treatment.
3. By mapping an interacting runandtumble model to one of beadsinurn, the manuscript presents a useful technical path in studies of interacting one dimensional active systems.
Weaknesses
The paper, although wellwritten at parts, is still too technical in its approach. While it is a judicious choice by the authors to maintain the flow and not bombard the text with jargon, I also feel that some further explanation of the mathematics would really help the article.
Report
The authors map a system of interacting restricted runandtumble particles to a system of beads in urns. This mapping allows for an analysis of the original system through Matrix Product Ansatz, and shows that certain models of interacting runandtumble particles cannot produce a sustained phase separated state in the steady state.
I believe the work is interesting and provides an alternate route (from the coarsegrained studies involving offlattice velocity jump processes) to study interacting active matter physics in onedimensional geometries. To my knowledge, the work is original and has the potential to lead to followup work by the community working in theoretical aspects of constrained active matter systems.
However, I also have the following concerns with the text and I feel addressing these will only enhance the readability of the manuscript for a wider audience (beyond the community working in exact analysis of interacting systems).
1. The analysis in page 3 can benefit from some physical explanations., for example, a termbyterm explanation of Eq. (3).
1a. Can we extract any physical insight from the choice of auxiliary matrices?
2. Are the simulations for obtaining F_r, F_l and F performed in the steady state? Then are the two urns (or particles) chosen at random? Or do the authors choose specific initial states?
3. What are the functional forms of \eta_{+} and \rho_{+} in terms of \gamma?
4. Can the authors comment on some transient properties? For example, if we had an initial state of a given density with all particles being placed on consecutive lattice sites and the rest empty, i.e., we start with a macrocluster. Is there a (heuristic) timescale for this cluster to break down into a state of globally homogeneous density? I expect the transient states to depend on the hopping and tumbling rates, as well as on the density of particles. Ca this line of thought give any insight for the nonexistence of MIPS in the model?
5. The authors assume that the hop rate u(m_k, m_{k+1}) depends on the number of particles in the departure site k and the arrival site k + 1. Can they elaborate on the insight behind this assumption? Should the hop rate also not depend on the internal states of the arrival and departure urns? Naively, one might assume the current within a domain to depend on the polarities of the particles at the end of the domain.
6. Can this restricted tumbling model show or ruleout microphase separation with formation of stable microclusters in the steady state?
7. Indeed the manuscript could benefit from citing further recent literature on interacting onedimensional runandtumble systems.
Requested changes
Kindly address the points mentioned in the report and correct typos, e.g., 'sate' > state ,'dependn' > depends .
Author: Pradeep Kumar Mohanty on 20221124 [id 3069]
(in reply to Report 3 on 20221020)A reply to Referee 3 is attached as a pdf file.
Attachment:
Report #1 by Anonymous (Referee 1) on 20221018 (Invited Report)
 Cite as: Anonymous, Report on arXiv:scipost_202209_00024v1, delivered 20221018, doi: 10.21468/SciPost.Report.5922
Report
Dear Editor,
in this Letter, the authors present a mapping from a onedimensional runandtumble dynamics with restricted tumbles to an urn model. They analyze the latter in the limit corresponding to two particles and show that, were this limit to faithfully describe the whole density range, phase separation would be impossible. They carry out a number of numerical simulations to support their claim that this limit is indeed qualitatively faithful.
The main appeal of the paper put forward by the authors is that questions have been raised "about the stability of [MIPS] states in 1D even though fluctuating hydrodynamic equations have predicted them [6,7]". References 6 & 7 correspond to early papers by Cates and Tailleur. However, Reference 6 states:
"In practice, of course, the very existence of a thermodynamic mapping in this (strictly 1D) system ensures that the bulk phase separation predicted by meanfield theory is replaced by Poissondistributed alternating domains of mean size".
I.e., there is no phase separation in 1d. I thus think that, contrary to the authors' claim, phase separation in 1d has never been considered possible. Furthermore, the authors consider a model in which at most one particle per site is allowed. It has been shown by Soto and Golestanian that this only leads to finitesize clusters (in 1d and 2d; See ref 14 and [Sep\'ulveda Soto Phys. Rev. E 94, 022603 2016]). I thus think that the results presented by the authors are rather in line with the existing literature, contrary to what the article suggests.
That said, the mapping to the urn model is interesting, and reminiscent of the mapping between the TASEP and the ZRP (See Section 2 of [MR Evans and T Hanney 2005 J. Phys. A: Math. Gen. 38 R195). Unfortunately, its analysis is very hard to follow and the algebra presented on page 3 will be impossible to understand for readers who have not worked on matrix ansatz and related models.
For these two reasons, I do not think that the short format chosen by the authors is appropriate to convey their otherwise interesting results: the main message is not novel enough that a short format is required and the latter makes the article very hard to read. I think this paper would have much more impact if it was extended and all the results were detailed and explained.
Finally, I do not think that Scipost Physics is the right venue for this article. The journal is described as publishing articles which "provide details on groundbreaking results obtained in any (sub)specialization of the field". I do not think that the current submission fits this type of description: it is a particular casestudy of a result that (I think) has been considered established since the early days of MIPS. That said, the results and the mapping are interesting and an extended, heavily improved, version could be considered for a more specialized publication.
More specific comments:
 In the introduction, the authors write "Theoretical investigations of this phenomenon have thus far concentrated on continuum models [811]". I think that Ref 8 deals with lattice models and would be better placed below with ref 14. The reference [Whitelam, Klymko, Mandal, "Phase separation and large deviations of lattice active matter.", The Journal of chemical physics 148.15 (2018): 154902] should also probably be cited.
 The second paragraph of the second column of the introduction should be fixed to reflect the content of the cited articles.
 Speaking about grand canonical ensemble in the introduction is surprising at first since ensemble equivalence for active matter has not been demonstrated. I think this is coming too early and that it would be much clearer later on, once the factorized steadystate has been put forward.
 "Nonexistence of MIPS transition in restricted tumbling model would imply that the same can not occur in any other RTP model in 1D". "The same" refers here to "Nonexistence of MIPS". I think the authors wanted to state the opposite: "Nonexistence of MIPS transition in restricted tumbling model would imply that MIPS cannot occur in any other RTP model in 1D".
 The "restricted" tumbling hypothesis is a crucial difference between this model and other models previously studied. Its definition and its discussion should come much earlier in the article. The role of the asymmetry (right neighbor vs any neighbor) should also be discussed.
 When the authors discuss the mapping to the beadsinurn model, they should comment on the relationship with the mapping from TASEP to ZRP.
 The beginning of the second column of page 2, when the current J(m) is discussed, is very hard to understand. This part should be clarified and detailed.
 The twourn section is not very clear: it implies $N=2$ but then an arbitrary $N$ comes back in Eq (2) and the N=2 case is only introduced later on in page 3. This should be completely rewritten: if the authors first want to do the general case, they should do so and only specify later on to the N=2 case. If they want to consider N=2 all along, they should do so. Furthermore, page 3 is amazingly dense and hard to read. I would strongly suggest detailing this part so that this can be understood by activematter readers who are not familiar with the matrixproduct ansatz.
 The authors should discuss the implication of the larg M limit, on page 3. This seems to imply a dilute limit which is not where MIPS is expected. Later on, they show the exponential distributions to survive at larger densities and I suggest that the authors clearly explain if this is a justification for using their results beyond the lowdensity limit.
 Figure 2: The authors should explain how the color code is related to the values of $\omega$.
Author: Pradeep Kumar Mohanty on 20221124 [id 3070]
(in reply to Report 1 on 20221018)A reply to Referee 1 is attached as apdf file.
Author: Pradeep Kumar Mohanty on 20221124 [id 3067]
(in reply to Report 4 on 20221028)A reply is attached as a pdf file.
Attachment:
referee_4.pdf