Difference between revisions of "User:Tohline/SSC/Structure/BiPolytropes/Analytic1 5"

From VistrailsWiki
Jump to navigation Jump to search
(→‎Murphy's Example Model Characteristics: Extend format of "analytic model characteristics" table)
 
(88 intermediate revisions by the same user not shown)
Line 3: Line 3:
{{LSU_HBook_header}}
{{LSU_HBook_header}}


Here we construct a [[User:Tohline/SSC/Structure/BiPolytropes#BiPolytropes|bipolytrope]] 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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|separately discussed, "mirror image" bipolytropes 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.  [On '''<font color="red">12 April 2015</font>''', J. E. Tohline wrote:  I became aware of the published discussions of this system by Murphy &#8212; and especially the work of [http://adsabs.harvard.edu/abs/1985PASAu...6..222M Murphy &amp; Fiedler (1985)] &#8212; (see itemization of [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Key_References|additional key references, below]]) in March of 2015 after searching the internet for previous analyses of radial oscillations in polytropes and, then, reading through [http://adsabs.harvard.edu/abs/2004ASSL..306.....H Horedt's (2004)] &sect;2.8.1 discussion of composite polytropes.]   
[[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 [[User:Tohline/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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#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 [http://adsabs.harvard.edu/abs/1983PASAu...5..175M J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)], and attributes of their physical structure were further discussed by [http://adsabs.harvard.edu/abs/1985PASAu...6..219M J. O. Murphy &amp; R. Fiedler (1985a, Proc. Astr. Soc. of Australia, 6, 219)]; [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Key_References|additional, closely related references are given below]]. In the discussion that follows, we will be heavily referencing Murphy's (1983) work.
 
==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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#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''


==Steps 2 &amp; 3==
==Steps 2 &amp; 3==
Based on the discussion [[User:Tohline/SSC/Structure/Polytropes#.3D_1_Polytrope|presented elsewhere of the structure of an isolated <math>n=1</math> polytrope]], the core of this bipolytrope will have the following properties:
Based on the discussion [[User:Tohline/SSC/Structure/Polytropes#.3D_1_Polytrope|presented elsewhere of the structure of an isolated <math>~n=1</math> polytrope]], the core of this bipolytrope will have the following properties:
<div align="center">
<div align="center">
<math>
<math>
Line 127: Line 115:
   </td>
   </td>
   <td align="left" colspan="2">
   <td align="left" colspan="2">
Setting <math>~n_c=1</math>, <math>~n_e=5</math>, and <math>~\phi_i = \phi_i ~~~~\Rightarrow</math>
Setting <math>~n_c=1</math> and <math>~n_e=5~~~~~\Rightarrow</math>
   </td>
   </td>
</tr>
</tr>
Line 218: Line 206:
</math>
</math>
</div>
</div>
As is explained in [[User:Tohline/SSC/Structure/BiPolytropes/MurphyUVplane#UV_Plane_Functions_as_Analyzed_by_Murphy_.281983.29|an accompanying discussion]], [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] 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 combine our standard <math>3^\mathrm{rd}</math> fitting expression with the last (''i.e.,'' the product) expression identified by Chandrasekhar.
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 [[User:Tohline/SSC/Structure/BiPolytropes/MurphyUVplane#UV_Plane_Functions_as_Analyzed_by_Murphy_.281983.29|an accompanying discussion]], [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] 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==
==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 should be sine-squared, as presented here.]]Following the work of [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] and of [http://adsabs.harvard.edu/abs/1985PASAu...6..219M Murphy &amp; Fiedler (1985a)], we will adopt for the envelope's structure the F-Type solution of the <math>~n=5</math> Lane-Emden function discovered by [http://adsabs.harvard.edu/abs/1962ApJ...136..680S S. Srivastava (1968, ApJ, 136, 680)] and described in an [[User:Tohline/SSC/Structure/Polytropes#Srivastava.27s_F-Type_Solution|accompanying discussion]], namely,
[[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 [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] and of [http://adsabs.harvard.edu/abs/1985PASAu...6..219M Murphy &amp; Fiedler (1985a)], we will adopt for the envelope's structure the F-Type solution of the <math>~n=5</math> Lane-Emden function discovered by [http://adsabs.harvard.edu/abs/1962ApJ...136..680S S. Srivastava (1968, ApJ, 136, 680)] and described in an [[User:Tohline/SSC/Structure/Polytropes#Srivastava.27s_F-Type_Solution|accompanying discussion]], namely,


<div align="center">
<div align="center">
Line 234: Line 222:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{B^{-1}\sin[\ln(A\eta)^{1/2})]}{\eta^{1/2}\{3-2\sin^2[\ln(A\eta)^{1/2}]\}^{1/2}} \, ,</math>
<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>
   </td>
</tr>
</tr>


</table>
</div>
specifically over the physically viable interval, <math>~e^{2\pi} \ge A\eta \ge \eta_\mathrm{crit} \equiv e^{2\tan^{-1}(1+2^{1/3})}  \, .</math>  The first derivative of this function is,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{d\phi}{d\eta}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 253: Line 234:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{B_0^{-1}\sin\Delta}{\eta^{1/2}(3-2\sin^2\Delta)^{1/2}} \, ,</math>
\frac{B^{-1}[3\cos\Delta-3\sin\Delta + 2\sin^3\Delta] }{2\eta^{3/2}(3-2\sin^2\Delta)^{3/2}} \, ,
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
where,
[[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">
<div align="center">
<math>~\Delta \equiv \ln(A\eta)^{1/2} = \ln A^{1/2} + \ln\eta^{1/2} \, .</math>
<math>~\Delta \equiv \ln(A_0\eta)^{1/2} = \frac{1}{2} (\ln A_0 + \ln\eta) \, .</math>
</div>
</div>


===First Constraint===
The first derivative of Srivastava's function is,
Calling upon Chandrasekhar's "product" expression, as just defined above &#8212; see also [[User:Tohline/SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar.27s_U_and_V_Functions|our accompanying discussion]] for a clarification of the "V<sub>5F</sub>" function notation &#8212; one fitting condition at the interface is,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~- \frac{2\xi_i }{ 3\theta_i } \biggl( \frac{d\theta}{d\xi} \biggr)_i \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>
<math>~\frac{d\phi}{d\eta}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 278: Line 256:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~2\biggl( V_{5F} \biggr)_i </math>
<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>
   </td>
</tr>
</tr>
</table>
</div>
As has been explained in the context of our [[User:Tohline/SSC/Structure/Polytropes#Srivastava.27s_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 [[User:Tohline/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>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<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>
   <td align="center">
   <td align="center">
Line 290: Line 279:
   </td>
   </td>
   <td align="left">
   <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>
<math>~\biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{1/2} \rho_e^{(1-n_e)/(2n_e)} \eta</math>
   </td>
   </td>
</tr>
</tr>
Line 296: Line 285:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow ~~~~\biggl( \frac{K_c}{3K_e} \biggr)\rho_e^{4/5}  </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 302: Line 291:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, .</math>
<math>~\biggl(\frac{\eta}{\xi}\biggr)^2 \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</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>
</div>


and, as in our [[User:Tohline/SSC/Structure/Polytropes#Example_Interval|separate discussion of the properties of Srivastava's function]], adopting the shorthand notation,
Referring back to the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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">
<div align="center">
<math>~y_i \equiv \tan\Delta_i \, ,</math>
<table border="0" cellpadding="5" align="center">
</div>


this interface condition becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\kappa_i</math>
<math>~\biggl(\frac{\eta}{\xi}\biggr)^2</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 330: Line 309:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{1}{3} \biggl[ \rho_0^{4}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-6} \theta^{4}_i\biggr]^{-1/5}  
\frac{3 - 2y_i^2(1+y_i^2)^{-1}- 3y_i^{-1} }{3-2y_i^2(1+y_i^2)^{-1}}  
\biggl[ \rho_0\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \phi_i^{-5} \biggr]^{4/5} </math>
</math>
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow ~~~~ \frac{\eta}{\xi}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 344: Line 322:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl(\frac{1}{3}\biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2} \, .</math>
\frac{3y_i(1+y_i^2)- 2y_i^3 -3(1+y_i^2)}{3y_i(1+y_i^2)-2y_i^3}  
</math>
   </td>
   </td>
</tr>
</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 [[User:Tohline/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>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<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>
   <td align="center">
   <td align="center">
Line 358: Line 342:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~2\biggl[ \frac{\eta (- d\phi/d\eta)}{\phi} \biggr]_i </math>
\frac{y_i^3 -3y_i^2 + 3y_i -3  }{3y_i+y_i^3}  
</math>
   </td>
   </td>
</tr>
</tr>
Line 366: Line 348:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow~~~~ \kappa_i(3y_i+y_i^3)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 372: Line 354:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<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>
y_i^3 -3y_i^2 + 3y_i -3 
</math>
   </td>
   </td>
</tr>
</tr>
Line 380: Line 360:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow~~~~  y_i^3(1-\kappa_i) -3 y_i^2 + 3(1-\kappa_i)y_i -3</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 386: Line 366:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, .</math>
0 \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</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">


<div align="center">
<table border="1" cellpadding="8" width="80%">
<tr>
<tr>
  <td align="right">
<math>~\kappa_i</math>
  </td>
   <td align="center">
   <td align="center">
<font color="red">'''ASIDE:'''</font> Analytic Solution of Cubic Equation
<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>
   </td>
</tr>
</tr>
<tr>
</table>
  <td align="left">
 
As is well known and documented &#8212; see, for example [http://mathworld.wolfram.com/CubicFormula.html Wolfram MathWorld] or [http://en.wikipedia.org/wiki/Cubic_function Wikipedia's discussion] of the topic &#8212; the roots of any cubic equation can be determined analytically.  In order to evaluate the root(s) of our particular cubic equation, we have drawn from the utilitarian [http://www.math.vanderbilt.edu/~schectex/courses/cubic/ online summary provided by Eric Schechter at Vanderbilt University].  For a cubic equation of the general form,
</td></tr>
<div align="center">
</table>
<math>~ay^3 + by^2 + cy + d = 0 \, ,</math>
</div>
a real root is given by the expression,
<div align="center">
<math>~
y = p + \{q + [q^2 + (r-p^2)^3]^{1/2}\}^{1/3} + \{q - [q^2 + (r-p^2)^3]^{1/2}\}^{1/3}
\, ,</math>
</div>
where,
<div align="center">
<math>~p \equiv -\frac{b}{3a} \, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
<math>~q \equiv \biggl[p^3 + \frac{bc-3ad}{6a^2} \biggr] \, ,</math>
&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp;
<math>~r=\frac{c}{3a} \, .</math>
</div>
</div>
In our particular case,
 
In a [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Analytic_Solution_of_Key_Interface_Relation|separate subsection of this chapter, below]], 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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Step_5:_Interface_Conditions|3<sup>rd</sup> interface condition, as detailed above]], states that,
<div align="center">
<div align="center">
<math>~a =(1-\kappa_i)\, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
<table border="0" cellpadding="5" align="center">
<math>~b =-3\, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
<math>~c = 3(1-\kappa_i) \, ,</math>
&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp;
<math>~d = - 3 \, .</math>
</div>
Hence,
<div align="center">
<math>~p = \frac{1}{(1-\kappa_i)} \, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
<math>~r=+1 \, ,</math>
&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp;
<math>~q = p^3 = \frac{1}{(1-\kappa_i)^3} \, ,</math>
</div>
which implies that the real root, <math>~y_\mathrm{root}</math>, is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~y_\mathrm{root}</math>
<math>~\frac{\eta_i}{\xi_i}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 447: Line 423:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2}</math>
p + \{p^3 + [p^6 + (1-p^2)^3]^{1/2}\}^{1/3} + \{p^3 - [p^6 + (1-p^2)^3]^{1/2}\}^{1/3}
</math>
   </td>
   </td>
</tr>
</tr>
Line 461: Line 435:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<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>
p + \{p^3 + p^3[1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} + \{p^3 - p^3[1 + (p^{-2}-1)^3]^{1/2}\}^{1/3}  
</math>
   </td>
   </td>
</tr>
</tr>
Line 469: Line 441:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow ~~~~ p^{-1} y_\mathrm{root}</math>
<math>~\Rightarrow ~~~~ B^2</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 475: Line 447:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<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>
1 + \{1 + [1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} + \{1 - [1 + (p^{-2}-1)^3]^{1/2}\}^{1/3}
\, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>


(There is also a pair of imaginary roots, but they are irrelevant in the context of our overarching astrophysical discussion.)
END EXTRACTION -->
  </td>
 
</tr>
===Second Constraint===
</table>
</div>


In summary, then,
====Obtained from Third Interface Condition====
* Once the location, <math>~\xi_i</math>, of the outer edge of the core has been specified, which determines <math>~\theta_i</math> and <math>~\theta^'_i</math> as well, the value of the parameters, <math>~\kappa_i</math> and <math>~p</math>, are known via the expressions,
Our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Step_5:_Interface_Conditions|3<sup>rd</sup> interface condition, as detailed above]], states that,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 497: Line 465:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\kappa_i \equiv - \frac{2\theta_i^' \xi_i}{3\theta_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)</math>
<math>~\frac{\eta_i}{\xi_i}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 503: Line 471:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{2}{3} \biggl(1 - \xi_i \cot\xi_i  \biggr)\biggl( \frac{\mu_e}{\mu_c} \biggr) \, ,</math>
<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2} \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
and,
</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">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~p \equiv (1-\kappa_i)^{-1} </math>
<math>~\eta_i</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 518: Line 490:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{3}{3-2(\mu_e/\mu_c)(1-\xi_i \cot\xi_i)}  \, .</math>
<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \xi_i \, ;</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
 
and, second, from the definition of Srivastava's function, <math>~\phi</math>, we deduce that the overall scaling parameter is,
* The value of <math>~y_\mathrm{root}</math> is determined from the just-derived solution to the governing cubic equation, which then gives the interface value of the envelope parameter (adjusted by a <math>~\pi m</math> phase angle, where <math>~m</math> is an, as yet unspecified, integer),
<div align="center">
<div align="center">
<math>~\Delta_i = \tan^{-1}(y_\mathrm{root}) +\pi m \, .</math>
<table border="0" cellpadding="5" align="center">
</div>
* This, in turn, can be interpreted as a coordinate root &#8212; which we shall refer to as <math>~\eta_\mathrm{root}</math> &#8212; via the expression,
<div align="center">
<math>~(A\eta)_\mathrm{root} = e^{2\Delta_i} = e^{2\pi m} \cdot e^{2\tan^{-1}(y_\mathrm{root})}  \, .</math>
</div>
 
The following log-log plot shows what value(s) of <math>~(A\eta)_\mathrm{root} </math> (vertical axis) result from a wide range of interface coordinate specifications, <math>~\xi_i</math> (horizontal axis), under the simplifying assumption that <math>~(\mu_e/\mu_c) = 1</math>; solid blue diamonds identify roots resulting from setting <math>~m=0</math> in the phase-shift specification, solid green triangles identify roots for which <math>~m=-1</math>, and solid purple squares identify roots for which <math>~m=-2</math>.  Notice that when <math>~\xi_i \ll 1</math> &#8212; and, hence, <math>~\kappa_i \ll 1</math> and <math>~p \rightarrow 1</math>,
<div align="center">
<math>~y_\mathrm{root} \rightarrow (1 + 2^{1/3})  ~~~~~~~
\Rightarrow ~~~~~~~ (A\eta)_\mathrm{root} \rightarrow e^{2\pi m} \cdot e^{2\tan^{-1}(1 + 2^{1/3})} = 10.05837e^{2\pi m} \, .</math>
</div>
Hence, as displayed in the figure, the curve traced by blue diamonds (m = 0) asymptotically approaches (to the left) <math>~10.05837</math>; the curve traced by green triangles (m = -1) asymptotically approaches <math>1.87834\times 10^{-2}</math>; and the curve traced by purple squares (m = -2) asymptotically approaches <math>3.50770\times 10^{-5}</math>.
 
 
<div align="center">
<table border="1" cellpadding="5">
<tr>
<tr>
   <td align="center" colspan="2">
   <td align="right">
Examination of F2 Function Discussed by 
<math>~B_0^2</math>
[http://adsabs.harvard.edu/abs/1983PASAu...5..175M J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)]
   </td>
   </td>
</tr>
<tr>
   <td align="center">
   <td align="center">
[[File:MurphyF2roots02.png|center|500px|F2 Roots by Murphy (1983)]]
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\sin^2\Delta_i}{\eta_i(3-2\sin^2\Delta_i)} </math>
   </td>
   </td>
  <td align="center">
</tr>
<table border="0" cellpadding="0">
 
<tr>
<tr>
  <td align="right">
&nbsp;
  </td>
   <td align="center">
   <td align="center">
Murphy's<p></p>
<math>~=</math>
Roots of<p></p>
F2 = 0<p></p>
Function<p></p>
----
   </td>
   </td>
</tr>
<tr>
   <td align="left">
   <td align="left">
<font face="Courier">3.5076E-05</font>
<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>
   </td>
</tr>
</tr>
</table>
</div>
Notice that, after the solution, <math>~\Delta_i(\kappa_i)</math>, of the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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>
<tr>
  <td align="right">
<math>~A_0\eta_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
   <td align="left">
   <td align="left">
<font face="Courier">1.8785E-02</font>
<math>~e^{2\Delta_i}</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ A_0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
   <td align="left">
   <td align="left">
<font face="Courier">4.2993E+00</font>
<math>~3^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \xi_i^{-1}  e^{2\Delta_i} \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
Throughout the envelope, therefore, the angle,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="left">
   <td align="right">
<font face="Courier">7.6882E+00</font>
<math>~\Delta</math>
   </td>
   </td>
</tr>
   <td align="center">
<tr>
<math>~=</math>
   <td align="left">
<font face="Courier">1.0913E+01</font>
   </td>
   </td>
</tr>
<tr>
   <td align="left">
   <td align="left">
<font face="Courier">1.4101E+01</font>
<math>~\ln(A_0\eta)^{1/2} = \frac{1}{2} \ln\biggl[ \xi \cdot \xi_i^{-1}  e^{2\Delta_i}  \biggr]
  </td>
= \Delta_i + \ln\biggl( \frac{\xi}{\xi_i} \biggr)^{1/2} \, .</math>
</tr>
</table>
   </td>
   </td>
</tr>
</tr>
Line 600: Line 573:
</div>
</div>


The red line in the figure shows where the coordinate value along the vertical axis equals the coordinate value along the horizontal axis.  Hence the points where this red line intersects the other curves in the figure can be interpreted as identifying solutions to the cubic equation for which <math>~A = 1</math>, and,
 
====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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Step_5:_Interface_Conditions|described above]], the U-constraint is an alternative interface condition that may be written as,
<div align="center">
<div align="center">
<math>~\eta_\mathrm{root} = \xi_i \, .</math>
<table border="0" cellpadding="5" align="center">
</div>
These points are also what [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] refers to as roots of his <math>~F2 = 0</math> function.  The values of ten successive coordinate roots of this function are listed in the right-hand column of Table 1 in [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)].  We have re-listed the values of six of these roots in a column immediately to the right of the above figure and have marked with orange circles the locations of these six points in the figure.  As they should, Murphy's identified roots lie precisely at the intersection of the red line with the other curves.
 
===Murphy's Example Model Characteristics===
 
As in the previous figure, the following log-log plot shows what value(s) of <math>~(A\eta)_\mathrm{root} </math> (vertical axis) result from a wide range of interface coordinate specifications, <math>~\xi_i</math> (horizontal axis), under the simplifying assumption that <math>~(\mu_e/\mu_c) = 1</math>.  The solid blue diamonds present exactly the same information as in the previous figure, that is, they identify roots resulting from setting <math>~m=0</math> in the phase-shift specification.  But here, focusing on a regime that highlights larger roots of the first constraint equation, solid black squares identify roots for which <math>~m=+1</math>, and solid maroon triangles identify roots for which <math>~m=+2</math>


<div align="center">
<table border="1" cellpadding="5">
<tr>
<tr>
   <td align="center" colspan="2">
   <td align="right">
Examination of Table 3 Data Presented in 
<math>~\frac{\xi_i \theta_i}{(-d\theta/d\xi)_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>
[http://adsabs.harvard.edu/abs/1983PASAu...5..175M J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)]
   </td>
   </td>
</tr>
<tr>
   <td align="center">
   <td align="center">
[[File:MurphyRootsTable3.png|center|500px|Table 3 Model Characteristics by Murphy (1983)]]
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
<table border="0" cellpadding="5">
<math>~\frac{\eta_i \phi_i^{5}}{(-d\phi/d\eta)_i} \, ,</math>
<tr>
  <td align="center" colspan="2" width="100%">
[[File:CommentButton02.png|center|100px|Comment by J. E. Tohline on 29 April 2015:  As is argued more fully below, it appears as though the two parameter values drawn directly from Murphy's Table 3 that are shown here in red contain typographical errors. ]]
Murphy's Table 3 Data
  </td>
</tr>
<tr>
  <td align="center" width="50%">
<math>~\zeta_J \rightarrow \xi_i</math><p></p>
----
  </td>
  <td align="center" width="50%">
<math>~\xi_J \rightarrow (A\eta)_\mathrm{root}</math><p></p>
----
   </td>
   </td>
</tr>
</tr>
</table>
</div>
which, in the particular case being examined here, becomes &#8212; again, see [[User:Tohline/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)] "U<sub>5F</sub>" and "U<sub>1E</sub>" function notations &#8212;
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<font face="Courier">3.2678E-02</font>
<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>
   <td align="right">
   <td align="left">
<font face="Courier" color="red">10.0164</font>
<math>~\biggl( U_\mathrm{5F} \biggr)_i </math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
<font face="Courier">0.8154</font>
&nbsp;
   </td>
   </td>
   <td align="right">
   <td align="center">
<font face="Courier">12.0083</font>
<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>
   </td>
</tr>
</tr>
</table>
</div>
Now, from our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#First_Constraint|discussion, above, of the first constraint]], we know that,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<font face="Courier">1.6598</font>
<math>~(3 - 2\sin^2\Delta_i - 3\cot\Delta_i)</math>
   </td>
   </td>
   <td align="right">
   <td align="center">
<font face="Courier">20.4312</font>
<math>~=</math>
   </td>
   </td>
</tr>
   <td align="left">
<tr>
<math>~(3-2\sin^2\Delta_i)\kappa_i \, .</math>
   <td align="right">
<font face="Courier" color="red">2.6914</font>
  </td>
  <td align="right">
<font face="Courier">33.0249</font>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
Hence, Chandrasekhar's U-constraint becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<font face="Courier">2.7302</font>
<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>
   <td align="right">
   <td align="left">
<font face="Courier">100<math>\pi</math></font>
<math>~
\frac{2B_0^{-4} \sin^4\Delta_i}{(3-2\sin^2\Delta_i)^2 \kappa_i}
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
<font face="Courier">3.1415</font>
<math>~\Rightarrow ~~~~ B_0^4 </math>
   </td>
   </td>
   <td align="right">
   <td align="center">
<font face="Courier">91275<math>\pi</math></font>
<math>~=</math>
   </td>
   </td>
</tr>
  <td align="left">
</table>
<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>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
which, as predicted, is identical to what we learned from the third interface condition, alone.


In his Table 3, [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] presents six pairs of "coordinate" roots, <math>~(\zeta_J, \xi_J)</math>, that not only satisfy the first constraint equation but that also, in his interpretation, identify physically viable model solutions.  The values of all six of these coordinate pairs have been re-listed immediately to the right of the log-log plot, in a separate panel of the above figure.  After mapping Murphy's coordinate notation to ours &#8212; specifically invoking the mapping <math>~(\zeta_J, \xi_J) ~\rightarrow~ (\xi_i, A\eta_\mathrm{root})</math> &#8212; the locations of all six models have been marked by orange circles on the plot.  Five of Murphy's models fall directly onto our generated curves, giving us considerable confidence that our analysis matches his in considerable detail.  As is explained more fully, below, we strongly suspect that the single Murphy model that does not lie on one of our generated curves &#8212; the one which, as published in Murphy's Table 3, has coordinates<math>~(\zeta_J, \xi_J) = (2.6914, 33.0249)</math> &#8212; is misplaced because of a typographical error in the published table.
===Comment on Murphy's Scalings===


[http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy's (1983)] 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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Murphy.27s_Example_Model_Characteristics|here]] and [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Murphy_and_Fiedler_.281985.29|here]] &#8212; 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)].


As a supplement to the preceding two figures, the following table shows in quantitative detail how we are able to ''analytically'' derive individual coordinate roots, <math>~(A\eta)_\mathrm{root}</math> &#8212; see the next-to-last column of the table &#8212; from ten different interface-coordinate locations, <math>~\xi_i = \zeta_J</math>, as specified in the first (left-most) column of the table.  The interface parameter, <math>~\kappa_i</math>, is determined directly from the specified value of <math>~\xi_i</math>; the parameter, <math>~p=(1-\kappa_i)^{-1}</math>; then, given <math>~p</math>, the root to the cubic equation and, in turn, the value of <math>~\Delta_i</math>, are determined; and, finally, a choice of the integer phase-shift parameter, <math>~m</math>, leads to the calculation of <math>~(A\eta)_\mathrm{root}</math>.  For comparison, the last (right-most) column of the table lists the corresponding coordinate, <math>~\xi_J</math>, as originally derived and published by [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)].  Coordinate pairs shown in the orange-colored panels of our table correspond to the six models whose characteristics are provided in Murphy's Table 3; coordinate pairs shown in the green-colored panels correspond to two of the ten "F2 = 0" roots identified in Murphy's Table 1; and the additional three coordinate pairs shown in the blue-colored panels have been drawn from Murphy's Table 2. 


In the lefthand column of the fourth page of his paper, [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] 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">
<div align="center">
<table border="1" cellpadding="7">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="center" rowspan="2">
   <td align="right">
Murphy (1983)
<math>~B</math>
 
Tables 1, 2, &amp; 3
   </td>
   </td>
   <td align="center" colspan="6">
   <td align="center">
Here &hellip; setting <math>\xi_i = \zeta_J</math> &hellip; and assuming <math>\mu_e/\mu_c = 1</math>
<math>~=</math>
   </td>
   </td>
   <td align="center" rowspan="2">
   <td align="left">
Murphy (1983)
<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>
 
Tables 1, 2, &amp; 3
  </td>
</tr>
<tr>
  <td align="center" colspan="6">
1<sup>st</sup> Constraint
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center">
   <td align="right">
<math>~\zeta_J</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\kappa_i</math>
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
<math>~p</math>
<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>
   </td>
  <td align="center">
</tr>
<math>~y_\mathrm{root}</math>
</table>
  </td>
</div>
   <td align="center">
Switching to our terminology, that is, setting,
<math>~\Delta_i = \tan^{-1}(y_\mathrm{root})</math>
<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>
   <td align="center">
   <td align="center">
<math>~m </math>
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
<math>~(A\eta)_\mathrm{root} = e^{2(\Delta_i + m\pi)}</math>
<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>
  <td align="center">
<math>~\xi_J</math>
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center" rowspan="4" bgcolor="#FFBB88">
   <td align="right">
3.2678E-02
&nbsp;
  </td>
  <td align="center" rowspan="4">
2.3732E-04
  </td>
  <td align="center" rowspan="4">
1.00023737
  </td>
  <td align="center" rowspan="4">
2.26083
  </td>
  <td align="center" rowspan="4">
1.15435
   </td>
   </td>
   <td align="center">
   <td align="center">
-3
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
6.55237E-08
<math>~\theta_i^{-1} \biggl[  \xi_i^{-1}  \biggl(\frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1} \biggr]^{1/2} \, .</math>
  </td>
  <td align="center" bgcolor="lightblue">
6.55E-08
   </td>
   </td>
</tr>
</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>
<tr>
   <td align="center">
   <td align="right">
-2
<math>~B</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
3.50874E-05
<math>~=</math>
   </td>
   </td>
   <td align="center" bgcolor="lightblue">
   <td align="left">
3.508E-05
<math>~\frac{B_0}{3^{1/4}} \biggl(  \frac{\mu_e}{\mu_c} \biggr)^{1/2} \theta_i^{-1} \, .</math>
   </td>
   </td>
</tr>
</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>
<tr>
   <td align="center">
   <td align="right">
-1
<math>~\pi</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
1.87890E-02
<math>~=</math>
   </td>
   </td>
   <td align="center" bgcolor="lightblue">
   <td align="left">
1.8789E-02
<math>~
\Delta_i + \ln\biggl( \frac{\xi_s}{\xi_i} \biggr)^{1/2}
</math>
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center">
   <td align="right">
0
<math>~\Rightarrow ~~~~ \ln\biggl( \frac{\xi_s}{\xi_i} \biggr)^{1/2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
10.0614
<math>~=</math>
   </td>
   </td>
   <td align="center" bgcolor="#FFBB88">
   <td align="left">
<font color="red">10.0164</font>
<math>~
\pi - \Delta_i
</math>
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center" bgcolor="#FFBB88">
   <td align="right">
0.8154
<math>~\Rightarrow ~~~~ \xi_s</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
0.15474
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
1.18307
<math>~
  </td>
\xi_i e^{2(\pi - \Delta_i )} \, .
  <td align="center">
</math>
2.93939
  </td>
  <td align="center">
1.24287
  </td>
  <td align="center">
0
  </td>
  <td align="center">
12.0101
  </td>
  <td align="center" bgcolor="#FFBB88">
12.0083
   </td>
   </td>
</tr>
</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>
<tr>
   <td align="center" bgcolor="#FFBB88">
   <td align="right">
1.6598
<math>~\eta_s \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-1}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
0.76541
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
4.26281
<math>~3^{-1/2} \xi_i e^{2(\pi - \Delta_i)} \, .</math>
  </td>
  <td align="center">
12.6311
  </td>
  <td align="center">
1.49179
  </td>
  <td align="center">
0
  </td>
  <td align="center">
19.7585
  </td>
  <td align="center" bgcolor="#FFBB88">
20.4312
   </td>
   </td>
</tr>
</tr>
</table>
</div>
==Analytic Solution of Key Interface Relation==
Returning to our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#KeyInterfaceRelation|previously derived]],
<div align="center" id="KeyInterfaceRelation">
<font color="#770000">'''Key Nonlinear Interface Relation'''</font>
<br />
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="center">
   <td align="right">
[&nbsp;2.0914&nbsp;]
<math>~\kappa_i</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
1.46609
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
-2.14552
<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, ,</math>
   </td>
   </td>
   <td align="center">
</tr>
-6.11816
</table>
 
</div>
 
and, as in our [[User:Tohline/SSC/Structure/Polytropes#Example_Interval|separate discussion of the properties of Srivastava's function]], adopting the shorthand notation,
<div align="center">
<math>~y_i \equiv \tan\Delta_i \, ,</math>
</div>
 
this key interface condition becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
   <td align="right">
<math>~\kappa_i</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
-1.40878
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
1
<math>~
  </td>
\frac{3 - 2y_i^2(1+y_i^2)^{-1}- 3y_i^{-1}  }{3-2y_i^2(1+y_i^2)^{-1}}
  <td align="center">
</math>
31.9964
  </td>
  <td align="center" rowspan="3" bgcolor="#FFBB88">
33.0249
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center">
   <td align="right">
[&nbsp;2.10961&nbsp;]
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
1.50744
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
-1.97067
<math>~
\frac{3y_i(1+y_i^2)- 2y_i^3 -3(1+y_i^2)}{3y_i(1+y_i^2)-2y_i^3}
</math>
   </td>
   </td>
   <td align="center">
</tr>
-5.56379
 
<tr>
   <td align="right">
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
-1.39296
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
1
<math>~
  </td>
\frac{y_i^3 -3y_i^2 + 3y_i -3  }{3y_i+y_i^3}
  <td align="center">
</math>
33.0249
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center" bgcolor="#FFBB88">
   <td align="right">
<font color="red">2.6914</font>
<math>~\Rightarrow~~~~ \kappa_i(3y_i+y_i^3)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
4.37925
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
-0.29592
<math>~
y_i^3 -3y_i^2 + 3y_i -3 
</math>
   </td>
   </td>
   <td align="center">
</tr>
-0.31486
 
<tr>
   <td align="right">
<math>~\Rightarrow~~~~  y_i^3(1-\kappa_i) -3 y_i^2 + 3(1-\kappa_i)y_i -3</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
-0.30503
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
1
<math>~
  </td>
0 \, .
  <td align="center">
</math>
290.943
   </td>
   </td>
</tr>
</tr>
</table>
</div>
This is a cubic equation whose solution, <math>~y_\mathrm{root}(\kappa_i)</math>, will also immediately supply the desired interface angle, <math>~\Delta_i</math>, and an interface coordinate root, <math>~(A_0\eta)_\mathrm{root}</math>. 


<div align="center">
<table border="1" cellpadding="8" width="80%">
<tr>
<tr>
  <td align="center" bgcolor="#FFBB88">
2.7302
  </td>
   <td align="center">
   <td align="center">
4.83853
<font color="red">'''ASIDE:'''</font>  Analytic Solution of Cubic Equation
   </td>
   </td>
   <td align="center">
</tr>
-0.26052
<tr>
  </td>
   <td align="left">
  <td align="center">
As is well known and documented &#8212; see, for example [http://mathworld.wolfram.com/CubicFormula.html Wolfram MathWorld] or [http://en.wikipedia.org/wiki/Cubic_function Wikipedia's discussion] of the topic &#8212; the roots of any cubic equation can be determined analytically.  In order to evaluate the root(s) of our particular cubic equation, we have drawn from the utilitarian [http://www.math.vanderbilt.edu/~schectex/courses/cubic/ online summary provided by Eric Schechter at Vanderbilt University].  For a cubic equation of the general form,
-0.27316
<div align="center">
  </td>
<math>~ay^3 + by^2 + cy + d = 0 \, ,</math>
  <td align="center">
</div>  
-0.26666
a real root is given by the expression,
  </td>
<div align="center">
  <td align="center">
<math>~
1
y = p + \{q + [q^2 + (r-p^2)^3]^{1/2}\}^{1/3} + \{q - [q^2 + (r-p^2)^3]^{1/2}\}^{1/3}
  </td>
\, ,</math>
  <td align="center">
</div>  
99.997<math>\pi</math>
where,
  </td>
<div align="center">
  <td align="center" bgcolor="#FFBB88">
<math>~p \equiv -\frac{b}{3a} \, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
100<math>\pi</math>
<math>~q \equiv \biggl[p^3 + \frac{bc-3ad}{6a^2} \biggr] \, ,</math>
  </td>
&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp;
</tr>
<math>~r=\frac{c}{3a} \, .</math>
 
</div>
<tr>
In our particular case,
  <td align="center" bgcolor="#FFBB88">
<div align="center">
3.1415
<math>~a =(1-\kappa_i)\, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
  </td>
<math>~b =-3\, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
   <td align="center">
<math>~c = 3(1-\kappa_i) \, ,</math>
22604.6
&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp;
<math>~d = - 3 \, .</math>
</div>
Hence,
<div align="center">
<math>~p = \frac{1}{(1-\kappa_i)} \, ,</math> &nbsp;&nbsp;&nbsp;&nbsp;
<math>~r=+1 \, ,</math>
&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp;
<math>~q = p^3 = \frac{1}{(1-\kappa_i)^3} \, ,</math>
</div>
which implies that the real root, <math>~y_\mathrm{root}</math>, is given by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~y_\mathrm{root}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
-4.42408E-05
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
-4.42408E-05
<math>~
p + \{p^3 + [p^6 + (1-p^2)^3]^{1/2}\}^{1/3} + \{p^3 - [p^6 + (1-p^2)^3]^{1/2}\}^{1/3}
</math>
   </td>
   </td>
   <td align="center">
</tr>
-4.42408E-05
 
<tr>
   <td align="right">
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
2
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
91268<math>\pi</math>
<math>~
  </td>
p + \{p^3 + p^3[1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} + \{p^3 - p^3[1 + (p^{-2}-1)^3]^{1/2}\}^{1/3}
  <td align="center" bgcolor="#FFBB88">
</math>
91275<math>\pi</math>
   </td>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="center" bgcolor="lightgreen">
   <td align="right">
4.2993
<math>~\Rightarrow ~~~~ p^{-1} y_\mathrm{root}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
-0.58961
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
0.62908
<math>~
1 + \{1 + [1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} + \{1 - [1 + (p^{-2}-1)^3]^{1/2}\}^{1/3}
\, .
</math>
   </td>
   </td>
  <td align="center">
</tr>
0.89355
</table>
</div>
 
(There is also a pair of imaginary roots, but they are irrelevant in the context of our overarching astrophysical discussion.)
   </td>
   </td>
   <td align="center">
</tr>
0.72924
</table>
</div>
 
In summary, then,
* Once the location, <math>~\xi_i</math>, of the outer edge of the core has been specified, which determines <math>~\theta_i</math> and <math>~\theta^'_i</math> as well, the value of the parameters, <math>~\kappa_i</math> and <math>~p</math>, are known via the expressions,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~\kappa_i \equiv - \frac{2\theta_i^' \xi_i}{3\theta_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
0
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
4.2994
<math>~\frac{2}{3} \biggl(1 - \xi_i \cot\xi_i  \biggr)\biggl( \frac{\mu_e}{\mu_c} \biggr) \, ,</math>
  </td>
  <td align="center" bgcolor="lightgreen">
4.2993
   </td>
   </td>
</tr>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="center" bgcolor="lightgreen">
   <td align="right">
7.6882
<math>~p \equiv (1-\kappa_i)^{-1} </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
-0.19091
<math>~=</math>
   </td>
   </td>
   <td align="center">
   <td align="left">
0.83969
<math>~\frac{3}{3-2(\mu_e/\mu_c)(1-\xi_i \cot\xi_i)}  \, .</math>
  </td>
  <td align="center">
1.6268
  </td>
  <td align="center">
1.01964
  </td>
  <td align="center">
0
  </td>
  <td align="center">
7.6851
  </td>
  <td align="center" bgcolor="lightgreen">
7.6882
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>


 
* The value of <math>~y_\mathrm{root}</math> is determined from the just-derived solution to the governing cubic equation, which then gives the interface value of the envelope parameter (adjusted by a <math>~\pi m</math> phase angle, where <math>~m</math> is an, as yet unspecified, integer),
 
<div align="center">
As this table and its two preceding figures illustrate, by and large, the agreement between our ''analytically'' determined coordinate roots and those derived &amp; published by [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] is excellent.  There are a couple of mismatches, but we attribute both to typesetting errors in the Murphy publication: 
<math>~\Delta_i = \tan^{-1}(y_\mathrm{root}) +\pi m \, .</math>
* The coordinate, <math>~\xi_J = 10.0164</math>, that has been displayed using a red font in the right-most column of our table is the value drawn directly from Murphy's Table 3.  It appears to us that two digits of this number were inadvertently transposed during publication and that the coordinate value should have appeared as <math>~\xi_J = 10.0614</math>.  We are convinced that, as published, this is a typographical error in Murphy's Table 3 because:  (a) the corrected value matches our derived value of that coordinate to more significant figures, and (b) the same model is also listed in Murphy's Table 2 with this corrected value.
</div>
* The coordinate, <math>~\zeta_J = 2.6914</math>, that has been displayed using a red font in the left-most column of our table is the value drawn directly from Murphy's Table 3.  As shown in our table, from this coordinate value we derive a root coordinate of <math>~(A\eta)_\mathrm{root} = 290.943</math>, which is almost a factor of ten larger than Murphy's published root coordinate value, <math>~\xi_J = 33.0249</math>.  Working backwards via an iterative procedure, we have determined that Murphy's listed root of <math>33.0249</math> can be obtained if the originating value of <math>~\zeta_J \leftrightarrow \xi_i</math> is set to <math>2.10961</math> (see the black-on-white panel entry in the left-most column of our table).  By replacing the "6" by a zero in Murphy's published coordinate parameter &#8212; that is, by assuming that there was a typographical error and Murphy's interface coordinate was ''actually'' <math>~\zeta_J = 2.0914</math>, we are able to derive a root coordinate value (31.9964) that is quite close to Murphy's published value (33.0249).
* This, in turn, can be interpreted as a coordinate root &#8212; which we shall refer to as <math>~\eta_\mathrm{root}</math> &#8212; via the expression,
 
 
If, following Murphy's lead, we now consider the dimensionless radial variable throughout the <math>~n=5</math> envelope, <math>~\eta</math>, to be simply a direct extension of the dimensionless radial variable, <math>~\xi</math>, that is used throughout the <math>~n=1</math> core &#8212; in which case, the interface between the core and the envelope is not just associated with <math>~\xi_i</math>, but also with <math>~\eta_\mathrm{root} = \xi_i</math> &#8212; then the so-called ''homology factor'', <math>~A</math>, in Srivastava's Lane-Emden function expression can be determined for each tabulated  ''interface'' coordinate pair via the ratio,
<div align="center">
<div align="center">
<math>~A = \frac{(A\eta)_\mathrm{root}}{\xi_i} \, .</math>
<math>~(A_0\eta)_\mathrm{root} = e^{2\Delta_i} = e^{2\pi m} \cdot e^{2\tan^{-1}(y_\mathrm{root})} \, .</math>
</div>
</div>
Furthermore, it is clear that the surface of the bipolytropic configuration is identified by the value of the radial coordinate, <math>~\xi_s = \eta_s</math>, at which the numerator of [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Step_6:__Envelope_Solution|Srivastava's function]], <math>~\phi</math>, drops to zero.  That is, the surface occurs where,
 


<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="1" cellpadding="8" width="80%">
 
<tr>
  <td align="center">
<font color="red">'''CAUTION:'''</font>  Solution Behavior When <math>~\kappa_i = 1</math> and <math>~\Delta_i = \pi/2 \, .</math>
  </td>
</tr>
<tr>
  <td align="left">
In building a ''sequence'' of bipolytropic configurations having <math>~(n_c, n_e) = (1, 5)</math>, it will make sense to steadily increase the value of the parameter that marks the edge of the core, <math>~\xi_i</math>, from zero &#8212; meaning no mass in the core &#8212; to its maximum allowed value, <math>~\pi</math> &#8212; meaning no mass in the envelope.  (A more complete discussion of physically viable parameter values is [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Physically_Viable_Parameter_Values|presented below]].)  As <math>~\xi_i</math> steadily increases from zero, for a while the parameter, <math>~\kappa_i</math>, will steadily increase from zero as well. As the value of <math>~\kappa_i</math> crosses through the value of "one," the associated parameter, <math>~p = (1-\kappa_i)^{-1}</math>, as well as the root of the cubic equation, <math>~y_\mathrm{root}</math>, will flip from positive infinity to negative infinity.  From the standpoint of the construction of physically realistic models, this does not pose a problem.  It simply reflects the fact that, when <math>~\kappa_i = 1</math>, the governing interface angle, <math>~\Delta_i = \pi/2</math> &#8212; or, allowing for the aforementioned phase shift, <math>~\Delta_i = \pi/2 + m\pi</math>. 
 
 
From a practical standpoint, however, it is useful to keep in mind that the root of the governing cubic equation will change abruptly from <math>~+ \infty</math> to <math>~- \infty</math> and its value will be quite sensitive to the choice of <math>~\xi_i</math> in the vicinity of <math>~\kappa_i = 1</math>.  Given the definition of the function, <math>~\kappa_i(\xi_i)</math>, this means that the abrupt transition will occur at an interface location, <math>~\xi_\mathrm{trans}</math>, whose value satisfies the condition,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Delta_s</math>
<math>~\frac{\xi_\mathrm{trans}}{\tan(\xi_\mathrm{trans})}  </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,077: Line 1,142:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\pi</math>
<math>~1-\frac{3}{2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \, .</math>
   </td>
   </td>
</tr>
</tr>


</table>
</div>
For the case, <math>~\mu_e/\mu_c = 1 \, ,</math> we have determined empirically (to a dozen significant digits), that,
<div align="center">
<math>~\xi_\mathrm{trans} = 1.836597203152 \, .</math>
</div>
Just for fun &#8212; and, again, for the case, <math>~\mu_e/\mu_c = 1 </math> &#8212;  we also have determined that <math>~\xi_\mathrm{trans}</math> satisfies the series,
<div align="center">
<math>~
0 = \sum_{n=0}^{\infty} (-1)^{n+1} \biggl(\frac{4n+3}{2n+1} \biggr) \frac{\xi_\mathrm{trans}^{2n}}{(2n)!}
\, .</math>
</div>
The black, vertical, long-dashed line in the following figure identifies where <math>~\xi_i = \xi_\mathrm{trans}</math>. 
Because the function <math>~\tan(\xi_\mathrm{trans})</math> is periodic, other values of <math>~\xi_\mathrm{trans}</math> will also satisfy the condition,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow ~~~~\xi_s = \eta_s</math>
<math>~\frac{\xi_\mathrm{trans}}{\tan(\xi_\mathrm{trans})}  </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,089: Line 1,175:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{A} \cdot e^{2\pi} \, .</math>
<math>~1-\frac{3}{2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
But the one whose value we have determined here is the only one that falls within the physically viable range of interface location values, <math>~0 \le \xi_i \le \pi </math>.


We can now construct a table of our ''analytically determined'' model characteristics that is patterned after and largely matches [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy's (1983)] Table 3, with only the overall scaling parameter, <math>~B</math>, yet to be determined.
  </td>
<div align="center">
</tr>
</table>
</div>
 
 
<span id="MurphyF2Function">The following log-log plot shows what value(s) of <math>~(A_0\eta)_\mathrm{root} </math> (vertical axis) result from a wide range of interface coordinate specifications, <math>~\xi_i</math> (horizontal axis), under the simplifying assumption that <math>~(\mu_e/\mu_c) = 1</math>: Solid blue diamonds identify roots resulting from setting <math>~m=0</math> in the phase-shift specification, solid green triangles identify roots for which <math>~m=-1</math>, and solid purple squares identify roots for which <math>~m=-2</math>. </span>  As has just been discussed in the context of a cautionary note, as <math>~\xi_i</math> is increased from zero for any given value of <math>~m</math> &#8212; that is, as one traverses the plot, moving from the lefthand edge toward the right along a curve of ''constant color'' &#8212; the solution, <math>~(A_0\eta)_\mathrm{root} </math>, makes an abrupt jump when <math>~\xi_i</math> crosses the value, <math>~\xi_\mathrm{trans}</math>.  This identifies the root of the cubic equation for which <math>~\Delta_i = \pi/2 + m\pi</math>.  One can continue to move in a smooth, continuous fashion along a single &#8212; but multi-colored &#8212; solution curve by letting <math>~m \rightarrow (m+1)</math> as <math>~\xi_i</math> crosses <math>~\xi_\mathrm{trans}</math>.
 
 
<div align="center">
<table border="1" cellpadding="5">
<table border="1" cellpadding="5">
<tr>
<tr>
  <td align="center" colspan="11">
Some Model Characteristics
  </td>
</tr>
<tr>
  <td align="center" rowspan="2">
Model
  </td>
   <td align="center" colspan="2">
   <td align="center" colspan="2">
Specified
Examination of F2 Function Discussed by 
  </td>
[http://adsabs.harvard.edu/abs/1983PASAu...5..175M J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)]
  <td align="center" colspan="4">
Analytically Determined Here
  </td>
  <td align="center" colspan="4">
&hellip; and 1<sup>st</sup> Constraint Implications<sup>&dagger;</sup>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center">
   <td align="center">
<math>~\xi_i</math>
[[File:MurphyF2roots05.png|center|500px|F2 Roots by Murphy (1983)]]
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~m</math>
<table border="0" cellpadding="0">
  </td>
<tr>
   <td align="center">
   <td align="center">
<math>~(A\eta)_\mathrm{root}</math>
Murphy's<p></p>
Roots of<p></p>
F2 = 0<p></p>
Function<p></p>
----
   </td>
   </td>
   <td align="center">
</tr>
<math>~A</math>
<tr>
   <td align="left">
<font face="Courier">3.5076E-05</font>
   </td>
   </td>
   <td align="center">
</tr>
<math>~B</math>
<tr>
   <td align="left">
<font face="Courier">1.8785E-02</font>
   </td>
   </td>
   <td align="center">
</tr>
<math>~\xi_s</math>
<tr>
   <td align="left">
<font face="Courier">4.2993E+00</font>
   </td>
   </td>
   <td align="center">
</tr>
<math>~\Delta_i = \ln(A\eta)_\mathrm{root}^{1/2}</math>
<tr>
   <td align="left">
<font face="Courier">7.6882E+00</font>
   </td>
   </td>
   <td align="center">
</tr>
<math>~V_{1E}</math>
<tr>
   <td align="left">
<font face="Courier">1.0913E+01</font>
   </td>
   </td>
   <td align="center">
</tr>
<math>~V_{5F}</math>
<tr>
   <td align="left">
<font face="Courier">1.4101E+01</font>
   </td>
   </td>
  <td align="center">
</tr>
<math>\frac{(n_c+1)V_{1E}}{(n_e+1)V_{5F}}</math>
</table>
   </td>
   </td>
</tr>
</tr>
</table>
</div>


<tr>
The red line in the figure shows where the coordinate value along the vertical axis equals the coordinate value along the horizontal axis.  Hence the points where this red line intersects the other curves in the figure can be interpreted as identifying solutions to the cubic equation for which,
  <td align="center">1</td>
<div align="center">
  <td align="right">3.2678E-02</td>
<math>~(A_0\eta)_\mathrm{root} = \xi_i \, .</math>
  <td align="center">0</td>
</div>
  <td align="right">10.0614</td>
 
  <td align="right">307.894</td>
The points of intersection that are highlighted in this figure are also what [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] refers to as roots of his <math>~F2 = 0</math> function.  The values of ten successive coordinate roots of this function are listed in the right-hand column of Table 1 in [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)]. We have re-listed the values of six of these roots in a column immediately to the right of the above figure and have marked with orange circles the locations of these six points in the figure.  As they should, Murphy's identified roots lie precisely at the intersection of the red line with the other curves.
  <td align="center">---</td>
 
  <td align="right">1.7392</td>
==Physically Viable Parameter Values==
  <td align="center">---</td>
===Interface Mapping===
  <td align="center">---</td>
 
  <td align="center">---</td>
As has already been stated in our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Steps_2_.26_3|above description of the <math>~n_c = 1</math> core]] of these bipolytropic configurations, the structure of the core will be physically viable as long as the radial location of the interface, <math>~\xi_i</math>, between the core and the envelope is positioned somewhere within the range,
  <td align="center">---</td>
<div align="center">
</tr>
<math>~0 \le \xi_i \le \pi \, .</math>
</div>
 
Similarly, our [[User:Tohline/SSC/Structure/Polytropes#Srivastava.27s_F-Type_Solution|discussion of Srivastava's Lane-Emden function]], which is being used to define the envelope of these bipolytropic configurations, makes it clear that the envelope will have a physically viable structure as long as the parameter, <math>~\Delta_i</math>, associated with the radial location of the interface is positioned somewhere within the range,
<div align="center">
<math>~\eta_\mathrm{crit}  < e^{2\Delta_i} < e^{2\pi} \, ,</math>
</div>
where,
<div align="center">
<math>\eta_\mathrm{crit} \equiv e^{2\tan^{-1}(1+2^{1/3})} = 10.05836783\, ,</math>
</div>
and, in keeping with the definition provided above,
<div align="center">
<math>\Delta_i = \tan^{-1}(y_\mathrm{root}) + m\pi \, .</math>
</div>
This coordinate range for the physically viable envelope can be rewritten as,
<div align="center">
<math>~\tan^{-1}(1+2^{1/3})  \le \Delta_i \le \pi \, .</math>
</div>
 
As is illustrated in the following figure, titled "Interface Mapping," our analytic solution defines a smooth, one-to-one mapping of the radial coordinate that defines the outer edge of the <math>~n_c=1</math> core, <math>~\xi_i</math>, to the parameter, <math>~\Delta_i</math>,  that defines the inner edge of the <math>~n_e = 5</math> envelope. 


<tr>
  <td align="center">2</td>
  <td align="right">0.8154</td>
  <td align="center">0</td>
  <td align="right">12.0101</td>
  <td align="right">14.729</td>
  <td align="center">---</td>
  <td align="right">36.356</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>


<div align="center">
<table border="1" cellpadding="0" align="center">
<tr>
<tr>
   <td align="center">3</td>
   <td align="center">
   <td align="right">1.6598</td>
[[File:InterfaceMap02.png|center|500px|Illustration of Interface Mapping]]
  <td align="center">0</td>
  </td>
  <td align="right">19.7584</td>
   <td align="center">
  <td align="right">11.9041</td>
<table border="0" cellpadding="5" align="center">
  <td align="center">---</td>
<tr>
  <td align="right">44.984</td>
<td align="center"><math>~\mathrm{Model}</math></td>
  <td align="center">---</td>
<td align="center"><math>~\xi_i</math></td>
  <td align="center">---</td>
<td align="center"><math>~m</math></td>
  <td align="center">---</td>
<td align="center"><math>~\Delta_i</math></td>
  <td align="center">---</td>
</tr>
<tr>
<td align="center">
----</td>
<td align="center">
----</td>
<td align="center">
----</td>
<td align="center">
----</td>
</tr>
<tr>
<td align="center">&hellip;</td>
<td align="center"><math>~0</math></td>
<td align="center"><math>~0</math></td>
<td align="center"><math>~\tan^{-1}(1+2^{1/3})</math></td>
</tr>
<tr>
<td align="center"><math>~(2)</math></td>
<td align="center"><math>~0.8154</math></td>
<td align="center"><math>~0</math></td>
<td align="center"><math>~1.24287</math></td>
</tr>
</tr>
<tr>
<tr>
  <td align="center">4</td>
<td align="center"><math>~(3)</math></td>
  <td align="right">2.0914</td>
<td align="center"><math>~1.6598</math></td>
  <td align="center">1</td>
<td align="center"><math>~0</math></td>
  <td align="right">31.9965</td>
<td align="center"><math>~1.49179</math></td>
  <td align="right">15.2991</td>
  <td align="center">---</td>
  <td align="right">35.0016</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
</tr>
<tr>
<tr>
  <td align="center">5</td>
<td align="center"><math>~(4)</math></td>
  <td align="right">2.7302</td>
<td align="center"><math>~2.0914</math></td>
  <td align="center">1</td>
<td align="center"><math>~1</math></td>
  <td align="right">99.997<math>\pi</math></td>
<td align="center"><math>~1.73281</math></td>
  <td align="right">115.064</td>
  <td align="center">---</td>
  <td align="right">4.6539</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
</tr>
<tr>
<tr>
  <td align="center">6</td>
<td align="center"><math>~(5)</math></td>
  <td align="right">3.1415</td>
<td align="center"><math>~2.7302</math></td>
  <td align="center">2</td>
<td align="center"><math>~1</math></td>
  <td align="right">91284<math>\pi</math></td>
<td align="center"><math>~2.87493</math></td>
  <td align="right">170.472<math>e^{2\pi}</math></td>
  <td align="center">---</td>
  <td align="right">3.1418</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
</tr>
</tr>
<tr>
<tr>
  <td align="left" colspan="11">
<td align="center">&hellip;</td>
<sup>&dagger;</sup>See an [[User:Tohline/SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar.27s_U_and_V_Functions|accompanying discussion]] for definitions of the functions, <math>~V_{1E}</math> and <math>~V_{5F} \, .</math>
<td align="center"><math>~\pi</math></td>
<td align="center"><math>~1</math></td>
<td align="center"><math>~\pi</math></td>
</tr>
</table>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
In the figure, the green curve traces the segment of the E-Type Lane-Emden function for the <math>~n_c = 1</math> polytropic core, <math>~\theta_{1E}</math> &#8212; that is, the segment of the "sinc" function &#8212; that covers the range <math>~[0,\pi]</math> for the interface parameter, <math>~\xi_i</math>;  the blue curve traces the segment of the F-Type Lane-Emden function for the <math>~n_e = 5</math> polytropic envelope, <math>~\phi_{5F}</math>, that covers this same range <math>~[0,\pi]</math>, but for the interface parameter, <math>~\Delta_i</math>.  (This blue curve is also displayed, and its properties described in more depth, in a [[User:Tohline/SSC/Structure/Polytropes#Example_Interval|separate discussion devoted to Srivastava's Lane-Emden function]] for <math>~n=5</math> polytropes.)  This entire displayed parameter range, <math>~[0,\pi]</math>, is associated with physically viable interface values for the core (green curve), but for the envelope (blue curve), only a subset of this range can be associated with physically viable interface values &#8212; namely, the range <math>~[\tan^{-1}(1+2^{1/3}),\pi]</math> over which <math>~\phi_{5F}</math> is positive but its slope is either negative or zero.
The purple dashed line segment labeled "(2)" in the above "Interface Mapping" figure intersects the green curve at the value, <math>~\xi_i = 0.8154</math>, and it intersects the blue curve at the interface value, <math>~\Delta_i = 1.24287</math>.  [These two numerical values are listed  in the row labeled "Model (2)" of the table that has been placed immediately to the right of the figure.]  This is intended to illustrate that an equilibrium bipolytropic configuration with <math>~(n_c, n_e) = (1, 5)</math> can be constructed by truncating the <math>~n_c = 1</math> core at a radius, <math>~\xi_i = 0.8154</math>, and matching it to an <math>~n_e = 5</math> envelope whose innermost radius is associated with the parameter, <math>~\Delta_i = 1.24287</math>.  (The corresponding, analytically determined values of the coefficients <math>~A_0</math> and <math>~B_0</math> dictate the manner in which the blue curve should be stretched both horizontally and vertically to complete a smooth attachment of the envelope to the core.)  In an analogous fashion, the dashed line segments labeled "(3)", "(4)", and "(5)" &#8212; and their corresponding coordinate values as listed in the accompanying table &#8212; illustrate how three additional equilibrium bipolytropic configurations with progressively larger cores can be constructed with the proper interface mapping.  We have specifically chosen to illustrate the interface-mapping of these four bipolytropic models because their cores are truncated at the same values of <math>~\xi_i</math> that [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] used for the models numbered "(2)", "(3)", "(4)", and "(5)" in his Table 3.  (See the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Murphy.27s_Example_Model_Characteristics|much more in-depth discussion, below]]; and note that the edge of the core is labeled by the parameter, <math>~\zeta_J</math>, rather than by <math>~\xi_i</math>, in Murphy's Table 3.)
For clarity we note that, the analytically determined values of <math>~\Delta_i</math> that are listed in the table that sits to the right of the above "Interface Mapping" figure have been shifted in phase by <math>~m\pi</math>, where the relevant value of the integer, <math>~m</math>, is also listed in the table.  In the figure, the two ''purple'' dashed-line segments are associated with models for which <math>~m=0</math>, while the two ''orange'' dashed-line segments are associated with models for which <math>~m=1</math>.
In his Table 3, [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] also tabulates the characteristics of two additional models:  The core of his "model 1" is truncated at a radius much less than <math>~\pi</math>, namely, <math>~\xi_i = 0.032678</math>; as a result, the bipolytropic configuration has a very small core and its structure is almost entirely that of an <math>~n = 5</math> polytrope.  At the other extreme, the core of his "model 6" is truncated at a radius that is almost, but not quite, equal to <math>~\pi</math>, namely, <math>~\xi_i = 3.1415</math>; the resulting bipolytropic configuration has a very tiny envelope and its structure is almost entirely an <math>~n = 1</math> polytrope.  Our analytic solution permits us to set <math>~\xi_i</math> to either of the two limiting values, <math>~0</math> or <math>~\pi</math>, and to show that, in these limits, <math>~\Delta_i</math> exactly equals, respectively, <math>~\tan^{-1}(1 + 2^{1/3})</math> and <math>~\pi</math>.  The mapping of the first of these two limits is illustrated by the red dashed line segment in the above figure; the second limit is illustrated simply by the intersection of the two curves at the coordinate location, <math>~\pi</math>.  It appears, therefore, that in these two limits Murphy's bipolytrope can be used, respectively, to define the structure of an isolated <math>~n=5</math> polytrope or an isolated <math>~n=1</math> polytrope.
'''<font color="red">ASIDE:</font>'''  It has previously been thought that Srivastava's Lane-Emden function, <math>~\phi_{5F}</math>, cannot be used on its own to define the structure of an ''isolated'' <math>~n=5</math> polytrope because the function's amplitude grows without bound and oscillates more and more rapidly between positive and negative values as the governing radial coordinate gets smaller and smaller.  From the results presented here, it now appears as though an isolated <math>~n=5</math> polytrope of this type ''can'' be constructed by letting the interface parameter <math>~\xi_i \rightarrow 0</math> &#8212; and, hence the alternative parameter <math>~\Delta_i \rightarrow \tan^{-1}(1+2^{1/3})</math> &#8212; in Murphy's bipolytrope.  It should be interesting to determine the values of the coefficients, <math>~A_0</math> and <math>~B_0</math>, that arise in this limit, and to examine in detail the structure of the complete <math>~n=5</math> model that results.  Does its radial density profile resemble &#8212; or, perhaps, exactly match &#8212; the radial density profile of the well-known ''isolated'' <math>~n=5</math> ploytrope?  Does the new structure have a finite radius as well as a finite central density?  In this context it is worth noting that, in two separate papers &#8212;  [http://adsabs.harvard.edu/abs/1980PASAu...4...37M Murphy (1980a)] and [http://adsabs.harvard.edu/abs/1981PASAu...4..205M Murphy (1981)] &#8212; Murphy has constructed and discussed the physical characteristics of equilibrium models that obey the <math>~n=5</math> polytropic equation of state all the way from the center to a surface which is of finite radius.  He accomplishes this by piecing together a core that is defined by the familiar, analytically specified, <math>~\phi_{5E}</math> Lane-Emden function and an envelope that is defined by Srivastava's <math>~\phi_{5F}</math> Lane-Emden function.
===Murphy's Example Model Characteristics===
As in an [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#MurphyF2Function|earlier figure]], the following log-log plot shows what value(s) of <math>~(A_0\eta)_\mathrm{root} </math> (vertical axis) result from a wide range of interface coordinate specifications, <math>~\xi_i</math> (horizontal axis), under the simplifying assumption that <math>~(\mu_e/\mu_c) = 1</math>.  The solid blue diamonds present exactly the same information as in the previous figure, that is, they identify roots resulting from setting <math>~m=0</math> in the phase-shift specification.  But here, focusing on a regime that highlights larger roots of the first constraint equation, solid black squares identify roots for which <math>~m=+1</math>, and solid maroon triangles identify roots for which <math>~m=+2</math>. 
In this version of the figure, colored dashed-line segments have been drawn to enclose two regions of parameter space in which solutions of the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#KeyInterfaceRelation|key nonlinear interface relation]] produce a [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Physically_Viable_Parameter_Values|physically viable model]] &#8212; that is, <math>~0 \le \xi_i \le \pi</math> and <math>~\tan^{-1}(1+2^{1/3}) \le \Delta_i \le \pi</math>.  The portion of the red-outlined rectangular region that lies to the left of the vertical, <math>~\xi_i = \xi_\mathrm{trans}</math> line accommodates some models that arise from specifying an <math>~m=0</math> phase shift  and the portion that lies to the right of this demarcation line accommodates some models that arise from specifying an <math>~m=+1</math> phase shift.  Similarly, the portion of the green-outlined rectangular region that lies to the left of the vertical, <math>~\xi_i = \xi_\mathrm{trans}</math> line accommodates additional models that arise from specifying an <math>~m=+1</math> phase shift  and the portion that lies to the right of this demarcation line accommodates models that arise from specifying an <math>~m=+2</math> phase shift.  Given the inherent periodicity of Srivastava's <math>~\phi_{5F}</math> function, it should be clear that we could equally well have identified other rectangular regions at both higher <math>~(m \ge 2)</math> and lower <math>~(m \le 0)</math> values of <math>~(A_0 \eta)_\mathrm{root}</math> in which additional physically viable models exist. It should also be clear that, because the figure presents a log-log plot and the lower limit of viable parameter values for the horizontal axis parameter, <math>~\xi_i</math>, is zero, the identified rectangular regions should be viewed as extending to the left, all the way to minus infinity.




<div align="center">
<div align="center">
<table border="1" cellpadding="0" align="center">
<table border="1" cellpadding="5">
<tr>
  <td align="center" colspan="2">
Examination of Table 3 Data Presented in 
[http://adsabs.harvard.edu/abs/1983PASAu...5..175M J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)]
  </td>
</tr>
<tr>
<tr>
   <td align="center">
   <td align="center">
Physically Viable Parameter Regime
[[File:MurphyPhysicallyViableTable3d.png|center|500px|Table 3 Model Characteristics by Murphy (1983)]]
   </td>
   </td>
</tr>
  <td align="center">
<tr><td align="center">
<table border="0" cellpadding="5">
[[File:MurphyPhysicallyViable.png|center|400px|Physically Viable Parameter Region]]
</td></tr>
</table>
</div>
 
===Second Constraint===
 
====Obtained from Third Interface Condition====
Now, our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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>
<tr>
   <td align="right">
   <td align="center" colspan="2" width="100%">
<math>~\frac{\eta_i}{\xi_i}</math>
[[File:CommentButton02.png|center|100px|Comment by J. E. Tohline on 29 April 2015:  As is argued more fully below, it appears as though the two parameter values drawn directly from Murphy's Table 3 that are shown here in red contain typographical errors. ]]
  </td>
Murphy's Table 3 Data
  <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>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="center" width="50%">
&nbsp;
<math>~\zeta_J \rightarrow \xi_i</math><p></p>
----
  </td>
  <td align="center" width="50%">
<math>~\xi_J \rightarrow (A_0\eta)_\mathrm{root}</math><p></p>
----
   </td>
   </td>
   <td align="center">
</tr>
<math>~=</math>
<tr>
   <td align="right">
<font face="Courier">3.2678E-02</font>
   </td>
   </td>
   <td align="left">
   <td align="right">
<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>
<font face="Courier" color="red">10.0164</font>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Rightarrow ~~~~ B^2</math>
<font face="Courier">0.8154</font>
   </td>
   </td>
   <td align="center">
   <td align="right">
<math>~=</math>
<font face="Courier">12.0083</font>
  </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>
   </td>
</tr>
</tr>
 
<tr>
</table>
</div>
 
====Obtained from Chandrasekhar's U-constraint====
We shall now demonstrate that the same expression for the scaling coefficient, <math>~B</math>, can alternatively be obtained from Chandrasekhar's "U" constraint after taking into account the result that already has been obtained from the "V" constraint.  As [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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">
   <td align="right">
<math>~\frac{\xi_i \theta_i}{(-d\theta/d\xi)_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>
<font face="Courier">1.6598</font>
   </td>
   </td>
   <td align="center">
   <td align="right">
<math>~=</math>
<font face="Courier">20.4312</font>
  </td>
  <td align="left">
<math>~\frac{\eta_i \phi_i^{5}}{(-d\phi/d\eta)_i} \, ,</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
which, in the particular case being examined here, becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{2\xi_i^2}{3\kappa_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)^2 </math>
<font face="Courier" color="red">2.6914</font>
   </td>
   </td>
   <td align="center">
   <td align="right">
<math>~=</math>
<font face="Courier">33.0249</font>
  </td>
</tr>
<tr>
  <td align="right">
<font face="Courier">2.7302</font>
   </td>
   </td>
   <td align="left">
   <td align="right">
<math>~\biggl( U_\mathrm{5F} \biggr)_i </math>
<font face="Courier">100<math>\pi</math></font>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<font face="Courier">3.1415</font>
   </td>
   </td>
   <td align="center">
   <td align="right">
<math>~=</math>
<font face="Courier">91275<math>\pi</math></font>
   </td>
   </td>
  <td align="left">
</tr>
<math>~
</table>
\frac{2B^{-4} \sin^4\Delta_i}{(3-2\sin^2\Delta_i)(3 - 2\sin^2\Delta_i - 3\cot\Delta_i)} \, .
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,352: Line 1,459:
</div>
</div>


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


In his Table 3, [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy (1983)] presents six pairs of "coordinate" roots, <math>~(\zeta_J, \xi_J)</math>, that not only satisfy the first constraint equation but that also, in his interpretation, identify physically viable model solutions.  The values of all six of these coordinate pairs have been re-listed immediately to the right of the log-log plot, in a separate panel of the above figure.  After mapping Murphy's coordinate notation to ours &#8212; specifically invoking the mapping <math>~(\zeta_J, \xi_J) ~\rightarrow~ (\xi_i, A_0\eta_\mathrm{root})</math> &#8212; the locations of all six models have been marked by orange circles on the plot.  Five of Murphy's models fall directly onto our generated curves and in regions marked as physically viable, giving us confidence that our analysis matches his in considerable detail.  As is explained more fully, below, we strongly suspect that the single Murphy model that does not lie on one of our generated curves &#8212; "Model 4" which, as published in Murphy's Table 3, has coordinates <math>~(\zeta_J, \xi_J) = (2.6914, 33.0249)</math> &#8212; is displaced, perhaps because of transcription errors in the published table.
As a supplement to the preceding two figures, we now display a table of our ''analytically determined'' model characteristics that is patterned after and largely matches [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy's (1983)] Table 3. The tabulated values have been derived assuming <math>~\mu_e/\mu_c = 1</math> throughout, that is, assuming that the mean molecular weights of the core and the envelope are identical.  Vertically, the table is divided into three segments:  The top portion contains six rows of numbers that result from our derivations; the middle portion displays six columns of numbers that have been extracted directly from Table 3 of Murphy (1983) and displayed here as image segments; and the bottom portion shows additional parameter values that we have generated in an effort to understand why our "Model 4" does not match the fourth model in Murphy's table.
<div align="center">
<table border="1" cellpadding="5">
<tr>
  <td align="center" colspan="9">
Some Model Characteristics<sup>&dagger;</sup> (assuming <math>\mu_e/\mu_c</math> = 1)
  </td>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="center" rowspan="2">
<math>~(3 - 2\sin^2\Delta_i - 3\cot\Delta_i)</math>
Model
   </td>
   </td>
   <td align="center">
   <td align="center" colspan="2">
<math>~=</math>
Specified
   </td>
   </td>
   <td align="left">
   <td align="center" colspan="6">
<math>~(3-2\sin^2\Delta_i)\kappa_i \, .</math>
Analytically Determined Here
   </td>
   </td>
</tr>
</tr>
</table>
</div>
Hence, Chandrasekhar's U-constraint becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="center">
<math>~\frac{2\xi_i^2}{3\kappa_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)^2 </math>
<math>~\xi_i</math>
  </td>
  <td align="center">
<math>~m</math>
  </td>
  <td align="center">
<math>~\Delta_i</math>
  </td>
  <td align="center">
<math>~(A_0\eta)_\mathrm{root}</math>
  </td>
  <td align="center">
<math>~\frac{A_0}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) </math>
  </td>
  <td align="center">
<math>~B_0</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\frac{B_0}{3^{1/4}} \biggl(\frac{\mu_e}{\mu_c}\biggr)^{1/2} \theta_i^{-1} </math>
   </td>
   </td>
  <td align="left">
   <td align="center">
<math>~
<math>~\xi_s</math>
\frac{2B^{-4} \sin^4\Delta_i}{(3-2\sin^2\Delta_i)^2 \kappa_i}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ B^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>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="right">
  <td align="center">1</td>
<math>~\Rightarrow ~~~~ B^2</math>
  <td align="right">3.2678E-02</td>
   </td>
  <td align="center">0</td>
   <td align="center">
  <td align="right">1.1544</td>
<math>~=</math>
  <td align="right">10.0614</td>
   </td>
  <td align="right">307.894</td>
   <td align="left">
  <td align="center">5.77929</td>
<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 align="center">4.39209</td>
   </td>
  <td align="right">1.7392</td>
</tr>
</tr>
</table>
 
</div>
<tr>
which, as predicted, is identical to what we learned from the third interface condition, alone.
  <td align="center">2</td>
  <td align="right">0.8154</td>
  <td align="center">0</td>
  <td align="right">1.2429</td>
  <td align="right">12.0101</td>
  <td align="right">14.729</td>
  <td align="center">1.25567</td>
  <td align="center">1.06865</td>
  <td align="right">36.356</td>
</tr>
 
<tr>
  <td align="center">3</td>
  <td align="right">1.6598</td>
  <td align="center">0</td>
  <td align="right">1.4918</td>
  <td align="right">19.7585</td>
  <td align="right">11.9041</td>
  <td align="center">1.01206</td>
  <td align="center">1.28146</td>
  <td align="right">44.984</td>
</tr>
 
<tr>
  <td align="center">4</td>
  <td align="right">2.6914</td>
  <td align="center">1</td>
  <td align="right">2.8365</td>
  <td align="right"><font color="red">290.9431</font></td>
  <td align="right"><font color="red">108.1010</font></td>
  <td align="center">0.14348</td>
  <td align="center"><font color="red">0.67430</font></td>
  <td align="right"><font color="red">4.9536</font></td>
</tr>
 
<tr>
  <td align="center">5</td>
  <td align="right">2.7302</td>
  <td align="center">1</td>
  <td align="right">2.8749</td>
  <td align="right">100.00<math>\pi</math></td>
  <td align="right">115.065</td>
  <td align="center">0.12408</td>
  <td align="center">0.64370</td>
  <td align="right">4.6539</td>
</tr>
 
<tr>
  <td align="center">6</td>
  <td align="right">3.1415</td>
  <td align="center">2</td>
  <td align="right">6.283141</td>
  <td align="right">91268<math>\pi</math></td>
  <td align="right">170.44223<math>e^{2\pi}</math></td>
  <td align="center">1.8966E-05</td>
  <td align="center">0.48862</td>
  <td align="right">3.141778</td>
</tr>
<tr>
  <td align="center" colspan="9">[[File:CommentButton02.png|right|50px|Comment by J. E. Tohline on 11 May 2015:
The heading on this last column of numbers extracted from Murphy's (1983) paper has a subscript "J" whereas the subscript should be "R" to indicate the configuration's total radius.]]<p></p>
Corresponding Values Extracted Directly from [http://adsabs.harvard.edu/abs/1983PASAu...5..175M Murphy's (1983)] Table 3
  </td>
</tr>
 
<tr>
<td align="center" colspan="9">[[Image:AAAwaiting01.png|600px|Columns extracted directly from Murphy (1983)]]</td>
<!--
  <td align="center">[[File:Table3CompleteCol1.png|center|50px|Column 1 from Table 3 of Murphy (1983)]]</td>
  <td align="right" colspan="2">[[File:Table3CompleteCol2.png|center|123px|Column 3 from Table 3 of Murphy (1983)]]</td>
  <td align="center">&nbsp;</td>
  <td align="right">[[File:Table3CompleteCol3.png|center|67px|Column 3 from Table 3 of Murphy (1983)]]</td>
  <td align="right">[[File:Table3CompleteCol4.png|center|117px|Column 4 from Table 3 of Murphy (1983)]]</td>
  <td align="center">&nbsp;</td>
  <td align="right">[[File:Table3CompleteCol5.png|center|71px|Column 5 from Table 3 of Murphy (1983)]]</td>
  <td align="right">[[File:Table3CompleteCol6.png|center|83px|Column 6 from Table 3 of Murphy (1983)]]</td>
-->
</tr>
<tr>
  <td align="center" colspan="9">
Attempting to Understand Murphy's "Model 4"
  </td>
</tr>
 
<tr>
  <td align="center">4b</td>
  <td align="right">0.048681</td>
  <td align="center">0</td>
  <td align="right">1.1545342</td>
  <td align="right">10.065043</td>
  <td align="right">206.75505</td>
  <td align="center">4.7358862</td>
  <td align="center" bgcolor="lightgreen">3.5999</td>
  <td align="right">2.589981</td>
</tr>
 
<tr>
  <td align="center">4c</td>
  <td align="right">1.14988</td>
  <td align="center">0</td>
  <td align="right">1.3234244</td>
  <td align="right">14.109506</td>
  <td align="right" bgcolor="lightgreen">12.2704</td>
  <td align="center">1.12444163</td>
  <td align="center">1.07640071</td>
  <td align="right" bgcolor="lightgreen">43.6409</td>
</tr>
 
<tr>
  <td align="center">4d</td>
  <td align="right">1.77027</td>
  <td align="center">0</td>
  <td align="right">1.5391618</td>
  <td align="right">21.721958</td>
  <td align="right" bgcolor="lightgreen">12.2704</td>
  <td align="center">0.9876639</td>
  <td align="center">1.355397</td>
  <td align="right" bgcolor="lightgreen">43.6409</td>
</tr>
 
<tr>
  <td align="center">4e</td>
  <td align="right">2.10961</td>
  <td align="center">1</td>
  <td align="right">1.7486314</td>
  <td align="right" bgcolor="lightgreen">33.0249</td>
  <td align="right">15.6545186</td>
  <td align="center">0.8651536</td>
  <td align="center">1.6157228</td>
  <td align="right">34.2068427</td>
</tr>
 
<tr>
  <td align="left" colspan="9">
<sup>&dagger;</sup>See an [[User:Tohline/SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar.27s_U_and_V_Functions|accompanying discussion]] for a complementary critique of Murphy's Table 3 models including an evlauation of the functions, <math>~U_{1E}</math>, <math>~V_{1E}</math>, <math>~U_{5F}</math>, and <math>~V_{5F} \, .</math>
  </td>
</tr>
 
</table>
</div>
 
 
Our six models (top portion of the above table) have been constructed by, first, specifying the location of the outer edge of the core, <math>~\xi_i</math> along with a corresponding phase shift parameter, <math>~m</math>, to match Murphy's chosen core-interface values &#8212; in his terminology, <math>~\zeta_J</math> (middle portion of the above table).  Once the value of <math>~\xi_i</math> has been specified for each model, the value of the parameter, <math>~p</math>, the root of the cubic equation, <math>~y_\mathrm{root}</math>, and all of the parameter values listed in the upper portion of the table were determined using the parameter relations derived above and/or detailed in the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Parameter_Values|table below]].  For example,
 
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~\Delta_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\tan^{-1}(y_\mathrm{root}) + m\pi \, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~(A_0\eta)_\mathrm{root}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~e^{2\Delta_i} \, ;</math>
  </td>
</tr>
</table>
</div>
and,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\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>
 
Comparing our tabulated numbers (upper portion of the table) with the numbers published by Murphy (middle portion of the table) adds strong support to the statement made in connection with the above figure, that is, that our analysis matches Murphy's in considerable detail.  But such a comparison also highlights the mismatch between Murphy's "Model 4" and ours &#8212; to aid in this particular comparison, the relevant numbers have been displayed in a red font in the upper portion of our table.  According to our analysis, a model with <math>~\xi_i = 2.6914</math> does not display characteristics that, in any fashion, match Murphy's model with the same core-interface value.
 
We initially suspected that the "Model 4" mismatch resulted from a typographical or transcription error in the value of <math>~\zeta_J</math> &#8212; in our terminology, <math>~\xi_i</math> &#8212; that appears in Murphy's published Table 3.  So we tried a variety of values of <math>~\xi_i</math>, looking for a set of generated characteristics that matched Murphy's published Model 4 characteristics.  We were unsuccessful in this effort.  The most tantalizing results from this search are displayed as cells with light-green backgrounds in the bottom segment of the above table:  (Model 4b) Setting <math>~\xi_i = 0.048681~(m = 0)</math> gives <math>~B=3.5999</math>, which matches Murphy's "Model 4" value of <math>~B</math>; (Model 4c) setting <math>~\xi_i = 1.14988~(m = 0)</math> gives values of <math>~A</math> and <math>~\zeta_J (\mathrm{i.e.,}~\xi_s)</math> that both match Murphy's values; (Model 4d) setting <math>~\xi_i = 1.77027~(m = 0)</math> also gives values of <math>~A</math> and <math>~\zeta_J</math> that both match Murphy's values; and (Model 4e) setting <math>~\xi_i = 2.10961~(m = 1)</math> gives a value of <math>~(A_0\eta)_\mathrm{root}</math> that matches Murphy's value of <math>~\xi_J</math>.  As we have been unable to find a value of <math>~\xi_i</math> whose associated parameter values match Murphy's "Model 4" values across the board, we remain at a loss to explain the mismatch between our "Model 4" and Murphy's.
 
==Step 8:  Envelope's Physical Profile==
 
Combining the set of expressions provided in the right-most column of the [[User:Tohline/SSC/Structure/BiPolytropes#TableSetup|accompanying bipolytrope "Table 1:  Setup"]] with the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Step_5:_Interface_Conditions|interface conditions already established, above]], we can now specify how the density, pressure, and integral mass vary throughout the envelope, that is, over the radial extent <math>\xi_i \le \xi \le \xi_s</math>.
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>~\rho</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\rho_e \phi^{5}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\rho_0 \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \phi^{5} \, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~P</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~K_e \rho_e^{6/5} \phi^{6}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~K_c \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  \biggr]^{6/5} \phi^{6}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~K_c \rho_0^2 \theta^{2}_i \phi^{6} \, ;</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{3K_e}{2\pi G} \biggr]^{3/2} \rho_e^{-1/5} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{3^3 \cdot (2^2\pi)^2}{2^3\pi^3 G^3} \biggr]^{1/2} \biggl[K_c^5 \rho_0^{4}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-6} \theta^{4}_i\biggr]^{3/10} 
\biggl[ \rho_0 \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i  \biggr]^{-1/5} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{2\cdot 3^3 }{\pi G^3} \biggr]^{1/2} K_c^{3/2} \rho_0 \theta_i \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2}
\biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr) \, .</math>
  </td>
</tr>
</table>
</div>
 
=Examples=
 
==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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#Step_4:__Throughout_the_core_.28.29|step &#35;4]] and [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1_5#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>
 
=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 [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#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=
=Related Discussions=

Latest revision as of 21:48, 11 March 2019

BiPolytrope with <math>n_c = 1</math> and <math>n_e=5</math>

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |
Comment by J. E. Tohline on 12 April 2015: I became aware of the published discussions of this system by Murphy (1983) and Murphy & 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) §2.8.1 discussion of composite polytropes.

Here we construct a 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 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 J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175), and attributes of their physical structure were further discussed by J. O. Murphy & R. Fiedler (1985a, Proc. Astr. Soc. of Australia, 6, 219); additional, closely related references are given below. In the discussion that follows, we will be heavily referencing Murphy's (1983) work.

Steps 2 & 3

Based on the discussion presented elsewhere of the structure of an isolated <math>~n=1</math> polytrope, the core of this bipolytrope will have the following properties:

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

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 (<math>0 \le \xi \le \xi_i</math>)

Specify: <math>~K_c</math> and <math>~\rho_0 ~\Rightarrow</math>

 

<math>~\rho</math>

<math>~=</math>

<math>\rho_0 \theta^{n_c}</math>

<math>~=</math>

<math>~\rho_0 \biggl( \frac{\sin\xi}{\xi} \biggr)</math>

<math>~P</math>

<math>~=</math>

<math>~K_c \rho_0^{1+1/n_c} \theta^{n_c + 1}</math>

<math>~=</math>

<math>~K_c \rho_0^{2} \biggl( \frac{\sin\xi}{\xi}\biggr)^{2}</math>

<math>~r</math>

<math>~=</math>

<math>~\biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{1/2} \rho_0^{(1-n_c)/(2n_c)} \xi</math>

<math>~=</math>

<math>~\biggl[ \frac{K_c}{2\pi G} \biggr]^{1/2} \xi</math>

<math>~M_r</math>

<math>~=</math>

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

<math>~=</math>

<math>~4\pi \biggl[ \frac{K_c}{2\pi G} \biggr]^{3/2} \rho_0 \biggl[\sin\xi - \xi \cos\xi \biggr]</math>

Step 5: Interface Conditions

 

Setting <math>~n_c=1</math> and <math>~n_e=5~~~~~\Rightarrow</math>

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

<math>~=</math>

<math>~\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{n_c}_i \phi_i^{-n_e}</math>

<math>~=</math>

<math>~\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \phi_i^{-5}</math>

<math>~\biggl( \frac{K_e}{K_c} \biggr) </math>

<math>~=</math>

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

<math>~=</math>

<math>~\biggl[\rho_0^{4}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-6} \theta^{4}_i\biggr]^{1/5}</math>

<math>~\frac{\eta_i}{\xi_i}</math>

<math>~=</math>

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

<math>~=</math>

<math>~\biggl( \frac{1}{3} \biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \phi_i^{-2}</math>

<math>~\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>

<math>~=</math>

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

<math>~=</math>

<math>~\biggl( \frac{1}{3} \biggr)^{1/2} \theta_i^{- 1} \biggl( \frac{d\theta}{d\xi} \biggr)_i \phi_i^3</math>


Alternative: In our introductory description of how to build a bipolytropic structure, we pointed out that, instead of employing these last two fitting conditions, Chandrasekhar [C67] 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,

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

and the product of the <math>3^\mathrm{rd}</math> and <math>4^\mathrm{th}</math> expressions, which in the present case generates,

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

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 an accompanying discussion, Murphy (1983) 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

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 Murphy (1983) and of Murphy & Fiedler (1985a), we will adopt for the envelope's structure the F-Type solution of the <math>~n=5</math> Lane-Emden function discovered by S. Srivastava (1968, ApJ, 136, 680) and described in an accompanying discussion, namely,

<math>~\phi</math>

<math>~=</math>

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

 

<math>~=</math>

<math>~\frac{B_0^{-1}\sin\Delta}{\eta^{1/2}(3-2\sin^2\Delta)^{1/2}} \, ,</math>

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 — the values of both will be determined presently from the interface conditions — and we have introduced the notation,

<math>~\Delta \equiv \ln(A_0\eta)^{1/2} = \frac{1}{2} (\ln A_0 + \ln\eta) \, .</math>

The first derivative of Srivastava's function is,

<math>~\frac{d\phi}{d\eta}</math>

<math>~=</math>

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

As has been explained in the context of our 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 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,

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

<math>~=</math>

<math>~\biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{1/2} \rho_e^{(1-n_e)/(2n_e)} \eta</math>

<math>~\Rightarrow ~~~~\biggl( \frac{K_c}{3K_e} \biggr)\rho_e^{4/5} </math>

<math>~=</math>

<math>~\biggl(\frac{\eta}{\xi}\biggr)^2 \, .</math>

Referring back to the 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,

<math>~\biggl(\frac{\eta}{\xi}\biggr)^2</math>

<math>~=</math>

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

<math>~\Rightarrow ~~~~ \frac{\eta}{\xi}</math>

<math>~=</math>

<math>~\biggl(\frac{1}{3}\biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2} \, .</math>

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 — not just at the interface — 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 — see also our accompanying discussion for elaboration on Murphy's (1983) "V5F" and "V1E" function notations — one fitting condition at the interface is,

<math>~- \frac{2\xi_i }{ 3\theta_i } \biggl( \frac{d\theta}{d\xi} \biggr)_i \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>

<math>~=</math>

<math>~2\biggl[ \frac{\eta (- d\phi/d\eta)}{\phi} \biggr]_i </math>

 

<math>~=</math>

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

 

<math>~=</math>

<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, .</math>

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,

<math>~ \kappa_i \equiv - \frac{2\theta_i^' \xi_i}{3\theta_i} \biggl( \frac{\mu_e}{\mu_c} \biggr) \, ,</math>

we will recast the first constraint into, what will henceforth be referred to as, the

Key Nonlinear Interface Relation

<math>~\kappa_i</math>

<math>~=</math>

<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, .</math>

In a separate subsection of this chapter, below, we present a closed-form analytic solution, <math>~\Delta_i(\kappa_i)</math>, to this nonlinear equation.


Second Constraint

Obtained from Third Interface Condition

Our 3rd interface condition, as detailed above, states that,

<math>~\frac{\eta_i}{\xi_i}</math>

<math>~=</math>

<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \phi_i^{-2} \, .</math>

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 3rd interface condition alone,

<math>~\eta_i</math>

<math>~=</math>

<math>~3^{-1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \xi_i \, ;</math>

and, second, from the definition of Srivastava's function, <math>~\phi</math>, we deduce that the overall scaling parameter is,

<math>~B_0^2</math>

<math>~=</math>

<math>~\frac{\sin^2\Delta_i}{\eta_i(3-2\sin^2\Delta_i)} </math>

 

<math>~=</math>

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

Notice that, after the solution, <math>~\Delta_i(\kappa_i)</math>, of the key nonlinear interface relation has been determined, the first of these two relations also permits us to write,

<math>~A_0\eta_i</math>

<math>~=</math>

<math>~e^{2\Delta_i}</math>

<math>~\Rightarrow ~~~~ A_0</math>

<math>~=</math>

<math>~3^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \xi_i^{-1} e^{2\Delta_i} \, .</math>

Throughout the envelope, therefore, the angle,

<math>~\Delta</math>

<math>~=</math>

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


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 described above, the U-constraint is an alternative interface condition that may be written as,

<math>~\frac{\xi_i \theta_i}{(-d\theta/d\xi)_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>

<math>~=</math>

<math>~\frac{\eta_i \phi_i^{5}}{(-d\phi/d\eta)_i} \, ,</math>

which, in the particular case being examined here, becomes — again, see our accompanying discussion for elaboration on Murphy's (1983) "U5F" and "U1E" function notations —

<math>~\frac{2\xi_i^2}{3\kappa_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)^2 </math>

<math>~=</math>

<math>~\biggl( U_\mathrm{5F} \biggr)_i </math>

 

<math>~=</math>

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

Now, from our discussion, above, of the first constraint, we know that,

<math>~(3 - 2\sin^2\Delta_i - 3\cot\Delta_i)</math>

<math>~=</math>

<math>~(3-2\sin^2\Delta_i)\kappa_i \, .</math>

Hence, Chandrasekhar's U-constraint becomes,

<math>~\frac{2\xi_i^2}{3\kappa_i} \biggl( \frac{\mu_e}{\mu_c}\biggr)^2 </math>

<math>~=</math>

<math>~ \frac{2B_0^{-4} \sin^4\Delta_i}{(3-2\sin^2\Delta_i)^2 \kappa_i} </math>

<math>~\Rightarrow ~~~~ B_0^4 </math>

<math>~=</math>

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

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

<math>~=</math>

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

which, as predicted, is identical to what we learned from the third interface condition, alone.

Comment on Murphy's Scalings

Murphy's (1983) 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,

<math>~A = \frac{\xi_J}{\zeta_J} \, ,</math>

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,

<math>~A = \frac{A_0 \eta_\mathrm{root}}{\xi_i} \, .</math>

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,

<math>~A = \frac{A_0 \eta_i}{\xi_i} \, .</math>

Recalling that <math>~\phi_i = 1</math>, we know from the interface conditions detailed above that,

<math>~\frac{\eta_i}{\xi_i} = \frac{1}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) \, .</math>

Hence Murphy's homology factor, <math>~A</math>, is related to our homology factor, <math>~A_0</math> via the expression,

<math>~A = \frac{A_0}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) \, .</math>

It is usually the value of this quantity, rather than simply our derived value of <math>~A_0</math>, that is tabulated below — both here and here — as we make quantitative comparisons between the characteristics of our derived models and those published by Murphy (1983) and by Murphy & Fiedler (1985a).


In the lefthand column of the fourth page of his paper, Murphy (1983) 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,

<math>~B</math>

<math>~=</math>

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

 

<math>~=</math>

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

Switching to our terminology, that is, setting,

<math>\ln\sqrt{A\zeta_J} \rightarrow \Delta_i</math>     and, as before,      <math>~\zeta_J \rightarrow \xi_i \, ,</math>

gives,

<math>~B</math>

<math>~=</math>

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

 

<math>~=</math>

<math>~\theta_i^{-1} \biggl[ \xi_i^{-1} \biggl(\frac{3}{\sin^2\Delta_i} - 2\biggr)^{-1} \biggr]^{1/2} \, .</math>

Hence, in terms of the definition of our scaling coefficient, <math>~B_0</math>, derived above, we have,

<math>~B</math>

<math>~=</math>

<math>~\frac{B_0}{3^{1/4}} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1/2} \theta_i^{-1} \, .</math>

As we make quantitative comparisons between the characteristics of our derived models and those published by Murphy (1983) and by Murphy & 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 — and, along with it, the envelope's density — drops to zero when,

<math>\Delta = \Delta_s \equiv \pi \, ,</math>

we know that the radius, <math>~\xi_s</math>, of the bipolytropic configuration is given by the expression,

<math>~\pi</math>

<math>~=</math>

<math>~ \Delta_i + \ln\biggl( \frac{\xi_s}{\xi_i} \biggr)^{1/2} </math>

<math>~\Rightarrow ~~~~ \ln\biggl( \frac{\xi_s}{\xi_i} \biggr)^{1/2}</math>

<math>~=</math>

<math>~ \pi - \Delta_i </math>

<math>~\Rightarrow ~~~~ \xi_s</math>

<math>~=</math>

<math>~ \xi_i e^{2(\pi - \Delta_i )} \, . </math>

In terms of the natural radial coordinate of the envelope, this is,

<math>~\eta_s \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-1}</math>

<math>~=</math>

<math>~3^{-1/2} \xi_i e^{2(\pi - \Delta_i)} \, .</math>

Analytic Solution of Key Interface Relation

Returning to our previously derived,

Key Nonlinear Interface Relation

<math>~\kappa_i</math>

<math>~=</math>

<math>~\frac{3 - 2\sin^2\Delta_i -3\cot\Delta_i}{(3-2\sin^2\Delta_i)} \, ,</math>

and, as in our separate discussion of the properties of Srivastava's function, adopting the shorthand notation,

<math>~y_i \equiv \tan\Delta_i \, ,</math>

this key interface condition becomes,

<math>~\kappa_i</math>

<math>~=</math>

<math>~ \frac{3 - 2y_i^2(1+y_i^2)^{-1}- 3y_i^{-1} }{3-2y_i^2(1+y_i^2)^{-1}} </math>

 

<math>~=</math>

<math>~ \frac{3y_i(1+y_i^2)- 2y_i^3 -3(1+y_i^2)}{3y_i(1+y_i^2)-2y_i^3} </math>

 

<math>~=</math>

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

<math>~\Rightarrow~~~~ \kappa_i(3y_i+y_i^3)</math>

<math>~=</math>

<math>~ y_i^3 -3y_i^2 + 3y_i -3 </math>

<math>~\Rightarrow~~~~ y_i^3(1-\kappa_i) -3 y_i^2 + 3(1-\kappa_i)y_i -3</math>

<math>~=</math>

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

This is a cubic equation whose solution, <math>~y_\mathrm{root}(\kappa_i)</math>, will also immediately supply the desired interface angle, <math>~\Delta_i</math>, and an interface coordinate root, <math>~(A_0\eta)_\mathrm{root}</math>.

ASIDE: Analytic Solution of Cubic Equation

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

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

a real root is given by the expression,

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

where,

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

In our particular case,

<math>~a =(1-\kappa_i)\, ,</math>      <math>~b =-3\, ,</math>      <math>~c = 3(1-\kappa_i) \, ,</math>      and      <math>~d = - 3 \, .</math>

Hence,

<math>~p = \frac{1}{(1-\kappa_i)} \, ,</math>      <math>~r=+1 \, ,</math>      and      <math>~q = p^3 = \frac{1}{(1-\kappa_i)^3} \, ,</math>

which implies that the real root, <math>~y_\mathrm{root}</math>, is given by the expression,

<math>~y_\mathrm{root}</math>

<math>~=</math>

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

 

<math>~=</math>

<math>~ p + \{p^3 + p^3[1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} + \{p^3 - p^3[1 + (p^{-2}-1)^3]^{1/2}\}^{1/3} </math>

<math>~\Rightarrow ~~~~ p^{-1} y_\mathrm{root}</math>

<math>~=</math>

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

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

In summary, then,

  • Once the location, <math>~\xi_i</math>, of the outer edge of the core has been specified, which determines <math>~\theta_i</math> and <math>~\theta^'_i</math> as well, the value of the parameters, <math>~\kappa_i</math> and <math>~p</math>, are known via the expressions,

<math>~\kappa_i \equiv - \frac{2\theta_i^' \xi_i}{3\theta_i} \biggl( \frac{\mu_e}{\mu_c} \biggr)</math>

<math>~=</math>

<math>~\frac{2}{3} \biggl(1 - \xi_i \cot\xi_i \biggr)\biggl( \frac{\mu_e}{\mu_c} \biggr) \, ,</math>

and,

<math>~p \equiv (1-\kappa_i)^{-1} </math>

<math>~=</math>

<math>~\frac{3}{3-2(\mu_e/\mu_c)(1-\xi_i \cot\xi_i)} \, .</math>

  • The value of <math>~y_\mathrm{root}</math> is determined from the just-derived solution to the governing cubic equation, which then gives the interface value of the envelope parameter (adjusted by a <math>~\pi m</math> phase angle, where <math>~m</math> is an, as yet unspecified, integer),

<math>~\Delta_i = \tan^{-1}(y_\mathrm{root}) +\pi m \, .</math>

  • This, in turn, can be interpreted as a coordinate root — which we shall refer to as <math>~\eta_\mathrm{root}</math> — via the expression,

<math>~(A_0\eta)_\mathrm{root} = e^{2\Delta_i} = e^{2\pi m} \cdot e^{2\tan^{-1}(y_\mathrm{root})} \, .</math>


CAUTION: Solution Behavior When <math>~\kappa_i = 1</math> and <math>~\Delta_i = \pi/2 \, .</math>

In building a sequence of bipolytropic configurations having <math>~(n_c, n_e) = (1, 5)</math>, it will make sense to steadily increase the value of the parameter that marks the edge of the core, <math>~\xi_i</math>, from zero — meaning no mass in the core — to its maximum allowed value, <math>~\pi</math> — meaning no mass in the envelope. (A more complete discussion of physically viable parameter values is presented below.) As <math>~\xi_i</math> steadily increases from zero, for a while the parameter, <math>~\kappa_i</math>, will steadily increase from zero as well. As the value of <math>~\kappa_i</math> crosses through the value of "one," the associated parameter, <math>~p = (1-\kappa_i)^{-1}</math>, as well as the root of the cubic equation, <math>~y_\mathrm{root}</math>, will flip from positive infinity to negative infinity. From the standpoint of the construction of physically realistic models, this does not pose a problem. It simply reflects the fact that, when <math>~\kappa_i = 1</math>, the governing interface angle, <math>~\Delta_i = \pi/2</math> — or, allowing for the aforementioned phase shift, <math>~\Delta_i = \pi/2 + m\pi</math>.


From a practical standpoint, however, it is useful to keep in mind that the root of the governing cubic equation will change abruptly from <math>~+ \infty</math> to <math>~- \infty</math> and its value will be quite sensitive to the choice of <math>~\xi_i</math> in the vicinity of <math>~\kappa_i = 1</math>. Given the definition of the function, <math>~\kappa_i(\xi_i)</math>, this means that the abrupt transition will occur at an interface location, <math>~\xi_\mathrm{trans}</math>, whose value satisfies the condition,

<math>~\frac{\xi_\mathrm{trans}}{\tan(\xi_\mathrm{trans})} </math>

<math>~=</math>

<math>~1-\frac{3}{2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \, .</math>

For the case, <math>~\mu_e/\mu_c = 1 \, ,</math> we have determined empirically (to a dozen significant digits), that,

<math>~\xi_\mathrm{trans} = 1.836597203152 \, .</math>

Just for fun — and, again, for the case, <math>~\mu_e/\mu_c = 1 </math> — we also have determined that <math>~\xi_\mathrm{trans}</math> satisfies the series,

<math>~ 0 = \sum_{n=0}^{\infty} (-1)^{n+1} \biggl(\frac{4n+3}{2n+1} \biggr) \frac{\xi_\mathrm{trans}^{2n}}{(2n)!} \, .</math>

The black, vertical, long-dashed line in the following figure identifies where <math>~\xi_i = \xi_\mathrm{trans}</math>.


Because the function <math>~\tan(\xi_\mathrm{trans})</math> is periodic, other values of <math>~\xi_\mathrm{trans}</math> will also satisfy the condition,

<math>~\frac{\xi_\mathrm{trans}}{\tan(\xi_\mathrm{trans})} </math>

<math>~=</math>

<math>~1-\frac{3}{2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \, .</math>

But the one whose value we have determined here is the only one that falls within the physically viable range of interface location values, <math>~0 \le \xi_i \le \pi </math>.


The following log-log plot shows what value(s) of <math>~(A_0\eta)_\mathrm{root} </math> (vertical axis) result from a wide range of interface coordinate specifications, <math>~\xi_i</math> (horizontal axis), under the simplifying assumption that <math>~(\mu_e/\mu_c) = 1</math>: Solid blue diamonds identify roots resulting from setting <math>~m=0</math> in the phase-shift specification, solid green triangles identify roots for which <math>~m=-1</math>, and solid purple squares identify roots for which <math>~m=-2</math>. As has just been discussed in the context of a cautionary note, as <math>~\xi_i</math> is increased from zero for any given value of <math>~m</math> — that is, as one traverses the plot, moving from the lefthand edge toward the right along a curve of constant color — the solution, <math>~(A_0\eta)_\mathrm{root} </math>, makes an abrupt jump when <math>~\xi_i</math> crosses the value, <math>~\xi_\mathrm{trans}</math>. This identifies the root of the cubic equation for which <math>~\Delta_i = \pi/2 + m\pi</math>. One can continue to move in a smooth, continuous fashion along a single — but multi-colored — solution curve by letting <math>~m \rightarrow (m+1)</math> as <math>~\xi_i</math> crosses <math>~\xi_\mathrm{trans}</math>.


Examination of F2 Function Discussed by J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)

F2 Roots by Murphy (1983)
Murphy's

Roots of

F2 = 0

Function


3.5076E-05

1.8785E-02

4.2993E+00

7.6882E+00

1.0913E+01

1.4101E+01

The red line in the figure shows where the coordinate value along the vertical axis equals the coordinate value along the horizontal axis. Hence the points where this red line intersects the other curves in the figure can be interpreted as identifying solutions to the cubic equation for which,

<math>~(A_0\eta)_\mathrm{root} = \xi_i \, .</math>

The points of intersection that are highlighted in this figure are also what Murphy (1983) refers to as roots of his <math>~F2 = 0</math> function. The values of ten successive coordinate roots of this function are listed in the right-hand column of Table 1 in Murphy (1983). We have re-listed the values of six of these roots in a column immediately to the right of the above figure and have marked with orange circles the locations of these six points in the figure. As they should, Murphy's identified roots lie precisely at the intersection of the red line with the other curves.

Physically Viable Parameter Values

Interface Mapping

As has already been stated in our above description of the <math>~n_c = 1</math> core of these bipolytropic configurations, the structure of the core will be physically viable as long as the radial location of the interface, <math>~\xi_i</math>, between the core and the envelope is positioned somewhere within the range,

<math>~0 \le \xi_i \le \pi \, .</math>

Similarly, our discussion of Srivastava's Lane-Emden function, which is being used to define the envelope of these bipolytropic configurations, makes it clear that the envelope will have a physically viable structure as long as the parameter, <math>~\Delta_i</math>, associated with the radial location of the interface is positioned somewhere within the range,

<math>~\eta_\mathrm{crit} < e^{2\Delta_i} < e^{2\pi} \, ,</math>

where,

<math>\eta_\mathrm{crit} \equiv e^{2\tan^{-1}(1+2^{1/3})} = 10.05836783\, ,</math>

and, in keeping with the definition provided above,

<math>\Delta_i = \tan^{-1}(y_\mathrm{root}) + m\pi \, .</math>

This coordinate range for the physically viable envelope can be rewritten as,

<math>~\tan^{-1}(1+2^{1/3}) \le \Delta_i \le \pi \, .</math>

As is illustrated in the following figure, titled "Interface Mapping," our analytic solution defines a smooth, one-to-one mapping of the radial coordinate that defines the outer edge of the <math>~n_c=1</math> core, <math>~\xi_i</math>, to the parameter, <math>~\Delta_i</math>, that defines the inner edge of the <math>~n_e = 5</math> envelope.


Illustration of Interface Mapping
<math>~\mathrm{Model}</math> <math>~\xi_i</math> <math>~m</math> <math>~\Delta_i</math>




<math>~0</math> <math>~0</math> <math>~\tan^{-1}(1+2^{1/3})</math>
<math>~(2)</math> <math>~0.8154</math> <math>~0</math> <math>~1.24287</math>
<math>~(3)</math> <math>~1.6598</math> <math>~0</math> <math>~1.49179</math>
<math>~(4)</math> <math>~2.0914</math> <math>~1</math> <math>~1.73281</math>
<math>~(5)</math> <math>~2.7302</math> <math>~1</math> <math>~2.87493</math>
<math>~\pi</math> <math>~1</math> <math>~\pi</math>


In the figure, the green curve traces the segment of the E-Type Lane-Emden function for the <math>~n_c = 1</math> polytropic core, <math>~\theta_{1E}</math> — that is, the segment of the "sinc" function — that covers the range <math>~[0,\pi]</math> for the interface parameter, <math>~\xi_i</math>; the blue curve traces the segment of the F-Type Lane-Emden function for the <math>~n_e = 5</math> polytropic envelope, <math>~\phi_{5F}</math>, that covers this same range <math>~[0,\pi]</math>, but for the interface parameter, <math>~\Delta_i</math>. (This blue curve is also displayed, and its properties described in more depth, in a separate discussion devoted to Srivastava's Lane-Emden function for <math>~n=5</math> polytropes.) This entire displayed parameter range, <math>~[0,\pi]</math>, is associated with physically viable interface values for the core (green curve), but for the envelope (blue curve), only a subset of this range can be associated with physically viable interface values — namely, the range <math>~[\tan^{-1}(1+2^{1/3}),\pi]</math> over which <math>~\phi_{5F}</math> is positive but its slope is either negative or zero.


The purple dashed line segment labeled "(2)" in the above "Interface Mapping" figure intersects the green curve at the value, <math>~\xi_i = 0.8154</math>, and it intersects the blue curve at the interface value, <math>~\Delta_i = 1.24287</math>. [These two numerical values are listed in the row labeled "Model (2)" of the table that has been placed immediately to the right of the figure.] This is intended to illustrate that an equilibrium bipolytropic configuration with <math>~(n_c, n_e) = (1, 5)</math> can be constructed by truncating the <math>~n_c = 1</math> core at a radius, <math>~\xi_i = 0.8154</math>, and matching it to an <math>~n_e = 5</math> envelope whose innermost radius is associated with the parameter, <math>~\Delta_i = 1.24287</math>. (The corresponding, analytically determined values of the coefficients <math>~A_0</math> and <math>~B_0</math> dictate the manner in which the blue curve should be stretched both horizontally and vertically to complete a smooth attachment of the envelope to the core.) In an analogous fashion, the dashed line segments labeled "(3)", "(4)", and "(5)" — and their corresponding coordinate values as listed in the accompanying table — illustrate how three additional equilibrium bipolytropic configurations with progressively larger cores can be constructed with the proper interface mapping. We have specifically chosen to illustrate the interface-mapping of these four bipolytropic models because their cores are truncated at the same values of <math>~\xi_i</math> that Murphy (1983) used for the models numbered "(2)", "(3)", "(4)", and "(5)" in his Table 3. (See the much more in-depth discussion, below; and note that the edge of the core is labeled by the parameter, <math>~\zeta_J</math>, rather than by <math>~\xi_i</math>, in Murphy's Table 3.)


For clarity we note that, the analytically determined values of <math>~\Delta_i</math> that are listed in the table that sits to the right of the above "Interface Mapping" figure have been shifted in phase by <math>~m\pi</math>, where the relevant value of the integer, <math>~m</math>, is also listed in the table. In the figure, the two purple dashed-line segments are associated with models for which <math>~m=0</math>, while the two orange dashed-line segments are associated with models for which <math>~m=1</math>.


In his Table 3, Murphy (1983) also tabulates the characteristics of two additional models: The core of his "model 1" is truncated at a radius much less than <math>~\pi</math>, namely, <math>~\xi_i = 0.032678</math>; as a result, the bipolytropic configuration has a very small core and its structure is almost entirely that of an <math>~n = 5</math> polytrope. At the other extreme, the core of his "model 6" is truncated at a radius that is almost, but not quite, equal to <math>~\pi</math>, namely, <math>~\xi_i = 3.1415</math>; the resulting bipolytropic configuration has a very tiny envelope and its structure is almost entirely an <math>~n = 1</math> polytrope. Our analytic solution permits us to set <math>~\xi_i</math> to either of the two limiting values, <math>~0</math> or <math>~\pi</math>, and to show that, in these limits, <math>~\Delta_i</math> exactly equals, respectively, <math>~\tan^{-1}(1 + 2^{1/3})</math> and <math>~\pi</math>. The mapping of the first of these two limits is illustrated by the red dashed line segment in the above figure; the second limit is illustrated simply by the intersection of the two curves at the coordinate location, <math>~\pi</math>. It appears, therefore, that in these two limits Murphy's bipolytrope can be used, respectively, to define the structure of an isolated <math>~n=5</math> polytrope or an isolated <math>~n=1</math> polytrope.


ASIDE: It has previously been thought that Srivastava's Lane-Emden function, <math>~\phi_{5F}</math>, cannot be used on its own to define the structure of an isolated <math>~n=5</math> polytrope because the function's amplitude grows without bound and oscillates more and more rapidly between positive and negative values as the governing radial coordinate gets smaller and smaller. From the results presented here, it now appears as though an isolated <math>~n=5</math> polytrope of this type can be constructed by letting the interface parameter <math>~\xi_i \rightarrow 0</math> — and, hence the alternative parameter <math>~\Delta_i \rightarrow \tan^{-1}(1+2^{1/3})</math> — in Murphy's bipolytrope. It should be interesting to determine the values of the coefficients, <math>~A_0</math> and <math>~B_0</math>, that arise in this limit, and to examine in detail the structure of the complete <math>~n=5</math> model that results. Does its radial density profile resemble — or, perhaps, exactly match — the radial density profile of the well-known isolated <math>~n=5</math> ploytrope? Does the new structure have a finite radius as well as a finite central density? In this context it is worth noting that, in two separate papers — Murphy (1980a) and Murphy (1981) — Murphy has constructed and discussed the physical characteristics of equilibrium models that obey the <math>~n=5</math> polytropic equation of state all the way from the center to a surface which is of finite radius. He accomplishes this by piecing together a core that is defined by the familiar, analytically specified, <math>~\phi_{5E}</math> Lane-Emden function and an envelope that is defined by Srivastava's <math>~\phi_{5F}</math> Lane-Emden function.

Murphy's Example Model Characteristics

As in an earlier figure, the following log-log plot shows what value(s) of <math>~(A_0\eta)_\mathrm{root} </math> (vertical axis) result from a wide range of interface coordinate specifications, <math>~\xi_i</math> (horizontal axis), under the simplifying assumption that <math>~(\mu_e/\mu_c) = 1</math>. The solid blue diamonds present exactly the same information as in the previous figure, that is, they identify roots resulting from setting <math>~m=0</math> in the phase-shift specification. But here, focusing on a regime that highlights larger roots of the first constraint equation, solid black squares identify roots for which <math>~m=+1</math>, and solid maroon triangles identify roots for which <math>~m=+2</math>.


In this version of the figure, colored dashed-line segments have been drawn to enclose two regions of parameter space in which solutions of the key nonlinear interface relation produce a physically viable model — that is, <math>~0 \le \xi_i \le \pi</math> and <math>~\tan^{-1}(1+2^{1/3}) \le \Delta_i \le \pi</math>. The portion of the red-outlined rectangular region that lies to the left of the vertical, <math>~\xi_i = \xi_\mathrm{trans}</math> line accommodates some models that arise from specifying an <math>~m=0</math> phase shift and the portion that lies to the right of this demarcation line accommodates some models that arise from specifying an <math>~m=+1</math> phase shift. Similarly, the portion of the green-outlined rectangular region that lies to the left of the vertical, <math>~\xi_i = \xi_\mathrm{trans}</math> line accommodates additional models that arise from specifying an <math>~m=+1</math> phase shift and the portion that lies to the right of this demarcation line accommodates models that arise from specifying an <math>~m=+2</math> phase shift. Given the inherent periodicity of Srivastava's <math>~\phi_{5F}</math> function, it should be clear that we could equally well have identified other rectangular regions at both higher <math>~(m \ge 2)</math> and lower <math>~(m \le 0)</math> values of <math>~(A_0 \eta)_\mathrm{root}</math> in which additional physically viable models exist. It should also be clear that, because the figure presents a log-log plot and the lower limit of viable parameter values for the horizontal axis parameter, <math>~\xi_i</math>, is zero, the identified rectangular regions should be viewed as extending to the left, all the way to minus infinity.


Examination of Table 3 Data Presented in J. O. Murphy (1983, Proc. Astr. Soc. of Australia, 5, 175)

Table 3 Model Characteristics by Murphy (1983)
Comment by J. E. Tohline on 29 April 2015: As is argued more fully below, it appears as though the two parameter values drawn directly from Murphy's Table 3 that are shown here in red contain typographical errors.

Murphy's Table 3 Data

<math>~\zeta_J \rightarrow \xi_i</math>


<math>~\xi_J \rightarrow (A_0\eta)_\mathrm{root}</math>


3.2678E-02

10.0164

0.8154

12.0083

1.6598

20.4312

2.6914

33.0249

2.7302

100<math>\pi</math>

3.1415

91275<math>\pi</math>


In his Table 3, Murphy (1983) presents six pairs of "coordinate" roots, <math>~(\zeta_J, \xi_J)</math>, that not only satisfy the first constraint equation but that also, in his interpretation, identify physically viable model solutions. The values of all six of these coordinate pairs have been re-listed immediately to the right of the log-log plot, in a separate panel of the above figure. After mapping Murphy's coordinate notation to ours — specifically invoking the mapping <math>~(\zeta_J, \xi_J) ~\rightarrow~ (\xi_i, A_0\eta_\mathrm{root})</math> — the locations of all six models have been marked by orange circles on the plot. Five of Murphy's models fall directly onto our generated curves and in regions marked as physically viable, giving us confidence that our analysis matches his in considerable detail. As is explained more fully, below, we strongly suspect that the single Murphy model that does not lie on one of our generated curves — "Model 4" which, as published in Murphy's Table 3, has coordinates <math>~(\zeta_J, \xi_J) = (2.6914, 33.0249)</math> — is displaced, perhaps because of transcription errors in the published table.


As a supplement to the preceding two figures, we now display a table of our analytically determined model characteristics that is patterned after and largely matches Murphy's (1983) Table 3. The tabulated values have been derived assuming <math>~\mu_e/\mu_c = 1</math> throughout, that is, assuming that the mean molecular weights of the core and the envelope are identical. Vertically, the table is divided into three segments: The top portion contains six rows of numbers that result from our derivations; the middle portion displays six columns of numbers that have been extracted directly from Table 3 of Murphy (1983) and displayed here as image segments; and the bottom portion shows additional parameter values that we have generated in an effort to understand why our "Model 4" does not match the fourth model in Murphy's table.


Some Model Characteristics (assuming <math>\mu_e/\mu_c</math> = 1)

Model

Specified

Analytically Determined Here

<math>~\xi_i</math>

<math>~m</math>

<math>~\Delta_i</math>

<math>~(A_0\eta)_\mathrm{root}</math>

<math>~\frac{A_0}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) </math>

<math>~B_0</math>

<math>~\frac{B_0}{3^{1/4}} \biggl(\frac{\mu_e}{\mu_c}\biggr)^{1/2} \theta_i^{-1} </math>

<math>~\xi_s</math>

1 3.2678E-02 0 1.1544 10.0614 307.894 5.77929 4.39209 1.7392
2 0.8154 0 1.2429 12.0101 14.729 1.25567 1.06865 36.356
3 1.6598 0 1.4918 19.7585 11.9041 1.01206 1.28146 44.984
4 2.6914 1 2.8365 290.9431 108.1010 0.14348 0.67430 4.9536
5 2.7302 1 2.8749 100.00<math>\pi</math> 115.065 0.12408 0.64370 4.6539
6 3.1415 2 6.283141 91268<math>\pi</math> 170.44223<math>e^{2\pi}</math> 1.8966E-05 0.48862 3.141778
Comment by J. E. Tohline on 11 May 2015: The heading on this last column of numbers extracted from Murphy's (1983) paper has a subscript "J" whereas the subscript should be "R" to indicate the configuration's total radius.

Corresponding Values Extracted Directly from Murphy's (1983) Table 3

Columns extracted directly from Murphy (1983)

Attempting to Understand Murphy's "Model 4"

4b 0.048681 0 1.1545342 10.065043 206.75505 4.7358862 3.5999 2.589981
4c 1.14988 0 1.3234244 14.109506 12.2704 1.12444163 1.07640071 43.6409
4d 1.77027 0 1.5391618 21.721958 12.2704 0.9876639 1.355397 43.6409
4e 2.10961 1 1.7486314 33.0249 15.6545186 0.8651536 1.6157228 34.2068427

See an accompanying discussion for a complementary critique of Murphy's Table 3 models including an evlauation of the functions, <math>~U_{1E}</math>, <math>~V_{1E}</math>, <math>~U_{5F}</math>, and <math>~V_{5F} \, .</math>


Our six models (top portion of the above table) have been constructed by, first, specifying the location of the outer edge of the core, <math>~\xi_i</math> along with a corresponding phase shift parameter, <math>~m</math>, to match Murphy's chosen core-interface values — in his terminology, <math>~\zeta_J</math> (middle portion of the above table). Once the value of <math>~\xi_i</math> has been specified for each model, the value of the parameter, <math>~p</math>, the root of the cubic equation, <math>~y_\mathrm{root}</math>, and all of the parameter values listed in the upper portion of the table were determined using the parameter relations derived above and/or detailed in the table below. For example,

<math>~\Delta_i</math>

<math>~=</math>

<math>~\tan^{-1}(y_\mathrm{root}) + m\pi \, ;</math>

<math>~(A_0\eta)_\mathrm{root}</math>

<math>~=</math>

<math>~e^{2\Delta_i} \, ;</math>

and,

<math>~\xi_s</math>

<math>~=</math>

<math>~\xi_i e^{2(\pi -\Delta_i)} \, .</math>

Comparing our tabulated numbers (upper portion of the table) with the numbers published by Murphy (middle portion of the table) adds strong support to the statement made in connection with the above figure, that is, that our analysis matches Murphy's in considerable detail. But such a comparison also highlights the mismatch between Murphy's "Model 4" and ours — to aid in this particular comparison, the relevant numbers have been displayed in a red font in the upper portion of our table. According to our analysis, a model with <math>~\xi_i = 2.6914</math> does not display characteristics that, in any fashion, match Murphy's model with the same core-interface value.

We initially suspected that the "Model 4" mismatch resulted from a typographical or transcription error in the value of <math>~\zeta_J</math> — in our terminology, <math>~\xi_i</math> — that appears in Murphy's published Table 3. So we tried a variety of values of <math>~\xi_i</math>, looking for a set of generated characteristics that matched Murphy's published Model 4 characteristics. We were unsuccessful in this effort. The most tantalizing results from this search are displayed as cells with light-green backgrounds in the bottom segment of the above table: (Model 4b) Setting <math>~\xi_i = 0.048681~(m = 0)</math> gives <math>~B=3.5999</math>, which matches Murphy's "Model 4" value of <math>~B</math>; (Model 4c) setting <math>~\xi_i = 1.14988~(m = 0)</math> gives values of <math>~A</math> and <math>~\zeta_J (\mathrm{i.e.,}~\xi_s)</math> that both match Murphy's values; (Model 4d) setting <math>~\xi_i = 1.77027~(m = 0)</math> also gives values of <math>~A</math> and <math>~\zeta_J</math> that both match Murphy's values; and (Model 4e) setting <math>~\xi_i = 2.10961~(m = 1)</math> gives a value of <math>~(A_0\eta)_\mathrm{root}</math> that matches Murphy's value of <math>~\xi_J</math>. As we have been unable to find a value of <math>~\xi_i</math> whose associated parameter values match Murphy's "Model 4" values across the board, we remain at a loss to explain the mismatch between our "Model 4" and Murphy's.

Step 8: Envelope's Physical Profile

Combining the set of expressions provided in the right-most column of the accompanying bipolytrope "Table 1: Setup" with the interface conditions already established, above, we can now specify how the density, pressure, and integral mass vary throughout the envelope, that is, over the radial extent <math>\xi_i \le \xi \le \xi_s</math>.

<math>~\rho</math>

<math>~=</math>

<math>~\rho_e \phi^{5}</math>

 

<math>~=</math>

<math>~\rho_0 \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \phi^{5} \, ;</math>

<math>~P</math>

<math>~=</math>

<math>~K_e \rho_e^{6/5} \phi^{6}</math>

 

<math>~=</math>

<math>~K_c \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 \biggr]^{6/5} \phi^{6}</math>

 

<math>~=</math>

<math>~K_c \rho_0^2 \theta^{2}_i \phi^{6} \, ;</math>

<math>~M_r</math>

<math>~=</math>

<math>~4\pi \biggl[ \frac{3K_e}{2\pi G} \biggr]^{3/2} \rho_e^{-1/5} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>

 

<math>~=</math>

<math>~\biggl[ \frac{3^3 \cdot (2^2\pi)^2}{2^3\pi^3 G^3} \biggr]^{1/2} \biggl[K_c^5 \rho_0^{4}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-6} \theta^{4}_i\biggr]^{3/10} \biggl[ \rho_0 \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i \biggr]^{-1/5} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>

 

<math>~=</math>

<math>~\biggl[ \frac{2\cdot 3^3 }{\pi G^3} \biggr]^{1/2} K_c^{3/2} \rho_0 \theta_i \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr) \, .</math>

Examples

Normalization

The dimensionless variables used in Tables 1 & 2 are defined as follows:

<math>~\rho^*</math>

<math>~\equiv</math>

<math>~\frac{\rho}{\rho_0}</math>

;    

<math>~r^*</math>

<math>~\equiv</math>

<math>~\frac{r}{(K_c/G)^{1/2}}</math>

<math>~P^*</math>

<math>~\equiv</math>

<math>~\frac{P}{K_c\rho_0^{2}}</math>

;    

<math>~M_r^*</math>

<math>~\equiv</math>

<math>~\frac{M_r}{\rho_0 (K_c/G)^{3/2}}</math>

<math>~H^*</math>

<math>~\equiv</math>

<math>~\frac{H}{K_c\rho_0}</math>

.    

 

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.


Properties of <math>~n_c=1</math>, <math>~n_e=5</math>, BiPolytrope Having Various Interface Locations, <math>~\xi_i</math>

Parameter

<math>~\xi_i</math>

<math>~\theta_i</math>

<math>~\frac{\sin\xi_i}{\xi_i} </math>

<math>~-\biggl(\frac{d\theta_i}{d\xi}\biggr)_i</math>

<math>~\frac{1}{\xi_i^2}(\sin\xi - \xi_i \cos\xi_i) </math>

<math>~r^*_\mathrm{core} \equiv r^*_i</math>

<math>~\biggl( \frac{1}{2\pi} \biggr)^{1/2} \xi_i</math>

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

<math>~\theta_i</math>

<math>P^*_i</math>

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

<math>H^*_i \biggr|_c = \frac{n_c+1}{n_e+1} \biggl( \frac{\mu_e}{\mu_c} \biggr) H^*_i \biggr|_e</math>

<math>~2\theta_i</math>

<math>M^*_\mathrm{core}</math>

<math>~\biggl( \frac{2}{\pi} \biggr)^{1/2} (\sin\xi_i - \xi_i \cos\xi_i)</math>

<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1}\eta_i</math>

<math>\frac{\xi_i}{\sqrt{3}}</math>

<math>~ -\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>

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

<math>~p</math>

<math>~\frac{3}{ 3-2(\mu_e/\mu_c)(1-\xi_i \cot\xi_i)}</math>

<math>~~y_\mathrm{root}</math>

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

<math>~\Delta_i - m\pi</math>

<math>~\tan^{-1}(y_\mathrm{root})</math>

<math>~(A_0\eta)_\mathrm{root}</math>

<math>~e^{2\Delta_i}</math>

<math>~\xi_s</math>

<math>~\xi_ie^{2(\pi-\Delta_i)}</math>

<math>~\biggl[ \frac{A_0}{3^{1/2}} \biggl( \frac{\mu_e}{\mu_c} \biggr) \biggr]</math>

<math>~\xi_i^{-1} e^{2\Delta_i}</math>

<math>~\biggl[\frac{B_0}{3^{1/4}} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1/2} \biggr]</math>

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

<math>~ -\biggl( \frac{d\phi}{d\eta} \biggr)_s \biggl( \frac{\mu_e}{\mu_c}\biggr)</math>

<math>~\frac{1}{ 2\xi_i e^{3(\pi-\Delta_i)}} \biggl( \frac{3}{\sin^2\Delta_i} - 2 \biggr)^{1/2}</math>

<math>~ R^*_s</math>

<math>~(2\pi)^{-1/2} \xi_i e^{2(\pi - \Delta_i)}</math>

<math>~ \biggl( \frac{\mu_e}{\mu_c}\biggr) M^*_\mathrm{tot}</math>

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


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 step #4 and step #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 2nd 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.

Table 2: Radial Profile of Various Physical Variables

Variable

Throughout the Core
<math>0 \le \xi \le \xi_i</math>

Throughout the Envelope
<math>\eta_i \le \eta \le \eta_s</math>

Plotted Profiles

<math>\xi_i = 0.5</math>

<math>\xi_i = 1.0</math>

<math>\xi_i = 3.0</math>

<math>~r^*</math>

<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \xi</math>

<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \biggl(\frac{3}{2\pi}\biggr)^{1/2}\eta</math>

 

<math>~\rho^*</math>

<math>\frac{\sin\xi}{\xi}</math>

<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta_i [\phi(\eta)]^5</math>

DenXi05.jpg
DenXi10.jpg
DenXi30.jpg

<math>~P^*</math>

<math>\biggl( \frac{\sin\xi}{\xi} \biggr)^2</math>

<math>\theta^{2}_i [\phi(\eta)]^{6}</math>

PresXi05.jpg
PresXi10.jpg
PresXi30.jpg

<math>~M_r^*</math>

<math>\biggl( \frac{2}{\pi}\biggr)^{1/2} (\sin\xi - \xi\cos\xi)</math>

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

MassXi05.jpg
MassXi10.jpg
MassXi30.jpg

In order to obtain the various envelope profiles, it is necessary to evaluate <math>~\phi(\eta)</math> and its first derivative using the information presented in Step 6, above.

Murphy and Fiedler (1985)

Table 1 from Murphy & Fiedler (1985, Proc. Astr. Soc. of Australia, 6, 219)

Murphy & Fiedler (1985) Table 1

Reproduction of Table 1 from MF85 Using Excel and Analytic Expressions Derived Here

Excel Regeneration of MF85 Table 1

Key References


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