Editing SSC/Stability/BiPolytrope00Old

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

=Radial Oscillations of a Zero-Zero Bipolytrope=
This is a chapter that summarizes [[SSC/Stability/BiPolytrope00Details#Radial_Oscillations_of_a_Zero-Zero_Bipolytrope|an accompanying, detailed derivation]].  
A newer version of this chapter, that takes into account a ''corrected'' interface matching condition provided by [http://adsabs.harvard.edu/abs/1958HDP....51..353L P. Ledoux &amp; Th. Walraven (1958)] is provided [[SSC/Stability/BiPolytrope00#Radial_Oscillations_of_a_Zero-Zero_Bipolytrope|here.]]

In a [[SSC/Structure/BiPolytropes/Analytic00#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|separate chapter on astrophysical interesting ''equilibrium structures'']], we have derived analytical expressions that define the equilibrium properties of bipolytropic configurations having <math>(n_c, n_e) = (0, 0)</math>, that is, bipolytropes in which both the core and the envelope are uniform in density, but the densities in the two regions are different from one another.  Letting <math>R</math> be the radius and <math>M_\mathrm{tot}</math> be the total mass of the bipolytrope, these configurations are fully defined once any two of the following three key parameters have been specified:  The envelope-to-core density ratio, <math>\rho_e/\rho_c</math>; the radial location of the envelope/core interface, <math>~q \equiv r_i/R</math>; and, the fractional mass that is contained within the core, <math>\nu \equiv M_\mathrm{core}/M_\mathrm{tot}</math>.  These three parameters are related to one another via the expression,
<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{q^3}{\nu} \biggl( \frac{1-\nu}{1-q^3} \biggr) \, .</math>
  </td>
</tr>
</table>
</div>

Equilibrium configurations can be constructed that have a wide range of parameter values; specifically,
<div align="center">
<math>~0 \le q \le 1 \, ;</math>
&nbsp; &nbsp; &nbsp; &nbsp; 
<math>~0 \le \nu \le 1 \, ;</math>
&nbsp; &nbsp; &nbsp; &nbsp; 
and,
&nbsp; &nbsp; &nbsp; &nbsp; 
<math>~0 \le \frac{\rho_e}{\rho_c} \le 1 \, .</math>
</div>
(We recognize from buoyancy arguments that any configuration in which the envelope density is larger than the core density will be Rayleigh-Taylor unstable, so we restrict our astrophysical discussion to structures for which <math>~\rho_e < \rho_c</math>.)


By employing the [[SSC/Perturbations#Spherically_Symmetric_Configurations_.28Stability_.E2.80.94_Part_II.29|linear stability analysis techniques described in an accompanying chapter]], we should, in principle, be able to identify a wide range of eigenvectors &#8212; that is, radial eigenfunctions and accompanying eigenfrequencies &#8212; that are associated with adiabatic radial oscillation modes in any one of these equilibrium, bipolytropic configurations.  Using numerical techniques, [http://adsabs.harvard.edu/abs/1985PASAu...6..222M Murphy &amp; Fiedler (1985)], for example, have carried out such an analysis of bipolytropic structures having <math>~(n_c, n_e) = (1,5)</math>.  A ''pair'' of  [[SSC/Perturbations#2ndOrderODE|linear adiabatic wave equations (LAWEs)]] must be solved &#8212; one tuned to accommodate the properties of the core and another tuned to accommodate the properties of the envelope &#8212;  then the pair of eigenfunctions must be matched smoothly at the radial location of the interface; the identified core- and envelope-eigenfrequencies must simultaneously match.  


After identifying the precise form of the LAWEs that apply to the case of <math>~(n_c, n_e) = (0,0)</math> bipolytropes, we discovered that, for a restricted range of key parameters, the pair of equations can both be solved ''analytically''. 

==Two Separate LAWEs==

In an [[User:Tohline/SSC/Perturbations#2ndOrderODE|accompanying discussion]], we derived the so-called,

<div align="center" id="2ndOrderODE">
<font color="#770000">'''Adiabatic Wave''' (or ''Radial Pulsation'') '''Equation'''</font><br />

{{User:Tohline/Math/EQ_RadialPulsation01}}
</div>

<!--
<div align="center" id="2ndOrderODE">
<font color="#770000">'''Adiabatic Wave Equation'''</font><br />

<math>
\frac{d^2x}{dr_0^2} + \biggl[\frac{4}{r_0} - \biggl(\frac{g_0 \rho_0}{P_0}\biggr) \biggr] \frac{dx}{dr_0} + \biggl(\frac{\rho_0}{\gamma_\mathrm{g} P_0} \biggr)\biggl[\omega^2 + (4 - 3\gamma_\mathrm{g})\frac{g_0}{r_0} \biggr]  x = 0 \, ,
</math>
</div>
-->
For both regions of the bipolytrope, we define the dimensionless (Lagrangian) radial coordinate,
<div align="center">
<math>~\xi  \equiv \frac{r_0}{r_i} \, .</math>
</div>
So, the interface is, by definition, located at <math>~\xi = 1</math>; and, the surface is necessarily at <math>~\xi = q^{-1}</math>.  As the material in the bipolytrope's core (envelope) is compressed/de-compressed during a radial oscillation, we will assume that heating/cooling occurs in a manner prescribed by an adiabat of index <math>~\gamma_c ~(\gamma_e)</math>; in general, <math>~\gamma_e \ne \gamma_c</math>.  For convenience, we will also adopt the frequently used shorthand "alpha" notation,
<div align="center">
<math>~\alpha_c \equiv 3 - \frac{4}{\gamma_c} \, ,</math>
&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp;
<math>~\alpha_e \equiv 3 - \frac{4}{\gamma_e} \, .</math>
</div>


===The Core's LAWE===
After adopting, for convenience, the function notation,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~g^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>
1  + \biggl(\frac{\rho_e}{\rho_c}\biggr)  \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_c} \biggr) \biggl( 1-q \biggr) + 
\frac{\rho_e}{\rho_c} \biggl(\frac{1}{q^2} - 1\biggr) \biggr]  \, ,
</math>
  </td>
</tr>
</table>
</div>
we [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Match_Prasad-like_Envelope_Eigenvector_to_the_Core_Eigenvector|have deduced]] that, for the core, the LAWE may be written in the form,

<div align="center">
<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>~
(1 - \eta^2)\frac{d^2x}{d\eta^2} +  
( 4 - 6\eta^2 )  \frac{1}{\eta} \cdot \frac{dx}{d\eta} 
+ \mathfrak{F}_\mathrm{core} x \, .
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<math>~\eta \equiv \frac{\xi}{g} \, ,</math>
&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp;
<math>~\mathfrak{F}_\mathrm{core} \equiv \frac{3\omega_\mathrm{core}^2}{2\pi G\gamma_c \rho_c} - 2\alpha_c\, .</math>
</div>
Not surprisingly, this is identical in form to the eigenvalue problem that was first presented &#8212; and solved analytically &#8212; by [[User:Tohline/SSC/UniformDensity#Setup_as_Presented_by_Sterne_.281937.29|Sterne (1937)]] in connection with his examination of radial oscillations in ''isolated'' uniform-density spheres.  As is demonstrated below, for the core of our zero-zero bipolytrope, we can in principle adopt any one of the [[User:Tohline/SSC/UniformDensity#Sterne.27s_General_Solution|polynomial eigenfunctions and corresponding eigenfrequencies]] derived by Sterne.  


===The Envelope's LAWE===

Subsequently, we also [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#More_General_Solution|have deduced]] that, for the envelope, the governing LAWE becomes,

<div align="center">
<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>~
\biggl[ 1 + \frac{(g^2-\mathcal{B}) \xi}{\mathcal{A}} - \mathcal{D} \xi^3\biggr] \frac{d^2x}{d\xi^2} 
+ \biggl\{ 3  + \frac{4(g^2-\mathcal{B}) \xi}{\mathcal{A}} - 6\mathcal{D} \xi^3 \biggr\} 
\frac{1}{\xi} \cdot \frac{dx}{d\xi} 
+ \biggl[ 
\mathcal{D}   \mathfrak{F}_\mathrm{env} \xi^3 -\alpha_e  
\biggr]\frac{x}{\xi^2} \, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathcal{A}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~2\biggl(\frac{\rho_e}{\rho_c}\biggr)  \biggl(1 - \frac{\rho_e}{\rho_c} \biggr)  \, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{B}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~1  + 2\biggl(\frac{\rho_e}{\rho_c}\biggr)  - 3\biggl(\frac{\rho_e}{\rho_c}\biggr)^2  
\, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{D}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{1}{\mathcal{A}}\biggl( \frac{\rho_e}{\rho_c}\biggr)^2  = \biggl[ \frac{\rho_e/\rho_c}{2(1-\rho_e/\rho_c)} \biggr]
\, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathfrak{F}_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{3\omega^2_\mathrm{env}}{2\pi G \gamma_e \rho_e}   - 2\alpha_e
\, .
</math>
  </td>
</tr>
</table>
</div>

<span id="KeyConstraint">We have been unable</span> to demonstrate that this governing equation can be solved analytically for ''arbitrary'' pairs of the key model parameters, <math>~q</math> and <math>~\rho_e/\rho_c</math>.  But, if we choose parameter value pairs that satisfy the constraint,
<div align="center">
<math>~g^2 = \mathcal{B} </math>
&nbsp; &nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math> &nbsp; &nbsp; &nbsp; &nbsp; 
<math>~g = \frac{1}{1+2q^3} \, ,</math>
&nbsp; &nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp; &nbsp; 
<math>~q^3 = \mathcal{D} =  \biggl[ \frac{\rho_e/\rho_c}{2(1-\rho_e/\rho_c)} \biggr] </math>
&nbsp; &nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math> &nbsp; &nbsp; &nbsp; &nbsp; 
<math>~\frac{\rho_e}{\rho_c} = \frac{2q^3}{1+2q^3} \, ,</math>
</div>


<div align="center">
<table border="1" cellpadding="8" align="center" width="70%">
<tr><td align="left">
<font color="red">'''WRONG!'''</font>  &nbsp; The expression that relates <math>~g^2</math> to <math>~q^3</math> should read,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~g^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1+8q^3}{(1+2q^3)^2} = \frac{1+8\mathcal{D} }{(1+2\mathcal{D})^2}</math>
  </td>
</tr>
</table>
</div>
</td></tr>
</table>
</div>


then the LAWE that is relevant to the envelope simplifies.  Specifically, it takes the form,

<div align="center">
<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>~
( 1  - q^3 \xi^3 ) \frac{d^2x}{d\xi^2} + ( 3  - 6q^3 \xi^3 ) \frac{1}{\xi} \cdot \frac{dx}{d\xi} 
+ 
\biggl[ q^3  \mathfrak{F}_\mathrm{env} \xi^3 -\alpha_e  
\biggr]\frac{x}{\xi^2} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{x}{\xi^2}\biggl\{
( 1  - q^3 \xi^3 ) 
\biggl[ \frac{d}{d\ln\xi} \biggl( \frac{d\ln x}{d\ln \xi} \biggr) - \biggl(  1 -  \frac{d\ln x}{d\ln \xi} \biggr)\cdot \frac{d\ln x}{d\ln \xi}\biggr]
+ ( 3  - 6q^3 \xi^3 )  \frac{d\ln x}{d\ln \xi} 
+ \biggl[ q^3  \mathfrak{F}_\mathrm{env} \xi^3 -\alpha_e  
\biggr] \biggr\}
\, .
</math>
  </td>
</tr>
</table>
</div>
Shortly after deriving this last expression, we realized that one possible solution is a simple power-law eigenfunction of the form,
<div align="center">
<math>~x=a_0 \xi^{c_0} \, ,</math>
</div>
where the (constant) exponent is one of the roots of the quadratic equation,
<div align="center">
<math>~c_0^2 + 2c_0 - \alpha_e = 0 \, ,</math>
&nbsp; &nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math> &nbsp; &nbsp; &nbsp; &nbsp; 
<math>~c_0 = -1 \pm \sqrt{1+\alpha_e} \, .</math>
</div>
This power-law eigenfunction must be paired with the associated, dimensionless eigenfrequency parameter,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{F}_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~c_0(c_0+5) = 3c_0 + \alpha_e</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{3\omega^2_\mathrm{env}}{2\pi G \gamma_e \rho_e}   </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 3(c_0 + \alpha_e) = 3[\alpha_e -1 \pm \sqrt{1+\alpha_e}] \, .</math>
  </td>
</tr>
</table>
</div>

Next, [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Eureka_Regarding_Prasad.27s_1948_Paper|we noticed]] the strong similarities between the mathematical properties of this eigenvalue problem and the one that was studied by [http://adsabs.harvard.edu/abs/1948MNRAS.108..414P C. Prasad (1948, MNRAS, 108, 414-416)] in connection with, what we now recognize to be, a closely related problem.  Drawing heavily from Prasad's analysis, we discovered an infinite number of eigenfunctions (each, a truncated polynomial expression) and associated eigenfrequencies that satisfy this governing envelope LAWE.  The eigenvectors associated with the lowest few modes are tabulated, below.

==Eigenvector==
===Core Segment===

<div align="center">
<table border="1" cellpadding="8">
<tr>
  <td align="center" colspan="1">
Mode
  </td>
  <td align="center" colspan="1">
Core Eigenfunction<br />
<br />
<math>
g^2 \equiv 1  + \biggl(\frac{\rho_e}{\rho_c}\biggr)  \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_c} \biggr) \biggl( 1-q \biggr) + 
\frac{\rho_e}{\rho_c} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] 
</math>
  </td>
  <td align="center">Core Eigenfrequency <br />
<math>~\frac{3\omega_\mathrm{core}^2}{2\pi \gamma_c G \rho_c} = 2[\alpha_c + j(2j+5)]</math>
  </td>
</tr>
<tr>
<td align="center">
<math>~j=0 </math>
</td>
<td align="left">
<math>~x_\mathrm{core} = a_0 </math>
</td>
<td align="center">
<math>~6-8/\gamma_c</math>
</td>
</tr>

<tr>
<td align="center">
<math>~j=1 </math>
</td>
<td align="left">
<math>~x_\mathrm{core} = a_0 \biggl[ 1 - \frac{7}{5}\biggr(\frac{\xi^2}{g^2}\biggr) \biggr]</math>
</td>
<td align="center">
<math>~20-8/\gamma_c</math>
</td>
</tr>

<tr>
<td align="center">
<math>~j=2 </math>
</td>
<td align="left">
<math>~x_\mathrm{core} = a_0 \biggl[ 1 - \frac{18}{5}\biggr(\frac{\xi^2}{g^2}\biggr) + \frac{99}{35}\biggr(\frac{\xi^2}{g^2}\biggr)^2 \biggr]</math>
</td>
<td align="center">
<math>~42-8/\gamma_c</math>
</td>
</tr>
</table>
</div>

===Envelope Segment===

<div align="center">
<table border="1" cellpadding="8">
<tr>
  <td align="center" colspan="1">
Mode
  </td>
  <td align="center" colspan="1">
Envelope Eigenfunction<br />
<br />
<math>~c_0 \equiv -1 \pm \sqrt{1+\alpha_e}</math>
  </td>
  <td align="center">Envelope Eigenfrequency<br />
<math>~\frac{3\omega_\mathrm{env}^2}{2\pi \gamma_e G \rho_e} = 3[\alpha_e + c_0(2\ell+1) + \ell(3\ell+5)]</math><br /><br /><math>~=2\alpha_e + (c_0+3\ell)(c_0+3\ell + 5)</math>
  </td>
</tr>
<tr>
<td align="center">
<math>~\ell=0 </math>
</td>
<td align="left">
<math>~x_\mathrm{env} = b_0 \xi^{c_0}</math>
</td>
<td align="center">
<math>~3[\alpha_e + c_0]</math>
</td>
</tr>

<tr>
<td align="center">
<math>~\ell=1 </math>
</td>
<td align="left">
<math>~x_\mathrm{env} = b_0 \xi^{c_0}\biggl\{1 + \biggl[  
\frac{c_0(c_0+5)-(c_0+3)(c_0+8)}{(c_0+3)(c_0+5) - \alpha_e}
\biggr](q\xi)^3  \biggr\}</math>
</td>
<td align="center">
<math>~3[\alpha_e + 3c_0 +8]</math>
</td>
</tr>

<tr>
<td align="center">
<math>~\ell=2 </math>
</td>
<td align="left">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x_\mathrm{env} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~b_0 \xi^{c_0}\biggl\{1 + \biggl[ \frac{c_0(c_0+5)-(c_0+6)(c_0+11)}{(c_0+3)(c_0+5) - \alpha_e}\biggr](q\xi)^3 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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]\biggl[ \frac{(c_0+3)(c_0+8)-(c_0+6)(c_0+11)}{(c_0+6)(c_0+8) - \alpha_e}\biggr](q\xi)^6 
\biggr\}
</math>
  </td>
</tr>
</table>


</td>
<td align="center">
<math>~3[\alpha_e + 5c_0 +22]</math>
</td>
</tr>
</table>
</div>

===Piecing Together===
Here we illustrate how the two segments of the eigenfunction can be successfully pieced together for the specific case of <math>~(\ell,j) = (2,1)</math>.

<span id="STEP1"><font color="red"><b>STEP 1:</b></font></span>  Choose a value of the adiabatic exponent for the envelope, <math>~\gamma_e</math>.  Then, the values of both <math>~\alpha_e</math> and <math>~c_0</math> are known as well; actually, because it is the root of a quadratic equation, <math>~c_0</math> can, in general, take on one of a ''pair'' of values.  We will elaborate on this further, below.

<span id="STEP2"><font color="red"><b>STEP 2:</b></font></span>  Acknowledging that the value of <math>~q</math> has yet to be determined, fix the value of the leading, overall scaling coefficient, <math>~b_0</math>, such that<sup>&dagger;</sup> <math>~x_\mathrm{env} = 1</math> at the interface, that is, at <math>~\xi = 1</math>.  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>
  <td align="right">
<math>~x_{\ell=2} |_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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>
  </td>
</tr>
</table>
</div>
where, the values of the newly introduced coefficients,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
  <td align="right">
<math>~A_{21}</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">
<math>~B_{21}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<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>
  </td>
</tr>
</table>
</div>
are also both known.

<table border="1" cellpadding="8" width="70%" align="center">
<tr><td align="center"><font color="red"><b>Wrong</b></font> expression for <math>~g^2</math></td></tr>
<tr><td align="left">
<font color="red"><b>STEP 3:</b></font>  Recognizing that this segment of the eigenfunction will only satisfy the envelope's LAWE if we restrict our discussion to equilibrium models for which <math>~g^2 = \mathcal{B} = (1+2q^3)^{-1}</math>, we must insert this same restriction on <math>~g^2</math> into the core's eigenfunction.  At the same time, we should fix the value of the leading, overall scaling coefficient, <math>~a_0</math>, such that<sup>&dagger;</sup> <math>~x_\mathrm{core} = 1</math> at the interface <math>~(\xi = 1)</math>.  For the case of <math>~j=1</math>, this means that, throughout the core, the eigenfunction is,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x_{j=1} |_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{5 -  7 (1+2q^3)^2 \xi^2}{5-7(1+2q^3)^2} \, .</math>
  </td>
</tr>
</table>
</div>
</td></tr></table>

<span id="STEP3"><font color="red"><b>STEP 3:</b></font></span>  Recognizing that this segment of the eigenfunction will only satisfy the envelope's LAWE if we restrict our discussion to equilibrium models for which <math>~g^2 = \mathcal{B} = [(1+8q^3)/(1+2q^3)^{2}]</math>, we must insert this same restriction on <math>~g^2</math> into the core's eigenfunction.  At the same time, we should fix the value of the leading, overall scaling coefficient, <math>~a_0</math>, such that<sup>&dagger;</sup> <math>~x_\mathrm{core} = 1</math> at the interface <math>~(\xi = 1)</math>.  For the case of <math>~j=1</math>, this means that, throughout the core, the eigenfunction is,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x_{j=1} |_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{5(1+8q^3) -  7 (1+2q^3)^2 \xi^2}{5(1+8q^3)-7(1+2q^3)^2} \, .</math>
  </td>
</tr>
</table>
</div>


<span id="STEP4"><font color="red"><b>STEP 4:</b></font></span>  Now we need to ensure that the first derivative with respect to the radial coordinate, <math>~\xi</math>, of both eigenfunction segments match at the interface <math>~(\xi = 1)</math>.  Making the notation substitution, <math>~\Chi \equiv q^3</math>, the mathematical relation that accomplishes this is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{14(1+2\Chi)^2}{7(1+2\Chi)^2 - 5(1+8\Chi)}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~
0
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2] [7(1+2\Chi)^2 - 5(1+8\Chi)]
- 14(1+2\Chi)^2 [1 + A_{21}\Chi + A_{21}B_{21}\Chi^2] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2] [2 - 12\Chi + 28\Chi^2 ]
- (14 + 56\Chi + 56 \Chi^2) [1 + A_{21}\Chi + A_{21}B_{21}\Chi^2] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2[c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2] 
-12\Chi [c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2] 
+ 28\Chi^2 [c_0 + (c_0 + 3)A_{21}\Chi + (c_0 + 6)A_{21}B_{21} \Chi^2] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
- 14 [1 + A_{21}\Chi + A_{21}B_{21}\Chi^2] 
- 56\Chi  [1 + A_{21}\Chi + A_{21}B_{21}\Chi^2] 
- 56 \Chi^2 [1 + A_{21}\Chi + A_{21}B_{21}\Chi^2] \, .
</math>
  </td>
</tr>
</table>
</div>

Or we have, equivalently,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~a\Chi^4 + b\Chi^3 + c\Chi^2 +d\Chi +e </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, ,</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~e</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~ 2c_0 - 14  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~d</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~2(c_0+3)A_{21} - 12c_0 - 14A_{21} - 56</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2(c_0-4)A_{21} - 4(14 + 3c_0) \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~c</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~2(c_0+6)A_{21}B_{21} - 12(c_0+3)A_{21} + 28c_0 - 14A_{21}B_{21} - 56A_{21} - 56  </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2(c_0 - 1)A_{21}B_{21} - 4(23 + 3c_0)A_{21} + 28(c_0 - 2)  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~b</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~-12(c_0+6)A_{21}B_{21} + 28(c_0+3)A_{21} - 56A_{21}B_{21} - 56A_{21} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- 4(3c_0+32)A_{21}B_{21} + 28(c_0+1)A_{21} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~a</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~28(c_0+4)A_{21}B_{21} \, .</math>
  </td>
</tr>
</table>

<table border="1" cellpadding="8" width="70%" align="center">
<tr><td align="center"><font color="red"><b>Wrong</b></font> expressions for quadratic equation coefficients, based on previous, wrong expression for <math>~g^2</math></td></tr>
<tr><td align="left">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~e</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~ 14- 2c_0  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~d</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~28(2-c_0) + (8-2c_0) A_{21} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~c</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~28(2-c_0) - 28(1+c_0) A_{21} + 2(1 - c_0)A_{21}B_{21}  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~b</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~-28A_{21}[(1+c_0) + (4+c_0)B_{21} ]  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~a</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~- 28 (4+c_0)A_{21}\cdot B_{21}  \, .</math>
  </td>
</tr>
</table>
</td></tr></table>

The physically relevant (real) root of this quartic equation in <math>~\Chi</math> &#8212; see our [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Quartic|accompanying detailed presentation]] &#8212; gives us the ''specific'' value of the dimensionless interface location, <math>~q</math>, for which both the values and the first derivatives of the two eigenfunction segments match at the interface.  As is shown by the plot displayed in the right-hand panel of Figure 1, we have found different values of <math>~q</math> for each choice (STEP 1) of <math>~\gamma_e</math> (or, equivalently, choice of <math>~\alpha_e</math>).  In this plot we have purposely flipped the horizontal axis so that the extreme left <math>~(\alpha_e = +3)</math> represents an incompressible <math>~(n = 0)</math> envelope, while the extreme right represents an isothermal <math>~(\gamma_e = 1)</math> envelope. 
<!--
The green and orange curves, respectively, show as well how the corresponding model parameters, <math>~\nu</math> and <math>~\rho_e/\rho_c</math>, vary with <math>~\alpha_e</math>.  The right-hand panel displays one example of an <math>~(\ell, j) = (2,1)</math> eigenfunction that simultaneously satisfies the LAWE of the core and the LAWE of the envelope, and matches smoothly at the interface.  This ''particular'' plotted solution corresponds to the case of <math>~\alpha_e = -0.35</math>, for which:  <math>~\gamma_e =1.1940299</math>; <math>~n_e = 5.1538462</math>; <math>~q = 0.7943853</math>; <math>~\nu =0.6675302</math>; and <math>~\rho_e/\rho_c =0.5006468</math>.
-->

<div align="center" id="Figure1">
<table border="1" cellpadding="8">
<tr><th align="center" colspan="4">Figure 1</th></tr>
<tr>
  <td align="center" colspan="3"><math>~\alpha_e = -0.35 \, ;~~~c_0 = \sqrt{1+\alpha_e} - 1</math></td>
  <td align="center" rowspan="18">

<table border="0" width="100%">
<tr><td align="center">
[[File:Quartic21Solution02Corrected.png|500px|quartic solution]]
</td></tr>
<tr><td align="left">''Directly Above:'' Plot shows for ''which'' equilibrium bipolytropic configurations with <math>~(n_c, n_e) = (0,0)</math> we are able to construct analytically prescribed eigenvectors for the radial oscillation mode, <math>~(\ell, j) = (2,1)</math>.  The top (blue), middle (green), and bottom (orange) curves show how <math>~q</math>, <math>~\nu</math>, and <math>~\rho_e/\rho_c</math> vary with the specified value of the envelope's adiabatic exponent over the full, physically reasonable range of the parameter, <math>~-1 \le \alpha_e \le 3</math>.  For the upper portion of each curve (dark blue, dark green, dark orange), the parameter, <math>~c_0</math>, is taken to be the "plus" root of its defining quadratic equation; the "minus" root defines <math>~c_0</math> along the lower portion of each curve (light blue, light green, light orange).  

''Upper-left Quadrant:'' An <math>~x(r_0/R)</math> plot showing the radial structure of the analytically prescribed eigenfunction for <math>~\alpha_e = -0.35</math> and <math>~c_0</math> (plus); its underlying, equilibrium model characteristics are identified by the black circular marker in the above plot.

''Lower-left Quadrant:'' The analytcially prescribed eigenfunction, <math>~x(r_0/R)</math>,  for <math>~\alpha_e = -0.9</math> and <math>~c_0</math> (minus); its underlying, equilibrium model characteristics are identified by the yellow circular marker in the above plot.

<sup>&dagger;</sup>Note that, as displayed here, the sign has been flipped on both <math>~x(r_0/R)</math> eigenfunctions so that, in practice, the amplitude at the interface is ''negative'' one, rather than positive one. Plotted in this way, we immediately recognize that both eigenfunctions are ''qualitatively'' similar to the <math>~j = 2</math> radial oscillation eigenfunction that [[User:Tohline/SSC/UniformDensity#Properties_of_Eigenfunction_Solutions|was derived by Sterne (1937) in the context of isolated, homogeneous spheres]].
</td></tr>
</table>

  </td>
</tr>
<tr>
  <td align="right"><math>~c_0</math> (plus):</td>
  <td align="center"><math>~-0.1937742</math></td>
  <td align="center" rowspan="8">
[[File:EigenfunctionP1Corrected.png|270px|quartic solution]]
  </td>
</tr>
<tr>
  <td align="right"><math>~\gamma_e</math>:</td>
  <td align="center"><math>~1.1940299</math></td>
</tr>
<tr>
  <td align="right"><math>~n_e</math>:</td>
  <td align="center"><math>~5.1538462</math></td>
</tr>
<tr>
  <td align="right"><math>~q</math>:</td>
  <td align="center"><math>~0.6840119</math></td>
</tr>
<tr>
  <td align="right"><math>~\nu</math>:</td>
  <td align="center"><math>~0.5466868</math></td>
</tr>
<tr>
  <td align="right"><math>~\rho_e/\rho_c</math>:</td>
  <td align="center"><math>~0.3902664</math></td>
</tr>
<tr>
  <td align="right"><math>~\alpha_c</math>:</td>
  <td align="center"><math>~+0.8326585</math></td>
</tr>
<tr>
  <td align="right"><math>~\gamma_c</math>:</td>
  <td align="center"><math>~+1.845579</math></td>
</tr>
<tr>
  <td align="center" colspan="3"><math>~\alpha_e = -0.9 \, ;~~~c_0 = -\sqrt{1+\alpha_e} - 1</math></td>
</tr>
<tr>
  <td align="right"><math>~c_0</math> (minus):</td>
  <td align="center"><math>~- 1.3162278</math></td>
  <td align="center" rowspan="8">
[[File:EigenfunctionM1Corrected.png|270px|quartic solution]]
  </td>
</tr>
<tr>
  <td align="right"><math>~\gamma_e</math>:</td>
  <td align="center"><math>~1.0256410</math></td>
</tr>
<tr>
  <td align="right"><math>~n_e</math>:</td>
  <td align="center"><math>~39</math></td>
</tr>
<tr>
  <td align="right"><math>~q</math>:</td>
  <td align="center"><math>~0.5728050</math></td>
</tr>
<tr>
  <td align="right"><math>~\nu</math>:</td>
  <td align="center"><math>~0.4586270</math></td>
</tr>
<tr>
  <td align="right"><math>~\rho_e/\rho_c</math>:</td>
  <td align="center"><math>~0.2731929</math></td>
</tr>
<tr>
  <td align="right"><math>~\alpha_c</math>:</td>
  <td align="center"><math>~-0.9595214</math></td>
</tr>
<tr>
  <td align="right"><math>~\gamma_c</math>:</td>
  <td align="center"><math>~+1.0102231</math></td>
</tr>
</table>
</div>


<span id="STEP5"><font color="red"><b>STEP 5:</b></font></span> Finally, for each choice of <math>~\gamma_e</math> &#8212; or, alternatively, <math>~\alpha_e</math> &#8212; the physically relevant value of the core's adiabatic exponent is set by demanding that the ''dimensional'' eigenfrequencies of the envelope and core precisely match.  That is, we demand that,
<table border="1" cellpadding="8" align="right">
<tr><th align="center" colspan="1">Figure 2</th></tr>
<tr>
  <td align="center">
[[File:AlphaVsAlpha21BothCorrected.png|350px|quartic solution]]
  </td>
</tr>
</table>
<div align="center">
<math>~\omega^2_\mathrm{env} = \omega^2_\mathrm{core} \, .</math>
</div>
From [[#Eigenvector|above]], we know that, for the core,
<div align="center">
<math>~3\omega^2_\mathrm{core}\biggr|_\mathrm{j=1} =  2\pi \gamma_c G \rho_c [ 20 - 8/\gamma_c] \, ;</math>
</div>
whereas, for the envelope,
<div align="center">
<math>~3\omega^2_\mathrm{env}\biggr|_\mathrm{\ell=2} =  2\pi \gamma_e G \rho_e [ 3(\alpha_e + 5c_0 + 22)] \, .</math>
</div>
By demanding that these frequencies be identical, we conclude that,
<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}{20} \biggl[ 8 + 3\gamma_e \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl(\alpha_e + 5c_0 + 22 \biggr)\biggr] \, .</math>
  </td>
</tr>
<!--
<tr>
  <td align="right">
<math>~ \Rightarrow ~~~ \gamma_c </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{20}\biggl\{ 8 + \biggl(\frac{12}{3-\alpha_e}\biggr) \biggl[\alpha_e + 5(\sqrt{1+\alpha_e}-1) + 22 \biggr]\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggr\}\, .</math>
  </td>
</tr>
-->
</table>
</div>

Figure 2 shows how the required value of <math>~\alpha_c</math> varies with the choice of <math>~\alpha_e</math>; here, both axes have been flipped in order to run from incompressible <math>~(\alpha = +3)</math> at the left/bottom, to isothermal <math>~(\alpha = -1)</math> at the right/top.  For the lower portion of the curve (red circular markers), the parameter, <math>~c_0</math>, is taken to be the "plus" root of its defining quadratic equation; the "minus" root defines <math>~c_0</math> along the upper portion of the curve (purple circular markers).  The diagonal dashed-black line identifies where <math>~\alpha_c = \alpha_e</math>; in models below and to the right of this line, the envelope is more compressible than is the core, whereas in models above and to the left of this line, the core is more compressible than the envelope.  

<!-- OMIT Evidently there is one model for which the <math>~(\ell,j) = (2,1)</math> eigenvector is analytically specifiable in which the envelope and core are equally compressible; it is the model with <math>~\gamma_c = \gamma_e \approx 1.13</math> that is identified by where the <math>~c_0</math> (minus) segment of the curve intersects the diagonal black-dashed line.
-->

The eigenfrequency that corresponds to the ''specific'' eigenfunction that is displayed in upper-left quadrant of Figure 1 is identified by the black circular marker in Figure 2; as is indicated by the row of numbers on the left in Figure 1, this model has,
<div align="center">
<math>~\gamma_c = 1.845579 </math>
&nbsp; &nbsp; &nbsp; <math>~\Rightarrow </math>&nbsp; &nbsp; &nbsp; 
<math>~\alpha_c = +0.8326535 \, . </math>
</div>
The yellow circular marker in Figure 2 identifies the model whose analytically prescribed, <math>~(\ell,j) = (2,1)</math> eigenfunction is displayed in the lower-left quadrant of Figure 1; it has,
<div align="center">
<math>~\gamma_c = 1.0102231 </math>
&nbsp; &nbsp; &nbsp; <math>~\Rightarrow </math>&nbsp; &nbsp; &nbsp; 
<math>~\alpha_c = -0.9595214 \, . </math>
</div>

==Examining Alignment with Surface Boundary Condition==

===Expectation===
As we have reviewed in [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Boundary_Condition|an accompanying discussion]], one astrophysically reasonable surface boundary condition provides a mathematical relationship between the logarithmic derivative of the eigenfunction with respect to the radius, in terms of the eigenfrequency as follows:


<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~ \frac{d\ln x}{d\ln \xi}\biggr|_{\xi = 1/q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{3\omega^2 }{4\pi G\rho_c \gamma_e} \biggl( \frac{\nu}{q^3}\biggr) - \biggl( 3 - \frac{4}{\gamma_e}\biggr) </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{3\omega^2 }{4\pi G\rho_c \gamma_e} \biggl( \frac{1+2q^3}{3q^3}\biggr) - \alpha_e </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{3} \biggl[\frac{3\omega^2 }{2\pi G\rho_c \gamma_e} \biggr] \biggl( \frac{\rho_c}{\rho_e}\biggr) - \alpha_e </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{3} \biggl[\mathfrak{F}_\mathrm{env} + 2\alpha_e \biggr]  - \alpha_e </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{3} \biggl[\mathfrak{F}_\mathrm{env} - \alpha_e \biggr]  \, . </math>
  </td>
</tr>
</table>
</div>

Now, according to our [[#Envelope_Segment|above-described envelope segment of the eigenfunction]], we established the analytic prescription,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{F}_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(c_0 + 3\ell)(c_0 + 3\ell+5) \, ,</math>
  </td>
</tr>
</table>
</div>
in which case the desired surface boundary condition is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~ 3 \cdot \frac{d\ln x}{d\ln \xi}\biggr|_{\xi = 1/q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ (c_0 + 3\ell)(c_0 + 3\ell+5) - \alpha_e  </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ [c_0^2 + c_0(6\ell + 5 ) + 3\ell(3\ell+5)] - (c_0^2 + 2c_0)</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 3[ c_0(2\ell + 1 ) + \ell(3\ell+5)]</math>
  </td>
</tr>
</table>
</div>

That is, we expect to find the following,
<div align="center" id="DesiredBoundaryCondition">
<table border="1" align="center"><tr><td align="center">

<table border="0" cellpadding="5" align="center">
<tr>
  <th align="center" colspan="3">
Desired Boundary Condition
  </th>
</tr>

<tr>
  <td align="right">
<math>~\frac{d\ln x}{d\ln \xi}\biggr|_{\xi = 1/q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ c_0(2\ell + 1 ) + \ell(3\ell+5) \, .</math>
  </td>
</tr>
</table>

</td></tr></table>
</div>

===Analytic2===
Continuing, from above, a discussion specifically of the case, <math>~\ell = 2</math>, the analytically specified envelope eigenfunction is,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x_{\ell=2} |_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi^{c_0}\biggl[ \frac{ 1 +  q^3 A_{\ell=2} \xi^{3} +  q^6 A_{\ell=2}B_{\ell=2}\xi^{6} }{ 1 +  q^3 A_{\ell=2}  +  q^6 A_{\ell=2}B_{\ell=2}}\biggr] \, ,
</math>
  </td>
</tr>
</table>
</div>
where, the values of the newly introduced coefficients,
<div align="center">
<table border="0" cellpadding="3" align="center">
<tr>
  <td align="right">
<math>~A_{\ell=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] = \frac{-2(2c_0+11)}{(2c_0+5)} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~B_{\ell=2}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{(c_0+3)(c_0+8) - (c_0 + 6)(c_0 + 11)}{(c_0 + 6)(c_0+8) - \alpha_e}\biggr] = \frac{-(c_0+7)}{2(c_0+4)} \, ,</math>
  </td>
</tr>
</table>
</div>
in which case,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{d\ln x}{d\ln \xi} = \frac{\xi}{x} \cdot \frac{dx}{d\xi} \biggl|_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\xi}{x} \biggl\{
c_0\xi^{c_0-1}\biggl[ \frac{ 1 +  q^3 A \xi^{3} +  q^6 AB\xi^{6} }{ 1 +  q^3 A  +  q^6 AB}\biggr] 
+ \xi^{c_0}\biggl[ \frac{ 3q^3 A \xi^{2} +  6q^6 AB\xi^{5} }{ 1 +  q^3 A  +  q^6 AB}\biggr] 
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0
+ \biggl[ \frac{ 3q^3 A \xi^{3} +  6q^6 AB\xi^{6} }{ 1 +  q^3 A\xi^3  +  q^6 AB\xi^6}\biggr] \, .
</math>
  </td>
</tr>
</table>
</div>
Hence, at the surface <math>~(\xi = 1/q)</math>, we find,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{d\ln x}{d\ln \xi} \biggl|_{\xi=1/q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0
+\biggl[ \frac{ 3 A  +  6 AB }{ 1 +  A  +  AB}\biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0
+\biggl[ \frac{ -12(2c_0+11)(c_0+4)  + 12(2c_0+11)(c_0+7) }{ 2(2c_0 + 5)(c_0+4) -  4(2c_0+11)(c_0+4)  +  2(2c_0+11) (c_0+7)}\biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0
+6\biggl[ \frac{ (2c_0^2 + 25c_0 + 77) -(2c_0^2 + 19c_0 +44)  }{ (2c_0^2 + 13c_0 + 20) -  2(2c_0^2 + 19c_0 + 44)  +  (2c_0^2 + 25c_0 + 77)}\biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
c_0
+6 \biggl[ \frac{ 6c_0 +33 }{ 9}\biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
5c_0 + 22 \, .
</math>
  </td>
</tr>
</table>
</div>

It is gratifying &#8212; although, somewhat surprising (to me!) &#8212; to find that this precisely matches the [[#DesiredBoundaryCondition|above-defined, desired boundary condition]] for the case of <math>~\ell = 2</math>.

===Duh!===

[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 4 February 2017:  This numerical determination of surface boundary conditions was carried out inside spreadsheet "FDflex22" of Excel file ''analyticeigenvectorcorrected.xlsx''.]]After also checking conformance with the expected boundary condition in [[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Check_Surface_Boundary_Condition|the case of analytic eigenfunctions having <math>~\ell = 3</math>]] and, separately (not shown), for numerically generated eigenfunctions having a wide range of oscillation frequencies, it dawned on us that the [[#DesiredBoundaryCondition|"desired" surface boundary condition]] may actually be a natural outcome of the envelope's LAWE.  

By constraining our discussion to models for which <math>~g^2 = \mathcal{B}</math> and <math>~\mathcal{D} = q^3</math>, the [[#The_Envelope.27s_LAWE|envelope's LAWE]] is,

<div align="center">
<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>~
\biggl[ 1  - q^3 \xi^3\biggr] \frac{d^2x}{d\xi^2} 
+ \biggl\{ 3   - 6q^3 \xi^3 \biggr\} 
\frac{1}{\xi} \cdot \frac{dx}{d\xi} 
+ \biggl[ 
q^3  \mathfrak{F}_\mathrm{env} \xi^3 -\alpha_e  
\biggr]\frac{x}{\xi^2} \, .
</math>
  </td>
</tr>
</table>
</div>
At the surface <math>~(\xi = 1/q)</math>, the coefficient of the second derivative term goes to zero, in which case the LAWE reduces in form to,

<div align="center">
<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>~
-\frac{3}{\xi} \cdot \frac{dx}{d\xi} 
+ \biggl[ 
\mathfrak{F}_\mathrm{env}  -\alpha_e  
\biggr]\frac{x}{\xi^2} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 3\cdot \frac{d\ln x}{d\ln \xi}\biggr|_\mathrm{surface} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\mathfrak{F}_\mathrm{env}  -\alpha_e  \, .
</math>
  </td>
</tr>
</table>
</div>

And this is ''precisely'' the condition that derives from the astrophysically reasonable boundary condition that we have [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Boundary_Condition|discussed separately]] and that has been [[#Expectation|reviewed, above]].

===Broader Analysis===

Let's, then, examine the behavior of the envelope's LAWE at the surface in the most general case &#8212; that is, when ''not'' constrained to <math>~g^2 = \mathcal{B}</math>.  First, we note that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~g^2 - \mathcal{B}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl(\frac{\rho_e}{\rho_c}\biggr)  \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_c} \biggr) \biggl( 1-q \biggr) + 
\frac{\rho_e}{\rho_c} \biggl(\frac{1}{q^2} - 1\biggr) \biggr]  -2\biggl(\frac{\rho_e}{\rho_c}\biggr)  + 3\biggl(\frac{\rho_e}{\rho_c}\biggr)^2
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl\{ 2 \biggl[1 - \biggl(\frac{\rho_e}{\rho_c} \biggr) -q  + q \biggl(\frac{\rho_e}{\rho_c} \biggr)\biggr] + 
\biggl( \frac{\rho_e}{\rho_c} \biggr)\frac{1}{q^2}  -2  + 2\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[ - 2q  + 2q \biggl(\frac{\rho_e}{\rho_c} \biggr) + 
\biggl( \frac{\rho_e}{\rho_c} \biggr)\frac{1}{q^2}    \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\frac{1}{q^2} \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) + 2q^3 \biggl(\frac{\rho_e}{\rho_c}  - 1  \biggr) \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

Hence, at the surface quite generally, the coefficient of the second derivative is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{1}{\mathcal{A}}\biggl[\mathcal{A} + (g^2 - \mathcal{B})\xi - \mathcal{A}\mathcal{D} \xi^3 \biggr]_{\xi=1/q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{\mathcal{A}}\biggl\{
2\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl(1-\frac{\rho_e}{\rho_c}\biggr)
+ \frac{1}{q^3} \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) + 2q^3 \biggl(\frac{\rho_e}{\rho_c}  - 1  \biggr) \biggr] 
- \frac{1}{q^3}\biggl(\frac{\rho_e}{\rho_c}\biggr)^2
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{\mathcal{A}} \biggl(\frac{\rho_e}{\rho_c}\biggr)\biggl\{
- 2\biggl(\frac{\rho_e}{\rho_c} -1\biggr)
+ \frac{1}{q^3}  \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) + 2q^3 \biggl(\frac{\rho_e}{\rho_c}  - 1  \biggr) \biggr] 
- \frac{1}{q^3}\biggl(\frac{\rho_e}{\rho_c}\biggr)
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
0 \, .
</math>
  </td>
</tr>
</table>
</div>

And, at the surface quite generally, the coefficient of the first derivative is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{1}{\mathcal{A}}\biggl[3\mathcal{A} + 4(g^2 - \mathcal{B})\xi - 6\mathcal{A}\mathcal{D} \xi^3 \biggr]_{\xi=1/q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{\mathcal{A}}\biggl\{
6\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl(1-\frac{\rho_e}{\rho_c}\biggr)
+ \frac{4}{q^3} \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) + 2q^3 \biggl(\frac{\rho_e}{\rho_c}  - 1  \biggr) \biggr] 
- \frac{6}{q^3}\biggl(\frac{\rho_e}{\rho_c}\biggr)^2
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2}{\mathcal{A}} \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl\{
-3\biggl(\frac{\rho_e}{\rho_c} - 1\biggr)
+ \frac{2}{q^3} \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) + 2q^3 \biggl(\frac{\rho_e}{\rho_c}  - 1  \biggr) \biggr] 
- \frac{3}{q^3}\biggl(\frac{\rho_e}{\rho_c}\biggr)
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2}{\mathcal{A}} \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[
\biggl(\frac{\rho_e}{\rho_c} - 1\biggr)
- \frac{1}{q^3}\biggl(\frac{\rho_e}{\rho_c}\biggr)
\biggr] \, .
</math>
  </td>
</tr>
</table>
</div>
Hence, at the surface quite generally, the envelope's LAWE becomes,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~
- \biggl[ 3\mathcal{A}  + 4(g^2-\mathcal{B}) \xi - 6\mathcal{A} \mathcal{D} \xi^3 \biggr]_{\xi=1/q} \frac{d\ln x}{d\ln\xi} \biggr|_\mathrm{surface}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\mathcal{A}\biggl[ 
\mathcal{D}   \mathfrak{F}_\mathrm{env} \xi^3 -\alpha_e  \biggr]_{\xi=1/q} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~
- 2\biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[
\biggl(\frac{\rho_e}{\rho_c} - 1\biggr)
- \frac{1}{q^3}\biggl(\frac{\rho_e}{\rho_c}\biggr)
\biggr]  \frac{d\ln x}{d\ln\xi} \biggr|_\mathrm{surface}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\rho_e}{\rho_c}\biggr)^2 \mathfrak{F}_\mathrm{env} \cdot \frac{1}{q^3} - 2\alpha_e \biggl(\frac{\rho_e}{\rho_c}\biggr)\biggl(1-  \frac{\rho_e}{\rho_c}\biggr)
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~
\biggl[2\biggl(\frac{\rho_e}{\rho_c}\biggr) +
2q^3\biggl(1 - \frac{\rho_e}{\rho_c} \biggr)
\biggr]  \frac{d\ln x}{d\ln\xi} \biggr|_\mathrm{surface}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\rho_e}{\rho_c}\biggr) \mathfrak{F}_\mathrm{env}  - 2q^3 \alpha_e \biggl(1-  \frac{\rho_e}{\rho_c}\biggr)
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~
\frac{d\ln x}{d\ln\xi} \biggr|_\mathrm{surface}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[2 + 2q^3\biggl(1 - \frac{\rho_e}{\rho_c} \biggr)\biggl(\frac{\rho_e}{\rho_c}\biggr)^{-1} \biggr]  ^{-1}
\biggl[ \mathfrak{F}_\mathrm{env}  - 2q^3 \alpha_e \biggl(1-  \frac{\rho_e}{\rho_c}\biggr)\biggl(\frac{\rho_e}{\rho_c}\biggr)^{-1} \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\mathfrak{F}_\mathrm{env} - \Kappa \alpha_e}{2+\Kappa} \, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<math>~\Kappa \equiv 2q^3\biggl(1 - \frac{\rho_e}{\rho_c} \biggr)\biggl(\frac{\rho_e}{\rho_c}\biggr)^{-1} \, .</math>
</div>

Notice that in the special case for which we have been able to identify analytically specifiable eigenvectors, namely, when
<div align="center">
<math>~g^2 = \mathcal{B} </math>&nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math>&nbsp; &nbsp; &nbsp; <math>~\Kappa = 1 \, ,</math>
</div>

this surface boundary condition simplifies to the ''expected'' expression,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~
\frac{d\ln x}{d\ln\xi} \biggr|_\mathrm{surface}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{3}
\biggl[ \mathfrak{F}_\mathrm{env}  - \alpha_e \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

Under what condition &#8212; other than when <math>~g^2=\mathcal{B}</math> &#8212; does the general expression generate the ''expected'' expression?  We need,

<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{1}{3}
\biggl[ \mathfrak{F}_\mathrm{env}  - \alpha_e \biggr] 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[2 + \Kappa \biggr]  ^{-1}
\biggl[ \mathfrak{F}_\mathrm{env}  - \alpha_e \Kappa \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ (2 + \Kappa )
\biggl[ \mathfrak{F}_\mathrm{env}  - \alpha_e \biggr] 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
3\biggl[ \mathfrak{F}_\mathrm{env}  - \alpha_e \Kappa \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
(2 + \Kappa ) \mathfrak{F}_\mathrm{env}  - (2 + \Kappa )  \alpha_e 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
3\mathfrak{F}_\mathrm{env}  - 3\Kappa \alpha_e 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
(\Kappa -1) \mathfrak{F}_\mathrm{env}  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-2(\Kappa-1 )  \alpha_e 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
\mathfrak{F}_\mathrm{env}  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-2\alpha_e \, .
</math>
  </td>
</tr>
</table>
</div>

But, given that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{F}_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{3\omega^2_\mathrm{env}}{2\pi G \gamma_e \rho_e}   - 2\alpha_e
\, ,
</math>
  </td>
</tr>
</table>
</div>

we see that the ''expected'' boundary condition will result only for <math>~\omega_\mathrm{env}^2 = 0</math>, that is, only for, <math>~\sigma_c^2 = 0</math>.  This is what we have been noticing as we have played with numerically generated eigenvectors:  When integrating from the center of the zero-zero bipolytrope, to its surface, the naturally resulting (first) derivative of the eigenfunction at the surface of the configuration matches the ''expected'' surface boundary condition &hellip;
* for all values of <math>\sigma_c^2</math>, when <math>~g^2= \mathcal{B}</math>, that is, when <math>~\Kappa=1</math>;
* only for <math>~\sigma_c^2 = 0</math> in all other configurations, that is, for all <math>~\Kappa \ne 1</math>.

What do we make of this?

==Five Mode Summary==

<div align="center">
<table border="1" cellpadding="5">
<tr>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (2,1)</math></td>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (2,2)</math></td>
  <td align="center" bgcolor="white">&nbsp;</td>
</tr>
<tr>
  <td align="center">[[File:Annotate21.png|300px|Log(amplitude) plot for (ell,j) = (2,1)]]</td>
  <td align="center">[[File:Annotate22.png|300px|Log(amplitude) plot for (ell,j) = (2,2)]]</td>
  <td align="center">&nbsp;</td>
</tr>
<tr>
  <td align="right">
<div align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 37.08874</math></div>
<font size="-2">[[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Illustration21|more details &hellip;]]</font>
  </td>
  <td align="right">
<div align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.95210</math></div>
<font size="-2">[[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Illustration22|more details &hellip;]]</font>
  </td>
  <td align="center">&nbsp;</td>
</tr>
<tr>
  <td align="center" bgcolor="lightblue" colspan="3">&nbsp;</td>
</tr>
<tr>
  <td align="center" bgcolor="white"><math>~(\ell, j) = (3,1)</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>
  <td align="center">[[File:Annotate31.png|300px|Log(amplitude) plot for (ell,j) = (3,1)]]</td>
  <td align="center">[[File:Annotate32.png|300px|Log(amplitude) plot for (ell,j) = (3,2)]]</td>
  <td align="center">[[File:Annotate33.png|300px|Log(amplitude) plot for (ell,j) = (3,3)]]</td>
</tr>
<tr>
  <td align="right">
<div align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 12.452545</math></div>
<font size="-2">[[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Illustration31|more details &hellip;]]</font>
  </td>
  <td align="right">
<div align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.05461</math></div>
<font size="-2">[[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Illustration32|more details &hellip;]]</font>
  </td>
  <td align="right">
<div align="center"><math>~\frac{3\omega^2}{2\pi G\rho_c} =  87.41594</math></div>
<font size="-2">[[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Illustration33|more details &hellip;]]</font>
  </td>
</tr>
</table>
</div>



<div align="center">
<table border="1" align="center" cellpadding="5">
<tr>
  <th align="center">Model</th>
  <th align="center"><math>~\ell</math></th>
  <th align="center"><math>~j</math></th>
  <th align="center"><math>~q</math></th>
  <th align="center"><math>~\gamma_e</math></th>
  <th align="center"><math>~\gamma_c</math></th>
  <th align="center"><math>~\nu</math></th>
  <th align="center"><math>~\frac{\rho_e}{\rho_c}</math></th>
  <th align="center"><math>~\alpha_e</math></th>
  <th align="center"><math>~\alpha_c</math></th>
  <th align="center"><math>~n_e</math></th>
  <th align="center"><math>~n_c</math></th>
  <th align="center"><math>~g^2</math></th>
  <th align="center"><math>~f</math></th>
  <th align="center"><math>~\sigma_\mathfrak{G}^2</math></th>
  <th align="center">Analytic<br /><math>~\sigma_c^2</math></th>
</tr>
<tr>
  <td align="center">Analytic21</td>
  <td align="center">2</td>
  <td align="center">1</td>
  <td align="center">0.684</td>
  <td align="center">1.194</td>
  <td align="center">1.846</td>
  <td align="center">0.547</td>
  <td align="center">0.390</td>
  <td align="center">-0.35</td>
  <td align="center">+0.833</td>
  <td align="center">5.15</td>
  <td align="center">1.18</td>
  <td align="center">1.324</td>
  <td align="center">2.542</td>
  <td align="center">+0.761</td>
  <td align="center" bgcolor="pink">28.9116</td>
</tr>
<tr>
  <td align="center">Analytic22</td>
  <td align="center">2</td>
  <td align="center">2</td>
  <td align="center">0.887</td>
  <td align="center">1.799</td>
  <td align="center">1.022</td>
  <td align="center">0.799</td>
  <td align="center">0.583</td>
  <td align="center">+0.776</td>
  <td align="center">-0.914</td>
  <td align="center">1.25</td>
  <td align="center">46</td>
  <td align="center">1.146</td>
  <td align="center">1.444</td>
  <td align="center">-0.878</td>
  <td align="center" bgcolor="pink">34.9155</td>
</tr>
<tr>
  <td align="center">Analytic31</td>
  <td align="center">3</td>
  <td align="center">1</td>
  <td align="center">0.406</td>
  <td align="center">1.180</td>
  <td align="center">1.009</td>
  <td align="center">0.378</td>
  <td align="center">0.118</td>
  <td align="center">-0.390</td>
  <td align="center">-0.964</td>
  <td align="center">5.56</td>
  <td align="center">111</td>
  <td align="center">1.194</td>
  <td align="center">3.568</td>
  <td align="center">-0.180</td>
  <td align="center" bgcolor="pink">12.1770</td>
</tr>
<tr>
  <td align="center">Analytic32</td>
  <td align="center">3</td>
  <td align="center">2</td>
  <td align="center">0.812</td>
  <td align="center">2.327</td>
  <td align="center">4.216</td>
  <td align="center">0.690</td>
  <td align="center">0.517</td>
  <td align="center">+1.281</td>
  <td align="center">+2.051</td>
  <td align="center">0.754</td>
  <td align="center">0.311</td>
  <td align="center">1.232</td>
  <td align="center">1.813</td>
  <td align="center">+9.654</td>
  <td align="center" bgcolor="pink">169.0733</td>
</tr>
</table>
</div>


<div align="center">
<table border="1" cellpadding="5" align="center">
<tr>
  <td align="center" colspan="6"><math>~\sigma_c^2</math><br />
Determined from Numerical Integration by Enforcing Boundary Condition (B.C.)<br />
<math>\frac{d \ln x_\mathrm{env}}{d\ln \xi}\biggr|_\mathrm{surface} = c_0(2\ell+1) + \ell(3\ell+5)</math>
  </td>
</tr>
<tr>
  <th align="center">Model</th>
  <td align="center"><math>~\ell=0</math><br />&nbsp;<br />'''B.C.:''' &nbsp; <math>~c_0</math></td>
  <td align="center"><math>~\ell=1</math><br />&nbsp;<br />'''B.C.:''' &nbsp; <math>~(3c_0+8)</math></td>
  <td align="center"><math>~\ell=2</math><br />&nbsp;<br />'''B.C.:''' &nbsp; <math>~(5c_0+22)</math></td>
  <td align="center"><math>~\ell=3</math><br />&nbsp;<br />'''B.C.:''' &nbsp; <math>~7(c_0+6)</math></td>
  <td align="center"><math>~\ell=4</math><br />&nbsp;<br />'''B.C.:''' &nbsp; <math>~(9c_0+68)</math></td>
</tr>
<tr>
  <td align="center">Analytic21</td>
  <td align="center">-0.76017962</td>
  <td align="center">+9.881793</td>
  <td align="center" bgcolor="pink">+28.9116</td>
  <td align="center">+56.32919</td>
  <td align="center">-</td>
</tr>
<tr>
  <td align="center">Analytic22</td>
  <td align="center">-4.890477</td>
  <td align="center">+5.5864115</td>
  <td align="center" bgcolor="pink">+34.91550</td>
  <td align="center">+83.09678</td>
  <td align="center">-</td>
</tr>
</table>
</div>

==Broad Application==

Here's one key lesson that can be drawn from our analytically specified oscillation modes.  As has been [[#Core_Segment|documented above]], the quantum number, <math>~j</math>, associated with the eigenvector of the core is related to the oscillation frequency, 
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\sigma_c^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{3\omega^2}{2\pi G \rho_c} \, ,</math>
  </td>
</tr>
</table>
</div>
via the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\sigma_c^2}{\gamma_c} - 2\alpha_c</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 2j(2j+5) </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~\sigma_c^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{8[\alpha_c + j(2j+5)]}{3-\alpha_c} \, .</math>
  </td>
</tr>
</table>
</div>
Also, as [[#Envelope_Segment|documented above]], the quantum number, <math>~\ell</math>, associated with the eigenvector of the envelope is related to the oscillation frequency via the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\sigma_c^2}{\gamma_e(\rho_e/\rho_c)} - 2\alpha_e</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(c_0 + 3\ell)(c_0 + 3\ell + 5) </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~\sigma_c^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\gamma_e \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[2\alpha_e + (c_0 + 3\ell)(c_0 + 3\ell + 5) \biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{\rho_e}{\rho_c}\biggr) \frac{4[2c_0^2 + 4c_0 + c_0^2 + c_0(6\ell+5) + 3\ell(3\ell+5) ]}{3-\alpha_e}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{\rho_e}{\rho_c}\biggr) \frac{12[c_0^2 + c_0(2\ell+3) + \ell(3\ell+5) ]}{3-\alpha_e} \, ,</math>
  </td>
</tr>
</table>
</div>
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 4 February 2017:  This expression for the density ratio is a necessary but not sufficient relationship.  According to each analytic solution, once the pair of quantum numbers <math>~(\ell,j)</math> has been specified, the two adiabatic indexes cannot be specified independently of one another.]]where, as a reminder, <math>~c_0 = [-1 \pm \sqrt{1+\alpha_e}]</math>.  Eliminating <math>~\sigma_c^2</math> between these two relations gives,
<div align="center">
<table border="1" cellpadding="8" align="center"><tr><td 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{2(3-\alpha_e) [\alpha_c + j(2j+5)]}{3(3-\alpha_c) [c_0^2 + c_0(2\ell+3) + \ell(3\ell+5) ]} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{(3-\alpha_e) \sigma_c^2}{12 [c_0^2 + c_0(2\ell+3) + \ell(3\ell+5) ]} \, .</math>
  </td>
</tr>
</table>
</td></tr></table>
</div>

===Case A===
Consider, first, the astrophysical system in which <math>~(n_c, n_e) = (5, 1)</math>.  Given that, in general, <math>~\alpha = (3-n)/(1+n)</math>, we therefore will be considering the system in which, <math>~(\alpha_c, \alpha_e) = (-\tfrac{1}{3}, 1) </math>.  For this system, we should be able to analytically specify eigenvectors having the properties specified in the following table.

<table border="1" cellpadding="5" align="center">
<tr>
  <td colspan="4" align="center">
Analytically Specifiable Eigenvectors for which<br />
<math>~(n_c, n_e) = (5, 1)</math> &nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math> &nbsp; &nbsp; &nbsp; <math>~(\alpha_c, \alpha_e) = (-\tfrac{1}{3}, 1)</math>
  </td>
</tr>
<tr>
  <td colspan="2">&nbsp;</td>
  <td colspan="2" align="center"><math>~\frac{\rho_e}{\rho_c}</math></td>
</tr>
<tr>
  <td align="center"><math>~\ell</math></td>
  <td align="center"><math>~c_0</math></td>
  <td colspan="1" align="center">
<math>~j=1</math> <br /><math>~(\sigma_c^2 = 16)</math>
  </td>
  <td colspan="1" align="center">
<math>~j=2</math> <br /><math>~(\sigma_c^2 = 42.4)</math>
  </td>
</tr>
<tr>
  <td align="center" rowspan="2">1</td>
  <td align="center"><math>~-1 + \sqrt{2}</math></td>
  <td align="center">0.26034953</td>
  <td align="center">n/a <math>~(q > 1)</math></td>
</tr>
<tr>
  <td align="center"><math>~-1 - \sqrt{2}</math></td>
  <td align="center">1.52 (n/a)</td>
  <td align="center">(n/a)</td>
</tr>
<tr>
  <td align="center" rowspan="2">2</td>
  <td align="center"><math>~-1 + \sqrt{2}</math></td>
  <td align="center">0.10636430</td>
  <td align="center">0.28186540</td>
</tr>
<tr>
  <td align="center"><math>~-1 - \sqrt{2}</math></td>
  <td align="center">0.24400066</td>
  <td align="center">0.64660175</td>
</tr>
<tr>
  <td align="center" rowspan="2">3</td>
  <td align="center"><math>~-1 + \sqrt{2}</math></td>
  <td align="center">0.05809795</td>
  <td align="center">0.15395957</td>
</tr>
<tr>
  <td align="center"><math>~-1 - \sqrt{2}</math></td>
  <td align="center">0.10216916</td>
  <td align="center">0.27074827</td>
</tr>
</table>

=Related Discussions=
* [[User:Tohline/Appendix/Ramblings/Additional_Analytically_Specified_Eigenvectors_for_Zero-Zero_Bipolytropes#Searching_for_Additional_Eigenvectors_of_Zero-Zero_Bipolytropes|Searching for Additional, Analytically Specified Eigenvectors of Zero-Zero Bipolytropes]]
* [[User:Tohline/Appendix/Ramblings/NumericallyDeterminedEigenvectors#Numerically_Determined_Eigenvectors_of_a_Zero-Zero_Bipolytrope|Numerically Determined Eigenvectors]]

<!-- OLD STUFF ...

<div align="center">
<math>~(\ell,j) = (2, 1)</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_e = 1.194</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_c = 2.254</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\frac{3\omega^2}{2\pi G\rho_c} = 37.08874</math><br />
&nbsp; &nbsp; <br />
<math>~(\ell,j) = (2, 2)</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_e =1.209 </math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_c = 1.046</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.95210</math><br />
&nbsp; &nbsp; <br />
<math>~(\ell,j) = (3, 1)</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_e = 1.344</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_c = 1.023</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\frac{3\omega^2}{2\pi G\rho_c} = 12.452545</math><br />
&nbsp; &nbsp; <br />
<math>~(\ell,j) = (3, 2)</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_e = 1.056</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_c = 1.025</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\frac{3\omega^2}{2\pi G\rho_c} = 35.05461</math><br />
&nbsp; &nbsp; <br />
<math>~(\ell,j) = (3, 3)</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_e = 1.840</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\gamma_c = 1.325</math>&nbsp; &nbsp; <math>~\biggl|</math>&nbsp; &nbsp; <math>~\frac{3\omega^2}{2\pi G\rho_c} = 87.41594</math><br />
&nbsp; &nbsp; <br />
</div>

<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>
-->

{{ SGFfooter }}