SciPost Submission Page
On the generality of symmetry breaking and dissipative freezing in quantum trajectories
by Joseph Tindall, Dieter Jaksch, Carlos Sánchez Muñoz
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Dieter Jaksch · Joseph Tindall 
Submission information  

Preprint Link:  https://arxiv.org/abs/2204.06585v2 (pdf) 
Date submitted:  20220622 11:38 
Submitted by:  Tindall, Joseph 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational 
Abstract
Recently, several studies involving open quantum systems which possess a strong symmetry have observed that every individual trajectory in the Monte Carlo unravelling of the master equation will dynamically select a specific symmetry sector to freeze into in the longtime limit. This phenomenon has been termed dissipative freezing, and in this paper we argue, by presenting several simple mathematical perspectives on the problem, that it is a general consequence of the presence of a strong symmetry in an open system with only a few exceptions. Using a number of example systems we illustrate these arguments, uncovering an explicit relationship between the spectral properties of the Liouvillian in offdiagonal symmetry sectors and the time it takes for freezing to occur. In the limiting case that eigenmodes with purely imaginary eigenvalues are manifest in these sectors, freezing fails to occur. Such modes indicate the preservation of information and coherences between symmetry sectors of the system and are associated with phenomena such as nonstationarity and synchronisation. The absence of freezing at the level of a single quantum trajectory provides a simple, computationally efficient way of identifying these traceless modes.
Author comments upon resubmission
Dear Editor and Referees,
We thank both referees for their thorough comments and reading of the manuscript. Following these comments, we have revised the manuscript substantially and believe it merits publication in Scipost Physics. Below, we first respond to the main criticisms of Referee 2, we then provided itemised responses to the individual comments of the referees. A list of the changes we have made to the manuscript is also provided for the resubmission.
The authors, Joseph Tindall, Dieter Jaksch and CarlosSanchez Munoz
Main Criticism from Referee 2 (Lack of precision and originality):
Here we address the referee’s criticism about the originality and precision of the manuscript.
Originality: We disagree that References [2426] communicated the picture we provide in Section 1 and 2 of our work. All these references considered a single microscopic model and numerically observed dissipative freezing (DF) (with Ref. 24 the only one to provide a mathematical proof and this was specific to the exceptional situation H = L). We instead make no assumptions about the details of the Liouvillians—other than its possession of a strong symmetry—and provide significant mathematical and physical insight into why freezing will emerge. We also identify a direct connection between rigorous mathematical results on ergodicity [4042] in open systems and dissipative freezing – something which Refs. [2426] all missed. We did not intend our work to constitute a rigorous mathematical proof of DF. The aim of this paper is, instead, to provide an insightful, comprehensive, mathematical and physical understanding of its origin, complemented with a range of concrete numerical examples. Alongside the intuition and comprehension it provides, such a perspective has allowed us to develop a number of original results without needing to rely on heavy mathematical artillery: exceptions to DF, a lower bound on the timescale on which DF occurs and a computationally efficient method for identifying traceless nondecaying modes in a Liouvillian. In our revised version of the manuscript we have made these points much clearer.
Precision: Both referees comments have significantly helped us improve the precision of the manuscript and the changes made are listed at the end of the document. These changes include mathematical definitions for certain quantities, rectifying grammatical and mathematical errors, avoiding discursive discussions, providing several elaborations when requested and deriving a rigorous lower bound on the timescale of dissipative freezing – which we demonstrate in our numerics. We believe that these changes have significantly improved the clarity and precision of the manuscript.
Report 1 Responses:
1) The timederivatives were in place to emphasize that the expression is `steady’ i.e. that the ensemble average of coefficients is 0 and unchanging in the longtime limit. The referee is correct, however, that this is superfluous and $m=0$ is sufficient because if the limit $t \rightarrow \infty$ converges asymptotically for $m=0$ then its derivative must vanish.
2) For a given trajectory (i.e. specific $i$) if freezing has occurred then there is only a single α for which the coefficients $c_{\alpha, \beta}^{(i)}(t)$ are not zero. Hence, for a given $i$, the product of coefficients for two different values of α is zero. We have now made this clear.
3) Yes, we did not intend our definition to imply you cannot have traceless nondecaying modes within a single symmetry subspace. We have renamed these as `intersector traceless modes’ in order to avoid suggesting these are the only possible traceless, nondecaying modes that can appear in a Liouvillian’s spectrum.
4) One is always able to introduce a new symmetry subspace to the system by expanding their Hamiltonian and jump operators via a direct sum. Physically this would correspond to considering an additional, independent open system on top of the existing one.
5) We now provide a mathematical definition of this average growth rate.
6) This time average is defined via an average between times $0$ and $\infty$. We have specified this and improved the readability here as the sentence was confusing in its original form.
7) The connection we were implying is that establishing an ergodic theorem with a single trajectory across multiple symmetry sectors is not possible due to the problem being over constrained in analogy with Eq. (10). In the interests of avoiding discursive discussion we have removed the reference to Eq. (10) here.
8) We have introduced explicit formulas for the quantities mentioned at the end of Section 2.3. With the addition of more explicit formulas for the freezetime and intersector spectral gap (see point 10), this should be easier to follow.
9) Yes, the first is a breakdown of freezing despite Eq. (8) holding whilst the second is a scenario in which Eq. (8) is no longer true. Section 2.5 has been reworded to make this clearer.
10 + 11) We have introduced an equation for the freeze time and intersector spectral gap. Our expression for the freezetime now follows from a new section where we derive a lowerbound on the timescale of freezing.
Report 2 Responses (under `Lack of Precision’)
1) The basis of steady states is at least $D$ dimensional, this is now stated explicitly.
2) We have provided additional discussion and references about the relationship between a stochastic unravelling of the GSKL equation and a continuous measurement process in an open system. We note that nonHermitian jump operators can result in a physical unravelling – for instance via photon counting in the environment.
3) The ‘Traceless’ mode exception completely covers the situation in which the intersector spectral gap vanishes, leading to Eq. (10) not being valid and freezing not occurring. The `Similar Symmetry Subspaces’ exception we detail is the only instance we are aware of where Eq. (10) can hold and the system does not undergo dissipative freezing. It is possible other examples where the structure of the two subspaces is sufficiently similar to prevent freezing exist. We have made this explicit, along with the fact we are, however, unaware of any such examples.
4) We thank the referee for spotting this error—we meant to say the projected Hamiltonian has the same form when focussing on the diagonal block where it is nonzero and have now made this mathematically precise, including the part about the unitary transform (change of basis).
List of changes
1) Introduction. We have mentioned the bound we derive on the freezetime and improved the wording to make the contributions of the paper to the field clearer.
2) Section 2.1. We have clarified that we mean the basis of steady states is at least D dimensional.
3) Section 2.2: We have expanded on what we mean by ‘generically’ in the definition of dissipative freezing.
Section 2.3: Recovering a Block Diagonal Steady State: The time derivatives have been removed from Eq. (10). The sentence following Eq. (10) has been made clearer. The discussion at the end of this subsection has been made more concise and clear.
4) Section 2.3: Long Products of Matrices: We have mathematically defined the average growth rate of the singular values of the propagation matrices in a given symmetry subspaces.
5) Section 2.3 Ergodicity and Quantum Systems Subject to Repeated Measurements: We have removed the discussion about Eq. (10) in this section to avoid being too discursive. We have mentioned that the timeaverage in the references is defined as between times $0$ and $\infty$.
6) New Section 2.4: The Timescales of Dissipative Freezing. This is a new section in which we derive a lower bound on the freezing time of trajectories. In this section we have also provided a mathematical definition for the intersector spectral gap, which plays a central role in the lower bound.
7) Section 2.5 (Previously Section 2.4): We have emphasized that Eqs. (8) and (10) are maintained in the first case whilst breaking down in the second, making the two exceptions very distinct. We have correlated the breakdown of freezing in terms of the bound derived in Sec 2.4. We have corrected a mistake in Eq. (15), including the unitary transform, and reworded the sentences following to provide clarity. We have also made it clear that other exceptions involving similar subspaces may exist beyond Eq. (15).
8) Section 3.1: We have provided a new and clearer definition (the qualitative results remain the same as before) of the average freezetime, connecting it to our new section 2.4 (on the timescales of freezing). We have mathematically defined the types of initial state used in the numerics to be precise.
9) Section 3.23.4 Numerical Results. We have redefined the freezetime and updated the plots based on this (the qualitative physics is unchanged). We have plotted the lower bound on freezetime in Figs 3c and 4c. We have significantly updated the discussions to reflect the earlier changes in the manuscript (e.g. new section and the earlier definition of quantities like the intersector spectral gap). This makes the discussion in this section much more precise.
10) Throughout: Grammatical errors pointed out by the referees have been fixed, included a naming error in the references.
Current status:
Reports on this Submission
Report 1 by Marco Cattaneo on 2022711 (Invited Report)
 Cite as: Marco Cattaneo, Report on arXiv:2204.06585v2, delivered 20220711, doi: 10.21468/SciPost.Report.5377
Report
The authors have modified the text according to my suggestions and I believe the presentation has considerably improved. In particular, the new section 2.4 provides us with a quantitative description of dissipative freezing that conveys the authors' message in a more clear way.
I would like to recommend the manuscript for publication, once the following minor typos have been fixed:
1) In eq.(10) the authors forgot to remove the text $m \in $ ...
2) In eqs. (21), (22) and the text below them there is a missing $'$ in $\beta$ in the definition of $\rho_R$.
Author: Joseph Tindall on 20220804 [id 2712]
(in reply to Report 1 by Marco Cattaneo on 20220711)We are pleased the referee recommends our work for publication and we have now fixed the typos they have pointed out.