Editing SSC/Structure/PolytropesEmbedded/n1 (section)

=Embedded Polytropic Spheres=
<table border="1" align="center" width="100%" colspan="8">
<tr>
  <td align="center" bgcolor="lightblue" width="25%"><br />[[SSC/Structure/PolytropesEmbedded|Part I: &nbsp; General Properties]]
&nbsp;
  </td>
  <td align="center" bgcolor="lightblue" width="25%"><br />[[SSC/Structure/PolytropesEmbedded/n1|Part II:&nbsp; Truncated Configurations with n = 1]]
&nbsp;
  </td>
  <td align="center" bgcolor="lightblue" width="25%"><br />[[SSC/Structure/PolytropesEmbedded/n5|Part III:&nbsp; Truncated Configurations with n = 5]]
&nbsp;
  </td>
  <td align="center" bgcolor="lightblue"><br />[[SSC/Structure/PolytropesEmbedded/Other|Part IV:&nbsp; Other Considerations]]
&nbsp;
  </td>
</tr>
</table>



==Truncated Configurations with n = 1==

Drawing from the [[SSC/Structure/Polytropes|earlier discussion of isolated polytropes]], we will reference various radial locations within the spherical configuration by the dimensionless radius,
<div align="center">
<math>
\xi \equiv \frac{r}{a_\mathrm{n=1}} ,
</math>
</div>
where,
<div align="center">
<math>
a_\mathrm{n=1} \equiv \biggl[\frac{1}{4\pi G}~ \biggl( \frac{H_c}{\rho_c} \biggr)_{n=1}\biggr]^{1/2} = \biggl[\frac{K}{2\pi G} \biggr]^{1/2} \, .
</math>
</div>
The solution to the Lane-Emden equation for <math>~n = 1</math> is,
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>
~\theta_1
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\frac{\sin\xi}{\xi} \, ,
</math>
  </td>
</tr>
</table>
</div>
hence,
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>
\frac{d\theta_1}{d\xi}
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \, .
</math>
  </td>
</tr>
</table>
</div>

<font color="darkblue">
===Review===
</font>

Again, from the [[SSC/Structure/Polytropes|earlier discussion]], we can describe the properties of an isolated, spherical {{Math/MP_PolytropicIndex}} = 1 polytrope as follows:
* <font color="red">Mass</font>:  
: In terms of the central density, <math>\rho_c</math>, and {{Math/MP_PolytropicConstant}}, the total mass is,
<div align="center">
<math>M = \frac{4}{\pi} \rho_c (\pi a_{n=1})^3 = 4\pi^2 \rho_c \biggl[\frac{K}{2\pi G} \biggr]^{3/2} = \rho_c \biggl[\frac{2\pi K^3}{G^3} \biggr]^{1/2}</math> ;
</div>

: and, expressed as a function of <math>M</math>, the mass that lies interior to the dimensionless radius <math>\xi</math> is,
<div align="center">
<math>\frac{M_\xi}{M} = \frac{1}{\pi} \biggl[ \sin\xi - \xi\cos\xi \biggr] \, ,~~~~~~\mathrm{for}~\pi \ge \xi \ge 0 \, .</math> 
</div>
: Hence,
<div align="center">
<math>M_\xi = \rho_c \biggl[\frac{2K^3}{\pi G^3} \biggr]^{1/2} \biggl[ \sin\xi - \xi\cos\xi \biggr] \, .</math> 
</div>

* <font color="red">Pressure</font>: 
: The central pressure of the configuration is,
<div align="center">
<math>P_c = \biggl[ \frac{G^3}{2\pi} \rho_c^4 M^2 \biggr]^{1/3} = \biggl[ \frac{G^3}{2\pi} \rho_c^6 \biggl(\frac{2\pi K^3}{G^3} \biggr) \biggr]^{1/3} = K\rho_c^2</math> ; 
</div>

: and, expressed in terms of the central pressure <math>P_c</math>, the variation with radius of the pressure is,
<div align="center">
<math>P_\xi= P_c \biggl[ \frac{\sin\xi}{\xi} \biggr]^2</math> .
</div>
: Hence,
<div align="center">
<math>P_\xi= K\rho_c^2 \biggl[ \frac{\sin\xi}{\xi} \biggr]^2</math> .
</div>

===Extension to Bounded Sphere===
Eliminating <math>\rho_c</math> between the last expression for <math>M_\xi</math> and the last expression for <math>P_\xi</math> gives,
<div align="center">
<math>P_\xi= \biggl[\frac{\pi}{2} \cdot \frac{G^3 M_\xi^2}{K^2} \biggr] \biggl[ \frac{\sin\xi}{\xi(\sin\xi - \xi \cos\xi )} \biggr]^2</math> .
</div>
Now, if we rip off an outer layer of the star down to some dimensionless radius <math>\xi_e < \pi</math>, the interior of the configuration that remains &#8212; containing mass <math>M_{\xi_e}</math> &#8212; should remain in equilibrium if we impose the appropriate amount of externally applied pressure <math>P_e = P_{\xi_e} </math> at that radius.  (This will work only for spherically symmetric configurations, as the gravitation acceleration at any location only depends on the mass contained inside that radius.)  If we rescale our solution such that the mass enclosed within <math>\xi_e</math> is the original total mass <math>M</math>, then the pressure that must be imposed by the external medium in which the configuration is embedded is,
<div align="center">
<math>P_e= \biggl[\frac{\pi}{2} \cdot \frac{G^3 M^2}{K^2} \biggr] \biggl[ \frac{\sin\xi_e}{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )} \biggr]^2</math> .
</div>
The associated equilibrium radius of this pressure-confined configuration is,
<div align="center">
<math>
R_\mathrm{eq} = \xi_e a_\mathrm{n=1} = \biggl[ \frac{K}{2\pi G} \biggr]^{1/2} \xi_e 
</math>
</div>

====Overlap with Whitworth's Presentation====
The solid green curve in the two top panels of Figure 1 shows how <math>R_\mathrm{eq}</math> varies with the applied external pressure <math>P_e</math> for this pressure-bounded <math>~n=1</math> model sequence.  In the top-right panel, following the lead of [http://adsabs.harvard.edu/abs/1981MNRAS.195..967W Whitworth] (1981, MNRAS, 195, 967) &#8212; for clarification, read the [[SSC/Structure/PolytropesASIDE1|accompanying ASIDE]] &#8212; these two quantities have been respectively normalized (or, "referenced") to,
<div align="center">
<math>
R_\mathrm{rf}\biggr|_\mathrm{n=1} \equiv \biggl( \frac{3^2 \cdot 5}{2^4 \pi} \biggr)^{1/2} \biggl(\frac{K}{G}\biggr)^{1/2} ~~~\Rightarrow ~~~ \frac{R_\mathrm{eq}}{R_\mathrm{rf}} = \biggl( \frac{2^3}{3^2 \cdot 5} \biggr)^{1/2} \xi_e \, ;
</math>
</div>
and,
<div align="center">
<math>
P_\mathrm{rf}\biggr|_\mathrm{n=1} \equiv \frac{2^6 \pi}{3^4 \cdot 5^3}  \biggl(\frac{G^3 M^2}{K^2}\biggr) ~~~\Rightarrow ~~~ \frac{P_e}{P_\mathrm{rf}} = \biggl( \frac{3^4 \cdot 5^3}{2^7} \biggr) \biggl[ \frac{\sin\xi_e}{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )} \biggr]^2 \, .
</math>
</div>
Note that this pair of mathematical expressions has been recorded to the immediate right of Whitworth's name in our [[SSC/Structure/PolytropesEmbedded#n1Summary|<math>~n=1</math> summary table]].
In the top-left panel of Figure 1, the solid green curve shows the identical sequence, but plotted as <math>~\log(p_a)</math> versus <math>~log(r_a)</math>, for easier comparison with Horedt's work.  The pair of mathematical expressions defining <math>~r_a(\xi_e)</math> and <math>~p_a(\xi_e)</math> has been recorded to the immediate right of Horedt's name in the same [[SSC/Structure/PolytropesEmbedded#n1Summary|summary table]].


<span id="WhitworthFig1b">
<div align="center">
<table border="2" cellpadding="8" width="85%">

<tr>
  <td align="center" colspan="2">
'''Figure 1:''' <font color="darkblue">
Equilibrium R-P Diagram &#8212; Referred to by [http://adsabs.harvard.edu/abs/1981PASJ...33..273K Kimura (1981)] as an "M<sub>1</sub> Sequence"
</font> 
  </td>
</tr>
<tr>
  <td align="left" colspan="2">
All of the plots shown in this figure illustrate how the equilibrium radius of a pressure-bounded polytrope varies with the applied external pressure.  In the right-hand column, the log-log plots display a normalized <math>~P_e</math> along the horizontal axis and a normalized <math>~R_\mathrm{eq}</math> along the horizontal axis; in the left-hand column, these axes are flipped, and a different normalization is used.  One primary intent of all the diagrams is to show that, for polytropic sequences having <math>~n > 3</math> (or, equivalently, sequences having <math>\gamma_g \equiv 1 + 1/n < 4/3),</math> no equilibrium models exist above some limiting external pressure.
  </td>
</tr>

<tr>
  <td align="center" bgcolor="white">
[[File:HoredtPlot2.png|250px|center|To be compared with Horedt (1970)]]
  </td>
  <td align="center" bgcolor="white">
[[File:WhitworthPlot2.png|250px|center|To be compared with Whitworth (1981)]]
  </td>
</tr>

<tr>
  <td align="center" bgcolor="white">
[[File:Horedt_PRdiagram0.png|250px|center|Horedt (1970) Figure 1]]
<!-- [[Image:AAAwaiting01.png|250px|center|Horedt (1970) Figure 1]] -->
  </td>
  <td align="center" bgcolor="white">
[[File:WhitworthFig1bCopy.jpg|300px|center|Whitworth (1981) Figure 1b]]
<!--[[Image:AAAwaiting01.png|300px|center|Whitworth (1981) Figure 1b]] -->
  </td>
</tr>

<tr>
  <td align="center" bgcolor="white">
[[File:Horedt_EmbeddedPolytrope.png|300px|center|Horedt (1970) Title Page]] 
<!--[[Image:AAAwaiting01.png|300px|center|Horedt (1970) Title Page]] -->
  </td>
  <td align="center" bgcolor="white">
[[File:Whitworth1981TitlePage0.png|200px|center|Whitworth (1981) Title Page]]
<!--[[Image:AAAwaiting01.png|200px|center|Whitworth (1981) Title Page]] -->
  </td>
</tr>

<tr>
  <td align="left" colspan="2">
''Bottom Left'' [reproduction of Figure 1 from [http://adsabs.harvard.edu/abs/1970MNRAS.151...81H Horedt (1970)]]:  All three displayed sequences &#8212; <math>~n=4</math> (<math>~\gamma_g = 1.25</math>), <math>~n=5</math> (<math>~\gamma_g = 1.20</math>), and <math>~n=\infty</math> (<math>~\gamma_g = 1</math>, hence, isothermal) &#8212; exhibit an upper limit for the bounding pressure.  Each sequence displays two segments &#8212; a solid segment and a dashed segment &#8212; indicating that, below the maximum allowed value of <math>~P_e</math>, it is possible to construct two (or more) equilibrium configurations; models lying along the solid segment of each displayed curve are expected to be dynamically stable while models lying along the dashed segments are unstable.

''Bottom Right'' [reproduction of Figure 1b from [http://adsabs.harvard.edu/abs/1981MNRAS.195..967W Whitworth (1981)]]:  Model sequences are shown for five different effective adiabatic indexes &#8212; <math>~\gamma_g = 1/3,~ 2/3,~ 1,~ 4/3,</math> and <math>~ 5/3</math> &#8212; corresponding, respectively, to polytropic indexes <math>~n = -2/3, -1/3, \infty, ~3/2, </math> and <math>~3</math>.  The three sequences having <math>~\gamma_g < 4/3</math> exhibit an upper limit for the bounding pressure.  Both the stable (solid) curve segment and the unstable (dashed) curve segment are drawn for the isothermal <math>~(\gamma_g = 1)</math> sequence, which is also displayed (as the <math>~n=\infty</math> sequence) in Horedt's diagram.

''Top'':  Plots that we have generated for direct comparison with Horedt's diagram (''left'') and with Whitworth's diagram (''right'').  Both plots display only the two sequences that are analytically prescribable:   <math>~n=1</math> (<math>~\gamma_g = 2</math>) and <math>~n=5</math> (<math>~\gamma_g = 1.20</math>).  Along the <math>~n=1</math> (green) sequence, stable equilibrium models can be constructed for all values of <math>~P_e</math>.  Along the <math>~n=5</math> sequence, equilibrium models only exist for values of <math>~P_e</math> less than the critical value, <math>~P_\mathrm{max} = (2^5\cdot 3^9/5^9) P_\mathrm{rf} = (3^{12}/2^{24}) P_\mathrm{Horedt}</math>; below this critical pressure, the sequence has two branches denoted by blue diamonds (stable models) and red squares (unstable models).
  </td>
</tr>
</table>
</div>
</span>

====Overlap with Stahler's Presentation====
We can invert the above expression for <math>~P_e(K,M)</math> to obtain the following expression for <math>~M(K,P_e)</math>:
<div align="center">
<math>~M= K \biggl[\frac{2}{\pi} \cdot \frac{P_e}{G^3} \biggr]^{1/2} \biggl[ \frac{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )}{\sin\xi_e} \biggr]</math> .
</div>
If, following Stahler's lead, we normalize this expression by <math>~M_\mathrm{SWS}</math> (evaluated for <math>~n=1</math>) and we normalize the above expression for <math>~R_\mathrm{eq}</math> by <math>~R_\mathrm{SWS}</math> (evaluated for <math>~n=1</math>), we obtain,
<div align="center">
<table border="0" cellpadding="3">

<tr>
  <td align="right">
<math>
\frac{M}{M_\mathrm{SWS}}
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
K \biggl[\frac{2}{\pi} \cdot \frac{P_e}{G^3} \biggr]^{1/2} \biggl[ \frac{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )}{\sin\xi_e} \biggr]
\biggl[ \biggl( \frac{G}{2} \biggr)^{3/2} K^{-1} P_\mathrm{ex}^{-1/2}  \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
(4\pi)^{-1/2} \biggl[ \frac{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )}{\sin\xi_e} \biggr] \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\frac{R_\mathrm{eq}}{R_\mathrm{SWS}}
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\biggl[ \frac{K}{2\pi G} \biggr]^{1/2} \xi_e \biggl[ \frac{G}{2K} \biggr]^{1/2} =
(4\pi)^{-1/2} \xi_e \, .
</math>
  </td>
</tr>
</table>
</div>

<span id="Stahler1983Fig17">
<div align="center">
<table border="2" cellpadding="8">
<tr>
  <td align="center" colspan="2">
'''Figure 2:''' <font color="darkblue">Equilibrium Mass-Radius Diagram </font> 
  </td>
</tr>
<tr>
  <td align="center" bgcolor="white">
[[File:Stahler1983TitlePage0.png|300px|center|Stahler (1983) Title Page]]
<!-- [[Image:AAAwaiting01.png|300px|center|Stahler (1983) Title Page]] -->
  </td>
  <td valign="top" width=350 rowspan="3">
''Top:'' A slightly edited reproduction of Figure 17 in association with Appendix B of  [http://adsabs.harvard.edu/abs/1983ApJ...268..165S Stahler] (1983, ApJ, 268, 165).  Stahler's figure caption reads, in part, "Mass-radius relation for bounded polytropes (schematic).  Each curve is labeled by the appropriate value or range" of {{Math/MP_PolytropicIndex}} &hellip; "As the cloud density increases from unity, all curves leave the origin with the same slope &hellip;"  


''Bottom:''  Curves depict the exact, analytically derived mass-radius relationship for truncated <math>~n = 1</math> (purple squares) and <math>~n = 5</math> (blue diamonds) polytropes that are embedded in an external medium of pressure <math>~P_e</math>; the relevant mathematical expressions are presented to the immediate right of Stahler's name in, respectively, our [[SSC/Structure/PolytropesEmbedded#n1Summary|<math>~n=1</math> summary table]] and our [[SSC/Structure/PolytropesEmbedded#n5Summary|<math>~n=5</math> summary table]].  As the dimensionless truncation radius, <math>~\xi_e</math>, increases steadily from zero, both curves exhibit very similar behavior up to <math>~M_n \equiv M/M_\mathrm{SWS} \approx 0.5</math>; thereafter the normalized mass and normalized radius continue to steadily increase along the <math>~n = 1</math> sequence, but the <math>~n = 5</math> sequence eventually bends back on itself, returning to the origin as <math>~\xi_e \rightarrow \infty</math>.  


''Comparison:'' The monotonic <math>P-R</math> behavior of the analytically derived solution for {{Math/MP_PolytropicIndex}} = 1 <math>(\gamma_g = 2)</math>, shown above, is consistent with the behavior of the numerically derived solutions presented by Whitworth for slightly lower values of <math>\gamma_g</math> = 5/3 and 4/3.  The analytically derived solution for {{Math/MP_PolytropicIndex}} = 5 <math>(\gamma_g = 6/5)</math> shows that, above some limiting pressure, no equilibrium configuration exists; this is consistent with the behavior of the numerically derived solutions presented by Whitworth for all values of <math>\gamma_g < 4/3 \, .</math>
 </td>
</tr>

<tr>
  <td align="center" bgcolor="white">
[[File:Stahler_MRdiagram1.png|300px|center|Stahler (1983) Figure 17 (edited)]]
<!-- [[Image:AAAwaiting01.png|300px|center|Stahler (1983) Figure 17 (edited)]] -->
  </td>
</tr>

<tr>
  <td align="center" bgcolor="white">
[[File:Stahler1983Comparison.png|300px|center|To be compared with Stahler (1983)]]
  </td>
</tr>
</table>
</div>
</span>

===Tabular Summary (n=1) ===
<span id="n1Summary">
<div align="center">
<table border="1" cellpadding="8" width="95%">
<tr>
  <th align="center" colspan="3">
Table 1: &nbsp;Properties of <math>~n=1</math> Polytropes Embedded in an External Medium of Pressure <math>~P_e</math> 
<br>
(and, accordingly, truncated at radius <math>~\xi_e</math>)
  </th>
</tr>
<tr>
  <td align="center" colspan="3">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>
~\theta_1 = \frac{\sin\xi_e}{\xi_e} 
</math>
  </td>

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

  <td align="right">
<math>
~\frac{d\theta_1}{d\xi} \biggr|_{\xi_e} = \frac{\cos\xi_e}{\xi_e} - \frac{\sin\xi_e}{\xi_e^2}
</math>
  </td>
</tr>
</table>
  </td>
</tr>

<tr>
  <td align="center" rowspan="1">
[http://adsabs.harvard.edu/abs/1970MNRAS.151...81H Horedt (1970)]
<br>for<br>
fixed <math>~(M,K_n)</math>
  </td>
  <td align="center">
<math>
~r_a = \frac{R_\mathrm{eq}}{R_\mathrm{Horedt}} = \xi_e
</math>
  </td>
  <td align="center">
<math>
~p_a = \frac{P_e}{P_\mathrm{Horedt}} = \biggl[ \frac{\sin\xi_e}{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )} \biggr]^2
</math>
  </td>
</tr>

<tr>
  <td align="center" rowspan="1">
[http://adsabs.harvard.edu/abs/1981MNRAS.195..967W Whitworth (1981)]
<br>for<br>
fixed <math>~(M,K_n)</math>
  </td>
  <td align="center">
<math>
\frac{R_\mathrm{eq}}{R_\mathrm{rf}} = \biggl( \frac{2^3}{3^2 \cdot 5} \biggr)^{1/2} \xi_e
</math>
  </td>
  <td align="center">
<math>
\frac{P_e}{P_\mathrm{rf}} = \biggl( \frac{3^4 \cdot 5^3}{2^7} \biggr) \biggl[ \frac{\sin\xi_e}{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )} \biggr]^2
</math>
  </td>
</tr>

<tr>
  <td align="center" rowspan="1">
[http://adsabs.harvard.edu/abs/1983ApJ...268..165S Stahler (1983)]
<br>for<br>
fixed <math>~(P_e,K_n)</math>
  </td>
  <td align="center">
<math>
\frac{R_\mathrm{eq}}{R_\mathrm{SWS}} = (4\pi)^{-1/2} \xi_e
</math>
  </td>
  <td align="center">
<math>
\frac{M}{M_\mathrm{SWS}} = (4\pi)^{-1/2} \biggl[ \frac{\xi_e(\sin\xi_e - \xi_e \cos\xi_e )}{\sin\xi_e} \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="left" colspan="3">
NOTE:  None of the analytic expressions for the dimensionless radius, pressure, or mass presented in this table explicitly appear in the referenced articles by Horedt, by Whitworth, or by Stahler but, as is discussed fully above, they are straightforwardly derivable from the more general relations that appear in these papers.  
  </td>
</tr>
</table>
</div>
</span>