Este OPE no es trivial, simplemente no tiene términos singulares.
Supongamos que hay un operador cuasi-primario $\mathcal{O}$ con peso $(h,\bar h)$ que aparece en el lado derecho. Podemos calcular el coeficiente con el que aparece al observar la función de tres puntos $$ \langle T(z)\bar T(\bar w)\mathcal{O}(x,\bar x)\rangle=\frac{f_{T\bar T\mathcal{O}}}{(z-w)^{2-h}(\bar z-\bar w)^{2-\bar h}(x-w)^{h-2}(\bar x-\bar w)^{\bar h+2}(z-x)^{2+h}(\bar z-\bar x)^{\bar h-2}}. $$ El lado derecho está fijado por la invariancia conforme global hasta el coeficiente $f_{T\bar T\mathcal{O}}$. Sin embargo, el lado izquierdo solo depende de $z$ y no de $\bar z$, así que debemos concluir que $\bar h=2$. De manera similar, debido a que solo depende de $\bar w$ y no de $w$, debemos concluir que $h=2$. Esto significa que no pueden aparecer términos singulares en el OPE porque necesariamente deben tener $h+\bar h<2$. Pero entonces podemos definir el operador $$ (T\bar T)(z,\bar z) \equiv T(z)\bar T(\bar z). $$ Es cuasi-primario y tiene dimensiones $(h,\bar h)=(2,2)$. De hecho, es el único cuasi-primario que aparece en el OPE. El OPE toma la forma simplemente $$ T(z)\bar T(\bar w) = (T\bar T)(z,\bar w)=\sum_{n=0}^\infty \frac{1}{n!}(z-w)^n\partial^n_w(T\bar T)(w,\bar w). $$
El operador $T\bar T$ de hecho puede ser definido en cualquier QFT 2D, no necesariamente conforme, pero el argumento es más sutil. Actualmente hay mucha investigación en las teorías que se obtienen al añadir $T\bar T$ a la lagrangiana. (Nota que esta es una deformación irrelevante). Busca "$T\bar T$-deformación".
Agregado: A pedido del OP, aquí hay un curso intensivo en OPEs en teorías de campo conformes. Dado que esta respuesta solo requiere invariancia conforme global, no discutiré las implicaciones de la simetría de Virasoro. Debido a esto, lo siguiente se aplica (con pequeñas modificaciones para acomodar el espin general) en CFT en $d\geq 2$. La simetría de Virasoro también conduce a modificaciones sencillas.
A continuación, $x_i$ denotan puntos espacio-temporales.
Cualquier CFT posee una expansión del producto de operadores que converge en estado de vacío. Es decir, $$ \mathcal{O}_1(x_1)\mathcal{O}_1(x_2)|0\rangle=\sum_i f_{\mathcal{O}_1\mathcal{O}_2\mathcal{O}_i} C_{12i}(x_1,x_2,x_3,\partial_{x_3})\mathcal{O}(x_3)|0\rangle. $$ El punto $x_3$ es en principio arbitrario y a menudo se toma como $x_3=x_2$. Aquí, el operador diferencial $C_{12i}(x_1,x_2,x_3,\partial_{x_3})$ está completamente fijado por la simetría conforme. Depende únicamente de los números cuánticos de los operadores $\mathcal{O}_1,\mathcal{O}_2,\mathcal{O}_i$. El coeficiente $f_{\mathcal{O}_1\mathcal{O}_2\mathcal{O}_i}$ no está fijado por la simetría conforme y representa la información dinámica de la teoría.
Esta expansión es exacta y convergente. A menudo se escribe omitiendo el estado de vacío $|0\rangle$. Esto se debe a que a menudo se utiliza dentro de funciones de correlación euclídeas, donde uno no necesariamente tiene que hablar de una cuantización particular. En funciones de correlación euclídeas, se interpreta el estado de vacío en cuantización radial alrededor del punto $x_3$. El OPE es aplicable en una función de correlación de $n$ puntos euclídea si existe una esfera alrededor de $x_3$ que solo contiene los operadores $\mathcal{O}_1,\mathcal{O}_2$ en $x_1$ y $x_2$ y ningún otro operador.
Uno puede calcular el coeficiente $f_{\mathcal{O}_1\mathcal{O}_2\mathcal{O}_i}$ observando la función de tres puntos $\langle\mathcal{O}_1\mathcal{O}_2\mathcal{O}_i\rangle$ y utilizando el OPE dentro de la función de tres puntos. Dado que las funciones de dos puntos se eligen canónicamente para ser diagonales $\langle\mathcal{O}_i\mathcal{O}_j\rangle\propto \delta_{i,j}$, tenemos $$ \langle\mathcal{O}_1(x_1)\mathcal{O}_2(x_2)\mathcal{O}_i(x_3)\rangle=f_{\mathcal{O}_1\mathcal{O}_2\mathcal{O}_i} C_{12i}(x_1,x_2,x'_3,\partial_{x'_3})\langle\mathcal{O}_i(x'_3)\mathcal{O}_i(x_3)\rangle. $$ Nuevamente, a menudo se usa $x'_3=x_2$. Dado que $C_{12i}(x_1,x_2,x'_3,\partial_{x'_3})\langle\mathcal{O}(x'_3)\mathcal{O}(x_3)\rangle$ está fijado por la simetría conforme y la normalización canónica de las funciones de dos puntos, el coeficiente $f_{\mathcal{O}_1\mathcal{O}_2\mathcal{O}_i}$ se calcula mediante las funciones de tres puntos. Sin embargo, este coeficiente aparece en el OPE, y el OPE es aplicable en todas las funciones de correlación de $n$ puntos, así que no hay falta de generalidad por la que el OP parece estar preocupado.