Fermionic Gaussian states: an introduction to numerical approaches

We really thanks the Referee 1 for their report and in-depth analysis of section 6.3. We implemented the suggested changes as follows:

• Added reference http://doi.org/10.1103/RevModPhys.80.517 and references therein as a review on the topic of entanglement in many-body systems.
• Added references https://doi.org/10.1088/1742-5468/2015/06/p06002 and https://doi.org/10.1088/1751-8121/aa6268 for the rainbow chain.

Sec.6.3:

• Added a section on the definition of MPS in the appendix A.2 and references https://doi.org/10.1016/j.aop.2014.06.013, https://doi.org/10.1103/RevModPhys.93.045003 and https://doi.org/10.1016/j.aop.2010.09.012.
• Fig. 10 and 11 are now referenced in the text (p.44). The captions have been corrected.
• The text has been updated to make it more clear, in particular we corrected to ""it makes $\vec{v}$ also an approximate eigenvalue of Λ" with the sentence "This makes $\vec{v}$ also an approximate eingenvector of Λ. In fact we know that the eigenvalues of Λ are 0 or 1, thus it is sufficient to extend the dimension of the eigenvector $\vec{v}$ adding to it N − l zeros to obtain an eigenvector of Λ." Changed "Since $\Lambda$ is a pure state" with "Since $\Lambda$ corresponds to a pure state". Changed "let $\Lambda$ be the ground state" with "Let $\Lambda$ be the correlation matrix of the ground state of a 1D local Hamiltonian."
• The notation has been uptated to make it consistent.
• The figures have been updated to make it clear that the label $m+1$ referers to the dimension of the block.
• The notation has been updated to make it consistent.
• We created new plots for bigger sizes and for a bigger value of $m$ as suggested by the referee. We updated the caption making clear that in such a situation one has an extensive scaling of the entropy with subsystem size.

Typos:

The typos have been corrected.

We gratefully thanks the Referee 2 for their careful and in-depth report. We implemented the changes suggested as follows:

1) We agree with the referee. We added a footnote with the definition of Hamming distance. 2) We agree with the referee. We moved the discussion where suggested. 3) We agree with the referee. We moved the box where suggested. 4) We thought about the possibility of defining the Pfaffian in the text. We believe that the Pfaffian is not so relevant to require its cumberson definition in the main text or in the appendix. Adding the definition of Pfaffian as a footnote is not going to help much either, as in order to define the Pfaffian one also have to introduce other elements such as the Levi-Civita symbol of groups of permutations. In the explicit examples following the introduction of the Pfaffian, some simple expressions of the Pfaffian are explicitly presented. 5) We agree with the referee. We wrote down the case for $N=2$ explicitly. 6) The variable $i$ was commented at the end of the box. Nevertheless, it was not discussed explicitly as an argument of the function. For this reason we added a comment about $i$ in the first sentence of the box as requested by the referee. 7) Added the minus sign. 8) We agree with the referee, we improved the position of the boxes. 9) Before the coding frame we wrote "in the following program we numerically compute the time evolution induced by a hopping Hamiltonian on a random translational invariant gaussian state with exponentially decaying correlation functions." The $\Gamma$ represents this state. 10) We agree with the referee. We expanded the caption of the table as requested. 11) We agree with the referee. We added more comments to the code. 12) Corrected the label of the figure. 13) We removed the word "experimental". 14) Added "Possible escape tricks exists. Perhaps one can slightly perturb the state and then compute the value of the energy after another step of the evolution, this can help in escaping from local minima. In general such tricks have to be adapted to the particular needs. " as requested. 15) Figure 11 is now called in the main text. 16) We corrected the figure numbering. 17) We agree with the referee. We merged figure 12 and 13. 18) There was a general problem with the tags of the figures. We corrected it. 19) We agree with the referee. We moved the box, as well as the code frames and the associated figures in the appropiate section. 20) We corrected the typos.

We really thanks the Referee 1 for their report and in-depth analysis of section 6.3. We implemented the changes suggested as follows:

+ Added reference http://doi.org/10.1103/RevModPhys.80.517 and references therein as a review on the topic of entanglement in many-body systems.
+ Added references https://doi.org/10.1088/1742-5468/2015/06/p06002 and https://doi.org/10.1088/1751-8121/aa6268 for the rainbow chain.

Sec.6.3:

+ Added a section on the definition of MPS in the appendix A.2 and references https://doi.org/10.1016/j.aop.2014.06.013, https://doi.org/10.1103/RevModPhys.93.045003 and https://doi.org/10.1016/j.aop.2010.09.012.
+ Fig. 10 and 11 are now referenced in the text (p.44). The captions have been corrected.
+ The text has been updated to make it more clear, in particular we corrected to ""it makes $\vec{v}$ also an approximate eigenvalue of Λ" with the sentence *"This makes $\vec{v}$ also an approximate eingenvector of Λ. In fact we know that the eigenvalues of Λ are 0 or 1, thus it is sufficient to extend the dimension of the eigenvector $\vec{v}$ adding to it N − l zeros to obtain an eigenvector of Λ."* Changed "Since $\Lambda$ is a pure state" with "Since $\Lambda$ corresponds to a pure state". Changed "let $\Lambda$ be the ground state" with "Let $\Lambda$ be the correlation matrix of the ground state of a 1D local Hamiltonian."
+ The notation has been uptated to make it consistent.
+ The figures have been updated to make it clear that the label $m+1$ referers to the dimension of the block.
+ The notation has been updated to make it consistent.
+ We created new plots for bigger sizes and for a bigger value of $m$ as suggested by the referee. We updated the caption making clear that in such a situation one has an extensive scaling of
the entropy with subsystem size.

Typos:

The typos have been corrected.

We really thanks the Referee 2 for their careful and in-depth report. We implemented the changes suggested as follows:

1) We agree with the referee. We added a footnote with the definition of Hamming distance.
2) We agree with the referee. We moved the discussion where suggested.
3) We agree with the referee. We moved the box where suggested.
4) We thought about the possibility of defining the Pfaffian in the text. We believe that the Pfaffian is not so relevant to require its cumberson definition in the main text or in the appendix. Adding the definition of Pfaffian as a footnote is not going to help much either, as in order to define the Pfaffian one also have to introduce other elements such as the Levi-Civita symbol of groups of permutations. In the explicit examples following the introduction of the Pfaffian, some simple expressions of the Pfaffian are explicitly presented.
5) We agree with the referee. We wrote down the case for $N=2$ explicitly.
6) The variable $i$ was commented at the end of the box. Nevertheless, it was not discussed explicitly as an argument of the function. For this reason we added a comment about $i$ in the first sentence of the box as requested by the referee.
7) Added the minus sign.
8) We agree with the referee, we improved the position of the boxes.
9) Before the coding frame we wrote "in the following program we numerically compute the time evolution induced by a hopping Hamiltonian on a random translational invariant gaussian state with exponentially decaying correlation functions." The $\Gamma$ represents this state.
10) We agree with the referee. We expanded the caption of the table as requested.
11) We agree with the referee. We added more comments to the code.
12) Corrected the label of the figure.
13) We removed the word "experimental".
14) Added "Possible escape tricks exists. Perhaps one can slightly perturb the state and then compute the value of the energy after another step of the evolution, this can help in escaping from local minima. In general such tricks have to be adapted to the particular needs. " as requested.
15) Figure 11 is now called in the main text.
16) We corrected the figure numbering.
17) We agree with the referee. We merged figure 12 and 13.
18) There was a general problem with the tags of the figures. We corrected it.
19) We agree with the referee. We moved the box, as well as the code frames and the associated figures in the appropiate section.
20) We corrected the typos.