SciPost Submission Page
Twospecies hardcore reversible cellular automaton: matrix ansatz for dynamics and nonequilibrium stationary state
by Marko Medenjak, Vladislav Popkov, Tomaž Prosen, Eric Ragoucy, Matthieu Vanicat
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Matthieu Vanicat 
Submission information  

Preprint Link:  https://arxiv.org/abs/1903.10590v1 (pdf) 
Date submitted:  20190328 01:00 
Submitted by:  Vanicat, Matthieu 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approach:  Theoretical 
Abstract
In this paper we study the statistical properties of a reversible cellular automaton in two outofequilibrium settings. In the first part we consider two instances of the initial value problem, corresponding to the inhomogeneous quench and the local quench. Our main result is an exact matrix product expression of the time evolution of the probability distribution, which we use to determine the time evolution of the density profiles analytically. In the second part we study the model on a finite lattice coupled with stochastic boundaries. Once again we derive an exact matrix product expression of the stationary distribution, as well as the particle current and density profiles in the stationary state. The exact expressions reveal the existence of different phases with either ballistic or diffusive transport depending on the boundary parameters.
Current status:
Reports on this Submission
Report #2 by Anonymous (Referee 2) on 201959 (Invited Report)
 Cite as: Anonymous, Report on arXiv:1903.10590v1, delivered 20190509, doi: 10.21468/SciPost.Report.942
Strengths
exact analytical results for interacting system
good interpretation with respect to general outofequilibrium theory
timely and subject of high current interest
Weaknesses
the system is very particular, universality of results / principles not clear yet
Report
In this paper, the authors study a new model of cellular automata, analysing various outofequilibrium situations. Analytic results are obtained both for the exact states emerging and for averages of various observables, which are accurate and rigorous. These are compares with explicit numerical evaluations, showing good agreement. The hydrodynamic limit of large space and time is taken, and interpreted within standard hydrodynamic theory. The model is a model for two species of interacting particles with hard core interaction.
The paper is precise, with exact analytical expressions for interesting physical quantities, and fits very well within the modern quests in nonequilibrium physics. A strength is that it combines analytical accuracy and physically interesting questions. I gives additional information that is of great use in deciphering nonequilibrium phenomena. I recommend publication, but I’d like to see one aspect clarified.
This is about the hydrodynamic interpretation. On page 8, the standard shock analysis is performed, solving the RankineHugoniot condition for the twocomponent fluid, and this is used to analyse shocks between separable states. The condition of equality of shock speed for the components (eq 22) is said to be on physical grounds. This is not so clear. In general, there is no problem in having many shocks at different velocities, and the Riemann problem for twocomponent fluids will usually display a shock and a rarefaction wave, or modifications of such two structures. I think what the authors are really doing here is establishing the conditions for a single shock to separate two currentcarrying state.
This is then used later, pages 1416, to interpret the result. But there are some questions. Wouldn’t it be possible to interpret all discontinuities in figure 3 using hydrodynamics? The sharp ones due to the light cone surely also have an interpretation; are they the normalmode velocities associated to the vacancy propagation? So these would be ``contact singularities”, and the associated modes would be essentially free (and of which there are two due to the oddeven separation), and freeness then explaining why they would not display diffusive spreading? While the middle one, with diffusive spreading, isn’t it also a contact discontinuity (i.e. a discontinuity parallel to the mode’s velocity  perhaps this is just a question of words, here referred to as marginally stable)? Overall, the fluid has three components. So don’t we indeed have in general three normal modes, and aren’t these the three velocities that we see emerging?
Finally, my understanding is that on both sides of the diffusive discontinuity, the authors expect the state to be essentially one of the currentcarrying states studied around eq 22. Can they see this analytically?
With all this, one would indeed have a standard picture for the Riemann problem of a 3component fluid, with 3 structures (here discontinuities), and in the present particular case, all being contact discontinuity (which is actually unusual, and point probably to the integrability of the model).
Requested changes
add discussion elements about hydrodynamic interpretation (see report), especially of first case studied (but also of other cases if possible), in order to possibly better extract general principles.
Report #1 by Anonymous (Referee 1) on 201951 (Invited Report)
 Cite as: Anonymous, Report on arXiv:1903.10590v1, delivered 20190501, doi: 10.21468/SciPost.Report.925
Strengths
1 Derives exact results for outofequilibrium transport
2 Simple approach based on matrix product ansatz, independent of the Bethe ansatz / generalized hydrodynamics technology
3 Rich physics, with features of both ballistic and diffusive transport
Weaknesses
1 The considered model has a deterministic dynamics, which might miss some aspects of interacting quantum integrable models
2 Writing and grammar should be improved
Report
The authors introduce a onedimensional twospecies model with deterministic dynamics, driven out of equilibrium from three different protocols : a local quench at the origin, an inhomogeneous quench corresponding to joining two halves of the system each in a different stationary state, and boundary driving.
The time evolution of the density profiles (or, respectively, their stationary distribution in the case of boundary driving) are determined analytically using a matrix product ansatz. From there the scaling limit (large time, long distances) is obtained by recasting the analytical expressions by contour integrals and evaluating them at their saddle point. This allows for the observation of nontrivial transport features, among which a diffusive dynamics and the propagation of ballistic fronts.
The issue of transport, whether ballistic, diffusive or subdiffusive in manybody systems has attracted a lot of attention over the last years. In this perspective, the present works is remarkable in that it derives exactly some interesting features of transport, without relying on the assumptions of generalized hydrodynamics. It therefore comes as an interesting alternative to the latter.
Some criticism can still be adressed. It may be argued that models such as the one considered here may not catch in full generality the richness of interacting quantum manybody system: the dynamics here is indeed purely deterministic, and it would be interesting if the authors could discuss a little what it might miss compared to nondeterministic models. For the latter, it seems reasonable to expect that no such simple matrix product solution could exist.
Last but not least, there is some room for improvement in the grammar and writing (see details in next paragraph).
Once these aspects and the one addressed below are taken in consideration, i recommend publication in JSTAT.
Requested changes
1 Many syntax or spelling errors, which need to be corrected. For instance, 2nd paragraph of first page, misspell of "analytically". Later, "in the context of ["the" missing] Totally Asymmetric Simple Exclusion Process". 3rd paragraph, "as a propagation of ["the" to be replaced by "a"] local observable". 3rd paragraph of section 2.1, "usefull" should be replaced by "useful", etc...
In many other places the writing seems a bit careless, and should be improved before publication.
2 In the introduction, since the TASEP is mentioned, it would be useful if the authors could recall in a few sentences the state of the art for the latter, as far as the computation of outofequilibrium density profiles is concerned.
3 The authors refer to a different approach to cellular automata, based on soliton counting. The deterministic dynamics considered here looks like it could allowed such methods for the present problems. Have these been considered by the authors? Do they fail where the matrix product approach works ? If so, this should probably be explained in the introduction.
4After equation (6) the author mention a connection with the YangBaxter equation, and "the fact that $R(0)=U$ implies the integrability of the model with periodic boundary conditions.". This Last sentence is in fact not completely clear to me : while integrable R matrices can indeed be used to build commuting sets of transfer matrices in different geometries (rowtorow, or diagonaltodiagonal), it is not clear whether the generator $\mathbb{U}$ for the dynamics considered here really belongs to this setup. It would correspond to a diagonaltodiagonal array of R matrices all carrying an infinite spectral parameter, which seems difficult to relate to the usual construction. Could the authors clarify whether the YangBaxter integrability plays any role here ?
5 in eqs (11) and (12), $p'$ has not been defined.
6 eq (13) is slightly confusing, as the labels $\alpha_{0,\pm}, \beta_{0,\pm}$ have already been used in eq. (10), but are a priori unrelated (even though in the particular case of boundary driving they turn out to be the same). Could this be fixed ?
7 The notation $\langle +0 \rangle_{12}$ (and other similar ones) should be explained.
8 before equuation (22) : "On [remove "the"] physical grounds"
9 section 2.2 : Could the authors briefly discuss what happens when the condition (23) is not met ? Does the shock smoothen and decay ?
10 section 3.3, after (62) : maybe a few words to explain what a "marginally stable" shock means in practice.
11 Section 4.1, after (64) : "This dependence is reminiscent [of] the defect..."
12. Figure 4 : is there a specific reason why the authors chose $\rho_+ = 0.01$, and not $0$ ? Would the plots look much different there ?