Hello,
I start from a certain lagrangian containing covariant derivative of a tensor field X[-a,-b], for example L=Sqrt[detG[]] CD[-a][X[-b,-c]] D[a][X[b,c]], where X is symmetric and a background field "no perturbations".
I am looking for the first and second perturbation order of this lagrangian with respect to the metric only and then to isolate the Tensor, Vector and Scalar sectors.
First I used DefMetric, SetSlicing, DefMetricField[h,dg,h] and DefProjectedTensor[X[-a,-b], TensorProperties->{"SymmetricTensor"},SpaceTimeOfDefinition->{"Background"}].
Then I can do Perturbed[L,1]//ExpandPerturbation etc..
But here I obtain an expression with both dg[1,a,b] and h[1,a,b], and only the 4D covariant Derivatives
This is a thing i don't understand.
But after that I would like to project perturbations on a flat FLRW background, in order to split my perturbations on the various sectors, but here i end up again with both dg^[1,a,b] and h^[1,a,b], prime derivatives and the 4d covariant derivatives of dg and h.
I suppose, if things had worked normally, i should have terms like X ' [a,b] X ' [-a,-b] + E' [a,b] E ' [-a,-b] X[c,d] X[-c,-d] etc...
Could you please indicate me how to get rid of that problem and have a good splitting.
Thank you for your help.
Charles