Difference between revisions of "User:Tohline/Appendix/Ramblings/Additional Analytically Specified Eigenvectors for Zero-Zero Bipolytropes"

From VistrailsWiki
Jump to navigation Jump to search
(→‎Solve Cubic Equation: Insert montage for model21)
 
(49 intermediate revisions by the same user not shown)
Line 7: Line 7:
In our [[User:Tohline/SSC/Stability/BiPolytrope0_0#Radial_Oscillations_of_a_Zero-Zero_Bipolytrope|accompanying summary]], we have demonstrated how analytically specified eigenvectors can be constructed for the mode labeled, <math>~(\ell, j) = (2,1)</math>.  This was done by specifying <math>~\gamma_e</math>, then solving a quartic equation for <math>~q</math>.  Shortly after completing this summary chapter, we noticed that an alternate approach may be to specify <math>~q</math>, then solve for <math>~\gamma_e</math>; and this path may be simpler because it may only involve solution of a quadratic equation. (Actually, we later have realized that the relevant equation is cubic, rather than quadratic.  This is nevertheless simpler than the quartic equation.)  If this proves to be the case, then it may also be possible to analytically construct eigenvectors of additional modes.  Let's see.
In our [[User:Tohline/SSC/Stability/BiPolytrope0_0#Radial_Oscillations_of_a_Zero-Zero_Bipolytrope|accompanying summary]], we have demonstrated how analytically specified eigenvectors can be constructed for the mode labeled, <math>~(\ell, j) = (2,1)</math>.  This was done by specifying <math>~\gamma_e</math>, then solving a quartic equation for <math>~q</math>.  Shortly after completing this summary chapter, we noticed that an alternate approach may be to specify <math>~q</math>, then solve for <math>~\gamma_e</math>; and this path may be simpler because it may only involve solution of a quadratic equation. (Actually, we later have realized that the relevant equation is cubic, rather than quadratic.  This is nevertheless simpler than the quartic equation.)  If this proves to be the case, then it may also be possible to analytically construct eigenvectors of additional modes.  Let's see.


==Seek Alternate Solution to Mode21 <math>~(\ell,j) = (2,1)</math>==


===Setup21===
<div align="center">
According to [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP4|STEP 4 in our accompanying summary discussion]], we need to solve the following "derivative matching" expression:
<table border="1" cellpadding="8" align="center" width="70%">
<tr><td align="left">
<font color="red">'''CAUTION!'''</font>  &nbsp; All of the "core" eigenfunction expressions are suspect because (I think!) they have been developed herein under the assumption that <math>~g^{-2} = (1+2q^3)^2</math>, whereas the relation should read,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 16: Line 17:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{14(1+2q^3)^2}{7(1+2q^3)^2 - 5}</math>
<math>~\frac{1}{g^2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 22: Line 23:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{(1+2q^3)^2}{1+8q^3} = \frac{(1+2\Chi)^2}{1+8\Chi }\, .</math>
\frac{c_0 + (c_0 + 3)A_{21}q^3 + (c_0 + 6)A_{21}B_{21} q^6}{1 + A_{21}q^3 + A_{21}B_{21}q^6} \, ,
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
</td></tr>
</table>
</div>
==Seek Alternate Solution to Mode21 <math>~(\ell,j) = (2,1)</math>==


where, recognizing that, <math>~\alpha_e = c_0(c_0+2) \, ,</math>
===Setup21===
According to [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP4|STEP 4 in our accompanying summary discussion]], we need to solve the following "derivative matching" expression:
<div align="center">
<div align="center">
<table border="0" cellpadding="3" align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~A_{21}</math>
<math>~
\frac{\gamma_c}{\gamma_e} \biggl[ \frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)} \biggr]
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{c_0(c_0+5) - (c_0 + 6)(c_0 + 11)}{(c_0 + 3)(c_0+5) - \alpha_e}\biggr] </math>
<math>~
c_0 + \frac{3A_{2}\Chi + 6A_{2}B_{2} \Chi^2}{1 + A_{2}\Chi + A_{2}B_{2}\Chi^2}
- 3\biggl(\frac{\gamma_c}{\gamma_e} -1\biggr)
</math>
   </td>
   </td>
</tr>
</tr>
Line 47: Line 60:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow~~~ \frac{\gamma_c}{\gamma_e} \biggl[\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)}\biggr] + 3\biggl(\frac{\gamma_c}{\gamma_e} -1\biggr) </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 53: Line 66:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{c_0^2 + 5c_0 - (c_0^2 + 17c_0 + 66)}{(c_0^2 + 8c_0 + 15) - (c_0^2+2c_0)}\biggr] </math>
<math>~
\frac{c_0 + (c_0 + 3)A_{2}\Chi + (c_0 + 6)A_{2}B_{2} \Chi^2}{1 + A_{2}\Chi + A_{2}B_{2}\Chi^2}  
</math>
   </td>
   </td>
</tr>
</tr>
 
</table>
</div>
where, recognizing that, <math>~\alpha_e = c_0(c_0+2) \, ,</math>
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
  <td align="right">
<math>~A_{2}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{c_0(c_0+5) - (c_0 + 6)(c_0 + 11)}{(c_0 + 3)(c_0+5) - \alpha_e}\biggr] </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{c_0^2 + 5c_0 - (c_0^2 + 17c_0 + 66)}{(c_0^2 + 8c_0 + 15) - (c_0^2+2c_0)}\biggr] </math>
  </td>
</tr>
 
<tr>
<tr>
   <td align="right">
   <td align="right">
Line 71: Line 114:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~B_{21}</math>
<math>~B_{2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 118: Line 161:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5} \, .</math>
<math>~\frac{\gamma_c}{\gamma_e} \biggl[\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)}\biggr] + 3\biggl(\frac{\gamma_c}{\gamma_e} -1\biggr)  \, .</math>
   </td>
   </td>
</tr>
</tr>
Line 637: Line 680:
<math>~c_0 = p - \mathfrak{r}^{1 / 3} \biggl[ \cos\biggl(\frac{\theta}{3}\biggr) - \sqrt{3} \sin\biggl(\frac{\theta}{3}\biggr) \biggr] \, .</math>
<math>~c_0 = p - \mathfrak{r}^{1 / 3} \biggl[ \cos\biggl(\frac{\theta}{3}\biggr) - \sqrt{3} \sin\biggl(\frac{\theta}{3}\biggr) \biggr] \, .</math>
</div>
</div>
In turn, for a given value of <math>~q</math>, the corresponding value of <math>~\alpha_e</math> is obtained via the relation, <math>~\alpha_e = c_0(c_0 + 2)</math>.
The right-hand panel of Figure 1 presents a plot of this <math>~c_0(q)</math> function; actually what has been plotted is the inverted relation, <math>~q(\alpha_e)</math>.  The open, red circular markers trace the portion of the function that provides physically viable solutions, in the sense that the corresponding value of <math>~\alpha_c</math> lies between the values, negative one and three; the filled, light=blue circular markers identify roots of the cubic equation that are not physically viable.


The light-blue diamonds in the right-hand panel of Figure 1 trace this <math>~c_0(q)</math> function, which is a solution to the above cubic equation.  In the left-hand panel, we re-display a plot that has been discussed in an [[User:Tohline/SSC/Stability/BiPolytrope0_0#Figure1|accompanying chapter]].  It contains a plot (blue markers) of the same function, but this time it is <math>~q(c_0)</math>, as determined from the root of a quartic equation.  In order to illustrate more clearly that the two curves are the same, we have plotted the quartic solution ''on top of'' the cubic solution in the right-hand panel.  The <math>~q(\alpha_e)</math> curve in the right-hand panel extends to smaller values of <math>~q</math> than does the corresponding curve in the left-hand panel.  The curve on the left has been purposely truncated because, in the context of its  [[User:Tohline/SSC/Stability/BiPolytrope0_0#Figure1|original presentation]], lower values of q would be associated with unrealistic values of the adiabatic exponent of the core.
In the left-hand panel of Figure 1, we re-display a plot that has been discussed in an [[User:Tohline/SSC/Stability/BiPolytrope0_0#Figure1|accompanying chapter]].  It contains a plot (blue markers) of the same same <math>~q(\alpha_e)</math> function but, this time, as determined from the root of a quartic equation.  In order to illustrate more clearly that the two curves are the same, we have plotted the quartic solution (small, purple circular markers) ''on top of'' the cubic solution in the right-hand panel.   


<div align="center">
<div align="center">
Line 645: Line 691:
</th></tr>
</th></tr>
<tr><td align="center">
<tr><td align="center">
[[File:QuarticSolution02.png|300px|quartic solution]]
[[File:Quartic21Solution02Corrected.png|300px|quartic solution]]
[[File:CubicSolution01.png|300px|quartic solution]]
[[File:Cubic21Solution01Corrected.png|300px|cubic solution]]
</td></tr>
</td></tr>
</table>
</table>
</div>
</div>


 
<!-- DELETE IMAGE (it appears immediately below)
<div align="center">
<div align="center">
[[File:Model21Montage2.png|750px|Montage of Stability Results for (ell,j) = (2,1) quantum numbers]]
[[File:Model21Montage2.png|750px|Montage of Stability Results for (ell,j) = (2,1) quantum numbers]]
</div>
</div>
-->


==Try Mode22 <math>~(\ell,j) = (2,2)</math>==
===Illustration21===
 
 
In this case we need to replace the core eigenfunction segment that was specified in <font color="red"><b>STEP 3</b></font> in our [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP3|separate discussion]] of mode  <math>~(\ell,j) = (2,1)</math> with the following:
<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center" colspan="2">
Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes<br />
 
<font color="red">Quantum Numbers:</font>&nbsp; <math>~(\ell,j) = (2,1)</math>
  </th>
</tr>
<tr>
  <th align="center" width="50%"> Analyzable Model Sequence</th>
  <th align="center" width="50%">One Example Eigenfunction</th>
</tr>
<tr>
  <td align="center" colspan="2">
[[File:Model21MontageCorrected.png|800px|Montage of Stability Results for (ell,j) = (2,1) quantum numbers]]
  </td>
</tr>
</table>
</div>
 
 
<font color="darkblue"><b>''Top-Left Panel:''</b></font>&nbsp; Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> &#8212; where <math>~\gamma_e</math> is the adiabatic exponent of the envelope &#8212; in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (2,1)</math>.  Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the ''core'' [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>.  The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
 
<font color="darkblue"><b>''Bottom-Left Panel:''</b></font> &nbsp; Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>.  Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top.  The core is ''more'' compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line.  The yellow circular marker identifies the same example model as it does in the top-left panel.
 
<font color="darkblue"><b>''Top-Right Panel:''</b></font>&nbsp; Displays &#8212; as a function of the fractional radius, <math>~r_0/R = q\xi</math> &#8212;  the analytically determined eigenfunction for the <math>~(\ell,j) = (2,1)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,
<div align="center">
<math>~\alpha_e =-0.35</math> &nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math> &nbsp; &nbsp; &nbsp; <math>~c_0~\mathrm{(plus)} = -0.1937742</math>
</div>
and, correspondingly (via our [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Quartic|separate solution of the governing ''quartic'' equation]]), <math>~(q,\alpha_c) = (0.6840119, +0.8326585)</math>.  (This same plot appears in [[User:Tohline/SSC/Stability/BiPolytrope0_0#Figure1|Figure 1 of our accompanying summary]].) Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction, namely,


<div align="center">
<div align="center">
Line 665: Line 742:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~x_{j=2} |_\mathrm{core}</math>
<math>~x_{j=1} |_\mathrm{core}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 672: Line 749:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{35 -  126 (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35 -  126 (1+2q^3)^2 + 99 (1+2q^3)^4 } \, .</math>
\frac{5 (1+8q^3) -  7 (1+2q^3)^2 \xi^2}{5 (1+8q^3)-7(1+2q^3)^2} \, ;</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
So, following the procedure outlined in <font color="red"><b>STEP 4</b></font> of our [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP3|separate discussion]], the requirement that the first derivatives match at the interface becomes,
and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction, namely,
 
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 684: Line 760:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{ -  252 (1+2q^3)^2 + 396 (1+2q^3)^4 }{35 -  126 (1+2q^3)^2 + 99 (1+2q^3)^4 }</math>
<math>~x_{\ell=2} |_\mathrm{env}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 691: Line 767:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{c_0 + (c_0 + 3)A_{21}q^3 + (c_0 + 6)A_{21}B_{21} q^6}{1 + A_{21}q^3 + A_{21}B_{21}q^6} \, .
\xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{21} \xi^{3} + q^6 A_{21}B_{21}\xi^{6} }{ 1 + q^3 A_{21} q^6 A_{21}B_{21} }\biggr] \, ,
</math>
</math>
   </td>
   </td>
Line 697: Line 773:
</table>
</table>
</div>
</div>
 
where the coefficients, <math>~A_{21}, B_{21}</math>, are as [[#Setup21|defined above]] in terms of the parameter, <math>~c_0</math>.  The corresponding eigenfrequency is, from the perspective of the core,
So, as described [[#Solve_Cubic_Equation|above]], in order to determine <math>~c_0</math> for any specified value of <math>~q</math> &#8212; or, equivalently, <math>~\Chi \equiv q^3</math> &#8212; we simply need to replace <math>~Q_{21}</math> with <math>~Q_{22}</math> in the definition of the coefficient of each term in the cubic expression, where,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 704: Line 779:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{22}</math>
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~20\gamma_c - 8 = \frac{8(7 + \alpha_c)}{3-\alpha_c} \approx 28.9115809\, ;</math>
\frac{ -  252 (1+2\Chi)^2  + 396 (1+2\Chi)^4 }{35 - 126 (1+2\Chi)^2 + 99 (1+2\Chi)^4 } \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
and, from the perspective of the envelope,
<div align="center">
<table border="0" cellpadding="5" align="center">


<div align="center">
<tr>
<table border="1" cellpadding="8">
  <td align="right">
<tr><td align="center" colspan="3">Figure 2:  Results for mode <math>~(\ell,j) = (2,2)</math></td></tr>
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
<tr><td align="center" rowspan="9">
  </td>
[[File:Cubic22Combo2.png|850px|Results for (ell,j) = (2,2)]]
  <td align="center">
</td>
<math>~=</math>
<td align="right"><math>~q</math></td>
  </td>
<td align="center"><math>~0.768375</math></td>
  <td align="left">
<math>~3\gamma_e\biggl[\alpha_e + 5c_0 + 22\biggr] \frac{\rho_e}{\rho_c} \approx 28.9115807\, ,</math>
  </td>
</tr>
</tr>
</table>
</div>
where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.3902664</math>.
<font color="darkblue"><b>''Bottom-Right Panel:''</b></font>&nbsp; The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale.  Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
<td align="right"><math>~c_0</math></td>
  <td align="right">
<td align="center"><math>~-0.16913291</math></td>
<math>~y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>
  </td>
</tr>
</tr>
<tr>
</table>
<td align="right"><math>~\gamma_e</math></td>
</div>
<td align="center"><math>~1.208583401</math></td>
where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the outer edge.  In this type of log-amplitude plot, the eigenfunction's various ''nodes'' &#8212; that is, radial locations where <math>~x</math> passes through zero &#8212; are highlighted; here, specifically, there is one node inside the core and one node resides in the envelope.  Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope.  Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.684)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.39)</math>.
</tr>
 
<tr>
NOTE:  As may be ascertained from our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes]], equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to ''any'' value, <math>~0 \le \rho_e/\rho_c \le 1</math>.  We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is [[User:Tohline/SSC/Stability/BiPolytrope0_0#KeyConstraint|limited by the constraint]],
<td align="right"><math>~\alpha_e</math></td>
<td align="center"><math>~-0.30965989</math></td>
</tr>
<tr>
<td align="right"><math>~n_e</math></td>
<td align="center"><math>~4.7942454</math></td>
</tr>
<tr>
<td align="right"><math>~\nu</math></td>
<td align="center"><math>~0.6357658</math></td>
</tr>
<tr>
<td align="right"><math>~\rho_e/\rho_c</math></td>
<td align="center"><math>~0.4756979</math></td>
</tr>
<tr>
<td align="right"><math>~\gamma_c</math></td>
<td align="center"><math>~1.0464786</math></td>
</tr>
<tr>
<td align="right"><math>~\alpha_c</math></td>
<td align="center"><math>~-0.822343</math></td>
</tr>
 
</table>
</div>
 
 
Finally, as explained in our [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Allow_Different_Adiabatic_Exponents|summary discussion]], in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 770: Line 834:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~ \gamma_c </math>
<math>~\frac{\rho_e}{\rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 776: Line 840:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{[ 6 + 2j(2j+5)] } \biggl\{ 8 +
<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>
\gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, .
</math>
   </td>
   </td>
</tr>
</tr>
Line 784: Line 846:
</div>
</div>


==Try Mode31 <math>~(\ell,j) = (3,1)</math>==
The black-dotted curve in the ''Bottom-Right Panel'' displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model &#8212; such as the example illustrated here &#8212; that has an analytically prescribable radial-oscillation eigenvector.
===Setup31===
 
==Try Mode22 <math>~(\ell,j) = (2,2)</math>==
 
===Setup22===
 
In this case we need to replace the core eigenfunction segment that was specified in <font color="red"><b>STEP 3</b></font> in our [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP3|separate discussion]] of mode  <math>~(\ell,j) = (2,1)</math> with the following:


In this case we need to replace the envelope eigenfunction segment that was specified in <font color="red"><b>STEP 2</b></font> in our [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP2|separate discussion]] of mode  <math>~(\ell,j) = (2,1)</math> with the following:
For the case of <math>~\ell=2</math>, this means that, throughout the envelope, the eigenfunction is,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 794: Line 859:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~x_{\ell=3} |_\mathrm{env}</math>
<math>~x_{j=2} |_\mathrm{core}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 801: Line 866:
   <td align="left">
   <td align="left">
<math>~
<math>~
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{31} \xi^{3} +  q^6 A_{31}B_{31}\xi^{6} + q^9 A_{31}B_{31}C_{31}\xi^{9} }{ 1 +  q^3 A_{31}  + q^6 A_{31}B_{31} + q^9 A_{31}B_{31}C_{31} }\biggr] \, ,
\frac{35(1+8q^3)^2 - 126 (1+8q^3) (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } </math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
where, the values of the newly introduced coefficients,
 
<!-- ORIGINAL WRONG CORE EIGENFUNCTION
<div align="center">
<div align="center">
<table border="0" cellpadding="3" align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~A_{31}</math>
<math>~x_{j=2} |_\mathrm{core}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{c_0(c_0+5) - (c_0 + 9)(c_0 + 14)}{(c_0 + 3)(c_0+5) - \alpha_e}\biggr] </math>
<math>~
\frac{35 -  126 (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35 -  126 (1+2q^3)^2 + 99 (1+2q^3)^4 } \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
ORIGINAL WRONG CORE EIGENFUNCTION -->
So, following the procedure outlined in <font color="red"><b>STEP 4</b></font> of our [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP3|separate discussion]], the requirement that the first derivatives match at the interface becomes,


<tr>
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\frac{ -  252 (1+8q^3) (1+2q^3)^2  + 396 (1+2q^3)^4 }{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 }</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 830: Line 905:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{c_0^2 + 5c_0 - (c_0^2 + 23c_0 + 126)}{(c_0^2 + 8c_0 + 15) - (c_0^2 + 2c_0)}\biggr] </math>
<math>~
\frac{c_0 + (c_0 + 3)A_{22}q^3 + (c_0 + 6)A_{22}B_{22} q^6}{1 + A_{22}q^3 + A_{22}B_{22}q^6} \, ,
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


where,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~A_{22}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 842: Line 924:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>
<math>~A_{21} = -\biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr) \, ,</math>
   </td>
   </td>
</tr>
</tr>
Line 848: Line 930:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~B_{31}</math>
<math>~B_{22}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{(c_0+3)(c_0+8) - (c_0 + 9)(c_0 + 14)}{(c_0 + 6)(c_0+8) - \alpha_e}\biggr] </math>
<math>~B_{21} = -\biggl( \frac{c_0 + 7 }{2c_0+8}\biggr) \, . </math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
So, as described [[#Solve_Cubic_Equation|above]], in order to determine <math>~c_0</math> for any specified value of <math>~q</math> &#8212; or, equivalently, <math>~\Chi \equiv q^3</math> &#8212; we simply need to replace <math>~Q_{21}</math> with <math>~Q_{22}</math> in the definition of the coefficient of each term in the cubic expression, where,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~Q_{22}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{(c_0^2 + 11c_0 + 24) - (c_0^2 + 23c_0 + 126)}{(c_0^2 + 14c_0 +48) - (c_0^2 + 2c_0) }\biggr] </math>
<math>~
\frac{ -  252 (1+8q^3) (1+2q^3)^2 + 396 (1+2q^3)^4 }{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 }\, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


<!-- REMOVE FIGURE2, WHICH WAS CONSTRUCTED FROM INCORRECT EXPRESSION FOR g^2
<div align="center">
<table border="1" cellpadding="8">
<tr><td align="center" colspan="3">Figure 2:  Results for mode <math>~(\ell,j) = (2,2)</math></td></tr>
<tr><td align="center" rowspan="9">
[[File:Cubic22Combo2.png|850px|Results for (ell,j) = (2,2)]]
</td>
<td align="right"><math>~q</math></td>
<td align="center"><math>~0.768375</math></td>
</tr>
<tr>
<tr>
  <td align="right">
<td align="right"><math>~c_0</math></td>
&nbsp;
<td align="center"><math>~-0.16913291</math></td>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
  </td>
</tr>
</tr>
<tr>
<tr>
  <td align="right">
<td align="right"><math>~\gamma_e</math></td>
<math>~C_{31}</math>
<td align="center"><math>~1.208583401</math></td>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{(c_0+6)(c_0+11) - (c_0 + 9)(c_0 + 14)}{(c_0 + 9)(c_0+14) - \alpha_e}\biggr] </math>
  </td>
</tr>
</tr>
 
<tr>
<td align="right"><math>~\alpha_e</math></td>
<td align="center"><math>~-0.30965989</math></td>
</tr>
<tr>
<td align="right"><math>~n_e</math></td>
<td align="center"><math>~4.7942454</math></td>
</tr>
<tr>
<td align="right"><math>~\nu</math></td>
<td align="center"><math>~0.6357658</math></td>
</tr>
<tr>
<td align="right"><math>~\rho_e/\rho_c</math></td>
<td align="center"><math>~0.4756979</math></td>
</tr>
<tr>
<td align="right"><math>~\gamma_c</math></td>
<td align="center"><math>~1.0464786</math></td>
</tr>
<tr>
<td align="right"><math>~\alpha_c</math></td>
<td align="center"><math>~-0.822343</math></td>
</tr>
 
</table>
</div>
END REMOVAL OF FIGURE2 -->
 
Finally, as explained in our [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Allow_Different_Adiabatic_Exponents|summary discussion]], in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~ \gamma_c </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 902: Line 1,021:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{(c_0^2 +17c_0 + 66) - (c_0^2 + 23 c_0 + 126)}{(c_0^2 + 23 c_0 + 126) (c_0^2 + 2c_0) }\biggr] </math>
<math>~\frac{1}{[ 6 + 2j(2j+5)] } \biggl\{ 8 +
\gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


<!-- DELETE because image has been moved immediately below
<div align="center">
[[File:Model22Montage2.png|750px|Montage of Stability Results for (ell,j) = (2,2) quantum numbers]]
</div>
-->
===Illustration22===
<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center" colspan="2">
Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes<br />
<font color="red">Quantum Numbers:</font>&nbsp; <math>~(\ell,j) = (2,2)</math>
  </th>
</tr>
<tr>
  <th align="center" width="50%"> Analyzable Model Sequence</th>
  <th align="center" width="50%">One Example Eigenfunction</th>
</tr>
<tr>
<tr>
  <td align="right">
   <td align="center" colspan="2">
&nbsp;
[[File:Model22MontageCorrected.png|800px|Montage of Stability Results for (ell,j) = (2,1) quantum numbers]]
  </td>
   <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6) } \, .</math>
   </td>
   </td>
</tr>
</tr>
Line 920: Line 1,058:
</div>
</div>


Then, after defining,
 
<font color="darkblue"><b>''Top-Left Panel:''</b></font>&nbsp; Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> &#8212; where <math>~\gamma_e</math> is the adiabatic exponent of the envelope &#8212; in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (2,2)</math>.  Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the ''core'' [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>.  The red marker that is farthest to the left identifies the model whose analytically determined eigenfunction is displayed, as an example, in the right-hand panels.
 
<font color="darkblue"><b>''Bottom-Left Panel:''</b></font> &nbsp; Plotted points (only 2, here!) show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>.  Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top.  The core is ''more'' compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line. 
 
<font color="darkblue"><b>''Top-Right Panel:''</b></font>&nbsp; Displays &#8212; as a function of the fractional radius, <math>~r_0/R = q\xi</math> &#8212;  the analytically determined eigenfunction for the <math>~(\ell,j) = (2,2)</math> mode in the model identified by the red circular marker that is farthest to the left in both left-hand panels, for which,
<div align="center">
<math>~q = \biggl[0.005 + 179\biggl( \frac{0.985}{200} \biggr) \biggr] = 0.886575</math>
</div>
and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (-2.332785, +0.7763158, -0.9146699)</math>.  Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction, namely,
 
<div align="center">
<div align="center">
<table border="0" cellpadding="3" align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{31}</math>
<math>~x_{j=2} |_\mathrm{core}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5} \, .</math>
<math>~
\frac{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } \, ;</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
 
and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction, namely,
the matching relation becomes,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 944: Line 1,092:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{31}</math>
<math>~x_{\ell=2} |_\mathrm{env}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 951: Line 1,099:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{c_0 + (c_0 + 3)A_{31}\Chi + (c_0 + 6)A_{31}B_{31} \Chi^2 + (c_0 + 9)A_{31}B_{31} C_{31}\Chi^3 }{1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 }
\xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{22} \xi^{3} + q^6 A_{22}B_{22}\xi^{6} }{ 1 + q^3 A_{22} + q^6 A_{22}B_{22} }\biggr] \, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
where the coefficients, <math>~A_{22}, B_{22}</math>, are as [[#Setup22|defined above]] in terms of the parameter, <math>~c_0</math>.  The corresponding eigenfrequency is, from the perspective of the core,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow~~~ 0</math>
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 964: Line 1,117:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~42\gamma_c - 8 = \frac{8(18 + \alpha_c)}{3-\alpha_c} \approx 34.915496 \, ;</math>
[c_0 + (c_0 + 3)A_{31}\Chi + (c_0 + 6)A_{31}B_{31} \Chi^2 + (c_0 + 9)A_{31}B_{31} C_{31}\Chi^3 ] - Q_{31}[1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 ]
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
and, from the perspective of the envelope,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 978: Line 1,134:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~3\gamma_e\biggl[\alpha_e + 5c_0 + 22\biggr] \frac{\rho_e}{\rho_c} \approx 34.915496 \, ,</math>
\biggl[c_0 - (c_0 + 3)\cdot \frac{6(c_0 + 7)}{2c_0 + 5 }  \cdot \Chi
+ (c_0 + 6) \cdot \frac{ 6(c_0 + 7)}{2c_0 + 5 }  \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \Chi^2
- (c_0 + 9)\cdot \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \frac{ 2( c_0 + 10)}{7( c_0 + 6)  } \cdot\Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.5822407</math>.
<font color="darkblue"><b>''Bottom-Right Panel:''</b></font>&nbsp; The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale.  Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~y</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>
- Q_{31} \biggl[1 - \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \Chi + \frac{ 6(c_0 + 7)}{2c_0 + 5 }  \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \Chi^2  
- \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \frac{ 2( c_0 + 10)}{7( c_0 + 6)  } \cdot \Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the center of the configuration.  In this type of log-amplitude plot, the eigenfunction's various ''nodes'' &#8212; that is, radial locations where <math>~x</math> passes through zero &#8212; are highlighted; here, specifically, there are two nodes inside the core and none in the envelope, although one of the nodes in the core lies just inside the core/envelope interface.  Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope.  Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.89)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.58)</math>.
NOTE:  As may be ascertained from our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes]], equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to ''any'' value, <math>~0 \le \rho_e/\rho_c \le 1</math>.  We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models.  Instead, our ability to derive analytically prescribed eigenvectors is [[User:Tohline/SSC/Stability/BiPolytrope0_0#KeyConstraint|limited by the constraint]],
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\frac{\rho_e}{\rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,009: Line 1,172:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>
\biggl[7c_0(2c_0 + 5)(c_0 +4)( c_0 + 6) - (c_0 + 3)(c_0 +4)7( c_0 + 6)\cdot 6(c_0 + 7) \cdot \Chi
+ (c_0 + 6) 7( c_0 + 6)\cdot  (c_0 + 7) \cdot  (6c_0 + 51) \cdot \Chi^2
- (c_0 + 9)\cdot  (c_0 + 7) \cdot  (6c_0 + 51)\cdot 2( c_0 + 10) \cdot\Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
The black-dotted curve in the ''Bottom-Right Panel'' displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model &#8212; such as the example illustrated here &#8212; that has an analytically prescribable radial-oscillation eigenvector.
==Try Mode31 <math>~(\ell,j) = (3,1)</math>==
===Setup31===
In this case we need to replace the envelope eigenfunction segment that was specified in <font color="red"><b>STEP 2</b></font> in our [[User:Tohline/SSC/Stability/BiPolytrope0_0#STEP2|separate discussion]] of mode  <math>~(\ell,j) = (2,1)</math> with the following:
For the case of <math>~\ell=2</math>, this means that, throughout the envelope, the eigenfunction is,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~x_{\ell=3} |_\mathrm{env}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
- Q_{31} \biggl[ 7(2c_0 + 5)(c_0 +4)( c_0 + 6)
\xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{31} \xi^{3} + q^6 A_{31}B_{31}\xi^{6} q^9 A_{31}B_{31}C_{31}\xi^{9} }{ 1 + q^3 A_{31}  + q^6 A_{31}B_{31} + q^9 A_{31}B_{31}C_{31} }\biggr] \, ,
- 42(c_0 + 7)(c_0 +4)( c_0 + 6)  \cdot \Chi
+ 7( c_0 + 6)(c_0 + 7) \cdot (6c_0 + 51) \cdot \Chi^2
- 2(c_0 + 7) \cdot (6c_0 + 51) ( c_0 + 10) \cdot \Chi^3 \biggr]
</math>
</math>
   </td>
   </td>
</tr>
</tr>
 
</table>
</div>
where, the values of the newly introduced coefficients,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~A_{31}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{c_0(c_0+5) - (c_0 + 9)(c_0 + 14)}{(c_0 + 3)(c_0+5) - \alpha_e}\biggr] </math>
\biggl[7c_0(2c_0 + 5)(c_0^2 + 10c_0 + 24) - 42 (c_0 + 3)(c_0^2 + 10c_0 + 24)(c_0 + 7) \cdot \Chi
+ 7(c_0 + 6) ( c_0 + 6)\cdot (6c_0^2+93c_0 + 357)\cdot \Chi^2
- 2(c_0 + 9)\cdot (6c_0^2+93c_0 + 357)( c_0 + 10) \cdot\Chi^3 \biggr]  
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,052: Line 1,220:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
&nbsp;
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{c_0^2 + 5c_0 - (c_0^2 + 23c_0 + 126)}{(c_0^2 + 8c_0 + 15) - (c_0^2 + 2c_0)}\biggr] </math>
- Q_{31} \biggl[ 7(2c_0 + 5)(c_0^2 + 10c_0 + 24)
- 42(c_0 + 7)(c_0^2 + 10c_0 + 24)   \cdot \Chi
+ 7( c_0 + 6)(6c_0^2+93c_0 + 357) \cdot \Chi^2
- 2(6c_0^2+93c_0 + 357)( c_0 + 10) \cdot \Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow ~~~ 0</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,081: Line 1,238:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>
\biggl[7c_0 [ 2c_0^3 + 25c_0^2  + 98c_0 + 120 ] - 42 (c_0^2 + 10c_0 + 21) (c_0^2 + 10c_0 + 24) \cdot \Chi
+ 7(c_0^2 + 12c_0 + 36) (6c_0^2+93c_0 + 357)\cdot \Chi^2
- 2(c_0^2 + 19c_0 + 90) (6c_0^2+93c_0 + 357) \cdot\Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,091: Line 1,244:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
<math>~B_{31}</math>
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{(c_0+3)(c_0+8) - (c_0 + 9)(c_0 + 14)}{(c_0 + 6)(c_0+8) - \alpha_e}\biggr] </math>
- Q_{31} \biggl[ 7 [ 2c_0^3 + 25c_0^2  + 98c_0 + 120 ]
- 42(c_0^3 + 17c_0^2 + 94c_0  + 168)   \cdot \Chi
+ 7( 6c_0^3 + 129c_0^2  + 915 c_0 + 2142) \cdot \Chi^2
- 2(6c_0^3+ 153c_0^2 + 1287c_0 + 3570) \cdot \Chi^3 \biggr] \, .
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,115: Line 1,262:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{(c_0^2 + 11c_0 + 24) - (c_0^2 + 23c_0 + 126)}{(c_0^2 + 14c_0 +48) - (c_0^2 + 2c_0) }\biggr] </math>
\biggl[7c_0 [ 2c_0^3 + 25c_0^2 + 98c_0 + 120 ] - 42 (c_0^4 + 20c_0^3 + 145c_0^2 + 450c_0 +504) \cdot \Chi
+ 7(6c_0^4 + 165c_0^3 +1689c_0^2 + 7632 c_0  + 12852 )\cdot \Chi^2
</math><!--Used to be 7596 c_0-->
   </td>
   </td>
</tr>
</tr>
Line 1,127: Line 1,271:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
- 2(6c_0^4 + 207c_0^3 + 2664c_0^2 + 15153c_0 + 32130) \cdot\Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,138: Line 1,280:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
<math>~C_{31}</math>
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{(c_0+6)(c_0+11) - (c_0 + 9)(c_0 + 14)}{(c_0 + 9)(c_0+11) - \alpha_e}\biggr] </math>
- Q_{31} \biggl[ 7 [ 2c_0^3 + 25c_0^2  + 98c_0 + 120 ]
- 42(c_0^3 + 17c_0^2 + 94c_0  + 168)   \cdot \Chi
+ 7( 6c_0^3 + 129c_0^2  + 915 c_0 + 2142) \cdot \Chi^2
- 2(6c_0^3+ 153c_0^2 + 1287c_0 + 3570) \cdot \Chi^3 \biggr] \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


===If the Envelope Follows a 4/3 Adiabat===
====Case of c<sub>0</sub> (plus)====
In this case, <math>~(\gamma_e, \alpha_e, c_0) = (\tfrac{4}{3}, 0,0)</math>, so the three coefficients,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~A_{31}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,169: Line 1,298:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } = - \frac{42}{5} \, ,</math>
<math>~\biggl[ \frac{(c_0^2 +17c_0 + 66) - (c_0^2 + 23 c_0 + 126)}{(c_0^2 + 20 c_0 + 99) - (c_0^2 + 2c_0) }\biggr] </math>
   </td>
   </td>
</tr>
</tr>
Line 1,175: Line 1,304:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~B_{31}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,181: Line 1,310:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } = - \frac{17}{8} \, ,</math>
<math>~\frac{- 2( c_0 + 10)}{3(2 c_0 + 11) } \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


Then, after defining,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~C_{31}</math>
<math>~Q_{31}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6) } = - \frac{10}{21} \, .</math>
<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)} \, ,</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
Hence, for this ''specific'' case, the derivative-matching constraint becomes,


the matching relation becomes,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 1,205: Line 1,340:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5}</math>
<math>~Q_{31}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,219: Line 1,354:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow~~~ 0</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,226: Line 1,361:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{  3A_{31} \Chi + 6A_{31} B_{31} \Chi^2 + 9 A_{31}B_{31} C_{31}\Chi^3 }{1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 }
[c_0 + (c_0 + 3)A_{31}\Chi + (c_0 + 6)A_{31}B_{31} \Chi^2 + (c_0 + 9)A_{31}B_{31} C_{31}\Chi^3 ] - Q_{31}[1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 ]
</math>
</math>
   </td>
   </td>
Line 1,233: Line 1,368:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow ~~~ \mathrm{LHS}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,239: Line 1,374:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\mathrm{RHS}
<math>~
\biggl[c_0 - (c_0 + 3)\cdot \frac{6(c_0 + 7)}{2c_0 + 5 } \cdot \Chi
+ (c_0 + 6) \cdot \frac{ 6(c_0 + 7)}{2c_0 + 5 }  \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \Chi^2
- (c_0 + 9)\cdot \frac{ 6(c_0 + 7)}{2c_0 + 5 }  \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \frac{ 2( c_0 + 10)}{3( 2c_0 + 11)  } \cdot\Chi^3 \biggr]
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
where,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\mathrm{LHS}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ 14(1+2\Chi)^2 \biggr] \biggl[ 2^3\cdot 3\cdot 5 \cdot 7  
<math>~
- 2^3\cdot 3 \cdot 7\cdot 42 \Chi + 3 \cdot 7 \cdot 42 \cdot 17  \Chi^2  
- Q_{31} \biggl[1 - \frac{ 6(c_0 + 7)}{2c_0 + 5 }  \cdot \Chi + \frac{ 6(c_0 + 7)}{2c_0 + 5 }  \cdot \frac{ (6c_0 + 51)}{ 6(c_0 +4) } \cdot \Chi^2  
- 42 \cdot 17 \cdot 10\Chi^3 \biggr] \, ,</math>
- \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{  (6c_0 + 51)}{ 6(c_0 +4) } \cdot \frac{ 2( c_0 + 10)}{3( 2c_0 + 11)  } \cdot \Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,266: Line 1,399:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\mathrm{RHS}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ 7(1+2\Chi)^2 - 5\biggr]
<math>~
\biggl[ - 2^3\cdot 3^2 \cdot 7 \cdot 42 \Chi + 6 \cdot 3 \cdot 7\cdot 42 \cdot 17  \Chi^2  
\biggl[3c_0(2c_0 + 5)(c_0 +4)( 2c_0 + 11) - (c_0 + 3)(c_0 +4)3( 2c_0 + 11)\cdot 6(c_0 + 7) \cdot \Chi  
- 9 \cdot 42 \cdot 17 \cdot 10 \Chi^3 \biggr] \, .
+ (c_0 + 6) 3( 2c_0 + 11)\cdot (c_0 + 7) \cdot (6c_0 + 51) \cdot \Chi^2  
- (c_0 + 9)\cdot (c_0 + 7) \cdot  (6c_0 + 51)\cdot 2( c_0 + 10) \cdot\Chi^3 \biggr]  
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


Empirically, what appears to work is:  <math>~\Chi = 0.276837296 ~~~\Rightarrow ~~~~q = 0.651740735 \, .</math>
<tr>
 
====Case of c<sub>0</sub> (minus)====
In this case, <math>~(\gamma_e, \alpha_e, c_0) = (\tfrac{4}{3}, 0,-2)</math>, so the three coefficients,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
   <td align="right">
   <td align="right">
<math>~A_{31}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } = - 30 \, ,</math>
<math>~
- Q_{31} \biggl[ 3(2c_0 + 5)(c_0 +4)( 2c_0 + 11)
- 18(c_0 + 7)(c_0 +4)( 2c_0 + 11)  \cdot \Chi
+ 3( 2c_0 + 11)(c_0 + 7)  \cdot (6c_0 + 51) \cdot \Chi^2
- 2(c_0 + 7) \cdot (6c_0 + 51) ( c_0 + 10) \cdot \Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,301: Line 1,432:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~B_{31}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,307: Line 1,438:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } = - \frac{13}{4} \, ,</math>
<math>~
\biggl[(6c_0^2 + 15c_0) (2c_0^2 + 19c_0 + 44) - 18(c_0^2 + 10c_0 + 21) (2c_0^2 + 19c_0 + 44) \cdot \Chi
+ 3 (2c_0^2 + 23c_0 +66) (6c_0^2 + 93c_0 + 357) \cdot \Chi^2
- 2(c_0^2 + 19c_0 + 90) (6c_0^2 + 93c_0 + 357) \cdot\Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,313: Line 1,448:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~C_{31}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6} = - \frac{4}{7} \, .</math>
<math>~
- Q_{31} \biggl[ 3(2c_0 + 5) (2c_0^2 + 19c_0 + 44)
- 18(c_0 + 7) (2c_0^2 + 19c_0 + 44)  \cdot \Chi
+ 3( 2c_0 + 11) (6c_0^2 + 93c_0 + 357\cdot \Chi^2
- 2 (6c_0^2 + 93c_0 + 357)  ( c_0 + 10) \cdot \Chi^3 \biggr]
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
Hence, for this ''specific'' case, the derivative-matching constraint becomes,


<div align="center">
<div align="center">
Line 1,331: Line 1,470:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5}</math>
<math>~\Rightarrow~~~ 0</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,338: Line 1,477:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{c_0 + (c_0 + 3)A_{31}\Chi + (c_0 + 6)A_{31}B_{31} \Chi^2 + (c_0 + 9)A_{31}B_{31} C_{31}\Chi^3 }{1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 }
\biggl[(12c_0^4 + 144c_0^3 + 549c_0^2 + 660c_0) - 18(2c_0^4 + 39c_0^3 + 276c_0^2 + 839c_0 + 924) \cdot \Chi  
+ 3 (12c_0^4 + 324c_0^3 + 3249c_0^2 + 14349c_0 + 23562) \cdot \Chi^2  
- 2(6c_0^4 + 207c_0^3 + 2664c_0^2 + 15153c_0 + 32130) \cdot\Chi^3 \biggr]
</math>
</math>
   </td>
   </td>
Line 1,348: Line 1,489:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{-2 + A_{31}\Chi + 4A_{31}B_{31} \Chi^2 + 7A_{31}B_{31} C_{31}\Chi^3 }{1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 }
- Q_{31} \biggl[ 3(4c_0^3 + 48c_0^2 + 183c_0 +220)
- 18(2c_0^3 + 33c_0^2 + 177c_0 + 308)  \cdot \Chi  
+ 3( 12c_0^3 + 252c_0^2 + 1737c_0 + 3927)  \cdot \Chi^2
- 2 (6c_0^3 + 153c_0^2 + 1287c_0 + 3570) \cdot \Chi^3 \biggr]
</math>
</math>
   </td>
   </td>
</tr>
</tr>


<tr>
</table>
   <td align="right">
</div>
<math>~\Rightarrow ~~~ \mathrm{LHS}</math>
 
   </td>
===Root of Quartic Equation===
   <td align="center">
 
To solve this equation analytically, we follow the [http://en.wikipedia.org/wiki/Quartic_function#Summary_of_Ferrari.27s_method Summary of Ferrari's method] that is presented in Wikipedia's discussion of the Quartic Function to identify the roots of an arbitrary quartic equation. 
 
<div align="center" id="Quartic">
<table border="1" cellpadding="8" width="85%">
<tr><td align="left">
 
First, we adopt the shorthand notation:
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~0</math>
   </td>
   <td align="center">
<math>~=</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\mathrm{RHS}
<math>~ac_0^4 + bc_0^3 + c c_0^2 +d c_0 +e \, ,</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
where, in our particular case,
where,


<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\mathrm{LHS}</math>
<math>~e</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,384: Line 1,538:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ 14(1+2\Chi)^2 \biggr] \biggl[ 2\cdot 7 - 2\cdot 7 \cdot 30\Chi + 7 \cdot 15\cdot 13 \Chi^2 -15 \cdot 13 \cdot 4 \Chi^3
<math>~ - 18\cdot 924\Chi +3\cdot 23562\Chi^2 -2\cdot 32130\Chi^3    -Q_{31}[ 3\cdot 220 - 18\cdot 308\Chi +3\cdot 3927\Chi^2 -2\cdot 3570\Chi^3  ]   \, ,</math>
  \biggr] \, ,</math>
   </td>
   </td>
</tr>
</tr>
Line 1,391: Line 1,544:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\mathrm{RHS}</math>
<math>~d</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,397: Line 1,550:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ 7(1+2\Chi)^2 - 5\biggr]
<math>~660 -18\cdot 839\Chi +3\cdot 14349\Chi^2 -2\cdot 15153\Chi^3  -Q_{31}[3\cdot 183 -18\cdot 177\Chi + 3\cdot 1737\Chi^2 -2\cdot 1287\Chi^3   ]   \, ,</math>
\biggl[ -2^2\cdot 7 (1+ 15 \Chi) + 4\cdot 7\cdot 15 \cdot 13 \Chi^2 - 7\cdot 15 \cdot 13 \cdot 4\Chi^3  
\biggr] \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


Empirically, what appears to work is:  <math>~\Chi = 0.063819021 ~~~\Rightarrow ~~~~q = 0.399622605  \, .</math>
<tr>
 
   <td align="right">
===Root of Quartic Equation===
<math>~c</math>
 
To solve this equation analytically, we follow the [http://en.wikipedia.org/wiki/Quartic_function#Summary_of_Ferrari.27s_method Summary of Ferrari's method] that is presented in Wikipedia's discussion of the Quartic Function to identify the roots of an arbitrary quartic equation. 
 
<div align="center" id="Quartic">
<table border="1" cellpadding="8" width="85%">
<tr><td align="left">
 
First, we adopt the shorthand notation:
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ac_0^4 + bc_0^3 + c c_0^2 +d c_0 +e \, ,</math>
  </td>
</tr>
</table>
where, in our particular case,
 
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~e</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,442: Line 1,562:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - 42\cdot 504\Chi +7\cdot 12852\Chi^2 -2\cdot 32130\Chi^3   -Q_{31}[ 7\cdot 120 - 42\cdot 168\Chi +7\cdot 2142\Chi^2 -2\cdot 3570\Chi^3  \, ,</math>
<math>~549 - 18\cdot 276\Chi +3\cdot 3249\Chi^2 -2\cdot 2664\Chi^3     -Q_{31}[ 3\cdot 48 - 18\cdot 33\Chi + 3\cdot 252\Chi^2 -2\cdot 153 \Chi^3     ] \, ,</math>
   </td>
   </td>
</tr>
</tr>
Line 1,448: Line 1,568:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~d</math>
<math>~b</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,454: Line 1,574:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~7\cdot 120 -42\cdot 450\Chi +7\cdot 7632\Chi^2 -2\cdot 15153\Chi^3   -Q_{31}[7\cdot 98 -42\cdot 94\Chi + 7\cdot 915\Chi^2 -2\cdot 1287\Chi^3   ]   \, ,</math>
<math>~144- 18\cdot 39\Chi +3\cdot 324\Chi^2 -2\cdot 207\Chi^3   - 12 Q_{31}[ 1 - 3\Chi +3\Chi^2 - \Chi^3   ]\, ,</math>
   </td>
   </td>
</tr>
</tr>
Line 1,460: Line 1,580:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~c</math>
<math>~a</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,466: Line 1,586:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~7\cdot 98 -42\cdot 145\Chi +7\cdot 1689\Chi^2 -2\cdot 2664\Chi^3     -Q_{31}[ 7\cdot 25 - 42\cdot 17\Chi + 7\cdot 129\Chi^2 -2\cdot 153 \Chi^3     ]  \, ,</math>
<math>~12( 1 - 3 \Chi +3\Chi^2 -\Chi^3)    \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>


<tr>
</td></tr>
  <td align="right">
<math>~b</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~7\cdot 25-42\cdot 20\Chi  +7\cdot 165\Chi^2 -2\cdot 207\Chi^3    -Q_{31}[ 7\cdot 2 - 42\Chi +7\cdot 6\Chi^2 -2\cdot 6\Chi^3    ]\, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~a</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~7\cdot 2 -42 \Chi +7\cdot 6\Chi^2 -2\cdot 6 \Chi^3    \, .</math>
  </td>
</tr>
</table>
</table>
</div>
<!-- OMITTED GAMMA 4/3 PARAGRAPH ...


It is prudent to check that this set of coefficients satisfies the quartic expression ''at least'' in the case of the [[#If_the_Envelope_Follows_a_4.2F3_Adiabat|two examples given above]] when <math>~\gamma_e = \tfrac{4}{3}</math>.  When [[#Case_of_c0_.28plus.29|c<sub>0</sub> (plus)]] is examined &#8212; that is, for <math>~c_0 = 0</math> &#8212; we should find that the coefficient, <math>~e</math>, by itself should be zero.  Indeed, it is zero to many significant digits when the empirically derived value of <math>~\Chi = 0.276837296 </math> is plugged into the expression for <math>~e</math>.  Similarly, the quartic expression is satisfied with <math>~c_0 = -2</math> if, as [[#Case_of_c0_.28minus.29|was derived empirically above when c<sub>0</sub> (minus) was examined]], if the value of <math>~\Chi = 0.063819021 </math> is used to determine the five separate coefficient values:  
It is prudent to check that this set of coefficients satisfies the quartic expression ''at least'' in the case of the [[#If_the_Envelope_Follows_a_4.2F3_Adiabat|two examples given above]] when <math>~\gamma_e = \tfrac{4}{3}</math>.  When [[#Case_of_c0_.28plus.29|c<sub>0</sub> (plus)]] is examined &#8212; that is, for <math>~c_0 = 0</math> &#8212; we should find that the coefficient, <math>~e</math>, by itself should be zero.  Indeed, it is zero to many significant digits when the empirically derived value of <math>~\Chi = 0.276837296 </math> is plugged into the expression for <math>~e</math>.  Similarly, the quartic expression is satisfied with <math>~c_0 = -2</math> if, as [[#Case_of_c0_.28minus.29|was derived empirically above when c<sub>0</sub> (minus) was examined]], if the value of <math>~\Chi = 0.063819021 </math> is used to determine the five separate coefficient values:  
<table border="0" align="center"><tr><td align="center"><math>~(a,b,c,d,e) = (+11.48754, +73.56560, - 262.9735, -2253.197, - 3049.777)\, .</math></td></tr></table>
<table border="0" align="center"><tr><td align="center"><math>~(a,b,c,d,e) = (+11.48754, +73.56560, - 262.9735, -2253.197, - 3049.777)\, .</math></td></tr></table>
 
END OMITTED PARAGRAPH -->
</td></tr>
</table>
</div>
 


<div align="center" id="Quartic2">
<div align="center" id="Quartic2">
Line 1,510: Line 1,607:
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
  <td align="center" colspan="4">&nbsp;</td>
  <td align="center">for &nbsp;&nbsp;<math>~~\Chi = 0.063819021</math></td>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
Line 1,523: Line 1,616:
   <td align="left">
   <td align="left">
<math>~ c^2 - 3bd + 12ae  \, ,</math>
<math>~ c^2 - 3bd + 12ae  \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~1.460152\times 10^{5}</math>
   </td>
   </td>
</tr>
</tr>
Line 1,540: Line 1,627:
   <td align="left">
   <td align="left">
<math>~ 2c^3 - 9bcd + 27b^2e + 27ad^2 - 72ace  \, ,</math>
<math>~ 2c^3 - 9bcd + 27b^2e + 27ad^2 - 72ace  \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~3.700597\times 10^{7}</math>
   </td>
   </td>
</tr>
</tr>
Line 1,557: Line 1,638:
   <td align="left">
   <td align="left">
<math>~ \frac{8ac - 3b^2}{8a^2}  \, ,</math>
<math>~ \frac{8ac - 3b^2}{8a^2}  \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~-3.827100\times 10^1</math>
   </td>
   </td>
</tr>
</tr>
Line 1,574: Line 1,649:
   <td align="left">
   <td align="left">
<math>~ \frac{b^3 - 4abc + 8a^2 d}{8a^3}  \, ,</math>
<math>~ \frac{b^3 - 4abc + 8a^2 d}{8a^3}  \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~-9.001430\times 10^1</math>
   </td>
   </td>
</tr>
</tr>
Line 1,591: Line 1,660:
   <td align="left">
   <td align="left">
<math>~ \frac{1}{2^{1 / 3}} \biggl[\Delta_1 + \sqrt{\Delta_1^2 - 4\Delta_0^3}  \biggr]^{1 / 3} \, ,</math>
<math>~ \frac{1}{2^{1 / 3}} \biggl[\Delta_1 + \sqrt{\Delta_1^2 - 4\Delta_0^3}  \biggr]^{1 / 3} \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
''complex''
   </td>
   </td>
</tr>
</tr>
Line 1,608: Line 1,671:
   <td align="left">
   <td align="left">
<math>~2\Delta_0^{1/2}\cos\biggl[ \frac{1}{3} \cos^{-1}\biggl( \frac{\Delta_1^2}{4\Delta_0^3}\biggr)^{1/2}\biggr] \, ,</math>
<math>~2\Delta_0^{1/2}\cos\biggl[ \frac{1}{3} \cos^{-1}\biggl( \frac{\Delta_1^2}{4\Delta_0^3}\biggr)^{1/2}\biggr] \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~7.006174\times 10^{2}</math>
   </td>
   </td>
</tr>
</tr>
Line 1,625: Line 1,682:
   <td align="left">
   <td align="left">
<math>~ \frac{1}{2} \biggl[ - \frac{2p}{3} + \frac{1}{3a}\biggl(\Kappa + \frac{\Delta_0}{\Kappa} \biggr)  \biggr]^{1 / 2} \, .</math>
<math>~ \frac{1}{2} \biggl[ - \frac{2p}{3} + \frac{1}{3a}\biggl(\Kappa + \frac{\Delta_0}{\Kappa} \biggr)  \biggr]^{1 / 2} \, .</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~3.385402</math>
   </td>
   </td>
</tr>
</tr>
Line 1,636: Line 1,687:
Then the four roots of the quartic equation are,
Then the four roots of the quartic equation are,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="center" colspan="4">&nbsp;</td>
  <td align="center">for &nbsp;&nbsp;<math>~~\Chi = 0.063819021</math></td>
</tr>


<tr>
<tr>
Line 1,651: Line 1,697:
   <td align="left">
   <td align="left">
<math>~-\frac{b}{4a} - S + \frac{1}{2}\biggl[ -4S^2 - 2p + \frac{\kappa}{S} \biggr]^{1 / 2}  \, ,</math>
<math>~-\frac{b}{4a} - S + \frac{1}{2}\biggl[ -4S^2 - 2p + \frac{\kappa}{S} \biggr]^{1 / 2}  \, ,</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~-3.972821</math>
   </td>
   </td>
</tr>
</tr>
Line 1,670: Line 1,710:
<math>~-\frac{b}{4a} - S - \frac{1}{2}\biggl[ -4S^2 - 2p + \frac{\kappa}{S} \biggr]^{1 / 2}  \, ,</math>
<math>~-\frac{b}{4a} - S - \frac{1}{2}\biggl[ -4S^2 - 2p + \frac{\kappa}{S} \biggr]^{1 / 2}  \, ,</math>
   </td>
   </td>
  <td align="right">
</tr>
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~-5.999955</math>
  </td>
</tr>


<tr>
<tr>
Line 1,687: Line 1,721:
   <td align="left">
   <td align="left">
<math>~-\frac{b}{4a} + S + \frac{1}{2}\biggl[ -4S^2 - 2p - \frac{\kappa}{S} \biggr]^{1 / 2}  \, .</math>
<math>~-\frac{b}{4a} + S + \frac{1}{2}\biggl[ -4S^2 - 2p - \frac{\kappa}{S} \biggr]^{1 / 2}  \, .</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~5.568830</math>
   </td>
   </td>
</tr>
</tr>
Line 1,704: Line 1,732:
   <td align="left">
   <td align="left">
<math>~-\frac{b}{4a} + S - \frac{1}{2}\biggl[ -4S^2 - 2p - \frac{\kappa}{S} \biggr]^{1 / 2}  \, .</math>
<math>~-\frac{b}{4a} + S - \frac{1}{2}\biggl[ -4S^2 - 2p - \frac{\kappa}{S} \biggr]^{1 / 2}  \, .</math>
  </td>
  <td align="right">
&nbsp; &nbsp; &nbsp; &hellip; &nbsp; &nbsp; &nbsp;
  </td>
  <td align="right">
<math>~-2.000000</math>
   </td>
   </td>
</tr>
</tr>
Line 1,916: Line 1,938:
</div>
</div>


==Try Mode32 <math>~(\ell,j) = (3,2)</math>==
===Illustration31===
 
<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center" colspan="2">
Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes<br />
 
<font color="red">Quantum Numbers:</font>&nbsp; <math>~(\ell,j) = (3,1)</math>
  </th>
</tr>
<tr>
  <th align="center" width="50%"> Analyzable Model Sequence</th>
  <th align="center" width="50%">One Example Eigenfunction</th>
</tr>
<tr>
  <td align="center" colspan="2">
[[File:Model31MontageCorrectedTwice.png|600px|Montage of Stability Results for (ell,j) = (3,1) quantum numbers]]
  </td>
</tr>
</table>
</div>
 
 
<font color="darkblue"><b>''Top-Left Panel:''</b></font>&nbsp; Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> &#8212; where <math>~\gamma_e</math> is the adiabatic exponent of the envelope &#8212; in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (3,1)</math>.  Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the ''core'' [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>.  The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
 
<font color="darkblue"><b>''Bottom-Left Panel:''</b></font> &nbsp; Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>.  Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top.  The core is ''more'' compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line.  The yellow circular marker identifies the same example model as it does in the top-left panel.


Here, we set the quantum numbers to <math>~(\ell,j) = (3,2)</math>.  In this case, [[#Try_Mode22|as when addessing Mode22, above]], the core eigenfunction segment should be,
<font color="darkblue"><b>''Top-Right Panel:''</b></font>&nbsp; Displays &#8212; as a function of the fractional radius, <math>~r_0/R = q\xi</math> &#8212;  the analytically determined eigenfunction for the <math>~(\ell,j) = (3,1)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,
<div align="center">
<math>~q = \biggl[0.01 + 40\biggl( \frac{0.98}{99} \biggr) \biggr] \approx 0.4059596</math>
</div>
and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (-1.7819827, -0.3885031, -0.9647648)</math>.  Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction, namely,


<div align="center">
<div align="center">
Line 1,925: Line 1,977:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~x_{j=2} |_\mathrm{core}</math>
<math>~x_{j=1} |_\mathrm{core}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,932: Line 1,984:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{35 -  126 (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35 -  126 (1+2q^3)^2 + 99 (1+2q^3)^4 } \, ,</math>
\frac{5(1+8q^3) -  7 (1+2q^3)^2 \xi^2}{5(1+8q^3)-7(1+2q^3)^2} \, ;</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
and, [[#Try_Mode31|as when addressing Mode31, above]], throughout the envelope, the eigenfunction should be,
and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction, namely,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 1,950: Line 2,002:
   <td align="left">
   <td align="left">
<math>~
<math>~
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{32} \xi^{3} +  q^6 A_{32}B_{32}\xi^{6} +  q^9 A_{32}B_{32}C_{32}\xi^{9} }{ 1 +  q^3 A_{32}  +  q^6 A_{32}B_{32} +  q^9 A_{32}B_{32}C_{32} }\biggr] \, ,
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{31} \xi^{3} +  q^6 A_{31}B_{31}\xi^{6} +  q^9 A_{31}B_{31}C_{31}\xi^{9} }{ 1 +  q^3 A_{31}  +  q^6 A_{31}B_{31} +  q^9 A_{31}B_{31}C_{31} }\biggr] \, ,
</math>
</math>
   </td>
   </td>
Line 1,956: Line 2,008:
</table>
</table>
</div>
</div>
where,
where the coefficients, <math>~A_{31}, B_{31}, C_{31}</math>, are as [[#Setup31|defined above]] in terms of the parameter, <math>~c_0</math>.  The corresponding eigenfrequency is, from the perspective of the core,
<div align="center">
<div align="center">
<table border="0" cellpadding="3" align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~A_{32}</math>
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,967: Line 2,020:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>
<math>~20\gamma_c - 8 = \frac{8(7 + \alpha_c)}{3-\alpha_c} \approx 12.17774\, ;</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
and, from the perspective of the envelope,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~B_{32}</math>
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,979: Line 2,037:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
<math>~\gamma_e\biggl[2\alpha_e + (c_0 + 9)(c_0+14)\biggr] \frac{\rho_e}{\rho_c} \approx 12.17774\, ,</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.118016</math>.
<font color="darkblue"><b>''Bottom-Right Panel:''</b></font>&nbsp; The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale.  Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~C_{32}</math>
<math>~y</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,991: Line 2,056:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6)  } \, .</math>
<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the outer edge.  In this type of log-amplitude plot, the eigenfunction's various ''nodes'' &#8212; that is, radial locations where <math>~x</math> passes through zero &#8212; are highlighted; here, specifically, there is one node inside the core and two nodes reside in the envelope.  Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope.  Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.41)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.12)</math>.


In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace <math>~Q_{31}</math> with,
NOTE:  As may be ascertained from our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes]], equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to ''any'' value, <math>~0 \le \rho_e/\rho_c \le 1</math>.  We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models.  Instead, our ability to derive analytically prescribed eigenvectors is [[User:Tohline/SSC/Stability/BiPolytrope0_0#KeyConstraint|limited by the constraint]],
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 2,003: Line 2,069:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{32}</math>
<math>~\frac{\rho_e}{\rho_c}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>
\frac{ -  252 (1+2\Chi)^2 + 396 (1+2\Chi)^4 }{35 -  126 (1+2\Chi)^2 + 99 (1+2\Chi)^4 } \, .
</math>
   </td>
   </td>
</tr>
</tr>
Line 2,017: Line 2,081:
</div>
</div>


The black-dotted curve in the ''Bottom-Right Panel'' displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model &#8212; such as the example illustrated here &#8212; that has an analytically prescribable radial-oscillation eigenvector.


===Check Surface Boundary Condition===


And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,
Let's determine analytically the logarithmic derivative of the envelope segment of the eigenfunction, at the surface <math>~(\xi = 1/q)</math> of the configuration.  Given that,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 2,025: Line 2,091:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~ \gamma_c </math>
<math>~x_{\ell=3} |_\mathrm{env}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,031: Line 2,097:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 +
<math>~
\gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\}
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{3} \xi^{3} +  q^6 A_{3}B_{3}\xi^{6} + q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3}  + q^6 A_{3}B_{3} + q^9 A_{3}B_{3}C_{3} }\biggr] \, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


where,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~A_{3}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,045: Line 2,116:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{42} \biggl\{ 8 +
<math>~\frac{- 6(c_0 + 7)}{(2c_0 + 5 )} \, ,</math>
\gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


<div align="center">
<table border="1" cellpadding="8">
<tr><td align="center" colspan="3">Figure 3:  Results for mode <math>~(\ell,j) = (3,2)</math></td></tr>
<tr><td align="center" rowspan="9">
[[File:Quartic32MontageSm.png|850px|Results for (ell,j) = (3,2)]]
</td>
<td align="right"><math>~q</math></td>
<td align="center"><math>~0.5940404</math></td>
</tr>
<tr>
<tr>
<td align="right"><math>~c_0</math></td>
  <td align="right">
<td align="center"><math>~-0.5381593</math></td>
<math>~B_{3}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
  </td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~\gamma_e</math></td>
  <td align="right">
<td align="center"><math>~1.0563279</math></td>
<math>~C_{3}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{- 2( c_0 + 10)}{3(2 c_0 + 11)  } \, ,</math>
  </td>
</tr>
</tr>
<tr>
<td align="right"><math>~\alpha_e</math></td>
<td align="center"><math>~-0.7867031</math></td>
</tr>
<tr>
<td align="right"><math>~n_e</math></td>
<td align="center"><math>~17.75</math></td>
</tr>
<tr>
<td align="right"><math>~\nu</math></td>
<td align="center"><math>~0.4730849</math></td>
</tr>
<tr>
<td align="right"><math>~\rho_e/\rho_c</math></td>
<td align="center"><math>~0.2954048</math></td>
</tr>
<tr>
<td align="right"><math>~\gamma_c</math></td>
<td align="center"><math>~1.0251098</math></td>
</tr>
<tr>
<td align="right"><math>~\alpha_c</math></td>
<td align="center"><math>~-0.9020212</math></td>
</tr>
</table>
</table>
</div>
</div>


==Try Mode33 <math>~(\ell,j) = (3,3)</math>==
we have,
 
Here, we set the quantum numbers to <math>~(\ell,j) = (3,3)</math>.  In this case, drawing from our summary of [[User:Tohline/SSC/UniformDensity#Sterne.27s_Presentation|Sterne's presentation]], the core eigenfunction segment should be,
 
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 2,108: Line 2,152:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~x_{j=3} |_\mathrm{core}</math>
<math>~\frac{ d\ln x_{\ell=3} }{d\ln \xi}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,114: Line 2,158:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{1}{x} \biggl\{
\frac{3\cdot 5 \cdot 7 -  33\cdot 3\cdot 7 (1+2q^3)^2 \xi^2 + 429\cdot 3 (1+2q^3)^4 \xi^4 - 143\cdot 5 (1+2q^3)^6 \xi^6}{3\cdot 5 \cdot 7 -  33\cdot 3\cdot 7 (1+2q^3)^
c_0\xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{3} \xi^{3} + q^6 A_{3}B_{3}\xi^{6} + q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3}  +  q^6 A_{3}B_{3} +  q^9 A_{3}B_{3}C_{3} }\biggr]
+ 429\cdot 3 (1+2q^3)^4 - 143\cdot 5 (1+2q^3)^6 }  
+ \xi^{c_0}\biggl[ \frac{ 3 q^3 A_{3} \xi^{3} + 6q^6 A_{3}B_{3}\xi^{6} + 9 q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3} + q^6 A_{3}B_{3} +  q^9 A_{3}B_{3}C_{3} }\biggr]
\, ,</math>
\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
</table>
 
</div>
and, [[#Try_Mode31|as when addressing Mode31, above]], throughout the envelope, the eigenfunction should be,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~x_{\ell=3} |_\mathrm{env}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,135: Line 2,175:
   <td align="left">
   <td align="left">
<math>~
<math>~
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{33} \xi^{3} +  q^6 A_{33}B_{33}\xi^{6} + q^9 A_{33}B_{33}C_{33}\xi^{9} }{ 1 +  q^3 A_{33} +  q^6 A_{33}B_{33} +  q^9 A_{33}B_{33}C_{33} }\biggr] \, ,
c_0
+ \biggl[ \frac{ 3 q^3 A_{3} \xi^{3} +  6q^6 A_{3}B_{3}\xi^{6} + 9 q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 +  q^3 A_{3} \xi^{3} +  q^6 A_{3}B_{3}\xi^{6} +  q^9 A_{3}B_{3}C_{3}\xi^{9} }\biggr] \, .
</math>
</math>
   </td>
   </td>
Line 2,141: Line 2,182:
</table>
</table>
</div>
</div>
where,
 
This means that,
<div align="center">
<div align="center">
<table border="0" cellpadding="3" align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~A_{33}</math>
<math>~\frac{ d\ln x_{\ell=3} }{d\ln \xi} \biggr|_{\xi=1/q}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,152: Line 2,195:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>
<math>~
c_0
+ \biggl[ \frac{ 3 A_{3}  +  6A_{3}B_{3} + 9 A_{3}B_{3}C_{3} }{ 1 +  A_{3}  + A_{3}B_{3}+ A_{3}B_{3}C_{3} }\biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 2,158: Line 2,204:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~B_{33}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,164: Line 2,210:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
<math>~
c_0
+ \biggl\{ \frac{ -3 [108(c_0 + 7)(c_0 +4)(2 c_0 + 11)] 
+  6[18(c_0 + 7) (6c_0 + 51)(2 c_0 + 11)] - 9 [ 12(c_0 + 7)(6c_0 + 51)( c_0 + 10) ] }{ 18(2c_0+5)(c_0+4)(2c_0+11) -  108(c_0 + 7)(c_0 +4)(2 c_0 + 11) 
+  18(c_0 + 7) (6c_0 + 51)(2 c_0 + 11) -  12(c_0 + 7)(6c_0 + 51)( c_0 + 10) }\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
Line 2,170: Line 2,221:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~C_{33}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,176: Line 2,227:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6)  } \, .</math>
<math>~
c_0
+ 108\biggl\{ \frac{ - 3(2c_0^3 + 33c_0^2 + 177c_0 + 308) 
+  (12c_0^3 + 252c_0^2 + 1737c_0 + 3927) - (6c_0^3 + 153c_0^2 + 1287c_0 + 3570) }{ 18(4c_0^3 + 48c_0^2 + 183c_0 +220) -  108(2c_0^3 + 33c_0^2 + 177c_0 + 308) 
+  18(12c_0^3 + 252c_0^2 + 1737c_0 + 3927) - 12(6c_0^3 + 153c_0^2 + 1287c_0 + 3570) }\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace <math>~Q_{31}</math> with,
<tr>
<div align="center">
  <td align="right">
<table border="0" cellpadding="5" align="center">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0
+ 108\biggl\{ \frac{0 c_0^3  + 0 c_0^2 -81 c_0 -567 }{ 0 c_0^3  + 0c_0^2 + 0c_0 -1458  }\biggr\}  
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0 + 6(c_0 + 7)
</math>
  </td>
</tr>


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{33}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{ -  2\cdot 33\cdot 3\cdot 7 (1+2q^3)^2 + 4\cdot 429\cdot 3 (1+2q^3)^4  - 6\cdot 143\cdot 5 (1+2q^3)^6 }{3\cdot 5 \cdot 7 -  33\cdot 3\cdot 7 (1+2q^3)^2 
7(c_0 + 6) \, ,
+ 429\cdot 3 (1+2q^3)^4  - 143\cdot 5 (1+2q^3)^6 }
\, .
</math>
</math>
   </td>
   </td>
Line 2,204: Line 2,281:
</div>
</div>


which precisely matches the desired surface boundary condition for <math>~\ell=3</math> that has been detailed in [[User:Tohline/SSC/Stability/BiPolytrope0_0#DesiredBoundaryCondition|an accompanying summary discussion]].


{{ LSU_WorkInProgress }}
==Try Mode32 <math>~(\ell,j) = (3,2)</math>==
===Setup32===
Here, we set the quantum numbers to <math>~(\ell,j) = (3,2)</math>.  In this case, [[#Try_Mode22|as when addessing Mode22, above]], the core eigenfunction segment should be,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x_{j=2} |_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } \, ;</math>
  </td>
</tr>
</table>
</div>
and, [[#Try_Mode31|as when addressing Mode31, above]], throughout the envelope, the eigenfunction should be,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x_{\ell=3} |_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{32} \xi^{3} +  q^6 A_{32}B_{32}\xi^{6} +  q^9 A_{32}B_{32}C_{32}\xi^{9} }{ 1 +  q^3 A_{32}  +  q^6 A_{32}B_{32} +  q^9 A_{32}B_{32}C_{32} }\biggr] \, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
  <td align="right">
<math>~A_{32}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~B_{32}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~C_{32}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{- 2( c_0 + 10)}{3( 2c_0 + 11)  } \, .</math>
  </td>
</tr>
</table>
</div>
In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace <math>~Q_{31}</math> with,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~Q_{32}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{ -  252(1+8q^3) (1+2\Chi)^2  + 396 (1+2\Chi)^4 }{35(1+8q^3)^2 -  126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } \, .
</math>
  </td>
</tr>
</table>
</div>


And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,
And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
 
 
<tr>
<tr>
   <td align="right">
  <td align="right">
<math>~ \gamma_c </math>
<math>~ \gamma_c </math>
   </td>
  </td>
   <td align="center">
  <td align="center">
<math>~=</math>
<math>~=</math>
   </td>
  </td>
   <td align="left">
  <td align="left">
<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 +
<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 +
\gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\}  
\gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\}
</math>
</math>
   </td>
  </td>
</tr>
</tr>
 
 
<tr>
<tr>
   <td align="right">
  <td align="right">
&nbsp;
&nbsp;
   </td>
  </td>
   <td align="center">
  <td align="center">
<math>~=</math>
<math>~=</math>
   </td>
  </td>
   <td align="left">
  <td align="left">
<math>~\frac{1}{72} \biggl\{ 8 +
<math>~\frac{1}{42} \biggl\{ 8 +
\gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, .
\gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, .
</math>
</math>
   </td>
  </td>
</tr>
</tr>
</table>
</table>
</div>
 
 
<!-- DELETE (moved below)
<div align="center">
[[File:Model32Montage1.png|750px|Montage of Stability Results for (ell,j) = (3,2) quantum numbers]]
</div>
-->
 
===Illustration32===
 
<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center" colspan="2">
Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes<br />
 
<font color="red">Quantum Numbers:</font>&nbsp; <math>~(\ell,j) = (3,2)</math>
  </th>
</tr>
<tr>
  <th align="center" width="50%"> Analyzable Model Sequence</th>
  <th align="center" width="50%">One Example Eigenfunction</th>
</tr>
<tr>
  <td align="center" colspan="2">
[[File:Model32MontageCorrected.png|600px|Montage of Stability Results for (ell,j) = (3,2) quantum numbers]]
  </td>
</tr>
</table>
</div>
 
 
<font color="darkblue"><b>''Top-Left Panel:''</b></font>&nbsp; Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> &#8212; where <math>~\gamma_e</math> is the adiabatic exponent of the envelope &#8212; in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (3,2)</math>.  Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the ''core'' [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>.  The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
 
<font color="darkblue"><b>''Bottom-Left Panel:''</b></font> &nbsp; Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>.  Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top.  The core is ''more'' compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line.  The yellow circular marker identifies the same example model as it does in the top-left panel.
 
<font color="darkblue"><b>''Top-Right Panel:''</b></font>&nbsp; Displays &#8212; as a function of the fractional radius, <math>~r_0/R = q\xi</math> &#8212;  the analytically determined eigenfunction for the <math>~(\ell,j) = (3,2)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,
<div align="center">
<math>~q = \biggl[0.01 +  185\biggl( \frac{0.98}{3\times 99} \biggr)\biggr] \approx 0.6204377</math>
</div>
and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (-2.3774068, + 0.8972496, +0.0231304)</math>.  Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction as defined by the expression for <math>~x_{j=2} |_\mathrm{core}</math>, given [[#Setup32|above]]; and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction as defined by the expression for <math>~x_{\ell=3} |_\mathrm{env}</math>, given [[#Setup32|above]].
 
The corresponding eigenfrequency is, from the perspective of the core,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~42\gamma_c - 8 = \frac{8(18 + \alpha_c)}{3-\alpha_c} \approx 48.4351\, ;</math>
  </td>
</tr>
</table>
</div>
and, from the perspective of the envelope,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\gamma_e\biggl[2\alpha_e + (c_0 + 9)(c_0+14)\biggr] \frac{\rho_e}{\rho_c} \approx 48.4351 \, ,</math>
  </td>
</tr>
</table>
</div>
where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.323257</math>.
 
<font color="darkblue"><b>''Bottom-Right Panel:''</b></font>&nbsp; The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale.  Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>
  </td>
</tr>
</table>
</div>
where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the enter of the configuration.  In this type of log-amplitude plot, the eigenfunction's various ''nodes'' &#8212; that is, radial locations where <math>~x</math> passes through zero &#8212; are highlighted; here, specifically, there is one node inside the core and one node resides in the envelope.  Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope.  Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.62)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.32)</math>.
 
NOTE:  As may be ascertained from our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes]], equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to ''any'' value, <math>~0 \le \rho_e/\rho_c \le 1</math>.  We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models.  Instead, our ability to derive analytically prescribed eigenvectors is [[User:Tohline/SSC/Stability/BiPolytrope0_0#KeyConstraint|limited by the constraint]],
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\rho_e}{\rho_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>
  </td>
</tr>
</table>
</div>
 
The black-dotted curve in the ''Bottom-Right Panel'' displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model &#8212; such as the example illustrated here &#8212; that has an analytically prescribable radial-oscillation eigenvector.
 
==Try Mode33 <math>~(\ell,j) = (3,3)</math>==
 
===Setup33===
 
Here, we set the quantum numbers to <math>~(\ell,j) = (3,3)</math>.  In this case, drawing from our summary of [[User:Tohline/SSC/UniformDensity#Sterne.27s_Presentation|Sterne's presentation]], the core eigenfunction segment should be,
 
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~x_{j=3} |_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{3\cdot 5 \cdot 7 -  33\cdot 3\cdot 7 (1+2q^3)^2 \xi^2 + 429\cdot 3 (1+2q^3)^4 \xi^4 - 143\cdot 5 (1+2q^3)^6 \xi^6}{3\cdot 5 \cdot 7 -  33\cdot 3\cdot 7 (1+2q^3)^2 
+ 429\cdot 3 (1+2q^3)^4  - 143\cdot 5 (1+2q^3)^6 }
\, ,</math>
  </td>
</tr>
</table>
</div>
and, [[#Try_Mode31|as when addressing Mode31, above]], throughout the envelope, the eigenfunction should be,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~x_{\ell=3} |_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{33} \xi^{3} +  q^6 A_{33}B_{33}\xi^{6} +  q^9 A_{33}B_{33}C_{33}\xi^{9} }{ 1 +  q^3 A_{33}  +  q^6 A_{33}B_{33} +  q^9 A_{33}B_{33}C_{33} }\biggr] \, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
  <td align="right">
<math>~A_{33}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~B_{33}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~C_{33}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6)  } \, .</math>
  </td>
</tr>
</table>
</div>
 
In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace <math>~Q_{31}</math> with,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~Q_{33}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{ -  2\cdot 33\cdot 3\cdot 7 (1+2q^3)^2 + 4\cdot 429\cdot 3 (1+2q^3)^4  - 6\cdot 143\cdot 5 (1+2q^3)^6 }{3\cdot 5 \cdot 7 -  33\cdot 3\cdot 7 (1+2q^3)^2 
+ 429\cdot 3 (1+2q^3)^4  - 143\cdot 5 (1+2q^3)^6 }
\, .
</math>
  </td>
</tr>
</table>
</div>
 
 
 
And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~ \gamma_c </math>
   </td>
   <td align="center">
<math>~=</math>
   </td>
   <td align="left">
<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 +
\gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\}  
</math>
   </td>
</tr>
 
<tr>
   <td align="right">
&nbsp;
   </td>
   <td align="center">
<math>~=</math>
   </td>
   <td align="left">
<math>~\frac{1}{72} \biggl\{ 8 +
\gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, .
</math>
   </td>
</tr>
</table>
</div>
 
 
<div align="center">
<table border="1" cellpadding="8">
<tr><td align="center" colspan="3">Figure 4:  Results for mode <math>~(\ell,j) = (3,3)</math></td></tr>
<tr><td align="center" rowspan="9">
[[|850px|Results for (ell,j) = (3,3)]]
</td>
<td align="right"><math>~q</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~c_0</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~\gamma_e</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~\alpha_e</math></td>
<td align="center"><math>~0. </math></td>
</tr>
<tr>
<td align="right"><math>~n_e</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~\nu</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~\rho_e/\rho_c</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~\gamma_c</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<td align="right"><math>~\alpha_c</math></td>
<td align="center"><math>~ </math></td>
</tr>
 
</table>
</div>
 
===Illustration33===
 
<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center" colspan="2">
Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes<br />
 
<font color="red">Quantum Numbers:</font>&nbsp; <math>~(\ell,j) = (3,3)</math>
  </th>
</tr>
<tr>
  <th align="center" width="50%"> Analyzable Model Sequence</th>
  <th align="center" width="50%">One Example Eigenfunction</th>
</tr>
<tr>
  <td align="center" colspan="2">
[[File:Model33Montage.png|800px|Montage of Stability Results for (ell,j) = (3,3) quantum numbers]]
  </td>
</tr>
</table>
</div>
 
 
<font color="darkblue"><b>''Top-Left Panel:''</b></font>&nbsp; Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> &#8212; where <math>~\gamma_e</math> is the adiabatic exponent of the envelope &#8212; in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (3,3)</math>.  Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the ''core'' [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>.  The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.
 
<font color="darkblue"><b>''Bottom-Left Panel:''</b></font> &nbsp; Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>.  Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top.  The core is ''more'' compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line.  The yellow circular marker identifies the same example model as it does in the top-left panel.
 
<font color="darkblue"><b>''Top-Right Panel:''</b></font>&nbsp; Displays &#8212; as a function of the fractional radius, <math>~r_0/R = q\xi</math> &#8212;  the analytically determined eigenfunction for the <math>~(\ell,j) = (3,2)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,
<div align="center">
<math>~q = 0.01 + 642\biggl( \frac{0.98}{990} \biggr) \approx 0.6455152 \, ,</math>
</div>
and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (+0.3511525, +0.8256130, -0.0183637)</math>.  Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction as defined by the expression for <math>~x_{j=3} |_\mathrm{core}</math>, given [[#Setup33|above]]; and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction as defined by the expression for <math>~x_{\ell=3} |_\mathrm{env}</math>, given [[#Setup33|above]].
 
The corresponding eigenfrequency is, from the perspective of the core,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~72\gamma_c - 8 = \frac{8(33 + \alpha_c)}{3-\alpha_c} \approx 87.41594 \, ;</math>
  </td>
</tr>
</table>
</div>
and, from the perspective of the envelope,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\gamma_e\biggl[2\alpha_e + (c_0 + 9)(c_0+14)\biggr] \frac{\rho_e}{\rho_c} \approx 87.41594 \, ,</math>
  </td>
</tr>
</table>
</div>
where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.3497877</math>.
 
<font color="darkblue"><b>''Bottom-Right Panel:''</b></font>&nbsp; The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale.  Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\tfrac{1}{10} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>
  </td>
</tr>
</table>
</div>
where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the outer edge.  In this type of log-amplitude plot, the eigenfunction's various ''nodes'' &#8212; that is, radial locations where <math>~x</math> passes through zero &#8212; are highlighted; here, specifically, there are three nodes inside the core and one node resides in the envelope.  Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope.  Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.65)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.35)</math>.
 
NOTE:  As may be ascertained from our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes]], equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to ''any'' value, <math>~0 \le \rho_e/\rho_c \le 1</math>.  We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models.  Instead, our ability to derive analytically prescribed eigenvectors is [[User:Tohline/SSC/Stability/BiPolytrope0_0#KeyConstraint|limited by the constraint]],
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\rho_e}{\rho_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>
  </td>
</tr>
</table>
</div>
 
The black-dotted curve in the ''Bottom-Right Panel'' displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model &#8212; such as the example illustrated here &#8212; that has an analytically prescribable radial-oscillation eigenvector.
 
==Summary Plots==
 
<!-- DELETE in favor of more complete table, immediately following
<div align="center">
<table border="1" align="center"><tr><td align="center">
[[File:DenRatioAmpsMontage01.png|900px|Analytic Eigenfunctions for zero-zero bipolytropes]]
</td></tr></table>
</div>
</div>
-->




<div align="center">
<div align="center">
<table border="1" cellpadding="8">
<table border="1" cellpadding="5">
<tr><td align="center" colspan="3">Figure 4:  Results for mode <math>~(\ell,j) = (3,3)</math></td></tr>
<tr><td align="center" rowspan="9">
[[|850px|Results for (ell,j) = (3,3)]]
</td>
<td align="right"><math>~q</math></td>
<td align="center"><math>~ </math></td>
</tr>
<tr>
<tr>
<td align="right"><math>~c_0</math></td>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (2,1)</math></td>
<td align="center"><math>~ </math></td>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (2,2)</math></td>
  <td align="center" bgcolor="white">&nbsp;</td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~\gamma_e</math></td>
  <td align="center">[[File:Annotate21.png|300px|Log(amplitude) plot for (ell,j) = (2,1)]]</td>
<td align="center"><math>~ </math></td>
  <td align="center">[[File:Annotate22.png|300px|Log(amplitude) plot for (ell,j) = (2,2)]]</td>
  <td align="center">&nbsp;</td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~\alpha_e</math></td>
  <td align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 37.08874</math></td>
<td align="center"><math>~0. </math></td>
  <td align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.95210</math></td>
  <td align="center">&nbsp;</td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~n_e</math></td>
  <td align="center" bgcolor="lightblue" colspan="3">&nbsp;</td>
<td align="center"><math>~ </math></td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~\nu</math></td>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (3,1)</math></td>
<td align="center"><math>~ </math></td>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (3,2)</math></td>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (3,3)</math></td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~\rho_e/\rho_c</math></td>
  <td align="center">[[File:Annotate31.png|300px|Log(amplitude) plot for (ell,j) = (3,1)]]</td>
<td align="center"><math>~ </math></td>
  <td align="center">[[File:Annotate32.png|300px|Log(amplitude) plot for (ell,j) = (3,2)]]</td>
</tr>
  <td align="center">[[File:Annotate33.png|300px|Log(amplitude) plot for (ell,j) = (3,3)]]</td>
<tr>
<td align="right"><math>~\gamma_c</math></td>
<td align="center"><math>~ </math></td>
</tr>
</tr>
<tr>
<tr>
<td align="right"><math>~\alpha_c</math></td>
  <td align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 12.452545</math></td>
<td align="center"><math>~ </math></td>
  <td align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.05461</math></td>
  <td align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 87.41594</math></td>
</tr>
</tr>
</table>
</table>
</div>
==Summary Plots==
<div align="center">
<table border="1" align="center"><tr><td align="center">
[[File:DenRatioAmpsMontage01.png|900px|Analytic Eigenfunctions for zero-zero bipolytropes]]
</td></tr></table>
</div>
</div>



Latest revision as of 22:16, 10 February 2017

Searching for Additional Eigenvectors of Zero-Zero Bipolytropes

This chapter is an extension of two accompanying discussions: The original discovery and detailed derivation; and the more readable, summary.

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

In our accompanying summary, we have demonstrated how analytically specified eigenvectors can be constructed for the mode labeled, <math>~(\ell, j) = (2,1)</math>. This was done by specifying <math>~\gamma_e</math>, then solving a quartic equation for <math>~q</math>. Shortly after completing this summary chapter, we noticed that an alternate approach may be to specify <math>~q</math>, then solve for <math>~\gamma_e</math>; and this path may be simpler because it may only involve solution of a quadratic equation. (Actually, we later have realized that the relevant equation is cubic, rather than quadratic. This is nevertheless simpler than the quartic equation.) If this proves to be the case, then it may also be possible to analytically construct eigenvectors of additional modes. Let's see.


CAUTION!   All of the "core" eigenfunction expressions are suspect because (I think!) they have been developed herein under the assumption that <math>~g^{-2} = (1+2q^3)^2</math>, whereas the relation should read,

<math>~\frac{1}{g^2}</math>

<math>~=</math>

<math>~\frac{(1+2q^3)^2}{1+8q^3} = \frac{(1+2\Chi)^2}{1+8\Chi }\, .</math>


Seek Alternate Solution to Mode21 <math>~(\ell,j) = (2,1)</math>

Setup21

According to STEP 4 in our accompanying summary discussion, we need to solve the following "derivative matching" expression:

<math>~ \frac{\gamma_c}{\gamma_e} \biggl[ \frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)} \biggr] </math>

<math>~=</math>

<math>~ c_0 + \frac{3A_{2}\Chi + 6A_{2}B_{2} \Chi^2}{1 + A_{2}\Chi + A_{2}B_{2}\Chi^2} - 3\biggl(\frac{\gamma_c}{\gamma_e} -1\biggr) </math>

<math>~\Rightarrow~~~ \frac{\gamma_c}{\gamma_e} \biggl[\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)}\biggr] + 3\biggl(\frac{\gamma_c}{\gamma_e} -1\biggr) </math>

<math>~=</math>

<math>~ \frac{c_0 + (c_0 + 3)A_{2}\Chi + (c_0 + 6)A_{2}B_{2} \Chi^2}{1 + A_{2}\Chi + A_{2}B_{2}\Chi^2} </math>

where, recognizing that, <math>~\alpha_e = c_0(c_0+2) \, ,</math>

<math>~A_{2}</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{c_0(c_0+5) - (c_0 + 6)(c_0 + 11)}{(c_0 + 3)(c_0+5) - \alpha_e}\biggr] </math>

 

<math>~=</math>

<math>~\biggl[ \frac{c_0^2 + 5c_0 - (c_0^2 + 17c_0 + 66)}{(c_0^2 + 8c_0 + 15) - (c_0^2+2c_0)}\biggr] </math>

 

<math>~=</math>

<math>~-\biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr) \, ,</math>

<math>~B_{2}</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{(c_0+3)(c_0+8) - (c_0 + 6)(c_0 + 11)}{(c_0 + 6)(c_0+8) - \alpha_e}\biggr] </math>

 

<math>~=</math>

<math>~\biggl[ \frac{(c_0^2 +11c_0 + 24) - (c_0^2 + 17c_0 + 66)}{(c_0^2+14c_0+48) - (c_0^2 + 2c_0)}\biggr] </math>

 

<math>~=</math>

<math>~-\biggl( \frac{c_0 + 7 }{2c_0+8}\biggr) \, . </math>

Here, we assume that <math>~\Chi \equiv q^3</math> is specified and seek the corresponding value of <math>~c_0</math>. Given that the LHS of this matching relation is known once <math>~\Chi</math> has been specified, in order to simplify notation we will also define,

<math>~Q</math>

<math>~\equiv</math>

<math>~\frac{\gamma_c}{\gamma_e} \biggl[\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)}\biggr] + 3\biggl(\frac{\gamma_c}{\gamma_e} -1\biggr) \, .</math>

Then the matching relation becomes,

<math>~Q</math>

<math>~=</math>

<math>~ \frac{c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2}{1 + A_{21}\Chi + A_{21}B_{21}\Chi^2} </math>

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

<math>~=</math>

<math>~[c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2 ] - Q[1 + A_{21}\Chi + A_{21}B_{21}\Chi^2 ] </math>

 

<math>~=</math>

<math>~\biggl[c_0 - (c_0 + 3)\biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr)\Chi + (c_0 + 6)\biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr)\biggl( \frac{c_0 + 7 }{2c_0+8}\biggr) \Chi^2 \biggr] - Q\biggl[1 - \biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr)\Chi + \biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr)\biggl( \frac{c_0 + 7 }{2c_0+8}\biggr)\Chi^2 \biggr] </math>

 

<math>~=</math>

<math>~\biggl[c_0(2c_0+5)(2c_0+8) - (c_0 + 3)(2c_0+8)( 4c_0 + 22)\Chi + (c_0 + 6)( 4c_0 + 22)(c_0 + 7 ) \Chi^2 \biggr] </math>

 

 

<math>~ - Q\biggl[(2c_0+5)(2c_0+8) - (2c_0+8)( 4c_0 + 22)\Chi + ( 4c_0 + 22)( c_0 + 7 )\Chi^2 \biggr] </math>

 

<math>~=</math>

<math>~\biggl\{c_0(4c_0^2 + 26c_0 + 40) + ( 4c_0 + 22)\Chi [(c_0^2 + 13c_0 + 42 ) \Chi- (2c_0^2 +14c_0 +24)] \biggr\} </math>

 

 

<math>~ - Q\biggl\{ (4c_0^2 + 26c_0 + 40) + ( 4c_0 + 22)\Chi [c_0(\Chi-2) + (7\Chi-8) ] \biggr\} </math>

 

<math>~=</math>

<math>~4c_0^3 +[ 26 - 4Q]c_0^2 + [40 - 26Q]c_0 - 40Q </math>

 

 

<math>~ + ( 4c_0 + 22)\Chi \biggl\{ [\Chi - 2]c_0^2 + [13\Chi -14 -Q (\Chi-2) ]c_0 + [42\Chi -Q (7\Chi-8) -24] \biggr\} </math>

 

<math>~=</math>

<math>~4c_0^3 +[ 26 - 4Q]c_0^2 + [40 - 26Q]c_0 - 40Q </math>

 

 

<math>~ + 4\Chi[\Chi - 2]c_0^3 + 4\Chi[13\Chi -14 -Q (\Chi-2) ]c_0^2 + 4\Chi[42\Chi -Q (7\Chi-8) -24]c_0 </math>

 

 

<math>~ + 22\Chi[\Chi - 2]c_0^2 + 22\Chi[13\Chi -14 -Q (\Chi-2) ]c_0 + 22\Chi[42\Chi -Q (7\Chi-8) -24] \, . </math>

Solve Cubic Equation

Here we draw from a separate discussion of solutions to a cubic equation.

Using <math>~y</math> in place of <math>~c_0</math>, this "derivative matching" relation can be written in the form of a standard cubic equation. Specifically,

<math>~a y^3 + b y^2 + c y + d</math>

<math>~=</math>

<math>~ 0 \, , </math>

where,

<math>~a</math>

<math>~\equiv</math>

<math>~ 4 + 4\Chi(\Chi - 2)\, , </math>

<math>~b</math>

<math>~\equiv</math>

<math>~ ( 26 - 4Q_{21}) + 4\Chi[ 13\Chi -14 - Q_{21} (\Chi-2) ] + 22\Chi (\Chi - 2)\, , </math>

<math>~c</math>

<math>~\equiv</math>

<math>~

[40  - 26 Q_{21}]+ 4\Chi[42\Chi  - Q_{21} (7\Chi-8) -24]  + 22\Chi[13\Chi  -14  - Q_{21}  (\Chi-2) ] \, ,

</math>

<math>~d</math>

<math>~\equiv</math>

<math>~ - 40 Q_{21}+ 22\Chi[42\Chi -Q_{21} (7\Chi-8) -24] \, . </math>

As is well known and documented — see, for example Wolfram MathWorld or Wikipedia's discussion of the topic — the roots of any cubic equation can be determined analytically. In order to evaluate the root(s) of our particular cubic equation, we have drawn from the utilitarian online summary provided by Eric Schechter at Vanderbilt University. For a cubic equation of the general form,

<math>~ay^3 + by^2 + cy + d = 0 \, ,</math>

a real root is given by the expression,

<math>~ y = p + \{z + [z^2 + (r-p^2)^3]^{1/2}\}^{1/3} + \{z - [z^2 + (r-p^2)^3]^{1/2}\}^{1/3} \, ,</math>

where,

<math>~p \equiv -\frac{b}{3a} \, ,</math>      <math>~z \equiv \biggl[p^3 + \frac{bc-3ad}{6a^2} \biggr] \, ,</math>      and      <math>~r=\frac{c}{3a} \, .</math>

(There is also a pair of imaginary roots, but they are irrelevant in the context of our overarching astrophysical discussion.)


Upon evaluation, we have found that the expression inside of the square root is negative over the region of parameter space that is of most physical interest. Hence, we need to call upon a separate discussion in which the cube root of complex numbers was discussed.


We'll shift to Wolfram's notation; specifically,

<math>~a_0 </math>

<math>~=</math>

<math>~\frac{d}{a}\, ,</math>

<math>~a_1 </math>

<math>~=</math>

<math>~\frac{c}{a}\, ,</math>

<math>~a_2 </math>

<math>~=</math>

<math>~\frac{b}{a}\, ,</math>

<math>~R </math>

<math>~\equiv</math>

<math>~\frac{3^2 a_2 a_1 - 3^3a_0 - 2a_2^3}{2\cdot 3^3}</math>

 

<math>~=</math>

<math>~\frac{bc - 3ad }{6 a^2} - \biggl( \frac{ b}{3 a} \biggr)^3 = z \, ,</math>

<math>~Q </math>

<math>~\equiv</math>

<math>~\frac{3 a_1 - a_2^2}{3^2}</math>

 

<math>~=</math>

<math>~\frac{ c }{3a} - \biggl(\frac{b}{3a}\biggr)^2 = r-p^2 \, .</math>

Then, after defining,

<math>~D \equiv Q^3 + R^2 = z^2 + (r-p^2)^3 \, ;</math>         <math>~S^3 \equiv R+ \sqrt{D} \, ;</math>         and         <math>~T^3 \equiv R- \sqrt{D} \, ;</math>

Wolfram states that the three roots of the cubic equation are (the first one being identical to the "real" root identified above),

<math>~y_1</math>

<math>~=</math>

<math>~p + (S + T) \, ,</math>

<math>~y_2</math>

<math>~=</math>

<math>~p - \frac{1}{2}(S + T) + \frac{1}{2} i \sqrt{3}(S-T) \, ,</math>

<math>~y_2</math>

<math>~=</math>

<math>~p - \frac{1}{2}(S + T) - \frac{1}{2} i \sqrt{3}(S-T) \, .</math>

Now, whenever <math>~D</math> is intrinsically negative, we need to treat both <math>~S^3</math> and <math>~T^3</math> as complex numbers. If we define,

<math>~\mathfrak{r} \equiv (R^2 + |D|)^{1 / 2} \, ,</math>

      and      

<math>~\theta \equiv \tan^{-1}\biggl[ \frac{\sqrt{|D|}}{R} \biggr] \, ,</math>

then we can write,

<math>~S^3 = \mathfrak{r} e^{+i\theta} \, ,</math>

      and      

<math>~T^3 = \mathfrak{r} e^{-i\theta} \, .</math>

As is explained in this online resource, both <math>~S</math> and <math>~T</math> must formally submit to three separate roots tagged by the integer index, <math>~(j=0,1,2)</math>. Working only with the <math>~j=0</math> root for both, we find that the above expressions for the three roots of our cubic equation become,

<math>~y_1</math>

<math>~=</math>

<math>~p + 2 \mathfrak{r}^{1 / 3} \cos\biggl(\frac{\theta}{3}\biggr) \, ,</math>

<math>~y_2</math>

<math>~=</math>

<math>~p - \mathfrak{r}^{1 / 3} \biggl[ \cos\biggl(\frac{\theta}{3}\biggr) + \sqrt{3} \sin\biggl(\frac{\theta}{3}\biggr) \biggr] \, ,</math>

<math>~y_3</math>

<math>~=</math>

<math>~p - \mathfrak{r}^{1 / 3} \biggl[ \cos\biggl(\frac{\theta}{3}\biggr) - \sqrt{3} \sin\biggl(\frac{\theta}{3}\biggr) \biggr] \, .</math>

We have deduced empirically that <math>~y_3</math> is the root that is physically relevant in our case. That is to say, for a given <math>~0 < q < 1</math>,

<math>~c_0 = p - \mathfrak{r}^{1 / 3} \biggl[ \cos\biggl(\frac{\theta}{3}\biggr) - \sqrt{3} \sin\biggl(\frac{\theta}{3}\biggr) \biggr] \, .</math>

In turn, for a given value of <math>~q</math>, the corresponding value of <math>~\alpha_e</math> is obtained via the relation, <math>~\alpha_e = c_0(c_0 + 2)</math>.

The right-hand panel of Figure 1 presents a plot of this <math>~c_0(q)</math> function; actually what has been plotted is the inverted relation, <math>~q(\alpha_e)</math>. The open, red circular markers trace the portion of the function that provides physically viable solutions, in the sense that the corresponding value of <math>~\alpha_c</math> lies between the values, negative one and three; the filled, light=blue circular markers identify roots of the cubic equation that are not physically viable.

In the left-hand panel of Figure 1, we re-display a plot that has been discussed in an accompanying chapter. It contains a plot (blue markers) of the same same <math>~q(\alpha_e)</math> function but, this time, as determined from the root of a quartic equation. In order to illustrate more clearly that the two curves are the same, we have plotted the quartic solution (small, purple circular markers) on top of the cubic solution in the right-hand panel.

Figure 1:  Comparing Roots to Quartic and Cubic Equations

quartic solution cubic solution


Illustration21

Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes

Quantum Numbers:  <math>~(\ell,j) = (2,1)</math>

Analyzable Model Sequence One Example Eigenfunction

Montage of Stability Results for (ell,j) = (2,1) quantum numbers


Top-Left Panel:  Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> — where <math>~\gamma_e</math> is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (2,1)</math>. Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>. The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.

Bottom-Left Panel:   Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>. Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line. The yellow circular marker identifies the same example model as it does in the top-left panel.

Top-Right Panel:  Displays — as a function of the fractional radius, <math>~r_0/R = q\xi</math> — the analytically determined eigenfunction for the <math>~(\ell,j) = (2,1)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,

<math>~\alpha_e =-0.35</math>       <math>~\Rightarrow</math>       <math>~c_0~\mathrm{(plus)} = -0.1937742</math>

and, correspondingly (via our separate solution of the governing quartic equation), <math>~(q,\alpha_c) = (0.6840119, +0.8326585)</math>. (This same plot appears in Figure 1 of our accompanying summary.) Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction, namely,

<math>~x_{j=1} |_\mathrm{core}</math>

<math>~=</math>

<math>~ \frac{5 (1+8q^3) - 7 (1+2q^3)^2 \xi^2}{5 (1+8q^3)-7(1+2q^3)^2} \, ;</math>

and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction, namely,

<math>~x_{\ell=2} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{21} \xi^{3} + q^6 A_{21}B_{21}\xi^{6} }{ 1 + q^3 A_{21} + q^6 A_{21}B_{21} }\biggr] \, , </math>

where the coefficients, <math>~A_{21}, B_{21}</math>, are as defined above in terms of the parameter, <math>~c_0</math>. The corresponding eigenfrequency is, from the perspective of the core,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~20\gamma_c - 8 = \frac{8(7 + \alpha_c)}{3-\alpha_c} \approx 28.9115809\, ;</math>

and, from the perspective of the envelope,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~3\gamma_e\biggl[\alpha_e + 5c_0 + 22\biggr] \frac{\rho_e}{\rho_c} \approx 28.9115807\, ,</math>

where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.3902664</math>.

Bottom-Right Panel:  The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,

<math>~y</math>

<math>~=</math>

<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>

where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the outer edge. In this type of log-amplitude plot, the eigenfunction's various nodes — that is, radial locations where <math>~x</math> passes through zero — are highlighted; here, specifically, there is one node inside the core and one node resides in the envelope. Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope. Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.684)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.39)</math>.

NOTE: As may be ascertained from our general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes, equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to any value, <math>~0 \le \rho_e/\rho_c \le 1</math>. We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,

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

<math>~=</math>

<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>

The black-dotted curve in the Bottom-Right Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radial-oscillation eigenvector.

Try Mode22 <math>~(\ell,j) = (2,2)</math>

Setup22

In this case we need to replace the core eigenfunction segment that was specified in STEP 3 in our separate discussion of mode <math>~(\ell,j) = (2,1)</math> with the following:

<math>~x_{j=2} |_\mathrm{core}</math>

<math>~=</math>

<math>~ \frac{35(1+8q^3)^2 - 126 (1+8q^3) (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } </math>


So, following the procedure outlined in STEP 4 of our separate discussion, the requirement that the first derivatives match at the interface becomes,

<math>~\frac{ - 252 (1+8q^3) (1+2q^3)^2 + 396 (1+2q^3)^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 }</math>

<math>~=</math>

<math>~ \frac{c_0 + (c_0 + 3)A_{22}q^3 + (c_0 + 6)A_{22}B_{22} q^6}{1 + A_{22}q^3 + A_{22}B_{22}q^6} \, , </math>

where,

<math>~A_{22}</math>

<math>~=</math>

<math>~A_{21} = -\biggl( \frac{ 4c_0 + 22}{2c_0 + 5}\biggr) \, ,</math>

<math>~B_{22}</math>

<math>~=</math>

<math>~B_{21} = -\biggl( \frac{c_0 + 7 }{2c_0+8}\biggr) \, . </math>

So, as described above, in order to determine <math>~c_0</math> for any specified value of <math>~q</math> — or, equivalently, <math>~\Chi \equiv q^3</math> — we simply need to replace <math>~Q_{21}</math> with <math>~Q_{22}</math> in the definition of the coefficient of each term in the cubic expression, where,

<math>~Q_{22}</math>

<math>~\equiv</math>

<math>~ \frac{ - 252 (1+8q^3) (1+2q^3)^2 + 396 (1+2q^3)^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 }\, . </math>


Finally, as explained in our summary discussion, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,

<math>~ \gamma_c </math>

<math>~=</math>

<math>~\frac{1}{[ 6 + 2j(2j+5)] } \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, . </math>


Illustration22

Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes

Quantum Numbers:  <math>~(\ell,j) = (2,2)</math>

Analyzable Model Sequence One Example Eigenfunction

Montage of Stability Results for (ell,j) = (2,1) quantum numbers


Top-Left Panel:  Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> — where <math>~\gamma_e</math> is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (2,2)</math>. Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>. The red marker that is farthest to the left identifies the model whose analytically determined eigenfunction is displayed, as an example, in the right-hand panels.

Bottom-Left Panel:   Plotted points (only 2, here!) show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>. Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line.

Top-Right Panel:  Displays — as a function of the fractional radius, <math>~r_0/R = q\xi</math> — the analytically determined eigenfunction for the <math>~(\ell,j) = (2,2)</math> mode in the model identified by the red circular marker that is farthest to the left in both left-hand panels, for which,

<math>~q = \biggl[0.005 + 179\biggl( \frac{0.985}{200} \biggr) \biggr] = 0.886575</math>

and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (-2.332785, +0.7763158, -0.9146699)</math>. Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction, namely,

<math>~x_{j=2} |_\mathrm{core}</math>

<math>~=</math>

<math>~ \frac{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } \, ;</math>

and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction, namely,

<math>~x_{\ell=2} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{22} \xi^{3} + q^6 A_{22}B_{22}\xi^{6} }{ 1 + q^3 A_{22} + q^6 A_{22}B_{22} }\biggr] \, , </math>

where the coefficients, <math>~A_{22}, B_{22}</math>, are as defined above in terms of the parameter, <math>~c_0</math>. The corresponding eigenfrequency is, from the perspective of the core,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~42\gamma_c - 8 = \frac{8(18 + \alpha_c)}{3-\alpha_c} \approx 34.915496 \, ;</math>

and, from the perspective of the envelope,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~3\gamma_e\biggl[\alpha_e + 5c_0 + 22\biggr] \frac{\rho_e}{\rho_c} \approx 34.915496 \, ,</math>

where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.5822407</math>.

Bottom-Right Panel:  The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,

<math>~y</math>

<math>~=</math>

<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>

where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the center of the configuration. In this type of log-amplitude plot, the eigenfunction's various nodes — that is, radial locations where <math>~x</math> passes through zero — are highlighted; here, specifically, there are two nodes inside the core and none in the envelope, although one of the nodes in the core lies just inside the core/envelope interface. Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope. Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.89)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.58)</math>.

NOTE: As may be ascertained from our general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes, equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to any value, <math>~0 \le \rho_e/\rho_c \le 1</math>. We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,

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

<math>~=</math>

<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>

The black-dotted curve in the Bottom-Right Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radial-oscillation eigenvector.

Try Mode31 <math>~(\ell,j) = (3,1)</math>

Setup31

In this case we need to replace the envelope eigenfunction segment that was specified in STEP 2 in our separate discussion of mode <math>~(\ell,j) = (2,1)</math> with the following: For the case of <math>~\ell=2</math>, this means that, throughout the envelope, the eigenfunction is,

<math>~x_{\ell=3} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{31} \xi^{3} + q^6 A_{31}B_{31}\xi^{6} + q^9 A_{31}B_{31}C_{31}\xi^{9} }{ 1 + q^3 A_{31} + q^6 A_{31}B_{31} + q^9 A_{31}B_{31}C_{31} }\biggr] \, , </math>

where, the values of the newly introduced coefficients,

<math>~A_{31}</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{c_0(c_0+5) - (c_0 + 9)(c_0 + 14)}{(c_0 + 3)(c_0+5) - \alpha_e}\biggr] </math>

 

<math>~=</math>

<math>~\biggl[ \frac{c_0^2 + 5c_0 - (c_0^2 + 23c_0 + 126)}{(c_0^2 + 8c_0 + 15) - (c_0^2 + 2c_0)}\biggr] </math>

 

<math>~=</math>

<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>

<math>~B_{31}</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{(c_0+3)(c_0+8) - (c_0 + 9)(c_0 + 14)}{(c_0 + 6)(c_0+8) - \alpha_e}\biggr] </math>

 

<math>~=</math>

<math>~\biggl[ \frac{(c_0^2 + 11c_0 + 24) - (c_0^2 + 23c_0 + 126)}{(c_0^2 + 14c_0 +48) - (c_0^2 + 2c_0) }\biggr] </math>

 

<math>~=</math>

<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>

<math>~C_{31}</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{(c_0+6)(c_0+11) - (c_0 + 9)(c_0 + 14)}{(c_0 + 9)(c_0+11) - \alpha_e}\biggr] </math>

 

<math>~=</math>

<math>~\biggl[ \frac{(c_0^2 +17c_0 + 66) - (c_0^2 + 23 c_0 + 126)}{(c_0^2 + 20 c_0 + 99) - (c_0^2 + 2c_0) }\biggr] </math>

 

<math>~=</math>

<math>~\frac{- 2( c_0 + 10)}{3(2 c_0 + 11) } \, .</math>

Then, after defining,

<math>~Q_{31}</math>

<math>~\equiv</math>

<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)} \, ,</math>

the matching relation becomes,

<math>~Q_{31}</math>

<math>~=</math>

<math>~ \frac{c_0 + (c_0 + 3)A_{31}\Chi + (c_0 + 6)A_{31}B_{31} \Chi^2 + (c_0 + 9)A_{31}B_{31} C_{31}\Chi^3 }{1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 } </math>

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

<math>~=</math>

<math>~ [c_0 + (c_0 + 3)A_{31}\Chi + (c_0 + 6)A_{31}B_{31} \Chi^2 + (c_0 + 9)A_{31}B_{31} C_{31}\Chi^3 ] - Q_{31}[1 + A_{31}\Chi + A_{31}B_{31}\Chi^2 + A_{31}B_{31}C_{31}\Chi^3 ] </math>

 

<math>~=</math>

<math>~ \biggl[c_0 - (c_0 + 3)\cdot \frac{6(c_0 + 7)}{2c_0 + 5 } \cdot \Chi + (c_0 + 6) \cdot \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{ (6c_0 + 51)}{ 6(c_0 +4) } \cdot \Chi^2 - (c_0 + 9)\cdot \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{ (6c_0 + 51)}{ 6(c_0 +4) } \cdot \frac{ 2( c_0 + 10)}{3( 2c_0 + 11) } \cdot\Chi^3 \biggr] </math>

 

 

<math>~ - Q_{31} \biggl[1 - \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \Chi + \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{ (6c_0 + 51)}{ 6(c_0 +4) } \cdot \Chi^2 - \frac{ 6(c_0 + 7)}{2c_0 + 5 } \cdot \frac{ (6c_0 + 51)}{ 6(c_0 +4) } \cdot \frac{ 2( c_0 + 10)}{3( 2c_0 + 11) } \cdot \Chi^3 \biggr] </math>

 

<math>~=</math>

<math>~ \biggl[3c_0(2c_0 + 5)(c_0 +4)( 2c_0 + 11) - (c_0 + 3)(c_0 +4)3( 2c_0 + 11)\cdot 6(c_0 + 7) \cdot \Chi + (c_0 + 6) 3( 2c_0 + 11)\cdot (c_0 + 7) \cdot (6c_0 + 51) \cdot \Chi^2 - (c_0 + 9)\cdot (c_0 + 7) \cdot (6c_0 + 51)\cdot 2( c_0 + 10) \cdot\Chi^3 \biggr] </math>

 

 

<math>~ - Q_{31} \biggl[ 3(2c_0 + 5)(c_0 +4)( 2c_0 + 11) - 18(c_0 + 7)(c_0 +4)( 2c_0 + 11) \cdot \Chi + 3( 2c_0 + 11)(c_0 + 7) \cdot (6c_0 + 51) \cdot \Chi^2 - 2(c_0 + 7) \cdot (6c_0 + 51) ( c_0 + 10) \cdot \Chi^3 \biggr] </math>

 

<math>~=</math>

<math>~ \biggl[(6c_0^2 + 15c_0) (2c_0^2 + 19c_0 + 44) - 18(c_0^2 + 10c_0 + 21) (2c_0^2 + 19c_0 + 44) \cdot \Chi + 3 (2c_0^2 + 23c_0 +66) (6c_0^2 + 93c_0 + 357) \cdot \Chi^2 - 2(c_0^2 + 19c_0 + 90) (6c_0^2 + 93c_0 + 357) \cdot\Chi^3 \biggr] </math>

 

 

<math>~ - Q_{31} \biggl[ 3(2c_0 + 5) (2c_0^2 + 19c_0 + 44) - 18(c_0 + 7) (2c_0^2 + 19c_0 + 44) \cdot \Chi + 3( 2c_0 + 11) (6c_0^2 + 93c_0 + 357) \cdot \Chi^2 - 2 (6c_0^2 + 93c_0 + 357) ( c_0 + 10) \cdot \Chi^3 \biggr] </math>

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

<math>~=</math>

<math>~ \biggl[(12c_0^4 + 144c_0^3 + 549c_0^2 + 660c_0) - 18(2c_0^4 + 39c_0^3 + 276c_0^2 + 839c_0 + 924) \cdot \Chi + 3 (12c_0^4 + 324c_0^3 + 3249c_0^2 + 14349c_0 + 23562) \cdot \Chi^2 - 2(6c_0^4 + 207c_0^3 + 2664c_0^2 + 15153c_0 + 32130) \cdot\Chi^3 \biggr] </math>

 

 

<math>~ - Q_{31} \biggl[ 3(4c_0^3 + 48c_0^2 + 183c_0 +220) - 18(2c_0^3 + 33c_0^2 + 177c_0 + 308) \cdot \Chi + 3( 12c_0^3 + 252c_0^2 + 1737c_0 + 3927) \cdot \Chi^2 - 2 (6c_0^3 + 153c_0^2 + 1287c_0 + 3570) \cdot \Chi^3 \biggr] </math>

Root of Quartic Equation

To solve this equation analytically, we follow the Summary of Ferrari's method that is presented in Wikipedia's discussion of the Quartic Function to identify the roots of an arbitrary quartic equation.

First, we adopt the shorthand notation:

<math>~0</math>

<math>~=</math>

<math>~ac_0^4 + bc_0^3 + c c_0^2 +d c_0 +e \, ,</math>

where, in our particular case,

<math>~e</math>

<math>~\equiv</math>

<math>~ - 18\cdot 924\Chi +3\cdot 23562\Chi^2 -2\cdot 32130\Chi^3 -Q_{31}[ 3\cdot 220 - 18\cdot 308\Chi +3\cdot 3927\Chi^2 -2\cdot 3570\Chi^3 ] \, ,</math>

<math>~d</math>

<math>~\equiv</math>

<math>~660 -18\cdot 839\Chi +3\cdot 14349\Chi^2 -2\cdot 15153\Chi^3 -Q_{31}[3\cdot 183 -18\cdot 177\Chi + 3\cdot 1737\Chi^2 -2\cdot 1287\Chi^3 ] \, ,</math>

<math>~c</math>

<math>~\equiv</math>

<math>~549 - 18\cdot 276\Chi +3\cdot 3249\Chi^2 -2\cdot 2664\Chi^3 -Q_{31}[ 3\cdot 48 - 18\cdot 33\Chi + 3\cdot 252\Chi^2 -2\cdot 153 \Chi^3 ] \, ,</math>

<math>~b</math>

<math>~\equiv</math>

<math>~144- 18\cdot 39\Chi +3\cdot 324\Chi^2 -2\cdot 207\Chi^3 - 12 Q_{31}[ 1 - 3\Chi +3\Chi^2 - \Chi^3 ]\, ,</math>

<math>~a</math>

<math>~\equiv</math>

<math>~12( 1 - 3 \Chi +3\Chi^2 -\Chi^3) \, .</math>

Now, define,

<math>~\Delta_0</math>

<math>~\equiv</math>

<math>~ c^2 - 3bd + 12ae \, ,</math>

<math>~\Delta_1</math>

<math>~\equiv</math>

<math>~ 2c^3 - 9bcd + 27b^2e + 27ad^2 - 72ace \, ,</math>

<math>~p</math>

<math>~\equiv</math>

<math>~ \frac{8ac - 3b^2}{8a^2} \, ,</math>

<math>~\kappa</math>

<math>~\equiv</math>

<math>~ \frac{b^3 - 4abc + 8a^2 d}{8a^3} \, ,</math>

<math>~\Kappa</math>

<math>~\equiv</math>

<math>~ \frac{1}{2^{1 / 3}} \biggl[\Delta_1 + \sqrt{\Delta_1^2 - 4\Delta_0^3} \biggr]^{1 / 3} \, ,</math>

(see below)   <math>~\biggl(\Kappa + \frac{\Delta_0}{\Kappa}\biggr)</math>

<math>~\equiv</math>

<math>~2\Delta_0^{1/2}\cos\biggl[ \frac{1}{3} \cos^{-1}\biggl( \frac{\Delta_1^2}{4\Delta_0^3}\biggr)^{1/2}\biggr] \, ,</math>

<math>~S</math>

<math>~\equiv</math>

<math>~ \frac{1}{2} \biggl[ - \frac{2p}{3} + \frac{1}{3a}\biggl(\Kappa + \frac{\Delta_0}{\Kappa} \biggr) \biggr]^{1 / 2} \, .</math>

Then the four roots of the quartic equation are,

<math>~(c_0)_{1}</math>

<math>~=</math>

<math>~-\frac{b}{4a} - S + \frac{1}{2}\biggl[ -4S^2 - 2p + \frac{\kappa}{S} \biggr]^{1 / 2} \, ,</math>

<math>~(c_0)_{2}</math>

<math>~=</math>

<math>~-\frac{b}{4a} - S - \frac{1}{2}\biggl[ -4S^2 - 2p + \frac{\kappa}{S} \biggr]^{1 / 2} \, ,</math>

<math>~(c_0)_{3}</math>

<math>~=</math>

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

<math>~(c_0)_{4}</math>

<math>~=</math>

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

It is this fourth root that interests us, here.


We have determined empirically that, in our specific case, the quantity,

<math>~\Delta_1^2 - 4\Delta_0^3</math>

is negative over the range of physically interesting values of <math>~\Chi</math>. Hence, the quantity, <math>~\Kappa^3</math>, is necessarily complex. Let's work carefully through a determination of <math>~\Kappa</math> and, by consequence, <math>~S</math>, in this situation.

<math>~2\Kappa^3</math>

<math>~=</math>

<math>~\Delta_1 + \sqrt{\Delta_1^2 - 4\Delta_0^3} </math>

 

<math>~=</math>

<math>~\Delta_1 + i \sqrt{4\Delta_0^3 - \Delta_1^2} </math>

 

<math>~=</math>

<math>~\Delta_1 + 2\Delta_0^{3/2} ~ i \biggl[ 1 - \Gamma^2 \biggr]^{1/2} </math>

<math>~\Rightarrow ~~~ \biggl[ \frac{\Kappa}{\Delta_0^{1/2}} \biggr]^3</math>

<math>~=</math>

<math>~\Gamma + i \sqrt{ 1 - \Gamma^2 } \, , </math>

where,

<math>\Gamma \equiv \biggl[ \frac{\Delta_1^2}{4\Delta_0^3}\biggr]^{1/2} \, .</math>

We therefore can state that, in the complex plane, the three roots <math>~(j=0,1,3)</math> of this expression are,

<math>~ \frac{\Kappa}{\Delta_0^{1/2}} </math>

<math>~=</math>

<math>~e^{i\theta_\Kappa/3} \cdot e^{i(2\pi j/3)} \, ,</math>

where,

<math>\theta_\Kappa \equiv \cos^{-1}\Gamma \, .</math>

Focusing on the simplest <math>~(j=0)</math> root, we have,

<math>~ \Kappa </math>

<math>~=</math>

<math>~\Delta_0^{1/2} e^{i\theta_\Kappa/3} </math>

<math>~ \Rightarrow ~~~ \biggl(\Kappa + \frac{\Delta_0}{\Kappa} \biggr)</math>

<math>~=</math>

<math>~\Delta_0^{1/2} e^{i\theta_\Kappa/3} + \Delta_0^{1/2} e^{- i\theta_\Kappa/3} </math>

 

<math>~=</math>

<math>~2\Delta_0^{1/2}\cos\biggl[ \frac{\cos^{-1}\Gamma}{3} \biggr] \, .</math>

Because this expression does not contain an imaginary component, we understand that <math>~S</math> is real.


Finally, as explained in our summary discussion, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,

<math>~ \gamma_c </math>

<math>~=</math>

<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{ 20} \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, , </math>

where, the last expression follows from plugging in the desired mode's quantum numbers, <math>~(\ell,j) = (3,1)</math>, and, again, using the relation,

<math>~\Chi = q^3</math>

<math>~=</math>

<math>~\frac{\rho_e/\rho_c}{2(1-\rho_e/\rho_c)}</math>

<math>~\Rightarrow ~~~ \frac{\rho_e}{\rho_c}</math>

<math>~=</math>

<math>~\frac{2\Chi}{1+2\Chi} \, .</math>

Illustration31

Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes

Quantum Numbers:  <math>~(\ell,j) = (3,1)</math>

Analyzable Model Sequence One Example Eigenfunction

Montage of Stability Results for (ell,j) = (3,1) quantum numbers


Top-Left Panel:  Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> — where <math>~\gamma_e</math> is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (3,1)</math>. Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>. The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.

Bottom-Left Panel:   Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>. Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line. The yellow circular marker identifies the same example model as it does in the top-left panel.

Top-Right Panel:  Displays — as a function of the fractional radius, <math>~r_0/R = q\xi</math> — the analytically determined eigenfunction for the <math>~(\ell,j) = (3,1)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,

<math>~q = \biggl[0.01 + 40\biggl( \frac{0.98}{99} \biggr) \biggr] \approx 0.4059596</math>

and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (-1.7819827, -0.3885031, -0.9647648)</math>. Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction, namely,

<math>~x_{j=1} |_\mathrm{core}</math>

<math>~=</math>

<math>~ \frac{5(1+8q^3) - 7 (1+2q^3)^2 \xi^2}{5(1+8q^3)-7(1+2q^3)^2} \, ;</math>

and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction, namely,

<math>~x_{\ell=3} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{31} \xi^{3} + q^6 A_{31}B_{31}\xi^{6} + q^9 A_{31}B_{31}C_{31}\xi^{9} }{ 1 + q^3 A_{31} + q^6 A_{31}B_{31} + q^9 A_{31}B_{31}C_{31} }\biggr] \, , </math>

where the coefficients, <math>~A_{31}, B_{31}, C_{31}</math>, are as defined above in terms of the parameter, <math>~c_0</math>. The corresponding eigenfrequency is, from the perspective of the core,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~20\gamma_c - 8 = \frac{8(7 + \alpha_c)}{3-\alpha_c} \approx 12.17774\, ;</math>

and, from the perspective of the envelope,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~\gamma_e\biggl[2\alpha_e + (c_0 + 9)(c_0+14)\biggr] \frac{\rho_e}{\rho_c} \approx 12.17774\, ,</math>

where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.118016</math>.

Bottom-Right Panel:  The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,

<math>~y</math>

<math>~=</math>

<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>

where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the outer edge. In this type of log-amplitude plot, the eigenfunction's various nodes — that is, radial locations where <math>~x</math> passes through zero — are highlighted; here, specifically, there is one node inside the core and two nodes reside in the envelope. Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope. Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.41)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.12)</math>.

NOTE: As may be ascertained from our general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes, equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to any value, <math>~0 \le \rho_e/\rho_c \le 1</math>. We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,

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

<math>~=</math>

<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>

The black-dotted curve in the Bottom-Right Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radial-oscillation eigenvector.

Check Surface Boundary Condition

Let's determine analytically the logarithmic derivative of the envelope segment of the eigenfunction, at the surface <math>~(\xi = 1/q)</math> of the configuration. Given that,

<math>~x_{\ell=3} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{3} \xi^{3} + q^6 A_{3}B_{3}\xi^{6} + q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3} + q^6 A_{3}B_{3} + q^9 A_{3}B_{3}C_{3} }\biggr] \, , </math>

where,

<math>~A_{3}</math>

<math>~=</math>

<math>~\frac{- 6(c_0 + 7)}{(2c_0 + 5 )} \, ,</math>

<math>~B_{3}</math>

<math>~=</math>

<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>

<math>~C_{3}</math>

<math>~=</math>

<math>~\frac{- 2( c_0 + 10)}{3(2 c_0 + 11) } \, ,</math>

we have,

<math>~\frac{ d\ln x_{\ell=3} }{d\ln \xi}</math>

<math>~=</math>

<math>~\frac{1}{x} \biggl\{ c_0\xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{3} \xi^{3} + q^6 A_{3}B_{3}\xi^{6} + q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3} + q^6 A_{3}B_{3} + q^9 A_{3}B_{3}C_{3} }\biggr] + \xi^{c_0}\biggl[ \frac{ 3 q^3 A_{3} \xi^{3} + 6q^6 A_{3}B_{3}\xi^{6} + 9 q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3} + q^6 A_{3}B_{3} + q^9 A_{3}B_{3}C_{3} }\biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ c_0 + \biggl[ \frac{ 3 q^3 A_{3} \xi^{3} + 6q^6 A_{3}B_{3}\xi^{6} + 9 q^9 A_{3}B_{3}C_{3}\xi^{9} }{ 1 + q^3 A_{3} \xi^{3} + q^6 A_{3}B_{3}\xi^{6} + q^9 A_{3}B_{3}C_{3}\xi^{9} }\biggr] \, . </math>

This means that,

<math>~\frac{ d\ln x_{\ell=3} }{d\ln \xi} \biggr|_{\xi=1/q}</math>

<math>~=</math>

<math>~ c_0 + \biggl[ \frac{ 3 A_{3} + 6A_{3}B_{3} + 9 A_{3}B_{3}C_{3} }{ 1 + A_{3} + A_{3}B_{3}+ A_{3}B_{3}C_{3} }\biggr] </math>

 

<math>~=</math>

<math>~ c_0 + \biggl\{ \frac{ -3 [108(c_0 + 7)(c_0 +4)(2 c_0 + 11)] + 6[18(c_0 + 7) (6c_0 + 51)(2 c_0 + 11)] - 9 [ 12(c_0 + 7)(6c_0 + 51)( c_0 + 10) ] }{ 18(2c_0+5)(c_0+4)(2c_0+11) - 108(c_0 + 7)(c_0 +4)(2 c_0 + 11) + 18(c_0 + 7) (6c_0 + 51)(2 c_0 + 11) - 12(c_0 + 7)(6c_0 + 51)( c_0 + 10) }\biggr\} </math>

 

<math>~=</math>

<math>~ c_0 + 108\biggl\{ \frac{ - 3(2c_0^3 + 33c_0^2 + 177c_0 + 308) + (12c_0^3 + 252c_0^2 + 1737c_0 + 3927) - (6c_0^3 + 153c_0^2 + 1287c_0 + 3570) }{ 18(4c_0^3 + 48c_0^2 + 183c_0 +220) - 108(2c_0^3 + 33c_0^2 + 177c_0 + 308) + 18(12c_0^3 + 252c_0^2 + 1737c_0 + 3927) - 12(6c_0^3 + 153c_0^2 + 1287c_0 + 3570) }\biggr\} </math>

 

<math>~=</math>

<math>~ c_0 + 108\biggl\{ \frac{0 c_0^3 + 0 c_0^2 -81 c_0 -567 }{ 0 c_0^3 + 0c_0^2 + 0c_0 -1458 }\biggr\} </math>

 

<math>~=</math>

<math>~ c_0 + 6(c_0 + 7) </math>

 

<math>~=</math>

<math>~ 7(c_0 + 6) \, , </math>

which precisely matches the desired surface boundary condition for <math>~\ell=3</math> that has been detailed in an accompanying summary discussion.



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   |


Try Mode32 <math>~(\ell,j) = (3,2)</math>

Setup32

Here, we set the quantum numbers to <math>~(\ell,j) = (3,2)</math>. In this case, as when addessing Mode22, above, the core eigenfunction segment should be,

<math>~x_{j=2} |_\mathrm{core}</math>

<math>~=</math>

<math>~ \frac{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 \xi^2 + 99 (1+2q^3)^4 \xi^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } \, ;</math>

and, as when addressing Mode31, above, throughout the envelope, the eigenfunction should be,

<math>~x_{\ell=3} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{32} \xi^{3} + q^6 A_{32}B_{32}\xi^{6} + q^9 A_{32}B_{32}C_{32}\xi^{9} }{ 1 + q^3 A_{32} + q^6 A_{32}B_{32} + q^9 A_{32}B_{32}C_{32} }\biggr] \, , </math>

where,

<math>~A_{32}</math>

<math>~=</math>

<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>

<math>~B_{32}</math>

<math>~=</math>

<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>

<math>~C_{32}</math>

<math>~=</math>

<math>~\frac{- 2( c_0 + 10)}{3( 2c_0 + 11) } \, .</math>

In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace <math>~Q_{31}</math> with,

<math>~Q_{32}</math>

<math>~\equiv</math>

<math>~ \frac{ - 252(1+8q^3) (1+2\Chi)^2 + 396 (1+2\Chi)^4 }{35(1+8q^3)^2 - 126(1+8q^3) (1+2q^3)^2 + 99 (1+2q^3)^4 } \, . </math>

And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,

<math>~ \gamma_c </math>

<math>~=</math>

<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{42} \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, . </math>


Illustration32

Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes

Quantum Numbers:  <math>~(\ell,j) = (3,2)</math>

Analyzable Model Sequence One Example Eigenfunction

Montage of Stability Results for (ell,j) = (3,2) quantum numbers


Top-Left Panel:  Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> — where <math>~\gamma_e</math> is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (3,2)</math>. Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>. The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.

Bottom-Left Panel:   Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>. Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line. The yellow circular marker identifies the same example model as it does in the top-left panel.

Top-Right Panel:  Displays — as a function of the fractional radius, <math>~r_0/R = q\xi</math> — the analytically determined eigenfunction for the <math>~(\ell,j) = (3,2)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,

<math>~q = \biggl[0.01 + 185\biggl( \frac{0.98}{3\times 99} \biggr)\biggr] \approx 0.6204377</math>

and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (-2.3774068, + 0.8972496, +0.0231304)</math>. Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction as defined by the expression for <math>~x_{j=2} |_\mathrm{core}</math>, given above; and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction as defined by the expression for <math>~x_{\ell=3} |_\mathrm{env}</math>, given above.

The corresponding eigenfrequency is, from the perspective of the core,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~42\gamma_c - 8 = \frac{8(18 + \alpha_c)}{3-\alpha_c} \approx 48.4351\, ;</math>

and, from the perspective of the envelope,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~\gamma_e\biggl[2\alpha_e + (c_0 + 9)(c_0+14)\biggr] \frac{\rho_e}{\rho_c} \approx 48.4351 \, ,</math>

where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.323257</math>.

Bottom-Right Panel:  The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,

<math>~y</math>

<math>~=</math>

<math>~\tfrac{1}{8} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>

where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the enter of the configuration. In this type of log-amplitude plot, the eigenfunction's various nodes — that is, radial locations where <math>~x</math> passes through zero — are highlighted; here, specifically, there is one node inside the core and one node resides in the envelope. Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope. Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.62)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.32)</math>.

NOTE: As may be ascertained from our general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes, equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to any value, <math>~0 \le \rho_e/\rho_c \le 1</math>. We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,

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

<math>~=</math>

<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>

The black-dotted curve in the Bottom-Right Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radial-oscillation eigenvector.

Try Mode33 <math>~(\ell,j) = (3,3)</math>

Setup33

Here, we set the quantum numbers to <math>~(\ell,j) = (3,3)</math>. In this case, drawing from our summary of Sterne's presentation, the core eigenfunction segment should be,

<math>~x_{j=3} |_\mathrm{core}</math>

<math>~=</math>

<math>~ \frac{3\cdot 5 \cdot 7 - 33\cdot 3\cdot 7 (1+2q^3)^2 \xi^2 + 429\cdot 3 (1+2q^3)^4 \xi^4 - 143\cdot 5 (1+2q^3)^6 \xi^6}{3\cdot 5 \cdot 7 - 33\cdot 3\cdot 7 (1+2q^3)^2 + 429\cdot 3 (1+2q^3)^4 - 143\cdot 5 (1+2q^3)^6 } \, ,</math>

and, as when addressing Mode31, above, throughout the envelope, the eigenfunction should be,

<math>~x_{\ell=3} |_\mathrm{env}</math>

<math>~=</math>

<math>~ \xi^{c_0}\biggl[ \frac{ 1 + q^3 A_{33} \xi^{3} + q^6 A_{33}B_{33}\xi^{6} + q^9 A_{33}B_{33}C_{33}\xi^{9} }{ 1 + q^3 A_{33} + q^6 A_{33}B_{33} + q^9 A_{33}B_{33}C_{33} }\biggr] \, , </math>

where,

<math>~A_{33}</math>

<math>~=</math>

<math>~\frac{- 6(c_0 + 7)}{2c_0 + 5 } \, ,</math>

<math>~B_{33}</math>

<math>~=</math>

<math>~\frac{ - (6c_0 + 51)}{ 6(c_0 +4) } \, ,</math>

<math>~C_{33}</math>

<math>~=</math>

<math>~\frac{- 2( c_0 + 10)}{7( c_0 + 6) } \, .</math>

In establishing the interface matching conditions, we should be able to follow the same steps as in our search for "31" modes, but replace <math>~Q_{31}</math> with,

<math>~Q_{33}</math>

<math>~\equiv</math>

<math>~ \frac{ - 2\cdot 33\cdot 3\cdot 7 (1+2q^3)^2 + 4\cdot 429\cdot 3 (1+2q^3)^4 - 6\cdot 143\cdot 5 (1+2q^3)^6 }{3\cdot 5 \cdot 7 - 33\cdot 3\cdot 7 (1+2q^3)^2 + 429\cdot 3 (1+2q^3)^4 - 143\cdot 5 (1+2q^3)^6 } \, . </math>


And, in order for the two (dimensional) eigenfrequencies to match, the adiabatic exponent for the core must be,

<math>~ \gamma_c </math>

<math>~=</math>

<math>~\frac{1}{ 6 + 2j(2j+5)] } \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell +5) \biggr]\frac{\rho_e}{\rho_c} \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{72} \biggl\{ 8 + \gamma_e \biggl[2\alpha_e + (c_0 + 9)(c_0 + 14) \biggr]\frac{\rho_e}{\rho_c} \biggr\} \, . </math>


Figure 4: Results for mode <math>~(\ell,j) = (3,3)</math>

[[|850px|Results for (ell,j) = (3,3)]]

<math>~q</math> <math>~ </math>
<math>~c_0</math> <math>~ </math>
<math>~\gamma_e</math> <math>~ </math>
<math>~\alpha_e</math> <math>~0. </math>
<math>~n_e</math> <math>~ </math>
<math>~\nu</math> <math>~ </math>
<math>~\rho_e/\rho_c</math> <math>~ </math>
<math>~\gamma_c</math> <math>~ </math>
<math>~\alpha_c</math> <math>~ </math>

Illustration33

Analytically Definable Eigenvectors in <math>~(n_c, n_e) = (0,0)</math> Bipolytropes

Quantum Numbers:  <math>~(\ell,j) = (3,3)</math>

Analyzable Model Sequence One Example Eigenfunction

Montage of Stability Results for (ell,j) = (3,3) quantum numbers


Top-Left Panel:  Plotted points show how the location of the core/envelope interface, <math>~q \equiv r_i/R</math>, varies with <math>~\alpha_e \equiv (3-4/\gamma_e)</math> — where <math>~\gamma_e</math> is the adiabatic exponent of the envelope — in equilibrium models that are amenable to analytic modal analysis for quantum numbers, <math>~(\ell,j) = (3,3)</math>. Red (alternatively, blue) markers identify models for which the corresponding value of the adiabatic exponent of the core [see bottom-left panel] falls inside (alternatively, outside) the physically viable range, namely, <math>~1 \le \gamma_c \le \infty</math>. The yellow circular marker identifies the model whose analytically determined eigenfunction is displayed on the right, as an example.

Bottom-Left Panel:   Plotted points show how <math>~\alpha_c \equiv (3-4/\gamma_c)</math> varies with <math>~\alpha_e</math> over the physically viable parameter range, <math>~-1 \le \alpha \le 3</math>. Both axes have been flipped so that incompressible models <math>~(\gamma = \infty)</math> lie on the left/bottom while isothermal models <math>~(\gamma =1)</math> lie on the right/top. The core is more compressible than the envelope in models that lie above and to the left of the black-dashed, diagonal line. The yellow circular marker identifies the same example model as it does in the top-left panel.

Top-Right Panel:  Displays — as a function of the fractional radius, <math>~r_0/R = q\xi</math> — the analytically determined eigenfunction for the <math>~(\ell,j) = (3,2)</math> mode in the model identified by the yellow circular marker in both left-hand panels, for which,

<math>~q = 0.01 + 642\biggl( \frac{0.98}{990} \biggr) \approx 0.6455152 \, ,</math>

and, correspondingly, <math>~(c_0, \alpha_e,\alpha_c) = (+0.3511525, +0.8256130, -0.0183637)</math>. Specifically, over the radial interval, <math>~0 \le \xi \le 1</math>, the green markers trace the core's contribution to the combined eigenfunction as defined by the expression for <math>~x_{j=3} |_\mathrm{core}</math>, given above; and, over the radial interval, <math>~1 \le \xi \le 1/q</math>, the purple markers trace the envelope's contribution to the combined eigenfunction as defined by the expression for <math>~x_{\ell=3} |_\mathrm{env}</math>, given above.

The corresponding eigenfrequency is, from the perspective of the core,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~72\gamma_c - 8 = \frac{8(33 + \alpha_c)}{3-\alpha_c} \approx 87.41594 \, ;</math>

and, from the perspective of the envelope,

<math>~\frac{3\omega^2}{2\pi G \rho_c}</math>

<math>~=</math>

<math>~\gamma_e\biggl[2\alpha_e + (c_0 + 9)(c_0+14)\biggr] \frac{\rho_e}{\rho_c} \approx 87.41594 \, ,</math>

where the relevant density ratio is, <math>~\rho_e/\rho_c = 2q^3/(1+2q^3) \approx 0.3497877</math>.

Bottom-Right Panel:  The green and purple markers present the same eigenfunction-amplitude information, <math>~x(r/R)</math>, as in the Top-Right panel, but on a logarithmic scale. Specifically, in this plot, the vertical displacement of the green and purple markers is given by the expression,

<math>~y</math>

<math>~=</math>

<math>~\tfrac{1}{10} \log_{10}[x^2 + \epsilon^2] + y_\mathrm{shift} \, ,</math>

where, for plotting purposes, we have used, <math>~\epsilon = 10^{-5}</math>, and have set <math>~y_\mathrm{shift}</math> to a value that ensures that <math>~y \approx 1</math> at the outer edge. In this type of log-amplitude plot, the eigenfunction's various nodes — that is, radial locations where <math>~x</math> passes through zero — are highlighted; here, specifically, there are three nodes inside the core and one node resides in the envelope. Using the vertical coordinate to represent, instead, the configuration's mass-density normalized to its central value, <math>~\rho/\rho_c</math>, the solid black line segments trace the unperturbed density distribution throughout this specific <math>~(n_c, n_e) = (0,0)</math> bipolytrope. Throughout the core, <math>~\rho/\rho_c = 1</math>; then, at the location of the interface <math>~(r_i/R = q \approx 0.65)</math>, the density abruptly drops to its envelope value <math>~(\rho/\rho_c = \rho_e/\rho_c \approx 0.35)</math>.

NOTE: As may be ascertained from our general discussion of the structural properties of <math>~(n_c, n_e) = (0,0)</math> bipolytropes, equilibrium "zero-zero" bipolytropes can be constructed with the envelope/core interface parameter set to any value across the range, <math>~0 \le q \le 1</math>; and for any chosen value of <math>~q</math>, the envelope/core density ratio can, in principle, be set to any value, <math>~0 \le \rho_e/\rho_c \le 1</math>. We have not, however, been able to analytically solve the relevant pair of linear-adiabatic wave equations (LAWEs) for this entire set of equilibrium models. Instead, our ability to derive analytically prescribed eigenvectors is limited by the constraint,

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

<math>~=</math>

<math>~\frac{2q^3}{1+2q^3} = \frac{2(r_i/R)^3}{1+2(r_i/R)^3}\, .</math>

The black-dotted curve in the Bottom-Right Panel displays the behavior of this constraint; accordingly, the step function depicted by the solid black line segments must necessarily drop from unity to a point on this black-dotted curve for any equilibrium model — such as the example illustrated here — that has an analytically prescribable radial-oscillation eigenvector.

Summary Plots

<math>~(\ell, j) = (2,1)</math> <math>~(\ell, j) = (2,2)</math>  
Log(amplitude) plot for (ell,j) = (2,1) Log(amplitude) plot for (ell,j) = (2,2)  
<math>~\frac{3\omega^2}{2\pi G\rho_c} = 37.08874</math> <math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.95210</math>  
 
<math>~(\ell, j) = (3,1)</math> <math>~(\ell, j) = (3,2)</math> <math>~(\ell, j) = (3,3)</math>
Log(amplitude) plot for (ell,j) = (3,1) Log(amplitude) plot for (ell,j) = (3,2) Log(amplitude) plot for (ell,j) = (3,3)
<math>~\frac{3\omega^2}{2\pi G\rho_c} = 12.452545</math> <math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.05461</math> <math>~\frac{3\omega^2}{2\pi G\rho_c} = 87.41594</math>

Related Discussions

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