User:Tohline/SSC/Structure/Other Analytic Ramblings

From VistrailsWiki
Jump to navigation Jump to search

Ramblings Regarding the Stability of Other Analytically Definable, Spherical Equilibrium Models

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

The material presented in this chapter was originally developed as a subsection of a chapter that discusses "Other" Analytic Equilibrium Models."



Work-in-progress.png

Material that appears after this point in our presentation is under development and therefore
may contain incorrect mathematical equations and/or physical misinterpretations.
|   Go Home   |


Generic Setup

Dividing the above, 2nd-order ODE through by the quantity, <math>~[R^2 (P_0/P_c)]</math>, gives,

<math> \frac{d^2x}{dr_0^2} + \biggl[\frac{4}{r_0} - \frac{1}{R}\biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \biggl(\frac{P_c}{P_0}\biggr)\biggr] \frac{dx}{dr_0} - \biggl[\frac{1}{R}\biggl(\frac{P_c}{P_0}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \frac{\alpha}{r_0} \biggr] x = - \frac{1}{R^2}\biggl(\frac{P_c}{P_0}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl[ \biggl( \frac{\tau_\mathrm{SSC}^2 \omega^2}{\gamma_g} \biggr) \biggr] x \, , </math>


which matches Prasad's (1949) equation (1), namely,

<math>~x^{' '} + \biggl[\frac{4}{r_0} - \frac{\mu(r_0) }{r_0}\biggr] x^{'} - \biggl[ \frac{\alpha \mu(r_0)}{r_0^2} \biggr] x</math>

<math>~=</math>

<math>~- \biggl(\frac{P_c}{P_0}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl[ \frac{n^2\rho_c}{\gamma_g P_c} \biggr] x \, ,</math>

where, primes indicate differentiation with respect to <math>~r_0</math>, and,

<math>~\mu(r_0) \equiv \frac{r_0}{R} \biggl(\frac{P_c}{P_0}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \, .</math>

(Note that Prasad's equation has the awkward units of inverse length-squared.) Regrouping terms in Prasad's governing equation, multiplying through by <math>~R^2</math> (to make the equation dimensionless), and now letting primes denote differentiation with respect to the dimensionless radial coordinate, <math>~\chi_0</math>, we quite generally can write the linear adiabatic wave equation as,

<math>~- \biggl(\frac{P_c}{P_0}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr) \sigma^2 x </math>

<math>~=</math>

<math>~\biggl[x^{' '} + \frac{4 x^'}{\chi_0}\biggr] - \frac{\mu(\chi_0)}{\chi_0} \biggl[ x^{'} + \frac{\alpha x}{\chi_0} \biggr]</math>

 

<math>~=</math>

<math>~\frac{1}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) - \frac{\mu(\chi_0)}{\chi_0} \biggl[\frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \biggr] \, . </math>

Defining,

<math>~A</math>

<math>~\equiv</math>

<math>~\biggl(\frac{P_0}{P_c}\biggr)\biggl(\frac{\rho_c}{\rho_0}\biggr) \, ,</math>

<math>~B</math>

<math>~\equiv</math>

<math>~\frac{A\mu(\chi_0)}{\chi_0} = \biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \, ,</math>

the governing equation becomes,

<math>~\sigma^2 x </math>

<math>~=</math>

<math>~\frac{B}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) -\frac{A}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) </math>

 

<math>~=</math>

<math>~B \biggl[ \frac{\alpha x}{\chi_0} + x^'\biggr] - A \biggl[ \frac{4x^'}{\chi_0}+ x^{' '} \biggr] \, . </math>

Notice that, because,

<math>~g_0 = - \frac{1}{\rho_0} ~\frac{dP_0}{dr_0} \, ,</math>

at every radial location throughout the configuration, it must also be true that, for any equilibrium configuration,

<math>~B</math>

<math>~=</math>

<math>~- \biggl(\frac{\rho_0}{\rho_c}\biggr)^{-1} \frac{d(P_0/P_c)}{d\chi_0}</math>

<math>~\Rightarrow ~~~~ \frac{B}{A}</math>

<math>~=</math>

<math>~- \frac{d}{d\chi_0}\biggl[\ln\biggl(\frac{P_0}{P_c}\biggr)\biggr] \, .</math>

The following table shows that this relationship holds for a collection of analytically described equilibrium structures.

Table 1a: Properties of Analytically Defined Equilibrium Structures
Model <math>~\frac{\rho_0}{\rho_c}</math> <math>~\frac{P_0}{P_c}</math> <math>~\frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr)</math>
Uniform-density <math>~1</math> <math>~1 - \chi_0^2</math> <math>~-2\chi_0</math>
Linear <math>~1-\chi_0</math> <math>~\tfrac{1}{5} (5 -24 \chi_0^2 + 28\chi_0^3 - 9\chi_0^4)</math> <math>~\tfrac{1}{5}[- 48\chi_0 + 84\chi_0^2 - 36\chi_0^3]</math>
Parabolic <math>~1-\chi_0^2</math> <math>~\tfrac{1}{2} (2 - 5\chi_0^2 + 4\chi_0^4 - \chi_0^6)</math> <math>~- 5\chi_0 + 8\chi_0^3 - 3\chi_0^5</math>
<math>~n=1</math> Polytrope <math>~\frac{\sin(\pi\chi_0)}{\pi\chi_0}</math> <math>~\biggl[\frac{\sin(\pi\chi_0)}{\pi\chi_0}\biggr]^2</math> <math>~\frac{2\sin(\pi\chi_0)}{(\pi^2\chi_0^3)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr]</math>


Table 1b: Properties of Analytically Defined Equilibrium Structures
Model <math>~\frac{\rho_0}{\rho_c}</math> <math>~B\equiv \frac{g_0}{g_\mathrm{SSC}}</math> <math>~A \equiv \biggl(\frac{P_0}{P_c}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr)^{-1}</math> <math>~\biggl(\frac{P_0}{P_c}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr)^{-2}</math>
Uniform-density <math>~1</math> <math>~2\chi_0</math> <math>~1 - \chi_0^2</math> <math>~1 - \chi_0^2</math>
Linear <math>~1-\chi_0</math> <math>~\tfrac{48}{5}(\chi_0 - \tfrac{3}{4}\chi_0^2)</math> <math>~\tfrac{1}{5} (1-\chi_0) (5 + 10\chi_0 - 9\chi_0^2)</math> <math>~(1 + 2\chi_0 - \tfrac{9}{5}\chi_0^2)</math>
Parabolic <math>~1-\chi_0^2</math> <math>~5\chi_0 - 3\chi_0^3</math> <math>~\tfrac{1}{2} (1-\chi_0^2) (2 - \chi_0^2)</math> <math>~(1 - \tfrac{1}{2} \chi_0^2)</math>
<math>~n=1</math> Polytrope <math>~\frac{\sin(\pi\chi_0)}{\pi\chi_0}</math> <math>~\frac{2}{\pi\chi_0^2}\biggl[ \sin(\pi\chi_0) - \pi\chi_0 \cos(\pi\chi_0) \biggr]</math> <math>~\frac{\sin(\pi\chi_0)}{\pi\chi_0}</math> <math>~1</math>

Leaning on this new expression for the ratio, <math>~B/A</math>, let's play with the form of the governing equation.

<math>~- \sigma^2 x </math>

<math>~=</math>

<math>~A \biggl\{ \frac{1}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) + \frac{d}{d\chi_0}\biggl[\ln\biggl(\frac{P_0}{P_c}\biggr)\biggr] \cdot \frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \biggr\} </math>

 

<math>~=</math>

<math>~A\biggl(\frac{P_0}{P_c}\biggr)^{-1} \biggl\{ \biggl(\frac{P_0}{P_c}\biggr)\frac{1}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) + \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr) \cdot \frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \biggr\} </math>

Polytropic Configurations

Let's compare this presentation of the LAWE to the form of the LAWE that has been derived specifically for polytropic equilibrium configurations, namely,

<math>~\frac{d^2x}{d\xi^2} + \biggl[\frac{4 - (n+1)V(\xi)}{\xi} \biggr] \frac{dx}{d\xi} + \biggl[\omega^2 \biggl(\frac{a_n^2 \rho_c }{\gamma_g P_c} \biggr) \frac{\theta_c}{\theta} - \biggl(3-\frac{4}{\gamma_g}\biggr) \cdot \frac{(n+1)V(x)}{\xi^2} \biggr] x </math>

<math>~=</math>

<math>0 \, ,</math>

where,

<math>~V(\xi)</math>

<math>~\equiv</math>

<math>~- \frac{\xi}{(\theta/\theta_c)} \frac{d (\theta/\theta_c)}{d\xi} = \frac{g_0}{a_n}\biggl(\frac{a_n^2\rho_0}{P_0}\biggr)\frac{\xi}{(n+1)} \, .</math>

[Note that <math>~\theta_c = 1</math> and, therefore for all practical purposes, it can be dropped. This notation was introduced in our separate discussion of the polytropic LAWE in order to make it clear how our derivations have overlapped earlier published work.] Regrouping terms, we have,

<math>~-\biggl[\omega^2 \biggl(\frac{a_n^2 \rho_c }{\gamma_g P_c} \biggr) \frac{\theta_c}{\theta}\biggr]x</math>

<math>~=</math>

<math>~\frac{d^2x}{d\xi^2} + \biggl[\frac{4 - (n+1)V(\xi)}{\xi} \biggr] \frac{dx}{d\xi} - \biggl[ \frac{\alpha (n+1)V(x)}{\xi^2} \biggr] x</math>

 

<math>~=</math>

<math>~\biggl[\frac{d^2x}{d\xi^2} + \biggl(\frac{4}{\xi}\biggr)\frac{dx}{d\xi} \biggr] -\biggl[\frac{(n+1)V(\xi)}{\xi} \biggr] \biggl[\frac{dx}{d\xi} + \frac{\alpha x}{\xi} \biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{\xi^4}\frac{d}{d\xi}\biggl(\xi^4 \frac{dx}{d\xi}\biggr) -\biggl[\frac{(n+1)V(\xi)}{\xi} \biggr] \biggl[\frac{1}{\xi^\alpha}\frac{d}{d\xi} \biggl(\xi^\alpha x\biggr)\biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{\xi^4}\frac{d}{d\xi}\biggl(\xi^4 \frac{dx}{d\xi}\biggr) +\biggl[(n+1)\frac{d\ln(\theta/\theta_c)}{d\xi} \biggr] \biggl[\frac{1}{\xi^\alpha}\frac{d}{d\xi} \biggl(\xi^\alpha x\biggr)\biggr] \, .</math>

Next we note that, written in terms of the traditional polytropic radial coordinate, <math>~\xi</math>, the fractional radius,

<math>~\chi_0 \equiv \frac{r_0}{R} = \frac{\xi}{\xi_1} = \frac{a_n \xi}{R} \, .</math>

Hence, multiplying the polytropic LAWE through by the quantity, <math>~(R/a_n)^2</math>, gives,

<math>~\frac{1}{\chi_0^4}\frac{d}{d\chi_0}\biggl(\chi_0^4 \frac{dx}{d\chi_0}\biggr) +\biggl[(n+1)\frac{d\ln(\theta/\theta_c)}{d\chi_0} \biggr] \biggl[\frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0} \biggl(\chi_0^\alpha x\biggr)\biggr] </math>

<math>~=</math>

<math>~-\biggl[\omega^2 \biggl(\frac{R^2 \rho_c }{\gamma_g P_c} \biggr) \frac{\theta_c}{\theta}\biggr]x</math>

 

<math>~=</math>

<math>~-\biggl(\frac{\theta_c}{\theta}\biggr)\sigma^2 x \, .</math>

Finally, noting that, for polytropic configurations,

<math>~\frac{\theta}{\theta_c}</math>

<math>~=</math>

<math>~\biggl( \frac{P_0}{P_c} \biggr)\biggl( \frac{\rho_0}{\rho_c} \biggr)^{-1} = \biggl( \frac{P_0}{P_c} \biggr)^{1/(n+1)} \, , </math>

we can rewrite the polytropic LAWE in the form,

<math>~\frac{1}{\chi_0^4}\frac{d}{d\chi_0}\biggl(\chi_0^4 \frac{dx}{d\chi_0}\biggr) +\biggl[\frac{d\ln(P_0/P_c)}{d\chi_0} \biggr] \biggl[\frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0} \biggl(\chi_0^\alpha x\biggr)\biggr] </math>

<math>~=</math>

<math>~-\biggl( \frac{P_0}{P_c} \biggr)^{-1}\biggl( \frac{\rho_0}{\rho_c} \biggr)\sigma^2 x \, ,</math>

which precisely matches the general expression for the LAWE presented at the end of our generic setup, directly above.

This seems to be a particularly insightful way to write the LAWE, as the only structural functions that appear explicitly are <math>~P_0(\chi_0)</math> and <math>~\rho_0(\chi_0)</math>. It appears as though the eigenfunctions that describe adiabatic radial pulsations do not explicitly depend a priori on the radial dependence of the equilibrium gravitational acceleration.

Examine Structural Pressure-Density Relation

Derivation

One striking property exhibited by the example configurations tabulated above is the structural relationship between the chosen function, <math>~\rho_0(\chi_0)</math>, and the corresponding radial pressure distribution, <math>~P_0(\chi_0)</math>, that is dictated by,

Hydrostatic Balance

<math>\frac{1}{\rho}\frac{dP}{dr} =- \frac{d\Phi}{dr} = - \frac{GM_r}{r^2} </math> ,

As has been detailed in the last column of Table 1b, in all four cases, the ratio, <math>~(P_0/P_c)(\rho_0/\rho_c)^{-2}</math>, is an analytically prescribed polynomial expression. That is, the pressure is "evenly divisible" by the square of the density. Let's examine how broadly reliable this behavior is. [Note that, for simplicity in typing, hereafter throughout this subsection we will drop the subscript zero and, rather than <math>~\chi_0</math>, we will use <math>~z \equiv r/R</math> to denote the dimensionless radial coordinate.]


Assume a mass-distribution given by the general quadratic function,

<math>~\frac{\rho}{\rho_c}</math>

<math>~=</math>

<math>~1 - a z - b z^2 \, ,</math>

where both coefficients, <math>~a</math> and <math>~b</math>, are positive.


ASIDE: Surface Location

The surface of the configuration will be defined by the radial location, <math>~z_s</math>, at which the density first goes to zero. If <math>~b = 0</math>, then the surface will be located at <math>~z_s = a^{-1}</math>; and if <math>~a = 0</math>, it will be located at <math>~z_s = b^{-1/2}</math>. More generally, however, the roots of the quadratic equation that results from setting <math>~\rho/\rho_c</math> to zero are,

<math>~z_\pm</math>

<math>~=</math>

<math>~- \frac{a}{2b}\biggl[1 \mp \biggl(1+\frac{4b}{a^2}\biggr)^{1/2} \biggr] \, .</math>

Because only the <math>~z_+</math> solution provides positive roots, we conclude that, when both <math>~a</math> and <math>~b</math> are nonzero, the radial coordinate of the surface is,

<math>~z_s = z_+</math>

<math>~=</math>

<math>~\frac{a}{2b}\biggl[\biggl(1+\frac{4b}{a^2}\biggr)^{1/2} - 1\biggr] \, .</math>

We acknowledge, as well, that the density profile can now be written in terms of these roots; specifically,

<math>~\frac{\rho}{\rho_c}</math>

<math>~=</math>

<math>~b(z_+ - z)(z - z_-) \, .</math>

In the discussion, below, it may be advantageous to adopt the following notation:

<math>~\ell^2 \equiv \frac{4b}{a^2} ~~~~~\Rightarrow ~~~~~ \ell = \frac{2b^{1/2}}{a} \, ,</math>

in which case,

<math>~az_s = \frac{2}{\ell^2}\biggl[\biggl(1+\ell^2\biggr)^{1/2} - 1\biggr]</math>         and         <math>~b^{1/2}z_s = \frac{1}{\ell}\biggl[\biggl(1+\ell^2\biggr)^{1/2} - 1\biggr] \, .</math>


This specified density profile implies a mass distribution,

<math>~M_r</math>

<math>~=</math>

<math>~4\pi R^3 \rho_c \int_0^z (1 - a z - b z^2)z^2 dz</math>

 

<math>~=</math>

<math>~4\pi R^3 \rho_c \biggl( \frac{z^3}{3} - \frac{a z^4}{4} - \frac{b z^5}{5} \biggr) \, .</math>

The hydrostatic balance condition therefore implies,

<math>~\frac{1}{R\rho_c} \biggl( \frac{\rho}{\rho_c} \biggr)^{-1} \frac{dP}{dz}</math>

<math>~=</math>

<math>~-\frac{G}{R^2 z^2} \biggl[ 4\pi R^3 \rho_c \biggl( \frac{z^3}{3} - \frac{a z^4}{4} - \frac{b z^5}{5} \biggr) \biggr]</math>

<math>~\Rightarrow ~~~~\biggl[ \frac{1}{4\pi G R^2 \rho_c^2 } \biggr] \frac{dP}{dz}</math>

<math>~=</math>

<math>~-\biggl( \frac{z}{3} - \frac{a z^2}{4} - \frac{b z^3}{5} \biggr)\biggl( \frac{\rho}{\rho_c} \biggr) </math>

 

<math>~=</math>

<math>~ - \biggl( \frac{z}{3} - \frac{a z^2}{4} - \frac{b z^3}{5} \biggr) + a \biggl( \frac{z^2}{3} - \frac{a z^3}{4} - \frac{b z^4}{5} \biggr) + b \biggl( \frac{z^3}{3} - \frac{a z^4}{4} - \frac{b z^5}{5} \biggr) </math>

 

<math>~=</math>

<math>~ -\biggl(\frac{1}{3}\biggr)z + z^2 \biggl( \frac{a}{4} + \frac{a}{3} \biggr) + z^3 \biggl( \frac{b}{5} - \frac{a^2}{4} + \frac{b}{3} \biggr) + z^4 \biggl( - \frac{ab}{5} - \frac{ab}{4} \biggr) - z^5 \biggl(\frac{b^2}{5} \biggr) </math>

 

<math>~=</math>

<math>~ -\biggl(\frac{1}{3}\biggr)z + z^2 \biggl( \frac{7a}{12} \biggr) + z^3 \biggl( \frac{8b}{15} - \frac{a^2}{4}\biggr) - z^4 \biggl( \frac{9ab}{20} \biggr) - z^5 \biggl(\frac{b^2}{5} \biggr) </math>

<math>\Rightarrow ~~~~\biggl( \frac{P}{P_n} \biggr)</math>

<math>~=</math>

<math>~\int \biggl[ -\biggl(\frac{1}{3}\biggr)z + z^2 \biggl( \frac{7a}{12} \biggr) + z^3 \biggl( \frac{8b}{15} - \frac{a^2}{4}\biggr) - z^4 \biggl( \frac{9ab}{20} \biggr) - z^5 \biggl(\frac{b^2}{5} \biggr) \biggr] dz </math>

 

<math>~=</math>

<math>~ \biggl[ -\biggl(\frac{1}{6}\biggr)z^2 + z^3 \biggl( \frac{7a}{36} \biggr) + z^4 \biggl( \frac{2b}{15} - \frac{a^2}{16}\biggr) - z^5 \biggl( \frac{9ab}{100} \biggr) - z^6 \biggl(\frac{b^2}{30} \biggr) \biggr] + C \, . </math>

where, <math>~C</math>, is an integration constant, and,

<math>~P_n \equiv 4\pi G \rho_c^2 R^2 \, .</math>

The integration constant — which also proves to be the normalized central pressure — is determined by ensuring that the pressure goes to zero at the surface of the configuration, <math>~z_s</math>. That is,

<math>~C = \frac{P_c}{P_n}</math>

<math>~=</math>

<math>~ - \biggl[ -\biggl(\frac{1}{6}\biggr)z_s^2 + z_s^3 \biggl( \frac{7a}{36} \biggr) + z_s^4 \biggl( \frac{2b}{15} - \frac{a^2}{16}\biggr) - z_s^5 \biggl( \frac{9ab}{100} \biggr) - z_s^6 \biggl(\frac{b^2}{30} \biggr) \biggr] \, . </math>

Hence, the pressure profile is,

<math>\frac{P}{P_n} </math>

<math>~=</math>

<math>~ \biggl(\frac{1}{6}\biggr)(z_s^2 - z^2) - \biggl( \frac{7a}{36} \biggr)(z_s^3 - z^3) - \biggl( \frac{2b}{15} - \frac{a^2}{16}\biggr) (z_s^4 - z^4) + \biggl( \frac{9ab}{100} \biggr)(z_s^5 - z^5) + \biggl(\frac{b^2}{30} \biggr)(z_s^6 - z^6) \, . </math>

Let's check this general expression against the specific cases described above.

Example1

First, let's set <math>~b=0</math>, but leave <math>~a</math> general. As described in the above ASIDE, this means that <math>~z_s=a^{-1}</math>. So the pressure distribution is,

<math>\frac{P}{P_n} </math>

<math>~=</math>

<math>~ \biggl(\frac{1}{6}\biggr)(z_s^2 - z^2) - \biggl( \frac{7a}{36} \biggr)(z_s^3 - z^3) + \biggl( \frac{a^2}{16}\biggr) (z_s^4 - z^4) </math>

 

<math>~=</math>

<math>~ \biggl(\frac{1}{6}\biggr)(a^{-2} - z^2) - \biggl( \frac{7a}{36} \biggr)(a^{-3} - z^3) + \biggl( \frac{a^2}{16}\biggr) (a^{-4} - z^4) </math>

 

<math>~=</math>

<math>~a^{-2} \biggl\{ \frac{1}{6}[1 - (az)^2] - \frac{7}{36} [1 - (az)^3] + \frac{1}{16} [1 - (az)^4] \biggr\} \, . </math>

And from the expression for the integration constant, we have,

<math>~\frac{P_c}{P_n}</math>

<math>~=</math>

<math>~\biggl[ \biggl(\frac{1}{6}\biggr)a^{-2} - a^{-3} \biggl( \frac{7a}{36} \biggr) + a^{-4} \biggl( \frac{a^2}{16}\biggr) \biggr] </math>

 

<math>~=</math>

<math>~a^{-2} \biggl[\frac{1}{6} - \frac{7}{36} + \frac{1}{16} \biggr] </math>

 

<math>~=</math>

<math>~\frac{5}{2^4\cdot 3^2 a^2} \, . </math>

Hence, dividing one expression by the other, we obtain,

<math>\frac{P}{P_c} </math>

<math>~=</math>

<math>~\frac{1}{5}\biggl\{ 24[1 - (az)^2] - 28 [1 - (az)^3] + 9 [1 - (az)^4] \biggr\} \, . </math>

Let's check to see if this "general linear" pressure distribution is evenly divisible by the square of the density distribution which, in this case, is,

<math>~\frac{\rho}{\rho_c}</math>

<math>~=</math>

<math>~1 - a z \, .</math>

Strategically rewriting the expression for the pressure distribution gives,

<math>\frac{P}{P_c} </math>

<math>~=</math>

<math>~\frac{1}{5}\biggl\{ 24[1 - (az)][1 + (az)] - 28 [1 - (az)][1 + (az) + (az)^2] + 9 [1 - (az)][1 + (az)][1 + (az)^2] \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{5} \biggl(\frac{\rho}{\rho_c}\biggr) \biggl\{ 24[1 + (az)] - 28 [1 + (az) + (az)^2] + 9 [1 + (az)][1 + (az)^2] \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{5} \biggl(\frac{\rho}{\rho_c}\biggr) \biggl\{5+ 24(az) - 28 [(az) + (az)^2] + 9 [(az) + (az)^2 + (az)^3] \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{5} \biggl(\frac{\rho}{\rho_c}\biggr) \biggl\{5+ 5(az) - 19 (az)^2 + 9 (az)^3 \biggr\} \, . </math>

And, as luck would have it, the expression inside the curly braces can be "divided evenly" by the quantity, <math>~[1-(az)]</math>, one more time. Specifically, the expression becomes,

<math>\frac{P}{P_c} </math>

<math>~=</math>

<math>~\frac{1}{5} \biggl(\frac{\rho}{\rho_c}\biggr)^2 [5+ 10(az) - 9 (az)^2 ] \, . </math>


Example2

First, let's set <math>~a=0</math>, but leave <math>~b</math> general. As described in the above ASIDE, this means that <math>~z_s=b^{-1/2}</math>. So the pressure distribution is,

<math>\frac{P}{P_n} </math>

<math>~=</math>

<math>~ \biggl(\frac{1}{6}\biggr)(z_s^2 - z^2) - \biggl( \frac{2b}{15} \biggr) (z_s^4 - z^4) + \biggl(\frac{b^2}{30} \biggr)(z_s^6 - z^6) </math>

 

<math>~=</math>

<math>~ \biggl(\frac{1}{6}\biggr)(b^{-1} - z^2) - \biggl( \frac{2b}{15} \biggr) (b^{-2} - z^4) + \biggl(\frac{b^2}{30} \biggr)(b^{-3} - z^6) </math>

 

<math>~=</math>

<math>~\frac{1}{b}\biggl\{ \biggl(\frac{1}{6}\biggr)[1 - bz^2] - \biggl( \frac{2}{15} \biggr) [1 - (bz^2)^2] + \biggl(\frac{1}{30} \biggr)[1 - (bz^2)^3] \biggr\} \, . </math>

And from the expression for the integration constant, we have,

<math>~\frac{P_c}{P_n}</math>

<math>~=</math>

<math>~ \biggl[ \biggl(\frac{1}{6}\biggr)z_s^2 - z_s^4 \biggl( \frac{2b}{15}\biggr) + z_s^6 \biggl(\frac{b^2}{30} \biggr) \biggr] </math>

 

<math>~=</math>

<math>~ \frac{1}{b} \biggl[ \frac{1}{6} - \frac{2}{15} + \frac{1}{30} \biggr] </math>

 

<math>~=</math>

<math>~ \frac{1}{3 \cdot 5b} \, . </math>

Hence, dividing one expression by the other, we obtain,

<math>\frac{P}{P_c} </math>

<math>~=</math>

<math>~\frac{1}{2}\biggl\{ 5 [1 - bz^2] - 4 [1 - (bz^2)^2] + [1 - (bz^2)^3] \biggr\} \, . </math>

Let's check to see if this "general parabolic" pressure distribution is evenly divisible by the square of the density distribution which, in this case, is,

<math>~\frac{\rho}{\rho_c}</math>

<math>~=</math>

<math>~1 - b z^2 \, .</math>

Strategically rewriting the expression for the pressure distribution gives,

<math>\frac{P}{P_c} </math>

<math>~=</math>

<math>~\frac{1}{2}\biggl\{ 5[1 - (bz^2)] - 4[1 - (bz^2)][1 + (bz^2)] + [1 - (bz^2)][1 + (bz^2) + (bz^2)^2] \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{2} \biggl( \frac{\rho}{\rho_c}\biggr) \biggl\{ 5 - 4 [1 + (bz^2)] + [1 + (bz^2) + (bz^2)^2] \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{2} \biggl( \frac{\rho}{\rho_c}\biggr) \biggl[ 2 - 3 (bz^2) + (bz^2)^2 \biggr] </math>

Again, as luck would have it, the expression inside the square brackets can be "divided evenly" by the quantity, <math>~[1-(bz^2)]</math>, one more time. Specifically, the expression becomes,

<math>\frac{P}{P_c} </math>

<math>~=</math>

<math>~\frac{1}{2} \biggl(\frac{\rho}{\rho_c}\biggr)^2 [2-(bz^2) ] \, . </math>

Example3

In the most general quadratic case, we should rewrite the pressure distribution as,

<math>\frac{P}{P_n} </math>

<math>~=</math>

<math>~ \frac{1}{2^4\cdot 3^2\cdot 5^2}\biggl\{ 2^3\cdot 3\cdot 5^2 (z_s^2 - z^2) - 2^2\cdot 5^2 \cdot 7 a(z_s^3 - z^3) - [2^5\cdot 3\cdot 5 b - 3^2\cdot 5^2 a^2] (z_s^4 - z^4) + 2^2\cdot 3^4 ab (z_s^5 - z^5) + 2^3\cdot 3\cdot 5b^2 (z_s^6 - z^6) \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{z_s^2}{2^4\cdot 3^2\cdot 5^2}\biggl\{ 600 (1 - \zeta^2) - 700 (a z_s) (1 - \zeta^3) - [480 (b z_s^2) - 225 (a z_s)^2] (1 - \zeta^4) + 324 (az_s)(b z_s^2) (1 - \zeta^5) + 120(b z_s^2)^2 (1 - \zeta^6) \biggr\} \, , </math>

where, <math>~\zeta \equiv z/z_s</math>. Similarly, let's rewrite the integration constant as,

<math>~\frac{P_c}{P_n}</math>

<math>~=</math>

<math>~ \frac{z_s^2}{2^4\cdot 3^2\cdot 5^2}\biggl\{ 600 - 700 (a z_s) - [480 (b z_s^2) - 225 (a z_s)^2] + 324 (az_s)(b z_s^2) + 120(b z_s^2)^2 \biggr\} \, . </math>

So the pressure can be written as,

<math>\biggl[ \frac{2^4\cdot 3^2\cdot 5^2}{z_s^2} \biggr] \biggl[ \frac{P}{P_n} \biggr]</math>

<math>~=</math>

<math>~ \biggl[ \frac{2^4\cdot 3^2\cdot 5^2}{z_s^2} \biggr] \biggl[ \frac{P_c}{P_n}\biggr] - 600 \zeta^2 + 700 (a z_s) \zeta^3 + [480 (b z_s^2) - 225 (a z_s)^2] \zeta^4 - 324 (az_s)(b z_s^2) \zeta^5 - 120(b z_s^2)^2 \zeta^6 \, . </math>

The question that remains to be answered is: Is this expression for the pressure distribution "evenly divisible" by the square (or even the first power) of the normalized density distribution which, as defined above for the general quadratic case, is,

<math>\frac{\rho}{\rho_c} = 1 - az - bz^2 = 1 - (az_s)\zeta - (bz_s^2)\zeta^2 \, .</math>

In attempting to answer this question, it may prove advantageous to refer back to the above ASIDE discussion of the roots of this quadratic function and, in particular, that,

<math>~az_s = \frac{2}{\ell^2}\biggl[\biggl(1+\ell^2\biggr)^{1/2} - 1\biggr]</math>         and         <math>~b^{1/2}z_s = \frac{1}{\ell}\biggl[\biggl(1+\ell^2\biggr)^{1/2} - 1\biggr] \, ,</math>

where, <math>~\ell^2 \equiv 4b/a^2</math>.


Uniform Density

In the case of a uniform-density configuration, the governing equation is,

<math>~\sigma^2 x </math>

<math>~=</math>

<math>~2\chi_0 \biggl[ \frac{\alpha x}{\chi_0} + x^'\biggr] - (1-\chi_0^2) \biggl[ \frac{4x^'}{\chi_0}+ x^{' '} \biggr] \, , </math>

where,

<math>~\sigma^2 \equiv \frac{\tau_\mathrm{SSC}^2 \omega^2}{\gamma_g} = \frac{6}{\gamma_g}\biggl[\frac{\omega^2}{4\pi G\bar\rho}\biggr] \, .</math>

The following individual mode analyses should be compared with the results found in our discussion of Sterne's general solution.

Mode 0

Try an eigenfunction of the form,

<math>x = a_0\, ,</math>

in which case,

<math>x^' = x^{' '} = 0 \, .</math>

In order for this to be a solution, we must have,

<math>~\sigma^2 a_0 </math>

<math>~=</math>

<math>~2\chi_0 \biggl[ \frac{\alpha a_0}{\chi_0} \biggr] </math>

<math>~\Rightarrow ~~~~ \frac{6}{\gamma_g}\biggl[\frac{\omega^2}{4\pi G\bar\rho}\biggr]</math>

<math>~=</math>

<math>~2\alpha = 2\biggl(3 - \frac{4}{\gamma_g}\biggr) </math>

<math>~\Rightarrow ~~~~ \frac{\omega^2}{4\pi G\bar\rho}</math>

<math>~=</math>

<math>~\gamma_g - \frac{4}{3}\, . </math>


Mode 2

Try an eigenfunction of the form,

<math>x = a_0 + a_2\chi_0^2 \, ,</math>

in which case,

<math>~x^' = 2 a_2\chi_0 </math>        and         <math>~x^{' '} = 2 a_2 \, . </math>

In order for this to be a solution, we must have,

<math>~\sigma^2 \biggl(a_0 + a_2\chi_0^2\biggr) </math>

<math>~=</math>

<math>~2 \biggl[ \alpha (a_0 + a_2\chi_0^2) + \chi_0 (2 a_2\chi_0 )\biggr] - (1-\chi_0^2) \biggl[ \frac{4(2 a_2\chi_0 )}{\chi_0}+ 2a_2 \biggr] </math>

 

<math>~=</math>

<math>~2\alpha a_0 + \chi_0^2[2a_2 (2+\alpha)] - 10a_2(1-\chi_0^2) </math>

<math>~\Rightarrow ~~~~ \sigma^2 a_0 - 2\alpha a_0 + 10a_2</math>

<math>~=</math>

<math>~\chi_0^2 [-\sigma^2 + 2 (2+\alpha) + 10 ]a_2 \, . </math>

Given that the coefficients on both sides of this expression must independently be zero, we have:

<math>~\sigma^2 </math>

<math>~=</math>

<math>~2 (2+\alpha) + 10</math>

<math>~\Rightarrow ~~~~ \frac{\omega^2}{4\pi G\bar\rho} </math>

<math>~=</math>

<math>~\frac{\gamma_g}{6}\biggl[14 +2\biggl(3-\frac{4}{\gamma_g} \biggr)\biggr]</math>

 

<math>~=</math>

<math>~\frac{\gamma_g}{6}\biggl[20 -\frac{8}{\gamma_g}\biggr] = \frac{1}{3}\biggl(10\gamma_g - 4\biggr) \, ,</math>

and

<math>~\frac{a_2}{a_0}</math>

<math>~=</math>

<math>~\frac{1}{10} \biggl[ 2\alpha - \sigma^2 \biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{10} \biggl\{ 2\alpha - [14+2\alpha) ] \biggr\} = - \frac{7}{5} \, .</math>


Parabolic Density Distribution

In the case of a parabolic density distribution, the governing equation is,

<math>~\sigma^2 x </math>

<math>~=</math>

<math>~(5\chi_0 - 3\chi_0^3)\biggl[ \frac{\alpha x}{\chi_0} + x^'\biggr] - \tfrac{1}{2} (1-\chi_0^2) (2 - \chi_0^2) \biggl[ \frac{4x^'}{\chi_0}+ x^{' '} \biggr] \, , </math>

where,

<math>~\sigma^2 \equiv \frac{\tau_\mathrm{SSC}^2 \omega^2}{\gamma_g} = \frac{15}{\gamma_g}\biggl[\frac{\omega^2}{4\pi G\bar\rho}\biggr] \, .</math>


First Trial

Try an eigenfunction of the form,

<math>x = (2-\chi_0^2)^{-1} (a + b\chi_0^2 + c\chi_0^4) \, ,</math>

in which case,

<math>~x^' </math>

<math>~=</math>

<math>~ (2-\chi_0^2)^{-1} (2 b\chi_0 + 4c\chi_0^3) + 2\chi_0 (2-\chi_0^2)^{-2} (a + b\chi_0^2 + c\chi_0^4) </math>

 

<math>~=</math>

<math>~ (2-\chi_0^2)^{-2}[(2-\chi_0^2) (2 b\chi_0 + 4c\chi_0^3) + 2\chi_0 (a + b\chi_0^2 + c\chi_0^4) ]</math>

 

<math>~=</math>

<math>~ 2\chi_0 (2-\chi_0^2)^{-2}[(2-\chi_0^2) (b + 2c\chi_0^2) + (a + b\chi_0^2 + c\chi_0^4) ]</math>

 

<math>~=</math>

<math>~ 2\chi_0 (2-\chi_0^2)^{-2}[ (2b+4c\chi_0^2-b\chi_0^2 -2c\chi_0^4) + (a + b\chi_0^2 + c\chi_0^4) ]</math>

 

<math>~=</math>

<math>~ 2\chi_0 (2-\chi_0^2)^{-2}[ (a + 2b)+ 4c\chi_0^2 -c \chi_0^4 ] \, ;</math>

and,

<math>~x^{' '} </math>

<math>~=</math>

<math>~(2-\chi_0^2)^{-1} (2 b + 12 c\chi_0^2) + 2\chi_0 (2-\chi_0^2)^{-2} (2 b\chi_0 + 4c\chi_0^3) + 2(2-\chi_0^2)^{-2} (a + b\chi_0^2 + c\chi_0^4) </math>

 

 

<math>~+ 2\chi_0 (a + b\chi_0^2 + c\chi_0^4) [-2(2-\chi_0^2)^{-3}(-2\chi_0)] + 2\chi_0 (2-\chi_0^2)^{-2} (2 b\chi_0 + 4c\chi_0^3) </math>

 

<math>~=</math>

<math>~(2-\chi_0^2)^{-3} [ 2(4-4\chi_0^2 + \chi_0^4) (b + 6 c\chi_0^2) + 4\chi_0^2 (2-\chi_0^2) (b + 2c\chi_0^2) + 2(2-\chi_0^2) (a + b\chi_0^2 + c\chi_0^4) </math>

 

 

<math>~+ 8\chi_0^2 (a + b\chi_0^2 + c\chi_0^4) + 4\chi_0^2 (2-\chi_0^2) (b + 2c\chi_0^2) ]</math>

 

<math>~=</math>

<math>~2(2-\chi_0^2)^{-3} \{ ( 4b+24c\chi_0^2 - 4b\chi_0^2-24c\chi_0^4 + b\chi_0^4 + 6c\chi_0^6 ) +(8b \chi_0^2+16c\chi_0^4 -4b\chi_0^4 - 8c\chi_0^6) </math>

 

 

<math>~ +(2a+2b\chi_0^2 + 2c\chi_0^4 -a\chi_0^2 - b\chi_0^4 - c\chi_0^6) + (4a\chi_0^2 + 4b\chi_0^4 + 4c\chi_0^6) \}</math>

 

<math>~=</math>

<math>~2(2-\chi_0^2)^{-3} \{ (4b+2a) + \chi_0^2(24c-4b+8b+2b-a+4a) + \chi_0^4(-24c+b+16c-4b+2c-b+4b) + \chi_0^6(6c-8c-c+4c) \} </math>

 

<math>~=</math>

<math>~2 (2-\chi_0^2)^{-3} \{(2a+4b) + \chi_0^2(3a + 6b + 24c) + \chi_0^4(-6c) + \chi_0^6(c)\} </math>


In order for this to be a solution, we must have,

<math>~\sigma^2 x </math>

<math>~=</math>

<math>~(5 - 3\chi_0^2)\biggl[ \alpha x + \chi_0 x^'\biggr] - (1-\chi_0^2) (2 - \chi_0^2) \biggl[ \frac{2x^'}{\chi_0}+ \frac{x^{' '}}{2} \biggr] </math>

<math>~\Rightarrow~~~~\sigma^2 (2-\chi_0^2)^{-1} (a + b\chi_0^2 + c\chi_0^4) </math>

<math>~=</math>

<math>~ (5 - 3\chi_0^2)\biggl[ \alpha (2-\chi_0^2)^{-1} (a + b\chi_0^2 + c\chi_0^4) + 2\chi_0^2 (2-\chi_0^2)^{-2}[ (a + 2b)+ 4c\chi_0^2 -c \chi_0^4 ]\biggr] </math>

 

 

<math>~ - (1-\chi_0^2) \biggl[ 4 (2-\chi_0^2)^{-1}[ (a + 2b)+ 4c\chi_0^2 -c \chi_0^4 ]+ (2-\chi_0^2)^{-2} \{(2a+4b) + \chi_0^2(3a + 6b + 24c) -6c \chi_0^4 + c\chi_0^6\} \biggr] \, . </math>

Multiplying through by <math>~(2-\chi_0^2)^2</math> gives,

<math>~\sigma^2 (2-\chi_0^2) (a + b\chi_0^2 + c\chi_0^4) </math>

<math>~=</math>

<math>~ (5 - 3\chi_0^2)\biggl\{ \alpha (2-\chi_0^2) (a + b\chi_0^2 + c\chi_0^4) + 2\chi_0^2 [ (a + 2b)+ 4c\chi_0^2 -c \chi_0^4 ]\biggr\} </math>

 

 

<math>~ - (1-\chi_0^2) \biggl\{ 4 (2-\chi_0^2)[ (a + 2b)+ 4c\chi_0^2 -c \chi_0^4 ]+ [(2a+4b) + \chi_0^2(3a + 6b + 24c) -6c \chi_0^4 + c\chi_0^6] \biggr\} </math>

<math>~\Rightarrow~~~~\sigma^2 [2a + \chi_0^2(2b-a) + \chi_0^4(2c -b) - c\chi_0^6]</math>

<math>~=</math>

<math>~ (5 - 3\chi_0^2)\biggl\{ \alpha [2a + \chi_0^2(2b-a) + \chi_0^4(2c -b) - c\chi_0^6]\biggr\} </math>

 

 

<math>~ + (5 - 3\chi_0^2)\biggl\{ (2a + 4b)\chi_0^2 + 8c\chi_0^4 - 2c \chi_0^6 \biggr\} </math>

 

 

<math>~ - (1-\chi_0^2) \biggl\{ [ 8(a + 2b)+ 32c\chi_0^2 - 8c \chi_0^4 ] - [ 4\chi_0^2(a + 2b)+ 16c\chi_0^4 -4c \chi_0^6 ]\biggr\} </math>

 

 

<math>~ - (1-\chi_0^2) \biggl\{ (2a+4b) + \chi_0^2(3a + 6b + 24c) -6c \chi_0^4 + c\chi_0^6 \biggr\} </math>

 

<math>~=</math>

<math>~ (5 - 3\chi_0^2)\biggl\{ 2a\alpha + \chi_0^2 [(2b-a)\alpha + (2a + 4b) ] + \chi_0^4[ (2c -b)\alpha + 8c] - (2+\alpha)c\chi_0^6 \biggr\} </math>

 

 

<math>~ - (1-\chi_0^2) \biggl\{ (10a + 20b)+ \chi_0^2[56c - a - 2b] +\chi_0^4 [-30c ]+ 5c \chi_0^6 \biggr\} </math>

So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~2a\sigma^2 - 10a\alpha +10a +20b =~a(2\sigma^2 - 10\alpha+10) + 20b</math>

<math>~\chi_0^2</math>   :   <math>~(2b-a)\sigma^2 -5[(2b-a)\alpha + (2a + 4b) ] +6a\alpha+[56c - a - 2b] - (10a + 20b)</math>

<math>=~a[-\sigma^2 -5(-\alpha + 2) +6\alpha-1-10] + b[2\sigma^2-10\alpha-20-2-20 ] + c[ 56]</math>

<math>=~a[-\sigma^2 + 11\alpha - 21] + b[2\sigma^2-10\alpha-42 ] + 56 c</math>

<math>~\chi_0^4</math>   :   <math>~(2c-b)\sigma^2 - 5[ (2c -b)\alpha + 8c] - 30c + 3[(2b-a)\alpha + (2a + 4b)] - [56c - a - 2b ] </math>

<math>=~a[-3\alpha+7] + b[-\sigma^2+5\alpha+6\alpha+12+2] + c[2\sigma^2-10\alpha-40-30-56]</math>

<math>=~a[7-3\alpha] + b[11\alpha-\sigma^2+14] + c[2\sigma^2-10\alpha-126]</math>

<math>~\chi_0^6</math>   :   <math>~- c\sigma^2 +(10+5\alpha)c +3[ (2c -b)\alpha + 8c] +5c +30c</math>

<math>=~b[-3\alpha] + c[-\sigma^2 + 10 + 5\alpha+6\alpha+24+35]</math>

<math>=~b[-3\alpha] + c[11\alpha -\sigma^2 + 69]</math>

<math>~\chi_0^8</math>   :  

<math>~- 3(2+\alpha)c -5c =~c[-11-3\alpha]</math>

Independent Investigation of Parabolic Distribution

In the specific case of a parabolic density distribution, the leading factor on the LHS is,

<math>~\frac{1}{A_\mathrm{parab}} \equiv \biggl(\frac{P_c}{P_0}\biggr)\biggl(\frac{\rho_0}{\rho_c}\biggr) </math>

<math>~=</math>

<math>~\frac{(1-\chi_0^2)}{(1-\chi_0^2)^2 (1-\tfrac{1}{2}\chi_0^2)} = \frac{1}{(1-\chi_0^2) (1-\tfrac{1}{2}\chi_0^2)} = \frac{2}{2 - 3\chi_0^2 + \chi_0^4} \, ,</math>

and the function appearing on the RHS is,

<math>~\mu(\chi_0)</math>

<math>~=</math>

<math>~\frac{\chi_0^2(1-\chi_0^2)(5-3\chi_0^2)}{(1-\chi_0^2)^2 (1-\tfrac{1}{2}\chi_0^2)} = \frac{\chi_0^2 (5-3\chi_0^2)}{A_\mathrm{parab}} \, .</math>

Multiplying the linear adiabatic wave equation through by <math>~A_\mathrm{parab}</math>, gives,

<math>~- \sigma^2 x </math>

<math>~=</math>

<math>~\frac{A_\mathrm{parab}}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) - \frac{B_\mathrm{parab} }{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \, , </math>

where,

<math> ~B_\mathrm{parab} \equiv \chi_0 (5-3\chi_0^2) \, . </math>

Now we note that,

<math>~\frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' \biggr]</math>

<math>~=</math>

<math>~ \frac{A_\mathrm{parab}}{\chi_0^4} \frac{d}{d\chi_0}\biggl[ \chi_0^4 x^' \biggr] + \chi_0^4 x^' \frac{d}{d\chi_0}\biggl[ \frac{A_\mathrm{parab}}{\chi_0^4} \biggr] </math>

<math>~\Rightarrow ~~~~\frac{A_\mathrm{parab}}{\chi_0^4} \frac{d}{d\chi_0}\biggl[ \chi_0^4 x^' \biggr] </math>

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' \biggr] - \chi_0^4 x^' \frac{d}{d\chi_0}\biggl[ \frac{A_\mathrm{parab}}{\chi_0^4} \biggr] </math>

 

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' \biggr] - \frac{\chi_0^4 x^'}{2} \frac{d}{d\chi_0}\biggl[ 1 - \frac{3}{\chi_0^2} + \frac{2}{\chi_0^4}\biggr] </math>

 

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' \biggr] - \chi_0^4 x^' \biggl[ \frac{3}{\chi_0^3} - \frac{4}{\chi_0^5}\biggr] </math>

 

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' \biggr] + \frac{x^'}{\chi_0} (4 - 3\chi_0^2 ) \, . </math>

Similarly we note that,

<math>~\frac{d}{d\chi_0}\biggl[ B_\mathrm{parab}~x \biggr]</math>

<math>~=</math>

<math>~ \frac{B_\mathrm{parab}}{\chi_0^\alpha} \frac{d}{d\chi_0}\biggl[ \chi_0^\alpha x \biggr] + \chi_0^\alpha x \frac{d}{d\chi_0}\biggl[ \frac{B_\mathrm{parab}}{\chi_0^\alpha} \biggr] </math>

<math>~\Rightarrow ~~~~ \frac{B_\mathrm{parab}}{\chi_0^\alpha} \frac{d}{d\chi_0}\biggl[ \chi_0^\alpha x \biggr] </math>

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ B_\mathrm{parab}~x \biggr] - \chi_0^\alpha x \frac{d}{d\chi_0}\biggl[ \frac{B_\mathrm{parab}}{\chi_0^\alpha} \biggr] </math>

 

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ B_\mathrm{parab}~x \biggr] - \chi_0^\alpha x \frac{d}{d\chi_0}\biggl[ 5\chi_0^{1-\alpha} -3 \chi_0^{3-\alpha}\biggr] </math>

 

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ B_\mathrm{parab}~x \biggr] - \chi_0^\alpha x \biggl[ 5(1-\alpha)\chi_0^{-\alpha} -3 (3-\alpha)\chi_0^{2-\alpha}\biggr] </math>

 

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[ B_\mathrm{parab}~x \biggr] - x \biggl[ 5(1-\alpha) -3 (3-\alpha)\chi_0^{2}\biggr] \, . </math>

Hence, the LAWE can be rewritten as,

<math>~- \sigma^2 x </math>

<math>~=</math>

<math>~\frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' \biggr] + \frac{x^'}{\chi_0} (4 - 3\chi_0^2 ) - \frac{d}{d\chi_0}\biggl[ B_\mathrm{parab}~x \biggr] + x \biggl[ 5(1-\alpha) -3 (3-\alpha)\chi_0^{2}\biggr] \, ; </math>

then multiplying through by <math>~\chi_0</math>, and rearranging terms gives,

<math>~- x \biggl\{ [ 5(1-\alpha)+ \sigma^2]\chi_0 -3 (3-\alpha)\chi_0^{3} \biggr\}</math>

<math>~=</math>

<math>~\chi_0 ~\frac{d}{d\chi_0}\biggl[ A_\mathrm{parab}~x^' - B_\mathrm{parab}~x\biggr] + x^'(4 - 3\chi_0^2 ) \, . </math>

Next, we note that,

<math>~ A_\mathrm{parab}~x^' </math>

<math>~=</math>

<math>~\frac{d}{d\chi_0} \biggl( A_\mathrm{parab}~x \biggr) - x \biggl[\frac{d}{d\chi_0}\biggl(A_\mathrm{parab}\biggr) \biggr]</math>

 

<math>~=</math>

<math>~\frac{d}{d\chi_0} \biggl( A_\mathrm{parab}~x \biggr) - \frac{x}{2} \biggl[\frac{d}{d\chi_0}\biggl(2 - 3\chi_0^2 + \chi_0^4\biggr) \biggr]</math>

 

<math>~=</math>

<math>~\frac{d}{d\chi_0} \biggl( A_\mathrm{parab}~x \biggr) + x (3\chi_0 - 2\chi_0^3 )</math>

<math>~\Rightarrow ~~~~A_\mathrm{parab}~x^' - B_\mathrm{parab}~x</math>

<math>~=</math>

<math>~\frac{d}{d\chi_0} \biggl( A_\mathrm{parab}~x \biggr) + \biggl[(3\chi_0 - 2\chi_0^3 ) - ( 5\chi_0 - 3\chi_0^3 )\biggr] x</math>

 

<math>~=</math>

<math>~\frac{d}{d\chi_0} \biggl( A_\mathrm{parab}~x \biggr) - (2\chi_0 - \chi_0^3 ) x \, .</math>

So, the LAWE becomes,

<math>~- x \biggl\{ [ 5(1-\alpha)+ \sigma^2]\chi_0 -3 (3-\alpha)\chi_0^{3} \biggr\}</math>

<math>~=</math>

<math>~\chi_0 ~\frac{d}{d\chi_0}\biggl[ \frac{d}{d\chi_0} \biggl( A_\mathrm{parab}~x \biggr) - (2\chi_0 - \chi_0^3 ) x\biggr] + x^'(4 - 3\chi_0^2 ) </math>

 

<math>~=</math>

<math>~ \chi_0 \frac{d^2}{d\chi_0^2} \biggl( A_\mathrm{parab}~x \biggr) - \chi_0 ~\frac{d}{d\chi_0}\biggl[ (2\chi_0 - \chi_0^3 ) x\biggr] + x^'(4 - 3\chi_0^2 ) </math>

 

<math>~=</math>

<math>~ \chi_0 \frac{d^2}{d\chi_0^2} \biggl( A_\mathrm{parab}~x \biggr) - \biggl\{ \frac{d}{d\chi_0}\biggl[ (2\chi_0^2 - \chi_0^4)x\biggr] - (2\chi_0 - \chi_0^3)x \biggr\} + \biggl\{ \frac{d}{d\chi_0}\biggl[ (4-3\chi_0^2)x\biggr] + 6\chi_0 x \biggr\} </math>

 

<math>~=</math>

<math>~ \chi_0 \frac{d^2}{d\chi_0^2} \biggl( A_\mathrm{parab}~x \biggr) + \frac{d}{d\chi_0}\biggl[(4-3\chi_0^2)x -(2\chi_0^2 - \chi_0^4)x\biggr] + (2\chi_0 - \chi_0^3)x + 6\chi_0 x </math>

 

<math>~=</math>

<math>~ \chi_0 \frac{d^2}{d\chi_0^2} \biggl( A_\mathrm{parab}~x \biggr) + \frac{d}{d\chi_0}\biggl[(4-5\chi_0^2 + \chi_0^4)x \biggr] + (8\chi_0 - \chi_0^3)x \, . </math>

Moving the last term on the RHS of this expression to the LHS, and factoring the polynomial coefficients of the terms inside of the first and second derivatives gives,

<math>~x \biggl\{ (5\alpha - 13 - \sigma^2)\chi_0 + (10 - 3\alpha)\chi_0^{3} \biggr\}</math>

<math>~=</math>

<math>~ \frac{\chi_0}{2} \frac{d^2}{d\chi_0^2} \biggl[ (1-\chi_0^2)(2-\chi_0^2)x \biggr] + \frac{d}{d\chi_0}\biggl[(1-\chi_0^2)(4-\chi_0^2)x \biggr] \, . </math>



First Trial (Same as Above)

Try an eigenfunction of the form,

<math>x = (2-\chi_0^2)^{-1} (a + b\chi_0^2 + c\chi_0^4) \, ,</math>

in which case,


LHS

<math>~=</math>

<math>~ \chi_0 (2-\chi_0^2)^{-2}\biggl\{ (a + b\chi_0^2 + c\chi_0^4) (2-\chi_0^2) [ (5\alpha - 13 - \sigma^2) + (10 - 3\alpha)\chi_0^{2}] \biggr\} </math>

 

<math>~=</math>

<math>~ \chi_0 (2-\chi_0^2)^{-2} (a + b\chi_0^2 + c\chi_0^4) \biggl\{ (10\alpha - 26 - 2\sigma^2) + (\sigma^2 -11\alpha + 33 )\chi_0^2 + (3\alpha -10)\chi_0^{4} \biggr\} </math>

 

<math>~=</math>

<math>~ \chi_0 (2-\chi_0^2)^{-2}\biggl\{ a(10\alpha - 26 - 2\sigma^2) + \chi_0^2[a(\sigma^2 -11\alpha + 33 ) + b(10\alpha - 26 - 2\sigma^2)] </math>

 

 

<math>~ +\chi_0^{4}[ a (3\alpha -10) + b(\sigma^2 -11\alpha + 33 ) + c(10\alpha - 26 - 2\sigma^2)] + \chi_0^{6}[ b (3\alpha -10) + c(\sigma^2 -11\alpha + 33 )] + c(3\alpha -10)\chi_0^{8} \biggr\} </math>


RHS (1st term)

<math>~=</math>

<math>~ \frac{\chi_0}{2} \frac{d^2}{d\chi_0^2} \biggl[ (1-\chi_0^2)(a + b\chi_0^2 + c\chi_0^4) \biggr] </math>

 

<math>~=</math>

<math>~ \frac{\chi_0}{2} \frac{d^2}{d\chi_0^2} \biggl[ a + (b-a)\chi_0^2 + (c-b)\chi_0^4 - c\chi_0^6\biggr] </math>

 

<math>~=</math>

<math>~ \frac{\chi_0}{2} \frac{d}{d\chi_0} \biggl[ 2(b-a)\chi_0 + 4(c-b)\chi_0^3 - 6c\chi_0^5\biggr] </math>

 

<math>~=</math>

<math>~ \frac{\chi_0}{2} [ 2(b-a) + 12(c-b)\chi_0^2 - 30c\chi_0^4] </math>

 

<math>~=</math>

<math>~ (b-a)\chi_0 + 6(c-b)\chi_0^3 - 15c\chi_0^5 </math>

 

<math>~=</math>

<math>~\chi_0 (2-\chi_0^2)^{-2}\biggl\{ (4-4\chi_0^2 + \chi_0^4) [(b-a) + 6(c-b)\chi_0^2 - 15c\chi_0^4 ] \biggr\} </math>

 

<math>~=</math>

<math>~\chi_0 (2-\chi_0^2)^{-2}\biggl\{ [(4b-4a) + (24c- 24b)\chi_0^2 - 60c\chi_0^4 ] + [(-4b + 4a)\chi_0^2 + (-24c + 24b)\chi_0^4 + 60c\chi_0^6 ] + [(b-a)\chi_0^4 + 6(c-b)\chi_0^6 - 15c\chi_0^8 ] \biggr\} </math>

 

<math>~=</math>

<math>~\chi_0 (2-\chi_0^2)^{-2}\biggl\{ (4b-4a) + [(24c- 24b) + (-4b + 4a)]\chi_0^2 +[ - 60c + (-24c + 24b) + (b-a)]\chi_0^4 + [60c + 6(c-b)]\chi_0^6 - 15c\chi_0^8 \biggr\} </math>

 

<math>~=</math>

<math>~\chi_0 (2-\chi_0^2)^{-2}\biggl\{ (4b-4a) + [ 24c- 28b + 4a]\chi_0^2 +[ - 84c + 25b -a ]\chi_0^4 + [66c -6b ]\chi_0^6 - 15c\chi_0^8 \biggr\} </math>


RHS (2nd term)

<math>~=</math>

<math>~ \frac{d}{d\chi_0}\biggl[(1-\chi_0^2)(4-\chi_0^2)(2-\chi_0^2)^{-1} (a + b\chi_0^2 + c\chi_0^4) \biggr] </math>

 

<math>~=</math>

<math>~ (2-\chi_0^2)^{-1} \frac{d}{d\chi_0}\biggl[(4-5\chi_0^2+\chi_0^4)(a + b\chi_0^2 + c\chi_0^4) \biggr] </math>

 

 

<math>~+ (4-5\chi_0^2+\chi_0^4)(a + b\chi_0^2 + c\chi_0^4) \frac{d}{d\chi_0}\biggl[(2-\chi_0^2)^{-1} \biggr] </math>

 

<math>~=</math>

<math>~ (2-\chi_0^2)^{-1} \frac{d}{d\chi_0}\biggl[4a + \chi_0^2(4b-5a) + \chi_0^4(4c-5b+a) + \chi_0^6(b-5c) + c\chi_0^8 \biggr] </math>

 

 

<math>~+ \biggl[4a + \chi_0^2(4b-5a) + \chi_0^4(4c-5b+a) + \chi_0^6(b-5c) + c\chi_0^8 \biggr] \frac{d}{d\chi_0}\biggl[(2-\chi_0^2)^{-1} \biggr] </math>

 

<math>~=</math>

<math>~ (2-\chi_0^2)^{-2}\biggl\{(2-\chi_0^2) \biggl[\chi_0 (8b-10a) + \chi_0^3 (16c-20b+4a) + \chi_0^5 (6b-30c) + 8c\chi_0^7 \biggr] </math>

 

 

<math>~+2\chi_0 \biggl[4a + \chi_0^2(4b-5a) + \chi_0^4(4c-5b+a) + \chi_0^6(b-5c) + c\chi_0^8 \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \chi_0 (2-\chi_0^2)^{-2}\biggl\{\biggl[(16b-20a) + \chi_0^2 (32c-40b+8a) + \chi_0^4 (12b-60c) + 16c\chi_0^6 \biggr] </math>

 

 

<math>~ - \biggl[\chi_0^2 (8b-10a) + \chi_0^4 (16c-20b+4a) + \chi_0^6 (6b-30c) + 8c\chi_0^8 \biggr] </math>

 

 

<math>~+ \biggl[8a + \chi_0^2 (8b-10a) + \chi_0^4 (8c-10b+2a) + \chi_0^6 (2b-10c) + 2c\chi_0^8 \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \chi_0 (2-\chi_0^2)^{-2}\biggl\{[(16b-20a) + 8a] + \chi_0^2 [(32c-40b+8a) + (10a-8b) + (8b-10a)] </math>

 

 

<math>~ + \chi_0^4 [(12b-60c)+ (20b-4a-16c) + (8c-10b+2a)] + \chi_0^6[16c + (30c-6b) + (2b-10c) ] - 6c\chi_0^8 \biggr\} </math>

 

<math>~=</math>

<math>~ \chi_0 (2-\chi_0^2)^{-2}\biggl\{[16b-12a] + \chi_0^2 [(32c-40b+8a) ] + \chi_0^4 [-2a + 22b -68c] + \chi_0^6[36c -4b ] - 6c\chi_0^8 \biggr\} </math>


So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~a(10\alpha - 26 - 2\sigma^2) - (4b-4a) - [16b-12a] ~=a(10\alpha - 10 - 2\sigma^2) - 20b</math>

<math>~\chi_0^2</math>   :   <math>~[a(\sigma^2 -11\alpha + 33 ) + b(10\alpha - 26 - 2\sigma^2)] - [ 24c- 28b + 4a]- [(32c-40b+8a) ]</math>

<math>= ~[a(\sigma^2 -11\alpha + 21 ) + b(10\alpha + 42 - 2\sigma^2)] - 56c</math>

<math>~\chi_0^4</math>   :   <math>~[ a (3\alpha -10) + b(\sigma^2 -11\alpha + 33 ) + c(10\alpha - 26 - 2\sigma^2)] - [ - 84c + 25b -a ] - [-2a + 22b -68c]</math>

<math>~=~[ a (3\alpha -7) + b(\sigma^2 -11\alpha -14 ) + c(10\alpha + 126 - 2\sigma^2)]</math>

<math>~\chi_0^6</math>   :   <math>~[ b (3\alpha -10) + c(\sigma^2 -11\alpha + 33 )] - [66c -6b ]-[36c -4b ]</math>

<math>~=~[ b (3\alpha ) + c(\sigma^2 -11\alpha - 69 )] </math>

<math>~\chi_0^8</math>   :  

<math>~c(3\alpha -10) + 15c + 6c = c[3\alpha+11]</math>

The expressions for the coefficients presented in this table exactly match the entire set of expressions derived earlier, except for the adopted sign convention — every term in this second derivation has the opposite sign to the corresponding term in the earlier derivation.

Conjecture

Returning to the generic formulation derived earlier, we have,

<math>~- \sigma^2 \biggl(\frac{\rho_0}{\rho_c}\biggr) x </math>

<math>~=</math>

<math>~\biggl(\frac{P_0}{P_c}\biggr)\frac{1}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) + \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr) \cdot \frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \, . </math>

Now, suppose that the expression on the RHS is of the form,

<math>~\mathrm{RHS}</math>

<math>~=</math>

<math>~u dv + v du \, ,</math>

where <math>~u \equiv (P_0/P_c)</math>? Then the function,

<math>~v \equiv \frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \, ,</math>

and the eigenvector, <math>~x</math>, must satisfy both of the relations:

Relation I

<math>~:</math>

<math>~- \sigma^2 \biggl(\frac{\rho_0}{\rho_c}\biggr) x = \frac{d}{d\chi_0}\biggl[ \biggl(\frac{P_0}{P_c}\biggr)\frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr)\biggr] </math>

Relation II

<math>~:</math>

<math>~\frac{d}{d\chi_0}\biggl[ \frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr)\biggr] = \frac{1}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) </math>

The validity (or not) of this conjecture can be tested against both configurations whose --- ABANDON!

Exploration

Compare LAWE to Hydrostatic Balance Condition

Returning to the generic formulation derived earlier, we have,

<math>~- \sigma^2 \biggl(\frac{\rho_0}{\rho_c}\biggr) x </math>

<math>~=</math>

<math>~\biggl(\frac{P_0}{P_c}\biggr)\frac{1}{\chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) + \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr) \cdot \frac{1}{\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \, . </math>

Dividing this entire expression through by <math>~(P_0/P_c)x</math> gives,

<math>~- \sigma^2 \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} </math>

<math>~=</math>

<math>~\frac{1}{x \chi_0^4} \frac{d}{d\chi_0}\biggl( \chi_0^4 x^' \biggr) + \frac{d}{d\chi_0}\biggl[ \ln \biggl(\frac{P_0}{P_c}\biggr)\biggr] \cdot \frac{1}{x\chi_0^\alpha}\frac{d}{d\chi_0}\biggl(\chi_0^\alpha x\biggr) \, . </math>

 

<math>~=</math>

<math>~\biggl( \frac{x^'}{x} \biggr) \cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \chi_0^4 x^' \biggr) \biggr] + \frac{d}{d\chi_0}\biggl[ \ln \biggl(\frac{P_0}{P_c}\biggr)\biggr] \cdot \frac{d}{d\chi_0}\biggl[\ln\biggl(\chi_0^\alpha x\biggr)\biggr] \, . </math>

Now, let's step aside from the LAWE and look directly at the differential relationship between the mass-density and the pressure, as dictated by combining the two principal governing relations, the

Hydrostatic Balance

<math>\frac{1}{\rho}\frac{dP}{dr} =- \frac{d\Phi}{dr} </math> ,

and,

Poisson Equation

<math>\frac{1}{r^2} \frac{d }{dr} \biggl( r^2 \frac{d \Phi}{dr} \biggr) = 4\pi G \rho </math> .

In combination, we have,

<math>~-4\pi G \rho_0 </math>

<math>~=</math>

<math>~\frac{1}{r_0^2}\frac{d}{dr_0}\biggl[ \frac{r_0^2}{\rho_0} \frac{dP_0}{dr_0}\biggr] </math>

<math>~\Rightarrow ~~~~ -4\pi G \rho_0 \biggl(\frac{R^2\rho_c}{P_c}\biggr)</math>

<math>~=</math>

<math>~\frac{1}{\chi_0^2}\frac{d}{d\chi_0}\biggl[ \frac{\chi_0^2}{(\rho_0/\rho_c)}\cdot \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr)\biggr] </math>

<math>~\Rightarrow ~~~~ -[4\pi G \rho_c \tau_\mathrm{SSC}^2] \biggl( \frac{\rho_0}{\rho_c} \biggr)</math>

<math>~=</math>

<math>~ \biggl( \frac{\rho_0}{\rho_c} \biggr)^{-1}\frac{1}{\chi_0^2}\frac{d}{d\chi_0}\biggl[ \chi_0^2 \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr)\biggr] + \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr) \cdot \frac{d}{d\chi_0}\biggl( \frac{\rho_0}{\rho_c}\biggr)^{-1} </math>

<math>~\Rightarrow ~~~~ -[4\pi G \rho_c \tau_\mathrm{SSC}^2] \biggl( \frac{\rho_0}{\rho_c} \biggr)^2\biggl(\frac{P_0}{P_c}\biggr)^{-1}</math>

<math>~=</math>

<math>~ \frac{1}{\chi_0^2 (P_0/P_c)}\frac{d}{d\chi_0}\biggl[ \chi_0^2 \frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr)\biggr] + \frac{d}{d\chi_0}\biggl[\ln\biggl(\frac{P_0}{P_c}\biggr)\biggr] \cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \frac{\rho_0}{\rho_c}\biggr)^{-1}\biggr] </math>

 

<math>~=</math>

<math>~\biggl(\frac{p^'}{p}\biggr) \frac{1}{\chi_0^2 p^'}\frac{d}{d\chi_0}\biggl[ \chi_0^2 p^'\biggr] + \frac{d\ln(p)}{d\chi_0}\cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \frac{\rho_0}{\rho_c}\biggr)^{-1}\biggr] </math>

 

<math>~=</math>

<math>~\frac{d\ln(p)}{d\chi_0}\cdot \frac{d}{d\chi_0}\biggl[ \ln(\chi_0^2 p^')\biggr] + \frac{d\ln(p)}{d\chi_0}\cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \frac{\rho_0}{\rho_c}\biggr)^{-1}\biggr] \, , </math>

where,

<math>~p^'</math>

<math>~\equiv</math>

<math>~\frac{d}{d\chi_0}\biggl(\frac{P_0}{P_c}\biggr) \, .</math>

Let's compare the form of this "equilibrium" relation with the form of the LAWE just constructed:

<math>~- [4\pi G \rho_c \tau_\mathrm{SSC}^2] \biggl( \frac{\rho_0}{\rho_c} \biggr)^2\biggl(\frac{P_0}{P_c}\biggr)^{-1}</math>

<math>~=</math>

<math>~\biggl( \frac{p^'}{p} \biggr) \cdot \frac{d}{d\chi_0}\biggl[ \ln(\chi_0^2 p^')\biggr] + \frac{d\ln(p)}{d\chi_0}\cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \frac{\rho_0}{\rho_c}\biggr)^{-1}\biggr] </math>

 

versus

 

<math>~- \sigma^2 \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} </math>

<math>~=</math>

<math>~\biggl( \frac{x^'}{x} \biggr) \cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \chi_0^4 x^' \biggr) \biggr] + \frac{d\ln(p)}{d\chi_0}\cdot \frac{d}{d\chi_0}\biggl[\ln\biggl(\chi_0^\alpha x\biggr)\biggr] </math>

I like this layout because it unveils similarities in the way the differential operators interact with the functions that describe the radial profiles of variables — specifically, the mass-density, the pressure, and the fractional radial displacement, <math>~x</math>, during pulsations. However, it is not yet obvious how best to translate between the two differential equations in order to aid in solving for the unknown variable, <math>~x(\chi_0)</math>.

Dabbling with Equilibrium Condition

In the meantime, I've found it instructive to play with the first of these two expressions to see how it might be restructured in order to most directly confirm that it is satisfied by the expressions presented in Table 1. Adopting the shorthand notation,

<math>~\Gamma \equiv 4\pi G\rho_c \tau_\mathrm{SSC}^2</math>         and         <math>~\varpi \equiv \frac{\rho_0}{\rho_c} \, ,</math>

and multiplying the "equilibrium" relation through by <math>~(-\varpi p)</math>, we have,

<math>~\Gamma \varpi^3</math>

<math>~=</math>

<math>~ - \varpi p^'\biggl\{ \frac{1}{\chi_0^2 p^'} \frac{d}{d\chi_0} (\chi_0^2 p^') - \frac{1}{\varpi}\frac{d\varpi}{d\chi_0} \biggr\} </math>

 

<math>~=</math>

<math>~ p^' \varpi^' - \varpi \frac{dp^'}{d\chi_0} -\frac{2\varpi p^'}{\chi_0} </math>

 

<math>~=</math>

<math>~ \frac{p^' }{\chi_0}\biggl[ \chi_0 \varpi^' - 2\varpi \biggr] - \varpi \frac{dp^'}{d\chi_0} \, ; </math>

or,

<math>~\Gamma \varpi^2</math>

<math>~=</math>

<math>~ \frac{p^' }{\chi_0 \varpi}\biggl[ \chi_0 \varpi^' - 2\varpi \biggr] - \frac{dp^'}{d\chi_0} \, . </math>

Specific Cases

Case 1 (Parabolic):

<math>~\varpi = 1 -\chi_0^2</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\varpi^' = -2\chi_0 </math>

<math>~p^' = -5\chi_0 + 8\chi_0^3 - 3\chi_0^5 = \chi_0(1-\chi_0^2)(-5+3\chi_0^2)</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\frac{p^'}{\chi_0 \varpi} = -5 +3\chi_0^2 \, .</math>

 

Also, note:

<math>~\frac{d(p^')}{d\chi_0} = -5 +24\chi_0^2 -15\chi_0^4 \, .</math>

For the parabolic case, therefore, the RHS of the "equilibrium" expression is,

RHS

<math>~=</math>

<math>~ (-5+3\chi_0^2)\biggl[ -2\chi_0^2 - 2(1-\chi_0^2) \biggr] - (-5 +24\chi_0^2 -15\chi_0^4) </math>

 

<math>~=</math>

<math>~ (10 - 6\chi_0^2) + (5 -24\chi_0^2 +15\chi_0^4) </math>

 

<math>~=</math>

<math>~15(1-2\chi_0^2+\chi_0^4) \, , </math>

which, indeed, matches the LHS of the "equilibrium" relation, if,

<math>~\Gamma = 15</math>        <math>~\Rightarrow</math>        <math>~\tau_\mathrm{SSC}^2 = \frac{15}{4\pi G \rho_c} \, .</math>

This has all worked satisfactorily because, as presented above, this is the correct value of <math>~\tau_\mathrm{SSC}^2</math> in the case of the parabolic density distribution.


Case 2 (Linear):

<math>~\varpi = 1 -\chi_0</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\varpi^' = -1 </math>

<math>~p^' = \tfrac{12}{5}[- 4\chi_0 + 7\chi_0^2 - 3\chi_0^3]</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\frac{p^'}{\chi_0 \varpi} = \tfrac{12}{5}(-4 +3\chi_0) \, .</math>

 

Also, note:

<math>~\frac{d(p^')}{d\chi_0} = \tfrac{12}{5}[- 4 + 14\chi_0 - 9\chi_0^2] \, .</math>

For the linear case, therefore, the RHS of the "equilibrium" expression is,

RHS

<math>~=</math>

<math>~ \tfrac{12}{5}(-4 +3\chi_0)\biggl[ -\chi_0 - 2(1-\chi_0) \biggr] - \tfrac{12}{5}(- 4 + 14\chi_0 - 9\chi_0^2) </math>

 

<math>~=</math>

<math>~ \tfrac{12}{5}\biggl[ (4 -3\chi_0)( 2-\chi_0 ) + (4 - 14\chi_0 + 9\chi_0^2) \biggr] </math>

 

<math>~=</math>

<math>~ \tfrac{12}{5}(12-24\chi_0 + 12\chi_0^2) </math>

 

<math>~=</math>

<math>~ \tfrac{2^4\cdot 3^2}{5}(1-2\chi_0 + \chi_0^2) \, , </math>

which, indeed, matches the LHS of the "equilibrium" relation, if,

<math>~\Gamma = \frac{2^4\cdot 3^2}{5}</math>        <math>~\Rightarrow</math>        <math>~\tau_\mathrm{SSC}^2 = \frac{2^2\cdot 3^2}{5\pi G \rho_c} \, .</math>

This has all worked satisfactorily because, as presented above, this is the correct value of <math>~\tau_\mathrm{SSC}^2</math> in the case of the linear density distribution.


Case 3 (n = 1 polytrope):

<math>~\varpi = \frac{\sin(\pi\chi_0)}{\pi\chi_0}</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\varpi^' = \frac{\cos(\pi\chi_0)}{\chi_0} - \frac{\sin(\pi\chi_0)}{\pi\chi_0^2}</math>

<math>~p^' = \frac{2\sin(\pi\chi_0)}{(\pi^2\chi_0^3)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr]</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\frac{p^'}{\chi_0 \varpi} = \frac{2}{(\pi\chi_0^3)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr] \, .</math>

Also, note:      <math>~\frac{d(p^')}{d\chi_0} </math>

<math>~=</math>

<math>~ \biggl[ \frac{2\pi \cdot \cos(\pi\chi_0)}{(\pi^2\chi_0^3)} - \frac{6\sin(\pi\chi_0)}{(\pi^2\chi_0^4)} \biggr] \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr] +\frac{2\sin(\pi\chi_0)}{(\pi^2\chi_0^3)} \biggl[ \pi\cos(\pi\chi_0) - \pi^2\chi_0 \sin(\pi\chi_0) - \pi \cos(\pi\chi_0) \biggr] </math>

 

<math>~=</math>

<math>~ \frac{1}{\pi^2 \chi_0^4} \biggl\{ \biggl[ 2\pi \chi_0\cdot \cos(\pi\chi_0) - 6\sin(\pi\chi_0) \biggr] \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr] +2 \chi_0 \sin(\pi\chi_0) \biggl[ - \pi^2\chi_0 \sin(\pi\chi_0) \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{1}{\pi^2 \chi_0^4} \biggl\{ 2\pi^2\chi_0^2\cos^2(\pi\chi_0) -8\pi\chi_0 \sin(\pi\chi_0) \cos(\pi\chi_0)+6\sin^2(\pi\chi_0) - 2 \pi^2 \chi_0^2 \sin^2(\pi\chi_0) \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{2}{\pi^2 \chi_0^4} \biggl\{3\sin^2(\pi\chi_0) -4\pi\chi_0 \sin(\pi\chi_0) \cos(\pi\chi_0)+\pi^2\chi_0^2\biggl[1 - 2\sin^2(\pi\chi_0) \biggr] \biggr\} \, . </math>

For the case of an n = 1 polytropic configuration, therefore, the equilibrium requirement is,

<math>~\Gamma \varpi^2</math>

<math>~=</math>

<math>~ \frac{p^' }{\chi_0 \varpi}\biggl[ \chi_0 \varpi^' - 2\varpi \biggr] - \frac{dp^'}{d\chi_0} </math>

 

<math>~=</math>

<math>~ \frac{2}{(\pi\chi_0^3)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr]\biggl[ \cos(\pi\chi_0) - \frac{\sin(\pi\chi_0)}{\pi\chi_0} - \frac{2\sin(\pi\chi_0)}{\pi\chi_0} \biggr] </math>

 

 

<math>~ - \frac{2}{\pi^2 \chi_0^4} \biggl\{3\sin^2(\pi\chi_0) -4\pi\chi_0 \sin(\pi\chi_0) \cos(\pi\chi_0)+\pi^2\chi_0^2\biggl[1 - 2\sin^2(\pi\chi_0) \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{2}{(\pi^2\chi_0^4)} \biggl\{ \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr]\biggl[\pi\chi_0 \cos(\pi\chi_0) - 3\sin(\pi\chi_0) \biggr] </math>

 

 

<math>~ -3\sin^2(\pi\chi_0) + 4\pi\chi_0 \sin(\pi\chi_0) \cos(\pi\chi_0) - \pi^2\chi_0^2\biggl[1 - 2\sin^2(\pi\chi_0) \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{2}{(\pi^2\chi_0^4)} \biggl\{3\sin^2(\pi\chi_0) - 4\pi\chi_0\sin(\pi\chi_0)\cos(\pi\chi_0) + (\pi\chi_0)^2 \biggl[1-\sin^2(\pi\chi_0) \biggr] </math>

 

 

<math>~ -3\sin^2(\pi\chi_0) + 4\pi\chi_0 \sin(\pi\chi_0) \cos(\pi\chi_0) - \pi^2\chi_0^2\biggl[1 - 2\sin^2(\pi\chi_0) \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{2}{(\pi^2\chi_0^4)} \biggl\{ (\pi\chi_0)^2 \sin^2(\pi\chi_0) \biggr\} </math>

 

<math>~=</math>

<math>~2\pi^2 \biggl[\frac{\sin(\pi\chi_0)}{\pi\chi_0} \biggr]^2 </math>

So, the equilibrium condition is satisfied if,

<math>~\Gamma = 2\pi^2</math>        <math>~\Rightarrow</math>        <math>~\tau_\mathrm{SSC}^2 = \frac{2\pi^2}{4\pi G \rho_c} = \frac{\pi}{2G \rho_c} \, .</math>

This has all worked satisfactorily because, as presented in a separate chapter discussion, this is the correct value of <math>~\tau_\mathrm{SSC}^2</math> in the case of an n = 1 polytropic configuration.

Dabbling with LAWE

Now, let's experiment with the LAWE as presented above, that is,

<math>~- \sigma^2 \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} </math>

<math>~=</math>

<math>~\biggl( \frac{x^'}{x} \biggr) \cdot \frac{d}{d\chi_0}\biggl[\ln\biggl( \chi_0^4 x^' \biggr) \biggr] + \frac{d\ln(p)}{d\chi_0}\cdot \frac{d}{d\chi_0}\biggl[\ln\biggl(\chi_0^\alpha x\biggr)\biggr] \, . </math>

After multiplying though by <math>~(-p)</math>, this expression may be written as,

<math>~-\sigma^2 \varpi</math>

<math>~=</math>

<math>~ \frac{p}{x}\biggl[ \frac{dx^'}{d\chi_0} + \frac{4x^'}{\chi_0}\biggr] + \frac{\alpha p^'}{\chi_0}\biggl[1 + \frac{\chi_0 x^'}{\alpha x}\biggr] </math>

 

<math>~=</math>

<math>~ \biggl(\frac{p}{x}\biggr) x^{' '} + p\biggl[ \frac{4}{\chi_0}\biggr]\frac{x^'}{x} + p^'\biggr[\frac{x^'}{x}\biggr] + \frac{\alpha p^'}{\chi_0} </math>

<math>~\Rightarrow ~~~~ -\biggl[\sigma^2 \varpi + \frac{\alpha p^'}{\chi_0} \biggr]x</math>

<math>~=</math>

<math>~ px^{' '} + \biggl[ \frac{4p}{\chi_0} + p^'\biggr]x^' </math>

<math>~\Rightarrow ~~~~ 0</math>

<math>~=</math>

<math>~ px^{' '} + \biggl[ 4p + \chi_0 p^'\biggr]\frac{x^'}{\chi_0} + \biggl[\sigma^2 \varpi + \frac{\alpha p^'}{\chi_0} \biggr]x \, . </math>

(We could have, perhaps, obtained this expression in a more direct fashion had we started directly from the form of the LAWE derived earlier.)

Specific Case Attempts
Uniform Density

Case 0 (Uniform density):

<math>~\varpi</math>

<math>~=</math>

<math>~1 \, ;</math>

<math>~p</math>

<math>~=</math>

<math>~1-\chi_0^2 \, ;</math>

<math>~\frac{\alpha p^'}{\chi_0}</math>

<math>~=</math>

<math>~- 2\alpha \, .</math>

For the uniform-density case, therefore, the the LAWE becomes,

<math>~-\sigma^2 </math>

<math>~=</math>

<math>~ \frac{(1-\chi_0^2)}{x}\biggl[ \frac{dx^'}{d\chi_0} + \frac{4x^'}{\chi_0}\biggr] -2\alpha \biggl[1 + \frac{\chi_0 x^'}{\alpha x}\biggr] </math>

<math>~\Rightarrow~~~~ (2\alpha -\sigma^2)x </math>

<math>~=</math>

<math>~ (1-\chi_0^2)\biggl[ \frac{dx^'}{d\chi_0} + \frac{4x^'}{\chi_0}\biggr] - 2\chi_0 x^' </math>

 

<math>~=</math>

<math>~ (1-\chi_0^2)\frac{dx^'}{d\chi_0} + (1-\chi_0^2)\biggl[ \frac{4x^'}{\chi_0}\biggr] - 2\chi_0 x^' </math>

 

<math>~=</math>

<math>~ (1-\chi_0^2)\frac{dx^'}{d\chi_0} + \frac{2x^'}{\chi_0}\biggl( 2 - 3\chi_0^2 \biggr) \, , </math>

where, as defined above,

<math>~\alpha \equiv 3 - \frac{4}{\gamma_g} \, .</math>

Mode 3

Try an eigenfunction of the form,

<math>x = a + b\chi_0^2 + c\chi_0^4 \, ,</math>

in which case,

<math>~\frac{2x^'}{\chi_0} = \frac{2}{\chi_0}(2 b\chi_0 +4c\chi_0^3) = 4b+8c\chi_0^2</math>         and         <math>~x^{' '} = 2 b + 12c\chi_0^2 \, . </math>

In order for this to be a solution, we must have,

<math>~(2\alpha -\sigma^2)( a + b\chi_0^2 + c\chi_0^4) </math>

<math>~=</math>

<math>~ (1-\chi_0^2)(2 b + 12c\chi_0^2 ) + ( 2 - 3\chi_0^2 )(4b+8c\chi_0^2 ) </math>

 

<math>~=</math>

<math>~ (2 b + 12c\chi_0^2) - \chi_0^2(2 b + 12c\chi_0^2 ) + 2(4b+8c\chi_0^2 ) - 3\chi_0^2 (4b+8c\chi_0^2 ) </math>

 

<math>~=</math>

<math>~10b + \chi_0^2(12c-2b+16c-12b) - \chi_0^4(12c + 24c) </math>

 

<math>~=</math>

<math>~10b + \chi_0^2(28c-14b) - \chi_0^4(36c) \, . </math>


So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~a\mathfrak{F} +10b</math>

<math>~\chi_0^2</math>   :  

<math>~b\mathfrak{F} -14b + 28c</math>

<math>~\chi_0^4</math>   :  

<math>~c[\mathfrak{F} -36]</math>

where, following Sterne's (1937) presentation,

<math>~\mathfrak{F} \equiv \sigma^2 - 2 \alpha \, .</math>

In order for all three of the coefficients to be zero, we must have:

First:     <math>~\mathfrak{F} = 36 \, ;</math>

Second:     <math>~22b = -28c ~~~~~\Rightarrow ~~~~~ c = - (11/14)b \, ;</math>

Third:     <math>~36a = -10b ~~~~~\Rightarrow ~~~~~ b = -(18/5)a \, .</math>

Hence, choosing <math>~a = 1</math> implies: <math>~b = -18/5</math>        and         <math>~ c = (11/7)(9/5) = +99/35 \, .</math> This precisely matches the "j = 2" mode identified by Sterne.


Parabolic

Case 1 (Parabolic):

<math>~\varpi = 1 -\chi_0^2</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\varpi^' = -2\chi_0 </math>

<math>~p = \tfrac{1}{2}(1 -\chi_0^2)^2 (2-\chi_0^2)</math>

 

 

<math>~p^' = -5\chi_0 + 8\chi_0^3 - 3\chi_0^5 = \chi_0(1-\chi_0^2)(-5+3\chi_0^2)</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\frac{\alpha p^'}{\chi_0} = \alpha (1-\chi_0^2)(-5+3\chi_0^2) \, .</math>


For the parabolic case, therefore, the the LAWE becomes,

<math>~0</math>

<math>~=</math>

<math>~ px^{' '} + \biggl[ 4p + \chi_0 p^'\biggr]\frac{x^'}{\chi_0} + \biggl[\sigma^2 \varpi + \frac{\alpha p^'}{\chi_0} \biggr]x </math>

 

<math>~=</math>

<math>~ \tfrac{1}{2}\varpi^2 (2-\chi_0^2)x^{' '} + \biggl[ 2\varpi^2 (2-\chi_0^2) + \chi_0^2\varpi(-5+3\chi_0^2)\biggr]\frac{x^'}{\chi_0} + \biggl[\sigma^2 \varpi + \alpha \varpi(-5+3\chi_0^2) \biggr]x </math>

 

<math>~=</math>

<math>~\frac{\varpi}{2}\biggl\{ \varpi (2-\chi_0^2)x^{' '} + \biggl[ 4\varpi (2-\chi_0^2) + \chi_0^2(-10+6\chi_0^2)\biggr]\frac{x^'}{\chi_0} + \biggl[\mathfrak{K}+6\alpha\chi_0^2 \biggr]x \biggr\} </math>

 

<math>~=</math>

<math>~\frac{\varpi}{2}\biggl\{ (2-3\chi_0^2 + \chi_0^4)x^{' '} + \biggl[ 8-22\chi_0^2 + 10\chi_0^4\biggr]\frac{x^'}{\chi_0} + \biggl[\mathfrak{K}+6\alpha\chi_0^2 \biggr]x \biggr\} </math>

where,

<math>~\mathfrak{K} \equiv 2(\sigma^2 - 5\alpha) \, .</math>

Mode Inverse

Try an eigenfunction of the form,

<math>x = (1 + a\chi_0^2)^{-\beta} = (1 + a\chi_0^2)^{-(\beta+2)} (1 + 2a\chi_0^2 + a^2\chi_0^4)\, ,</math>

in which case,

<math>~x^'</math>

<math>~=</math>

<math>~ - 2a\beta \chi_0 (1 + a\chi_0^2)^{-\beta-1} \, ; </math>

 

<math>~=</math>

<math>~ - 2a\beta (\chi_0 + a\chi_0^3)(1 + a\chi_0^2)^{-(\beta+2)} \, ; </math>

<math>~x^{' '}</math>

<math>~=</math>

<math>~ - 2a\beta (1 + a\chi_0^2)^{-\beta-1} - 2a\beta \chi_0\biggl[ -2a(\beta+1) \chi_0(1 + a\chi_0^2)^{-\beta-2}\biggr] </math>

 

<math>~=</math>

<math>~ - 2a\beta (1 + a\chi_0^2)^{-(\beta+2)}\biggl[(1 + a\chi_0^2) -2a(\beta+1) \chi_0^2\biggr] </math>

 

<math>~=</math>

<math>~ - 2a\beta (1 + a\chi_0^2)^{-(\beta+2)}\biggl[1 - a(2\beta+1)\chi_0^2 \biggr] </math>

In order for this to be a solution, we must have,

<math>~0</math>

<math>~=</math>

<math>~ (2-3\chi_0^2 + \chi_0^4)x^{' '} + \biggl[ 8-22\chi_0^2 + 10\chi_0^4\biggr]\frac{x^'}{\chi_0} + \biggl[\mathfrak{K}+6\alpha\chi_0^2 \biggr]x </math>

 

<math>~=</math>

<math>~ - 2a\beta \biggl[1 - a(2\beta+1)\chi_0^2 \biggr](2-3\chi_0^2 + \chi_0^4) - 2a\beta (1 + a\chi_0^2) \biggl[ 8-22\chi_0^2 + 10\chi_0^4\biggr] + \biggl[\mathfrak{K}+6\alpha\chi_0^2 \biggr](1 + 2a\chi_0^2 + a^2\chi_0^4) </math>

 

<math>~=</math>

<math>~ 2a\beta (-2+3\chi_0^2 - \chi_0^4) + 2a^2\beta (2\beta+1) (2\chi_0^2-3\chi_0^4 + \chi_0^6) </math>

 

 

<math>~ + 2a\beta \biggl[ -8+22\chi_0^2 - 10\chi_0^4\biggr] + 2a^2\beta \biggl[ -8\chi_0^2+22\chi_0^4 - 10\chi_0^6\biggr] </math>

 

 

<math>~ + \mathfrak{K} (1 + 2a\chi_0^2 + a^2\chi_0^4) + 6\alpha (\chi_0^2 + 2a\chi_0^4 + a^2\chi_0^6) </math>


So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~\mathfrak{K} - 20a\beta </math>

<math>~\chi_0^2</math>   :  

<math>~50a\beta + 4a^2\beta (2\beta-3) + 6\alpha + 2a\mathfrak{K}</math>

<math>~\chi_0^4</math>   :  

<math>~- 22a\beta + 2a^2\beta (19-6\beta) + 12a\alpha + a^2\mathfrak{K}</math>

<math>~\chi_0^6</math>   :  

<math>~2a^2\beta (2\beta-9) + 6a^2\alpha</math>


Mode 3P

Try an eigenfunction of the form,

<math>x = a + b\chi_0^2 + c\chi_0^4 \, ,</math>

in which case,

<math>~\frac{x^'}{\chi_0} = \frac{1}{\chi_0}(2 b\chi_0 +4c\chi_0^3) = 2b+4c\chi_0^2</math>         and         <math>~x^{' '} = 2 b + 12c\chi_0^2 \, . </math>

In order for this to be a solution, we must have,

<math>~0</math>

<math>~=</math>

<math>~\frac{\varpi}{2}\biggl\{ (2-3\chi_0^2 + \chi_0^4)(2 b + 12c\chi_0^2 ) + (8-22\chi_0^2 + 10\chi_0^4 )(2b+4c\chi_0^2 ) + (\mathfrak{K}+6\alpha\chi_0^2 )( a + b\chi_0^2 + c\chi_0^4 ) \biggr\} </math>

 

<math>~=</math>

<math>~\frac{\varpi}{2}\biggl\{ (4b+24c\chi_0^2 - 6b\chi_0^2 -36c\chi_0^4 + 2b\chi_0^4 + 12c\chi_0^6) + (16b + 32c\chi_0^2 - 44b\chi_0^2 - 88c\chi_0^4 + 20b\chi_0^4 + 40c\chi_0^6 ) </math>

 

 

<math>~+ (a\mathfrak{K} + b\mathfrak{K}\chi_0^2 + c\mathfrak{K}\chi_0^4 + 6a\alpha\chi_0^2 + 6b\alpha\chi_0^4 + 6c\alpha\chi_0^6 ) \biggr\} </math>

 

<math>~=</math>

<math>~\frac{\varpi}{2}\biggl[ (20b + a\mathfrak{K}) \chi_0^0 +(24c - 6b + 32c - 44b + b\mathfrak{K} + 6a\alpha)\chi_0^2 + (-36c+2b -88c+20b +c\mathfrak{K} + 6b\alpha) \chi_0^4 + (12c + 40c + 6c\alpha )\chi_0^6 \biggr] </math>


So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~20b + a\mathfrak{K}</math>

<math>~\chi_0^2</math>   :  

<math>~56c + (\mathfrak{K}- 50)b + 6a\alpha</math>

<math>~\chi_0^4</math>   :  

<math>~b(22+6\alpha) + c(\mathfrak{K}-124)</math>

<math>~\chi_0^6</math>   :  

<math>~c(52 + 6\alpha)</math>


This is disappointing, as it does not result in nonzero coefficient values.

Polytrope

Case 3 (n = 1 polytrope):

<math>~\varpi = \frac{\sin(\pi\chi_0)}{\pi\chi_0}</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\varpi^' = \frac{\cos(\pi\chi_0)}{\chi_0} - \frac{\sin(\pi\chi_0)}{\pi\chi_0^2}</math>

<math>~p = \biggl[\frac{\sin(\pi\chi_0)}{\pi\chi_0}\biggr]^2 = \varpi^2</math>

 

 

<math>~p^' = \frac{2\varpi}{(\pi \chi_0^2)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr]</math>

      <math>~~~~\Rightarrow~~~~</math>     

<math>~\frac{\alpha p^'}{\chi_0 } = \frac{2\alpha \varpi}{(\pi \chi_0^3)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr] \, .</math>


For the n = 1 polytropic case, therefore, the the LAWE becomes,

<math>~0</math>

<math>~=</math>

<math>~ px^{' '} + \biggl[ 4p + \chi_0 p^'\biggr]\frac{x^'}{\chi_0} + \biggl[\sigma^2 \varpi + \frac{\alpha p^'}{\chi_0} \biggr]x </math>

 

<math>~=</math>

<math>~ \biggl[\frac{\sin(\pi\chi_0)}{\pi\chi_0}\biggr] \varpi x^{' '} + \biggl\{ 4\biggl[\frac{\sin(\pi\chi_0)}{\pi\chi_0}\biggr] + \frac{2}{(\pi \chi_0)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr]\biggr\}\frac{\varpi x^'}{\chi_0} </math>

 

 

<math>~ + \biggl\{\sigma^2 + \frac{2\alpha }{(\pi \chi_0^3)} \biggl[ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) \biggr] \biggr\}\varpi x </math>

 

<math>~=</math>

<math>~ \frac{\varpi}{\pi \chi_0^2} \biggl\{ [\sin(\pi\chi_0)] \chi_0 x^{' '} + [ 2\sin(\pi\chi_0) + 2\pi\chi_0 \cos(\pi\chi_0) ] x^' + \{(\pi \chi_0^2)\sigma^2 + 2\alpha \chi_0^{-1} [ \pi\chi_0 \cos(\pi\chi_0) - \sin(\pi\chi_0) ] \} x \biggr\} </math>


Exploration2

Let's begin with the LAWE written in the following form (see, for example, the related context discussion):

<math>~- \sigma^2 \mathcal{G}_\sigma</math>

<math>~=</math>

<math>~ \biggl(\frac{P}{\rho }\biggr)\frac{1}{x^4} \frac{d}{dx}\biggl( x^4 \mathcal{G}_\sigma^' \biggr) + \biggl(\frac{P^'}{\rho }\biggr) \frac{1}{x^\alpha}\frac{d}{dx}\biggl(x^\alpha \mathcal{G}_\sigma\biggr) \, . </math>

One advantage of beginning with this construction is that — as the following table shows — it might be reasonable to expect in general that both pre-factors — <math>~(P/\rho )</math> and <math>~(P^'/\rho )</math> — will have a relatively simple mathematical form. In addition, however, it appears as though both terms on the RHS want to be logarithmic derivatives.

Properties of Analytically Defined Astrophysical Structures
Model <math>~\rho(x)</math> <math>~\biggl[\frac{P(x)}{\rho(x)}\biggr]</math> <math>~\biggl[ \frac{P^'(x)}{\rho(x)} \biggr]</math>
Uniform-density <math>~1</math> <math>~1 - x^2</math> <math>~-2x</math>
Linear <math>~1-x</math> <math>~(1-x)(1 + 2x - \tfrac{9}{5}x^2)</math> <math>~-\tfrac{12}{5}x (4-3x)</math>
Parabolic <math>~1-x^2</math> <math>~(1-x^2)(1 - \tfrac{1}{2} x^2)</math> <math>~-x (5-3x^2)</math>
<math>~n=1</math> Polytrope <math>~\frac{\sin }{ x}</math> <math>~\frac{\sin x}{x}</math> <math>~\frac{2}{x} \biggl[ \cos x - \frac{\sin x}{x} \biggr] </math>

New Form of LAWE

Multiplying the LAWE through by <math>~[\rho/(P\mathcal{G}_\sigma)]</math> gives,

<math>~- \sigma^2 \biggl(\frac{\rho}{P}\biggr)</math>

<math>~=</math>

<math>~ \biggl(\frac{\mathcal{G}_\sigma^'}{\mathcal{G}_\sigma }\biggr)\frac{1}{(x^4 \mathcal{G}_\sigma^')} \frac{d}{dx}\biggl( x^4 \mathcal{G}_\sigma^' \biggr) + \biggl(\frac{P^'}{P}\biggr) \frac{1}{(x^\alpha \mathcal{G}_\sigma)}\frac{d}{dx}\biggl(x^\alpha \mathcal{G}_\sigma\biggr) </math>

 

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{dx } \cdot \frac{d\ln( x^4 \mathcal{G}_\sigma^' )}{dx} + \frac{d\ln P}{dx } \cdot \frac{d\ln( x^\alpha \mathcal{G}_\sigma )}{dx} </math>

 

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{dx } \biggl[ \frac{d\ln( x^4 )}{dx} + \frac{d\ln( \mathcal{G}_\sigma^' )}{dx} \biggr] + \frac{d\ln P}{dx } \biggl[ \frac{d\ln( x^\alpha )}{dx} + \frac{d\ln( \mathcal{G}_\sigma )}{dx} \biggr] </math>

 

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{dx } \biggl[ \frac{d\ln( x^4 )}{dx} + \frac{d\ln( \mathcal{G}_\sigma^' )}{dx} + \frac{d\ln P}{dx } \biggr] + \frac{d\ln P}{dx } \biggl[ \frac{d\ln( x^\alpha )}{dx} \biggr] </math>

<math>~\Rightarrow ~~~~~ - \biggl[ \sigma^2 \biggl(\frac{\rho}{P}\biggr) + \frac{d\ln P}{dx } \cdot \frac{d\ln( x^\alpha )}{dx} \biggr]</math>

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{dx } \biggl[ \frac{d\ln( x^4 P \mathcal{G}_\sigma^' )}{dx} \biggr] </math>

<math>~\Rightarrow ~~~~~ - \frac{\mathcal{G}_\sigma}{P}\biggl[ \sigma^2 \rho + \frac{\alpha P^'}{x} \biggr]</math>

<math>~=</math>

<math>~ \frac{1}{x^4 P } \biggl[ \frac{d( x^4 P \mathcal{G}_\sigma^' )}{dx} \biggr] </math>

<math>~\Rightarrow ~~~~~ \frac{d( x^4 P \mathcal{G}_\sigma^' )}{dx} </math>

<math>~=</math>

<math>~- \mathcal{G}_\sigma \biggl[ \sigma^2 x^4 \rho + \alpha x^3 P^' \biggr] </math>

 

<math>~=</math>

<math>~-x^4 \rho \mathcal{G}_\sigma \biggl[ \sigma^2 + \frac{\alpha P^'}{x\rho} \biggr] \, .</math>

Trial Logarithmic Eigenfunction

Defining,

<math>\mathcal{F}(x) \equiv \biggl[ \sigma^2 + \frac{\alpha P^'}{x\rho} \biggr] \, ,</math>

the LAWE becomes,

<math>~ \frac{d( x^4 P \mathcal{G}_\sigma^' )}{dx} </math>

<math>~=</math>

<math>~-x^4 \rho \mathcal{F} \mathcal{G}_\sigma \, .</math>

First Try

Let's try an eigenvector of the form,

<math>~\mathcal{G}_\sigma </math>

<math>~=</math>

<math>~A(x) \ln P + B(x) \, ,</math>

in which case,

<math>~\mathcal{G}_\sigma^' </math>

<math>~=</math>

<math>~A^' \ln P + \frac{A\cdot P^'}{P}+ B^' </math>

<math>~\Rightarrow ~~~~ x^4 P \mathcal{G}_\sigma^' </math>

<math>~=</math>

<math>~x^4\biggl[ PA^' \ln P + A\cdot P^'+ PB^' \biggr]</math>

<math>~\Rightarrow ~~~~\mathrm{LHS} ~~\equiv \frac{d}{dx}\biggl[x^4 P \mathcal{G}_\sigma^' \biggr]</math>

<math>~=</math>

<math>~ x^4\biggl[ P'A^' \ln P + PA^{' '} \ln P + PA^' \biggl(\frac{P^'}{P}\biggr) + A'\cdot P^'+ PB^{' '} + A\cdot P^{' '}+ P'B^' \biggr] </math>

 

 

<math>~ + 4x^3\biggl[ PA^' \ln P + A\cdot P^'+ PB^' \biggr] </math>

 

<math>~=</math>

<math>~x^4 \biggl\{\ln P \biggl[P'A^' + PA^{' '} + \frac{4}{x} \cdot PA^' \biggr] + \biggl[ A^' P^' + A'\cdot P^'+ PB^{' '} + A\cdot P^{' '}+ P'B^' + \frac{4}{x} \biggl( A\cdot P^'+ PB^' \biggr) \biggr]\biggr\} \, . </math>

Now, in order for this expression to match the RHS of the LAWE, we must have, first of all,

<math>~\mathcal{F} A</math>

<math>~=</math>

<math>~-~\frac{1}{\rho}\biggl[ P'A^' + PA^{' '} + \frac{4}{x} \cdot PA^'\biggr] \, ;</math>

and, second,

<math>~\mathcal{F} B</math>

<math>~=</math>

<math>~-~\frac{1}{\rho}\biggl[A^' P^' + A'\cdot P^'+ PB^{' '} + A\cdot P^{' '}+ P'B^' + \frac{4}{x} \biggl( A\cdot P^'+ PB^' \biggr) \biggr] \, .</math>

Case 1 (Parabolic):

<math>\mathcal{F}(x) \equiv \biggl[ \sigma^2 - \alpha(5 - 3x^2) \biggr] = f_0 + f_2 x^2 \, ,</math>

where,

<math>f_0 \equiv \sigma^2 - 5\alpha</math>         and         <math>f_2 \equiv 3\alpha \, .</math>

Also, the first condition is,

<math>~- 2( f_0 + f_2 x^2 )A</math>

<math>~=</math>

<math>~2( -5 + 3x^2 ) xA^' + (2-3x^2 + x^4) \biggl[A^{' '} + \frac{4A^'}{x} \biggr] \, .</math>

So, if we adopt a polynomial expression for the function, <math>~A(x)</math>, of the form,

<math>~A(x) = a_0 + a_2 x^2 + a_4 x^4 \, ,</math>


<math>~\Rightarrow ~~~~ A^' = 2a_2 x + 4a_4 x^3</math>         and         <math>A^{' '} = 2a_2 + 12a_4 x^2 \, ,</math>

the condition becomes,

<math>~- 2( f_0 + f_2 x^2 )(a_0 + a_2 x^2 + a_4 x^4)</math>

<math>~=</math>

<math>~( -10 + 6x^2 ) (2a_2 x^2 + 4a_4 x^4) + (2-3x^2 + x^4) \biggl[(2a_2 + 12a_4 x^2) + 4 (2a_2 + 4a_4 x^2) \biggr] </math>

<math>~\Rightarrow ~~~~ ( f_0 + f_2 x^2 )(a_0 + a_2 x^2 + a_4 x^4)</math>

<math>~=</math>

<math>~( 10 - 6x^2 ) (a_2 x^2 + 2 a_4 x^4) - (2-3x^2 + x^4) (5a_2 + 14a_4 x^2) </math>

So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~f_0 a_0 + 10 a_2</math>

<math>~\chi_0^2</math>   :  

<math>~f_0 a_2 + f_2 a_0 -10a_2 + 28a^4 - 15a_2</math>

<math>~\chi_0^4</math>   :  

<math>~f_0 a_4 + f_2 a_2 -20a_4 +6a_2 -42a_4 + 5a_2</math>

<math>~\chi_0^6</math>   :  

<math>~f_2 a_4 +12a_4 +14a_4</math>

This does not seem to work.

Another Trial

Start with a form of the LAWE found midway through the above derivation, namely,

<math>~- \biggl[ \sigma^2 \biggl(\frac{\rho}{P}\biggr) + \frac{d\ln P}{dx } \cdot \frac{d\ln( x^\alpha )}{dx} \biggr]</math>

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{dx } \biggl[ \frac{d\ln( x^4 P \mathcal{G}_\sigma^' )}{dx} \biggr] \, ; </math>

and multiplying through by <math>~x^2</math> gives,

<math>~- \biggl[ \sigma^2 \biggl(\frac{x^2 \rho}{P}\biggr) + \alpha \cdot \frac{d\ln P}{d\ln x } \biggr]</math>

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{d\ln x } \biggl[ \frac{d\ln( x^4 P \mathcal{G}_\sigma^' )}{d\ln x} \biggr] </math>

<math>~\Rightarrow ~~~~ -\frac{d\ln P}{d\ln x } \biggl[ \sigma^2 \biggl(\frac{P^'}{x \rho}\biggr)^{-1} + \alpha \biggr]</math>

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{d\ln x } \biggl[ \frac{d\ln( \mathcal{G}_\sigma^' )}{d\ln x} + \frac{d\ln( P )}{d\ln x} + \frac{d\ln( x^4 )}{d\ln x} \biggr] </math>

<math>~\Rightarrow ~~~~ -\frac{d\ln P}{d\ln x } \biggl[\alpha + \sigma^2 \biggl(\frac{P^'}{x \rho}\biggr)^{-1} + \frac{d\ln \mathcal{G}_\sigma}{d\ln x } \biggr]</math>

<math>~=</math>

<math>~ \frac{d\ln \mathcal{G}_\sigma}{d\ln x } \biggl[ \frac{d\ln( \mathcal{G}_\sigma^' )}{d\ln x} +4 \biggr] \, . </math>

Note that,

<math>~\frac{d\ln \mathcal{G}_\sigma}{d\ln x } </math>

<math>~=</math>

<math>~\frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} \, ;</math>

<math>~\frac{d\ln \mathcal{G}_\sigma}{d\ln x } \cdot \frac{d\ln( \mathcal{G}_\sigma^' )}{d\ln x}</math>

<math>~=</math>

<math>~\frac{x^2 \mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma} \, .</math>

Consider Parabolic Case

In the case of a parabolic density distribution, the LAWE becomes,

<math>~\frac{2x^2(5-3x^2)}{(1-x^2)(2-x^2)} \biggl[\alpha - \sigma^2 \biggl(5-3x^2\biggr)^{-1} + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} \biggr]</math>

<math>~=</math>

<math>~ \biggl(\frac{x^2 \mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +4 \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} </math>

<math>~\Rightarrow ~~~~ \frac{2}{(1-x^2)(2-x^2)} \biggl[ \biggl( \alpha + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}\biggr)(5-3x^2) -\sigma^2 \biggr]</math>

<math>~=</math>

<math>~ \biggl(\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +\frac{4}{x^2} \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} </math>


Let's try,

<math>~\mathcal{G}_\sigma</math>

<math>~=</math>

<math>~(a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^m \, ,</math>

which implies,

<math>~\mathcal{G}_\sigma^'</math>

<math>~=</math>

<math>~n(a_0 + a_2x^2)^{n-1}(2a_2x) \cdot (b_0 + b_2x^2)^m +m (a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^{m-1}(2b_2x)</math>

<math>~\Rightarrow ~~~~ \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}</math>

<math>~=</math>

<math>~n(a_0 + a_2x^2)^{-1}(2a_2x^2) +m (b_0 + b_2x^2)^{-1}(2b_2x^2) </math>

 

<math>~=</math>

<math>~\frac{2x^2}{(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl[ n a_2 (b_0 + b_2x^2) +mb_2 (a_0 + a_2x^2) \biggr] </math>

 

<math>~=</math>

<math>~\frac{2x^2}{(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl[ (n a_2 b_0 + mb_2 a_0) +(na_2 b_2+ mb_2 a_2)x^2\biggr] \, ,</math>

and,

<math>~\mathcal{G}_\sigma^{' '}</math>

<math>~=</math>

<math>~n m (a_0 + a_2x^2)^{n-1}(2a_2x) \cdot (b_0 + b_2x^2)^{m-1}(2b_2x) + n(a_0 + a_2x^2)^{n-1}(2a_2) \cdot (b_0 + b_2x^2)^m + n(n-1)(a_0 + a_2x^2)^{n-2}(2a_2x)^2 \cdot (b_0 + b_2x^2)^m </math>

 

 

<math>~+m n(a_0 + a_2x^2)^{n-1}(2a_2x) \cdot (b_0 + b_2x^2)^{m-1}(2b_2x) +m (a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^{m-1}(2b_2) +m(m-1) (a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^{m-2}(2b_2x)^2</math>

<math>~\Rightarrow ~~~~ \frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}</math>

<math>~=</math>

<math>~8n m a_2b_2 x^2 (a_0 + a_2x^2)^{-1}\cdot (b_0 + b_2x^2)^{-1} + n2a_2 (a_0 + a_2x^2)^{-1} + n(n-1)4a_2^2 x^2 (a_0 + a_2x^2)^{-2} +m2b_2 (b_0 + b_2x^2)^{-1} +m(m-1)4 b_2^2 x^2 (b_0 + b_2x^2)^{-2}</math>

 

<math>~=</math>

<math>~\frac{2n a_2(b_0 + b_2x^2) + 2m b_2 (a_0 + a_2x^2)}{ (a_0 + a_2x^2)(b_0 + b_2x^2)} + \biggl[ \frac{4n(n-1) a_2^2 }{ (a_0 + a_2x^2)^{2}} + \frac{8n m a_2b_2}{ (a_0 + a_2x^2)(b_0 + b_2x^2)}+ \frac{4m(m-1) b_2^2 }{(b_0 + b_2x^2)^{2}} \biggr]x^2 </math>

So, we have for the LAWE:

LHS

<math>~=</math>

<math>~ \frac{2}{(1-x^2)(2-x^2)} \biggl[ \biggl( \alpha + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}\biggr)(5-3x^2) -\sigma^2 \biggr]</math>

 

<math>~=</math>

<math>~ \frac{2}{(1-x^2)(2-x^2)(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl\{ \biggl[ \alpha(a_0 + a_2x^2) (b_0 + b_2x^2) + 2x^2(n a_2 b_0 + mb_2 a_0) + 2x^4 (na_2 b_2+ mb_2 a_2) \biggr](5-3x^2) </math>

 

 

<math>~ -\sigma^2 (a_0 + a_2x^2) (b_0 + b_2x^2) \biggr\} \, ;</math>

RHS

<math>~=</math>

<math>~ \biggl(\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +\frac{4}{x^2} \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} </math>

 

<math>~=</math>

<math>~\frac{2n a_2(b_0 + b_2x^2) + 2m b_2 (a_0 + a_2x^2)}{ (a_0 + a_2x^2)(b_0 + b_2x^2)} + \biggl[ \frac{4n(n-1) a_2^2 }{ (a_0 + a_2x^2)^{2}} + \frac{8n m a_2b_2}{ (a_0 + a_2x^2)(b_0 + b_2x^2)}+ \frac{4m(m-1) b_2^2 }{(b_0 + b_2x^2)^{2}} \biggr]x^2 </math>

 

 

<math>~ + \frac{8}{(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl[ (n a_2 b_0 + mb_2 a_0) +(na_2 b_2+ mb_2 a_2)x^2\biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{(a_0 + a_2x^2)(b_0 + b_2x^2)} \biggl\{ 2n a_2(b_0 + b_2x^2) + 2m b_2 (a_0 + a_2x^2) + 8(n a_2 b_0 + mb_2 a_0) + 8(na_2 b_2+ mb_2 a_2)x^2 </math>

 

 

<math>~ + \biggl[8n m a_2b_2+ \frac{4n(n-1) a_2^2(b_0 + b_2x^2) }{ (a_0 + a_2x^2)} + \frac{4m(m-1) b_2^2(a_0 + a_2x^2) }{(b_0 + b_2x^2)} \biggr]x^2 \biggr\} \, . </math>

Putting these together gives,

<math>~ 0 </math>

<math>~=</math>

<math>~ \biggl[ \alpha(a_0 + a_2x^2) (b_0 + b_2x^2) + 2x^2(n a_2 b_0 + mb_2 a_0) + 2x^4 (na_2 b_2+ mb_2 a_2) \biggr](5-3x^2) -\sigma^2 (a_0 + a_2x^2) (b_0 + b_2x^2) </math>

 

 

<math>~ - \biggl[ n a_2(b_0 + b_2x^2) + m b_2 (a_0 + a_2x^2) + 4(n a_2 b_0 + mb_2 a_0) + 4(na_2 b_2+ mb_2 a_2)x^2+ 4n m a_2b_2x^2 \biggr](1-x^2)(2-x^2) </math>

 

 

<math>~ - \frac{(1-x^2)(2-x^2)}{ (a_0 + a_2x^2)(b_0 + b_2x^2)}\biggl[2n(n-1) a_2^2(b_0 + b_2x^2)^2 + 2m(m-1) b_2^2(a_0 + a_2x^2)^2 \biggr]x^2 \, . </math>


First Guess

Now, if we are very lucky, we will find that,

<math>~(a_0 + a_2x^2) = (1-x^2)</math>      <math>~\Rightarrow</math>      <math>~a_0 = 1</math>      and      <math>~a_2 = -1</math>;

and, simultaneously,

<math>~(b_0 + b_2x^2) = (2-x^2)</math>      <math>~\Rightarrow</math>      <math>~b_0 = 2</math>      and      <math>~b_2 = -1</math>.

In this case, the fractional coefficient in the last term of the LAWE will become unity and the LAWE becomes,

<math>~ 0 </math>

<math>~=</math>

<math>~ \biggl[ \alpha(1-x^2) (2-x^2) - 2x^2(2n + m) + 2x^4 (n+ m) \biggr](5-3x^2) </math>

 

 

<math>~ - \biggl[\sigma^2 - n (2-x^2) - m (1-x^2) - 4(2n + m) + 4(n + m )x^2+ 4n m x^2 \biggr](1-x^2)(2-x^2) </math>

 

 

<math>~ - \biggl[2n(n-1) (2-x^2)^2 + 2m(m-1) (1-x^2)^2 \biggr]x^2 </math>

 

<math>~=</math>

<math>~ \biggl[ 2\alpha - x^2(3\alpha +4n + 2m) + x^4 (\alpha + 2n+ 2m) \biggr](5-3x^2) </math>

 

 

<math>~ + \biggl[(-\sigma^2 + 10n + 5m ) - x^2( 5n+5m+4nm )\biggr] (2-3x^2+x^4) </math>

 

 

<math>~ - \biggl\{ [8n(n-1) +2m(m-1) ] - [8n(n-1) + 4m(m-1) ] x^2 + [ 2n(n-1)+ 2m(m-1) ] x^4 \biggr\} x^2 \, . </math>

So, the coefficients of each even power of <math>~\chi_0^n</math> are:

<math>~\chi_0^0</math>   :  

<math>~10\alpha + 2(-\sigma^2 + 10n + 5m )</math>

<math>~\chi_0^2</math>   :  

<math>~-6\alpha - 5(3\alpha +4n + 2m) -2( 5n+5m+4nm ) - 3(-\sigma^2 + 10n + 5m ) - [8n(n-1) +2m(m-1) ]</math>

<math>~\chi_0^4</math>   :  

<math>~5(\alpha + 2n+ 2m)+3(3\alpha +4n + 2m) + (-\sigma^2 + 10n + 5m ) +3( 5n+5m+4nm ) + [8n(n-1) + 4m(m-1) ] </math>

<math>~\chi_0^6</math>   :  

<math>~-3(\alpha + 2n+ 2m) - ( 5n+5m+4nm ) - [ 2n(n-1)+ 2m(m-1) ] </math>


Now let's begin simplification. The <math>~x^0</math> coefficient implies,

<math>~(-\sigma^2 + 10n + 5m )=-5\alpha \, .</math>

Hence,

<math>~\chi_0^2</math>   :  

<math>~-6\alpha - 5( 4n + 2m) -2( 5n+5m+4nm ) - [8n(n-1) +2m(m-1) ]</math>

<math>~\chi_0^4</math>   :  

<math>~9\alpha + 5(2n+ 2m)+3(4n + 2m) +3( 5n+5m+4nm ) + [8n(n-1) + 4m(m-1) ] </math>

<math>~\chi_0^6</math>   :  

<math>~-3\alpha -3(2n+ 2m) - ( 5n+5m+4nm ) - [ 2n(n-1)+ 2m(m-1) ] </math>

Using the <math>~x^6</math> coefficient to define <math>~\alpha</math>, that is, setting,

<math>~3\alpha = \{-3(2n+ 2m) - ( 5n+5m+4nm ) - [ 2n(n-1)+ 2m(m-1) ]\} \, ,</math>

means that the other two coefficient expressions are,

<math>~\chi_0^2</math>   :  

<math>~6(2n+ 2m) + 2 ( 5n+5m+4nm ) + 4 [ n(n-1)+ m(m-1) ] - 5( 4n + 2m) -2( 5n+5m+4nm ) - [8n(n-1) +2m(m-1) ]</math>

<math>~\chi_0^4</math>   :  

<math>~-9(2n+ 2m) - 3( 5n+5m+4nm ) - 6[ n(n-1)+ m(m-1) ] + 5(2n+ 2m)+3(4n + 2m) +3( 5n+5m+4nm ) + [8n(n-1) + 4m(m-1) ] </math>

The only question remaining is, what pair of values for <math>~(n, m)</math> result in both of these expressions going to zero? Simplifying the first expression gives,

<math>~0</math>

<math>~=</math>

<math>~6(2n+ 2m) + 2 ( 5n+5m+4nm ) + 4 [ n(n-1)+ m(m-1) ] - 5( 4n + 2m) -2( 5n+5m+4nm ) - [8n(n-1) +2m(m-1) ]</math>

 

<math>~=</math>

<math>~2m -8n -4n(n-1) + 2m(m-1) </math>

 

<math>~=</math>

<math>~ 2[m^2 - 2n(n+1)] \, . </math>

Simplifying the second expression gives,

<math>~0</math>

<math>~=</math>

<math>~-9(2n+ 2m) - 3( 5n+5m+4nm ) - 6[ n(n-1)+ m(m-1) ] + 5(2n+ 2m)+3(4n + 2m) +3( 5n+5m+4nm ) + [8n(n-1) + 4m(m-1) ] </math>

 

<math>~=</math>

<math>~4n - 2m + 2n(n-1) - 2m(m-1) </math>

 

<math>~=</math>

<math>~2[n(n+1) -m^2] </math>

Okay. Because it is not possible for both of these last two constraints to be simultaneously satisfied, I conclude that this last, specific eigenfunction guess is incorrect.


Second Guess

Let's try again, keeping the same values of the <math>~b_0</math> and <math>~b_2</math> — that is,

<math>~(b_0 + b_2x^2) = (2-x^2)</math>      <math>~\Rightarrow</math>      <math>~b_0 = 2</math>      and      <math>~b_2 = -1</math>

— but leaving the values of <math>~a_0</math> and <math>~a_2</math> unspecified. In this case, the LAWE becomes,

<math>~ 0 </math>

<math>~=</math>

<math>~ \biggl[ \alpha(a_0 + a_2x^2) (2 - x^2) + 2x^2(2n a_2 - m a_0) - 2x^4 (na_2 + m a_2) \biggr](5-3x^2)(a_0 + a_2x^2) -\sigma^2 (a_0 + a_2x^2)^2 (2 - x^2) </math>

 

 

<math>~ - \biggl[ n a_2(2 - x^2) -m (a_0 + a_2x^2) + 4(2n a_2 - m a_0) - 4(na_2 + m a_2)x^2 - 4n m a_2 x^2 \biggr](1-x^2)(2-x^2)(a_0 + a_2x^2) </math>

 

 

<math>~ - \biggl[2n(n-1) a_2^2(2 - x^2)^2 + 2m(m-1) (a_0 + a_2x^2)^2 \biggr]x^2 (1-x^2) </math>

 

<math>~=</math>

<math>~ \biggl\{ \alpha [2a_0 ] + x^2[(4n a_2 - 2m a_0) + \alpha (2a_2-a_0) ] - x^4 [(2na_2 + 2ma_2 ) + a_2\alpha ]\biggr\} [ 5a_0 + (5a_2-3a_0)x^2 -3a_2x^4] -\sigma^2 [ 2a_0 + (2a_2-a_0)x^2 - a_2x^4 ] (a_0+a_2x^2) </math>

 

 

<math>~ + \biggl[ ( 5m a_0 - 10n a_2) + (4n m a_2 + 5na_2 + 5m a_2)x^2 \biggr] (1-x^2)[ 2a_0 + (2a_2-a_0)x^2 - a_2x^4 ] </math>

 

 

<math>~ - \biggl\{ [ 8n(n-1) a_2^2 + 2m(m-1)a_0^2 ] + [ -8n(n-1) a_2^2 + 4m(m-1)a_0 a_2 ]x^2 + [ 2n(n-1) a_2^2 + 2m(m-1)a_2^2 ]x^4 \biggr\} x^2 (1-x^2) </math>

 

<math>~=</math>

<math>~ \biggl\{ \alpha [2a_0 ] + x^2[(4n a_2 - 2m a_0) + \alpha (2a_2-a_0) ] - x^4 [(2na_2 + 2ma_2 ) + a_2\alpha ]\biggr\} [ 5a_0 + (5a_2-3a_0)x^2 -3a_2x^4] </math>

 

 

<math>~ -\sigma^2\biggl\{ 2a_0^2 + [2a_0a_2 + a_0(2a_2-a_0)]x^2 +[a_2(2a_2-a_0) -a_0a_2]x^4 - a_2^2 x^6 \biggr\} </math>

 

 

<math>~ + \biggl\{ [ ( 5m a_0 - 10n a_2) ] + [(4n m a_2 + 5na_2 + 5m a_2)- ( 5m a_0 - 10n a_2) ]x^2 - [ 4n m a_2 + 5na_2 + 5m a_2 ]x^4 \biggr\} [ 2a_0 + (2a_2-a_0)x^2 - a_2x^4 ] </math>

 

 

<math>~ - \biggl\{ [ 8n(n-1) a_2^2 + 2m(m-1)a_0^2 ]x^2 + [ -8n(n-1) a_2^2 + 4m(m-1)a_0 a_2 ]x^4 + [ 2n(n-1) a_2^2 + 2m(m-1)a_2^2 ]x^6 \biggr\} (1-x^2) \, . </math>


So, the coefficients of each even power of <math>~x^n</math> are:

<math>~x^0</math>   :  

<math>~ 10a_0^2 \alpha - 2a_0^2\sigma^2 + 2a_0[ ( 5m a_0 - 10n a_2) ] </math>

<math>~x^2</math>   :  

<math>~5a_0[(4n a_2 - 2m a_0) + \alpha (2a_2-a_0) ] + \alpha [2a_0 ](5a_2-3a_0)-\sigma^2[2a_0a_2 + a_0(2a_2-a_0)]

</math>

<math>~+ 2a_0[(4n m a_2 + 5na_2 + 5m a_2)- ( 5m a_0 - 10n a_2) ] + (2a_2-a_0)[ ( 5m a_0 - 10n a_2) ] - [ 8n(n-1) a_2^2 + 2m(m-1)a_0^2 ] </math>

<math>~x^4</math>   :  

<math>~ - 5a_0[(2na_2 + 2ma_2 ) + a_2\alpha ] + (5a_2-3a_0)[(4n a_2 - 2m a_0) + \alpha (2a_2-a_0) ] - 6a_0 a_2\alpha -\sigma^2[a_2(2a_2-a_0) -a_0a_2]

</math>

<math>~- 2a_0[ 4n m a_2 + 5na_2 + 5m a_2 ] + (2a_2-a_0)[(4n m a_2 + 5na_2 + 5m a_2)- ( 5m a_0 - 10n a_2) ] - a_2(5ma_0 - 10na_2) </math>

<math>~-[ -8n(n-1) a_2^2 + 4m(m-1)a_0 a_2 ] + [ 8n(n-1) a_2^2 + 2m(m-1)a_0^2 ]

</math>

<math>~x^6</math>   :  

<math>~ -3a_2[(4n a_2 - 2m a_0) + \alpha (2a_2-a_0) ] - (5a_2-3a_0)[(2na_2 + 2ma_2 ) + a_2\alpha ] +\sigma^2 a_2^2

</math>

<math>~- (2a_2-a_0)[ 4n m a_2 + 5na_2 + 5m a_2 ] - a_2[(4n m a_2 + 5na_2 + 5m a_2)- ( 5m a_0 - 10n a_2) ] </math>

<math>~- [ 2n(n-1) a_2^2 + 2m(m-1)a_2^2 ] + [ -8n(n-1) a_2^2 + 4m(m-1)a_0 a_2 ]

</math>

<math>~x^8</math>   :  

<math>~ 3a_2 a_2\alpha + a_2[ 4n m a_2 + 11na_2 + 11m a_2 ] + [ 2n(n-1) a_2^2 + 2m(m-1)a_2^2 ] </math>


After simplification:

<math>~x^0</math>   :  

<math>~ 10a_0^2 \alpha - 2a_0^2\sigma^2 + 10m a_0^2 - 20n a_0a_2 </math>

<math>~x^2</math>   :  

<math>~\alpha (20a_0a_2-11a_0^2) -\sigma^2[4a_0a_2 -a_0^2]

</math>

<math>~+ 60na_0a_2 -20na_2^2 + 20m a_0a_2 -25m a_0^2 + 8n m a_0a_2 - [ 8n(n-1) a_2^2 + 2m(m-1)a_0^2 ] </math>

<math>~x^4</math>   :  

<math>~ \alpha (10a_2^2 - 22a_0a_2+3a_0^2) -\sigma^2 (2a_2^2 -2a_0a_2) -47n a_0a_2 + 60n a_2^2 - 50ma_0a_2 + 11m a_0^2 + 10m a_2^2-12n m a_0a_2 + 8n m a_2^2

</math>

<math>~+ 16n(n-1) a_2^2 - 4m(m-1)a_0 a_2 + 2m(m-1)a_0^2 </math>

<math>~x^6</math>   :  

<math>~ \alpha (-a_2^2) +\sigma^2 a_2^2 -47n a_2^2 + 11na_0a_2+ 22 m a_0a_2 -25ma_2^2 -12n m a_2^2 + 4n m a_0a_2

</math>

<math>~-10n(n-1) a_2^2 - 2m(m-1)a_2^2 + 4m(m-1)a_0 a_2 </math>

<math>~x^8</math>   :  

<math>~ \{ 3\alpha + [ 4n m + 11n + 11m ] + [ 2n(n-1) + 2m(m-1) ]\}a_2^2 </math>



Third Guess

Let's try again, keeping the same values of the <math>~b_0</math> and <math>~b_2</math> — that is,

<math>~(b_0 + b_2x^2) = (2-x^2)</math>      <math>~\Rightarrow</math>      <math>~b_0 = 2</math>      and      <math>~b_2 = -1</math>

— but leaving the values of <math>~a_0</math> and <math>~a_2</math> unspecified. In this case, the LAWE becomes,

<math>~ 0 </math>

<math>~=</math>

<math>~ \biggl[ \alpha(a_0 + a_2x^2) (2 - x^2) + 2x^2(2n a_2 - m a_0) - 2a_2 x^4 (n + m ) \biggr](5-3x^2) -\sigma^2 (a_0 + a_2x^2) (2 - x^2) </math>

 

 

<math>~ + \biggl[- n a_2(2 - x^2) + m (a_0 + a_2x^2) - 4(n a_2 2 - m a_0) + 4(na_2 + m a_2)x^2 + 4n m a_2 x^2 \biggr](1-x^2)(2-x^2) </math>

 

 

<math>~ - \frac{1}{ (a_0 + a_2x^2)}\biggl[2n(n-1) a_2^2(2 - x^2)^2 + 2m(m-1) (a_0 + a_2x^2)^2 \biggr](1-x^2)x^2 </math>

 

<math>~=</math>

<math>~\biggl\{[2a_0\alpha ] + [\alpha(- a_0 +2a_2) + (4n a_2 - 2m a_0) ]x^2 + [-a_2\alpha - 2a_2 (n + m ) ]x^4 \biggr\} (5-3x^2) -\sigma^2 [2a_0 + (-a_0 + 2a_2)x^2 -a_2x^4] </math>

 

 

<math>~ + \biggl[( -10na_2 + 5ma_0 ) + (5ma_2 + 5na_2 + 4n m a_2 )x^2 \biggr](2-3x^2 +x^4) </math>

 

 

<math>~ - \frac{1}{ (a_0 + a_2x^2)}\{ [ 8n(n-1)a_2^2 +2m(m-1)a_0^2 ] + [ -8n(n-1)a_2^2 +4m(m-1)a_0 a_2 ]x^2 + [ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ]x^4 \} (x^2-x^4) </math>

 

<math>~=</math>

<math>~\biggl\{ [10a_0\alpha - 2a_0\sigma^2 ] + [5\alpha(- a_0 +2a_2) + 5(4n a_2 - 2m a_0) -6a_0\alpha + (a_0 - 2a_2)\sigma^2 ]x^2 </math>

 

 

<math>~ + [-5a_2\alpha - 10a_2 (n + m ) -3\alpha(- a_0 +2a_2) -3 (4n a_2 - 2m a_0) + a_2\sigma^2]x^4 + [3a_2\alpha +6a_2 (n + m ) ]x^6 \biggr\} (a_0 + a_2x^2) </math>

 

 

<math>~ +\biggl\{ [2( -10na_2 + 5ma_0 )] + [2 (5ma_2 + 5na_2 + 4n m a_2 ) -3( -10na_2 + 5ma_0 )]x^2 </math>

 

 

<math>~ + [ -3 (5ma_2 + 5na_2 + 4n m a_2 ) + ( -10na_2 + 5ma_0 )]x^4 + (5ma_2 + 5na_2 + 4n m a_2 )x^6 \biggr\} (a_0 + a_2x^2) </math>

 

 

<math>~ - \{ [ 8n(n-1)a_2^2 +2m(m-1)a_0^2 ] + [ -8n(n-1)a_2^2 +4m(m-1)a_0 a_2 ]x^2 + [ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ]x^4 \} (x^2-x^4) </math>

 

<math>~=</math>

<math>~\biggl\{ [10a_0\alpha - 2a_0\sigma^2 ] + [\alpha(- 11a_0 +10a_2) + (20n a_2 - 10m a_0) + (a_0 - 2a_2)\sigma^2 ]x^2 </math>

 

 

<math>~ + [\alpha(3a_0 - 11a_2) + (-22n a_2 +6m a_0 - 10a_2 m) + a_2\sigma^2]x^4 + [3a_2\alpha +6a_2 (n + m ) ]x^6 \biggr\} (a_0 + a_2x^2) </math>

 

 

<math>~ + \biggl\{ [ -20na_2 + 10ma_0 ] + [10ma_2 + 40na_2 + 8n m a_2 -15ma_0 ]x^2 </math>

 

 

<math>~ + [ -15ma_2 - 25na_2 -12 n m a_2 + 5ma_0 ]x^4 + [5ma_2 + 5na_2 + 4n m a_2 ]x^6 \biggr\} (a_0 + a_2x^2) </math>

 

 

<math>~ - \{ [ 8n(n-1)a_2^2 +2m(m-1)a_0^2 ]x^2 + [ -8n(n-1)a_2^2 +4m(m-1)a_0 a_2 ]x^4 + [ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ]x^6 \} (1-x^2) </math>

So, the coefficients of each even power of <math>~x^n</math> are:

<math>~x^0</math>   :  

<math>~a_0[10a_0\alpha - 2a_0\sigma^2 ] + a_0 [ -20na_2 + 10ma_0 ]</math>

<math>~x^2</math>   :  

<math>~a_0 [\alpha(- 11a_0 +10a_2) + (20n a_2 - 10m a_0) + (a_0 - 2a_2)\sigma^2 ] + a_2[10a_0\alpha - 2a_0\sigma^2 ] + a_0[10ma_2 + 40na_2 + 8n m a_2 -15ma_0 ]

</math>

<math> + a_2[ -20na_2 + 10ma_0 ] - [ 8n(n-1)a_2^2 +2m(m-1)a_0^2 ] </math>

<math>~x^4</math>   :  

<math>~ a_0 [\alpha(3a_0 - 11a_2) + (-22n a_2 +6m a_0 - 10a_2 m) + a_2\sigma^2] + a_2[\alpha(- 11a_0 +10a_2) + (20n a_2 - 10m a_0) + (a_0 - 2a_2)\sigma^2 ]

</math>

<math>~ + a_0 [ -15ma_2 - 25na_2 -12 n m a_2 + 5ma_0 ] + a_2[10ma_2 + 40na_2 + 8n m a_2 -15ma_0 ] </math>

<math> ~ - [ -8n(n-1)a_2^2 +4m(m-1)a_0 a_2 ] + [ 8n(n-1)a_2^2 +2m(m-1)a_0^2 ]

</math>

<math>~x^6</math>   :  

<math>~ a_0[3a_2\alpha +6a_2 (n + m ) ] + a_2[\alpha(3a_0 - 11a_2) + (-22n a_2 +6m a_0 - 10a_2 m) + a_2\sigma^2]

</math>

<math>~ a_0[5ma_2 + 5na_2 + 4n m a_2 ] + a_2[ -15ma_2 - 25na_2 -12 n m a_2 + 5ma_0 ] </math>

<math>~ -[ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ] + [ -8n(n-1)a_2^2 +4m(m-1)a_0 a_2 ]

</math>

<math>~x^8</math>   :  

<math>~ a_2[3a_2\alpha +6a_2 (n + m ) ] + a_2[5ma_2 + 5na_2 + 4n m a_2 ] + [ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ] </math>


After simplification:

<math>~x^0</math>   :  

<math>~ \alpha(10a_0^2) + \sigma^2(- 2a_0^2) -20n a_0a_2 + 10ma_0^2 </math>

<math>~x^2</math>   :  

<math>~\alpha(- 11a_0^2 +20 a_0a_2) + \sigma^2(a_0^2 - 4 a_0a_2) + 60n a_0a_2 -20na_2^2- 25m a_0^2 + 20m a_0a_2 + 8n m a_0a_2

</math>

<math> - [ 8n(n-1)a_2^2 + 2m(m-1)a_0^2 ] </math>

<math>~x^4</math>   :  

<math>~ \alpha(10a_2^2 - 22 a_0a_2 +3a_0^2) - \sigma^2(2a_2^2- 2a_0a_2 ) - 27n a_0a_2+ 60n a_2^2 - 50m a_0a_2 +11m a_0^2+ 10ma_2^2-12 n m a_0a_2 + 8n m a_2^2

</math>

<math> ~ + 16n(n-1)a_2^2 - 4m(m-1)a_0 a_2 + 2m(m-1)a_0^2 </math>

<math>~x^6</math>   :  

<math>~ a_0[3a_2\alpha +6a_2 (n + m ) ] + a_2[\alpha(3a_0 - 11a_2) + (-22n a_2 +6m a_0 - 10a_2 m) + a_2\sigma^2]

</math>

<math>~ a_0[5ma_2 + 5na_2 + 4n m a_2 ] + a_2[ -15ma_2 - 25na_2 -12 n m a_2 + 5ma_0 ] </math>

<math>~ -[ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ] + [ -8n(n-1)a_2^2 +4m(m-1)a_0 a_2 ]

</math>

<math>~x^8</math>   :  

<math>~ a_2[3a_2\alpha +6a_2 (n + m ) ] + a_2[5ma_2 + 5na_2 + 4n m a_2 ] + [ 2n(n-1)a_2^2 +2m(m-1)a_2^2 ] </math>


Whitworth's (1981) Isothermal Free-Energy Surface

© 2014 - 2021 by Joel E. Tohline
|   H_Book Home   |   YouTube   |
Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS |
Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation