The Metric Matrix: How Can I Invert a Non-Diagonal Matrix?

  • Context: Graduate 
  • Thread starter Thread starter MathematicalPhysicist
  • Start date Start date
  • Tags Tags
    Metric
Click For Summary

Discussion Overview

The discussion revolves around the challenges of inverting a non-diagonal metric tensor in the context of general relativity, specifically focusing on the metric ##ds^2=dudv+F(y,z)du^2+dy^2+dz^2##. Participants explore the implications of this inversion for calculating Christoffel symbols, Ricci tensor, and Ricci scalar, while considering various methods and approaches to tackle the problem.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant expresses the need to calculate Christoffel symbols and Ricci tensor but is unsure how to find the inverse of a non-diagonal metric without inverting it directly.
  • Another participant suggests that the metric has a straightforward inverse and questions the necessity of avoiding direct inversion.
  • A participant provides the matrix representation of the metric and notes that inverting it requires linear algebra techniques.
  • Some participants propose alternative methods for computing the Levi-Civita connection coefficients, suggesting that using geodesic equations may simplify the process.
  • There are discussions about the specific non-zero Christoffel symbols and Ricci tensor components, with some participants providing their calculations and others questioning or refining these results.
  • One participant mentions using a Maple package for calculations and seeks assistance with its implementation.
  • Another participant confirms the vanishing of the Ricci scalar based on their calculations.
  • Discrepancies arise regarding the values of certain Christoffel symbols, with participants offering different constants and methods for their derivation.

Areas of Agreement / Disagreement

Participants express differing views on the necessity and complexity of inverting the metric directly. While some agree that the inversion is straightforward, others suggest alternative methods that may be less tedious. There is also no consensus on the exact values of certain Christoffel symbols, indicating ongoing debate and refinement of calculations.

Contextual Notes

Some calculations and assumptions regarding the metric and its properties remain unresolved, particularly concerning the specific forms of the Christoffel symbols and the Ricci tensor components. The discussion highlights the dependence on definitions and the potential for different interpretations of the results.

Who May Find This Useful

This discussion may be useful for students and researchers in general relativity, particularly those interested in the mathematical techniques for handling non-diagonal metrics and the computation of geometric quantities in curved spacetime.

MathematicalPhysicist
Science Advisor
Gold Member
Messages
4,662
Reaction score
372
I have the metric ##ds^2=dudv+F(y,z)du^2+dy^2+dz^2##.

I have shown that for ##F(y,z)=0## it's Minkowski metric, but for ##F(y,z)\ne 0 ## I want to calculate Christoffel Symbols, Ricci tensor and scalar, the problem is that the metric ##g_{\mu\nu}## is not diagonal, and I need to find ##g^{\mu\nu}## in order to calculate christoffel and the rest, but how can I solve this without inverting the metric?

Thanks!
 
Physics news on Phys.org
Why do you need to do it without inverting the metric? The metric has a straight-forward inverse.
 
Well we have a matrix: {0,1,0,0},{1,F(y,z),0,0},{0,0,1,0},{0,0,0,1}

So it's not diagonal, so inverting it you need to use the algorithm from Linear Algebra.
 
Ah, yes, the inverse is easy, didn't notice it.
:-D
 
So? Is that a problem?

EDIT: Still, if it wasn't easy I see no problem in using linear algebra. That's why you learned it - to use it.
 
Also, I strongly suggest against using the typical expression in terms of the metric to compute the connection coefficients of the Levi-Civita connection. Instead, use the geodesic equations obtained from finding the optima of
$$
\int g_{\mu\nu}\dot x^\mu \dot x^\nu ds.
$$
It is the same maths, but the bookkeeping is way easier.
 
weirdoguy said:
So? Is that a problem?

EDIT: Still, if it wasn't easy I see no problem in using linear algebra. That's why you learned it - to use it.
I know, but if there are less tedious methods just in case...
 
  • Like
Likes   Reactions: weirdoguy
BTW, does someone know how to implement the calculation of Christoffel symbol, Ricci tensor and Ricci scalar in maple?

I just want to check if my calculations are right, I got:

The only Christoffel symbols that don't vanish are: ##\Gamma_{uu}^y=-F'_y(y,z)/2## and ##\Gamma_{uu}^z=-F'_z(y,z)/2##; and the Ricci tensor the components that don't vanish are:
$$R^y_{uyu}=-F''_{yy}/2$$
$$R^z_{uzu}=-F''_{zz}/2$$
$$R^y_{uzu}=R^z_{uyu}=-F''_{yz}/2$$

I get that Ricci scalar vanishes since ##R=R^{\alpha}_{\beta\mu\alpha}g^{\beta\mu}=0##, am I right?

I don't understand how to use this package in maple 2017, I asked for help in a suitable forum for maple, but didn't get yet an answer.
 
Did you get the point of #6? It really is the most straight-forward way of computing the Levi-Civita connection coefficients.

MathematicalPhysicist said:
The only Christoffel symbols that don't vanish are: ##\Gamma_{uu}^y=-F'_y(y,z)/2## and ##\Gamma_{uu}^z=-F'_z(y,z)/2##;

I get more non-zero Christoffel symbols. In particular, look for ##\Gamma^v_{uy} = \Gamma^v_{yu}## and ##\Gamma^v_{uz} = \Gamma^v_{zu}##. Can you show your work?

and the Ricci tensor the components that don't vanish are:
$$R^y_{uyu}=-F''_{yy}/2$$
$$R^z_{uzu}=-F''_{zz}/2$$
$$R^y_{uzu}=R^z_{uyu}=-F''_{yz}/2$$
That's the Riemann tensor. The Ricci tensor is ##R_{\mu\nu} = R^{\gamma}_{\mu\gamma\nu}##. There are also some non-zero elements of the form ##R^v_{abc}##, where ##abc## are all different from ##v##, but they do not affect the computation of the Ricci scalar or the Ricci tensor since they never correspond to one of the lower indices and therefore disappear in the trace. Therefore:

I get that Ricci scalar vanishes since ##R=R^{\alpha}_{\beta\mu\alpha}g^{\beta\mu}=0##, am I right?
Yes. I get a vanishing Ricci scalar as well.
 
  • #10
Yes, I get that ##\Gamma_{uy}^v=1/4F_y'## and ##\Gamma^v_{uz}=1/4F_z'##, I think this about sums it for Christoffel symbol.

I can scan my work, I just inverted the metric with an LA algorithm (where you augment ##g_{\mu\nu}## with the identity matrix and then get on the left the identity matrix and on the right the inverse matrix).
 
  • #11
MathematicalPhysicist said:
BTW, does someone know how to implement the calculation of Christoffel symbol, Ricci tensor and Ricci scalar in maple?

Some years ago I used a Maple package called GRTensor II which I ran under the symbolic toolbox in Matlab.

Cheers
 
  • #12
MathematicalPhysicist said:
Yes, I get that ##\Gamma_{uy}^v=1/4F_y'## and ##\Gamma^v_{uz}=1/4F_z'##, I think this about sums it for Christoffel symbol.
I got a different constant. Using the method described in #6 your work for the Christoffel symbols should not be longer than a few lines (one for the variation with respect to each coordinate). The relevant variations are the variations with respect to ##u## and ##v## of
$$
\int \underbrace{[\dot u \dot v + F \dot u^2 + \dot y^2 + \dot z^2]}_{\equiv \mathscr L} ds.
$$
 
  • #13
MathematicalPhysicist said:
Well we have a matrix: {0,1,0,0},{1,F(y,z),0,0},{0,0,1,0},{0,0,0,1}

So it's not diagonal, so inverting it you need to use the algorithm from Linear Algebra.
Its easy to invert$$\left(
\begin{array}{cc|c|c}
0 & \tfrac12 & 0 & 0 \\
\tfrac12 & F & 0 & 0 \\
\hline 0 & 0 & 1 & 0 \\
\hline 0 & 0 & 0 & 1
\end{array}\right)
$$provided you remember how to invert$$
\begin{pmatrix}
0 & \tfrac12 \\
\tfrac12 & F
\end{pmatrix}
$$
 
  • Like
Likes   Reactions: MathematicalPhysicist

Similar threads

  • · Replies 62 ·
3
Replies
62
Views
7K
  • · Replies 24 ·
Replies
24
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
Replies
1
Views
2K
  • · Replies 76 ·
3
Replies
76
Views
4K
  • · Replies 6 ·
Replies
6
Views
3K