Editing SSC/Structure/BiPolytropes/MurphyUVplane (section)

==Critique of Murphy's Model Characteristics==


{{ SGFworkInProgress }}


As a supplement to the preceding two figures &#8212; [[SSC/Structure/BiPolytropes/Analytic15/Pt2#MurphyF2Function|one, here]] and [[SSC/Structure/BiPolytropes/Analytic15/Pt2#SecondMurphyF2Function|the other, here]] &#8212; the following table shows in quantitative detail how we are able to ''analytically'' derive individual coordinate roots, <math>(A_0\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_0\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 {{ Murphy83a }}.  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.  


<div align="center">
<table border="1" cellpadding="7">
<tr>
  <td align="center" rowspan="2">
{{ Murphy83a }}<br />
Tables 1, 2, &amp; 3
  </td>
  <td align="center" colspan="6">
Here &hellip; setting <math>\xi_i = \zeta_J</math> &hellip; and assuming <math>\mu_e/\mu_c = 1</math>
  </td>
  <td align="center" rowspan="2">
{{ Murphy83a }}<br />
Tables 1, 2, &amp; 3
  </td>
</tr>
<tr>
  <td align="center" colspan="6">
1<sup>st</sup> Constraint
  </td>
</tr>

<tr>
  <td align="center">
<math>\zeta_J</math>
  </td>
  <td align="center">
<math>\kappa_i</math>
  </td>
  <td align="center">
<math>p</math>
  </td>
  <td align="center">
<math>y_\mathrm{root}</math>
  </td>
  <td align="center">
<math>(\Delta_i-m\pi)</math><p></p><math>~= \tan^{-1}(y_\mathrm{root})</math>
  </td>
  <td align="center">
<math>m </math>
  </td>
  <td align="center">
<math>(A_0\eta)_\mathrm{root} = e^{2\Delta_i}</math>
  </td>
  <td align="center">
<math>\xi_J</math>
  </td>
</tr>

<tr>
  <td align="center" rowspan="4" bgcolor="#FFBB88">
3.2678E-02
  </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 align="center">
-3
  </td>
  <td align="center">
6.55237E-08
  </td>
  <td align="center" bgcolor="lightblue">
6.55E-08
  </td>
</tr>

<tr>
  <td align="center">
-2
  </td>
  <td align="center">
3.50874E-05
  </td>
  <td align="center" bgcolor="lightblue">
3.508E-05
  </td>
</tr>

<tr>
  <td align="center">
-1
  </td>
  <td align="center">
1.87890E-02
  </td>
  <td align="center" bgcolor="lightblue">
1.8789E-02
  </td>
</tr>

<tr>
  <td align="center">
0
  </td>
  <td align="center">
10.0614
  </td>
  <td align="center" bgcolor="#FFBB88">
<font color="red">10.0164</font>
  </td>
</tr>

<tr>
  <td align="center" bgcolor="#FFBB88">
0.8154
  </td>
  <td align="center">
0.15474
  </td>
  <td align="center">
1.18307
  </td>
  <td align="center">
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>
</tr>

<tr>
  <td align="center" bgcolor="#FFBB88">
1.6598
  </td>
  <td align="center">
0.76541
  </td>
  <td align="center">
4.26281
  </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>
</tr>

<tr>
  <td align="center">
[&nbsp;2.0914&nbsp;]
  </td>
  <td align="center">
1.46609
  </td>
  <td align="center">
-2.14552
  </td>
  <td align="center">
-6.11816
  </td>
  <td align="center">
-1.40878
  </td>
  <td align="center">
1
  </td>
  <td align="center">
31.9964
  </td>
  <td align="center" rowspan="3" bgcolor="#FFBB88">
33.0249
  </td>
</tr>

<tr>
  <td align="center">
[&nbsp;2.10961&nbsp;]
  </td>
  <td align="center">
1.50744
  </td>
  <td align="center">
-1.97067
  </td>
  <td align="center">
-5.56379
  </td>
  <td align="center">
-1.39296
  </td>
  <td align="center">
1
  </td>
  <td align="center">
33.0249
  </td>
</tr>

<tr>
  <td align="center" bgcolor="#FFBB88">
<font color="red">2.6914</font>
  </td>
  <td align="center">
4.37925
  </td>
  <td align="center">
-0.29592
  </td>
  <td align="center">
-0.31486
  </td>
  <td align="center">
-0.30503
  </td>
  <td align="center">
1
  </td>
  <td align="center">
290.943
  </td>
</tr>

<tr>
  <td align="center" bgcolor="#FFBB88">
2.7302
  </td>
  <td align="center">
4.83853
  </td>
  <td align="center">
-0.26052
  </td>
  <td align="center">
-0.27316
  </td>
  <td align="center">
-0.26666
  </td>
  <td align="center">
1
  </td>
  <td align="center">
99.997<math>\pi</math>
  </td>
  <td align="center" bgcolor="#FFBB88">
100<math>\pi</math>
  </td>
</tr>

<tr>
  <td align="center" bgcolor="#FFBB88">
3.1415
  </td>
  <td align="center">
22604.6
  </td>
  <td align="center">
-4.42408E-05
  </td>
  <td align="center">
-4.42408E-05
  </td>
  <td align="center">
-4.42408E-05
  </td>
  <td align="center">
2
  </td>
  <td align="center">
91268<math>\pi</math>
  </td>
  <td align="center" bgcolor="#FFBB88">
91275<math>\pi</math>
  </td>
</tr>

<tr>
  <td align="center" bgcolor="lightgreen">
4.2993
  </td>
  <td align="center">
-0.58961
  </td>
  <td align="center">
0.62908
  </td>
  <td align="center">
0.89355
  </td>
  <td align="center">
0.72924
  </td>
  <td align="center">
0
  </td>
  <td align="center">
4.2994
  </td>
  <td align="center" bgcolor="lightgreen">
4.2993
  </td>
</tr>

<tr>
  <td align="center" bgcolor="lightgreen">
7.6882
  </td>
  <td align="center">
-0.19091
  </td>
  <td align="center">
0.83969
  </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>
</tr>

</table>
</div>



As this table and the log-log plots presented in two preceding figures illustrate &#8212; [[SSC/Structure/BiPolytropes/Analytic15/Pt2#MurphyF2Function|one, here]] and [[SSC/Structure/BiPolytropes/Analytic15/Pt2#SecondMurphyF2Function|the other, here]] &#8212; by and large, the agreement between our ''analytically'' determined coordinate roots and those derived &amp; published by {{ Murphy83a }} is excellent.  There are a couple of mismatches, but we attribute both to typesetting errors in the Murphy publication:  
* 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.
* 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_0\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).

<!--  OMIT FOLLOWING SECTION

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">
<math>~A = \frac{(A\eta)_\mathrm{root}}{\xi_i} \, .</math>
</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 [[SSC/Structure/BiPolytropes/Analytic15#Step_6:__Envelope_Solution|Srivastava's function]], <math>\phi</math>, drops to zero.  That is, the surface occurs where,

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

<tr>
  <td align="right">
<math>~\Delta_s</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\pi</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~~\xi_s = \eta_s</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{A} \cdot e^{2\pi} \, .</math>
  </td>
</tr>
</table>
</div>

STOP OMISSION -->


<div align="center">
<table border="1" cellpadding="5" align="center" width="90%">
<tr>
  <td align="center" colspan="4">
Model Characteristics from Table 3 of &hellip;<br />{{ Murphy83afigure }}
  </td>
  <td align="center" colspan="4">
Implications
  </td>
</tr>
<tr>
  <td align="center">
Model
  </td>
  <td align="center">
<math>\xi = \zeta_J</math>
  </td>
  <td align="center">
<math>(A\xi) = \xi_J</math>
  </td>
  <td align="center">
<math>A = \xi_J/\zeta_J</math>
  </td>
  <td align="center">
<math>\Delta = \ln(A\xi)^{1/2}</math>
  </td>
  <td align="center">
<math>V_{1E}</math>
  </td>
  <td align="center">
<math>V_{5F}</math>
  </td>
  <td align="center">
<math>\frac{(n_c+1)V_{1E}}{(n_e+1)V_{5F}}</math>
  </td>
</tr>
<tr>
  <td align="center">
1
  </td>
  <td align="right">
0.032678
  </td>
  <td align="right">
<font color="red">10.0614</font>
  </td>
  <td align="right">
307.895
  </td>
  <td align="right">
1.1544
  </td>
  <td align="right">
3.5598E-04
  </td>
  <td align="right">
1.1871E-04
  </td>
  <td align="center">
1.000
  </td>
</tr>

<tr>
  <td align="center">
2
  </td>
  <td align="right">
0.8154
  </td>
  <td align="right">
12.0083
  </td>
  <td align="right">
14.727
  </td>
  <td align="right">
1.2428
  </td>
  <td align="right">
2.3212E-01
  </td>
  <td align="right">
7.7302E-02
  </td>
  <td align="center">
1.001
  </td>
</tr>
<tr>
  <td align="center">
3
  </td>
  <td align="right">
1.6598
  </td>
  <td align="right">
20.4312
  </td>
  <td align="right">
12.309
  </td>
  <td align="right">
1.5085
  </td>
  <td align="right">
1.1481
  </td>
  <td align="right">
0.40720
  </td>
  <td align="center">
0.940
  </td>
</tr>

<tr>
  <td align="center">
4
  </td>
  <td align="right">
<font color="red">2.6914</font>
  </td>
  <td align="right">
33.0249
  </td>
  <td align="right">
12.2704
  </td>
  <td align="right">
1.7486
  </td>
  <td align="right">
6.5689
  </td>
  <td align="right">
0.75371
  </td>
  <td align="center">
2.905
  </td>
</tr>

<tr>
  <td align="center">
5
  </td>
  <td align="right">
2.7302
  </td>
  <td align="right">
314.159
  </td>
  <td align="right">
115.068
  </td>
  <td align="right">
2.8750
  </td>
  <td align="right">
7.2578
  </td>
  <td align="right">
2.4193
  </td>
  <td align="center">
1.000
  </td>
</tr>

<tr>
  <td align="center">
6
  </td>
  <td align="right">
3.1415
  </td>
  <td align="right">
2.8675E+05
  </td>
  <td align="right">
9.1278E+04
  </td>
  <td align="right">
6.2832
  </td>
  <td align="right">
3.3907E+04
  </td>
  <td align="right">
2.9313E+04
  </td>
  <td align="center">
0.386
  </td>
</tr>

<tr>
  <td align="left" colspan="8">
<font color="red">NOTE:</font> As is [[SSC/Structure/BiPolytropes/Analytic15#Murphy.27s_Example_Model_Characteristics|explained in an accompanying discussion]], we suspect that the two numbers drawn from Murphy's Table 3 that are displayed here in a red font contain typographical errors.
  </td>
</tr>

</table>
</div>



<div align="center">
<table border="1" cellpadding="5">
<tr>
  <td align="center" colspan="14">
Some Model Characteristics (assuming <math>\mu_e/\mu_c</math> = 1)
  </td>
</tr>
<tr>
  <td align="center" rowspan="2">
Model
  </td>
  <td align="center" colspan="2">
Specified
  </td>
  <td align="center" colspan="6">
Analytically Determined Here
  </td>
  <td align="center" colspan="5">
&hellip; and Constraint Implications<sup>&dagger;</sup>
  </td>
</tr>

<tr>
  <td align="center">
<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 align="center">
<math>~\frac{B_0}{3^{1/4}} \biggl(\frac{\mu_e}{\mu_c}\biggr)^{1/2} \theta_i^{-1} </math>
  </td>
  <td align="center">
<math>~\xi_s</math>
  </td>
  <td align="center">
<math>~V_{1E}</math>
  </td>
  <td align="center">
<math>~V_{5F}</math>
  </td>
  <td align="center">
<math>\frac{(n_c+1)V_{1E}}{(n_e+1)V_{5F}}</math>
  </td>
  <td align="center">
<math>~U_{1E}</math>
  </td>
  <td align="center">
<math>\frac{U_{1E}}{U_{5F}}</math>
  </td>
</tr>

<tr>
  <td align="center">1</td>
  <td align="right">3.2678E-02</td>
  <td align="center">0</td>
  <td align="right">1.1544</td>
  <td align="right">10.0614</td>
  <td align="right">307.894</td>
  <td align="center">5.77929</td>
  <td align="center">4.39209</td>
  <td align="right">1.7392</td>
  <td align="right">3.5598E-04</td>
  <td align="right">1.1866E-4</td>
  <td align="center">1.000</td>
  <td align="right">2.9998</td>
  <td align="center">1.000</td>
</tr>

<tr>
  <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>
  <td align="right">0.23212</td>
  <td align="right">0.077372</td>
  <td align="center">1.000</td>
  <td align="right">2.8644</td>
  <td align="center">1.000</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>
  <td align="right">1.14812</td>
  <td align="right">0.38271</td>
  <td align="center">1.000</td>
  <td align="right">2.3995</td>
  <td align="center">1.000</td>
</tr>

<tr>
  <td align="center">4</td>
  <td align="right">2.0914</td>
  <td align="center">1</td>
  <td align="right">1.7328</td>
  <td align="right">31.9964</td>
  <td align="right">15.2990</td>
  <td align="center">0.87563</td>
  <td align="center">1.60398</td>
  <td align="right">35.0016</td>
  <td align="right">2.19913</td>
  <td align="right">0.73304</td>
  <td align="center">1.000</td>
  <td align="right">1.98894</td>
  <td align="center">1.000</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>
  <td align="right">7.25779</td>
  <td align="right">2.41926</td>
  <td align="center">1.000</td>
  <td align="right">1.02703</td>
  <td align="center">1.000</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>
  <td align="right">3.3907E+04</td>
  <td align="right">1.1302E+04</td>
  <td align="center">1.000</td>
  <td align="right">2.9106E-04</td>
  <td align="center">1.000</td>
</tr>
<tr>
  <td align="center" colspan="9">
Corresponding Values Extracted Directly from Table 3 of &hellip;<br />{{ Murphy83afigure }}.
  </td>
  <td align="center" colspan="5">
&nbsp;
  </td> 
</tr>

<tr>
  <td align="center"><!-- [[File:Table3CompleteCol1.png|center|50px|Column 1 from Table 3 of Murphy (1983)]] -->
1<br />2<br />3<br />4<br />5<br />6
  </td>
  <td align="left" colspan="2"><!-- [[File:Table3CompleteCol2.png|center|123px|Column 3 from Table 3 of Murphy (1983)]] -->
3.2678 &times; 10<sup>-2</sup><br />0.8154<br />1.6598<br />2.6914<br />2.7302<br />3.1415
  </td>
  <td align="center">&nbsp;</td>
  <td align="left"><!-- [[File:Table3CompleteCol3.png|center|67px|Column 3 from Table 3 of Murphy (1983)]] -->
10.0164<br />12.0083<br />20.4312<br />33.0249<br />100&pi;<br />91275&pi;
  </td>
  <td align="right"><!-- [[File:Table3CompleteCol4.png|center|117px|Column 4 from Table 3 of Murphy (1983)]] -->
307.8945<br />14.7269<br />12.3093<br />12.2704<br />115.0678<br />170.4515e<sup>2&pi;</sup>
  </td>
  <td align="center">&nbsp;</td>
  <td align="right"><!-- [[File:Table3CompleteCol5.png|center|71px|Column 5 from Table 3 of Murphy (1983)]] -->
4.3921<br />1.0686<br />1.2859<br />3.5999<br />0.6437<br />0.3357
  </td>
  <td align="right"><!-- [[File:Table3CompleteCol6.png|center|83px|Column 6 from Table 3 of Murphy (1983)]] -->
1.7391<br />36.3615<br />43.5030<br />43.6409<br />4.4654<br />3.1416
  </td>
  <td align="left" colspan="5">[[File:CommentButton02.png|left|100px|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.]]</td>
</tr>

<tr>
  <td align="left" colspan="14">
<sup>&dagger;</sup>See an [[SSC/Structure/BiPolytropes/MurphyUVplane#Chandrasekhar.27s_U_and_V_Functions|accompanying discussion]] for definitions of the functions, <math>U_{1E}</math>, <math>V_{1E}</math>, <math>U_{5F}</math>, and <math>V_{5F} \, .</math>
  </td>
</tr>

</table>
</div>