[guided]The key point is that the commutator product is designed so that all first-order terms cancel and the first nonzero term is the matrix commutator. We work in a fixed submultiplicative matrix norm $\|\cdot\|_*$ on $M(n,\mathbb C)$ so that convergence and product estimates are controlled in one finite-dimensional norm.
Define
\begin{align*}
A:=[X,Y]=XY-YX.
\end{align*}
For $t$ near $0$, the matrix exponential series gives
\begin{align*}
\exp(tX)=I+tX+\frac{t^2}{2}X^2+t^3R_X(t),
\end{align*}
where $R_X(t)$ is bounded near $0$. This boundedness follows from absolute convergence of the exponential [power series](/page/Power%20Series) in the chosen submultiplicative norm: the tail beginning at degree $3$, divided by $t^3$, is locally bounded. The same statement holds with $Y$, $-X$, and $-Y$ in place of $X$.
Multiplying the four second-order expansions gives
\begin{align*}
\exp(tX)\exp(tY)\exp(-tX)\exp(-tY)=I+t^2(XY-YX)+t^3R(t),
\end{align*}
where $R(t)$ is bounded near $0$. Let us indicate the cancellation. The degree-$1$ part is
\begin{align*}
tX+tY-tX-tY=0.
\end{align*}
The degree-$2$ part consists of the intrinsic second-order terms and the pairwise products of first-order terms. After collecting them, the terms $X^2/2$, $Y^2/2$, $X^2/2$, and $Y^2/2$ cancel against the corresponding products coming from inverse factors, and the surviving mixed contribution is exactly
\begin{align*}
XY-YX=[X,Y].
\end{align*}
All remaining terms have degree at least $3$ in $t$, and boundedness of the exponential tails packages them as $t^3R(t)$ with $R$ bounded near $0$.
Now fix $u\geq 0$ and put
\begin{align*}
t_m=\sqrt{\frac{u}{m}}.
\end{align*}
For all sufficiently large $m$, $t_m$ lies in the neighbourhood where $R$ is bounded. If $K>0$ satisfies $\|R(t)\|_*\leq K$ there, then
\begin{align*}
C_m=I+t_m^2A+t_m^3R(t_m)=I+\frac{u}{m}A+E_m,
\end{align*}
where
\begin{align*}
E_m:=t_m^3R(t_m)
\end{align*}
and
\begin{align*}
\|E_m\|_*\leq K\left(\frac{u}{m}\right)^{3/2}.
\end{align*}
Multiplying by $m$ gives
\begin{align*}
m(C_m-I)=uA+mE_m.
\end{align*}
Since
\begin{align*}
\|mE_m\|_*\leq Ku^{3/2}m^{-1/2}\to 0,
\end{align*}
we obtain
\begin{align*}
m(C_m-I)\to uA.
\end{align*}
The matrix exponential limit says that if $D_m\in M(n,\mathbb C)$ satisfies $m(D_m-I)\to B$, then $D_m^m\to\exp(B)$. Applying this with $D_m=C_m$ and $B=uA$ yields
\begin{align*}
C_m^m\to\exp(uA)=\exp(u[X,Y]).
\end{align*}[/guided]