Multiple Quantum Mpemba Effect: Exceptional Points and Oscillations: Appendix D: QMPE in Region (e)

cover
5 Mar 2024

This paper is available on arxiv under CC 4.0 license.

Authors:

(1) Amit Kumar Chatterjee, Yukawa Institute for Theoretical Physics, Kyoto University & Department of Physics, Ramakrishna Mission Vidyamandira;

(2) Satoshi Takada, Department of Mechanical Systems Engineering and Institute of Engineering;

(3) Hisao Hayakawa, Yukawa Institute for Theoretical Physics, Kyoto University.

Appendix D: QMPE in region (e): 3rd order EP

The region (e) in the (d,˜ Γ) plane of Hatano’s model is composed of a single point

FIG. 15. The figure (a) illustrates the two intersection time solutions t˜± [Eq. (D8)] in region (e) as functions of ˜dI and ˜dII.For comparatively high values of these control parameters, both solutions exist, leading to double QMPE. Parameters used are d˜ = 2√ 2, Γ =˜ Γ˜I = Γ˜II = 6√ 3. The figure (b) shows a surface from the three dimensional plot in figure (a). It demonstrates the intersection times [Eq. (D8)] by fixing ˜dI and varying ˜dII. We observe single QMPE as well as double QMPE depending on the parameter choices. Depending on the sign of ( ˜dI − ˜dII), the two solutions interchange their branches. Parameters used are ˜d = 2√ 2, Γ =˜ Γ˜I = Γ˜II = 6√ 3, ˜dI = 8.0.

It is a third order EP where all the non-zero eigenvalues become equal and their eigenvectors coalesce. Since this is a single point with fixed values of ˜d and Γ, the right eigenvectors are explicitly given by

Of course |r1i and hℓ1| have same expressions both in Eqs. (3) and (D2). The left eigenvectors are

) The non-zero eigenvalue of the Lindbladian [Eq. (2)] in region (e) is −iλ = −i4√ 3. Importantly, as a consequence of the 3rd order EP, we can only achieve an approximately diagonal form i.e. Jordan normal form LbJn leading to the equation

Due to the off-diagonal nonzero defects in Eq. (D5), the time evolution of density matrix elements and other observables pick up additional algebraic terms in time.

FIG. 16. The figure (a) shows an example of double QMPE in the ground state density matrix element in region (e). The main figure shows the first intersection between two copies I and II whereas the inset shows their second intersection that is weaker in magnitude. The figure (b) shows QMPE in average energy in region (e). The main figure shows that the energy difference between two copies I and II intersect zero two times, resulting in double QMPE. The insets show the intersections in the time evolution of the individual energies of the two copies. Parameters used are d˜= 2√ 2, Γ =˜ Γ˜I = Γ˜II = 6√ 3, d˜I = 5.0, d˜II = 7.0.

1. Density matrix elements

In region (e), the density matrix elements have the expressions below

In Fig. 15(a), we present existence of t˜+ and t˜− as function of the control parameters ˜dI and ˜dII using the fixed dissipation protocol. The solutions’ surfaces in Fig. 15(a) reveal that there are regions where only one of the solutions t˜+ or t˜− exist leading to single QMPE, and in other regions both the solutions exist resulting in double QMPE. Particularly, lower values of ˜dI and ˜dII correspond to single intersection time solution and therefore single QMPE. To understand the behaviors of t˜+,− more clearly, we fix ˜dI and display t˜+,− as functions of ˜dII in Fig. 15(b). The

FIG. 17. The figure (a) illustrates QMPE in von Neumann entropy in region (e). The difference of entropies between two copies I and II intersect zero twice, leading to double QMPE. The inset shows the second intersection that is weaker in magnitude with respect to the first intersection shown in the main figure. Parameters used in figure (a) are same as in Fig. 16. The figure (b) illustrates an example of thermal QMPE where the intersection between initially hotter and colder copies occur after the initial stage of negative temperature. The inset of figure (b) shows the full time evolution including the initial negative temperatures. Parameters used are ˜d = 2√ 2, Γ = 6 ˜ √ 3, Γ˜I = 14.0, Γ˜II = 15.0, ˜dI = 6.0, ˜dII = 3.0.

Fig. 15(b) shows that the two branches t˜+ and t˜− interchange as ˜dI < ˜dII. More specifically, for ˜dI > ˜dII we have t˜+ > t˜− and for ˜dI < ˜dII we have t˜+ < t˜−. The underlying reason is that the discriminant p γ 2 1 − 4γ2γ0 in Eq. (D8) is proportional to sign of ( ˜dI − ˜dII). Thus the two intersection times interchange their roles by changing the sign accompanying the discriminant as we go from ˜dI > ˜dII to ˜dI < ˜dII. The Fig. 16(a) demonstrates an explicit example of double QMPE in ground state probability in region (e). Like other regions, here also the second intersection (shown in the inset of Fig. 16(a)) is much feeble in comparison to the first intersection.

2. Average energy

For the same parameters used to obtain QMPE for ρgg(t) in Fig. 16(a), we would like to examine the possibility of QMPE in average energy E(t). Indeed, we observe in Fig. 16(b) that energy also exhibits double QMPE, marked by the twice intersection of zero by the energy difference ∆E(t) between the two copies. The insets of Fig. 16(b) shows separately the two intersections in the time evolution of the energies for the two copies. The second intersection is much weaker in magnitude compared to the first intersection just like the case of ground state probability. Note that the first intersection between the energies of I and II occur after they overshoot their steady state value which they reach during some early time of the relaxation process.

3. von Neumann entropy

Here we investigate the von Neumann entropy SvN(t), starting from two copies I and II. We use same parameters that have produced double QMPE both in ground state probability [Fig. 16(a)] and average energy [Fig. 16(b)]. In fact, the Fig. 17(a) shows that the entropy also exhibits double QMPE indicated by twice zero intersections (the second one shown in inset of Fig. 17(a)) of the entropy difference ∆SvN(t) between the two copies.

4. Temperature

In this section we study the thermal QMPE in region (e). We observe thermal QMPE in region (e) by using variable dissipation protocol in Fig. 17(b). There we observe single intersection and reversal of identities (i.e. hotter becoming colder and vice versa) occur after the negative temperatures have happened. So, such temperatures can be considered as well-defined and therefore this intersection in Fig. 17(b) corresponds to single QMPE in temperature in region (e).