Part 2: How I Implemented the Stellar Halos for all .Sci Galaxies
Well, the actual implementation of stellar halos for all 10000+ celestia.Sci galaxies requires the use of all
methods I have previously exposed in a tutorial fashion:
Following the literature, I have assumed for simplicity that the stellar halo properties measured explicitly for the Andromeda galaxy (M 31) are approximately universal throughout. Moreover (in a 3rd part), I use a supposedly universal stellar Luminosity function as extracted from nearby Milkyway stars (including HIPPARCOS stars) and as parametrized in E) above. Stellar Halo Properties (see A) above)
The galactic halo properties were taken from the FREERef 1
= S. Courteau et al., APJ, 739 20 (2011)
, The Luminosity Profile and Structural Parameters of the Andromeda Galaxy
and a second even more recent APJ paper, specializing on the surface brightness profile of the stellar M 31 halo,Ref 2
= Karoline M. Gilbert et al. 2012 ApJ 760 76 (2012)
, Global Properties of M31's Stellar Halo from the Splash Survey: I. Surface Brightness Profile
, unfortunately requires a costly subscription. However, here is the FREE arXive version
The shape of the skyplane projected halo is essentially circular
, corresponding to a spherical
shape in 3D. The authors find that the following inverse power law function fits the measured surface brightness of the halo as function of the projected radial distance R from the center:
R* = 30 kpc and a "core radius" parameter a_h = 5.2 kpc
. From the analysis, a surprisingly large halo radius R_halo >= 175 kpc
(2/3 of M 31's virial radius!) is deduced. The power index α was determined as α = 1.1 ± 0.1
which is very close to a drop off expected from pure gravity ~ 1/R^2, i.e. α = 1.0.
Here is a summary display of halo measurements and a display of the above fit from Ref 2 (arXiv, Fig. 10).
halo_data_M31.jpg [ 30.29 KiB | Viewed 2377 times ]
the data points are color-coded by the field's position angle, counted from M 31's minor axis. It is evident that the data is consistent with a circular shape
of the skyplane projected halo!3D Deprojection of the Measured Surface Brightness Profile (see B) above)
The next important step is to deproject
the observed circular skyplane projected surface brightness distribution into a proper distribution referring to 3D space
, by following the exposed mathematics in B) above.
The constant R* = 30 kpc above does not enter the present discussion, since we can equivalently rewrite the fit function somewhat simplified as
in terms of a new reference flux constant i0.
As a quick orientation about the geometries involved, let me recall my drawing from the general discussion B):
[click on image for a bigger size]
sphere3d.jpg [ 59.93 KiB | Viewed 2377 times ]
Keep in mind that R and r denote the projected and true halo radius, respectively, as apparent from the figure!
According to B), the radial behaviour of the halo's 3D surface brightness is obtained as solution of an Abel integral equation that takes a simple integral form in terms of the fitted 2D projected expression i(R):
The result reduces again to a simple inverse power law that decreases a bit faster, however, like 1/r^(2*alpha+1) for large r! Note the Γ function being a generalization of the familiar factorial n! = 1*2*3*....*n to non-integer argument.
Next, it is convenient to introduce 2 new parameters and a dimensionless, fractional
via the transformation
Obviously, η takes values in (0,1)
. The second dimensionless parameter ratio
expresses the relative
central core size in terms of a ratio of the core radius and the total halo radius r_halo.
The result now takes the simpler form
As it should be, the 3D surface brightness is seen to vanish for η = 1 or equivalently for r = r_halo. Random Halo Star Generation (see C, D) above)
While in 3D, the random generation of the isotropic
angular dependence involves of course just a uniform
(i.e. constant!) Probability Density Function (PDF
) (see C) above), the task to properly generate the star density as function of the distance r from center is a much less trivial affair (see D) above).
For the random halo star generation we conveniently use spherical polar coordinates
and express the differential volume
element dV in this familiar form
dV= r^2 dr dΏ = r_halo^3 * η^2 dη dΏ
For the random generation of the isotropic angular distribution, we again take a look at the differential solid angle element dΏ
as in C) above
dΏ = sinΦ dΦ dθ = - d(cosΦ
treating correspondingly θ and cosΦ (NOT Φ) as uniformly distributed random variables.
Next let us focus our attention on the radial dependence.
Note first of all, that our deprojected radial surface brightness distribution i_3D above has the meaning of a volume density (i.e. [#stars / volume]!). Thus, in order to normalize its integral to 1, as needed for a proper PDF, we now need to integrate over the spherical halo volume
, like so
Note the resulting factor of η^2 which will suppress the halo star density towards small η and the factor 4π from the trivial integration over the solid angle.
Hence the properly normalized halo star density that will serve as our PDF, takes the following form
Since the normalization integral over i_3D is rather long, I displayed it in unevaluated form. Of course (by construction),
Here is a plot of our PDF that may well be more instructive:
You clearly see that our PDF with its strong depletion towards η -> 0, its sharp peak and its vanishing for η ->1 is very far from a constant, uniform PDF. This is why need to use a combination of the appropriate two methods described in D).
The first step consists in finding a suitable majorant function
, simple enough in order to allow using the Inverse Transform Method
described in D). Remember:
The majorant function needs also to be normalized to a unit area (just like the exact PDF), but after multiplication with a suitable constant C has to enclose our exact PDF entirely! The closer the shape of the majorant is to the exact PDF, the higher the efficiency!
In our case, finding the best majorant with these properties is relatively easy.
Apart from a constant normalization factor, it looks like so:
where the majorant property is based on the obvious inequality:
and thus (for α > 0, w>0)
By comparing the exact PDF with this majorant, you now can prove the majorant property in 2 lines
. Indeed, here is a comparison plot which clearly illustrates the claim:
Look at the legend below the plot. The constant C = 1.214911633 > 1
and correspondingly we get a very high efficiency (100/C)
i.e. for our case
resulting in a sufficiently simple elementary function on the right!
Therefore, we may easily solve the inversion equation (see D) above!) U(η) = u
for η with the result
This equation based on the majorant PDF allows ingenously to use uniform random values u in (0,1) as input from your built-in random generator routine and get output values η in (0,1) that are distributed according to our majorant function. Since our majorant was already pretty close in shape to the exact PDF, we only need to apply relatively small subsequent corrections to the generated sample of η values by means of the Acceptance - Rejection Method by John Von Neumann
discussed in D).
All it takes is to throw another uniform random value uu in (0,1) and compare it with the
if uu <= f_by_ch
accept this particular η, else reject η and try with a new η!
Here is now the final plot of success. I prepared a histogram of binned η values of 73728 halo stars
sampled from the .Sci
code by means of the described algorithm. Then I overlaid the previous plot with the properly normalized C * majorant (red) and our exact PDF (green). The histogram shape and normalization agrees perfectly with the exact (green) PDF for all η in (0,1)!
This quite long discussion has so far just covered the random generation of halo star positions
in 3D space
. There is another long section to come (after some rest
) that deals with the random generation of the halo star magnitudes and colors
, in agreement with observation!
Please let me know whether the presented level was hopelessly too high, or whether people with some math/astrophysics background were getting out something from this report. The material is certainly not meant as a Newbie tutorial
. It rather should provide a concise report about some scientific-level aspects of celestia.Sci