User:Tohline/SSC/Virial/FormFactors
Structural Form Factors
|
|---|
| | Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
As has been defined in a companion, introductory discussion, three key dimensionless structural form factors are:
|
<math>~\mathfrak{f}_M </math> |
<math>~\equiv</math> |
<math>~ \int_0^1 3\biggl[ \frac{\rho(x)}{\rho_c}\biggr] x^2 dx \, ,</math> |
|
<math>~\mathfrak{f}_W</math> |
<math>~\equiv</math> |
<math>~ 3\cdot 5 \int_0^1 \biggl\{ \int_0^x \biggl[ \frac{\rho(x)}{\rho_c}\biggr] x^2 dx \biggr\} \biggl[ \frac{\rho(x)}{\rho_c}\biggr] x dx\, ,</math> |
|
<math>~\mathfrak{f}_A</math> |
<math>~\equiv</math> |
<math>~ \int_0^1 3\biggl[ \frac{P(x)}{P_c}\biggr] x^2 dx \, ,</math> |
where, <math>~x \equiv r/R_\mathrm{limit}</math>.
One Detailed Example (n = 5)
Here we derive detailed expressions for the above subset of structural form factors in the case of spherically symmetric configurations that obey an <math>~n=5</math> polytropic equation of state. The hope is that this will illustrate, in a clear and helpful manner, how the task of calculating form factors is to be carried out, in practice; and, in particular, to provide one nontrivial example for which analytic expressions are derivable. This should simplify the task of debugging numerical algorithms that are designed to calculate structural form factors for more general cases that cannot be derived analytically. The limits of integration will be specified in a general enough fashion that the resulting expressions can be applied, not only to the structures of isolated polytropes, but to pressure-truncated polytropes that are embedded in a hot, tenuous external medium and to the cores of bipolytropes.
Key Foundations
We use the following normalizations, as drawn from our more general introductory discussion:
|
Adopted Normalizations <math>~(n=5; ~\gamma=6/5)</math> | |||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
|||||||||||||||||||
|
Note that the following relations also hold: <math>~E_\mathrm{norm} = P_\mathrm{norm} R_\mathrm{norm}^3 = \frac{G M_\mathrm{tot}^2}{ R_\mathrm{norm}} = \biggl( \frac{3}{4\pi} \biggr) M_\mathrm{tot} c_\mathrm{norm}^2</math> | |||||||||||||||||||
As is detailed in our accompanying discussion of bipolytropes — see also our discussion of the properties of isolated polytropes — in terms of the dimensionless Lane-Emden coordinate, <math>~\xi \equiv r/a_{5}</math>, where,
<math> a_{5} =\biggr[ \frac{3K}{2\pi G} \biggr]^{1/2} \rho_0^{-2/5} \, , </math>
the radial profile of various physical variables is as follows:
|
<math>~\frac{r}{[K^{1/2}/(G^{1/2}\rho_0^{2/5})]}</math> |
<math>~=</math> |
<math>~\biggl( \frac{3}{2\pi} \biggr)^{1/2} \xi \, ,</math> |
|
<math>~\frac{\rho}{\rho_0}</math> |
<math>~=</math> |
<math>~\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2} \, ,</math> |
|
<math>~\frac{P}{K\rho_0^{6/5}}</math> |
<math>~=</math> |
<math>~\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3} \, ,</math> |
|
<math>~\frac{M_r}{[K^{3/2}/(G^{3/2}\rho_0^{1/5})]}</math> |
<math>~=</math> |
<math>~\biggl( \frac{2\cdot 3}{\pi } \biggr)^{1/2} \biggl[ \xi^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3/2} \biggr] \, .</math> |
Notice that, in these expressions, the central density, <math>~\rho_0</math>, has been used instead of <math>~M_\mathrm{tot}</math> to normalize the relevant physical variables. We can switch from one normalization to the other by realizing that — see, again, our accompanying discussion — in isolated <math>~n=5</math> polytropes, the total mass is given by the expression,
<math>M_\mathrm{tot} = \biggr[ \frac{2\cdot 3^4 K^3}{\pi G^3} \biggr]^{1/2} \rho_0^{-1/5} ~~~~\Rightarrow ~~~~ \rho_0^{1/5} = \biggr[ \frac{2\cdot 3^4 K^3}{\pi G^3} \biggr]^{1/2} M_\mathrm{tot}^{-1} \, .</math>
Employing this mapping to switch to our "preferred" adopted normalizations, as defined in the above boxed-in table, the four radial profiles become,
|
<math>~\frac{r}{R_\mathrm{norm}}</math> |
<math>~=</math> |
<math>~ \biggl( \frac{\pi}{2\cdot 3^4} \biggr) \biggl( \frac{3}{2\pi} \biggr)^{1/2} \xi = \biggl( \frac{\pi}{2^3\cdot 3^7} \biggr)^{1/2} \xi \, ,</math> |
|
<math>~\frac{\rho}{\rho_\mathrm{norm}}</math> |
<math>~=</math> |
<math>~\biggl( \frac{2^3\cdot 3^6}{\pi} \biggr)^{3/2} \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2} \, ,</math> |
|
<math>~\frac{P}{P_\mathrm{norm}}</math> |
<math>~=</math> |
<math>~\biggl( \frac{2\cdot 3^4}{\pi} \biggr)^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3} \, ,</math> |
|
<math>~\frac{M_r}{M_\mathrm{tot}}</math> |
<math>~=</math> |
<math>~ \biggl( \frac{\pi}{2\cdot 3^4} \biggr)^{1/2} \biggl( \frac{2\cdot 3}{\pi } \biggr)^{1/2} \biggl[ \xi^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3/2} \biggr] = \biggl[\frac{\xi^2}{3}\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-1}\biggr]^{3/2} \, .</math> |
|
|
|---|
|
© 2014 - 2021 by Joel E. Tohline |