Lo demostraré por cada $X:M \rightarrow TM$ campo vectorial, podemos asociar un operador diferencial $X:C^{ \infty }(M) \rightarrow C^{ \infty }(M)$ .
A estas alturas probablemente sepas que podemos trabajar en coordenadas $(U, \mathcal {X})$ , $U \underset {open} \subset M$ , $ \mathcal {X}:U \rightarrow \mathbb {R}^n$ y así podemos trabajar con el campo vectorial $X$ escrito localmente como $X= \sum_ {i=1}^nX_i \partial x_i$ . Es decir, estamos escribiendo $ \dfrac { \partial }{ \partial x_i}(p)= \partial x_i(p)$ por conveniencia y borrando el $p$ para evitar la anotación pesada. También recuerde que $\{ \partial x_1(p),..., \partial x_n (p)\}$ es una base para $T_p M$ .
Ahora, dado este campo vectorial, define un mapeo $ G:\mathfrak {X}^ \infty (M) \rightarrow \mathcal {D}^ \infty (M)$ desde el espacio de los campos vectoriales $ \mathfrak {X}^ \infty (M)$ sobre $M$ al conjunto de "operadores diferenciales direccionales $ \mathcal {D}^ \infty (M)$ (Acabo de inventar esto) de mapeos $ \mathcal {C}^ \infty (M) \rightarrow \mathcal {C}^ \infty (M)$ de tal manera que $$ G(X)(f)= \sum_ {i-1}^n X_i \dfrac { \partial f}{ \partial x_i} $$ que está en $C^ \infty (M)$ (recuerde que cada $X_i$ está en $C^ \infty (M)$ ). Básicamente, $G$ es sólo el "cambio de punto de vista" que permite ver un campo vectorial $X$ como un derivado direccional $X$ .
Ahora, dado un operador diferencial direccional, se puede definir un campo vectorial que está asociado a él. Espero que ya se hayan dado cuenta ahora de que la notación para la base de $T_p M$ no es accidental y está motivada por esta asociación entre los campos vectoriales y la diferenciación direccional.
Si necesitas una referencia más fuerte, intenta buscar la "Geometría de Riemann" de Do Carmo. También estoy tomando este curso este semestre y este libro está siendo de gran ayuda. Si tienes alguna pregunta, no dudes en preguntar.