Metrics and coordinate transforms

[chroot]

Hi all, I'm trying to solve Exercise 1.4.3 in Foster & Nightingale's "A Short Course in General Relativity."

The question essentially provides the metric for Euclidean space in spherical coordinates and the matrix representing the coordinate transformation from spherical to cylindrical coordinates. The question prompts the deduction of the metric for Euclidean space in cylindrical coordinates using this information and the rule

g_{i^\prime j^\prime}<br /> = U^k_{i^\prime}<br /> U^l_{j^\prime}<br /> g_{kl}<br />

where primed coordinates are cylindrical, and unprimed coordinates are spherical.

The matrices involved are

1) the metric of Euclidean space in spherical coordinates

[g_{ij}] =<br /> \left[<br /> \begin{array}{ccc}<br /> 1 &amp; 0 &amp; 0\\<br /> 0 &amp; r^2 &amp; 0\\<br /> 0 &amp; 0 &amp; r^2 \sin^2 \theta<br /> \end{array}<br /> \right]<br />

2) the coordinate transformation matrix between spherical and cylindrical coordinates:

[U^{i^\prime}_j] =<br /> \left[<br /> \begin{array}{ccc}<br /> \sin \theta &amp; r \cos \theta &amp; 0\\<br /> 0 &amp; 0 &amp; 1\\<br /> \cos \theta &amp; -r \sin \theta &amp; 0<br /> \end{array}<br /> \right]<br />

Now, it seems to me that, in matrix notation:

<br /> [g_{i^\prime j^\prime}] = <br /> <u> [g_{ij}] [U^T]<br /> </u>

But this doesn't seem to be correct. I should end up with the metric for Euclidean space in cylindrical coordinates:

<br /> g_{i^\prime j^\prime} = <br /> \left[<br /> \begin{array}{ccc}<br /> 1 &amp; 0 &amp; 0\\<br /> 0 &amp; r^2 \sin^2 \theta &amp; 0\\<br /> 0 &amp; 0 &amp; 1<br /> \end{array}<br /> \right]<br /> <br /> =<br /> <br /> \left[<br /> \begin{array}{ccc}<br /> 1 &amp; 0 &amp; 0\\<br /> 0 &amp; \rho^2 &amp; 0\\<br /> 0 &amp; 0 &amp; 1<br /> \end{array}<br /> \right]<br />

If my thinking *is* correct, somewhere I'm just getting confused. Can anyone help me figure out what's going on? Thanks in advance.

- Warren

 

[chroot]

Er wait...

[g_{ij}] = [g_{ji}]

via matrix transposition, yes?

So then

[U^{i^\prime}_j] = [U^j_{i^\prime}]

via matrix inversion, not transposition, yes?

Oooops... I still am not getting the correct answer, but I suspect I should have written

<br /> [g_{i^\prime j^\prime}] = <br /> <u> [g_{ij}] [U^{-1}]<br /> </u>

- Warren

 

[chroot]

Erm, no one can help me?

- Warren

 

[Hurkyl]

Your coordinate change matrix doesn't look right...

 

[chroot]

Hurkyl,

Any idea what it should be then? That came right out of my book. :-/

- Warren

 

[lethe]

Originally posted by chroot
Hurkyl,

Any idea what it should be then? That came right out of my book. :-/

- Warren

I agree with Hurkyl. the matrix that i got was
\left (\begin{array}{ccc}<br /> \frac{r}{\sqrt{r^2+z^2}} &amp; \frac{z}{\sqrt{r^2+z^2}} &amp; 0\\<br /> \frac{z}{\sqrt{r^2+z^2}} &amp; \frac{-r}{\sqrt{z^2+r^2}} &amp; 0\\<br /> 0 &amp; 0 &amp; 1<br /> \end{array}<br /> \right)<br />

 

[lethe]

oh, actually, the thing i wrote is the same thing that you had... its OK

anyway. the point of the story is you don t conjugate to transform the metric. you only conjugate to transform operators

 

[lethe]

the transformation rule is
<br /> g&#039;_{\mu\nu}=\frac{\partial x^\rho}{\partial x&#039;^\mu}\frac{\partial x^\sigma}{\partial x&#039;^\nu} g_{\rho\sigma}<br />

 

[selfAdjoint]

Courtesy of Lillian Lieber, "Co, low, primes below". That is, covariant indices are subscript, and the transformation derivatives are with respect to the new, primed coordinates. And of course contravariant indices are just the opposite.

 

[lethe]

Originally posted by chroot

<br /> [g_{i^\prime j^\prime}] = <br /> <u> [g_{ij}] [U^{-1}]<br /> </u>

- Warren

a minor point, but i would like to correct your notation. \mbox{$ U^i_j$} is a matrix element. to represent the matrix that these components make up you can either put the single component in brackets, getting something like \mbox{$\left [ U^i_j\right ] $} (which is what you have done for the metric) or you can use \mbox{$ U $} to stand for the whole tensor. but to put \mbox{$ U $} inside the brackets is meaningless.

 

[chroot]

Originally posted by lethe
the transformation rule is
<br /> g&#039;_{\mu\nu}=\frac{\partial x^\rho}{\partial x&#039;^\mu}\frac{\partial x^\sigma}{\partial x&#039;^\nu} g_{\rho\sigma}<br />

Yeah, I'm aware. And in matrix notation, doesn't that end up being

[g_{i^\prime j^\prime}] = <br /> [U^{i^\prime}_k] [g_{kl}] [U^l_{j^\prime}]

So that the matrix multiplication is correct? Maybe I'm being naive, but I'm just flipping the U's upside down (inverting the corresponding matrix) and reordering terms so that each multiplication has the same index for the columns of the first matrix as the rows of the second.

Is this not the right way to approach setting up this kind of a multiplication?

Is there some deeper, darker secret I need to know to properly evaluate

g_{i^\prime j^\prime}<br /> = U^k_{i^\prime}<br /> U^l_{j^\prime}<br /> g_{kl}

using matrices? The tensor forms make perfect sense to me, but using them as matrices seems to be a stumbling block to me.

- Warren

 

[lethe]

Originally posted by chroot
Yeah, I'm aware. And in matrix notation, doesn't that end up being

[g_{i^\prime j^\prime}] = <br /> [U^{i^\prime}_k] [g_{kl}] [U^l_{j^\prime}]

So that the matrix multiplication is correct?

no, matrix multiplication is not correct. matrix multiplication is the rule for multiplying a (1,1) tensor with another (1,1) tensor.

when you change coordinates, a (1,1) tensor gets conjugated, but not a (0,2) tensor. the metric is a (0,2) tensor.

forget matrices. thinking of tensors as matrices is a barrier to understanding what a tensor is.

also, the formula you wrote above, i can tell just by looking at it that it is not correct, since the indices don t match on both sides of the equation.

 

[chroot]

Originally posted by lethe
when you change coordinates, a (1,1) tensor gets conjugated, but not a (0,2) tensor. the metric is a (0,2) tensor.

Conjugation = matrix transposition? Inversion?

forget matrices. thinking of tensors as matrices is a barrier to understanding what a tensor is.

I can't forget matrices. I'm trying to understand this exercise, which asks me to deal with matrices. Further, how can you actually compute something with tensors if you don't consider their components?

- Warren

 

[lethe]

Originally posted by chroot
Conjugation = matrix transposition? Inversion?

conjugation means ABA^-1. here, B is your tensor, and A is the change of coordinates matrix. this is the operation you were trying to do to the metric, which i claim is only valid for (1,1) tensors.

I can't forget matrices. I'm trying to understand this exercise, which asks me to deal with matrices. Further, how can you actually compute something with tensors if you don't consider their components?

- Warren

alright, keep your matrices. i don t find them so useful, but perhaps you can get more use from them than i do.

 

[lethe]

OK, so let my spherical coordinates be r\mbox{, }\theta\mbox{, and } \phi, and my cylindrical coordinates be \rho, \psi, z

\frac{\partial x^\text{sph}}{\partial x^\text{cyl}} = \left(\begin{array}{ccc}<br /> \frac{\rho}{\sqrt{\rho^2+z^2}} &amp; \frac{z}{\sqrt{\rho^2+z^2}} &amp; 0\\<br /> \frac{z}{\rho^2+z^2} &amp; \frac{-\rho}{\rho^2+z^2} &amp; 0\\<br /> 0 &amp; 0 &amp; 1<br /> \end{array}\right)<br />
where the first column is differentiation with respect to &rho;, the second with respect to z, and the third with respect to &psi;

the metric in spherical is

<br /> g^\text{sph} = \left(\begin{array}{ccc}<br /> 1 &amp; 0 &amp; 0\\<br /> 0 &amp; r^2 &amp; 0\\<br /> 0 &amp; 0 &amp; r^2\sin^2\theta<br /> \end{array}\right)<br />

so using the transformation formula, i can get the elements of g^cyl

i will only calculate the diagonal components, since the basis is orthogonal.

<br /> \begin{align*}<br /> g^\text{cyl}_{\rho\rho}&amp;=\frac{\partial r}{\partial \rho}\frac{\partial r}{\partial \rho}g^\text{sph}_{rr}+\frac{\partial\theta}{\partial\rho}\frac{\partial\theta}{\partial\rho}g^\text{sph}_{\theta\theta}+\frac{\partial\phi}{\partial\rho}\frac{\partial \phi}{\partial\rho}g^\text{sph}_{\phi\phi} \\<br /> =&amp; \frac{\rho^2}{\rho^2+z^2} + \frac{z^2}{(\rho^2+z^2)^2}(\rho^2+z^2)<br /> = 1<br /> \end{align*}<br />

ok...

<br /> g^\text{cyl}_{\psi\psi}=g^\text{sph}_{\phi\phi}=\rho^2<br />

and
<br /> g^\text{cyl}_{zz}=\frac{z^2}{\rho^2+z^2}+\frac{\rho^2}{(\rho^2+z^2)^2}(\rho^2+z^2)=1<br />

and thus, in cylindrical coordinates, i have

<br /> ds^2=d\rho^2+\rho^2d\psi^2+dz^2<br />

voil&agrave;