First of all this is not a quantum gravity problem and to be fair not a problem. We take a torsion free theory because torsion is not observed as far as I know. Furthermore you do not need Dirac fermions in order that torsion is not vanishing! Also you can have fermions and arrange them such that torsion vanishes! Simply by introducing the cotorsion in your vierbein connection and write your new connection as
˜ωabμ=ωabμ−Kabμ
where
Kρμν=(Tρμν+Tρμν+Tρνμ)/2
Of course C is the torsion tensor defined as Tρμν=2Γρ[μν]. Then ˜ωabμ only depends on the vierbein eμa as
˜ωabμ=eaρ∂[eebρ]+12eaρebσecμ∂[σecρ]−(a↔b)
The cotorsion does not contribute to the Lorentz transformations at all ( of Ωμ actually, too much background I need to write to make it precise). It turns out that you can choose ωabμ such that the torsion vanishes in a gravitational theory that contains spinorial fields.
Background:
But to the point, on why fermions induce torsion in GR. First of all let us remember that GR is invariant under the big group GL(4,R) that people sometimes call it "symmetry of the world" (yes I have heard it indeed). The not necessarily symmetric Christoffel connection Γ then belongs to its associated algebra which is gl(2,R). Good. One important thing we require that any realistic theory of gravity has is the metric compatibility condition, i.e. ∇λgμν=0, something you will most certainly agree upon. Remember that this covariant derivative involves a lot of Christoffel symbols and it acts on a vector field Vμ as ∇νVμ=∂νVμΓμλνVλ. Another thing we have to remember is what is the infinitesimal length in GR. This is given by the usual ds2=gμνdxμdxν but we can also use the famous "vierbeins" eμa which form an orthonormal frame that transforms under both GL(4,R) and the connected component of the Lorentz group SO(3,1). Then the metric can be written in therms of the vierbeins as gμν=ημνeμaeνb. Now, this metric is transforming as a usual two-ranrk tensor under GL(4,R) and as a scalar under SO(3,1). In correspondence with the connection of the metric we also need a connection for this new orthtonormal frame. This is given by the connection ωabμ. It turns out that its possible to define the covariant derivative of any field representation of SO(3,1) using ωabμ, for example for a collection of scalar fields we get that the covariant derivative of acts as DμϕI=∂μϕI+(Ωμ)IJϕJ
. Here, we have introduced a new tensor Ωμ≡−i2ωabμSab. Of course Sab are the usual spinor matrices which are given by Sab=Σab/2=i4[γa,γb] and γa are the usual Dirac γ-matrices that form a Clifford algebra. In any case we will be intrested in the spinor representation of the Lorentz group and how the covariant derivative acts on it: Dμψ=∂μψ−i4ωabμΣabψ.
The curious thing is on the way the Ωμ transforms. It turns out that in order to transform properly the new connections ωμab must obey
ωabμ→ω′abμ=ωabμ+λacωcbμ+λbcωcaμ−∂μλab
where λab=−λba are the usual transformation parameters (the ones we find in the exponential). Very good, now that we have established the basics lets go to the main point. The covariant derivative ∇μ needs to be generalized in order to incorporate the new connections. The new covariant derivative can be written as Dν and acts on the vierbeins as
Dνeaμ=∂νeaμ+Γμρνeaρ+ωabνeμb=0
and this condition determines the form of the spin connection ωabμ=gρνeaρ∇μebν with thisoption not being unique though. How does torsion come into the game? Well, let us see what happens with the Riemann curvature tensor first. The curvature tensor should transform as a GL(4,R) tensor, as usual, and as a second rank Lorentz tensor now. This can be written in terms of the new connection as
Rabμν=∂μωabν+ωacνωbνc−∂νωabμ−ωacνωbμc
simply by projecting out 2 components of a 4-rank tensor as following Rρσμν=eρaebσRabμν. If you replace this to into (Curv) you will get back the usual Riemann tensor. Now note that we have everything we need in order to formulate a gravitational theory that includes Lorentz spinors (which transform as world scalars). The catch is, as we said in the beginning, that the Christoffel symbols might not be symmetric. The antisymmetric part of the Christoffel symbol transforms as a GL(4,R) tensor and is written as
Tσμν=2Γσ[μν]
Then, if the right hand side is non-vanishing the Einstein tensor also becomes assymetric and as a result the stress-energy-momentum tensor.