Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Math / Wedge World / tensor wedge doc / Support doc files

volume element REVIEWED

DOCX · 48.4 KB
Open DOCX file

Working notes by Phil dated May 2016, marked reviewed in red with all items moved to Appendix F. They start from the 2-form surface area factor K(t)^2 as a sum of squared 2x2 minors, check in Maple that it equals det(RTR), then prove the general theorem for a tall m x n matrix using multiindex and permutation-sum arguments. They also link the result to k-form integration and the tangent-space volume element.

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
The Volume Element Issue PhL 5.3.16 Reviewing in red 5/7/16. All items in this doc have been dealt with, no need ever to reread. I did learn various things here for the first time, all written up now in Appendix F. 5.16.16: Includes the sum of squared minors theorem and implications, all in Appendix F. My first shot at this is in the doc "Measure in Sja" where I look at various possibilities but reach no conclusions. One conjecture presented there is this: The volume element of a dV of the tangent space TxM is given by dt1dt2....dtn [ this is now proven in (F.5.2) ] So I have written up a 2-form type integration and I do get this result for the "volume element" dA = K(t) dt1dt2 where K(t)2 = det2 + det2 + det2 Is it possible that this expression can be written K(t)2 = det(RTR) ? Let's first just try this in Maple to see if by any chance there is a match. The R matrix in this case n = 2 and m = 3 is a 3 x 2 matrix. Here are the Maple results: So the answer is YES: det(RTR) IS another way to write the K2(t) factor! So now I have to learn why this works! And show that it works for any tall m x n R matrix! Note that RTR is an nxn matrix. So consider (T is the matrix transpose, NOT the covariant transpose), K(t)2 = det([RTR]**) = ΣP (-1)S(P) [RTR]1P(1) [RTR]2P(2) ... [RTR]nP(n) = ΣP (-1)S(P) [RTR]1P(1) [RTR]2P(2) ... [RTR]nP(n) = ΣP (-1)S(P)Σi1,i2...in=1m (RT)1i1Ri1P(1)(RT)2i2Ri2P(2) ... (RT)nikRikP(n) = ΣP (-1)S(P)Σi1,i2...in=1m Ri11Ri1P(1)Ri22Ri2P(2) ...RinnRikP(n) Now approach from the other direction in order to find the right threading, and I will stick with this simple example at first, K(t)2 = det2 + det2 + det2 = [ ΣP (-1)S(P) R1P(1)R2P(2)]2 + [ ΣP (-1)S(P) R1P(1)R3P(2)]2 + [ ΣP (-1)S(P) R2P(1)R3P(2)]2 There must be some easier way. Here is what I want to prove: Theorem: Given a tall m x n R matrix, then sum of squares of all minors of width n = det(RTR) // proven now in F.6 Proof: I need some way to write out the LHS. I need notation to indicate a minor. How about this: minorI = minor(ri1, ri2, ....rin) where I is an ordered multiindex! This would enumerate all the minors! Then the claim of the theorem is this ; Σ'I [minorI]2 = det(RTR) In the above, I can think of the rs as n rows selected from the R matrix. I can treat those rows as columns of a different matrix which is just the transpose of the matrix of rows, and the notation above that minorI = det(ri1, ri2, ....rin) is saying the minor is the determinant of these n column vectors. Conjectured Lemma: Σ'I [minorI]2 = (1/n!) ΣI [minorI]2 // part of F.6 Proof: ΣI [minorI]2 = Σi1,i2...in=1m [minorI]2 = Σi1≠i2≠...≠in [minorI]2 = (Σi1<i2<..<in + n!-1 other orderings) [minorI]2 = (ΣP [ΣP(i)<P(i)<...<P(i)]) [minorI]2 = Σi<i<...<i [ΣP fP(i)P(i)...P(i)] . by (A.9.1) with fii...in = [minorI]2 But fii...in = [minorI]2 is a totally symmetric function of the indices since row swaps don't affect the squared determinant. This = Σi<i<...<i [ΣP [minorI]2] = Σi<i<...<i (n!) [minorI]2 = n! ΣI [minorI]2 QED Meanwhile, I know that minorI = ΣP(-1)S(P) Ri1P(1) Ri2P(2) ....RinP(n) = det(RIZ) [minorI]2 = ΣP(-1)S(P) Ri1P(1) Ri2P(2) ....RinP(n) ΣP'(-1)S(P') Ri1P'(1) Ri2P'(2) ....RinP'(n) Then Σ'I [minorI]2 = Σ1≤1i<i2..<in<m ΣP(-1)S(P) Ri1P(1) Ri2P(2) ....RinP(n) ΣP'(-1)S(P') Ri1P'(1) Ri2P'(2) ....RinP'(n) Then at least I have "an expression" for this object. Each term has 2n powers of R! Meanwhile, det(RTR) = ΣP(-1)S(P) (RTR)1P(1) (RTR)2P(2) ....(RTR)nP(n) and this also has 2n powers of R, so they stand a chance of being equal. Let's now write this thing out in gory detail, where I assume the "matrix transpose" and not the covariant transpose. det(RTR) = ΣP(-1)S(P) Σi1,i2...in=1m Ri11Ri1P(1)Ri22Ri2P(2) ...RinnRikP(n) How might I state the desired theorem in multiindex notation? Σ'I [minorI]2 =Σ'IΣP(-1)S(P) Ri1P(1) Ri2P(2) ....RinP(n) ΣP'(-1)S(P') Ri1P'(1) Ri2P'(2) ....RinP'(n) =Σ'IΣP(-1)S(P) RIP(Z) ΣP'(-1)S(P') RIP'(Z) = ΣP(-1)S(P)ΣP'(-1)S(P') Σ'I RIP(Z) RIP'(Z) det(RTR) = ΣP(-1)S(P) Σi1,i2...in=1m Ri11Ri1P(1)Ri22Ri2P(2) ...RinnRikP(n) = ΣP(-1)S(P) ΣI Ri11Ri1P(1)Ri22Ri2P(2) ...RinnRikP(n) = ΣP(-1)S(P) ΣI (Ri11Ri22...Rinn) (Ri1P(1)Ri2P(2) ...RikP(n)) = ΣP(-1)S(P) ΣI RIZ RIP(Z) So here then is what I want to show, expressed in multiindex notation, ΣP(-1)S(P)ΣP'(-1)S(P') Σ'I RIP'(Z) RIP(Z) = ΣP(-1)S(P) ΣI RIZ RIP(Z) Plan A. I think the rearrangement theorem will be called into action here. [ correct ] Consider the left side, LHS = Σ'I ΣP(-1)S(P)RIP(Z) ΣP'(-1)S(P')RIP'(Z) I can replace P' by P'Q and make no difference, so this says LHS = Σ'I ΣP(-1)S(P)RIP(Z) [ΣP'(-1)S(P'Q)RIP'Q(Z)] Now since P is a fixed quantity in the the last bracket, why not take Q = ?? I don't see a way to do anything here. And if I could do something, I would end up with the freed-up ΣP creating n! which I don't want. So maybe this is not the path. Plan B: Maybe this theorem will prove helpful, [ no, was not helpful ] Fact : ΣI TIJ x^I = Σ'I det(TIJ) x^I if TIJ has factored form (A.8.36) ΣI TJI x^I = Σ'I det(TJI) x^I if TJI has factored form Let's back up and try this new form: Σ'I [minorI]2 = Σ'I[ΣP(-1)S(P) Ri1P(1) Ri2P(2) ....RinP(n)[ minorI] ] = Σ'I[ΣP(-1)S(P) Ri1P(1) Ri2P(2) ....RinP(n) det(RIZ) ] = Σ'I ΣP(-1)S(P) RIP(Z) det(RIZ) Then I need to show that ΣP(-1)S(P) Σ'I RIP(Z) det(RIZ) = ΣP(-1)S(P) ΣI RIP(Z) RIZ Plan C. Go back to minorI = minor(ri1, ri2, ....rin) I = ordered multiindex This was just one way to write the minor. Another way might be this Σ'I [minorI]2 = (1/n!) ΣI [minorI]2 which I have now proven above. Then the theorem I need to prove is this: (1/n!) ΣP(-1)S(P)ΣP'(-1)S(P') ΣI RIP'(Z) RIP(Z) = ΣP(-1)S(P) ΣI RIZ RIP(Z) ? and now I have (1/n!) ready to cancel any n! I can develop by playing with perm sums! How about rewriting the above as (1/n!) ΣP(-1)S(P)ΣP'(-1)S(P') ΣI (RT)P'(Z)I RIP(Z) = ΣP(-1)S(P) ΣI (RT)ZI RIP(Z) ? Now do the ΣI normal matrix multiplications to get (1/n!) ΣP(-1)S(P)ΣP'(-1)S(P') (RTR)P'(Z)P(Z) = ΣP(-1)S(P) (RTR)ZP(Z) ? I may be close here! Let RTR = A just to simplify (1/n!) ΣP(-1)S(P)ΣP'(-1)S(P')AP'(Z)P(Z) = ΣP(-1)S(P) AZP(Z) ? Now in the P' sum on the left replace P' by QP' to get (1/n!) ΣP(-1)S(P)ΣP'(-1)S(QP')AQP'(Z)P(Z) = ΣP(-1)S(P) AZP(Z) ? Now since we are inside the ΣP sum, select Q = P to get (1/n!) ΣP(-1)S(P)ΣP'(-1)S(PP')APP'(Z)P(Z) = ΣP(-1)S(P) AZP(Z) ? Now consider that AIJ = Ai1j1Ai2j2..... so it does have "factored form". Then I can say APP'(Z)P(Z) = AP'(Z)Z using (A.8.31) which leaves me with (1/n!) ΣP(-1)S(P)ΣP'(-1)S(PP')AP'(Z)Z = ΣP(-1)S(P) AZP(Z) ? Now swap dummy names P↔P' on the left (1/n!) ΣP'(-1)S(P')ΣP(-1)S(P'P)AP(Z)Z = ΣP(-1)S(P) AZP(Z) ? (1/n!) ΣP' ΣP(-1)S(P)AP(Z)Z = ΣP(-1)S(P) AZP(Z) ? ΣP(-1)S(P)AP(Z)Z = ΣP(-1)S(P) AZP(Z) ? But in this last statement, both sides are representations of det(A)!!! I think I have done it! So let's now pretend I have a solid proof of this theorem: Σ'I [minorI]2 = det(RTR) this plan is mostly what I did in F.6 I have seen now this works out properly in the case of m = 3 and r = 2 which is a special tall R matrix. How did the left side arise in my 2-form case of Section 10.13? It seems a bit specialized there as I look at it The general case for a k-form says ∫φ αx = ∫φ Σ'I fI(x) dxi ^ dxi ^ ... ^ dxi // αx = Σ'I fI(x) dx^I ≡ ∫[0,1]k Σ'J gJ(t) dtj ^ dtj ^ ... ^ dtj // first definition (pull back) ≡ (!Syntax Error, I!Syntax Error, I ... !Syntax Error, I) Σ'J gJ(t) dtjdtj ... dtj // second definition where gJ(t) = Σ'I fI(φ(t)) det(RIJ) and x = φ(t) , R = (Dφ) . (10.12.1) Write all together as ∫φ αx = (!Syntax Error, I!Syntax Error, I ... !Syntax Error, I)Σ'J Σ'I fI(φ(t)) det(RIJ)dtjdtj ... dtj Now det(RIJ) is in fact a minor of R and RIJ is a particular k x k submatrix of R. The measure here seems to be measure = Σ'J det(RIJ)dtjdtj ... dtj In making a connection to theorem Σ'I [minorI]2 = det(RTR), I think k = n is required, because only then are you dealing with an n-dimensional TxM volume element, so let's now restrict to k = n. In that case, ∫φ αx = (!Syntax Error, I!Syntax Error, I ... !Syntax Error, I)Σ'J Σ'I fI(φ(t)) det(RIJ)dtjdtj ... dtjn where det(RIJ) is a full-width minor of R which is an m x n matrix. I think there is only one term in the J' sum so we rewrite the above as ∫φ αx = (!Syntax Error, I!Syntax Error, I ... !Syntax Error, I) Σ'I fI(φ(t)) det(RIZ) dt1dt2 ... dtn Now det(RIZ) = minorI and we get ∫φ αx = (!Syntax Error, I!Syntax Error, I ... !Syntax Error, I) Σ'I fI(φ(t)) [minorI dt1dt2 ... dtn] The theorem says Σ'I [minorI]2 = det(RTR) = (1/n!) ΣI [minorI]2 Look at the general-case 2-form formulas, ∫φ αx = !Syntax Error, I!Syntax Error, I G(t1,t2) dt1dt2 G(t) = Σ'I fI(φ(t)) Now think of this as the dot product of two vectors each having (m,2) components, so G(t) = f(φ(t)) = f(φ(t)) det(RIZ) = f(φ(t)) minor = Σ'I fI(φ(t)) minorI = ( f(φ(t)) minor-hat ) | minor | And then we would say | minor |2 = Σ'I [minorI]2 = det(RTR) and THAT is now it gets into the works. Suppose you had f(φ(t)) = | f(φ(t)) | minor-hat . Then in the case you would have G(t) = | f(φ(t)) | | minor | = | f(φ(t)) | and finally ∫φ αx = !Syntax Error, I!Syntax Error, I | f(φ(t)) | dt1dt2 // this idea now in F.7 Asynchronous Claim: Volume of tangent space is dV = dt1dt2...dtn Is this true? TBC. Continuing 5.5.16