Editing SSC/Structure/BiPolytropes/Analytic15

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

=BiPolytrope with n<sub>c</sub> = 1 and n<sub>e</sub> = 5=
{| class="PGEclass" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|- 
! style="height: 125px; width: 125px; background-color:white;" |
<font size="-1">[[H_BookTiledMenu#MoreModels|<b>Murphy (1983)<br /><br />Analytic</b>]]<br />(n<sub>c</sub>, n<sub>e</sub>) = (1, 5)</font>
|}
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 12 April 2015:  I became aware of the published discussions of this system by Murphy (1983) and Murphy &amp; Fiedler (1985b) in March of 2015 after searching the internet for previous analyses of radial oscillations in polytropes and, then, reading through Horedt's (2004) &sect;2.8.1 discussion of composite polytropes.]]Here we construct a [[SSC/Structure/BiPolytropes#BiPolytropes|system of bipolytropic configurations]] in which the core has an <math>~n_c=1</math> polytropic index and the envelope has an <math>~n_e=5</math> polytropic index.  As in the case of our [[SSC/Structure/BiPolytropes/Analytic51#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|separately discussed, "mirror image" bipolytropic configurations having <math>~(n_c, n_e) = (5, 1)</math>]], this system is particularly interesting because the entire structure can be described by closed-form, analytic expressions.  Bipolytropes of this type were first constructed by {{ Murphy83afull }}, and attributes of their physical structure were further discussed by {{ MF85afull }}; [[SSC/Structure/BiPolytropes/Analytic15#Key_References|additional, closely related references are given below]].  In the discussion that follows, we will be heavily referencing {{ Murphy83a }} &#8212; hereafter, {{ Murphy83ahereafter }}.
&nbsp;<br />
<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>

==Steps 2 &amp; 3==
Based on the discussion [[SSC/Structure/Polytropes/Analytic#n_=_1_Polytrope|presented elsewhere of the structure of an isolated n = 1 polytrope]], the core of this bipolytrope will have the following properties:
<div align="center">
<math>
\theta(\xi) = \frac{\sin\xi}{\xi} ~~~~\Rightarrow ~~~~  \theta_i = \frac{\sin\xi_i}{\xi_i} ;
</math>

<math>
\frac{d\theta}{d\xi} = \biggl[ \frac{\cos\xi}{\xi}- \frac{\sin\xi}{\xi^2}\biggr] ~~~~\Rightarrow ~~~~  
\biggl(\frac{d\theta}{d\xi}\biggr)_i = \biggl[\frac{\cos\xi_i}{\xi_i}- \frac{\sin\xi_i}{\xi_i^2} \biggr] \, .
</math>
</div>
The first zero of the function <math>~\theta(\xi)</math> and, hence, the surface of the corresponding isolated <math>~n=1</math> polytrope is located at <math>~\xi_s = \pi</math>.  Hence, the interface between the core and the envelope can be positioned anywhere within the range, <math>~0 < \xi_i < \pi</math>.

==Step 4:  Throughout the core (0 &le; &xi; &le; &xi;<sub>i</sub>)==

<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="center" colspan="3">
Specify:  <math>~K_c</math> and <math>~\rho_0 ~\Rightarrow</math>
  </td>
  <td colspan="2">
&nbsp;
  </td>
</tr>
<tr>
  <td align="right">
<math>~\rho</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>\rho_0 \theta^{n_c}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\rho_0 \biggl( \frac{\sin\xi}{\xi} \biggr)</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~P</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~K_c \rho_0^{1+1/n_c} \theta^{n_c + 1}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~K_c \rho_0^{2} \biggl( \frac{\sin\xi}{\xi}\biggr)^{2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~r</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{1/2} \rho_0^{(1-n_c)/(2n_c)} \xi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{K_c}{2\pi G} \biggr]^{1/2} \xi</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~M_r</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~4\pi \biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{3/2} \rho_0^{(3-n_c)/(2n_c)} \biggl(-\xi^2 \frac{d\theta}{d\xi} \biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~4\pi \biggl[ \frac{K_c}{2\pi G} \biggr]^{3/2} \rho_0 \biggl[\sin\xi  - \xi \cos\xi \biggr]</math>
  </td>
</tr>
</table>

</div>

==Step 5: Interface Conditions==

<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td colspan="3">
&nbsp;
  </td>
  <td align="left" colspan="2">
Setting <math>~n_c=1</math> and <math>~n_e=5~~~~~\Rightarrow</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\frac{\rho_e}{\rho_0}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{n_c}_i \phi_i^{-n_e}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \phi_i^{-5}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\biggl( \frac{K_e}{K_c} \biggr) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\rho_0^{1/n_c - 1/n_e}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-(1+1/n_e)} \theta^{1 - n_c/n_e}_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[\rho_0^{4}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-6} \theta^{4}_i\biggr]^{1/5}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\frac{\eta_i}{\xi_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{n_c + 1}{n_e+1} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \theta_i^{(n_c-1)/2} \phi_i^{(1-n_e)/2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{1}{3} \biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \phi_i^{-2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{n_c + 1}{n_e + 1} \biggr]^{1/2} \theta_i^{- (n_c + 1)/2} \phi_i^{(n_e+1)/2} \biggl( \frac{d\theta}{d\xi} \biggr)_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{1}{3} \biggr)^{1/2} \theta_i^{- 1} \biggl( \frac{d\theta}{d\xi} \biggr)_i \phi_i^3</math>
  </td>
</tr>
</table>
</div>

<span id="Alternative"><font color="red">Alternative:</font>  In our</span> [[SSC/Structure/BiPolytropes#UVplane|introductory description of how to build a bipolytropic structure]], we pointed out that, instead of employing these last two fitting conditions, Chandrasekhar [[Appendix/References#C67|[<b><font color="red">C67</font></b>]]] found it useful to employ, instead, the ratio of the <math>3^\mathrm{rd}</math> to <math>4^\mathrm{th}</math> expressions, which in the present case produces,
<div align="center">
<math>
\frac{\eta_i \phi_i^{5}}{(d\phi/d\eta)_i} = \frac{\xi_i \theta_i}{(d\theta/d\xi)_i} \biggl( \frac{\mu_e}{\mu_c}\biggr) \, ,
</math>
</div>
and the product of the <math>3^\mathrm{rd}</math> and <math>4^\mathrm{th}</math> expressions, which in the present case generates,
<div align="center">
<math>
\frac{3\eta_i (d\phi/d\eta)_i}{ \phi_i } = \frac{\xi_i (d\theta/d\xi)_i}{ \theta_i } \biggl( \frac{\mu_e}{\mu_c}\biggr) \, .
</math>
</div>
In what follows we will sometimes refer to the first of these two expressions as Chandrasekhar's "U-constraint" and we will sometimes refer to the second as Chandrasekhar's "V-constraint."  As is explained in [[SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar's_U_and_V_Functions|an accompanying discussion]], {{ Murphy83a }} followed Chandrasekhar's lead and extracted fitting conditions from this last pair of expressions.  In seeking the most compact analytic solution, we have found it advantageous to invoke our standard <math>3^\mathrm{rd}</math> fitting expression in tandem with the Chandrasekhar's V-constraint.

==Step 6:  Envelope Solution==

[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 20 April 2015:  There is a type-setting error in this function expression as published in the upper left-hand column of the second page of the article by Murphy (1983); the sine function in the denominator should be sine-squared, as presented here.]]Following the work of {{ Murphy83a }} and of {{ MF85a }}, we will adopt for the envelope's structure the F-Type solution of the n = 5 Lane-Emden function discovered by {{ Srivastava62full }} and described in an [[SSC/Structure/Polytropes/Analytic#Srivastava's_F-Type_Solution|accompanying discussion]], namely,

<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\phi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{B_0^{-1}\sin[\ln(A_0\eta)^{1/2})]}{\eta^{1/2}\{3-2\sin^2[\ln(A_0\eta)^{1/2}]\}^{1/2}} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{B_0^{-1}\sin\Delta}{\eta^{1/2}(3-2\sin^2\Delta)^{1/2}} \, ,</math>
  </td>
</tr>

</table>
</div>
[[File:CommentButton02.png|right|100px|Note that our homology factor and scaling coefficient serve virtually the same roles as the homology factor, A, and scaling coefficient, B, used by Murphy (1983), but they are not mathematically identical so we have added a subscript "0" to highlight the distinctions.]]
where <math>~A_0</math> is a "homology factor" and <math>~B_0</math> is an overall scaling coefficient &#8212; the values of both will be determined presently from the interface conditions &#8212; and we have introduced the notation,
<div align="center">
<math>~\Delta \equiv \ln(A_0\eta)^{1/2} = \frac{1}{2} (\ln A_0 + \ln\eta) \, .</math>
</div>

The first derivative of Srivastava's function is,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{d\phi}{d\eta}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{B_0^{-1}[3\cos\Delta-3\sin\Delta + 2\sin^3\Delta] }{2\eta^{3/2}(3-2\sin^2\Delta)^{3/2}}  \, .
</math>
  </td>
</tr>
</table>
</div>

As has been explained in the context of our [[SSC/Structure/Polytropes/Analytic#Srivastava's_F-Type_Solution|more general discussion of Srivastava's function]], if we ignore, for the moment, the additional "<math>m\pi</math>" phase shift that can be attached to a determination of the angle, <math>~\Delta</math>, the physically viable interval for the dimensionless radial coordinate is, <math>~e^{2\pi} \ge A_0\eta \ge \eta_\mathrm{crit} \equiv e^{2\tan^{-1}(1+2^{1/3})}  \, .</math>  


For this bipolytropic configuration, it is worth emphasizing how the dimensionless radial coordinate of the <math>~n_e = 5</math> envelope, <math>~\eta</math>, is related to the dimensionless radial coordinate of the <math>~n_c = 1</math> core, <math>~\xi</math>.  Referring to the general setup procedure for constructing any bipolytropic configuration that has been [[SSC/Structure/BiPolytropes#TableSetup|presented in tabular form in a separate discussion]], it is clear that in order for the radial coordinate, <math>~r</math>, to carry a consistent meaning throughout the model, we must have,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~r = \biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{1/2} \rho_0^{(1-n_c)/(2n_c)} \xi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{1/2} \rho_e^{(1-n_e)/(2n_e)} \eta</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~\biggl( \frac{K_c}{3K_e} \biggr)\rho_e^{4/5}  </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{\eta}{\xi}\biggr)^2 \, .</math>
  </td>
</tr>
</table>
</div>

Referring back to the [[SSC/Structure/BiPolytropes/Analytic15#Step_5:_Interface_Conditions|already established interface conditions, above]], to relate <math>~\rho_e</math> to <math>~\rho_0</math>, and to re-express the ratio, <math>~K_e/K_c</math>, we therefore have,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl(\frac{\eta}{\xi}\biggr)^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{3} \biggl[ \rho_0^{4}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-6} \theta^{4}_i\biggr]^{-1/5} 
\biggl[ \rho_0\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \phi_i^{-5} \biggr]^{4/5}  </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ \frac{\eta}{\xi}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{1}{3}\biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2} \, .</math>
  </td>
</tr>
</table>
</div>
While this result is not a surprise because the right-hand-side is the same expression that was presented, above, as the interface condition for the ratio, <math>~\eta_i/\xi_i</math>, it is nevertheless useful because it shows that the same relation works throughout the system &#8212; not just at the interface &#8212; and it clearly defines how we can swap back and forth between the two dimensionless radial coordinates when examining the structure and characteristics of this composite bipolytropic structure.

===First Constraint===
Calling upon Chandrasekhar's V-constraint, as just defined above &#8212; see also [[SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar.27s_U_and_V_Functions|our accompanying discussion]] for elaboration on [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy's (1983)] "V<sub>5F</sub>" and "V<sub>1E</sub>" function notations &#8212; one fitting condition at the interface is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~- \frac{2\xi_i }{ 3\theta_i } \biggl( \frac{d\theta}{d\xi} \biggr)_i \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2\biggl[ \frac{\eta (- d\phi/d\eta)}{\phi} \biggr]_i </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{[3\sin\Delta_i - 2\sin^3\Delta_i -3\cos\Delta_i ] }{\sin\Delta_i (3-2\sin^2\Delta_i)}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, .</math>
  </td>
</tr>
</table>
</div>
The left-hand side of this expression is inherently positive over the physically viable radial coordinate range, <math>~0 \ge \xi_i \ge \pi</math> and its value is known once the radial coordinate of the edge of the core has been specified.  So, defining the interface parameter,
<div align="center">
<math>~
\kappa_i \equiv - \frac{2\theta_i^' \xi_i}{3\theta_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)
\, ,</math>
</div>
we will recast the first constraint into, what will henceforth be referred to as, the
<div align="center" id="KeyInterfaceRelation">
<table border="1" cellpadding="8" align="center">
<tr><td align="center">

<font color="#770000">'''Key Nonlinear Interface Relation'''</font>
<br />
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\kappa_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, .</math>
  </td>
</tr>
</table>

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

In a [[SSC/Structure/BiPolytropes/Analytic15/Pt2#Analytic_Solution_of_Key_Interface_Relation|separate subsection (Part II) of this chapter]], we present a closed-form analytic solution, <math>\Delta_i(\kappa_i)</math>, to this nonlinear equation.


<!-- BEGIN EXTRACTION

===Second Constraint (Earliest Try)===

====Obtained from Third Interface Condition====
Now, our [[SSC/Structure/BiPolytropes/Analytic15#Step_5:_Interface_Conditions|3<sup>rd</sup> interface condition, as detailed above]], states that,
<div align="center">
<table border="0" cellpadding="5" align="center">

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)  \biggl[ \frac{\eta^{1/2} (3-2\sin^2\Delta_i)^{1/2}}{B^{-1}\sin\Delta_i} \biggr]^{2} </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ B^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\sqrt{3}}{\xi_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1}  \biggl( \frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1} \, .  </math>
  </td>
</tr>

</table>
</div>

END EXTRACTION -->

===Second Constraint===

====Obtained from Third Interface Condition====
Our [[SSC/Structure/BiPolytropes/Analytic15#Step_5:_Interface_Conditions|3<sup>rd</sup> interface condition, as detailed above]], states that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\eta_i}{\xi_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2} \, .</math>
  </td>
</tr>

</table>
</div>

If we now choose to normalize the interface amplitude such that, <math>~\phi_i = 1</math>, then this condition establishes two relations:  First, from the 3<sup>rd</sup> interface condition alone,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\eta_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \xi_i \, ;</math>
  </td>
</tr>

</table>
</div>
and, second, from the definition of Srivastava's function, <math>~\phi</math>, we deduce that the overall scaling parameter is,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~B_0^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\sin^2\Delta_i}{\eta_i(3-2\sin^2\Delta_i)} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\sqrt{3}}{\xi_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1}  \biggl( \frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1} \, .  </math>
  </td>
</tr>

</table>
</div>
Notice that, after the solution, <math>~\Delta_i(\kappa_i)</math>, of the [[SSC/Structure/BiPolytropes/Analytic15#KeyInterfaceRelation|key nonlinear interface relation]] has been determined, the first of these two relations also permits us to write,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~A_0\eta_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~e^{2\Delta_i}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ A_0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \xi_i^{-1}  e^{2\Delta_i} \, .</math>
  </td>
</tr>
</table>
</div>
Throughout the envelope, therefore, the angle,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\Delta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\ln(A_0\eta)^{1/2} = \frac{1}{2} \ln\biggl[ \xi \cdot \xi_i^{-1}  e^{2\Delta_i}  \biggr]
= \Delta_i + \ln\biggl( \frac{\xi}{\xi_i} \biggr)^{1/2} \, .</math>
  </td>
</tr>
</table>
</div>


====Obtained from Chandrasekhar's U-constraint====
We shall now demonstrate that the same expression for the scaling coefficient, <math>~B_0</math>, can alternatively be obtained from Chandrasekhar's U-constraint, without assuming that <math>~\phi_i = 1</math>, after taking into account the result that already has been obtained from the V-constraint.  As [[SSC/Structure/BiPolytropes/Analytic15#Step_5:_Interface_Conditions|described above]], the U-constraint is an alternative interface condition that may be written as,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\xi_i \theta_i}{(-d\theta/d\xi)_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\eta_i \phi_i^{5}}{(-d\phi/d\eta)_i} \, ,</math>
  </td>
</tr>
</table>
</div>
which, in the particular case being examined here, becomes &#8212; again, see [[SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar.27s_U_and_V_Functions|our accompanying discussion]] for elaboration on the "U<sub>5F</sub>" and "U<sub>1E</sub>" function notations used by {{ Murphy83a }} &#8212; 
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{2\xi_i^2}{3\kappa_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( U_\mathrm{5F} \biggr)_i </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{2B_0^{-4} \sin^4\Delta_i}{(3-2\sin^2\Delta_i)(3 - 2\sin^2\Delta_i - 3\cot\Delta_i)} \, .
</math>
  </td>
</tr>
</table>
</div>

Now, from our [[SSC/Structure/BiPolytropes/Analytic15#First_Constraint|discussion, above, of the first constraint]], we know that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~(3 - 2\sin^2\Delta_i - 3\cot\Delta_i)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(3-2\sin^2\Delta_i)\kappa_i \, .</math>
  </td>
</tr>
</table>
</div>

Hence, Chandrasekhar's U-constraint becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">

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

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

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ B_0^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\sqrt{3}}{\xi_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1}  \biggl( \frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1}  \, ,</math>
  </td>
</tr>
</table>
</div>
which, as predicted, is identical to what we learned from the third interface condition, alone.

===Comment on Murphy's Scalings===

The {{ Murphy83a }} derivations also include an homology factor, <math>~A</math>, and an overall scaling factor, <math>~B</math>, but they are calculated differently from our <math>~A_0</math> and <math>~B_0</math>.  In the righthand column of the third page of his paper, Murphy states that,
<div align="center">
<math>~A = \frac{\xi_J}{\zeta_J} \, ,</math>
</div>
which, when translated into our notation <math>~(\zeta_J \rightarrow \xi_i </math> and <math>~\xi_J \rightarrow A_0\eta_\mathrm{root})</math> gives,
<div align="center">
<math>~A = \frac{A_0 \eta_\mathrm{root}}{\xi_i} \, .</math>
</div>
Now, in our derivation, <math>~\eta_\mathrm{root}</math> is synonymous with the location of the envelope interface, <math>~\eta_i</math>, as expressed in terms of the dimensionless radial coordinate associated with Srivastava's Lane-Emden function, so we can equally well state that, 
<div align="center">
<math>~A = \frac{A_0 \eta_i}{\xi_i} \, .</math>
</div>
Recalling that <math>~\phi_i = 1</math>, we know from the [[SSC/Structure/BiPolytropes/Analytic15#Step_5:_Interface_Conditions|interface conditions detailed above]] that,
<div align="center">
<math>~\frac{\eta_i}{\xi_i} = \frac{1}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) \, .</math>
</div>
Hence Murphy's homology factor, <math>~A</math>, is related to our homology factor, <math>~A_0</math> via the expression,
<div align="center">
<math>~A =  \frac{A_0}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) \, .</math>
</div>
It is usually the value of this quantity, rather than simply our derived value of <math>~A_0</math>, that is tabulated below &#8212; both [[SSC/Structure/BiPolytropes/Analytic15#Murphy.27s_Example_Model_Characteristics|here]] and [[SSC/Structure/BiPolytropes/Analytic15#Murphy_and_Fiedler_.281985.29|here]] &#8212; as we make quantitative comparisons between the characteristics of our derived models and those published by {{ Murphy83a }} and by {{ MF85a }}.


In the lefthand column of the fourth page of his paper, {{ Murphy83a }} defines the coefficient <math>B</math> in such a way that the ''value'' of the envelope function, <math>\phi_{5F}</math>, equals the ''value'' of the core function, <math>\theta_{1E}</math>, at the interface.  Specifically, he sets,

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

<tr>
  <td align="right">
<math>~B</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{\zeta_J}{\sin\zeta_J} \biggr] \biggl[  \frac{A^{1/2} \sin(\ln\sqrt{A\zeta_J}) }{(A\zeta_J)^{1/2} \{2 + \cos[\ln(A\zeta_J)]\}^{1/2}}\biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{\zeta_J}{\sin\zeta_J} \biggr] \biggl[  \frac{\sin(\ln\sqrt{A\zeta_J}) }{\zeta_J^{1/2} \{3 - 2\sin^2(\ln\sqrt{A\zeta_J}) \}^{1/2}}\biggr] \, .</math>
  </td>
</tr>
</table>
</div>
Switching to our terminology, that is, setting, 
<div align="center">
<math>\ln\sqrt{A\zeta_J} \rightarrow \Delta_i</math>&nbsp; &nbsp;&nbsp; and, as before, &nbsp; &nbsp;&nbsp; 
<math>~\zeta_J \rightarrow \xi_i \, ,</math>
</div>
gives,

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

<tr>
  <td align="right">
<math>~B</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{\xi_i}{\sin\xi_i} \biggr] \biggl[  \frac{\sin\Delta_i }{\xi_i^{1/2}  (3 - 2\sin^2\Delta_i )^{1/2}}\biggr] </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\theta_i^{-1} \biggl[  \xi_i^{-1}  \biggl(\frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1} \biggr]^{1/2} \, .</math>
  </td>
</tr>
</table>
</div>
Hence, in terms of the definition of ''our'' scaling coefficient, <math>~B_0</math>, derived above, we have,

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

<tr>
  <td align="right">
<math>~B</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{B_0}{3^{1/4}} \biggl(  \frac{\mu_e}{\mu_c} \biggr)^{1/2} \theta_i^{-1} \, .</math>
  </td>
</tr>
</table>
</div>

As we make quantitative comparisons between the characteristics of our derived models and those published by [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] and by  [http://adsabs.harvard.edu/abs/1985PASAu...6..219M Murphy &amp; Fiedler (1985a)], below, we usually will tabulate the value of this quantity, 
rather than simply our derived value of <math>~B_0</math>.

==Step 7:  Identifying the Surface==

Because Shrivastava's function &#8212; and, along with it, the envelope's density &#8212; drops to zero when,
<div align="center">
<math>\Delta = \Delta_s \equiv \pi \, ,</math>
</div>
we know that the radius, <math>~\xi_s</math>, of the bipolytropic configuration is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\pi</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\Delta_i + \ln\biggl( \frac{\xi_s}{\xi_i} \biggr)^{1/2} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ \ln\biggl( \frac{\xi_s}{\xi_i} \biggr)^{1/2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\pi - \Delta_i 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ \xi_s</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\xi_i e^{2(\pi - \Delta_i )} \, .
</math>
  </td>
</tr>
</table>
</div>

In terms of the natural radial coordinate of the envelope, this is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\eta_s \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-1}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3^{-1/2} \xi_i e^{2(\pi - \Delta_i)} \, .</math>
  </td>
</tr>
</table>
</div>

=Key References=
* [http://adsabs.harvard.edu/abs/1962ApJ...136..680S S. Srivastava (1968, ApJ, 136, 680)] ''A New Solution of the Lane-Emden Equation of Index n = 5''
* [http://adsabs.harvard.edu/abs/1978AuJPh..31..115B H. A. Buchdahl (1978, Australian Journal of Physics, 31, 115)]:  ''Remark on the Polytrope of Index 5'' &#8212; the result of this work by Buchdahl has been [[SSC/Structure/BiPolytropes/Analytic51#Buchdahl1978|highlighted inside our discussion of bipolytropes with <math>~(n_c, n_e) = (5, 1)</math>]].
* [http://adsabs.harvard.edu/abs/1980PASAu...4...37M J. O. Murphy (1980a, Proc. Astr. Soc. of Australia, 4, 37)]:  ''A Finite Radius Solution for the Polytrope Index 5''
* [http://adsabs.harvard.edu/abs/1980PASAu...4...41M J. O. Murphy (1980b, Proc. Astr. Soc. of Australia, 4, 41)]:  ''On the F-Type and M-Type Solutions of the Lane-Emden Equation''
* [http://adsabs.harvard.edu/abs/1981PASAu...4..205M J. O. Murphy (1981, Proc. Astr. Soc. of Australia, 4, 205)]:  ''Physical Characteristics of a Polytrope Index 5 with Finite Radius''
* [http://adsabs.harvard.edu/abs/1982PASAu...4..376M J. O. Murphy (1982, Proc. Astr. Soc. of Australia, 4, 376)]:  ''A Sequence of E-Type Composite Analytical Solutions of the Lane-Emden Equation''
* [http://adsabs.harvard.edu/abs/1983AuJPh..36..453M J. O. Murphy (1983, Australian Journal of Physics, 36, 453)]:  ''Structure of a Sequence of Two-Zone Polytropic Stellar Models with Indices 0 and 1''
* [http://adsabs.harvard.edu/abs/1983PASAu...5..175M J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)]:  ''Composite and Analytical Solutions of the Lane-Emden Equation with Polytropic Indices n = 1 and n = 5''
* [http://adsabs.harvard.edu/abs/1985PASAu...6..219M J. O. Murphy &amp; R. Fiedler (1985a, Proc. Astr. Soc. of Australia, 6, 219)]:  ''Physical Structure of a Sequence of Two-Zone Polytropic Stellar Models''
* [http://adsabs.harvard.edu/abs/1985PASAu...6..222M J. O. Murphy &amp; R. Fiedler (1985b, Proc. Astr. Soc. of Australia, 6, 222)]:  ''Radial Pulsations and Vibrational Stability of a Sequence of Two-Zone Polytropic Stellar Models''

=Related Discussions=
* [[SSC/Structure/PolytropesEmbedded#Embedded_Polytropic_Spheres|Polytropes emdeded in an external medium]]
* [[SSC/Structure/BiPolytropes#BiPolytropes|Constructing BiPolytropes]]

{{ SGFfooter }}