Home / Math and Physics Files / Math / Curvilinear Systems / Tensor Doc and Support / Files related to May 2015 update
How does the determinant of a tensor transform
DOCX · 207.5 KB
Open DOCX file
Phil's exploratory notes dated 4.10.15 in the tensor document update files. Part 1 works the N=2 case using the permutation tensor and finds det(M') = J^2 g^-1 det(M), with no clean result. Later parts use a covariant form of the determinant with two permutation tensors, compare conventions with an outside PDF, and tabulate the weights of det for mixed index placements, concluding weights of 0, -2 and +2.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
How does the determinant of a tensor transform? PhL 4.10.15
Part 1.
Just look at the N = 2 case. We have
det(M) = εabM1aM2b
det(M') ≡ ε'abM'1aM'2b
= εab (R1a' Rab' Ma'b')(R2a" Rbb" Ma"b") ε'ab = εab Weinberg
= (εab Rab'Rbb") (R1a'R2a"Ma'b'Ma"b")
Define K as follows,
εabRab'Rbb" = K εb'b"
because the thing is AS and my theorem says it has that form. Then
εabRa1Rb2 = K ε12 = K
K = εabRa1Rb2 = det(Rij) = det(Sji) = det(S) = J
So we end up with
det(M') = J εb'b" (R1a'R2a"Ma'b'Ma"b")
= J R1a'R2a" (εb'b"Ma'b'Ma"b")
Now study this last object. It is AS on a' and a", so write
(εb'b"Ma'b'Ma"b") = Q εa'a"
Then set a'=1 and a" = 2,
(εb'b"M1b'M2b") = Q ε12= Q g ε12 = Qg // Weinberg; (D.5.10)
so that
det(M) = Q g => Q = g-1 det(M)
We then have
det(M') = J R1a'R2a" (εb'b"Ma'b'Ma"b")
= J R1a'R2a"Q εa'a"
= J g-1 det(M) R1a'R2a" εa'a"
= J g-1 det(M) det(Rij)
= J g-1 det(M) det(Sji)
= J g-1 det(M) J
= J2 g-1 det(M)
Now try
g'α det(M') = J2 g'α g-1 det(M)
= J2 g'α g-1-α gαdet(M)
= J2 g'α(1/g)1+α gαdet(M)
Would like to see 1+α = α which means 1 = 0, so cannot find α to make this work!
So I have not clean answer, and this is a question I have never resolved!
If g = 1, then we have
det(M') = J2 det(M)
and then det(M) is a scalar density of weight - 2. But I don't have a result for the general case!
Part 2.
Comments on Part 1. I think my starting point det(M) = εabM1aM2b is not really covariant, so is a poor starting point. An idea from Part 3 below is that you could write
det(M) = (1/2!) εabεa'b' Ma'a Mb'b
I have not proven this yet, but it I think I will, and I will add it somewhere to tensor doc.
Given this covariant definition of det(M), I can just read off that
det(M') = J2 det(M) scalar density of weight - 2
Then I am completely done and I can just ignore the efforts of Part 1 above.
Part 3.
Search the web. I find in an obscure PDF the following ("Project Report pdf" I saved it)
But he also gives this comment
I am now scanning through this well-written article. Here is how basis vectors transform
`
something new to me
but looks like my general App D notation.
So in his notation, εi means my εabc...x with i being "an index array". σ is a permutation operator I would say, like my P. And the sign is like my (-1)p.
Next comes his definition:
In my notation I have
εabc... = det(gij) εabc... = g εabc... // relating all down to all up
ε'abc... = det(g'ij) ε'abc... = g' ε'abc... // = g' εabc... (D.5.10)
which would say in his notation εi = g εi. He then makes a different underlying definition compared to Weinberg. I now continue in his PDF
I think this would translate into
det(Aij) = (1/N!) εabc.. εa'b'c'... Aaa'Abb'......
compared to my thing
det(Aij) = εa'b'c'... A1a'A2b'......
and I can imagine these would agree. His notation is more covariant.
Assume this is true:
det(Aij) = (1/N!) εabc.. εa'b'c'... Aaa'Abb'......
Then if A were a normal tensor, I would claim that det(Aij) had weight = -2, since each of my ε tensors has weight -1 at least in my world. This agrees with my Part 1 result at least for g = 1, maybe I made a mistake?
He then claims these strange results
strange
He then does tensor densities like this
For me, J =
the Jacobian ≡ J(x') ≡ det(S(x')) = det(∂xi/∂x'k) = 1/det(R(x(x')) = 1/ det(∂x'i/∂xk) . (5.12.6)
so I would translate this something like
Ti'j' = J-W Ri'i Sjj' Tij = J-W Ri'i Rj'j Tij
and this last form agrees with my M'ab = Raa' Rbb' Ma'b' notation. So far I am still with the author.
Now comes his claim quoted above.
Let's now study his following proof. But it is based on his strange claims above.
Alternate Path to Evaluate the transformation property of det(M).
What would happen in my Appendix D if I made his basic assumption that εi = εi ?
That is to say
εabc.. = εabc...
My sections D.1 D.2 D.3 are unchanged since ε has not come up yet.
In D.4 I make this assumption
det(Rij) = R1a' R2b' ... εa'b'c'..
and I guess I now have to think about that. Here ε is the permutation tensor in Weinberg. I think I can write this as
det(Rij) = (1/N!) RAa' RBb' ... εa'b'c'.. εAB...
and then it is seen to be a scalar. That fact is concealed in the first form, but I think the first form is accurate. But, then the question is this: what is the meaning of εAB... ? For Gootvilig εi= εi, but for Weinberg we have instead εAB.. = g εAB... . So maybe I would have to have some g or g' factor out front.
Let's start again.
det(Rij) = R1a' R2b' ... εa'b'c'.. ε... = Weinberg permutation tensor
This is a "mechanical " definition of det (R) because it is what you get actually doing the det!
I think I could show that
det(Rij) = (1/N!) RAa' RBb' ... εa'b'c'.. εAB...
where now there are two permutation tensors. Let's try to prove this right here.
det(Rij) = (1/N!) RAa' RBb' ... εa'b'c'.. εAB...
= (1/N!) εAB...{ εa'b'c'.. RAa' RBb' ... }
= (1/N!) εAB...{ RA1 RB2 ... + all signed permutations of 2nd index }
= (1/N!) εAB... ΣP p P2 (RA1 RB2 ...)
Now define
QAB..X ≡ ΣP p P2(RA1 RB2 ...RXN) (D.10.1)
Argue that this is AS as in App D text, and therefore
QAB..X = K εABC..
Then we have
Q123.. = K ε123.. = K
so then
K = Q123.. = ΣP p P2(R11 R22 ...RNN)
But this is the mech def of the determinant! Thus
K = det(Rij)
and we then have
QAB..X = K εABC.. = det(Rij)εABC..
Going back we then find the thing we started with
det(Rij) = (1/N!) RAa' RBb' ... εa'b'c'.. εAB...
= (1/N!) εAB... QAB..X
= (1/N!) εAB... det(Rij)εABC..
= det(Rij) (1/N!) [ εAB... εAB...]
= det(Rij) (1/N!) [ εAB... ]2
But ΣAB [ εAB... ]2 gets a contribution of +1 each time AB... is a permutation of 12...N. There are N! such permutations, so ΣAB [ εAB... ]2 = N! This then gives
det(Rij) = det(Rij) (1/N!) N! = det(Rij)
and I have finally verified what I thought was true. So we know this is true:
det(Rij) = (1/N!) RAa' RBb' ... εa'b'c'.. εAB...
Now what comes next regarding covariance?? I would use my Weinberg result
εabc... = det(gij) εabc... = g εabc... // relating all down to all up
ε'abc... = det(g'ij) ε'abc... = g' ε'abc... // = g' εabc... (D.5.10)
so that then
εAB... = (1/g) εAB...
and then I would write
det(Rij) = (1/N!) (1/g) εa'b'c'.. εABC.. RAa' RBb' ...
Now suppose R is M, a true tensor, so all the above leads to
det(Mij) = (1/N!) (1/g) εa'b'c'.. εABC.. MAa' MBb' ...
That factor of g really is sitting there if you define det(Mij) in its mechanical sense! So write
g det(Mij) = (1/N!) εa'b'c'.. εABC.. MAa' MBb' ...
I then argue from my usual rules:
g det(Mij) transforms as a scalar density of weight - 2 // if M is a tensor
Now I also know that g' = J2g so g is a scalar density of weight - 2. This would seem then do imply that the object det(Mij) is a true scalar. This is the first time I have ever arrived at this conclusion!
What about other tilts?
Mij = giaMaj => det(Mij) = det(gij)det(Mij) = g det(Mij)
W = -2 W = -2 W=0
Mij = gjkMik => det(Mij) = det(gij)det(Mij) = (1/g) det(Mij)
W = 0 W = +2 W=-2
Mij = gikMkj => det(Mij) = det(gij)det(Mij) = (1/g) det(Mij)
W = 2 W = +2 W=0
So here is a Big Summary of all this det of tensor stuff
det(Mij) W = 0
det(Mij) = g det(Mij) W = -2
det(Mij) = (1/g) det(Mij) W = 0
det(Mij) = (1/g) det(Mij) W = 2
Now rewrite again
det(Mij) = g det(Mij)
det(Mij) = (1/g) det(Mij) = (1/g) g det(Mij) = det(Mij)
det(Mij) = (1/g) det(Mij) = (1/g) det(Mij)
Summarize yet again
det(Mij) = det(Mij) W = 0
det(Mij) = g det(Mij) W = -2
det(Mij) = (1/g) det(Mij) W = +2
Now here is a little check:
g ≡ det(gij) from Ch 5 so weight (g) = -2 agrees
1/g = det(gij) so weight = +2
Everything is finally consistent.
This is an allowed matrix mult form, so I claim that
det(Mij) = det(giaMaj) = det(gij) det(Mij) = g det(Mij).
Then the weight of det(Mij) is -2 + 0 = -2. So far then
det(Mij) = scalar density with weight W = 0
det(Mij) = scalar density with weight W = -2
Next, consider
Mij = giaMaj => det(Mij) = det(gij) det(Mij) = (1/g) det(Mij)
W = 0 W=2 W=-2
The fourth form:
Mij = gjkMik = Mik gkj => det(Mij) = det(Mij) det(gij) = det(Mij)(1/g)
W = 0 W=-2 W= 2
So here I can fill out my chart
First:
det(Mij) = g det(Mij)
det(Mij)
What happens if we make the Gootvilig convention that εabc = εabc ?
In Section D.4 I get down to this general result
ε'abc.. = Kεabc... = J-(W+1) εabc... (D.4.7)
but then I don't know what to do with ε'abc.. !!!
Here are results I end up with the Weinberg convention:
ε'abc.. = εabc... // assumption
εabc... = g εabc...
ε'abc... = g' ε'abc...
εabc.. = gaa'gbb'..... εabc..
ε'abc.. = g'aa'g'bb'..... ε'abc.. // assumption that ε is like any other tensor object
Now how would I "get to" the Loot convention?
εLabc... ≡ εabc...
εLabc... ≡ (1/g)εabc.. = (1/g) g εabc... = εabc...
Then I would have that εLabc... transforms with weight -1, but εLabc... would have weight 2 - 1 = +1
and that must be his method.
Question 1: What about G = det(g) ? Here g is a true rank-2 tensor. So by my claim above, G should be a normal scalar, but in fact it is a scalar density of weight - 2. What have I done wrong? Maybe this is the way out : G = det(gij) and not det(gij).
7. Jacobian and related
Here are translations of a few equations from Section 5.12:
(5.12.2) det(S) = 1/det(R) → det(Sij) = 1/ det(Rij) (7.5.18)
(5.12.6) J = det(S) → J = det(Sij) (7.5.19)
(5.12.12) g = det() → g = det(gij)
g' = det(') → g' = det(g'ij) (7.5.20)
(5.12.14) 1/g = det(g) → 1/g = det(gij)
1/g' = det(g') → 1/g' = det(g'ij) (7.5.21)
(5.12.14) g' = J2 g → g' = J2 g
|J| = → |J| = (7.5.22)