Consider a covariance function (positive definite function) on $\mathbb{Z}$:
$$
\gamma(k)=(1+|k|)^{-\alpha},\quad \alpha>0.
$$
It is guaranteed to be positive definite by Polya's criterion (monotonicity and convexity). Consider matrices
$$
\Sigma_1=[\gamma(i-j)]_{1\le i,j\le m},\quad \Sigma_2=[\gamma(k+i-j)]_{1\le i,j\le m},\quad m>0,~k\ge 0.
$$
**Show that all entries of $\Sigma_1^{-1}\Sigma_2$ are nonnegative.**

*Remark: I checked it numerically for different $k$, $m$ and $\alpha$. A weaker statement involving large enough $k\asymp m^{1+\epsilon}$, where $\epsilon>0$ is arbitrarily small, is also of great interest. A further question is to identify general conditions on a positive definite $\gamma(k)$ which fulfills the preceding requirement.*

*Remark:* The problem can be reinterpreted in terms of linear prediction of stationary time series $\{X(n)\}$ whose covariance function is $\gamma(k)$.
It is the same as to show that the linear predictor $\hat{X}(k)=\sum_{j=1}^m c_jX(j)$ which minimizes the mean square error $E |X(k)- \hat{X}(k)|^2$ satisfies $c_j\ge 0$.