I am trying to understand the application of the operator product expansion to calculate the radially ordered product in the complex plain of Tzz(z)∂wXρ(w) which should result in
⟨R(Tzz(z)∂wXρ(w))⟩=−l2s1(z−w)2∂wXρ(w)−l2s1(z−w)∂2zXρ(z)+⋯
but embarassingly I encounter a stumbling block right at the beginning. After inserting Tzz(z)≐:ημν∂zXμ∂zXν: one has
⟨R(Tzz(z)∂wXρ(w))⟩=R(:ημν∂zXμ(z)∂zXν(z):∂wXρ(w))
which can obviously be further expanded to
...=ημν⟨∂zXμ(z)∂wXρ(w)⟩∂zXν(z)+ημν⟨∂zXν(z)∂wXρ(w)⟩∂zXμ(z)
It is exactly this last step I dont understand. If this initial stumbling block is removed, I can understand the rest of the derivation, so can somebody help me remove it?
To generalize a bit, it seems I do not yet properly understand how such expressions involving normal and radial (time) ordered products are generally evaluated. So if somebody could give me a more general hint about this, I would probably be able to see how the last expression in my particular example is obtained.