Difference between revisions of "User:Tohline/SSC/VirialStability"

From VistrailsWiki
Jump to navigation Jump to search
(→‎Strategy: Create table of chosen dimensional scaling parameters)
 
(16 intermediate revisions by the same user not shown)
Line 3: Line 3:
=Virial Stability of BiPolytropes=
=Virial Stability of BiPolytropes=
{{LSU_HBook_header}}
{{LSU_HBook_header}}
==Discussion with Kundan Kadam==
See particularly … [https://ui.adsabs.harvard.edu/abs/2016MNRAS.462.2237K/abstract K. Kadam ''et al.'' (2016, MNRAS, 462, pp. 2237-2245)]
The comments provided inside the following box are for Kundan Kadam, providing feedback on a derivation he sent to me ( J. E. Tohline) on Monday, 27 January 2014.
<table align="center" border="1" width="95%" cellpadding="5">
<tr><td align="left">
[<font color="red">31 January 2014</font>] I agree with the following derived expressions; they all already appear in my presentation elsewhere on this page:
<div align="center">
Throughout the envelope &hellip; &nbsp; &nbsp; &nbsp;
<math>
M_r = \frac{4\pi}{3} \biggl( r^3\rho_e + r_i^3\rho_c - r_i^3\rho_e\biggr) \, ;
</math>
<math>
M_\mathrm{tot} = \frac{4\pi}{3}\rho_c \biggl( \frac{q^3 R^3}{\nu} \biggr) \, ;
</math>
<math>
\nu = \biggl[1 + \biggl(\frac{\rho_e}{\rho_c}\biggr)(q^{-3} - 1) \biggr]^{-1} \, .
</math>
</div>
During the recent December holiday break, I also used the hydrostatic balance relation to derive expressions for the pressure throughout the two-component (uniform density core + uniform density envelope) model, as you have done.  I agree completely with your derivation of the expressions for &#8212; and relationships between &#8212; <math>P_{cc}</math>, <math>P_{ci}</math>, and <math>P_{ei}</math>.  In addition, I chose to normalize all the pressures to,
<div align="center">
<math>
P_\mathrm{norm} \equiv \biggl(G^3 \rho_c^4 M_\mathrm{tot}^2 \biggr)^{1/3} \, .
</math>
</div>
Letting an asterisk superscript denote normalized pressures, that is, <math>P^* \equiv P/P_\mathrm{norm}</math>, my derived expressions are:
<div align="center">
<math>
P^*_{cc} = P^*_{ci} + \biggr[\frac{\pi \nu^2}{6} \biggr]^{1/3} \, ,
</math>
</div>
and,
<div align="center">
<math>
P^*_{ei} = P^*_{ci} = \biggr[\frac{\pi \nu^2}{6} \biggr]^{1/3} \biggl( \frac{\rho_e}{\rho_c}\biggr) \biggl[ \frac{\rho_e}{\rho_c}
\biggl(\frac{1}{q^2}-1\biggr) + 2 \biggl(1- \frac{\rho_e}{\rho_c}\biggr)(1-q) \biggr] \, .
</math>
</div>
<br />
<div align="center">
<hr width="85%" />
</div>
<br />
[<font color="red">4 February 2014</font>] In deriving the above analytic expression for the pressure, we have effectively determined the detailed [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|equilibrium structure of a bipolytrope that has <math>n_c = 0</math> and <math>n_e = 0</math>]] without ever using the virial theorem.  This is an excellent accomplishment for a variety of reasons. 
<ins>First</ins>, we can now compare and contrast these newly derived expressions with the analogous expressions that have been derived earlier for a [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|bipolytrope with  <math>n_c = 5</math> and <math>n_e = 1</math>]].  For example, at the end of [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#Step_7:__Surface_Boundary_Condition|Step #7 in the new derivation]], we find the following expressions for the equilibrium radius and total mass of a bipolytrope with <math>(n_c, n_e) = (0, 0)</math> expressed in terms of the chosen central density, <math>\rho_0</math>, and central pressure, <math>P_0</math>, of the configuration:
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\biggl[ \frac{G\rho_0^2}{P_0} \biggr]^{1/2} R</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp;
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \frac{\sqrt{3}}{q g} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>
\biggl[ \frac{G^3\rho_0^4}{P_0^3} \biggr]^{1/2} M_\mathrm{tot}</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp;
  </td>
  <td align="left">
<math>\biggl( \frac{2}{\pi} \biggr)^{1/2} \frac{\sqrt{3}}{\nu g^3} \, ,</math>
  </td>
</tr>
</table>
</div>
where the function (<font color="red">note that, for a few days, there was a sign error in this definition of <math>g</math></font>),
<div align="center">
<math>
g(q, \rho_e/\rho_0) \equiv \biggl\{ 1  +
\biggl(\frac{\rho_e}{\rho_0}\biggr)  \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl(1- q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] \biggr\}^{1/2}\, .
</math>
</div>
We will also find it useful to combine these two expressions to eliminate direct reference to the central density, <math>\rho_0</math>, obtaining,
<div align="center">
<table border="0">
<tr>
  <td align="right">
<math>
\biggl[ \frac{R^4}{GM_\mathrm{tot}^2} \biggr] P_0</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp;
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2^3\pi} \biggr) \frac{3\nu^2 g^2}{q^4} \, .</math>
  </td>
</tr>
</table>
</div>
For comparison, referring back to the "Normalization" discussion and table of "Parameter Values" provided in our [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Examples|example bipolytrope with <math>(n_c, n_e) = (5, 1)</math>]], we have,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\biggl( \frac{\mu_e}{\mu_c} \biggr) \biggl[ \frac{G\rho_0^2}{(K_c \rho_0^{6/5})} \biggr]^{1/2} R</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp;
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \frac{\eta_s}{\theta_i^2} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>
\biggl( \frac{\mu_e}{\mu_c} \biggr) ^2 \biggl[ \frac{G^3\rho_0^4}{(K_c \rho_0^{6/5})^3} \biggr]^{1/2} M_\mathrm{tot}</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp;
  </td>
  <td align="left">
<math>\biggl( \frac{2}{\pi} \biggr)^{1/2} \frac{A\eta_s}{\theta_i} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>
\biggl[ \frac{R^4}{GM^2_\mathrm{tot}} \biggr]^{1/2} K_c \rho_0^{6/5}</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp;
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2^3\pi} \biggr) \frac{\eta_s^2}{A^2\theta_i^6} \, ,</math>
  </td>
</tr>
</table>
</div>
where expressions for the various functions, <math>A, \eta_s, \theta_i</math>, are also provided in the table of "Parameter Values."  The dimensional normalizations are clearly the same in both types of bipolytropes because, in the latter case, <math>P_0 = K_c \rho_0^{6/5}</math>.
<ins>Second</ins>, a comparison between these two types of bipolytropes may help clarify how a discontinuous jump in the mean molecular weight should be handled in the more general virial analysis.  Usually <math>\rho_e/\rho_c</math> is set equal to <math>\mu_e/\mu_c</math> at the interface, but if the core and envelope are both treated as having uniform densities, setting the ratio of densities to the ratio of mean molecular weights seems to overconstrain the problem.  I haven't figured out yet how to handle this, but maybe this model comparison will help.
<ins>Third</ins>, when conducting the virial analysis, it should now be clear how to evaluate the system's free energy.  When the system is treated as being composed of a central uniform-density spherical core surrounded by a uniform-density envelope, we have already derived an analytic expression for the total gravitational potential energy, <math>W</math>.  The above derivation of the pressure distribution throughout such a configuration &#8212; that is, throughout a bipolytrope with <math>(n_c, n_e) = (0, 0)</math> &#8212; now lets us derive an analytic expression for the thermal energy of the core, <math>S_\mathrm{core}</math>, and an analytic expression for the thermal energy of the envelope, <math>S_\mathrm{env}</math>.  Whether focusing on the core or the envelope, start with an appreciation that,
<div align="center">
<math>
dS = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) dm = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) \biggl( \frac{4\pi}{3} \biggr) \rho r^2 dr
= 2\pi Pr^2 dr \, .
</math>
</div>
Hence,
<div align="center">
<math>
S_\mathrm{core} = 2\pi \int_0^{r_i} [P(r)]_\mathrm{core} r^2 dr; ~~~~~\mathrm{and}~~~~~
S_\mathrm{env} = 2\pi \int_{r_i}^R [P(r)]_\mathrm{env} r^2 dr \, .
</math>
</div>
Then, for either volume segment, the relationship between the internal energy, <math>U</math>, and the thermal energy, <math>S</math>, is,
<div align="center">
<math>
U = \biggl[ \frac{2}{3(\gamma - 1)} \biggr] S \, .
</math>
</div>
My initial derivation gives,
<table border="1" align="center" width="85%">
<tr><td align="left">
<font color="darkblue">
Do not use these expressions as they are flawed.  The correct expressions follow in the comments dated 12 February 2014.
</font>
</td></tr>
<tr><td>
<math>
S_\mathrm{core} = 2\pi P_0 R^3 \biggl[q^3 - q^5 \biggl( \frac{2\pi G\rho_0^2 R^2}{5P_0} \biggr) \biggr] \, ;
</math>
</td></tr>
<tr><td>
<math>
S_\mathrm{env} = \frac{4\pi}{3} R^3 \biggl\{ \frac{3P_i}{2} (1-q^3) + \pi G \rho_0^2 R^2 \biggl(\frac{\rho_e}{\rho_c} \biggr)
\biggl[  \frac{1}{6}(6q^3 - q^2 - 5q^5) - \frac{1}{10}\biggl(\frac{\rho_e}{\rho_c} \biggr) (2 + 10q^3 - 5q^2-7q^5)      \biggr] \biggr\} \, .
</math>
</td></tr>
</table>
<br />
<div align="center">
<hr width="85%" />
</div>
<br />
[<font color="red">12 February 2014</font>] In an [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#Thermal_Energy_Content|accompanying discussion of the thermal energy content of an <math>(n_c, n_e) = (0, 0)</math> bipolytrope]], I've re-derived and cross-checked the expressions for <math>S_\mathrm{core}</math> and <math>S_\mathrm{env}</math>.  I am quite confident that the correct expressions are:
<div align="center">
<table border="0">
<tr>
  <td align="right">
<math>~S_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>\biggl( \frac{4\pi}{5} \biggr) R^3 q^5 \biggl (\frac{5P_i}{2q^2} + \Pi \biggr)  \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~S_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\frac{\pi R^3}{5} \biggl\{ 10 P_i  (1-q^3) 
+ 10 \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggl[-2q^2 + 3q^3 - q^5  \biggr]
  + 6 \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr)^2 \biggl[ -1 + 5q^2 -5q^3 +  q^5  \biggr] \biggr\} \, ,
</math>
  </td>
</tr>
</table>
</div>
where, <math>P_i</math> is the pressure at the interface,
<div align="center">
<math>
\Pi \equiv \frac{3}{2^3 \pi}
\biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) \, ,
</math>
</div>
and the relationship between the central pressure and the pressure at the interface is,
<div align="center">
<math>~P_0 = P_i + \Pi q^2 \, .</math>
</div>
<br />
<div align="center">
<hr width="85%" />
</div>
<br />
[<font color="red">7 February 2014</font>] More succinctly, this <math>(n_c, n_e) = (0,0)</math> bipolytrope is very interesting because it provides three independent ways of deriving the proper equilibrium structure.  If we work through the detailed derivation in all three cases, we should have a very firm foundation from which to perform the stability analysis of a broad range of bipolytropic configurations. The three derivations are &hellip;
<ol>
<li>Detailed force balance (Kundan has already done this):  Use the hydrostatic balance equation to determine how the pressure varies with radius throughout the two-component (core-envelope) structure.  For a given choice of <math>P_0, \rho_0, \nu,</math> and <math>q</math>, the radius and, hence, total mass of the configuration is determined by the surface boundary condition, that is, <math>P = 0</math> at <math>r = R</math>.  These results are summarized above.
<li>Virial equilibrium:  The structures derived from demanding detailed force balance should also be in virial equilibrium.  (The virial equilibrium condition is less demanding than detailed force balance because it is a global statement of energy balance and is not concerned with internal structural details. Virial equilibrium is a necessary but not sufficient condition; the force-balance solution gives the more detailed information but must, itself, also define a structure that has virial energy balance.)  It is well known that, for the simplest spherically symmetric stars, the statement of virial balance is,
<div align="center">
<math>
2S + W = 0 \, .
</math>
</div>
Since we have analytic expressions for both <math>S</math> and <math>W</math>, we should check to see if they sum appropriately to zero as expected from this simple statement of virial equilibrium.  [Note inserted <font color="red">12 February 2014</font>:  An [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#Virial_Equilibrium|accompanying discussion of the virial equilibrium of this particular bipolytrope]] shows that the relation does hold precisely.]
</ol>
<ol start="3">
<li>Extrema in the Free Energy:  We should also expect to be able to define the free energy of this bipolytrope by the expression,
<div align="center">
<math>
\mathfrak{G} = U + W \, .
</math>
</div>
Then, for a given choice of <math>P_0, \rho_0, \nu,</math> and <math>q</math> (perhaps we should choose <math>M_\mathrm{tot}</math> instead of <math>P_0</math>), the equilibrium radius should be defined by the radius at which,
<div align="center">
<math>
\frac{\partial \mathfrak{G}}{\partial R} = 0 \, .
</math>
</div>
We should check to make sure that this is also true.
</ol>
Items #2 and #3 are the next things I will be working on.
<br />
<div align="center">
<hr width="85%" />
</div>
<br />
[<font color="red">17 February 2014; by J. E. Tohline</font>] Using derivatives of the free energy function, <math>~\mathfrak{G}</math>, I have finished deriving the equilibrium radius for bipolytropic configurations where the density of the core <math>~(\rho_0)</math> and the density of the envelope <math>~(\rho_e)</math> are both assumed to be uniform.  I have also finished deriving the condition for the dynamical stability of such configurations assuming that, as the equilibrium structure undergoes a radial perturbation, the core compresses along a <math>~\gamma_c</math> adiabat while the envelope compresses along a <math>~\gamma_e</math> adiabat.  The [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#Equilibrium_Condition|equilibrium radius is derived here]] and the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic0_0#Stability_Condition|stability condition is described here]].
<br />
<div align="center">
<hr width="85%" />
</div>
<br />
[<font color="red">25 February 2014; by J. E. Tohline</font>] I have extended the free-energy-based stability analysis to models of <math>~(n_c, n_e) = (5, 1)</math> bipolytropes.  As is [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Expression_for_Free_Energy|detailed here]], I have been able to derive analytic, closed form expressions for the gravitational potential energy and internal thermal energies of these configurations.  As is [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Equilibrium_Condition|shown here]], an examination of the first derivative of the free energy identifies equilibrium configurations that exactly match those derived from the earlier detailed force-balance derivation; and, finally, as [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Stability_Condition|shown here (see especially Figure 3)]], the second derivative of the free energy expression allows us to determine which equilibrium models are stable and which are dynamically unstable.  Hooray!
</td></tr>
</table>
{{LSU_WorkInProgress}}


=BiPolytrope Structural Relations=
=BiPolytrope Structural Relations=
Line 290: Line 596:


=Strategy=
=Strategy=
==Scaling Parameters==
We want to vary the total radius, <math>R \,</math>, of the configuration and look for extrema in the free energy, while holding the following parameters fixed:  <math>M_\mathrm{tot} \,</math>, <math>c_s^2 \,</math> (or <math>K_c \,</math>), <math>\nu \,</math>, and <math>q \,</math>.  So we need to rewrite the expressions for <math>W \,</math> and <math>U \,</math> such that everything is constant except for <math>R \,</math>, or <math>\chi \equiv R/R_0 \,</math>.  And, in order to put everything explicitly in terms of the fixed parameters just specified, let's go ahead and define the length, time, and density against which all dimensional quantities will ultimately be scaled.
We want to vary the total radius, <math>R \,</math>, of the configuration and look for extrema in the free energy, while holding the following parameters fixed:  <math>M_\mathrm{tot} \,</math>, <math>c_s^2 \,</math> (or <math>K_c \,</math>), <math>\nu \,</math>, and <math>q \,</math>.  So we need to rewrite the expressions for <math>W \,</math> and <math>U \,</math> such that everything is constant except for <math>R \,</math>, or <math>\chi \equiv R/R_0 \,</math>.  And, in order to put everything explicitly in terms of the fixed parameters just specified, let's go ahead and define the length, time, and density against which all dimensional quantities will ultimately be scaled.


Line 350: Line 657:




==Energy Expressions==
==Rescaled Energy Expressions==
The expression for <math>W</math> shown above needs only a slight modification to put it in the appropriate form, namely,
===Isothermal Core===
In the case of an isothermal core, the expression for <math>W</math> shown above needs only a slight modification to put it in the appropriate form, namely,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


Line 418: Line 726:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>\mathfrak{G}</math>
<math>\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{M_\mathrm{tot} c_s^2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 425: Line 733:
   <td align="left">
   <td align="left">
<math>
<math>
  - \frac{3}{5}  M_\mathrm{tot} c_s^2 \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \chi^{-1}  
  - \frac{3}{5}  \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \chi^{-1}  
+ \nu M_\mathrm{tot} c_s^2 [\ln(\nu/q^3) -3 \ln\chi]  + (1-\nu) n_e M_\mathrm{tot} K_e \rho_0^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \chi^{-3/n_e}
+ \nu [\ln(\nu/q^3) -3 \ln\chi]  + (1-\nu) n_e \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \chi^{-3/n_e}
</math>
</math>
   </td>
   </td>
Line 440: Line 748:
   <td align="left">
   <td align="left">
<math>
<math>
\mathfrak{G}_0 - A \chi^{-1} - C \ln\chi  + D \chi^{-3/n_e} \, ,
\mathfrak{G}_0^* - A \chi^{-1} - C \ln\chi  + D \chi^{-3/n_e} \, ,
</math>
</math>
   </td>
   </td>
Line 451: Line 759:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>\mathfrak{G}_0 \,</math>
<math>\mathfrak{G}_0^* \,</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 458: Line 766:
   <td align="left">
   <td align="left">
<math>
<math>
\nu M_\mathrm{tot} c_s^2 \ln(\nu/q^3)  \, ,
\nu \ln(\nu/q^3)  \, ,
</math>
</math>
   </td>
   </td>
Line 472: Line 780:
   <td align="left">
   <td align="left">
<math>
<math>
\frac{3}{5} M_\mathrm{tot} c_s^2 \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, ,
\frac{3}{5} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, ,
</math>
</math>
   </td>
   </td>
Line 486: Line 794:
   <td align="left">
   <td align="left">
<math>
<math>
3 \nu M_\mathrm{tot} c_s^2 \, ,
3 \nu  \, ,
</math>
</math>
   </td>
   </td>
Line 500: Line 808:
   <td align="left">
   <td align="left">
<math>
<math>
(1-\nu) n_e M_\mathrm{tot} K_e \rho_0^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e}  
(1-\nu) n_e \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e}  
= (1-\nu) n_e M_\mathrm{tot} K_e \biggl[  \frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr)  \biggr]^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \, .
= (1-\nu) n_e \biggl( \frac{K_e}{c_s^2} \biggr) \biggl[  \frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr)  \biggr]^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \, .
</math>
  </td>
</tr>
</table>
 
===Polytropic Core===
In the case of a polytropic core,
 
<table border="0" align="center" cellpadding="5">
<tr>
  <td align="right">
<math>\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{M_\mathrm{tot} K_c \rho_0^{1/n_c}}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
- A \chi^{-1}  + B_c \chi^{-3/n_c} + B_e \chi^{-3/n_e} \, ,
</math>
  </td>
</tr>
 
</table>
where,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>A \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
\frac{3}{5} \biggl[ \frac{GM_\mathrm{tot}}{R_0 K_c \rho_0^{1/n_c}} \biggr] \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) =
\frac{3}{5} \biggl( \frac{4\pi}{3} \biggr)^{1/n_c} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>B_c \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_c \nu \biggl( \frac{\nu}{q^3}\biggr)^{1/n_c} \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>B_e \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_e (1-\nu) \biggl[ \biggl( \frac{K_e }{K_c}\biggr)\rho_0^{1/n_e - 1/n_c} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e}
\, .
</math>
</math>
   </td>
   </td>
Line 518: Line 893:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>\biggl[ \frac{P_c}{(\rho_c/\mu_c)} \biggr]_E</math>
<math>\biggl[ \frac{P_c}{(\rho_c/\mu_c)} \biggr]_E \, .</math>
  </td>
</tr>
</table>
Note that this reflects the same physical condition as the constraint that is placed on the enthalpy at the interface when we analyze the detailed structure of <math>n_c = 5</math>, <math>n_e=1</math> bipolytropes (see the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Parameter_Values|middle of Table 1 in the accompanying discussion]]), namely,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\biggl[ \frac{\mu_e H_e}{(n_e+1)}\biggr]_i</math>
  </td>
  <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl[ \frac{\mu_c H_c}{(n_c+1)} \biggr]_i \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>


In the case of an isothermal core, this implies,
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ c_s^2}</math>
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ c_s^2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 536: Line 929:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e \rho_0^{1/n_e}}{ c_s^2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 542: Line 935:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>\rho_0^{-1/n_e} \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>
<math> \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>
   </td>
   </td>
</tr>
</tr>
Line 555: Line 948:
   <td align="left">
   <td align="left">
<math>
<math>
\biggl\frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr)   \biggr]^{-1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e}
\biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e}
\biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e}
= \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \chi_E^{3/n_e}
</math>
</math>
  </td>
</tr>
</table>
<table border="1" cellpadding="5" align="center">
<tr>
  <td colspan="3" align="center">
Not Necessarily Useful
   </td>
   </td>
</tr>
</tr>
Line 597: Line 998:
   </td>
   </td>
</tr>
</tr>
</table>
Therefore, the coefficient <math>D\,</math> becomes,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>D \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_e (1-\nu) \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e}
= n_e (1-\nu) \biggl( \frac{\mu_c }{\mu_e}\biggr) \chi_E^{3/n_e} \, .
</math>
  </td>
</tr>
</table>
In the case of a polytropic core, this implies,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ K_c}</math>
  </td>
  <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\rho_0^{1/n_c} \biggl[\rho_e^{-1/n_e}\biggr]_E</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl( \frac{K_e }{ K_c} \biggr) \rho_0^{(1/n_e-1/n_c)}</math>
  </td>
  <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math> \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e}
= \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \chi_E^{3/n_e}
</math>
  </td>
</tr>
</table>
Therefore, the coefficient <math>B_e\,</math> becomes,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>B_e \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_e (1-\nu) \biggl( \frac{\mu_c }{\mu_e}\biggr) \chi_E^{3/n_e} \, .
</math>
  </td>
</tr>
</table>
</table>



Latest revision as of 21:50, 17 November 2020


Virial Stability of BiPolytropes

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

Discussion with Kundan Kadam

See particularly … K. Kadam et al. (2016, MNRAS, 462, pp. 2237-2245)

The comments provided inside the following box are for Kundan Kadam, providing feedback on a derivation he sent to me ( J. E. Tohline) on Monday, 27 January 2014.

[31 January 2014] I agree with the following derived expressions; they all already appear in my presentation elsewhere on this page:

Throughout the envelope …       <math> M_r = \frac{4\pi}{3} \biggl( r^3\rho_e + r_i^3\rho_c - r_i^3\rho_e\biggr) \, ; </math>

<math> M_\mathrm{tot} = \frac{4\pi}{3}\rho_c \biggl( \frac{q^3 R^3}{\nu} \biggr) \, ; </math>

<math> \nu = \biggl[1 + \biggl(\frac{\rho_e}{\rho_c}\biggr)(q^{-3} - 1) \biggr]^{-1} \, . </math>

During the recent December holiday break, I also used the hydrostatic balance relation to derive expressions for the pressure throughout the two-component (uniform density core + uniform density envelope) model, as you have done. I agree completely with your derivation of the expressions for — and relationships between — <math>P_{cc}</math>, <math>P_{ci}</math>, and <math>P_{ei}</math>. In addition, I chose to normalize all the pressures to,

<math> P_\mathrm{norm} \equiv \biggl(G^3 \rho_c^4 M_\mathrm{tot}^2 \biggr)^{1/3} \, . </math>

Letting an asterisk superscript denote normalized pressures, that is, <math>P^* \equiv P/P_\mathrm{norm}</math>, my derived expressions are:

<math> P^*_{cc} = P^*_{ci} + \biggr[\frac{\pi \nu^2}{6} \biggr]^{1/3} \, , </math>

and,

<math> P^*_{ei} = P^*_{ci} = \biggr[\frac{\pi \nu^2}{6} \biggr]^{1/3} \biggl( \frac{\rho_e}{\rho_c}\biggr) \biggl[ \frac{\rho_e}{\rho_c} \biggl(\frac{1}{q^2}-1\biggr) + 2 \biggl(1- \frac{\rho_e}{\rho_c}\biggr)(1-q) \biggr] \, . </math>




[4 February 2014] In deriving the above analytic expression for the pressure, we have effectively determined the detailed equilibrium structure of a bipolytrope that has <math>n_c = 0</math> and <math>n_e = 0</math> without ever using the virial theorem. This is an excellent accomplishment for a variety of reasons.

First, we can now compare and contrast these newly derived expressions with the analogous expressions that have been derived earlier for a bipolytrope with <math>n_c = 5</math> and <math>n_e = 1</math>. For example, at the end of Step #7 in the new derivation, we find the following expressions for the equilibrium radius and total mass of a bipolytrope with <math>(n_c, n_e) = (0, 0)</math> expressed in terms of the chosen central density, <math>\rho_0</math>, and central pressure, <math>P_0</math>, of the configuration:

<math> \biggl[ \frac{G\rho_0^2}{P_0} \biggr]^{1/2} R</math>

  <math>~=</math> 

<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \frac{\sqrt{3}}{q g} \, ,</math>

<math> \biggl[ \frac{G^3\rho_0^4}{P_0^3} \biggr]^{1/2} M_\mathrm{tot}</math>

  <math>~=</math> 

<math>\biggl( \frac{2}{\pi} \biggr)^{1/2} \frac{\sqrt{3}}{\nu g^3} \, ,</math>

where the function (note that, for a few days, there was a sign error in this definition of <math>g</math>),

<math> g(q, \rho_e/\rho_0) \equiv \biggl\{ 1 + \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl(1- q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] \biggr\}^{1/2}\, . </math>

We will also find it useful to combine these two expressions to eliminate direct reference to the central density, <math>\rho_0</math>, obtaining,

<math> \biggl[ \frac{R^4}{GM_\mathrm{tot}^2} \biggr] P_0</math>

  <math>~=</math> 

<math>\biggl( \frac{1}{2^3\pi} \biggr) \frac{3\nu^2 g^2}{q^4} \, .</math>

For comparison, referring back to the "Normalization" discussion and table of "Parameter Values" provided in our example bipolytrope with <math>(n_c, n_e) = (5, 1)</math>, we have,

<math> \biggl( \frac{\mu_e}{\mu_c} \biggr) \biggl[ \frac{G\rho_0^2}{(K_c \rho_0^{6/5})} \biggr]^{1/2} R</math>

  <math>~=</math> 

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

<math> \biggl( \frac{\mu_e}{\mu_c} \biggr) ^2 \biggl[ \frac{G^3\rho_0^4}{(K_c \rho_0^{6/5})^3} \biggr]^{1/2} M_\mathrm{tot}</math>

  <math>~=</math> 

<math>\biggl( \frac{2}{\pi} \biggr)^{1/2} \frac{A\eta_s}{\theta_i} \, ,</math>

<math> \biggl[ \frac{R^4}{GM^2_\mathrm{tot}} \biggr]^{1/2} K_c \rho_0^{6/5}</math>

  <math>~=</math> 

<math>\biggl( \frac{1}{2^3\pi} \biggr) \frac{\eta_s^2}{A^2\theta_i^6} \, ,</math>

where expressions for the various functions, <math>A, \eta_s, \theta_i</math>, are also provided in the table of "Parameter Values." The dimensional normalizations are clearly the same in both types of bipolytropes because, in the latter case, <math>P_0 = K_c \rho_0^{6/5}</math>.

Second, a comparison between these two types of bipolytropes may help clarify how a discontinuous jump in the mean molecular weight should be handled in the more general virial analysis. Usually <math>\rho_e/\rho_c</math> is set equal to <math>\mu_e/\mu_c</math> at the interface, but if the core and envelope are both treated as having uniform densities, setting the ratio of densities to the ratio of mean molecular weights seems to overconstrain the problem. I haven't figured out yet how to handle this, but maybe this model comparison will help.

Third, when conducting the virial analysis, it should now be clear how to evaluate the system's free energy. When the system is treated as being composed of a central uniform-density spherical core surrounded by a uniform-density envelope, we have already derived an analytic expression for the total gravitational potential energy, <math>W</math>. The above derivation of the pressure distribution throughout such a configuration — that is, throughout a bipolytrope with <math>(n_c, n_e) = (0, 0)</math> — now lets us derive an analytic expression for the thermal energy of the core, <math>S_\mathrm{core}</math>, and an analytic expression for the thermal energy of the envelope, <math>S_\mathrm{env}</math>. Whether focusing on the core or the envelope, start with an appreciation that,

<math> dS = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) dm = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) \biggl( \frac{4\pi}{3} \biggr) \rho r^2 dr = 2\pi Pr^2 dr \, . </math>

Hence,

<math> S_\mathrm{core} = 2\pi \int_0^{r_i} [P(r)]_\mathrm{core} r^2 dr; ~~~~~\mathrm{and}~~~~~ S_\mathrm{env} = 2\pi \int_{r_i}^R [P(r)]_\mathrm{env} r^2 dr \, . </math>

Then, for either volume segment, the relationship between the internal energy, <math>U</math>, and the thermal energy, <math>S</math>, is,

<math> U = \biggl[ \frac{2}{3(\gamma - 1)} \biggr] S \, . </math>

My initial derivation gives,

Do not use these expressions as they are flawed. The correct expressions follow in the comments dated 12 February 2014.

<math> S_\mathrm{core} = 2\pi P_0 R^3 \biggl[q^3 - q^5 \biggl( \frac{2\pi G\rho_0^2 R^2}{5P_0} \biggr) \biggr] \, ; </math>

<math> S_\mathrm{env} = \frac{4\pi}{3} R^3 \biggl\{ \frac{3P_i}{2} (1-q^3) + \pi G \rho_0^2 R^2 \biggl(\frac{\rho_e}{\rho_c} \biggr) \biggl[ \frac{1}{6}(6q^3 - q^2 - 5q^5) - \frac{1}{10}\biggl(\frac{\rho_e}{\rho_c} \biggr) (2 + 10q^3 - 5q^2-7q^5) \biggr] \biggr\} \, . </math>




[12 February 2014] In an accompanying discussion of the thermal energy content of an <math>(n_c, n_e) = (0, 0)</math> bipolytrope, I've re-derived and cross-checked the expressions for <math>S_\mathrm{core}</math> and <math>S_\mathrm{env}</math>. I am quite confident that the correct expressions are:

<math>~S_\mathrm{core}</math>

<math>~=~</math>

<math>\biggl( \frac{4\pi}{5} \biggr) R^3 q^5 \biggl (\frac{5P_i}{2q^2} + \Pi \biggr) \, ,</math>

<math>~S_\mathrm{env}</math>

<math>~=~</math>

<math> \frac{\pi R^3}{5} \biggl\{ 10 P_i (1-q^3) + 10 \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggl[-2q^2 + 3q^3 - q^5 \biggr]

 + 6 \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr)^2 \biggl[ -1 + 5q^2 -5q^3 +  q^5   \biggr] \biggr\} \, ,

</math>

where, <math>P_i</math> is the pressure at the interface,

<math> \Pi \equiv \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) \, , </math>

and the relationship between the central pressure and the pressure at the interface is,

<math>~P_0 = P_i + \Pi q^2 \, .</math>




[7 February 2014] More succinctly, this <math>(n_c, n_e) = (0,0)</math> bipolytrope is very interesting because it provides three independent ways of deriving the proper equilibrium structure. If we work through the detailed derivation in all three cases, we should have a very firm foundation from which to perform the stability analysis of a broad range of bipolytropic configurations. The three derivations are …

  1. Detailed force balance (Kundan has already done this): Use the hydrostatic balance equation to determine how the pressure varies with radius throughout the two-component (core-envelope) structure. For a given choice of <math>P_0, \rho_0, \nu,</math> and <math>q</math>, the radius and, hence, total mass of the configuration is determined by the surface boundary condition, that is, <math>P = 0</math> at <math>r = R</math>. These results are summarized above.
  2. Virial equilibrium: The structures derived from demanding detailed force balance should also be in virial equilibrium. (The virial equilibrium condition is less demanding than detailed force balance because it is a global statement of energy balance and is not concerned with internal structural details. Virial equilibrium is a necessary but not sufficient condition; the force-balance solution gives the more detailed information but must, itself, also define a structure that has virial energy balance.) It is well known that, for the simplest spherically symmetric stars, the statement of virial balance is,

    <math> 2S + W = 0 \, . </math>

    Since we have analytic expressions for both <math>S</math> and <math>W</math>, we should check to see if they sum appropriately to zero as expected from this simple statement of virial equilibrium. [Note inserted 12 February 2014: An accompanying discussion of the virial equilibrium of this particular bipolytrope shows that the relation does hold precisely.]

  1. Extrema in the Free Energy: We should also expect to be able to define the free energy of this bipolytrope by the expression,

    <math> \mathfrak{G} = U + W \, . </math>

    Then, for a given choice of <math>P_0, \rho_0, \nu,</math> and <math>q</math> (perhaps we should choose <math>M_\mathrm{tot}</math> instead of <math>P_0</math>), the equilibrium radius should be defined by the radius at which,

    <math> \frac{\partial \mathfrak{G}}{\partial R} = 0 \, . </math>

    We should check to make sure that this is also true.

Items #2 and #3 are the next things I will be working on.




[17 February 2014; by J. E. Tohline] Using derivatives of the free energy function, <math>~\mathfrak{G}</math>, I have finished deriving the equilibrium radius for bipolytropic configurations where the density of the core <math>~(\rho_0)</math> and the density of the envelope <math>~(\rho_e)</math> are both assumed to be uniform. I have also finished deriving the condition for the dynamical stability of such configurations assuming that, as the equilibrium structure undergoes a radial perturbation, the core compresses along a <math>~\gamma_c</math> adiabat while the envelope compresses along a <math>~\gamma_e</math> adiabat. The equilibrium radius is derived here and the stability condition is described here.




[25 February 2014; by J. E. Tohline] I have extended the free-energy-based stability analysis to models of <math>~(n_c, n_e) = (5, 1)</math> bipolytropes. As is detailed here, I have been able to derive analytic, closed form expressions for the gravitational potential energy and internal thermal energies of these configurations. As is shown here, an examination of the first derivative of the free energy identifies equilibrium configurations that exactly match those derived from the earlier detailed force-balance derivation; and, finally, as shown here (see especially Figure 3), the second derivative of the free energy expression allows us to determine which equilibrium models are stable and which are dynamically unstable. Hooray!



Work-in-progress.png

Material that appears after this point in our presentation is under development and therefore
may contain incorrect mathematical equations and/or physical misinterpretations.
|   Go Home   |


BiPolytrope Structural Relations

[Following a discussion that Tohline had with Kundan Kadam on 3 July 2013, we have decided to carry out a virial equilibrium and stability analysis of nonrotating bipolytropes.]

We will adopt the following approach:

  • Properties of the core <math>\cdots</math>
    • Uniform density, <math>\rho_c</math>;
    • Polytropic constant, <math>K_c</math>, and polytropic index, <math>n_c</math>;
    • Surface of the core at <math>r_i</math>;
  • Properties of the envelope <math>\cdots</math>
    • Uniform density, <math>\rho_e</math>;
    • Polytropic constant, <math>K_e</math>, and polytropic index, <math>n_e</math>;
    • Base of the core at <math>r_i</math> and surface at <math>R</math>.

Use the dimensionless radius,

<math>\xi \equiv \frac{r}{r_i}</math>.

Then, <math>\xi_i = 1</math> and <math>\xi_s \equiv R/r_i</math>.

Expressions for Mass

Inside the core, the expression for the mass interior to any radius, <math>0 \le \xi \le 1</math>, is,

<math>M_\xi = \frac{4\pi}{3} \rho_c r_i^3 \xi^3</math> .

The expression for the mass interior to any position within the envelope, <math>1 \le \xi \le \xi_s</math>, is,

<math>M_\xi = \frac{4\pi}{3} r_i^3 \biggl[\rho_c + \rho_e(\xi^3 - 1) \biggr]</math> .

Hence, the mass of the core, the mass of the envelope, and the total mass are, respectively,

<math> M_\mathrm{core} = \frac{4\pi}{3} \rho_c r_i^3 </math> ;

<math> M_\mathrm{env} = \frac{4\pi}{3} r_i^3 \biggl[\rho_e (\xi_s^3 - 1) \biggr] </math> ;

<math> M_\mathrm{tot} = \frac{4\pi}{3} r_i^3 \biggl[\rho_c + \rho_e(\xi_s^3 - 1) \biggr] </math> .

Following the work of Schönberg & Chandrasekhar (1942) — see our accompanying discussion — we are seeking equilibrium configurations in the <math>\nu - q</math> plane where,

<math>\nu \equiv \frac{M_\mathrm{core}}{M_\mathrm{tot}} </math>     and     <math>q \equiv \frac{r_i}{R} = \frac{1}{\xi_s}</math>.

So we can combine the above expressions to obtain,

<math>\frac{\rho_e}{\rho_c} = \frac{M_\mathrm{env}}{M_\mathrm{core}} (\xi_s^3 - 1)^{-1} = \biggl[ \frac{1-\nu}{\nu}\biggr] (\xi_s^3 - 1)^{-1} = \frac{q^3}{\nu}\biggl( \frac{1 - \nu}{1- q^3} \biggr) \, , </math>

or,

<math>\nu = \biggl[ 1 + \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^3} - 1\biggr) \biggr]^{-1}</math> .

It is worth noting that exactly the same result arises from an examination of the analytically definable, structural properties of bipolytropes having <math>n_c = 5</math> and <math>n_e = 1</math>. That is, the ratio of the average density in the envelope to the average density in the core is,

<math> \frac{\bar{\rho}_e}{\bar{\rho}_c} = \frac{(M^*_\mathrm{tot} - M^*_\mathrm{core}) / [(R^*)^3 - (r^*_\mathrm{core})^3]}{M^*_\mathrm{core} / (r^*_\mathrm{core})^3 } = \frac{q^3 (1- \nu)}{\nu (1-q^3)} \, . </math>

This is, of course, at it should be.


It is worth noting that, because <math>\bar\rho \equiv 3M_\mathrm{tot}/(4\pi R^3) \,</math>, we can write,

<math>\frac{\rho_c}{\bar\rho} = \frac{\nu}{q^3} \, ,</math>    and     <math>\frac{\rho_e}{\bar\rho} = \frac{1-\nu}{1-q^3} \, ,</math>

which is consistent with the above expression for the ratio, <math>\rho_e/\rho_c \, .</math> The following figure shows how <math>\nu \,</math> varies with <math>q \,</math> for various choices of the mass density ratio, <math>\rho_e/\rho_c \,</math>. It illustrates that, for a given core-to-total mass ratio, <math>\nu \,</math>, the relative location of the interface radius, <math>q \,</math>, can vary between zero and one, but each value of <math>q \,</math> reflects a different ratio of envelope-to-core mass density.

Nu versus Q

Energy Expressions

The gravitational potential energy of the bipolytropic configuration is obtained by integrating over the following differential energy contribution,

<math>dW = - \biggl( \frac{GM_r}{r} \biggr) dm</math> .

Hence,

<math>W = W_\mathrm{core} + W_\mathrm{env}</math>

<math> = - G \biggl\{ \int_0^{r_i} \biggl( \frac{M_r}{r} \biggr) 4\pi r^2 \rho_c dr + \int^R_{r_i} \biggl( \frac{M_r}{r} \biggr) 4\pi r^2 \rho_e dr \biggr\} </math>

 

<math> = - G \biggl\{ \int_0^1 \biggl( \frac{4\pi }{3} \rho_c r_i^3 \xi^3 \biggr) 4\pi r_i^2 \rho_c \xi d\xi + \int_1^{\xi_s} \frac{4\pi}{3} \rho_c r_i^3 \biggl[ 1 + \frac{\rho_e}{\rho_c}(\xi^3 - 1) \biggr] 4\pi r_i^2 \rho_e \xi d\xi \biggr\} </math>

 

<math> = - \frac{3GM^2_\mathrm{core}}{r_i} \biggl\{ \int_0^1 \xi^4 d\xi + \int_1^{\xi_s} \biggl[ 1 + \frac{\rho_e}{\rho_c}(\xi^3 - 1) \biggr] \biggl( \frac{\rho_e}{\rho_c} \biggr) \xi d\xi \biggr\} </math>

 

<math> = - \frac{3GM^2_\mathrm{core}}{r_i} \biggl\{ \frac{1}{5} + \biggl( \frac{\rho_e}{\rho_c} \biggr) \int_1^{\xi_s} \xi d\xi + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \int_1^{\xi_s} (\xi^3 - 1) \xi d\xi \biggr\} </math>

 

<math> = - \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) 3\nu^2 \xi_s \biggl\{ \frac{1}{5} + \frac{1}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \frac{1}{5}(\xi_s^5 - 1) - \frac{1}{2}(\xi_s^2-1) \biggr] \biggr\} </math>

 

<math> = - \biggl( \frac{3GM^2_\mathrm{tot}}{5R_0} \biggr) \biggl( \frac{R}{R_0}\biggr)^{-1} \nu^2 \xi_s f(\nu,q) \, , </math>

where <math>R_0</math> is an, as yet unspecified, normalization radius, and

<math> f(\nu,q) \equiv 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \, . </math>

I like the form of this expression. The leading term, which scales as <math>R^{-1}</math>, encapsulates the behavior of the gravitational potential energy for a given choice of the internal structure, namely, a given choice of <math>\xi_s</math>, <math>\nu</math>, and density ratio <math>(\rho_e/\rho_c)</math>. Actually, only two internal structural parameters need to be specified — <math>\nu</math> and <math>\xi_s</math> (or, <math>q</math>). From these two, the expression shown above allows the determination of <math>(\rho_e/\rho_c)</math>.

Sanity Check: Uniform Density Configuration

<math> M_\mathrm{core} = \frac{4\pi}{3} r_i^3 \rho_c \, ; </math>

<math> M_\mathrm{env} = \frac{4\pi}{3} r_i^3 \rho_e (\xi_s^3 -1) \, ; </math>

<math> M_\mathrm{tot} = \frac{4\pi}{3} r_i^3 \rho_c \biggl[ 1 + \frac{\rho_e}{\rho_c} (\xi_s^3 -1) \biggr] \, . </math>

If <math>\rho_e/\rho_c = 1</math>, then,

<math> \frac{M_\mathrm{env}}{M_\mathrm{core}} = \biggl(\frac{1}{\nu} - 1\biggr) = (\xi_s^3 -1) = \biggl(\frac{1}{q^3}-1\biggr) ~~~~\Rightarrow ~~~~ \nu = q^3 \, . </math>

The gravitational potential energy is,

<math>W\biggr|_{\rho_e/\rho_c = 1}</math>

<math> = - \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) 3\nu^2 \xi_s \biggl\{ \frac{1}{5} + \frac{1}{2} (\xi_s^2 - 1) + \biggl[ \frac{1}{5}(\xi_s^5 - 1) - \frac{1}{2}(\xi_s^2-1) \biggr] \biggr\} </math>

 

<math> = - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr)\nu^2 \xi_s \biggl\{ \frac{1}{5} + \biggl[ \frac{1}{5}(\xi_s^5 - 1) \biggr] \biggr\} </math>

 

<math> = - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr)\nu^2 \xi_s^6 = - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) \biggl( \frac{\nu}{q^3}\biggr)^2 = - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) \, . </math>

Drawing on expressions developed in our introductory discussion of the virial equation, the internal energy of the bipolytropic configuration is,

<math> U = U_\mathrm{core} + U_\mathrm{env} </math>

<math>=</math>

<math> \biggl\{ M_\mathrm{core} \biggl[ (1 - \delta_{\infty n_c}) n_c K_c \rho_c^{1/n_c} + \delta_{\infty n_c} c_s^2 \ln(\rho_c/\rho_0) \biggr] + n_e M_\mathrm{env} K_e \rho_e^{1/n_e} \biggr\} \, , </math>

where <math>\rho_0</math> is an, as yet unspecified, normalization density, and we have allowed for either an isothermal (<math>\delta_{\infty n_c} = 1</math>) or an adiabatic (<math>\delta_{\infty n_c} = 0</math>) core.


Strategy

Scaling Parameters

We want to vary the total radius, <math>R \,</math>, of the configuration and look for extrema in the free energy, while holding the following parameters fixed: <math>M_\mathrm{tot} \,</math>, <math>c_s^2 \,</math> (or <math>K_c \,</math>), <math>\nu \,</math>, and <math>q \,</math>. So we need to rewrite the expressions for <math>W \,</math> and <math>U \,</math> such that everything is constant except for <math>R \,</math>, or <math>\chi \equiv R/R_0 \,</math>. And, in order to put everything explicitly in terms of the fixed parameters just specified, let's go ahead and define the length, time, and density against which all dimensional quantities will ultimately be scaled.

Chosen Scaling Parameters

 

Polytropic Core

Isothermal Core

<math>R_0 \,</math>

<math>\biggl[ \biggl( \frac{K_c}{G} \biggr)^{n_c} M_\mathrm{tot}^{1-n_c} \biggr]^{1/(3-n_c)} \,</math>

<math>\frac{GM_\mathrm{tot}}{c_s^2} \,</math>

<math>t_0 \,</math>

<math>\biggl[ \biggl( \frac{K_c^3}{M_\mathrm{tot}^2} \biggr)^{n_c} G^{-(2n_c+3)} \biggr]^{1/(6-2n_c)} \,</math>

<math>\frac{GM_\mathrm{tot}}{c_s^3} \,</math>

<math>\rho_0 \equiv \frac{3M_\mathrm{tot}}{4\pi R_0^3} \,</math>

<math>\frac{3}{4\pi} \biggl[ \frac{G^3 M_\mathrm{tot}^2}{K_c^3} \biggr]^{n_c/(3-n_c)} \,</math>

<math>\frac{3}{4\pi} \biggl[ \frac{c_s^6}{G^3 M_\mathrm{tot}^2} \biggr] \,</math>


Rescaled Energy Expressions

Isothermal Core

In the case of an isothermal core, the expression for <math>W</math> shown above needs only a slight modification to put it in the appropriate form, namely,

<math>W \,</math>

<math> = - \frac{3}{5} M_\mathrm{tot} c_s^2 \chi^{-1} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, , </math>

In the case of an isothermal core, the expression for the total internal energy may be rewritten as,

<math> \frac{U}{M_\mathrm{tot}} \, </math>

<math>= \,</math>

<math> \nu c_s^2 \ln(\rho_c/\rho_0) + (1-\nu) n_e K_e \rho_e^{1/n_e} </math>

 

<math>=</math>

<math> \nu c_s^2 [\ln(\rho_c/\bar{\rho}) + \ln(\bar\rho/\rho_0)] + (1-\nu) n_e K_e \rho_0^{1/n_e} \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr) \biggr]^{1/n_e} </math>

 

<math>=</math>

<math> \nu c_s^2 [\ln(\nu/q^3) -3 \ln\chi] + (1-\nu) n_e K_e \rho_0^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \chi^{-3/n_e} \, . </math>

Hence (in the case of an isothermal core),

<math>\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{M_\mathrm{tot} c_s^2}</math>

<math>=</math>

<math>

- \frac{3}{5}  \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \chi^{-1} 

+ \nu [\ln(\nu/q^3) -3 \ln\chi] + (1-\nu) n_e \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \chi^{-3/n_e} </math>

 

<math>=</math>

<math> \mathfrak{G}_0^* - A \chi^{-1} - C \ln\chi + D \chi^{-3/n_e} \, , </math>

where,

<math>\mathfrak{G}_0^* \,</math>

<math>\equiv \,</math>

<math> \nu \ln(\nu/q^3) \, , </math>

<math>A \,</math>

<math>\equiv \,</math>

<math> \frac{3}{5} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, , </math>

<math>C \,</math>

<math>\equiv \,</math>

<math> 3 \nu \, , </math>

<math>D \, </math>

<math>\equiv \,</math>

<math> (1-\nu) n_e \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} = (1-\nu) n_e \biggl( \frac{K_e}{c_s^2} \biggr) \biggl[ \frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr) \biggr]^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \, . </math>

Polytropic Core

In the case of a polytropic core,

<math>\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{M_\mathrm{tot} K_c \rho_0^{1/n_c}}</math>

<math>=</math>

<math> - A \chi^{-1} + B_c \chi^{-3/n_c} + B_e \chi^{-3/n_e} \, , </math>

where,

<math>A \,</math>

<math>\equiv \,</math>

<math> \frac{3}{5} \biggl[ \frac{GM_\mathrm{tot}}{R_0 K_c \rho_0^{1/n_c}} \biggr] \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) = \frac{3}{5} \biggl( \frac{4\pi}{3} \biggr)^{1/n_c} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, , </math>

<math>B_c \,</math>

<math>\equiv \,</math>

<math> n_c \nu \biggl( \frac{\nu}{q^3}\biggr)^{1/n_c} \, , </math>

<math>B_e \, </math>

<math>\equiv \,</math>

<math> n_e (1-\nu) \biggl[ \biggl( \frac{K_e }{K_c}\biggr)\rho_0^{1/n_e - 1/n_c} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \, . </math>

Temperature Across the Interface

In order to ensure that the temperature of the envelope is the same as the temperature of the core when there is a drop in the mean molecular weight at the interface, we need to have,

<math>\biggl[ \frac{P_e}{(\rho_e/\mu_e)}\biggr]_E</math>

<math>= \, </math>

<math>\biggl[ \frac{P_c}{(\rho_c/\mu_c)} \biggr]_E \, .</math>

Note that this reflects the same physical condition as the constraint that is placed on the enthalpy at the interface when we analyze the detailed structure of <math>n_c = 5</math>, <math>n_e=1</math> bipolytropes (see the middle of Table 1 in the accompanying discussion), namely,

<math>\biggl[ \frac{\mu_e H_e}{(n_e+1)}\biggr]_i</math>

<math>= \, </math>

<math>\biggl[ \frac{\mu_c H_c}{(n_c+1)} \biggr]_i \, .</math>


In the case of an isothermal core, this implies,

<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ c_s^2}</math>

<math>= \, </math>

<math>\biggl[\rho_e^{-1/n_e}\biggr]_E</math>

<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e \rho_0^{1/n_e}}{ c_s^2}</math>

<math>= \, </math>

<math> \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>

 

<math>= \, </math>

<math> \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e} = \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \chi_E^{3/n_e} </math>

Not Necessarily Useful

<math>\Rightarrow ~~~~~ \chi_E^{3/n_e}</math>

<math>= \, </math>

<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ c_s^2} \biggl[ \frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr) \biggr]^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e}</math>

 

<math>= \, </math>

<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl[ \frac{D}{(1-\nu)n_e M_\mathrm{tot} c_s^2} \biggr]</math>

 

<math>= \, </math>

<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl[ \frac{3\nu}{(1-\nu)n_e} \biggl(\frac{D}{C}\biggr)\biggr]</math>

Therefore, the coefficient <math>D\,</math> becomes,

<math>D \, </math>

<math>\equiv \,</math>

<math> n_e (1-\nu) \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} = n_e (1-\nu) \biggl( \frac{\mu_c }{\mu_e}\biggr) \chi_E^{3/n_e} \, . </math>


In the case of a polytropic core, this implies,

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

<math>= \, </math>

<math>\rho_0^{1/n_c} \biggl[\rho_e^{-1/n_e}\biggr]_E</math>

<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl( \frac{K_e }{ K_c} \biggr) \rho_0^{(1/n_e-1/n_c)}</math>

<math>= \, </math>

<math> \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>

 

<math>= \, </math>

<math> \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e} = \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \chi_E^{3/n_e} </math>

Therefore, the coefficient <math>B_e\,</math> becomes,

<math>B_e \, </math>

<math>\equiv \,</math>

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


Pressure Across the Interface

We will relate <math>K_e</math> to <math>K_c</math> by demanding that initially the pressure is identical in both layers. The relevant algebraic relation will depend on whether the core is isothermal (<math>n_c = \infty</math>), or whether it has a finite polytropic index and therefore adjusts adiabatically to compressions or expansions.

Isothermal Core

Guided by the interface conditions presented in Table 2 of our accompanying discussion of the structure of bipolytropes, the condition for pressure balance in the case of an isothermal core should be,

<math>\frac{c_s^2}{K_e \rho_0^{1/n_e}} = \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\rho_e}{\rho_c} \biggr)^{1+1/n_e} = \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, .</math>

Actually, in the full structural solution,

<math>\frac{c_s^2}{K_e \rho_0^{1/n_e}}</math>

<math>= </math>

<math> \biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e} e^{\psi_i} \biggl( \frac{\rho_e}{\rho_c} \biggr)^{1+1/n_e} \phi_i^{1+n_e} </math>

 

<math>= </math>

<math> \biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e} e^{\psi_i} \biggl[ \frac{\mu_e}{\mu_c} e^{-\psi_i} \phi_i^{-{n_e}}\biggr]^{1+1/n_e} \phi_i^{1+n_e} </math>

 

<math>= </math>

<math> \biggl( \frac{\rho_c}{\rho_0} e^{-\psi_i} \biggr)^{1/n_e} \biggl[ \frac{\mu_e}{\mu_c} \biggr]^{1+1/n_e} \, . </math>

Adiabatic Core

Guided by the interface conditions presented in Table 2 of our accompanying discussion of the structure of bipolytropes, the condition for pressure balance in the case of polytropic core should be,

<math>\biggl(\frac{K_c}{K_e}\biggr) \rho_0^{1/n_c -1/n_e} = \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\rho_e}{\rho_c} \biggr)^{1+1/n_e} = \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, .</math>

Actually, in the full structural solution,

<math>\biggl(\frac{K_c}{K_e}\biggr) \rho_0^{1/n_c -1/n_e} </math>

<math>=</math>

<math> \rho_0^{1/n_c -1/n_e} \rho_c^{-(1+1/n_c)} \theta_i^{-(1+n_c)} \rho_e^{1+1/n_e} \phi_i^{1+n_e} </math>

 

<math>=</math>

<math> \biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e - 1/n_c} \theta_i^{-(1+n_c)} \biggl( \frac{\rho_e}{\rho_c}\biggr)^{1+1/n_e} \phi_i^{1+n_e} </math>

 

<math>=</math>

<math> \biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e - 1/n_c} \biggl[ \frac{\mu_e}{\mu_c} \theta_i^{n_c}

\biggr]^{1+1/n_e} \theta_i^{-(1+n_c)}

</math>

 

<math>=</math>

<math> \biggl( \frac{\rho_c}{\rho_0} \theta_i^{n_c}\biggr)^{1/n_e - 1/n_c} \biggl[ \frac{\mu_e}{\mu_c}

\biggr]^{1+1/n_e} \, .

</math>

Virial Analysis

Isothermal Core Equilibrium Condition

For a given set of fixed coefficients in the free energy expression, equilibria are identified by setting <math>\partial\mathfrak{G}/\partial\chi = 0</math>. Generally,

<math>\frac{\partial\mathfrak{G}}{\partial\chi}</math>

<math>= \,</math>

<math> A\chi^{-2} - C\chi^{-1} - \biggl(\frac{3}{n_e}\biggr)D\chi^{-(1+3/n_e)} \, . </math>

So, for the case of <math>n_e = 3/2</math>, the equilibrium radius is given by the condition,

<math>A\chi_E^{-2}</math>

<math>= \,</math>

<math>

C\chi_E^{-1} + 2D\chi_E^{-3}

</math>

<math> \Rightarrow ~~~~ C\chi_E^{2} - A\chi_E + 2D </math>

<math>= \,</math>

<math>0 \,</math>

<math> \Rightarrow ~~~~ \chi_E </math>

<math>= \,</math>

<math> \frac{1}{2C} \biggl[ A \pm \biggl( A^2 - 8DC \biggr)^{1/2}\biggr] </math>

<math> \Rightarrow ~~~~ \chi_E </math>

<math>= \,</math>

<math> \frac{A}{2C} \biggl[ 1 \pm \biggl( 1 - \frac{8DC}{A^2} \biggr)^{1/2}\biggr] </math>

 

<math>= \,</math>

<math> \chi_0 [ 1 \pm \sqrt{1-\Pi} ] \, , </math>

where,

<math>\chi_0 \,</math>

<math>\equiv \,</math>

<math> \frac{A}{2C} \, , </math>

    and    

<math>\Pi \,</math>

<math>\equiv \,</math>

<math> \frac{8DC}{A^2} \, . </math>

Combined Equilibrium Constraints

But, from the condition on the temperature at the interface we also need,

<math>\chi_E \,</math>

<math>=\,</math>

<math> \biggl(\frac{\mu_e}{\mu_c} \biggr)^{1/2} \biggl[ \frac{2\nu}{(1-\nu)} \biggl( \frac{D}{C} \biggr)\biggr]^{1/2} </math>

 

<math>=\,</math>

<math> \chi_0 \biggl(\frac{\mu_e}{\mu_c} \biggr)^{1/2} \biggl[ \biggl(\frac{\nu}{1-\nu}\biggr) \Pi \biggr]^{1/2} \, . </math>

Hence, the two conditions combined imply,

<math>1 \pm \sqrt{1-\Pi}</math>

<math>=\,</math>

<math> \lambda^{1/2} \Pi^{1/2} \, , </math>

where,

<math>\lambda \,</math>

<math>\equiv\,</math>

<math> \biggl(\frac{\mu_e}{\mu_c} \biggr) \biggl(\frac{\nu}{1-\nu}\biggr) \, . </math>

This, in turn, implies,

<math>1-\Pi \, </math>

<math>=\,</math>

<math> (\lambda^{1/2} \Pi^{1/2} -1 )^2\, </math>

 

<math>=\,</math>

<math> \lambda\Pi -2(\lambda\Pi)^{1/2} + 1 \, </math>

<math>\Rightarrow ~~~~2(\lambda\Pi)^{1/2} \,</math>

<math>=\,</math>

<math> \Pi(1+\lambda) \, </math>

<math>\Rightarrow ~~~~\Pi \,</math>

<math>=\,</math>

<math> \frac{4\lambda}{(1+\lambda)^2} \, . </math>

But notice, as well, that in order to have real roots of the equilibrium condition, we need <math>\Pi \le 1</math>. This means,

<math>4\lambda \,</math>

<math>\le \,</math>

<math> (1+\lambda)^2 \, </math>

<math>\Rightarrow ~~~~ \lambda^2 - 2\lambda + 1 \,</math>

<math>\ge \,</math>

<math> 0 \, . </math>

<math>\Rightarrow ~~~~ (\lambda - 1)^2 \,</math>

<math>\ge \,</math>

<math> 0 \, </math>

<math>\Rightarrow ~~~~ \lambda \,</math>

<math>\ge \,</math>

<math> 1 \, </math>

<math>\Rightarrow ~~~~ \frac{\mu_e}{\mu_c} \,</math>

<math>\ge \,</math>

<math> \frac{1}{\nu} - 1 \, </math>

<math>\Rightarrow ~~~~ \nu \,</math>

<math>\ge \,</math>

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

Free Energy Expression

To within an additive constant, the free energy may now be written as,

<math> \mathfrak{G} = W + U = - A \chi^{-1} + B_e \chi^{-3/n_e} + (1-\delta_{\infty n_c}) B_c \chi^{-3/n_c} - \delta_{\infty n_c} \biggl[ B_I \ln\chi - \frac{1}{3}B_I \ln \biggl( \frac{\rho_c|_0}{\rho_0} \biggr) \biggr] \, , </math>

where, <math>\chi \equiv R/R_0</math> and,

<math> A \, </math>

<math>=\,</math>

<math> \biggl( \frac{3GM^2_\mathrm{tot}}{5R_0} \biggr) \nu^2 \xi_s \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \biggr\} = \biggl( \frac{3GM^2_\mathrm{tot}}{5R_0} \biggr) \frac{\nu^2}{q} f(\nu, q) \, , </math>

<math>B_e \,</math>

<math>=\,</math>

<math> M_\mathrm{tot}(1-\nu) n_eK_e (\rho_e|_0)^{1/n_e} = M_\mathrm{tot} (1-\nu) K_e\rho_0^{1/n_e} n_e \biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{1/n_e} \, , </math>

<math>B_c \,</math>

<math>= \,</math>

<math> M_\mathrm{tot} \nu n_c K_c (\rho_c|_0)^{1/n_c} = M_\mathrm{tot} K_e\rho_0^{1/n_e} \nu n_c \biggl[\biggl(\frac{K_c}{K_e}\biggr) \rho_0^{1/n_c - 1/n_e} \biggr] \biggl( \frac{\rho_c|_0}{\rho_0}\biggr)^{1/n_c} = M_\mathrm{tot} K_e\rho_0^{1/n_e} \nu n_c \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \biggr] \biggl( \frac{\rho_c|_0}{\rho_0}\biggr)^{1/n_c} \, , </math>

<math>B_I\,</math>

<math>=\,</math>

<math> 3 M_\mathrm{tot} c_s^2 \nu = M_\mathrm{tot} K_e\rho_0^{1/n_e} \biggl[ \frac{c_s^2}{K_e\rho_0^{1/n_e} } \biggr] 3\nu = M_\mathrm{tot} K_e\rho_0^{1/n_e} \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e}\biggr] 3\nu \, . </math>

Derivatives of Free Energy

<math> \frac{\partial\mathfrak{G}}{\partial \chi} = A \chi^{-2} -\frac{3}{n_e} B_e \chi^{-(1+3/n_e)} - (1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi^{-(1+3/n_c)} - \delta_{\infty n_c} B_I \chi^{-1} \, ; </math>

<math> \frac{\partial^2\mathfrak{G}}{\partial \chi^2} = -2 A \chi^{-3} + \frac{3}{n_e} \biggl(1+\frac{3}{n_e}\biggr) B_e \chi^{-(2+3/n_e)} + (1-\delta_{\infty n_c}) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi^{-(2+3/n_c)} + \delta_{\infty n_c} B_I \chi^{-2} \, . </math>

Equilibrium Condition

We obtain the equilibrium radius, <math>\chi_E</math>, when <math>\partial\mathfrak{G}/\partial\chi = 0</math>. Hence, the relation governing the equilibrium radius is,

<math> A \chi_E^{-2} </math>

<math>=</math>

<math> (1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi_E^{-1- 3/n_c} +\delta_{\infty n_c} B_I \chi_E^{-1} +\frac{3}{n_e} B_e \chi_E^{-1-3/n_e} </math>

<math> \Rightarrow ~~~~~ \frac{n_e A}{3B_e} </math>

<math>=</math>

<math> (1-\delta_{\infty n_c}) \frac{n_e B_c}{n_c B_e} \chi_E^{1- 3/n_c} +\delta_{\infty n_c} \frac{n_e B_I}{3B_e} \chi_E + \chi_E^{1-3/n_e} </math>

<math> \Rightarrow ~~~~~ \chi_E^{1-3/n_e} </math>

<math>=</math>

<math> \alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \, , </math>

where,

<math>\alpha \,</math>

<math>\equiv \,</math>

<math> \frac{n_e A}{3B_e} = \biggl[ \frac{GM_\mathrm{tot}}{5R_0 K_e (\rho_e|_0)^{1/n_e} } \biggr] \frac{\nu^2 f(\nu, q)}{q} = \biggl[ \frac{GM_\mathrm{tot}}{5R_0 K_e \rho_0^{1/n_e} } \biggr]\biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{-1/n_e} \frac{\nu^2 f(\nu, q)}{q}

\, ;

</math>

<math>\beta \,</math>

<math>\equiv \,</math>

<math>\frac{n_e B_c}{n_c B_e} = \nu \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\rho_c|_0}{\rho_0}\biggr)^{1/n_c} \biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{-1/n_e} \biggr] \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, ; </math>

<math>\beta_I \,</math>

<math>\equiv \,</math>

<math> \frac{n_e B_I}{3B_e} = \frac{c_s^2 \nu}{(1-\nu)K_e (\rho_e|_0)^{1/n_e}} = \nu \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{-1/n_e} \biggr] \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, . </math>

Isothermal Core

In the case of an isothermal core (<math>\delta_{\infty n_c} = 1</math>),

<math> A \chi_E^{-2} </math>

<math>=</math>

<math> B_I \chi_E^{-1} +\frac{3}{n_e} B_e \chi_E^{-1-3/n_e} </math>

<math> \Rightarrow ~~~~ B_I \chi_E^{3/n_e} - A \chi_E^{3/n_e-1} +\frac{3}{n_e} B_e </math>

<math>=\,</math>

<math> 0 \, . </math>

Or, alternatively, via a dimensionless treatment,

<math> \chi_E^{1-3/n_e} </math>

<math>=</math>

<math> \alpha - \beta_I \chi_E \, . </math>

Adiabatic Core

In the case of an adiabatice core (<math>\delta_{\infty n_c} = 0</math>),

<math> \chi_E^{1-3/n_e} </math>

<math>=</math>

<math> \alpha - \beta \chi_E^{1- 3/n_c} \, . </math>

Stability

At this equilibrium radius, the second derivative of the free energy has the value,

<math> \chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E </math>

<math> = </math>

<math> -2 A\biggl( \frac{n_e}{3B_e} \biggr) + (1-\delta_{\infty n_c}) \biggl( \frac{n_e}{3B_e} \biggr) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi_E^{1-3/n_c} + \delta_{\infty n_c} \biggl( \frac{n_e}{3B_e} \biggr) B_I \chi_E + \frac{3}{n_e}\biggl( \frac{n_e}{3B_e} \biggr) \biggl(1+\frac{3}{n_e}\biggr) B_e \chi_E^{1-3/n_e} </math>

 

<math> = </math>

<math> -2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} + \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \chi_E^{1-3/n_e} \, , </math>

which, when combined with the condition for equilibrium gives,

<math> \chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E </math>

<math> = </math>

<math> -2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} + \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \biggl[ \alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \biggr] </math>

 

<math> = </math>

<math> \alpha \biggl(\frac{3}{n_e}-1\biggr) + (1-\delta_{\infty n_c}) \beta \biggl(\frac{3}{n_c}-\frac{3}{n_e}\biggr) \chi_E^{1-3/n_c} - \delta_{\infty n_c} \beta_I \biggl(\frac{3}{n_e}\biggr)\chi_E </math>

<math> \Rightarrow ~~~~ \chi_E^3 \biggl( \frac{n_e^2}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E </math>

<math> = </math>

<math> \alpha (3-n_e) - (1-\delta_{\infty n_c}) 3\beta \biggl(1 - \frac{n_e}{n_c} \biggr) \chi_E^{1-3/n_c} - \delta_{\infty n_c} 3\beta_I \chi_E \, . </math>

The equilibrium configuration is stable as long as this second derivative is positive.

Isothermal Core

Hence, for a bipolytrope with an isothermal core (<math>\delta_{\infty n_c} = 1</math>), the configuration is stable as long as,

<math> \chi_E < \frac{\alpha (3-n_e)}{3\beta_I} \, . </math>

Adiabatic Core

In the adiabatic case (<math>\delta_{\infty n_c} = 0</math>), the configuration is stable as long as,

<math> \chi_E^{1-3/n_c} < \frac{\alpha n_c (3-n_e)}{3\beta (n_c-n_e)} \, . </math>


Examples

Isothermal Core with <math>n_e=3/2</math>

Consider the case examined by Schönberg & Chandrasekhar (1942), that is, the case of an isothermal core and an envelope with <math>n_e = 3/2</math>. The equilibrium radius is given by the expression,

<math> \chi_E^{-1} </math>

<math>=</math>

<math> \alpha - \beta_I \chi_E \, </math>

<math> \Rightarrow ~~~~~ \beta_I \chi_E^{2} -\alpha \chi_E + 1 </math>

<math>=</math>

<math> 0 \, </math>

<math> \Rightarrow ~~~~~ \chi_E </math>

<math>=</math>

<math> \frac{1}{2\beta_I} \biggl[\alpha \pm \biggr(\alpha^2-4\beta_I \biggr)^{1/2}\biggr] \, </math>

 

<math>=</math>

<math> \frac{\alpha}{2\beta_I} \biggl[1 \pm \biggr(1-\frac{4\beta_I}{\alpha^2} \biggr)^{1/2}\biggr] \, . </math>

And the system is stable when,

<math> \chi_E < \chi_0 \equiv \frac{\alpha }{2\beta_I} \, . </math>

A couple of physical attributes are now clear:

  • Physical configurations only exist for <math>(4\beta_I/\alpha^2) \le 1</math>.
  • For each value of <math>(4\beta_I/\alpha^2) < 1 \,</math>, there are two equilibrium configurations, given by the <math>\pm</math> roots of the quadratic equation for <math>\chi_E</math>; the "negative" branch is stable but the "positive" branch is unstable.

Note that,

<math> \chi_0 \equiv \frac{\alpha }{2\beta_I} = \biggl( \frac{GM_\mathrm{tot}}{10R_0 c_s^2} \biggr) \frac{\nu f(\nu,q)}{q} ~~~\Rightarrow ~~~ \biggl( \frac{10R_0 c_s^2}{GM_\mathrm{tot}} \biggr) = \frac{\nu f(\nu,q)}{q \chi_0} \, , </math>

and,

<math> \frac{4\beta_I}{\alpha^2 } = \biggl( \frac{10 R_0 c_s^2 }{GM_\mathrm{tot}} \biggr)^2 \biggl[ \frac{K_e (\rho_e|_0)^{1/n_e}}{c_s^2} \biggr] \frac{q^2(1-\nu)}{\nu^3 f^2(\nu,q)} = \biggl[ \frac{K_e (\rho_e|_0)^{1/n_e}}{c_s^2} \biggr] \frac{q^2(1-\nu)}{\nu^3 f^2(\nu,q)} \biggl[\frac{\nu f(\nu,q)}{q \chi_0} \biggr]^2 = \biggl[ \frac{\mu_c}{\mu_e} \biggr] \biggl( \frac{1}{\nu} - 1 \biggr) \, . </math>

But, this last expression must be less than or equal to unity, which implies,

<math> \frac{1}{\nu} \le 1 + \frac{\mu_e}{\mu_c} ~~~\Rightarrow ~~~ \nu \ge \biggl(1 + \frac{\mu_e}{\mu_c} \biggr)^{-1} </math> This doesn't seem to have the correct behavior, for example, the smaller values of <math>\nu</math> should be the stable ones, so there must be a mistake in the derivation.

Adiabatic Core with <math>n_c = 5</math> and <math>n_e=1</math>

Consider the case with an analytical structure derived by Eagleton, Faulkner, and Cannon (1998, MNRAS, 298, 831), that is, the case of an adiabatic core having <math>n_c=5</math> and an envelope with <math>n_e = 1</math>. The equilibrium radius is,

Old and Probably Irrelevant Discussion

Summary Expressions (New)

In the above derivations, we have adopted the notation,

<math> \rho_\mathrm{norm} \equiv \frac{3M_\mathrm{tot}}{4\pi R_0^3} \, . </math>

Now, guided by the earlier discussion of pressure-bounded isothermal spheres, we choose the following normalization energy and radius:

<math> E_0 = 3M_\mathrm{tot} c_s^2 </math>         and         <math> R_0 = \frac{GM_\mathrm{tot}}{5c_s^2} \, . </math>

Also, by analogy, it is useful to define the dimensionless parameter,

<math> \Pi_I \equiv \frac{K_e \rho_\mathrm{norm}^{1/n_e}}{c_s^2} = \frac{K_e}{c_s^2} \biggl[ \frac{3M_\mathrm{tot}}{4\pi R_0^3} \biggr]^{1/n_e} = \biggl( \frac{3\cdot 5^3}{4\pi} \biggr)^{1/n_e} \frac{K_e}{c_s^2} \biggl[ \frac{c_s^6}{G^3 M_\mathrm{tot}^2} \biggr]^{1/n_e} \, . </math>

(It is worth noting that if we set <math>n_e = -1</math>, the dimensionless parameter <math>\Pi_I</math> becomes identical to the parameter <math>\Pi</math> as defined in the context of our discussion of the Bonnor-Ebert sphere. But in order to complete the analogy with the Bonnor-Ebert sphere discussion, we would also need to change the sign on the last term in the above expression for the free energy because in the earlier discussion the external pressure was an external, confining condition whereas here it is included as an internal energy of the system.)


Relevant Expressions for Isothermal Core

<math> \frac{\rho_e}{\rho_c}= \frac{\mu_e}{\mu_c} </math>

<math> \frac{q^3}{\nu}\biggl( \frac{1 - \nu}{1- q^3} \biggr) </math>

 

<math> \chi \equiv \frac{R}{R_0} </math>

<math> \Pi_I^{n_e/3} [\nu^{-n_e}\ (1-\nu)^{n_e+1} ]^{1/3} q^{n_e} (1-q^3)^{-(n_e+1)/3} </math>

 

<math> \chi^3 = \biggl(\frac{R}{R_0} \biggr)^{3} </math>

<math> \Pi_I^{n_e} \biggl( \frac{\rho_e}{\rho_c}\biggr)^{n_e + 1} \biggl[ q^3 + \biggl( \frac{\rho_e}{\rho_c} \biggr)(1-q^3) \biggr]^{-1} </math>

 

<math> \frac{A}{E_0} </math>

<math> \nu^2 \xi_s \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \biggr\} </math>

<math>= \biggl( \frac{\nu^2}{q} \biggr) \biggl\{ f_A(\nu,q) \biggr\}</math>

<math> \frac{B_e}{E_0} </math>

<math> \Pi_I \biggl( \frac{n_e}{3} \biggr) (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e} (\xi_s^3 - 1)^{-1/n_e} </math>

 

<math> \frac{B_I}{E_0} </math>

<math>\nu </math>

 

<math> \frac{\mathfrak{G}}{E_0} </math>

<math> - \frac{A}{E_0} \chi^{-1} - \frac{B_I}{E_0} \ln\chi + \frac{B_e}{E_0} \chi^{-3/n_e} </math>

<math> = \nu \biggl[ - \biggl( \frac{\nu}{q} \biggr) \frac{f_A(\nu,q)}{\chi} - \ln\chi + \frac{n_e}{3} \biggl(\frac{1}{q^3}-1\biggr)\biggr] </math>

[On 8 November 2013, J. E. Tohline wrote: I just confirmed that the simpler expression for the normalized total free energy, <math>\mathfrak{G}/E_0</math>, matches the more complicated version. I don't like the result because the third term in the free energy -- the one contributed by the internal energy of the envelope -- is independent of the radius of the configuration, <math>\chi</math>; it works out this way because my expression for <math>\Pi_I</math> has a radial dependence that exactly cancels out the explicit radial dependence that appears in the more complicated expression. But maybe it's okay after all because this expression is intended to show how the free energy varies across the <math>(q,\nu)</math> plane, and the effect of <math>\Pi_I</math> appears implicitly through the specification of <math>\chi</math>, or visa versa.]


Summary Expressions (Old)

In the above derivations, we have adopted the notation,

<math> \rho_\mathrm{norm} \equiv \frac{3M_\mathrm{tot}}{4\pi R_0^3} \, . </math>

Now, guided by the dimensional aspects of the various coefficients in the free energy expression, we choose the following normalization energy and radius:

<math> E_0 = M_\mathrm{tot} K_e \rho_\mathrm{norm}^{1/n_e} </math>         and         <math> R_0 = \frac{GM_\mathrm{tot}}{K_e \rho_\mathrm{norm}^{1/n_e}} \, . </math>

When combined with the expression for <math>\rho_\mathrm{norm}</math>, these become,

<math> E_0 = \biggl[ \biggl( \frac{4\pi}{3} \biggr) G^3 M_\mathrm{tot}^{5-n_e} K_e^{-n_e}\biggr]^{1/(3-n_e)} </math>         and         <math> R_0 = \biggl[ \frac{3M_\mathrm{tot}}{4\pi} \biggl( \frac{K_e}{GM_\mathrm{tot}} \biggr)^{n_e} \biggr]^{1/(3-n_e)} \, . </math>

So, the primary scales are determined after specifying two parameters: <math>M_\mathrm{tot}</math> and <math>K_e</math>. We also obtain,

<math> \kappa_I \equiv \frac{c_s^2}{K_e \rho_\mathrm{norm}^{1/n_e}} = \frac{M_\mathrm{tot} c_s^2}{E_0} = c_s^2 \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{K_e^{n_e}}{G^3 M_\mathrm{tot}^2} \biggr]^{1/(3-n_e)} \, . </math>


Relevant Expressions for Isothermal Core

<math> \frac{\rho_e}{\rho_c} </math>

<math> \frac{q^3}{\nu}\biggl( \frac{1 - \nu}{1- q^3} \biggr) </math>

<math> \chi \equiv \frac{R}{R_0} </math>

<math> q^{n_e} (1-q^3)^{-(n_e+1)/3} \biggl[\nu^{-n_e}\ (1-\nu)^{n_e+1} \kappa_I^{-n_e} \biggr]^{1/3} </math>

<math> \frac{A}{E_0} </math>

<math> \biggl( \frac{3}{5} \biggr)\nu^2 \xi_s \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \biggr\} </math>

<math> \frac{B_e}{E_0} </math>

<math> n_e(1 - \nu)^{1+1/n_e} \xi_s^{3/n_e} (\xi_s^3 - 1)^{-1/n_e} </math>

<math> \frac{B_I}{E_0} </math>

<math> 3 \kappa_I \nu </math>

<math> \frac{\mathfrak{G}}{E_0} </math>

<math> - \frac{A}{E_0} \chi^{-1} - \frac{B_I}{E_0} \ln\chi + \frac{B_e}{E_0} \chi^{-3/n_e} </math>

Subsequently, we will also find it useful to have expressions for the following coefficient ratios:

<math>\frac{n_e A}{3 B_e}</math>

<math>=</math>

<math> \biggl( \frac{GM^2_\mathrm{tot}}{5R_0} \biggr) \nu^2 \xi_s \biggl\{ f(\rho_e/\rho_c , \xi_s) \biggr\} \biggl[K_e M_\mathrm{tot} \rho_\mathrm{norm}^{1/n_e} (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e} (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1} </math>

 

<math>=</math>

<math> \biggl( \frac{GM_\mathrm{tot}}{5K_e \rho_\mathrm{norm}^{1/n_e} R_0} \biggr) \nu^2 \xi_s \biggl\{ f(\rho_e/\rho_c , \xi_s) \biggr\} (1 - \nu)^{-(1+1/n_e)} q^{3/n_e} \biggl(\frac{1}{q^3} - 1 \biggr)^{1/n_e} </math>

 

<math>=</math>

<math> \biggl( \frac{GM_\mathrm{tot}}{5K_e \rho_\mathrm{norm}^{1/n_e} R_0} \biggr) \nu^2 q^{-1} [ (1 - \nu)^{-(n_e+1)} (1 - q^3 ) ]^{1/n_e} \biggl\{ f(\rho_e/\rho_c , \xi_s) \biggr\} \, ; </math>


<math>\frac{n_e B_I}{3 B_e}</math>

<math>=</math>

<math> M_\mathrm{tot} c_s^2 \nu \biggl[K_e M_\mathrm{tot} \rho_\mathrm{norm}^{1/n_e} (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e} (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1} </math>

 

<math>=</math>

<math> \kappa_I \nu (1 - \nu)^{-(1+1/n_e)} q^{3/n_e} \biggl(\frac{1}{q^3} - 1 \biggr)^{1/n_e} </math>

 

<math>=</math>

<math> \kappa_I \nu \biggl[ (1 - \nu)^{-(n_e+1)} (1-q^{3}) \biggr]^{1/n_e} \, ; </math>


<math>\frac{n_e B_c}{n_c B_e}</math>

<math>=</math>

<math> \biggl[K_c M_\mathrm{tot} \rho_\mathrm{norm}^{1/n_c} \nu^{1+1/n_c}\xi_s^{3/n_c} \biggr] \biggl[K_e M_\mathrm{tot} \rho_\mathrm{norm}^{1/n_e} (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e} (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1} </math>

 

<math>=</math>

<math> \kappa_A \nu^{1+1/n_c}\xi_s^{3/n_c} \biggl[ (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e} (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1} </math>

 

<math>=</math>

<math> \kappa_A \biggl[\nu^{1+1/n_c} (1 - \nu)^{-(1+1/n_e)} \biggr] \biggl[ \xi_s^{3(1/n_c-1/n_e)} (\xi_s^3 - 1)^{1/n_e}\biggr] </math>

 

<math>=</math>

<math> \kappa_A \biggl[\nu^{1+1/n_c} (1 - \nu)^{-(1+1/n_e)} \biggr] \biggl[ q^{3(1/n_e-1/n_c)} \biggl( \frac{1}{q^3} - 1\biggr)^{1/n_e}\biggr] </math>

 

<math>=</math>

<math> \kappa_A \biggl[\nu^{1+1/n_c} (1 - \nu)^{-(1+1/n_e)} \biggr] \biggl[ q^{-3n_e/n_c} (1-q^3) \biggr]^{1/n_e} </math>

 

<math>=</math>

<math> \kappa_A \biggl[\nu^{1+1/n_c} q^{-3/n_c} \biggr] \biggl[ (1 - \nu)^{-(n_e+1)} (1-q^3) \biggr]^{1/n_e} \, . </math>

Derivatives of Free Energy

<math> \frac{\partial\mathfrak{G}}{\partial \chi} = A \chi^{-2} -(1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi^{-(1+3/n_c)} - \delta_{\infty n_c} B_I \chi^{-1} -\frac{3}{n_e} B_e \chi^{-(1+3/n_e)} \, ; </math>

<math> \frac{\partial^2\mathfrak{G}}{\partial \chi^2} = -2 A \chi^{-3} + (1-\delta_{\infty n_c}) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi^{-(2+3/n_c)} + \delta_{\infty n_c} B_I \chi^{-2} + \frac{3}{n_e} \biggl(1+\frac{3}{n_e}\biggr) B_e \chi^{-(2+3/n_e)} \, . </math>

Equilibrium Condition

We obtain the equilibrium radius, <math>\chi_E</math>, when <math>\partial\mathfrak{G}/\partial\chi = 0</math>. Hence, the relation governing the equilibrium radius is,

<math> A \chi_E^{-2} </math>

<math>=</math>

<math> (1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi_E^{-1- 3/n_c)} +\delta_{\infty n_c} B_I \chi_E^{-1} +\frac{3}{n_e} B_e \chi_E^{-1-3/n_e)} </math>

<math> \Rightarrow ~~~~~ \frac{n_e A}{3B_e} </math>

<math>=</math>

<math> (1-\delta_{\infty n_c}) \frac{n_e B_c}{n_c B_e} \chi_E^{1- 3/n_c} +\delta_{\infty n_c} \frac{n_e B_I}{3B_e} \chi_E + \chi_E^{1-3/n_e} </math>

<math> \Rightarrow ~~~~~ \chi_E^{1-3/n_e} </math>

<math>=</math>

<math> \alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \, , </math>

where,

<math>\alpha \equiv \frac{n_e A}{3B_e} \, ; ~~~ \beta \equiv \frac{n_e B_c}{n_c B_e} \, ; ~~~ \beta_I \equiv \frac{n_e B_I}{3B_e} \, .</math>

Stability

At this equilibrium radius, the second derivative of the free energy has the value,

<math> \chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E </math>

<math> = </math>

<math> -2 A\biggl( \frac{n_e}{3B_e} \biggr) + (1-\delta_{\infty n_c}) \biggl( \frac{n_e}{3B_e} \biggr) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi_E^{1-3/n_c} + \delta_{\infty n_c} \biggl( \frac{n_e}{3B_e} \biggr) B_I \chi_E + \frac{3}{n_e}\biggl( \frac{n_e}{3B_e} \biggr) \biggl(1+\frac{3}{n_e}\biggr) B_e \chi_E^{1-3/n_e} </math>

 

<math> = </math>

<math> -2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} + \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \chi_E^{1-3/n_e} \, , </math>

which, when combined with the condition for equilibrium gives,

<math> \chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E </math>

<math> = </math>

<math> -2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} + \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \biggl[ \alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \biggr] </math>

 

<math> = </math>

<math> \alpha \biggl(\frac{3}{n_e}-1\biggr) + (1-\delta_{\infty n_c}) \beta \biggl(\frac{3}{n_c}-\frac{3}{n_e}\biggr) \chi_E^{1-3/n_c} - \delta_{\infty n_c} \beta_I \biggl(\frac{3}{n_e}\biggr)\chi_E </math>

<math> \Rightarrow ~~~~ \chi_E^3 \biggl( \frac{n_e^2}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E </math>

<math> = </math>

<math> \alpha (3-n_e) - (1-\delta_{\infty n_c}) 3\beta \biggl(1 - \frac{n_e}{n_c} \biggr) \chi_E^{1-3/n_c} - \delta_{\infty n_c} 3\beta_I \chi_E \, . </math>

Finally, the equilibrium configuration is stable as long as this second derivative is positive. Hence, for a bipolytrope with an isothermal core (<math>\delta_{\infty n_c} = 1</math>), the configuration is stable as long as,

<math> \chi_E < \frac{\alpha (3-n_e)}{3\beta_I} \, . </math>

In the adiabatic case (<math>\delta_{\infty n_c} = 0</math>), the configuration is stable as long as,

<math> \chi_E^{1-3/n_c} < \frac{\alpha n_c (3-n_e)}{3\beta (n_c-n_e)} \, . </math>


Examples

Isothermal Core with <math>n=3/2</math> Envelope

When the core is isothermal and <math>n_e = 3/2</math>, the equilibrium condition is:

<math> \chi_E^{-1} = \alpha - \beta_I \chi_E \, , </math>

<math> \Rightarrow ~~~~ \beta_I \chi_E^2 - \alpha \chi_E + 1 = 0 \, , </math>

<math> \Rightarrow ~~~~ \chi_E = \frac{1}{2\beta_I} \biggl[ \alpha \pm \sqrt{\alpha^2 - 4\beta_I} \biggr] = \frac{\alpha}{2\beta_I} \biggl[ 1 \pm \sqrt{1 - \frac{4\beta_I}{\alpha^2}} \biggr] \, . </math>

At the same time, the condition for stability is,

<math> \chi_E < \frac{\alpha}{2\beta_I} \, . </math>

Isothermal Core with <math>n=1</math> Envelope

When the core is isothermal and <math>n_e = 1</math>, the equilibrium condition is:

<math> \chi_E^{-2} = \alpha - \beta_I \chi_E \, , </math>

<math> \Rightarrow ~~~~ \chi_E^3 - \frac{\alpha}{\beta_I} \chi_E^2 + \frac{1}{\beta_I} = 0 \, . </math>

(We need to solve this cubic equation.)

At the same time, the condition for stability is,

<math> \chi_E < \frac{2\alpha}{3\beta_I} \, . </math>

Old (and probably incorrect) cases

Envelope with <math>n=3/2</math>

If we choose an <math>n_e = 3/2</math> envelope, we obtain stability for,

<math> \chi_E < \frac{2\beta}{\alpha}\, . </math>

In this case, the equilibrium radius condition is,

<math> \chi_E^2 - \alpha \chi_E + \beta =0 </math>

<math> \Rightarrow ~~~~ \chi_E = \frac{1}{2}\biggl[\alpha \pm \biggl( \alpha^2 -4\beta \biggr)^{1/2} \biggr] = \frac{\alpha}{2}\biggl[1 \pm \biggl( 1 -\frac{4\beta}{\alpha^2} \biggr)^{1/2} \biggr] </math>


Envelope with <math>n=1</math>

If, instead, we choose an <math>n_e = 1</math> envelope, we obtain stability for,

<math> \chi_E < \sqrt{\frac{3\beta}{\alpha} }\, . </math>

In this case, the equilibrium radius condition is,

<math> \alpha = \chi_E + \beta \chi_E^{-2} \, , </math>

<math> \Rightarrow ~~~~ \chi_E^3 - \alpha \chi_E^2 + \beta = 0 \, . </math>


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