SciPost Submission Page
Chaotic Correlation Functions with Complex Fermions
by Ritabrata Bhattacharya, Dileep P. Jatkar, Arnab Kundu
|As Contributors:||Ritabrata Bhattacharya · Dileep Jatkar|
|Date submitted:||2020-09-21 17:50|
|Submitted by:||Bhattacharya, Ritabrata|
|Submitted to:||SciPost Physics|
We study correlation functions in the complex fermion SYK model. We focus, specifically, on the h = 2 mode which explicitly breaks conformal invariance and exhibits the chaotic behaviour. We explicitly compute fermion six-point function and extract the corresponding six-point OTOC which exhibits an exponential growth. Following the program of Gross-Rosenhaus, we estimate the triple short time limit of the six point function. Unlike the conformal modes with high values of h, the h = 2 mode has contact interaction dominating over the planar in the large q limit.
Author comments upon resubmission
List of changes
We thank both the Referees for carefully going through our manuscript and raising extremely constructive comments and criticisms. We apologise for late resubmission, which was partly due to Covid19 related lockdown and poor internet connectivity.
In this revised version, we have attempted our best to address these issues. Below, we enlist the specific changes that have been incorporated in the revised version.
List of changes:
1. We agree with the Referee’s comment and we have removed the corresponding remark.
2. We agree completely with the Referee on this point. Our aim was to provide an impression of the structure. Therefore, we have added a footnote clarifying further the notational abuse that we have made use of, for simplicity. In the explicit calculation, we have taken care of this shift inside the calculations.
Also, we have further added comments regarding the Referee’s comment on “more” out-of-time order vs “less” out-of-time order.
3. We have attempted a more self-contained description in section 2 and 3. The updates include further explanation following eqns (2.5) and (2.6), (3.1), (3.2) and diagram of contact and planar contributions.
4. We have completely re-written the earlier subsection 4.4 into a new section 5. In this section, we have provided many explicit details of the calculations, as well as improved on the figures and related details. We hope this better explains the results and addresses the questions raised by the Referee.
1. We have added a short paragraph in the Conclusion, adding a comment in connection of the ref 1712.04963. We thank the Referee for pointing out this reference to us.
2. We have added a footnote demonstrating explicitly how the combinatoric factor of 90 arises in the eqn (3.3).
3. We have stated the goal clearly, at the end of eqn (3.3).
4. The paragraph following eqn (4.15) is edited to clarify the issues. We have also explicitly included the reference  here.
5. Added Appendix B with the relevant and explicit details of the various functions.
6. Clarified that the 6-point correlator decreases with increasing \mu\beta.
7. We have clarified first two points summarising results of the computation of six point correlation function. Hopefully this addresses the referee’s concern.
Submission & Refereeing History
You are currently on this page
Reports on this Submission
Anonymous Report 2 on 2020-10-20 Invited Report
The authors have accommodated the changes suggested in the previous reports. I find that the manuscript is a lot more streamlined than the previous version. I recommend it for publication after some rather minor corrections listed below.
1- Equation of curly-J on page 4 is incorrect by a square-root.
2- The organisation of the paper at the end of the section 1 needs an update because of the inclusion of an additional section 5.
3- there seems to be a typo above equation 2.7 “They are antisymmetric, respectively symmetric under …”
4- In figure 1, the vertex between \tau_1, \tau_5 and \tau_3 is not correctly depicted.
5- There are (6C2 x 4C2)/3! options. The division by 3! Is because interchanging the vertices among themselves doesn’t give new diagrams. Moreover, the more pressing issue is that there are index contractions between 1-2, 3-4 and 5-6 labeled fermions. So unless I am missing something, these additional diagrams are suppressed in N counting.
Anonymous Report 1 on 2020-10-11 Invited Report
I thank the authors for improving the manuscript and explaining their methods. Unfortunately I still cannot recommend the paper for publication without substantial improvements.
In sections 4.3 and 5 it is still unclear to me what is being done. It would help if the authors could define exactly what 6-point OTOC they are computing (in Lorentzian signature). After (5.1) it is said that the six insertions are ordered in Euclidean time (t) as 0<t1<t3<t5<t2<t4<t6<2pi. This Euclidean ordering determines the order in which the operators appear on the time folded contour. It seems the authors never address the Lorentzian time insertions. Whether the analytic continuation of the 6-point function with this specific Euclidean ordering gives a 6-point OTOC, or a TOC, or something in between that can be represented on fewer than 3 timefolds, depends on the Lorentzian times as well. If I guess that the authors have in mind Lorentzian times (T) of the form T1=T2>T3=T4>T5=T6 with pairwise identical operators, then the Euclidean ordering described after eq. (5.1) is actually not a 6-point OTOC. It can be represented on a contour with only two timefolds instead of three. I would therefore not expect this to display the features of a true 6-point OTOC.
Relatedly, as explained around (5.1), the authors seem to study the 6-point function only as a function of t1+t2, holding all other times as well as t1-t2 fixed. Ignoring the fact that these aren't even Lorentzian times, there remains another issue: defining the 6-point Lyapunov exponent as growth w.r.t. a single time difference seems unnatural (of course, depending on the configuration). But in any case I would think that a good definition of 6-point Lyapunov growth should know something about all six insertion times.
As a concrete suggestion, I would recommend computing a 6-point function with Euclidean times 0<t1<t3<t2<t5<t4<t6<2pi and Lorentzian times as above. This is a true 6-point OTOC (i.e., it cannot be drawn on a contour with less than 3 timefolds). Comparing it to the configuration described above would be interesting. This goes hand in hand with a comparison to reference  which the authors have added in response to the other referee: in the updated discussion section it is claimed that results found here are consistent with the "maximally braided" 6-point OTOC of . I don't think this is correct for two reasons:
(i) As explained above, the OTOC that is being computed here does not seem to be truly "6-point OTO", so it is quite different from the one studied in .
(ii) Reference  claims that the Lyapunov exponent in the 6-point function is unchanged compared to the 4-point function, which seems to contradict the main result of the present paper (and therefore should be clearly addressed by the authors).
Regarding the numerics: thank you for the additional explanations. I understand the method now. However, I am not convinced that it is correct. The authors compute a purely Euclidean 6-point function (which depends on 6 times), fit its dependence on a single time coordinate to a sinusoidal and expect to read off the Lyapunov exponent from that. Several questions:
(i) It is well known that numerical analytic continuation is quite subtle (essentially because it is not unique and hence ill-defined). For example, at the end of Appendix G of ref.  a method is described for performing such an analytic continuation in case of 2-point functions: this numerical algorithm is much more subtle than what the authors do here. If my concerns are unwarranted, I apologize, but do encourage some explanation as to why none of these issues matter here.
(ii) More importantly: for the reason I explained above, knowing only the Euclidean ordering of the operators in the 6-point function does not uniquely determine the real-time OTOC at all. It might not even be an OTOC, depending on what Lorentzian times one chooses. For this reason I am suspicious about obtaining the Lyapunov exponent without doing any real-time calculation at all, or even specifying the Lorentzian insertion times.
To summarize: I appreciate the more detailed explanations. However, after understanding better how the 6-point OTOC was computed, I am now worried about the correctness of this calculation. The main reason for my worry is that the authors never seem to do any calculation in Lorentzian time. As explained above, having only the Euclidean 6-point function is not enough to even specify what 6-point OTOC one is talking about. It might be that the authors have a very specific Lorentzian configuration in mind and their simple method of analytic continuation actually computes it correctly. But if that is the case, I am not able to see it from the current manuscript.