Out-of-equilibrium Eigenstate Thermalization Hypothesis

The paper is certainly interesting and perhaps even overdue as since early days of ETH there were many discussions of relation between ETH and relaxation of observables. It was also realized that the key are the correlations between the matrix elements of observables and the overlaps of the initial state with the eigenstates but to my knowledge this is the first serious work in this direction.

I think the paper is well written and well organized. In addition it contains various numerical and analytical tests of consistency of results so there is no question in my mind that it fully deserves publication.

Below I want to list some questions and suggestions, which are optional, and the authors might (or might not) want to address them.

1. The authors effectively study cross-correlations between the matrix elements of the i) projector operator to the initial state $P=|\Psi\rangle\langle \Psi|$, which, as they note, is the rank one operator and has special properties, and ii) the matrix elements of some observable $A$. I think it would be nice first to formulate similar cross-correlation function between two observables say $A$ and $B$ and define the analogue of the function $f(\omega) g^\ast (\omega)$. When $A=B$ this would be $|f(\omega)|^2$, which is studied a lot in the literature, but generally $f\neq g$. After formulating this cross-correlation for arbitrary operators it would be easier to understand the subtleties peculiar to $ B=P$.

2. From Eq. (5) it follows that $|\Psi_{ij}|^2=\Psi_{ii}\Psi_{jj}$ , which seems to suggest that there are roughly $D^2$ independent phases in $\Psi_{ij}$ but only $D$ independent amplitudes. I wonder if this constraint plays any role in defining $\tilde R_{ij}$. Perhaps this is what allows the authors to derive Eq. (13), but maybe the authors could say in words what is important. I clearly can extend Eq. (13) to add products with more terms and more constraints on $\tilde R$. Do they follow automatically from (13) or are there more constraints?
3. Related to the previous question. Can I use this machinery to study objects like $\langle \psi | A(t) A(t') |\psi\rangle$ or one must introduce more functions like in higher odder ETH the authors developed. Maybe a general comment of how one can use the formalism for practical calculations. By the way, I am not sure that "out of equilibrium ETH" is a good title. ETH is still equilibrium to my taste as the authors use equilibrium eigenstates, just as I wrote, the important point is that one of the operators is the projector operator. If we e.g. use a projector to a mixed state, for example a pure polarized state for 10 sites in the middle and the identity matrix or some thermal state for the rest. For which length of the partial projector the suggested ETH will become equilibrium then?
4, I wonder if the authors can comment on how their formalism connects with the Kubo linear response when the state $|\Psi\rangle$ is a perturbed eigenstate of the Hamiltonian with say the operator $\epsilon A$ or more generally $\epsilon B$. In this case $\langle \Psi| A(t) |\Psi\rangle$ should reduce to a standard integral of a two-point function and hence be expressed through the standard ETH. This question connects to my question 1.

5. In Fig. 2 there seem to be a good convergence to the thermodynamic limit at large frequencies as perhaps expected. Can the authors make some claims about say short time relaxation. Does this convergence differ qualitatively from convergence of the cross-correlation (or the spectral function) of normal observables? Does the new ETH formalism allow them to reduce finite size effects?

6. The easiest one :). There is a typo right below Eq. (28).

Publish (surpasses expectations and criteria for this Journal; among top 10%)