Editing SSC/Structure/BiPolytropes/Analytic15/Pt3 (section)

=Examples=
<table border="1" align="center" width="100%" colspan="8">
<tr>
  <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Structure/BiPolytropes/Analytic15|Part I: &nbsp; Steps 2 thru 7]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Structure/BiPolytropes/Analytic15/Pt2|Part II:&nbsp; Analytic Solution of Interface Relation]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Structure/BiPolytropes/Analytic15/Pt3|III:&nbsp; Modeling]]<br />&nbsp;</td>
  <td align="center" bgcolor="lightblue"><br />[[SSC/Structure/BiPolytropes/MurphyUVplane|IV:&nbsp; Murphy's UV Plane]]<br />&nbsp;</td>
</tr>
</table>

==Normalization==
The dimensionless variables used in Tables 1 &amp; 2 are defined as follows:
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>~\rho^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{\rho}{\rho_0}</math>
  </td>

  <td align="center">; &nbsp;&nbsp;&nbsp;</td>

  <td align="right">
<math>~r^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{r}{(K_c/G)^{1/2}}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~P^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{P}{K_c\rho_0^{2}}</math>
  </td>

  <td align="center">; &nbsp;&nbsp;&nbsp;</td>

  <td align="right">
<math>~M_r^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{M_r}{\rho_0 (K_c/G)^{3/2}}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~H^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{H}{K_c\rho_0}</math>
  </td>

  <td align="center">. &nbsp;&nbsp;&nbsp;</td>

  <td align="right" colspan="3">
&nbsp;
  </td>
</tr>

</table>
</div>

==Parameter Values==

The <math>2^\mathrm{nd}</math> column of Table 1 catalogues the analytic expressions that define various parameters and physical properties (as identified, respectively, in column 1) of the <math>~n_c=1</math>, <math>~n_e=5</math> bipolytrope.  


<!--  BEGIN TABLE OF PARAMETERS --->
<div align="center">
<table border="1" cellpadding="5">
<tr>
  <td align="center" colspan="2">
'''Properties of <math>~n_c=1</math>, <math>~n_e=5</math>, BiPolytrope Having Various Interface Locations, <math>~\xi_i</math>'''
  </td>
</tr>

<tr>
  <td align="center">
Parameter
  </td>
  <td align="center">
<math>~\xi_i</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\theta_i</math>
  </td>
  <td align="center">
<math>~\frac{\sin\xi_i}{\xi_i} </math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~-\biggl(\frac{d\theta_i}{d\xi}\biggr)_i</math>
  </td>
  <td align="center">
<math>~\frac{1}{\xi_i^2}(\sin\xi - \xi_i \cos\xi_i) </math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~r^*_\mathrm{core} \equiv r^*_i</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{1}{2\pi} \biggr)^{1/2} \xi_i</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\rho^*_i \biggr|_c = \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \rho^*_i \biggr|_e</math>
  </td>
  <td align="center">
<math>~\theta_i</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>P^*_i</math>
  </td>
  <td align="center">
<math>~\theta_i^2</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>H^*_i \biggr|_c = \frac{n_c+1}{n_e+1} \biggl( \frac{\mu_e}{\mu_c} \biggr) H^*_i \biggr|_e</math>
  </td>
  <td align="center">
<math>~2\theta_i</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>M^*_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{2}{\pi} \biggr)^{1/2} (\sin\xi_i - \xi_i \cos\xi_i)</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1}\eta_i</math>
  </td>
  <td align="center">
<math>\frac{\xi_i}{\sqrt{3}}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~ -\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>
  </td>
  <td align="center">
<math>~\frac{1}{\sqrt{3} \theta_i}\biggl( - \frac{d\theta}{d\xi} \biggr)_i = \frac{1}{\sqrt{3}}\biggl( \frac{1}{\xi_i}- \cot\xi_i\biggr)</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~p</math>
  </td>
  <td align="center">
<math>~\frac{3}{ 3-2(\mu_e/\mu_c)(1-\xi_i \cot\xi_i)}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~~y_\mathrm{root}</math>
  </td>
  <td align="center">
<math>~p\biggl( 1 + \{1 + [1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} + \{1 - [1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} \biggr)</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\Delta_i - m\pi</math>
  </td>
  <td align="center">
<math>~\tan^{-1}(y_\mathrm{root})</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~(A_0\eta)_\mathrm{root}</math>
  </td>
  <td align="center">
<math>~e^{2\Delta_i}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\xi_s</math>
  </td>
  <td align="center">
<math>~\xi_ie^{2(\pi-\Delta_i)}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\biggl[ \frac{A_0}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) \biggr]</math>
  </td>
  <td align="center">
<math>~\xi_i^{-1} e^{2\Delta_i}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\biggl[\frac{B_0}{3^{1/4}} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1/2} \biggr]</math>
  </td>
  <td align="center">
<math>~\biggl[ \xi_i \biggl( \frac{3}{\sin^2\Delta_i} - 2\biggr) \biggr]^{-1/2}
= \theta_i \biggl[ \frac{\xi_i^{1/2}}{\sin\xi_i}\biggl( \frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1/2} \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~ -\biggl( \frac{d\phi}{d\eta} \biggr)_s \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>
  </td>
  <td align="center">
<math>~\frac{1}{ 2\xi_i e^{3(\pi-\Delta_i)}} \biggl( \frac{3}{\sin^2\Delta_i} - 2 \biggr)^{1/2}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~ R^*_s</math>
  </td>
  <td align="center">
<math>~(2\pi)^{-1/2} \xi_i e^{2(\pi - \Delta_i)}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~ \biggl( \frac{\mu_e}{\mu_c}\biggr) M^*_\mathrm{tot}</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{3}{2\pi} \biggr)^{1/2} \sin\xi_i \biggl( \frac{3}{\sin^2\Delta_i} - 2 \biggr)^{1/2} e^{(\pi - \Delta_i)}</math>
  </td>
</tr>
</table>
</div>


<!-- OMIT THIS FIGURE
<div align="center">
<table border="1" cellpadding="0" align="center">
<tr>
  <td align="center">
Physically Viable Parameter Regime
  </td>
</tr>
<tr><td align="center">
[[File:MurphyPhysicallyViable.png|center|400px|Physically Viable Parameter Region]]
</td></tr>
</table>
</div>
STOP OMISSION -->

==Profile==

Once the values of the key set of parameters have been determined as illustrated in the preceding formula table, the radial profile of various physical variables can be determined throughout the bipolytrope as detailed in [[SSC/Structure/BiPolytropes/Analytic15#Step_4:__Throughout_the_core_.28.29|step &#35;4]] and [[SSC/Structure/BiPolytropes/Analytic15#Step_8:__Envelope.27s_Physical_Profile|step &#35;8]], above.  The following table summarizes the mathematical expressions that define the profile throughout the core (column 2) and throughout the envelope (column 3) of the normalized mass density, <math>~\rho^*(r^*)</math>, the normalized gas pressure, <math>~P^*(r^*)</math>, and the normalized mass interior to <math>~r^*</math>, <math>~M_r^*(r^*)</math>.  For all profiles, the relevant normalized radial coordinate is <math>~r^*</math>, as defined in the 2<sup>nd</sup> row of the table.  Graphical illustrations of these resulting profiles can be viewed by clicking on the thumbnail images posted in the last few columns of the table.  

<div align="center">
<b>Table 2:  Radial Profile of Various Physical Variables</b>
<table border="1" cellpadding="6">
<tr>
  <td align="center" rowspan="2">
Variable
  </td>
  <td align="center" rowspan="2">
Throughout the Core<br>
<math>0 \le \xi \le \xi_i</math>
  </td>
  <td align="center" rowspan="2">
Throughout the Envelope<sup>&dagger;</sup><br>
<math>\eta_i \le \eta \le \eta_s</math>
  </td>

  <td align="center" colspan="3">
Plotted Profiles
  </td>
</tr>

<tr>
  <td align="center">
<math>\xi_i = 0.5</math>
  </td>
  <td align="center">
<math>\xi_i = 1.0</math>
  </td>
  <td align="center">
<math>\xi_i = 3.0</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~r^*</math>
  </td>
  <td align="center">
<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \xi</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \biggl(\frac{3}{2\pi}\biggr)^{1/2}\eta</math>
  </td>

  <td align="center" colspan="3">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>~\rho^*</math>
  </td>
  <td align="center">
<math>\frac{\sin\xi}{\xi}</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i [\phi(\eta)]^5</math>
  </td>

  <td align="center">
<!--  [[File:PlotDensity_xi_0.5.jpg|thumb|75px]] -->
[[Image:DenXi05.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:DenXi10.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:DenXi30.jpg|thumb|75px]]
  </td>
</tr>

<tr>
  <td align="center">
<math>~P^*</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\sin\xi}{\xi} \biggr)^2</math>
  </td>
  <td align="center">
<math>\theta^{2}_i [\phi(\eta)]^{6}</math>
  </td>

  <td align="center">
<!-- [[File:PlotPressure_xi_0.5.jpg|thumb|75px]] -->
[[Image:PresXi05.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:PresXi10.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:PresXi30.jpg|thumb|75px]]
  </td>
</tr>

<tr>
  <td align="center">
<math>~M_r^*</math>
  </td>
  <td align="center">
<math>\biggl( \frac{2}{\pi}\biggr)^{1/2} (\sin\xi - \xi\cos\xi)</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \biggl( \frac{2\cdot 3^3 }{\pi} \biggr)^{1/2} \theta_i 
\biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr) </math>
  </td>

  <td align="center">
<!-- [[File:PlotPressure_xi_0.5.jpg|thumb|75px]] -->
[[Image:MassXi05.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:MassXi10.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:MassXi30.jpg|thumb|75px]]
  </td>
</tr>

<tr>
  <td align="left" colspan="6">
<sup>&dagger;</sup>In order to obtain the various envelope profiles, it is necessary to evaluate <math>~\phi(\eta)</math> and its first derivative using the information [[SSC/Structure/BiPolytropes/Analytic15#Step_6:__Envelope_Solution|presented in Step 6, above]].
  </td>
</tr>
</table>
</div>

==Murphy and Fiedler (1985)==


<div align="center">
<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" bgcolor="lightblue">
Table 1 from [http://adsabs.harvard.edu/abs/1985PASAu...6..219M Murphy &amp; Fiedler (1985, Proc. Astr. Soc. of Australia, 6, 219)]
  </td>
</tr>
<tr>
  <td align="center">
<!-- [[File:MF85Table1.png|center|600px|Murphy &amp; Fiedler (1985) Table 1]] -->
[[Image:AAAwaiting02.png|center|300px|Murphy &amp; Fiedler (1985) Table 1]]
  </td>
</tr>
<tr>
  <td align="center" bgcolor="lightblue">
Reproduction of Table 1 from MF85 Using Excel and Analytic Expressions Derived Here
  </td>
</tr>
<tr>
  <td align="center">
[[File:MF85Table1byTohline02.png|center|600px|Excel Regeneration of MF85 Table 1]]
  </td>
</tr>
</table>
</div>