Consider the Lie group $GL(n, \mathbb{R})$. Since $GL(n, \mathbb{R})$ is an open subset of the space $M_{n,n}(\mathbb{R})$ of $n \times n$ matrices, we can identify the tangent space (Lie algebra) $T_{e}GL(n, \mathbb{R})$ with $M_{n,n}(\mathbb{R})$, where $e$ is the identity matrix. I'm trying to understand the following:
For all $x, y \in M_{n,n}(\mathbb{R}) = T_{Id}GL(n, \mathbb{R})$ we have $[x, y] = xy -yx$ where $[\cdot, \cdot]$ is the Lie algebra bracket.
According to the question Is the Lie bracket of a Lie Algebra of a Matrix Lie Group always the commutator?, this seems to be the case for matrix Lie groups in general.
So far, I've worked out the following: Let $x \in M_{n,n}(\mathbb{R}) = T_{Id}GL(n, \mathbb{R})$ and let $X \in \Gamma^{inv}(TG)$ be the left invariant vector field corresponding to $x$, i.e. with $X(e) = x$. By openness, we can choose $\epsilon > 0$ such that the curve $(-\epsilon, \epsilon) \ni t \mapsto e + tx$ through the identity is contained in $GL(n, \mathbb{R})$. For $g \in GL(n, \mathbb{R})$ let $L_g$ be the left multiplication by $g$ in $GL(n, \mathbb{R})$. Using left invariance of $X$ calculate: $$X(g) = ((L_g)_{*}X)(g) = D_eL_g(X(e)) = D_eL_g(x) = \frac{d}{dt}|_{t=0}g(e+tx) = gx$$ Now, the above implies for $f \in C^{\infty}(GL(n, \mathbb{R}))$ and $g \in GL(n, \mathbb{R})$: $$X(g)(f) = (D_gf)(gx)$$ where we identify the tangent vector $X(g)$ with the corresponding derivation at $g$.
Now how do I calculate $[x,y] \stackrel{def.}{=} [X, Y](e) \stackrel{\text{WHY?}}{=}xy-yx$ from this?
I've now been able to do the calculation myself with a bit of inspiration from freeRmodule's answer:
Using the calculations presented in the question, we calculate the adjoint operator for the Lie group $G := GL(n, \mathbb{R})$. For $g \in G$ and $x \in M_{n,n}(\mathbb{R})=T_{Id}G$ with corresponding left invariant vector field $X$ we have: $$Ad(g)(x) = ((R_{g^{-1}})_{*}X)(Id) = D_gR_{g^{-1}}(\overbrace{X(g)}^{=gx}) = gxg^{-1}$$
By a statement which holds for Lie groups in general we then calculate the Lie bracket of $x$ and $y$ (noting that the Lie group exponential $exp_G$ on $G$ is simply the matrix exponential $Exp$): \begin{align*} [x,y] &= \frac{d}{dt}|_{t=0} Ad(exp_G(tx))(y) = \frac{d}{dt}|_{t=0}Exp(tx)yExp(-tx) \\ &= \left( \frac{d}{dt}|_{t=0}Exp(tx) \right)\cdot y \cdot Id + Id \cdot y \cdot \left( \frac{d}{dt}|_{t=0}Exp(-tx) \right) \\ &= xy -yx \end{align*}