Editing Jaycall/KillingVectorApproach

__FORCETOC__ 
<!-- __NOTOC__ will force TOC off -->

=Characteristic Vector Approach to Integrals of Motion Problem=

Motivated by our recent work in relativistic fluid dynamics, one insightful approach to the integrals of motion problem involves rewriting the equations of motion in terms of an arbitrary vector field, which we're currently calling a ''characteristic vector''.  The benefit here is that we can produce a mathematical condition that tells us what combintations of variables should be grouped together in order to produce conserved quantities.  The details follow.

==Newtonian Equations of Motion==

Using index summation notation (wherein repeated indices are summed over), it is possible to write the Newtonian equations of motion for a test particle in a fixed, static potential in an arbitrary orthogonal coordinate system.  The equation describing motion in the <math>i</math>-dimension is:

<div align="center">
<math>
\frac{d}{dt} \left( m \ {h_i}^2 \dot{\lambda}_i \right) = m \ {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k - m \ \partial_i \Phi
</math>
</div>

where <math>m</math> is the mass of the test particle, the <math>\lambda</math>s are the chosen coordinates, and the <math>\dot{\lambda}</math>s are the coordinate velocities (or the total time-derivatives of the coordinates), the <math>h</math>s are the metric scale factors associated with each of the coordinates, <math>\Phi</math> is the gravitational potential, and the <math>\Gamma</math>s are the Christoffel symbols. 

==Metric Components, Orthogonality Condition, and Scale Factors==

For any coordinate system it is possible to express the components of the metric in terms of the partial derivatives defining the transformation from any other coordinate system.  Written in terms of the transformation from cylindrical coordinates, each metric component is

<div align="center">
<math>
g_{ij} = \left( \partial_i R \right) \left( \partial_j R \right) + \left( \partial_i z \right) \left( \partial_j z \right) + R^2 \left( \partial_i \phi \right) \left( \partial_j \phi \right) .
</math>
</div>

The orthogonality condition is that the metric be diagonal.  That is, each off-diagonal component of the metric must equal zero.

<div align="center">
<math>
\left( \partial_i R \right) \left( \partial_j R \right) + \left( \partial_i z \right) \left( \partial_j z \right) + R^2 \left( \partial_i \phi \right) \left( \partial_j \phi \right) = 0_{ij}, \ \ \ \ \ i \neq j .
</math>
</div>

If the coordinates should be orthogonal, then scale factors can also be defined in terms of the diagonal components of the metric such that <math>g_{ii} = {h_i}^2</math>.  This results in the following definition of the scale factors if the orthogonality condition is met.

<div align="center">
<math>
h_i = \left[ \left( \partial_i R \right)^2 + \left( \partial_i z \right)^2 + R^2 \left( \partial_i \phi \right)^2 \right]^{1/2} .
</math>
</div>

==Christoffel Symbols==
For an orthogonal coordinate system, the Christoffel symbols can be written (depending on which, if any, indices are repeated) as

<table align="center" border="1" cellpadding="5">
<tr>
  <th align="center" colspan="3">
<font color="darkblue">
Christoffel Symbols<br />
</font>
(for orthogonal coordinate systems)
  </th>
</tr>
<tr>
  <td align="center">
<math>
\Gamma^k_{ij} = 0
</math>
  </td>
  <td align="center">
for
  </td>
  <td align="center">
<math>
i \ne j \ne k
</math>
  </td>
</tr>
<tr>
  <td align="center">
<math>
\Gamma^k_{ii} = - \frac{h_i}{h_k} \frac{\partial_k h_i}{h_k}
</math>
  </td>
  <td align="center">
for
  </td>
  <td align="center">
<math>
i \ne k
</math>
  </td>
</tr>
<tr>
  <td align="center">
<math>
\Gamma^i_{ij} = \Gamma^i_{ji} = \frac{\partial_j h_i}{h_i}
</math>
  </td>
  <td align="center" colspan="2">
otherwise<br/>
<font color="blue" size="-2">
(Mistake in original expression corrected <br/>06/10/2010, [[User_talk:Jaycall#Killing_Vector_Approach|as discussed]].)
</font>
  </td>
</tr>
<tr>
  <td align="center">
<math>
\Gamma^i_{ii} = \frac{\partial_i h_i}{h_i}
</math>
  </td>
  <td align="center" colspan="2">
otherwise
  </td>
</tr>
</table>

==Identification of Conserved Quantity==

If, through some miracle, the right-hand side (RHS) of the equation of motion should equal zero, then the combination of variables inside the parentheses must by definition be a conserved quantity.  But good luck guessing a coordinate system in which this miracle will occur!  

===Introduce Characteristic Vector===
The trick is to introduce an arbitrary vector field, which will be used to build a location-dependent weighted linear combination of the equations of motion.  Since the equations of motion associated with each of the coordinates collectively form a single vector equation, we will just form an inner product of our characteristic vector with each side of the vector equation of motion.  The result is

<span id="CV.01"><table align="right" border="1" cellpadding="10" width="10%">
<tr><th><font color="darkblue">CV.01</font></th></tr>
</table></span>
<div align="center">
<math>
C_i \frac{d}{dt} \left( m \ {h_i}^2 \dot{\lambda}_i \right) = m \ {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k C_i - m \ C_i \ \partial_i \Phi
</math>
</div>

Next, we bring <math>C_i</math> inside the total time-derivative on the left-hand side (LHS).  This produces an additional term, which we promptly move over to the RHS and include as part of what is commonly referred to as the ''source'' (since it's the source of any change in the quantity in parentheses).

<div align="center">
<math>
\frac{d}{dt} \left( m \ {h_i}^2 \dot{\lambda}_i C_i \right) = m \ {h_i}^2 \dot{\lambda}_i \dot{C}_i + m \ {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k C_i - m \ C_i \ \partial_i \Phi
</math>
</div>

By doing this, we have formed a new ''conservative'' quantity; that is, a new quantity in the parentheses which will be conserved if and only if the source is zero.  The utility is in the fact that we should be able to force the source associated with this new conservative quantity to go to zero by choosing the right characteristic vector.  Once we find the right characteristic vector, we'll be able to use it directly to build the conserved quantity <math>m \ {h_i}^2 \dot{\lambda}_i C_i</math>.

===Strategic Choice of Characteristic Vector===
So we're looking for a vector <math>\vec{C}</math> such that

<div align="center">
<math>
{h_i}^2 \dot{\lambda}_i \dot{C}_i + {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k C_i - C_i \ \partial_i \Phi = 0 .
</math>
</div>

The good news is that this is a single scalar equation constraining the three independent components of <math>\vec{C}</math>, so several families of solutions should exist.  Consequently, we can afford to be choosey!  

In the problem we're exploring using [[User:Tohline/Appendix/Ramblings/T3Integrals|T3 coordinates]], we already know a conserved quantity associated with the azimuthal coordinate &#8212; angular momentum.  We're looking for an additional (independent) conserved quantity associated with the two radial coordinates.  So it makes sense to look for one of the solutions that has <math>C_3 = 0</math>.  Furthermore, we'd like to avoid dealing with the unknown potential function (which varies only with <math>\lambda_1</math>) as much as possible, so as long as we're being choosey, let's look for a solution that has <math>C_1 = 0</math>.  We now have three conditions on three components.  The condition constraining <math>C_2</math> is

<div align="center">
<math>
{h_2}^2 \dot{\lambda}_2 \dot{C}_2 + {h_k}^2 \Gamma^k_{2j} \dot{\lambda}_j \dot{\lambda}_k C_2 = 0 .
</math>
</div>

This can be rewritten more simply as

<div align="center">
<math>
\dot{C}_2 = - \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) C_2 .
</math>
</div>

Perhaps the best approach to solving this condition is separation of variables.

<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="right">
<math>
\int \frac{dC_2}{C_2} 
</math>
  </td>
  <td align="center">
<math>
= 
</math>
  </td>
  <td align="left">
<math>
- \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow 
</math>
  </td>
  <td align="right">
<math>
\ln C_2 
</math>
  </td>
  <td align="center">
<math>
= 
</math>
  </td>
  <td align="left">
<math>
- \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow 
</math>
  </td>
  <td align="right">
<math>
C_2
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\exp \left\{ - \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt \right\} .
</math>
  </td>
</tr>

</table>

The final necessary step will be to write the quantity in parentheses as the exact derivative of some other quantity, call it <math>\Xi</math>.  The long-sought conserved quantity will then be <math> m {h_2}^2 \dot{\lambda}_2 \exp \left( - \Xi \right) </math> , where

<div align="center">
<math>
\Xi \equiv \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt .
</math>
</div>


====Question from Joel====
Is the following logic correct?  

Suppose we examine the term in which <math>j=2</math> and <math>k=1</math>.  The integral becomes,

<div align="center">
<math>
\Xi\biggr|_{j=2,k=1} = \int \left( \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_1}{\dot{\lambda}_2} \right) dt = \int \left( \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{22}  \right) d\lambda_1
</math>
</div>

Given that the scale factors and the Christoffel symbols are expressible entirely in terms of <math>\lambda_1</math> and <math>\lambda_2</math>, one could imagine a situation &#8212; for example in the quadratic case of <math>q^2=2</math> &#8212; in which this integral could be completed analytically.  (I presume that wherever <math>\lambda_2</math> appears inside this integral, it can be treated as a constant because the two coordinates are independent of one another.)

====Response from Jay====

Yes, I believe you've worked this term out correctly, but as you can see in what follows, it cancels out with the other cross term.  Writing out each of the terms in <math>\Xi</math> gives,

<div align="center">
<math>
\frac{d\Xi}{dt} = \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{21} \frac{\dot{\lambda}_1 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{23} \frac{\dot{\lambda}_3 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \Gamma^2_{21} \frac{\dot{\lambda}_1 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \Gamma^2_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \Gamma^2_{23} \frac{\dot{\lambda}_3 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \Gamma^3_{21} \frac{\dot{\lambda}_1 \dot{\lambda}_3}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \Gamma^3_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_3}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \Gamma^3_{23} \frac{\dot{\lambda}_3 \dot{\lambda}_3}{\dot{\lambda}_2} .
</math>
</div>

Plugging in values for the Christoffel symbols leads to the expression

<div align="center">
<math>
\frac{d\Xi}{dt} = \frac{{h_1}^2}{{h_2}^2} \frac{\partial_2 h_1}{h_1} \frac{\dot{\lambda}_1 \dot{\lambda}_1}{\dot{\lambda}_2} - \frac{{h_1}^2}{{h_2}^2} \frac{h_2}{h_1} \frac{\partial_1 h_2}{h_1} \frac{\dot{\lambda}_2 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \frac{\partial_1 h_2}{h_2} \frac{\dot{\lambda}_1 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \frac{\partial_2 h_2}{h_2} \frac{\dot{\lambda}_2 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \frac{\partial_2 h_3}{h_3} \frac{\dot{\lambda}_3 \dot{\lambda}_3}{\dot{\lambda}_2} .
</math>
</div>

And limiting our interest to motion within the meridional plane (setting <math>\dot{\lambda}_3 = 0</math>) and simplifying

<span id="CV.02"><table align="right" border="1" cellpadding="10" width="10%">
<tr><th><font color="darkblue">CV.02</font></th></tr>
</table></span>
<table align="center" border="1" cellpadding="5">
<tr>
  <td align="right">
<math>
\frac{d\Xi}{dt} \equiv - \frac{d\ln C_2}{dt}
</math>
  </td>
  <td align="center">
<math>
= 
</math>
  </td>
  <td align="left">
<math>
\frac{h_1}{h_2} \frac{\partial_2 h_1}{h_2} \frac{{\dot{\lambda}_1}^2}{\dot{\lambda}_2} - \cancel{\frac{\partial_1 h_2}{h_2} \dot{\lambda}_1} + \cancel{\frac{\partial_1 h_2}{h_2} \dot{\lambda}_1} + \frac{\partial_2 h_2}{h_2} \dot{\lambda}_2 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
= 
</math>
  </td>
  <td align="left">
<math>
\biggl[ \biggl(\frac{h_1 \dot{\lambda}_1}{h_2 \dot{\lambda}_2}\biggr)^2 \frac{\partial \ln h_1}{\partial\ln\lambda_2}  + \frac{\partial \ln h_2}{\partial \ln\lambda_2}\biggr] \frac{d \ln{\lambda}_2}{dt} 
</math>
  </td>
</tr>

<tr>
  <td align="left" colspan="3">
<font color="red"><b>NOTE:</b> Sign error fixed on 15 July 2010</font>.  Specifically, <math>d\Xi/dt \equiv d\ln C_2/dt</math> changed to <math>d\Xi/dt \equiv - d\ln C_2/dt</math>.
  </td>
</tr>
</table>

=Application to T3 Coordinates=

On a [[User:Tohline/Appendix/Ramblings/T3CharacteristicVector|separate page]], we have applied this "characteristic vector" approach specifically to T3 coordinates.

=Related Talk Session=
==Joel's First Talk Message to Jay==
Jay:  In the email message that you sent to me today, you indicated that you had added the 3rd component to the dx/dt equation. I don't see that modification in the current T3Coordinates page. --[[User:Tohline|Tohline]] 15:05, 29 May 2010 (MDT)

:Joel:  It's at the bottom of the section entitled "Time-Derivative of Position and Velocity Vectors".  Under the history tab, can you see the edits that I made?  It was the very first edit and was fairly minor.  I only added a couple terms. --[[User:Jaycall|Jaycall]] 20:17, 29 May 2010 (MDT)

::The mediaWiki "talk" page recommends that when you reply to a query, you should not start a new subsection but, rather, just indent (using one or more colons) immediately following the query.  Also, it recommends that you "sign" each "talk" message.  You do this by typing 2 dashes followed by 4 tildes!  I'm going to edit your "reply" (and this additional one from me) to put it in this recommended format.  But you have to "sign" your own remark. --[[User:Tohline|Tohline]] 16:23, 29 May 2010 (MDT)

::By the way, I ''can'' see the edits that you made to add the 3rd component to the dx/dt equation.  I did not spot the changes earlier, but via the history "diff" function, I'm able to see the edits clearly. --[[User:Tohline|Tohline]] 16:41, 29 May 2010 (MDT)

==Killing Vector Approach==
Thanks for writing out all the terms in the expression for <math>d\Xi/dt</math> in your "[[User:Jaycall/KillingVectorApproach|Killing Vector Approach]]" discussion.  On the one hand, I'm happy that the term I mentioned cancels with another one because that means we have fewer terms to integrate.  On the other hand, it was one term that I thought we might actually be able to manipulate analytically.  Of the two terms that remain, one is relatively simple &#8212; the one that is proportional to <math>\dot{\lambda}_2</math> &#8212; but the first term is a bear!  That will take some thinking! --[[User:Tohline|Tohline]] 10:50, 30 May 2010 (MDT)

:I know it seems like the wrong direction, but unless we can express this summary integral entirely in terms of <math>\lambda_1</math>, <math>\lambda_2</math>, <math>\dot{\lambda}_1</math> and <math>\dot{\lambda}_2</math>, I think that in order to make progress we're going to have to write things out in terms of <math>\varpi</math>, <math>z</math>, <math>\dot{\varpi}</math> and <math>\dot{z}</math>.  Of course, that will make the expression much messier, but I think it will be very difficult to recognize this quantity as an exact derivative unless everything's in terms of the same coordinates. --[[User:Jaycall|Jaycall]] 12:55, 30 May 2010 (MDT)

::I agree.  It is for similar reasons that I am planning on focusing on the case where <math>q^2=2</math>.  At least in this special case we can invert the coordinate expressions to give <math>\varpi</math> and <math>z</math> in terms of the <math>\lambda</math>s.  Incidentally, I have added a link to your "Killing Vector" page as well as a link to this associated "Talk" page on the [[User:Tohline/Appendix/Ramblings|page where I itemize my various "ramblings"]]. --[[User:Tohline|Tohline]] 16:16, 30 May 2010 (MDT)

:::Good point.  That would be a good reason to focus on one specific case.  I haven't looked closely at the inverted coordinate transformation for the <math>q^2=2</math> case for several days.  Since the relation is quadratic, is it obvious which root should be used?

::The proper physical root was obvious to me when I performed the coordinate inversion in the case of [[User:Tohline/Appendix/Ramblings/T1Coordinates#Coordinate_Inversion|T1 Coordinates]], so I presume it will be obvious for the T3 Coordinate system.  But the inversion needs to be redone for the case of T3 Coordinates; do you want to take care of this and type up the result?  In the (quadratic) case of <math>q^2=2</math>, it may be as simple as replacing <math>\chi_2</math> with <math>1/\lambda_2</math>. --[[User:Tohline|Tohline]] 11:21, 31 May 2010 (MDT)

:::Sure. --[[User:Jaycall|Jaycall]] 13:34, 31 May 2010 (MDT)

'''<font color="red">Mistake?</font>''' Please check the next to last row in the [[User:Jaycall/KillingVectorApproach#Christoffel_Symbols|tabulated expressions for Christoffel symbols]]; shouldn't the <math>i</math> and <math>j</math> indexes be swapped? --[[User:Tohline|Tohline]] 09:39, 6 June 2010 (MDT)

:You are correct.  Good eye.--[[User:Jaycall|Jaycall]] 09:42, 9 June 2010 (MDT)

::Actually, it wasn't my good eye.  I discovered this from my inelegant brute-force derivations.  In particular, when I was trying to derive the same form of the [[User:Tohline/Appendix/Ramblings/T3CharacteristicVector#Two_Views_of_Equation_of_Motion|equation of motion from two different points of view]], I had to plug in some of your Christoffel symbol expressions.  At first I could not get my EOM to match yours; it was while trying to understand this mismatch that I realized that your ''general'' expression for <math>\Gamma^i_{ij}</math> did not match the ''specific'' expression for <math>\Gamma^2_{21}</math> that you typed on the "T3 Coordinates" page.  In my effort to get the EOMs to match, I decided that your ''general'' expression was the incorrect one. --[[User:Tohline|Tohline]] 14:58, 9 June 2010 (MDT)

:::For the record, on 10 June 2010, I modified the [[User:Jaycall/KillingVectorApproach#Christoffel_Symbols|table that contains the general Christoffel symbol expressions]] in order to correct the order of these indexes.--[[User:Tohline|Tohline]] 18:55, 11 June 2010 (MDT)

==Logarithmic Derivatives of T3 Scale Factors==

Check out my new subsection (under T3 Coordinates) entitled [[User:Tohline/Appendix/Ramblings/T3Integrals#Logarithmic_Derivatives_of_Scale_Factors|Logarithmic Derivatives of Scale Factors]].  First, see if you agree that there is a mistake (typo) in one of your tables of partial derivatives.  Specifically, I think that the correct expression is:
<div align="center">
<math>
\frac{\partial z}{\partial\lambda_2} = - (q^2-1)\frac{\varpi^2 z \ell^2}{\lambda_2} .
</math>
</div>
Second, see if you agree with my derived expressions for <math>\partial\ln h_i/\partial\ln\lambda_j</math>. --[[User:Tohline|Tohline]] 22:28, 31 May 2010 (MDT)

:You're certainly right about the sign error.  I have corrected the expression in the relevant table.  --[[User:Jaycall|Jaycall]] 14:28, 1 June 2010 (MDT)

:I have not yet been able to confirm your expressions for logarithmic derivatives of the scale factors.  I'm not sure what approach you took in deriving them, but since I had already calculated partials of the scale factors (<math>\partial_i h_j</math> and so forth), I derived a little trick to help simplify the work.  Do you agree that

<div align="center">
<math>
\frac{\partial \ln h_i}{\partial \ln \lambda_j} \equiv \frac{\lambda_j}{h_i} \partial_j h_i ?
</math>
</div>

::Yes, this is precisely how I defined the logarithmic derivatives.  And I used your tables of partial derivatives (e.g., <math>\partial\varpi/\partial\lambda_j</math> and <math>\partial z/\partial\lambda_j</math>) to obtain <math>\partial_i h_j</math> and so forth. --[[User:Tohline|Tohline]] 17:30, 1 June 2010 (MDT)

'''<font color="red">Mistake!</font>''' Jay:  I have found a mistake in my original derivation of the logarithmic derivatives of the <math>h_1</math> scale factor.  Perhaps now they will match yours. Check it out and let me know! --[[User:Tohline|Tohline]] 12:41, 4 June 2010 (MDT)

::Joel:  I am finally able to confirm your derivation of the logarithmic derivatives of the scale factors.  I came at them from a different angle with a fresh head, and I am now confident that what you have is correct.  Next I will look into the general case.  --[[User:Jaycall|Jaycall]] 14:29, 10 July 2010 (MDT)

==Two Views of Equation of Motion==

Jay: I have created a [[User:Tohline/Appendix/Ramblings/T3CharacteristicVector|new wiki page]] to explicitly analyze how your Characteristic Vector formalism applies to the T3 Coordinate system.  On this page is a subsection entitled "[[User:Tohline/Appendix/Ramblings/T3CharacteristicVector#Two_Views_of_Equation_of_Motion|Two Views of Equation of Motion]]" in which I have shown that the <math>2^\mathrm{nd}</math> component of the equation of motion can be written in the following form,
<div align="center">
<math>\frac{d(h_2^2 \dot{\lambda}_2)}{dt} = \biggl( h_1 \frac{\partial h_1}{\partial \lambda_2} \biggr) \dot{\lambda}_1^2 + \biggl( h_2 \frac{\partial h_2}{\partial \lambda_2} \biggr) \dot{\lambda}_2^2</math> ,
</div>
assuming that <math>\partial\Phi/\partial\lambda_2 = 0</math>.  More to the point, I have shown that I can derive this form of this component of the equation of motion using either (A) your Christoffel symbol formalism, or (B) the more classical formalism that I have been using, which I obtained from Appendix 1.B of Binney &amp; Tremaine (1987) and Morse &amp; Feshbach (1953).  

Now, I shouldn't be shocked that both derivations lead to the same form of the equation &#8212; after all, we're doing physics aren't we?  But here is what surprised me:  The Christoffel symbol formalism produces this form of the equation without ever assuming ''a particular coordinate system'' whereas, in order to derive the equation using the more classical formalism, I plugged in some relations between the partial derivatives of the scale factors that ''I thought'' were specific to T3 coordinates.  This leads me to ask, "How many of the T3 Coordinate relations [[User:Tohline/Appendix/Ramblings/T3Integrals|derived earlier]] are generic relations that apply to any orthogonal, axisymmetric coordinate system, and how many are unique to the T3 Coordinate System?"  I'm particularly interested in knowing how generalizable the relations are that appear inside boxes labeled "T3 Coordinates" under the subsections I have entitled, "[[User:Tohline/Appendix/Ramblings/T3Integrals#Derived_Identity_for_T3_Coordinates|Derived Identity for T3 Corodinates]]" and "[[User:Tohline/Appendix/Ramblings/T3Integrals#Logarithmic_Derivatives_of_Scale_Factors|Logarithmic Derivatives of Scale Factors]]". Can you answer this? --[[User:Tohline|Tohline]] 10:46, 7 June 2010 (MDT)

:Joel:  Yes, this is the power of using a covariant formulation.  The only assumptions we made were that the coordinates were orthogonal and axisymmetric.  I believe that the relation under "[[User:Tohline/Appendix/Ramblings/T3Integrals#Derived_Identity_for_T3_Coordinates|Derived Identity for T3 Coordinates]]" would apply to any coordinate system meeting these criteria.  I will double check, though.  (How did you determine that the position vector could be written <math>\vec{x} = \hat{e}_1 (h_1 \lambda_1) + \hat{e}_2 (h_2 \lambda_2)</math>?)

:The logarithmic derivatives of the scale factors cannot be general in their current form due to the appearance of factors of <math>q</math>, but I suspect they can be generalized and molded into a form that closely resembles what you have.  I will work on this, too. --[[User:Jaycall|Jaycall]] 10:39, 9 June 2010 (MDT)

::Jay: You asked how I determined the expression for the [[User:Tohline/Appendix/Ramblings/T3Integrals#Definition|position vector]].  I don't have my detailed notes with me right now, but I can outline the steps.  In Morse &amp; Feshbach's classical presentation of orthogonal curvilinear coordinates, they define a set of so-called ''direction cosines,'' <math>\gamma_{ij}</math>.  (I don't have the generalized definition with me at this time, but they are expressible in terms of derivatives of the scale factors with respect to the various coordinates; I'm certain ''direction cosines'' are definable in terms of Christoffel symbols.) Each of the Cartesian unit vectors (<math>\hat{\imath},\hat{\jmath},\hat{k}</math>) can be written in terms of products of these direction cosines and the three unit vectors of your chosen curvilinear coordinate system (<math>\hat{e}_1, \hat{e}_2, \hat{e}_3</math>).  According to Morse &amp; Feshbach, for example, <math>\hat{\imath} = (\hat{e}_1\gamma_{11} + \hat{e}_2\gamma_{21} + \hat{e}_3\gamma_{31})</math>.  So, all I did to obtain an expression for the position vector in T3 coordinates was define <math>\vec{x} = (\hat\imath x + \hat\jmath y + \hat{k}z)</math>, then replace each Cartesian unit vector by its expression in terms of products of the T3-coordinate unit vectors &amp; the appropriate direct cosine, then gather together and simplify terms.  Does this make sense? --[[User:Tohline|Tohline]] 14:59, 10 June 2010 (MDT)

:::Yes, it does.  That's very helpful. --[[User:Jaycall|Jaycall]] 09:49, 11 June 2010 (MDT)

::::Jay:  For the record, I have put together a wiki page that discusses the definition and utility of [[User:Tohline/Appendix/Ramblings/DirectionCosines|Direction Cosines]].--[[User:Tohline|Tohline]] 18:47, 3 July 2010 (MDT)

Joel, I finally figured out what was going on this afternoon with the strange conserved quantities we were getting, like <math>{h_2}^4 \dot{\lambda_2}^2</math> and, consequently, <math>\lambda_1 \lambda_2</math>.  There ''is'' a mistake in the derivation stemming all the way back from Eq. CV.02 on my page on the [[User:Jaycall/KillingVectorApproach|characteristic vector approach]].  My derivation of the brute force condition on <math>C_2</math> was correct, but you can see from above that <math>\Xi</math> does not equal <math>\ln C_2</math>, it equals <math>- \ln C_2</math>.  Consequently, the conclusion is no longer that <math>C_2</math> must equal <math>{h_2}^2 \dot{\lambda_2}</math>, but rather <math>\left( {h_2}^2 \dot{\lambda_2} \right)^{-1}</math>.  So in fact you ''did'' find a conserved quantity--<math>1</math>!   : )  I guess it was a tautology after all; your intuition was right.  --[[User:Jaycall|Jaycall]] 18:50, 14 July 2010 (MDT)

:Jay, very good catch!  Looking back through the history of that wiki page, it is clear that I am the one who introduced this sign error.  I mistakingly set <math>\Xi = \ln C_2</math> in equation [[User:Jaycall/KillingVectorApproach#CV.02|CV.02]] at the same time that I "Added a couple of right-justified equation numbers."  I'll march through the relevent wiki pages today and attempt to correct this mistake wherever it appears. --[[User:Tohline|Tohline]] 08:30, 15 July 2010 (MDT)

::Joel, I'm sorry yesterday's excitement didn't pan out.  It would have been really neat.  Well, at least I learned something.  On my drive home yesterday I was thinking about what it meant for there to be a constraint that involved only the coordinates (and not their rates of change).  We were uncomfortable with this idea, but now I understand why in general no conserved quantity can take this form.  If there were a constraint on only the coordinates, then this would limit a particle to a curve or a surface...a subspace of physical space (as opposed to a subspace of phase space).  Yet it would definitely be possible to get a particle to leave the subspace by giving it an initial velocity in any direction normal to the subspace.  Any conserved quantity must definitely involve the coordinate velocities (unless you're artificially restricting the motion). --[[User:Jaycall|Jaycall]] 14:15, 15 July 2010 (MDT)

:::Good observation. --[[User:Tohline|Tohline]] 16:28, 15 July 2010 (MDT)

==Special Quadratic Case==
Jay: On a [[User:Tohline/Appendix/Ramblings/T3Integrals/QuadraticCase#Special_Case_.28Quadratic.29|new accompanying wiki page]] I have begun to investigate how our system behaves in the special case of <math>q^2 = 2</math>.  I have added a pointer/link to this new page &#8212; and a separate pointer/link to the parallel investigation that you have told me you are conducting &#8212; on the page containing a table of contents to my [[User:Tohline/Appendix/Ramblings|various research ramblings]].  I have not yet actually gone to the page that details your analysis because, at least for a while, I prefer for our trains of thought to remain independent of one another.  However, in an effort to simplify the comparison that we ultimately will need to make between our independent analyses, I strongly recommend that you write your expressions for the scale factors, etc. ''for the special case of <math>q^2=2</math>'' in terms of a variable, <math>\Lambda</math>, defined such that,
<div align="center">
<math>
\Lambda \equiv \biggl[1 + \biggl(\frac{2\lambda_1}{\lambda_2}\biggr)^2\biggr]^{1/2} = \cosh(2\Zeta) .
</math>
</div>
By the way, if you don't agree that <math>\sqrt{1+(2\lambda_1/\lambda_2)^2}</math> equals <math>\cosh(2\Zeta)</math>, please let me know!--[[User:Tohline|Tohline]] 19:49, 11 June 2010 (MDT)

:: Joel: I have checked and agree that in the <math>q^2=2</math> case only, the above relation among <math>\Lambda</math>, <math>Z</math> and the T3 coordinates is correct. --[[User:Jaycall|Jaycall]] 14:58, 10 July 2010 (MDT)


{{ SGFfooter }}