Is this Markov chain of Gaussian matrix products $G_1 G_2 \dots G_m$ ergodic? Consider Markov chain $\{X_t\}_{t\in N}\subseteq R^{n\times n}$ defined by $X_{t} = X_0 G_1 \dots G_t$ where $G_i$'s are iid Gaussian matrices $G_1,\dots,G_t\sim N(0,1/n)^{n\times n}$, and $X_0$ is some deterministic matrix with full rank fixed scale, $\|X_0\|_F^2=n$ and $\operatorname{rank}(X_0)=n$. Is this chain ergodic? Put in other words, assuming that initial matrices $X_0$ and $Y_0$ share the some "nice" properties (eg. similar scale and full rank), can we find a coupling(*) of two copies of this chain $\{X_t\}_{t\in N}$ and $\{Y_t\}_{t\in N}$ such that $P(X_t\neq Y_t) \le \alpha^t$ for some $\alpha<1$?
Intuitively this coupling should be plausible: Since Gaussian matrices $G_i$'s are rotation invariant, we can couple the $X_t$ and $Y_t$ by optimising for rotation. But I come short of making these intuitions any more rigorous. I've looked up literature on random matrices but can't find anything relevant.
(*) More elaborately, can we design $\{(G_t,W_t)\}_{t\in N}$ where both $G_1,\dots, $ and $W_1,\dots$ are independent, both jointly are coupled to lead to mixing of the two chains?
 A: You are dealing here with the products of invertible matrices, and the resulting Markov chain is known as a random walk on the corresponding group $GL(n,\mathbb R)$. A qualitative asymptotic "boundary" theory of such products was created by Furstenberg in the early 60's (see his 1963 papers "A Poisson formula for semi-simple Lie groups" and "Noncommuting random products"). He mostly talks about $SL(n,\mathbb R)$ and more general semi-simple Lie groups, but it doesn't make much difference.
I will return to Furstenberg's theory in a moment, but let me first say a couple of words about the "coupling" conditions you mention, in the "rawest" form. When talking about a general Markov chain one should in principle distinguish its mixing and ergodicity. Mixing is equivalent to the triviality of the tail $\sigma$-algebra of the chain or to the asymptotic independence of the time $n$ distributions of initial states:
$$
\| (\delta_x - \delta_y) P^n \| \to 0 \qquad \forall x,y \;,
$$
where $\|\cdot\|$ is the total variation, and $P$ is the transition operator of the Markov chain. In a similar way, ergodicity is equivalent to the triviality of the shift invariant $\sigma$-algebra in the path space (i.e., to the absence of non-constant bounded harmonic functions), or to the independence of the Cesaro averages of time $n$ distributions of initial states:
$$
\frac1n \left\| \sum_{k=1}^n (\delta_x - \delta_y) P^k) \right\| \to 0 \qquad \forall x,y \;.
$$
For the random walks that you consider mixing and ergodicity are actually equivalent. It follows from so-called 0-2 laws that also provide further information about the above differences.
Now, one of the consequences of Furstenbrg's theory is that random walks on matrix groups can only be ergodic ($\equiv$ mixing) in very special cases. Gaussian random walks on the linear group are not among them. In particular, in your setup
$$
\lim_n \| (\delta_x - \delta_y) P^n \| \neq 0
$$
for any two initial matrices $x,y\in GL(n,\mathbb R)$ unless $x$ and $y$ are proportional.
