Signs of Kravchuk matrix asymptotically produce a large circular region with hyperbolic sinks. Why? The Kravchuk matrix of dimension $n+1$ is such that its entries satisfy $$K_{i,j}^{(n+1)}=[x^i](1+x)^{n-j}(1-x)^j\quad\forall0\le i,j\le n.$$ It enjoys properties such as involution and has various expressions involving combinatorial sums.
Figure. Signs of the transpose of $K_{i,j}^{(1000)}$: Positive entries filled orange and negative entries filled white.
 

We observe that as odd $n\to\infty$, there is a large circular region, inside of which lie numerous hyperbolic sinks of various sizes.

A few weeks ago I posted this observation on MSE, and recently got a detailed answer explaining the circular shape. Briefly, using the theory of $_2F_1$-hypergeometric functions, we find that the circular parametrisation $$(1-2\alpha)^2+(1-2\beta)^2=1$$ with $\alpha=i/n$ and $\beta=j/n$ causes equation (2.5) of Paris (2013) to be reduced to the boundary line $z_\ast^+=-1$.
However, the question of why the hyperbolic sinks appear, and their fractal-like appearances remain unanswered. It also appears that the largest sinks lie on the line $i=j$.
Is there an explanation for this behaviour?
P.S. I found only one publication with a plot of the circular shape (Figure 7). Note that the sinks do not appear there since $n$ is too small.
 A: In Section 2 of our 1996 paper "Local statistics for random domino tilings of the Aztec diamond" (https://arxiv.org/abs/math/0008243), Jim Propp, Noam Elkies, and I analyzed this behavior inside the circle, and the same techniques work outside. The short answer is that you can write the coefficient of $x^i$ in $(1+x)^{n-j} (1-x)^j$ as a contour integral and compute asymptotics using the saddle point method. There are two critical points, and a phase mismatch between them give the Moiré effect shown in the picture. What we analyze in Proposition 4 in the paper is not actually your $K_{i,j}$, but rather $K_{i,j}K_{j,i}$. However, $K_{i,j}$ and $K_{j,i}$ differ by some simple factor (involving factorials), so this is essentially the same as analyzing $K_{i,j}$ in isolation.
I don't know of a simple description of the Moiré effect. It amounts to the factor of $\cos^2 \Phi(\ell,m;n)$ in Proposition 4 (but for your case of just $K_{i,j}$, the cosine factor wouldn't be squared, so it would change sign). We don't write out what $\Phi$ is explicitly in our paper, but it's determined by the computations in the proof, and we give a few properties in lemmas later in the section. For our purposes, all we needed to know about this $\cos^2 \Phi$ factor was that it averaged to $1/2$ in certain sums, and that could be proved with some exponential sum estimates.
My memory of these formulas is that they were rather complicated and not so easy to interpret nicely. However, there may well be a conceptually clearer description, which would be really interesting. It's not something we explored carefully, since for us the Moiré effect was more of an obstacle than an object of study itself.
