Difference between revisions of "User:Tohline/SSC/Stability/n3PolytropeLAWE"

From VistrailsWiki
Jump to navigation Jump to search
 
(41 intermediate revisions by the same user not shown)
Line 4: Line 4:
{{LSU_HBook_header}}
{{LSU_HBook_header}}


==Groundwork==
==Background==
 
===Our Formulation of the Problem===


In an [[User:Tohline/SSC/Perturbations#2ndOrderODE|accompanying discussion]], we derived the so-called,
In an [[User:Tohline/SSC/Perturbations#2ndOrderODE|accompanying discussion]], we derived the so-called,
Line 70: Line 72:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~- \frac{3\omega^2}{4\pi G \rho_c}  \, .</math>
<math>~\frac{3\omega^2}{2\pi G \rho_c}  \, .</math>
   </td>
   </td>
</tr>
</tr>
Line 76: Line 78:
</div>
</div>


Here we perform a numerical integration of the governing LAWE for <math>~n=3</math> polytropes.  We can directly compare our results with [http://adsabs.harvard.edu/abs/1941ApJ....94..245S Schwarzschild's (1941)] published work on "Overtone Pulsations for the Standard [Stellar] Model."
In this chapter we carry out a numerical integration of this governing LAWE for <math>~n=3</math> polytropes.  The results are presented [[#Numerical_Integration|below]].
 
===Schwarzschild (1941)===
We can directly compare our results with [http://adsabs.harvard.edu/abs/1941ApJ....94..245S Schwarzschild's (1941)] published work on "Overtone Pulsations for the Standard [Stellar] Model." To begin with, it is straightforward to demonstrate that the last form of the LAWE, provided above, matches equation (2) from [[User:Tohline/SSC/Perturbations#Schwarzschild_.281941.29|Schwarzschild (1941)]], if <math>~n</math> is set to 3 &#8212; see the boxed-in excerpt, immediately below.  Note as well that Schwarzschild's dimensionless oscillation frequency &#8212; defined in his equation (1) and which we will label, <math>~\omega_\mathrm{Sch}</math> &#8212; is related to our dimensionless frequency via the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\sigma_c^2</math>
  </td>
  <td align="center">
<math>~~\leftrightarrow~~</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{3\gamma_g}{2} \biggr) \omega_\mathrm{Sch}^2 \, .</math>
  </td>
</tr>
</table>
</div>
[http://adsabs.harvard.edu/abs/1941ApJ....94..245S Schwarzschild (1941)] numerically integrated the LAWE for <math>~n=3</math> polytropic spheres to find eigenvectors (''i.e.,'' the spatially discrete eigenfunction and corresponding eigenfrequency) for five separate oscillation modes (the fundamental mode, plus the 1<sup>st</sup>, 2<sup>nd</sup>, 3<sup>rd</sup>, and 4<sup>th</sup> overtones) for models having four different adopted adiabatic indexes <math>~\gamma_g = \tfrac{4}{3}, \tfrac{10}{7}, \tfrac{20}{13}, \tfrac{5}{3})</math>. 
 
 
<div align="center">
<table border="1" cellpadding="5">
<tr><td align="center">
Paragraph extracted from [http://adsabs.harvard.edu/abs/1941ApJ....94..245S M. Schwarzschild (1941)]<p></p>
"''Overtone Pulsations for the Standard Model''"<p></p>
ApJ, vol. 94, pp. 245 - 252 &copy; American Astronomical Society
</td></tr>
<tr><td>
[[File:Schwarzschild1941.png|700px|center|Schwarzschild (1941, ApJ, 94, 245)]]
</td></tr>
<tr><td align="left">
&nbsp; &nbsp; &nbsp; <sup>3</sup>A. S. Eddington (1930), [https://archive.org/details/TheInternalConstitutionOfTheStars ''The Internal Constitution of the Stars''], pp. 188 and 192.
</td></tr>
</table>
</div>
 
Drawing from our [[User:Tohline/SSC/Perturbations#Ensure_Finite-Amplitude_Fluctuations|discussion of the historical treatment of boundary conditions]], we presume that Schwarzschild imposed the following constraint at the surface:
 
<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}{\gamma_g} \biggl( 4 - 3\gamma_g + \frac{\omega^2 R^3}{GM_\mathrm{tot}}\biggr) </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{3\omega^2 R^3}{4\pi G \gamma \bar\rho} - \alpha \biggr] </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2} \biggl[ \frac{\sigma_c^2}{\gamma } \biggl(\frac{\rho_c}{\bar\rho}\biggr) -2 \alpha \biggr] </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2} \biggl\{ \biggl[ \mathfrak{F} + 2\alpha \biggr] \biggl(\frac{\rho_c}{\bar\rho}\biggr) -2 \alpha \biggr\} \, .</math>
  </td>
</tr>
</table>
</div>
 
Recognizing from an [[User:Tohline/SSC/Structure/Polytropes#Horedt2004|accompanying tabulation]] that, for <math>~n=3</math> polytropes,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\rho_c}{\bar\rho}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~54.18248 \, ,</math>
  </td>
</tr>
</table>
</div>
we presume that the surface boundary condition imposed by Schwarzschild was,
 
<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>~27.09124 ( \mathfrak{F} + 2\alpha )  - \alpha  \, .</math>
  </td>
</tr>
</table>
</div>
 
Our Table 1 catalogs the eigenfrequencies that Schwarzschild determined (drawn from ''his'' Table 1) for these twenty different models/modes.
 
<table border="1" cellpadding="8" align="center">
<tr>
  <th colspan="6" align="center">
<font size="+1">Table 1:</font> &nbsp; From Table 1 of [http://adsabs.harvard.edu/abs/1941ApJ....94..245S M. Schwarzschild (1941)]
  </th>
</tr>
<tr>
  <td align="center" rowspan="2">Mode</td>
  <td align="center"><math>~\alpha = 0.0</math><p></p><math>~(\gamma_g = 4/3)</math></td>
  <td align="center"><math>~\alpha = 0.2</math><p></p><math>~(\gamma_g = 10/7)</math></td>
  <td align="center" colspan="2"><math>~\alpha = 0.4</math><p></p><math>~(\gamma_g = 20/13)</math></td>
  <td align="center"><math>~\alpha = 0.6</math><p></p><math>~(\gamma_g = 5/3)</math></td>
</tr>
<tr>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
  <td align="center"><math>~\mathfrak{F} = \biggl[\frac{3\omega_\mathrm{Sch}^2}{2} - 2\alpha \biggr]</math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
</tr>
<tr>
  <td align="center">0</td>
  <td align="center">0.00000</td>
  <td align="center">0.05882</td>
  <td align="center">0.10391</td>
  <td align="center" bgcolor="pink">-0.64414</td>
  <td align="center">0.13670</td>
</tr>
<tr>
  <td align="center">1</td>
  <td align="center">0.16643</td>
  <td align="center">0.19139</td>
  <td align="center">0.21998</td>
  <td align="center" bgcolor="pink">-0.47003</td>
  <td align="center">0.25090</td>
</tr>
<tr>
  <td align="center">2</td>
  <td align="center">0.3392</td>
  <td align="center">0.3648</td>
  <td align="center">0.3920</td>
  <td align="center" bgcolor="pink">-0.2120</td>
  <td align="center">0.4209</td>
</tr>
<tr>
  <td align="center">3</td>
  <td align="center">0.5600</td>
  <td align="center">0.5863</td>
  <td align="center">0.6136</td>
  <td align="center" bgcolor="pink">+0.1204</td>
  <td align="center">0.6420</td>
</tr>
<tr>
  <td align="center">4</td>
  <td align="center">0.8283</td>
  <td align="center">0.8554</td>
  <td align="center">0.8832</td>
  <td align="center" bgcolor="pink">+0.5248</td>
  <td align="center">0.9117</td>
</tr>
</table>
 
[http://adsabs.harvard.edu/abs/1941ApJ....94..245S Schwarzschild (1941)] also documented the radial structure of the eigenfunction that is associated with each of these twenty model/mode eigenfrequencies.  Each column of his Table 4, except the first, presents numerical values of the amplitude of a specific model/mode at 84 discrete radial locations throughout the n = 3 polytrope; the first column of the table lists the corresponding radial coordinate, <math>~\xi</math>.  Focusing on the model that he analyzed assuming <math>~\alpha = 0.4</math>, we have typed his five columns of data into an Excel spreadsheet and have used this data to generate the pair of plots displayed, below, in Figure 1.  The left-hand panel displays the eigenfunction amplitude versus radius, <math>~x(\xi)</math>, for the fundamental mode as well as for the first four overtones; it essentially replicates Figure 1 from [http://adsabs.harvard.edu/abs/1941ApJ....94..245S Schwarzschild (1941)].  The right-hand panel displays the same data, but as a semi-log plot; specifically, it displays <math>~y(\xi)</math>, where,
<div align="center">
<math>~y \equiv \frac{1}{2} \log_{10}[x^2 + 10^{-8}] \, .</math>
</div>
Each sharp valley in this semi-log plot highlights the location of a node in the corresponding eigenfunction, that is, it identifies where <math>~x(\xi)</math> crosses through zero.
 
 
<table border="1" cellpadding="5" align="center">
<tr>
  <th align="center"><font size="+1">Figure 1:</font>&nbsp;
Schwarzschild's Eigenfunctions for an n = 3 Polytrope with &alpha; = 0.4 (&gamma; = 20/13)
  </th>
</tr>
<tr>
<td align="center">[[File:Schwarzschild1941Combined3small.png|800 px|Schwarzschild (1941) eigenfunctions]]</td>
</tr>
</table>


==Numerical Integration==
==Numerical Integration==


It can be shown straightforwardly that this matches the LAWE used by [[User:Tohline/SSC/Perturbations#Schwarzschild_.281941.29|Schwarzschild (1941)]], if <math>~n</math> is set to 3.  Here we use the [[User:Tohline/SSC/Stability/Polytropes#Numerical_Integration_from_the_Center.2C_Outward|finite-difference algorithm described separately]] to integrate the discretized LAWE from the center of the polytropic configuration, outward to its surface, which in this case &#8212; see, for example, [[User:Tohline/SSC/Structure/Polytropes#Horedt2004|p. 77 of Horedt (2004)]] &#8212; is located at the polytropic-coordinate location,
===From the Core to the Surface===
 
Here we use the [[User:Tohline/SSC/Stability/Polytropes#Numerical_Integration_from_the_Center.2C_Outward|finite-difference algorithm described separately]] to integrate the discretized LAWE from the center of the polytropic configuration, outward to its surface, which in this case &#8212; see, for example, [[User:Tohline/SSC/Structure/Polytropes#Horedt2004|p. 77 of Horedt (2004)]] &#8212; is located at the polytropic-coordinate location,
<div align="center">
<div align="center">
<math>~\xi_\mathrm{max} = 6.89684862 \, .</math>
<math>~\xi_\mathrm{max} = 6.89684862 \, .</math>
Line 96: Line 310:
   </li>
   </li>
   <li>Set the eigenfunction to unity at the center <math>~(\xi_0 = 0)</math> of the configuration, that is, set <math>~x_0 = 1</math>.</li>
   <li>Set the eigenfunction to unity at the center <math>~(\xi_0 = 0)</math> of the configuration, that is, set <math>~x_0 = 1</math>.</li>
   <li>Determine the value of the eigenfunction at the first grid ''line'' away from the center &#8212; having coordinate location, <math>~\xi_1 = \Delta_\xi </math> &#8212; via the expression,
   <li>Determine the value of the eigenfunction at the first grid ''line'' away from the center &#8212; having coordinate location, <math>~\xi_1 = \Delta_\xi </math> &#8212; via the [[User:Tohline/Appendix/Ramblings/PowerSeriesExpressions#PolytropicDisplacement|derived power-series expression]],
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 111: Line 325:
   <td align="left">
   <td align="left">
<math>~
<math>~
x_0 \biggl[ 1 - \frac{\Delta_\xi^2 (n+1) \mathfrak{F}}{12} \biggr] \, .</math>
x_0 \biggl[ 1 - \frac{\Delta_\xi^2 (n+1) \mathfrak{F}}{60} \biggr] \, .</math>
   </td>
   </td>
</tr>
</tr>
Line 141: Line 355:
</ul>
</ul>


Now, in searching for an appropriate boundary condition at the surface of the configuration, it will be useful to tabulate, not only the value of the eigenfunction at the surface, <math>~x_N</math>, but also its logarithmic derivative.  A finite-difference expression of this logarithmic derivative that is consistent with the above-described finite-difference algorithm, is,
We divided our model into <math>~N = 200</math> radial zones and, using this algorithm, integrated the LAWE from the center of the configuration to the surface, for <math>~\alpha = 0.4</math>, and approximately 40 different chosen values of the frequency parameter across the range, <math>~-0.7 \le \mathfrak{F} \le + 0.3</math>.  The radial displacement functions resulting from these integrations are presented in Figure 2 as an animation sequence.  The specified value of <math>~\mathfrak{F}</math> is displayed at the top of each animation frame, and the resulting displacement function, <math>~x(r/R)</math>, is traced by the small, red circular markers in each frame.
 
<div align="center" id="SchwarzschildMovie">
<table border="0" cellpadding="5" align="center">
<tr>
  <th colspan="1" align="center"><font size="+1">Figure 2:</font>&nbsp; Numerically Determined Eigenfunctions for Various <math>~\mathfrak{F}</math></th>
  <td align="center" rowspan="2">
<table border="1" cellpadding="5">
<tr>
  <th align="center" colspan="3"><font size="+1">Table 2</font></th>
</tr>
<tr>
  <td align="center" rowspan="2">Mode</td>
  <td align="center" colspan="1">Match Schwarzschild</td>
  <td align="center" colspan="1">Match B.C.</td>
</tr>
<tr>
  <td align="center"><math>~\mathfrak{F}</math></td>
  <td align="center"><math>~\mathfrak{F}</math></td>
</tr>
<tr>
  <td align="center">0</td>
  <td align="center">-0.644131578154</td>
  <td align="center">-0.644131577959</td>
</tr>
<tr>
  <td align="center">1</td>
  <td align="center">-0.47013976423</td>
  <td align="center">-0.47013975308</td>
</tr>
<tr>
  <td align="center">2</td>
  <td align="center">-0.2121284391</td>
  <td align="center">-0.2121282667</td>
</tr>
<tr>
  <td align="center">3</td>
  <td align="center">+0.1202565375</td>
  <td align="center">+0.120257856</td>
</tr>
</table>
</td>
</tr>
<tr>
<td align="center" bgcolor="white">
[[File:Schwarzschild1941movie.gif|Eigenfunctions for Standard Model]]
</td>
</tr>
</table>
</div>
 
Each frame of the Figure 2 animation also displays, as smooth solid curves, the radial eigenfunctions that Schwarzschild (1941) obtained for the fundamental mode (blue curve) and the first three overtone modes (green, purple, &amp; orange curves, repectively) for his model with <math>~\alpha = 0.4</math>.  These are the same curves that appear in the left-hand panel of Figure 1, but here the displacement amplitude has been renormalized such that <math>~x_0 = 1</math>, and, along the horizontal axis, the radial location is marked in terms of the fractional radius, <math>~r/R \equiv \xi/\xi_\mathrm{max}</math>.  In our examination of this model, as we approached each specific value of a modal eigenfrequency identified by Schwarzschild &#8212; see the frequencies highlighted in pink in our Table 1 &#8212; we fine-tuned our choice of the eigenfrequency in order to find a displacement function whose surface ''amplitude'' matched, to a high level of precision, the surface amplitude associated with Schwarzschild's corresponding published eigenfunction.  The column of our Table 2 whose heading is "Match Schwarzschild" identifies &#8212; to at least 10 digits precision &#8212; the frequency choice that was required in order for these surface ''amplitudes'' to match in each case. 
 
It is gratifying to see that our resulting frequencies match well the values published by Schwarzschild (as highlighted in pink, above).  But this does not satisfactorily explain why, among the entire range of displacement functions displayed (in red) in the Figure 2 animation, Schwarzschild labeled these specific ones as the eigenmodes.  As we shall now demonstrate, his eigenmode identifications resulted from the imposition of a specific, physically justified constraint on the ''slope'', rather than the value, of the displacement function at the surface of the configuration.  (See also our separate brief answer to the question, "[[User:Tohline/SSC/Perturbations#What_Makes_This_an_Eigenvalue_Problem.3F|What makes this an eigenvalue problem?]]".)
 
===Surface Boundary Condition===
 
As was stated, [[#Schwarzschild_.281941.29|above]], we presume that as Schwarzschild searched for natural modes of oscillation in isolated, <math>~n=3</math> polytropes, he imposed the following boundary condition at the surface of the configuration:
 
<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>~27.09124 ( \mathfrak{F} + 2\alpha ) - \alpha  \, .</math>
  </td>
</tr>
</table>
</div>
 
In order to duplicate his findings, then, we need to fine tune our specification of the oscillation frequency such that the resulting displacement function presents this behavior at the surface of our model. A finite-difference expression of this logarithmic derivative that is consistent with the above-described finite-difference algorithm, is,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 175: Line 465:
+ x_{N-1} \biggl[\frac{4\Delta_\xi \theta_{N}}{\xi_\mathrm{max}} - \Delta_\xi (n+1)(- \theta^')_{N} - 2\theta_{N}\biggr]  
+ x_{N-1} \biggl[\frac{4\Delta_\xi \theta_{N}}{\xi_\mathrm{max}} - \Delta_\xi (n+1)(- \theta^')_{N} - 2\theta_{N}\biggr]  
\, .</math>
\, .</math>
  </td>
</tr>
</table>
</div>
We added to our numerical algorithm a step that evaluates, in this manner, the logarithmic derivative of the displacement function at the surface of our polytropic configuration.  The eigenfrequencies that generated displacement functions with this surface behavior are listed for four separate modes in the column of Table 2 titled, "Match B.C."  In every case the values agree to at least five decimal places with the "Match Schwarzschild" eigenfrequencies.  We conclude, therefore, that it ''was'' the implementation of this surface boundary condition that permitted Schwarzschild to quantitatively identify the properties of the eigenvectors associated with natural radial modes of oscillation in <math>~n=3</math> polytropes.
===Our Results===
<table border="1" cellpadding="8" align="center">
<tr>
  <th colspan="9" align="center">
<font size="+1">Table 3:</font> &nbsp; Our Results (to be compared w/ Table 1, above)
  </th>
</tr>
<tr>
  <td align="center" rowspan="2">Mode</td>
  <td align="center" colspan="2"><math>~\alpha = 0.0</math><p></p><math>~(\gamma_g = 4/3)</math></td>
  <td align="center" colspan="2"><math>~\alpha = 0.2</math><p></p><math>~(\gamma_g = 10/7)</math></td>
  <td align="center" colspan="2"><math>~\alpha = 0.4</math><p></p><math>~(\gamma_g = 20/13)</math></td>
  <td align="center" colspan="2"><math>~\alpha = 0.6</math><p></p><math>~(\gamma_g = 5/3)</math></td>
</tr>
<tr>
  <td align="center"><math>~\mathfrak{F}</math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2 = \frac{2}{3}\biggl(\mathfrak{F}+2\alpha  \biggr)</math>
  <td align="center"><math>~\mathfrak{F}</math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
  <td align="center"><math>~\mathfrak{F} </math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
  <td align="center"><math>~\mathfrak{F}</math>
  <td align="center"><math>~\omega_\mathrm{Sch}^2</math>
</tr>
<tr>
  <td align="center">0</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">-0.311782342981</td>
  <td align="center">0.058812</td>
  <td align="center" bgcolor="pink">-0.644131577959</td>
  <td align="center">0.103912</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
<tr>
  <td align="center">1</td>
  <td align="center">+0.24946512002</td>
  <td align="center">0.166310</td>
  <td align="center">-0.113086698932</td>
  <td align="center">0.191276</td>
  <td align="center" bgcolor="pink">-0.47013975308</td>
  <td align="center">0.219907</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
<tr>
  <td align="center">2</td>
  <td align="center">+0.50882623652</td>
  <td align="center">0.339217</td>
  <td align="center">+0.14705874055</td>
  <td align="center">0.364706</td>
  <td align="center" bgcolor="pink">-0.2121282667</td>
  <td align="center">0.391914</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
<tr>
  <td align="center">3</td>
  <td align="center">+0.83977118</td>
  <td align="center">0.559847</td>
  <td align="center">+0.479241829</td>
  <td align="center">0.586161</td>
  <td align="center" bgcolor="pink">+0.120257856</td>
  <td align="center">0.613505</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
<tr>
  <td align="center">4</td>
  <td align="center">+1.24253191</td>
  <td align="center">0.828355</td>
  <td align="center">+0.8832297</td>
  <td align="center">0.855486</td>
  <td align="center" bgcolor="pink">+0.52498863</td>
  <td align="center">0.883326</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
</table>
=Truncated n = 3 Polytropes=
We understand that, for ''isolated'' <math>~n=3</math> polytropic spheres, the value of the adiabatic exponent for which the configuration is marginally unstable is <math>~\gamma_g = 4/3</math>, which is equivalent to, <math>~\alpha \equiv (3-4/\gamma_g) = 0</math>.  This critical condition is identified by examining when the oscillation frequency of the ''fundamental'' mode goes to zero.  Let's use our numerical integration tool to determine what this critical value of the adiabatic exponent is for ''truncated'', <math>~n=3</math> polytropes.  We will accomplish this as follows:
* At various truncation radii, <math>~0 < \xi_\mathrm{surf}/\xi_\mathrm{max} < 1</math> &hellip;
** Force <math>~\mathfrak{F} = -2\alpha</math>;
** Iterate on the choice of <math>~\alpha</math> until the displacement function with no radial nodes (''i.e.,'' the ''fundamental'' mode) satisfies the surface boundary condition of <math>~(d\ln x/d\ln \xi)_\mathrm{surf} = -3</math>, to a desired level of accuracy.
The following table shows the values of <math>~\alpha_\mathrm{crit}</math> &#8212; and associated values of <math>~\gamma_\mathrm{crit}</math> &#8212; that we obtained for nine different values of <math>~\xi_\mathrm{surf}/\xi_\mathrm{max}</math>; in each case, iterations were continued until the desired surface boundary condition was satisfied to six significant digits.
<div align="center">
<table border="1" cellpadding="5" align="center">
<tr>
  <td align="center" colspan="6">Pressure-Truncated n = 3 Polytropes<br /><math>~N_\mathrm{zones} = 200</math>
</tr>
<tr>
  <td align="center">Edge Zone</td>
  <td align="center"><math>~\frac{\xi_\mathrm{surf}}{\xi_\mathrm{max}}</math></td>
  <td align="center"><math>~\alpha_\mathrm{crit}</math></td>
  <td align="center"><math>~\mathfrak{F}</math></td>
  <td align="center">Surface B.C</td>
  <td align="center"><math>~\gamma_\mathrm{crit}</math></td>
</tr>
<tr>
  <td align="center">200</td>
  <td align="center">1.00</td>
  <td align="center">0</td>
  <td align="center">0</td>
  <td align="center">---</td>
  <td align="center"><math>~\tfrac{4}{3}</math></td>
</tr>
<tr>
  <td align="center">180</td>
  <td align="center">0.90</td>
  <td align="center">-0.000201541</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">1.333244</td>
</tr>
<tr>
  <td align="center">160</td>
  <td align="center">0.8</td>
  <td align="center">-0.00327575</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">1.331879</td>
</tr>
<tr>
  <td align="center">150</td>
  <td align="center">0.75</td>
  <td align="center">-0.00808603</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">1.329749</td>
</tr>
<tr>
  <td align="center">120</td>
  <td align="center">0.60</td>
  <td align="center">-0.0576031</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">1.308214</td>
</tr>
<tr>
  <td align="center">100</td>
  <td align="center">0.5</td>
  <td align="center">-0.159111</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">1.266179</td>
</tr>
<tr>
  <td align="center">80</td>
  <td align="center">0.40</td>
  <td align="center">-0.405712</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">1.174497</td>
</tr>
<tr>
  <td align="center">50</td>
  <td align="center">0.25</td>
  <td align="center">-1.74909</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">0.842266</td>
</tr>
<tr>
  <td align="center">20</td>
  <td align="center">0.10</td>
  <td align="center">-14.6648</td>
  <td align="center"><math>~-2\alpha_\mathrm{crit}</math></td>
  <td align="center">-3.00000</td>
  <td align="center">0.226439</td>
</tr>
</table>
</div>
=Analytic Inquiry=
<font color="red"><b>NOTE</b> (from J. E. Tohline in April, 2017):</font> &nbsp; The following subsections present some exploratory ideas that were pursued while I was searching for analytic solutions to the polytropic LAWE.  For the most part this material has been superseded by a [[User:Tohline/SSC/Stability/InstabilityOnsetOverview#Polytropic|separate discussion]] in which we describe the desired analytic solution, which we discovered in March, 2017.
==Fundamental-Mode, Homentropic Oscillations==
The LAWE, [[#Background|presented above]], that is relevant to polytropic spheres, may be rewritten as,
<div align="center">
{{ User:Tohline/Math/EQ_RadialPulsation02 }}
</div>
Now, if we assume that oscillations occur adiabatically with an adiabatic index that is consistent with the chosen polytropic index &#8212; that is to say,
<div align="center">
<math>~\gamma_g = \frac{n+1}{n}</math>&nbsp; &nbsp; &nbsp; <math>~\Rightarrow</math> &nbsp; &nbsp; &nbsp; <math>~\alpha =  \frac{3-n}{n+1} \, ,</math>
</div>
in which case the configuration remains homentropic as it oscillates &#8212; and if we look only for a (marginally unstable) configuration that has <math>~\sigma_c^2 = 0</math>, then the relevant 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>~
\frac{d^2x}{d\xi^2} + \biggl[ 4 - (n+1) Q \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} 
- (3-n)  Q \frac{x}{\xi^2} \, .
</math>
  </td>
</tr>
</table>
</div>
==Specific case of n = 3 Polytropes==
===Homologous Collapse===
If we examine only an <math>~n=3</math> polytropic configuration, then the last term disappears.  This means that ''in this very special case'', a perfectly valid solution to the LAWE is <math>~x = \mathrm{constant}</math>.  This is presumably the eigenfunction that Schwarzschild deduced; the fundamental-mode "oscillations" are perfectly homologous. Given that the model is marginally unstable, an ensuing dynamical collapse will presumably begin in a perfectly homologous fashion.  This is precisely the type of "free-fall" collapse that was discussed and modeled by [[User:Tohline/Apps/GoldreichWeber80#Homologous_Solution|Goldreich &amp; Weber (1980)]].
===Another Potential Option===
We have wondered whether, in this very special case, one or more additional fundamental-mode eigenfunction(s) might satisfy the governing LAWE.  Here is a relevant line of arguments, beginning with the LAWE for the n = 3 polytropic sphere.
<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{d^2x}{d\xi^2} + 4( 1 +  Q ) \frac{1}{\xi} \cdot \frac{dx}{d\xi} 
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{dy}{d\xi} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- 4( 1 +  Q ) \frac{y}{\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{1}{4}\cdot \frac{d\ln y}{d\ln\xi} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- ( 1 +  Q ) \, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<math>~y \equiv \frac{dx}{d\xi} \, .</math>
</div>
But, by definition, the function <math>~Q(\xi)</math> is a logarithmic derivative of the Lane-Emden function.  Hence, we also can write,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~-1
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{1}{4}\cdot \frac{d\ln y}{d\ln\xi} + \frac{d\ln\theta}{d\ln\xi}</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{d\ln (\theta y^{1/4})}{d\ln\xi} </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~- d\ln\xi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ d\ln (\theta y^{1/4}) \, . </math>
  </td>
</tr>
</table>
</div>
Integrating this equation once gives,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\ln(\theta y^{1/4}) + \ln\xi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ln (c_0)</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \xi \theta y^{1/4} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~c_0</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\frac{dx}{d\xi} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{c_0}{\xi\theta}\biggr)^4 \, .</math>
  </td>
</tr>
</table>
</div>
Referring to the [[User:Tohline/Appendix/Ramblings/PowerSeriesExpressions#PolytropicLaneEmden|power-series expansion of the polytropic Lane-Emden function]], <math>~\theta(\xi)</math>, about the configuration's center, we see that the product, <math>~\xi\theta</math>, goes to zero as the first power of <math>~\xi</math>.  This means that the right-hand side of this last differential equation blows up at the center.  This, therefore, does not appear to provide a physically viable avenue by which to identify an alternative fundamental-mode eigenfunction.
==Play With Form of LAWE==
===Logarithmic Derivative Rewrite===
We have noticed that the LAWE that governs the eigenfunction associated with the fundamental mode of the marginally unstable model (FMMUM),
<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{d^2x}{d\xi^2} + \biggl[ 4 - (n+1) Q \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} 
- (3-n)  Q \frac{x}{\xi^2} \, ,
</math>
  </td>
</tr>
</table>
</div>
may be rewritten entirely as an expression that relates the logarithmic derivatives of <math>~x, \xi,</math> and <math>~\theta</math>.  Multiplying through by <math>~\xi^2/x</math>, then drawing on a differential relation that has been [[User:Tohline/SSC/Stability/BiPolytrope0_0Details#Idea_Involving_Logarithmic_Derivatives|derived in a separate context]], namely,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\xi^2}{x} \cdot \frac{d^2x}{d\xi^2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr]
+ \biggl[  \frac{d\ln x}{d\ln \xi} -1 \biggr]\cdot \frac{d\ln x}{d\ln \xi} \, ,
</math>
</td>
</tr>
</table>
</div>
this LAWE associated with the FMMUM 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>~
\frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr]
+ \biggl[  \frac{d\ln x}{d\ln \xi} -1 \biggr]\cdot \frac{d\ln x}{d\ln \xi} + \biggl[ 4 - (n+1) Q \biggr] \frac{d\ln x}{d\ln\xi}  - (3-n)  Q
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr]
+ \frac{d\ln x}{d\ln \xi} \cdot \biggl[  \frac{d\ln x}{d\ln \xi} + 3 - (n+1) Q \biggr]  - (3-n)  Q
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr]
+ \frac{d\ln x}{d\ln \xi} \cdot \biggl\{  \frac{d\ln [x \xi^3 \theta^{(n+1)}] }{d\ln \xi} \biggr\}  + \frac{d\ln \theta^{(3-n) }}{d\ln\xi} \, .
</math>
  </td>
</tr>
</table>
</div>
I'm not sure if anyone else has previously appreciated that the "fundamental mode" polytropic LAWE can be written in this form.  I'm even less sure that this form sheds light on its solution.
Play a little more &hellip; &nbsp; Start by letting, <math>~A \equiv  [x \xi^3 \theta^{(n+1)}] \, ,</math> in which case we have,
<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{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr]
+ \frac{d\ln x}{d\ln \xi} \cdot \biggl\{  \frac{1}{A} \cdot \frac{dA}{d\ln \xi} \biggr\}  + \frac{d\ln \theta^{(3-n) }}{d\ln\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{A} \cdot \frac{d}{d\ln\xi} \biggl[ A \cdot \frac{d\ln x}{d\ln \xi} \biggr]  + \frac{d\ln \theta^{(3-n) }}{d\ln\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
\frac{d}{d\xi} \biggl[ A \cdot \frac{d\ln x}{d\ln \xi} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{A}{\theta^{(3-n)}} \cdot \frac{d\theta^{(3-n) }}{d\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
\frac{d}{d\xi} \biggl[ x \xi^3 \theta^{(n+1)} \cdot \frac{d\ln x}{d\ln \xi} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{x \xi^3 \theta^{(n+1)}}{\theta^{(3-n)}} \cdot \frac{d\theta^{(3-n) }}{d\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
\frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \cdot \frac{dx}{d\xi} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- x \xi^3 \theta^{2(n-1)} \cdot \frac{d\theta^{(3-n) }}{d\xi} \, .
</math>
  </td>
</tr>
</table>
</div>
===One Feeble Guess===
Now, what if, <math>~x \equiv [\xi^{-2}\theta^{-n}]</math> &nbsp; ?
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\Rightarrow~~~
\xi^2 \theta^{(n+1)} \cdot \frac{d(\xi^{-2}\theta^{-n})}{d\xi}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{d\theta}{d\xi} - \frac{1}{\xi^{2}\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi} \, ,
</math>
  </td>
</tr>
</table>
</div>
in which case, the LAWE becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\frac{d}{d\xi} \biggl[ \xi^2 \frac{d\theta}{d\xi} - \frac{1}{\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi}\biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- [\xi^{-2}\theta^{-n}] \xi^3 \theta^{2(n-1)} \cdot \frac{d\theta^{(3-n) }}{d\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
\frac{d}{d\xi} \biggl[ \xi^2 \frac{d\theta}{d\xi} \biggr] - \frac{d}{d\xi} \biggl[\frac{1}{\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi}\biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \xi \theta^{n-2} \cdot \frac{d\theta^{(3-n) }}{d\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
-\xi^2\theta^n - \frac{d}{d\xi} \biggl[\frac{1}{\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi}\biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \xi  (3-n)\frac{d\theta }{d\xi}
</math>
  </td>
</tr>
</table>
</div>
===Behavior for Known n=5 Solution===
We know that the [[User:Tohline/SSC/Stability/n5PolytropeLAWE#Eureka_Moment|FMMUM for pressure-truncated, n = 5 polytropic configurations]] takes the form,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~x </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~1- f_n(\xi) \, ,</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<math>~f_5(\xi) = \frac{\xi^2}{15} \, .</math>
</div>
The governing LAWE therefore gives,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \cdot \frac{dx}{d\xi} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- x \xi^3 \theta^{2(n-1)} \cdot \frac{d\theta^{(3-n) }}{d\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- (3-n)x \xi^3 \theta^{n} \cdot \frac{d\theta }{d\xi}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
(3-n)(1-f_n) \xi^3 \theta^{n} \cdot \frac{d\theta }{d\xi}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \cdot \frac{df_n}{d\xi} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi^4 \theta^{(n+1)} \cdot\frac{d^2f_n}{d\xi^2}
+ \frac{df_n}{d\xi} \cdot \frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi^4 \theta^{(n+1)} \cdot\frac{d^2f_n}{d\xi^2}
+ \frac{df_n}{d\xi} \cdot \biggl\{
4\xi^3 \theta^{(n+1)} + (n+1)\xi^4 \theta^n \frac{d\theta}{d\xi}
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
(3-n)(1-f_n) \frac{d\theta }{d\xi}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi \theta \cdot\frac{d^2f_n}{d\xi^2}
+ \frac{df_n}{d\xi} \cdot \biggl\{
4 \theta + (n+1)\xi \frac{d\theta}{d\xi}
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
\frac{1}{\theta} \frac{d\theta }{d\xi} \biggl[ (3-n)(1-f_n) - (n+1)\xi \frac{df_n}{d\xi} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi \cdot\frac{d^2f_n}{d\xi^2}
+ 4 \cdot \frac{df_n}{d\xi} \, .
</math>
  </td>
</tr>
</table>
</div>
Let's check to see whether the known <math>~f_5(\xi)</math> function properly satisfies this last ODE when n = 5.
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
-\xi (3+\xi^2)^{-1}  \biggl[ -2 \biggl(1-\frac{\xi^2}{15} \biggr) - \frac{12\xi^2}{15} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2\xi}{15} + \frac{8\xi}{15}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~
2 + \frac{2\xi^2}{3}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2}{3} (3+\xi^2) \, ,
</math>&nbsp; &nbsp; &nbsp; Yes!
  </td>
</tr>
</table>
</div>
Given that,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\theta_5^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{3}{3+\xi^2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\xi^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3\biggl[ \frac{1}{\theta_5^2} - 1 \biggr] \, ,</math>
  </td>
</tr>
</table>
</div>
we could ''presume'' that, when defined in terms of <math>~\theta_5</math>, the defining function,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~f_5(\theta_5) = \frac{1}{n}\biggl[ \frac{1}{\theta_5^2} - 1 \biggr]</math>
  </td>
  <td align="center">
&nbsp; &nbsp;<math>~\Rightarrow~</math>&nbsp; &nbsp;
  </td>
  <td align="left">
<math>~
\frac{df_5}{d\xi} = - \frac{2}{n\theta_5^3} \cdot \frac{d\theta_5}{d\xi}
</math>
  </td>
  <td align="center">
&nbsp; &nbsp;<math>~\Rightarrow~</math>&nbsp; &nbsp;
  </td>
  <td align="left">
<math>~
\frac{d^2f_5}{d\xi^2} = - \frac{2}{n\theta_5^3} \cdot \frac{d^2\theta_5}{d\xi^2} + \frac{6}{n\theta_5^4} \cdot \biggl( \frac{d\theta_5}{d\xi} \biggr)^2 \, .
</math>
  </td>
</tr>
</table>
</div>
In this case, 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>~
\xi \cdot\frac{d^2f_n}{d\xi^2}
- \frac{(3-n)}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ 1-f_n \biggr]
+ \biggl[ 4 + \frac{(n+1)\xi}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggr] \biggl[ \frac{df_n}{d\xi} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi \biggl\{ - \frac{2}{n\theta_5^3} \cdot \frac{d^2\theta_5}{d\xi^2} + \frac{6}{n\theta_5^4} \cdot \biggl( \frac{d\theta_5}{d\xi} \biggr)^2\biggr\}
- \biggl[ 4 + \frac{(n+1)\xi}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggr] \biggl[ \frac{2}{n\theta_5^3} \cdot \frac{d\theta_5}{d\xi}\biggr]
- \frac{(3-n)}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ \frac{(n + 1)\theta_5^2 -1  }{n\theta_5^2} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{n\theta_5^3}\biggl\{
\biggl[ 2\xi \cdot \frac{d^2\theta_5}{d\xi^2} - \frac{6\xi}{\theta_5} \cdot \biggl( \frac{d\theta_5}{d\xi} \biggr)^2\biggr]
+ 2\biggl[ 4 + \frac{(n+1)\xi}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggr]  \frac{d\theta_5}{d\xi}
+ (3-n) \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ (n + 1)\theta_5^2 -1 \biggr]
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{n\theta_5^3}\biggl\{
2\xi \cdot \frac{d^2\theta_5}{d\xi^2}
+ (n-2) \frac{2\xi}{\theta_5}\biggl( \frac{d\theta_5}{d\xi} \biggr)^2
+ \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ (3-n) (n + 1)\theta_5^2 + 5+ n  \biggr] 
\biggr\} \, .
</math>
  </td>
</tr>
</table>
</div>
Now, from the polytropic Lane-Emden equation, we also know that,
<div align="center">
{{ User:Tohline/Math/EQ_SSLaneEmden01 }}
</div>
That is,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~2\xi \biggl(\frac{d\theta_5}{d\xi}\biggr) + \xi^2 \cdot \frac{d^2\theta_5}{d\xi^2} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \xi^2 \theta_5^n</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~  \frac{d^2\theta_5}{d\xi^2} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \theta_5^n - \frac{2}{\xi} \biggl(\frac{d\theta_5}{d\xi}\biggr)</math>
  </td>
</tr>
</table>
</div>
So, the 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>~
(n-2) \frac{2\xi}{\theta_5}\biggl( \frac{d\theta_5}{d\xi} \biggr)^2
+ (n+1)\biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ (3-n) \theta_5^2 + 1 \biggr] 
- 2\xi \theta_5^n \, .
</math>
  </td>
</tr>
</table>
</div>
Again, let's check to see if the case of n=5 works &hellip;
<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( \frac{d\theta_5}{d\xi} \biggr)^2
+ \frac{\theta_5}{\xi} \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ 1 - 2 \theta_5^2 \biggr] 
- \frac{\theta_5^{6}}{3}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
3\xi^2(3+\xi^2)^{-3}
-3(3+\xi^2)^{-2}\biggl[ 1 - 6(3+\xi^2)^{-1} \biggr] 
- 3^2(3+\xi^2)^{-3}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(3+\xi^2)^{-3}\biggl\{
3\xi^2 -3\biggl[ (3+\xi^2)  -6 \biggr]  - 3^2 \biggr\}
</math> &nbsp; &nbsp; Yes!
   </td>
   </td>
</tr>
</tr>
Line 187: Line 1,518:
** n = 3:&nbsp; [[User:Tohline/SSC/Stability/n3PolytropeLAWE|Numerical Solution]] to compare with [http://adsabs.harvard.edu/abs/1941ApJ....94..245S M. Schwarzschild (1941)]
** n = 3:&nbsp; [[User:Tohline/SSC/Stability/n3PolytropeLAWE|Numerical Solution]] to compare with [http://adsabs.harvard.edu/abs/1941ApJ....94..245S M. Schwarzschild (1941)]
** n = 5:&nbsp; [[User:Tohline/SSC/Stability/n5PolytropeLAWE|Attempt at Formulating an Analytic Solution]]
** n = 5:&nbsp; [[User:Tohline/SSC/Stability/n5PolytropeLAWE|Attempt at Formulating an Analytic Solution]]
{{LSU_HBook_footer}}

Latest revision as of 16:39, 11 April 2017

Radial Oscillations of n = 3 Polytropic Spheres

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

Background

Our Formulation of the Problem

In an accompanying discussion, we derived the so-called,

Adiabatic Wave (or Radial Pulsation) Equation

LSU Key.png

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

whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations. Because this widely used form of the radial pulsation equation is not dimensionless but, rather, has units of inverse length-squared, we have found it useful to also recast it in the following dimensionless form:

<math> \frac{d^2x}{d\chi_0^2} + \biggl[\frac{4}{\chi_0} - \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} \biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \biggr] \frac{dx}{d\chi_0} + \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} \biggl(\frac{1}{\gamma_\mathrm{g}} \biggr)\biggl[\tau_\mathrm{SSC}^2 \omega^2 + (4 - 3\gamma_\mathrm{g})\biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \frac{1}{\chi_0} \biggr] x = 0 , </math>

where,

<math>~g_\mathrm{SSC} \equiv \frac{P_c}{R\rho_c} \, ,</math>       and       <math>~\tau_\mathrm{SSC} \equiv \biggl[\frac{R^2 \rho_c}{P_c}\biggr]^{1/2} \, .</math>

In a separate discussion, we showed that specifically for isolated, polytropic configurations, this linear adiabatic wave equation (LAWE) can be rewritten as,

<math>~0 </math>

<math>~=</math>

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

 

<math>~=</math>

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

where we have adopted the dimensionless frequency notation,

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

<math>~\equiv</math>

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

In this chapter we carry out a numerical integration of this governing LAWE for <math>~n=3</math> polytropes. The results are presented below.

Schwarzschild (1941)

We can directly compare our results with Schwarzschild's (1941) published work on "Overtone Pulsations for the Standard [Stellar] Model." To begin with, it is straightforward to demonstrate that the last form of the LAWE, provided above, matches equation (2) from Schwarzschild (1941), if <math>~n</math> is set to 3 — see the boxed-in excerpt, immediately below. Note as well that Schwarzschild's dimensionless oscillation frequency — defined in his equation (1) and which we will label, <math>~\omega_\mathrm{Sch}</math> — is related to our dimensionless frequency via the expression,

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

<math>~~\leftrightarrow~~</math>

<math>~\biggl( \frac{3\gamma_g}{2} \biggr) \omega_\mathrm{Sch}^2 \, .</math>

Schwarzschild (1941) numerically integrated the LAWE for <math>~n=3</math> polytropic spheres to find eigenvectors (i.e., the spatially discrete eigenfunction and corresponding eigenfrequency) for five separate oscillation modes (the fundamental mode, plus the 1st, 2nd, 3rd, and 4th overtones) for models having four different adopted adiabatic indexes <math>~\gamma_g = \tfrac{4}{3}, \tfrac{10}{7}, \tfrac{20}{13}, \tfrac{5}{3})</math>.


Paragraph extracted from M. Schwarzschild (1941)

"Overtone Pulsations for the Standard Model"

ApJ, vol. 94, pp. 245 - 252 © American Astronomical Society

Schwarzschild (1941, ApJ, 94, 245)

      3A. S. Eddington (1930), The Internal Constitution of the Stars, pp. 188 and 192.

Drawing from our discussion of the historical treatment of boundary conditions, we presume that Schwarzschild imposed the following constraint at the surface:

<math>~\frac{d\ln x}{d\ln \xi}\biggr|_\mathrm{surface}</math>

<math>~=</math>

<math>~\frac{1}{\gamma_g} \biggl( 4 - 3\gamma_g + \frac{\omega^2 R^3}{GM_\mathrm{tot}}\biggr) </math>

 

<math>~=</math>

<math>~\biggl[ \frac{3\omega^2 R^3}{4\pi G \gamma \bar\rho} - \alpha \biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{2} \biggl[ \frac{\sigma_c^2}{\gamma } \biggl(\frac{\rho_c}{\bar\rho}\biggr) -2 \alpha \biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{2} \biggl\{ \biggl[ \mathfrak{F} + 2\alpha \biggr] \biggl(\frac{\rho_c}{\bar\rho}\biggr) -2 \alpha \biggr\} \, .</math>

Recognizing from an accompanying tabulation that, for <math>~n=3</math> polytropes,

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

<math>~\approx</math>

<math>~54.18248 \, ,</math>

we presume that the surface boundary condition imposed by Schwarzschild was,

<math>~\frac{d\ln x}{d\ln \xi}\biggr|_\mathrm{surface}</math>

<math>~=</math>

<math>~27.09124 ( \mathfrak{F} + 2\alpha ) - \alpha \, .</math>

Our Table 1 catalogs the eigenfrequencies that Schwarzschild determined (drawn from his Table 1) for these twenty different models/modes.

Table 1:   From Table 1 of M. Schwarzschild (1941)

Mode <math>~\alpha = 0.0</math>

<math>~(\gamma_g = 4/3)</math>
<math>~\alpha = 0.2</math>

<math>~(\gamma_g = 10/7)</math>
<math>~\alpha = 0.4</math>

<math>~(\gamma_g = 20/13)</math>
<math>~\alpha = 0.6</math>

<math>~(\gamma_g = 5/3)</math>
<math>~\omega_\mathrm{Sch}^2</math> <math>~\omega_\mathrm{Sch}^2</math> <math>~\omega_\mathrm{Sch}^2</math> <math>~\mathfrak{F} = \biggl[\frac{3\omega_\mathrm{Sch}^2}{2} - 2\alpha \biggr]</math> <math>~\omega_\mathrm{Sch}^2</math>
0 0.00000 0.05882 0.10391 -0.64414 0.13670
1 0.16643 0.19139 0.21998 -0.47003 0.25090
2 0.3392 0.3648 0.3920 -0.2120 0.4209
3 0.5600 0.5863 0.6136 +0.1204 0.6420
4 0.8283 0.8554 0.8832 +0.5248 0.9117

Schwarzschild (1941) also documented the radial structure of the eigenfunction that is associated with each of these twenty model/mode eigenfrequencies. Each column of his Table 4, except the first, presents numerical values of the amplitude of a specific model/mode at 84 discrete radial locations throughout the n = 3 polytrope; the first column of the table lists the corresponding radial coordinate, <math>~\xi</math>. Focusing on the model that he analyzed assuming <math>~\alpha = 0.4</math>, we have typed his five columns of data into an Excel spreadsheet and have used this data to generate the pair of plots displayed, below, in Figure 1. The left-hand panel displays the eigenfunction amplitude versus radius, <math>~x(\xi)</math>, for the fundamental mode as well as for the first four overtones; it essentially replicates Figure 1 from Schwarzschild (1941). The right-hand panel displays the same data, but as a semi-log plot; specifically, it displays <math>~y(\xi)</math>, where,

<math>~y \equiv \frac{1}{2} \log_{10}[x^2 + 10^{-8}] \, .</math>

Each sharp valley in this semi-log plot highlights the location of a node in the corresponding eigenfunction, that is, it identifies where <math>~x(\xi)</math> crosses through zero.


Figure 1: 

Schwarzschild's Eigenfunctions for an n = 3 Polytrope with α = 0.4 (γ = 20/13)

Schwarzschild (1941) eigenfunctions

Numerical Integration

From the Core to the Surface

Here we use the finite-difference algorithm described separately to integrate the discretized LAWE from the center of the polytropic configuration, outward to its surface, which in this case — see, for example, p. 77 of Horedt (2004) — is located at the polytropic-coordinate location,

<math>~\xi_\mathrm{max} = 6.89684862 \, .</math>

It is assumed, at the outset, that we have in hand an appropriately discretized description of the unperturbed, equilibrium properties of an <math>~n=3</math> polytrope; specifically, at each radial grid line, we have tabulated values of the radial coordinate, <math>~0 \le \xi_i \le \xi_\mathrm{max}</math>, the Lane-Emden function, <math>~\theta_i</math>, and its first radial derivative, <math>~\theta_i'</math>.

The algorithm is as follows (substitute <math>~n=3</math> everywhere):

  • Establish an equally spaced radial-coordinate grid containing <math>~N</math> grid zones (and, accordingly, <math>~N+1</math> grid lines), in which case the grid-spacing parameter, <math>~\Delta_\xi \equiv \xi_\mathrm{max}/N</math>.
  • Specify a value of the adiabatic exponent, <math>~\gamma</math>, which, in turn, determines the value of the parameter, <math>~\alpha \equiv (3-4/\gamma) \, .</math>
  • Choose a value for the (square of the) dimensionless oscillation frequency, <math>~\sigma_c^2</math>, which we will accomplish by assigning a value to the parameter,

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

  • Set the eigenfunction to unity at the center <math>~(\xi_0 = 0)</math> of the configuration, that is, set <math>~x_0 = 1</math>.
  • Determine the value of the eigenfunction at the first grid line away from the center — having coordinate location, <math>~\xi_1 = \Delta_\xi </math> — via the derived power-series expression,

    <math>~ x_1 </math>

    <math>~=</math>

    <math>~ x_0 \biggl[ 1 - \frac{\Delta_\xi^2 (n+1) \mathfrak{F}}{60} \biggr] \, .</math>

  • At all other grid lines, <math>~i=2,N</math>, determine the value of the eigenfunction, <math>~x_i</math>, via the expression,

    <math>~x_i \biggl[2\theta_{i-1} +\frac{4\Delta_\xi \theta_{i-1}}{\xi_{i-1}} - \Delta_\xi (n+1)(- \theta^')_{i-1}\biggr] </math>

    <math>~=</math>

    <math>~ x_{i-1}\biggl\{4\theta_{i-1} - \frac{\Delta_\xi^2(n+1)}{3}\biggl[ \mathfrak{F}+2\alpha - 2\alpha \biggl(- \frac{3\theta^'}{\xi}\biggr)_{i-1} \biggr] \biggr\} + x_{i-2} \biggl[\frac{4\Delta_\xi \theta_{i-1}}{\xi_{i-1}} - \Delta_\xi (n+1)(- \theta^')_{i-1} - 2\theta_{i-1}\biggr] \, .</math>

We divided our model into <math>~N = 200</math> radial zones and, using this algorithm, integrated the LAWE from the center of the configuration to the surface, for <math>~\alpha = 0.4</math>, and approximately 40 different chosen values of the frequency parameter across the range, <math>~-0.7 \le \mathfrak{F} \le + 0.3</math>. The radial displacement functions resulting from these integrations are presented in Figure 2 as an animation sequence. The specified value of <math>~\mathfrak{F}</math> is displayed at the top of each animation frame, and the resulting displacement function, <math>~x(r/R)</math>, is traced by the small, red circular markers in each frame.

Figure 2:  Numerically Determined Eigenfunctions for Various <math>~\mathfrak{F}</math>
Table 2
Mode Match Schwarzschild Match B.C.
<math>~\mathfrak{F}</math> <math>~\mathfrak{F}</math>
0 -0.644131578154 -0.644131577959
1 -0.47013976423 -0.47013975308
2 -0.2121284391 -0.2121282667
3 +0.1202565375 +0.120257856

Eigenfunctions for Standard Model

Each frame of the Figure 2 animation also displays, as smooth solid curves, the radial eigenfunctions that Schwarzschild (1941) obtained for the fundamental mode (blue curve) and the first three overtone modes (green, purple, & orange curves, repectively) for his model with <math>~\alpha = 0.4</math>. These are the same curves that appear in the left-hand panel of Figure 1, but here the displacement amplitude has been renormalized such that <math>~x_0 = 1</math>, and, along the horizontal axis, the radial location is marked in terms of the fractional radius, <math>~r/R \equiv \xi/\xi_\mathrm{max}</math>. In our examination of this model, as we approached each specific value of a modal eigenfrequency identified by Schwarzschild — see the frequencies highlighted in pink in our Table 1 — we fine-tuned our choice of the eigenfrequency in order to find a displacement function whose surface amplitude matched, to a high level of precision, the surface amplitude associated with Schwarzschild's corresponding published eigenfunction. The column of our Table 2 whose heading is "Match Schwarzschild" identifies — to at least 10 digits precision — the frequency choice that was required in order for these surface amplitudes to match in each case.

It is gratifying to see that our resulting frequencies match well the values published by Schwarzschild (as highlighted in pink, above). But this does not satisfactorily explain why, among the entire range of displacement functions displayed (in red) in the Figure 2 animation, Schwarzschild labeled these specific ones as the eigenmodes. As we shall now demonstrate, his eigenmode identifications resulted from the imposition of a specific, physically justified constraint on the slope, rather than the value, of the displacement function at the surface of the configuration. (See also our separate brief answer to the question, "What makes this an eigenvalue problem?".)

Surface Boundary Condition

As was stated, above, we presume that as Schwarzschild searched for natural modes of oscillation in isolated, <math>~n=3</math> polytropes, he imposed the following boundary condition at the surface of the configuration:

<math>~\frac{d\ln x}{d\ln \xi}\biggr|_\mathrm{surface}</math>

<math>~=</math>

<math>~27.09124 ( \mathfrak{F} + 2\alpha ) - \alpha \, .</math>

In order to duplicate his findings, then, we need to fine tune our specification of the oscillation frequency such that the resulting displacement function presents this behavior at the surface of our model. A finite-difference expression of this logarithmic derivative that is consistent with the above-described finite-difference algorithm, is,

<math>~\frac{d\ln x}{d\ln \xi} \biggr|_\mathrm{surface}</math>

<math>~\approx</math>

<math>~\frac{\xi_\mathrm{max}}{x_N} \biggl[ \frac{x_{N+1}-x_{N-1}}{2\Delta_\xi} \biggr] \, .</math>

Everything is known here, except for the quantity, <math>~x_{N+1}</math>, which can be evaluated using the last expression in our algorithm one more time to, in effect, evaluate the eigenfunction just outside the surface. That is, we obtain <math>~x_{N+1}</math> and, in turn, obtain a value for the logarithmic derivative at the surface, via the expression,

<math>~x_{N+1} \biggl[2\theta_{N} +\frac{4\Delta_\xi \theta_{N}}{\xi_\mathrm{max}} - \Delta_\xi (n+1)(- \theta^')_{N}\biggr] </math>

<math>~=</math>

<math>~ x_{N}\biggl\{4\theta_{N} - \frac{\Delta_\xi^2(n+1)}{3}\biggl[ \mathfrak{F}+2\alpha - 2\alpha \biggl(- \frac{3\theta^'}{\xi}\biggr)_{N} \biggr] \biggr\} + x_{N-1} \biggl[\frac{4\Delta_\xi \theta_{N}}{\xi_\mathrm{max}} - \Delta_\xi (n+1)(- \theta^')_{N} - 2\theta_{N}\biggr] \, .</math>

We added to our numerical algorithm a step that evaluates, in this manner, the logarithmic derivative of the displacement function at the surface of our polytropic configuration. The eigenfrequencies that generated displacement functions with this surface behavior are listed for four separate modes in the column of Table 2 titled, "Match B.C." In every case the values agree to at least five decimal places with the "Match Schwarzschild" eigenfrequencies. We conclude, therefore, that it was the implementation of this surface boundary condition that permitted Schwarzschild to quantitatively identify the properties of the eigenvectors associated with natural radial modes of oscillation in <math>~n=3</math> polytropes.

Our Results

Table 3:   Our Results (to be compared w/ Table 1, above)

Mode <math>~\alpha = 0.0</math>

<math>~(\gamma_g = 4/3)</math>
<math>~\alpha = 0.2</math>

<math>~(\gamma_g = 10/7)</math>
<math>~\alpha = 0.4</math>

<math>~(\gamma_g = 20/13)</math>
<math>~\alpha = 0.6</math>

<math>~(\gamma_g = 5/3)</math>
<math>~\mathfrak{F}</math> <math>~\omega_\mathrm{Sch}^2 = \frac{2}{3}\biggl(\mathfrak{F}+2\alpha \biggr)</math> <math>~\mathfrak{F}</math> <math>~\omega_\mathrm{Sch}^2</math> <math>~\mathfrak{F} </math> <math>~\omega_\mathrm{Sch}^2</math> <math>~\mathfrak{F}</math> <math>~\omega_\mathrm{Sch}^2</math>
0 --- --- -0.311782342981 0.058812 -0.644131577959 0.103912 --- ---
1 +0.24946512002 0.166310 -0.113086698932 0.191276 -0.47013975308 0.219907 --- ---
2 +0.50882623652 0.339217 +0.14705874055 0.364706 -0.2121282667 0.391914 --- ---
3 +0.83977118 0.559847 +0.479241829 0.586161 +0.120257856 0.613505 --- ---
4 +1.24253191 0.828355 +0.8832297 0.855486 +0.52498863 0.883326 --- ---

Truncated n = 3 Polytropes

We understand that, for isolated <math>~n=3</math> polytropic spheres, the value of the adiabatic exponent for which the configuration is marginally unstable is <math>~\gamma_g = 4/3</math>, which is equivalent to, <math>~\alpha \equiv (3-4/\gamma_g) = 0</math>. This critical condition is identified by examining when the oscillation frequency of the fundamental mode goes to zero. Let's use our numerical integration tool to determine what this critical value of the adiabatic exponent is for truncated, <math>~n=3</math> polytropes. We will accomplish this as follows:

  • At various truncation radii, <math>~0 < \xi_\mathrm{surf}/\xi_\mathrm{max} < 1</math> …
    • Force <math>~\mathfrak{F} = -2\alpha</math>;
    • Iterate on the choice of <math>~\alpha</math> until the displacement function with no radial nodes (i.e., the fundamental mode) satisfies the surface boundary condition of <math>~(d\ln x/d\ln \xi)_\mathrm{surf} = -3</math>, to a desired level of accuracy.

The following table shows the values of <math>~\alpha_\mathrm{crit}</math> — and associated values of <math>~\gamma_\mathrm{crit}</math> — that we obtained for nine different values of <math>~\xi_\mathrm{surf}/\xi_\mathrm{max}</math>; in each case, iterations were continued until the desired surface boundary condition was satisfied to six significant digits.

Pressure-Truncated n = 3 Polytropes
<math>~N_\mathrm{zones} = 200</math>
Edge Zone <math>~\frac{\xi_\mathrm{surf}}{\xi_\mathrm{max}}</math> <math>~\alpha_\mathrm{crit}</math> <math>~\mathfrak{F}</math> Surface B.C <math>~\gamma_\mathrm{crit}</math>
200 1.00 0 0 --- <math>~\tfrac{4}{3}</math>
180 0.90 -0.000201541 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 1.333244
160 0.8 -0.00327575 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 1.331879
150 0.75 -0.00808603 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 1.329749
120 0.60 -0.0576031 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 1.308214
100 0.5 -0.159111 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 1.266179
80 0.40 -0.405712 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 1.174497
50 0.25 -1.74909 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 0.842266
20 0.10 -14.6648 <math>~-2\alpha_\mathrm{crit}</math> -3.00000 0.226439

Analytic Inquiry

NOTE (from J. E. Tohline in April, 2017):   The following subsections present some exploratory ideas that were pursued while I was searching for analytic solutions to the polytropic LAWE. For the most part this material has been superseded by a separate discussion in which we describe the desired analytic solution, which we discovered in March, 2017.

Fundamental-Mode, Homentropic Oscillations

The LAWE, presented above, that is relevant to polytropic spheres, may be rewritten as,

LSU Key.png

<math>~0 = \frac{d^2x}{d\xi^2} + \biggl[ 4 - (n+1) Q \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} + (n+1) \biggl[ \biggl( \frac{\sigma_c^2}{6\gamma_g } \biggr) \frac{\xi^2}{\theta} - \alpha Q\biggr] \frac{x}{\xi^2} </math>

where:    <math>~Q(\xi) \equiv - \frac{d\ln\theta}{d\ln\xi} \, ,</math>    <math>~\sigma_c^2 \equiv \frac{3\omega^2}{2\pi G\rho_c} \, ,</math>     and,     <math>~\alpha \equiv \biggl(3 - \frac{4}{\gamma_\mathrm{g}}\biggr)</math>

Now, if we assume that oscillations occur adiabatically with an adiabatic index that is consistent with the chosen polytropic index — that is to say,

<math>~\gamma_g = \frac{n+1}{n}</math>      <math>~\Rightarrow</math>       <math>~\alpha = \frac{3-n}{n+1} \, ,</math>

in which case the configuration remains homentropic as it oscillates — and if we look only for a (marginally unstable) configuration that has <math>~\sigma_c^2 = 0</math>, then the relevant LAWE is,

<math>~0 </math>

<math>~=</math>

<math>~ \frac{d^2x}{d\xi^2} + \biggl[ 4 - (n+1) Q \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} - (3-n) Q \frac{x}{\xi^2} \, . </math>

Specific case of n = 3 Polytropes

Homologous Collapse

If we examine only an <math>~n=3</math> polytropic configuration, then the last term disappears. This means that in this very special case, a perfectly valid solution to the LAWE is <math>~x = \mathrm{constant}</math>. This is presumably the eigenfunction that Schwarzschild deduced; the fundamental-mode "oscillations" are perfectly homologous. Given that the model is marginally unstable, an ensuing dynamical collapse will presumably begin in a perfectly homologous fashion. This is precisely the type of "free-fall" collapse that was discussed and modeled by Goldreich & Weber (1980).

Another Potential Option

We have wondered whether, in this very special case, one or more additional fundamental-mode eigenfunction(s) might satisfy the governing LAWE. Here is a relevant line of arguments, beginning with the LAWE for the n = 3 polytropic sphere.

<math>~0 </math>

<math>~=</math>

<math>~ \frac{d^2x}{d\xi^2} + 4( 1 + Q ) \frac{1}{\xi} \cdot \frac{dx}{d\xi} </math>

<math>~\Rightarrow ~~~ \frac{dy}{d\xi} </math>

<math>~=</math>

<math>~ - 4( 1 + Q ) \frac{y}{\xi} </math>

<math>~\Rightarrow ~~~ \frac{1}{4}\cdot \frac{d\ln y}{d\ln\xi} </math>

<math>~=</math>

<math>~ - ( 1 + Q ) \, , </math>

where,

<math>~y \equiv \frac{dx}{d\xi} \, .</math>

But, by definition, the function <math>~Q(\xi)</math> is a logarithmic derivative of the Lane-Emden function. Hence, we also can write,

<math>~-1 </math>

<math>~=</math>

<math>~ \frac{1}{4}\cdot \frac{d\ln y}{d\ln\xi} + \frac{d\ln\theta}{d\ln\xi}</math>

 

<math>~=</math>

<math>~ \frac{d\ln (\theta y^{1/4})}{d\ln\xi} </math>

<math>~\Rightarrow~~~- d\ln\xi</math>

<math>~=</math>

<math>~ d\ln (\theta y^{1/4}) \, . </math>

Integrating this equation once gives,

<math>~\ln(\theta y^{1/4}) + \ln\xi</math>

<math>~=</math>

<math>~\ln (c_0)</math>

<math>~\Rightarrow ~~~ \xi \theta y^{1/4} </math>

<math>~=</math>

<math>~c_0</math>

<math>~\Rightarrow ~~~\frac{dx}{d\xi} </math>

<math>~=</math>

<math>~\biggl( \frac{c_0}{\xi\theta}\biggr)^4 \, .</math>

Referring to the power-series expansion of the polytropic Lane-Emden function, <math>~\theta(\xi)</math>, about the configuration's center, we see that the product, <math>~\xi\theta</math>, goes to zero as the first power of <math>~\xi</math>. This means that the right-hand side of this last differential equation blows up at the center. This, therefore, does not appear to provide a physically viable avenue by which to identify an alternative fundamental-mode eigenfunction.

Play With Form of LAWE

Logarithmic Derivative Rewrite

We have noticed that the LAWE that governs the eigenfunction associated with the fundamental mode of the marginally unstable model (FMMUM),

<math>~0 </math>

<math>~=</math>

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

may be rewritten entirely as an expression that relates the logarithmic derivatives of <math>~x, \xi,</math> and <math>~\theta</math>. Multiplying through by <math>~\xi^2/x</math>, then drawing on a differential relation that has been derived in a separate context, namely,

<math>~\frac{\xi^2}{x} \cdot \frac{d^2x}{d\xi^2} </math>

<math>~=</math>

<math>~ \frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr] + \biggl[ \frac{d\ln x}{d\ln \xi} -1 \biggr]\cdot \frac{d\ln x}{d\ln \xi} \, , </math>

this LAWE associated with the FMMUM becomes,

<math>~0 </math>

<math>~=</math>

<math>~ \frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr] + \biggl[ \frac{d\ln x}{d\ln \xi} -1 \biggr]\cdot \frac{d\ln x}{d\ln \xi} + \biggl[ 4 - (n+1) Q \biggr] \frac{d\ln x}{d\ln\xi} - (3-n) Q </math>

 

<math>~=</math>

<math>~ \frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr] + \frac{d\ln x}{d\ln \xi} \cdot \biggl[ \frac{d\ln x}{d\ln \xi} + 3 - (n+1) Q \biggr] - (3-n) Q </math>

 

<math>~=</math>

<math>~ \frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr] + \frac{d\ln x}{d\ln \xi} \cdot \biggl\{ \frac{d\ln [x \xi^3 \theta^{(n+1)}] }{d\ln \xi} \biggr\} + \frac{d\ln \theta^{(3-n) }}{d\ln\xi} \, . </math>

I'm not sure if anyone else has previously appreciated that the "fundamental mode" polytropic LAWE can be written in this form. I'm even less sure that this form sheds light on its solution.

Play a little more …   Start by letting, <math>~A \equiv [x \xi^3 \theta^{(n+1)}] \, ,</math> in which case we have,

<math>~0 </math>

<math>~=</math>

<math>~ \frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr] + \frac{d\ln x}{d\ln \xi} \cdot \biggl\{ \frac{1}{A} \cdot \frac{dA}{d\ln \xi} \biggr\} + \frac{d\ln \theta^{(3-n) }}{d\ln\xi} </math>

 

<math>~=</math>

<math>~ \frac{1}{A} \cdot \frac{d}{d\ln\xi} \biggl[ A \cdot \frac{d\ln x}{d\ln \xi} \biggr] + \frac{d\ln \theta^{(3-n) }}{d\ln\xi} </math>

<math>~\Rightarrow~~~ \frac{d}{d\xi} \biggl[ A \cdot \frac{d\ln x}{d\ln \xi} \biggr] </math>

<math>~=</math>

<math>~ - \frac{A}{\theta^{(3-n)}} \cdot \frac{d\theta^{(3-n) }}{d\xi} </math>

<math>~\Rightarrow~~~ \frac{d}{d\xi} \biggl[ x \xi^3 \theta^{(n+1)} \cdot \frac{d\ln x}{d\ln \xi} \biggr] </math>

<math>~=</math>

<math>~ - \frac{x \xi^3 \theta^{(n+1)}}{\theta^{(3-n)}} \cdot \frac{d\theta^{(3-n) }}{d\xi} </math>

<math>~\Rightarrow~~~ \frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \cdot \frac{dx}{d\xi} \biggr] </math>

<math>~=</math>

<math>~ - x \xi^3 \theta^{2(n-1)} \cdot \frac{d\theta^{(3-n) }}{d\xi} \, . </math>

One Feeble Guess

Now, what if, <math>~x \equiv [\xi^{-2}\theta^{-n}]</math>   ?

<math>~\Rightarrow~~~ \xi^2 \theta^{(n+1)} \cdot \frac{d(\xi^{-2}\theta^{-n})}{d\xi} </math>

<math>~=</math>

<math>~ \frac{d\theta}{d\xi} - \frac{1}{\xi^{2}\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi} \, , </math>

in which case, the LAWE becomes,


<math>~ \frac{d}{d\xi} \biggl[ \xi^2 \frac{d\theta}{d\xi} - \frac{1}{\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi}\biggr] </math>

<math>~=</math>

<math>~ - [\xi^{-2}\theta^{-n}] \xi^3 \theta^{2(n-1)} \cdot \frac{d\theta^{(3-n) }}{d\xi} </math>

<math>~\Rightarrow~~~ \frac{d}{d\xi} \biggl[ \xi^2 \frac{d\theta}{d\xi} \biggr] - \frac{d}{d\xi} \biggl[\frac{1}{\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi}\biggr] </math>

<math>~=</math>

<math>~ - \xi \theta^{n-2} \cdot \frac{d\theta^{(3-n) }}{d\xi} </math>

<math>~\Rightarrow~~~ -\xi^2\theta^n - \frac{d}{d\xi} \biggl[\frac{1}{\theta^{n}} \cdot \frac{d(\xi^2 \theta^{(n+1)} )}{d\xi}\biggr] </math>

<math>~=</math>

<math>~ - \xi (3-n)\frac{d\theta }{d\xi} </math>

Behavior for Known n=5 Solution

We know that the FMMUM for pressure-truncated, n = 5 polytropic configurations takes the form,

<math>~x </math>

<math>~=</math>

<math>~1- f_n(\xi) \, ,</math>

where,

<math>~f_5(\xi) = \frac{\xi^2}{15} \, .</math>

The governing LAWE therefore gives,

<math>~ \frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \cdot \frac{dx}{d\xi} \biggr] </math>

<math>~=</math>

<math>~ - x \xi^3 \theta^{2(n-1)} \cdot \frac{d\theta^{(3-n) }}{d\xi} </math>

 

<math>~=</math>

<math>~ - (3-n)x \xi^3 \theta^{n} \cdot \frac{d\theta }{d\xi} </math>

<math>~\Rightarrow~~~ (3-n)(1-f_n) \xi^3 \theta^{n} \cdot \frac{d\theta }{d\xi} </math>

<math>~=</math>

<math>~ \frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \cdot \frac{df_n}{d\xi} \biggr] </math>

 

<math>~=</math>

<math>~ \xi^4 \theta^{(n+1)} \cdot\frac{d^2f_n}{d\xi^2} + \frac{df_n}{d\xi} \cdot \frac{d}{d\xi} \biggl[ \xi^4 \theta^{(n+1)} \biggr] </math>

 

<math>~=</math>

<math>~ \xi^4 \theta^{(n+1)} \cdot\frac{d^2f_n}{d\xi^2} + \frac{df_n}{d\xi} \cdot \biggl\{ 4\xi^3 \theta^{(n+1)} + (n+1)\xi^4 \theta^n \frac{d\theta}{d\xi} \biggr\} </math>

<math>~\Rightarrow~~~ (3-n)(1-f_n) \frac{d\theta }{d\xi} </math>

<math>~=</math>

<math>~ \xi \theta \cdot\frac{d^2f_n}{d\xi^2} + \frac{df_n}{d\xi} \cdot \biggl\{ 4 \theta + (n+1)\xi \frac{d\theta}{d\xi} \biggr\} </math>

<math>~\Rightarrow~~~ \frac{1}{\theta} \frac{d\theta }{d\xi} \biggl[ (3-n)(1-f_n) - (n+1)\xi \frac{df_n}{d\xi} \biggr] </math>

<math>~=</math>

<math>~ \xi \cdot\frac{d^2f_n}{d\xi^2} + 4 \cdot \frac{df_n}{d\xi} \, . </math>

Let's check to see whether the known <math>~f_5(\xi)</math> function properly satisfies this last ODE when n = 5.

<math>~ -\xi (3+\xi^2)^{-1} \biggl[ -2 \biggl(1-\frac{\xi^2}{15} \biggr) - \frac{12\xi^2}{15} \biggr] </math>

<math>~=</math>

<math>~ \frac{2\xi}{15} + \frac{8\xi}{15} </math>

<math>~\Rightarrow~~~ 2 + \frac{2\xi^2}{3} </math>

<math>~=</math>

<math>~ \frac{2}{3} (3+\xi^2) \, , </math>      Yes!

Given that,

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

<math>~=</math>

<math>~\frac{3}{3+\xi^2}</math>

<math>~\Rightarrow ~~~\xi^2 </math>

<math>~=</math>

<math>~3\biggl[ \frac{1}{\theta_5^2} - 1 \biggr] \, ,</math>

we could presume that, when defined in terms of <math>~\theta_5</math>, the defining function,

<math>~f_5(\theta_5) = \frac{1}{n}\biggl[ \frac{1}{\theta_5^2} - 1 \biggr]</math>

   <math>~\Rightarrow~</math>   

<math>~ \frac{df_5}{d\xi} = - \frac{2}{n\theta_5^3} \cdot \frac{d\theta_5}{d\xi} </math>

   <math>~\Rightarrow~</math>   

<math>~ \frac{d^2f_5}{d\xi^2} = - \frac{2}{n\theta_5^3} \cdot \frac{d^2\theta_5}{d\xi^2} + \frac{6}{n\theta_5^4} \cdot \biggl( \frac{d\theta_5}{d\xi} \biggr)^2 \, . </math>

In this case, the governing LAWE becomes,

<math>~0 </math>

<math>~=</math>

<math>~ \xi \cdot\frac{d^2f_n}{d\xi^2} - \frac{(3-n)}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ 1-f_n \biggr] + \biggl[ 4 + \frac{(n+1)\xi}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggr] \biggl[ \frac{df_n}{d\xi} \biggr] </math>

 

<math>~=</math>

<math>~ \xi \biggl\{ - \frac{2}{n\theta_5^3} \cdot \frac{d^2\theta_5}{d\xi^2} + \frac{6}{n\theta_5^4} \cdot \biggl( \frac{d\theta_5}{d\xi} \biggr)^2\biggr\} - \biggl[ 4 + \frac{(n+1)\xi}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggr] \biggl[ \frac{2}{n\theta_5^3} \cdot \frac{d\theta_5}{d\xi}\biggr] - \frac{(3-n)}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ \frac{(n + 1)\theta_5^2 -1 }{n\theta_5^2} \biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{n\theta_5^3}\biggl\{ \biggl[ 2\xi \cdot \frac{d^2\theta_5}{d\xi^2} - \frac{6\xi}{\theta_5} \cdot \biggl( \frac{d\theta_5}{d\xi} \biggr)^2\biggr] + 2\biggl[ 4 + \frac{(n+1)\xi}{\theta} \biggl( \frac{d\theta }{d\xi} \biggr) \biggr] \frac{d\theta_5}{d\xi} + (3-n) \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ (n + 1)\theta_5^2 -1 \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~\frac{1}{n\theta_5^3}\biggl\{ 2\xi \cdot \frac{d^2\theta_5}{d\xi^2} + (n-2) \frac{2\xi}{\theta_5}\biggl( \frac{d\theta_5}{d\xi} \biggr)^2 + \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ (3-n) (n + 1)\theta_5^2 + 5+ n \biggr] \biggr\} \, . </math>

Now, from the polytropic Lane-Emden equation, we also know that,

LSU Key.png

<math>~\frac{1}{\xi^2} \frac{d}{d\xi}\biggl( \xi^2 \frac{d\Theta_H}{d\xi} \biggr) = - \Theta_H^n</math>

That is,

<math>~2\xi \biggl(\frac{d\theta_5}{d\xi}\biggr) + \xi^2 \cdot \frac{d^2\theta_5}{d\xi^2} </math>

<math>~=</math>

<math>~- \xi^2 \theta_5^n</math>

<math>~\Rightarrow~~~ \frac{d^2\theta_5}{d\xi^2} </math>

<math>~=</math>

<math>~- \theta_5^n - \frac{2}{\xi} \biggl(\frac{d\theta_5}{d\xi}\biggr)</math>

So, the LAWE becomes,

<math>~0</math>

<math>~=</math>

<math>~ (n-2) \frac{2\xi}{\theta_5}\biggl( \frac{d\theta_5}{d\xi} \biggr)^2 + (n+1)\biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ (3-n) \theta_5^2 + 1 \biggr] - 2\xi \theta_5^n \, . </math>

Again, let's check to see if the case of n=5 works …

<math>~0</math>

<math>~=</math>

<math>~ \biggl( \frac{d\theta_5}{d\xi} \biggr)^2 + \frac{\theta_5}{\xi} \biggl( \frac{d\theta }{d\xi} \biggr) \biggl[ 1 - 2 \theta_5^2 \biggr] - \frac{\theta_5^{6}}{3} </math>

 

<math>~=</math>

<math>~ 3\xi^2(3+\xi^2)^{-3} -3(3+\xi^2)^{-2}\biggl[ 1 - 6(3+\xi^2)^{-1} \biggr] - 3^2(3+\xi^2)^{-3} </math>

 

<math>~=</math>

<math>~(3+\xi^2)^{-3}\biggl\{ 3\xi^2 -3\biggl[ (3+\xi^2) -6 \biggr] - 3^2 \biggr\} </math>     Yes!

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