Issue 
A&A
Volume 635, March 2020



Article Number  A20  
Number of page(s)  8  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201937202  
Published online  28 February 2020 
A precise analytical approximation for the deprojection of the Sérsic profile
Institut d’Astrophysique de Paris (UMR 7095: CNRS & Sorbonne Université), 98 Bis Bd Arago, 75014 Paris, France
email: vitral@iap.fr
Received:
27
November
2019
Accepted:
8
January
2020
The Sérsic model shows a close fit to the surface brightness (or surface density) profiles of elliptical galaxies and galaxy bulges, and possibly also those of dwarf spheroidal galaxies and globular clusters. The deprojected density and mass profiles are important for many astrophysical applications, in particular for massorbit modeling of these systems. However, the exact deprojection formula for the Sérsic model employs special functions that are not available in most computer languages. We show that all previous analytical approximations to the 3D density profile are imprecise at low Sérsic index (n ≲ 1.5). We derived a more precise analytical approximation to the deprojected Sérsic density and mass profiles by fitting twodimensional tenthorder polynomials to the residuals of the analytical approximations by Lima Neto et al. (1999, MNRAS, 309, 481; LGM) for these profiles, relative to the numerical estimates. Our LGMbased polynomial fits have typical relative precision better than 0.2% for both density and mass profiles, for Sérsic indices 0.5 ≤ n ≤ 10 and radii 0.001 < r/R_{e} < 1000. Our approximation is much more precise than previously published approximations (except, in some models, for a few discrete values of the index). An appendix compares the deprojected Sérsic profiles with those of other popular simple models.
Key words: methods: numerical / galaxies: elliptical and lenticular, cD / galaxies: structure / globular clusters: general / galaxies: bulges
© E. Vitral and G. A. Mamon 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The Sérsic model (Sérsic 1963; Sersic 1968) is the generalization of the R^{1/4} law (de Vaucouleurs 1948) to describe the surface brightness profiles of elliptical galaxies (Caon et al. 1993) and the bulges of spiral galaxies (Andredakis et al. 1995). It has also been used to describe the surface density profiles of nuclear star clusters (Graham & Spitler 2009), resolved dwarf spheroidal galaxies (Battaglia et al. 2006), and globular clusters (Barmby et al. 2007).
The surface (mass or number) density (or equivalently surface brightness) of the Sérsic model is
where R is the projected distance to the source center, R_{e} is the effective radius containing half of the projected luminosity, n is the Sérsic index, and Σ_{0} is the central surface density. The term b_{n} is a function of n, obtained by solving the equation:
where is the lower incomplete gamma function.
Since the Sérsic model accurately represents astronomical objects viewed in projection, it is important to know its corresponding threedimensional (3D) density and mass profiles. These serve as a reference for comparison with other possible observational tracers, as well as to dark matter. Moreover, the 3D density profile is required for modeling the kinematics of spherical structures because it appears in the Jeans equation of local dynamical equilibrium. Since the Jeans equation also contains the total mass profile, the 3D mass profiles of stellar components are required to estimate the dark matter mass profile of elliptical and dwarf spheroidal galaxies.
Many authors assume that simple threedimensional models resemble Sérsic models for certain values of the Sérsic index: It is often assumed that massive ellipticals and spiral bulges are well represented by the Hernquist (1990) model (e.g., Widrow & Dubinski 2005). On the other hand, dwarf spheroidal galaxies are often described with the Plummer (1911) model (e.g., Muñoz et al. 2018, who also tried Sérsic and other models), while ultra diffuse galaxies have been described with the Einasto (Einasto 1965; Navarro et al. 2004)^{1} model (Nusser 2019). Łokas & Mamon (2001) noted that the projected Navarro et al. (1996, hereafter NFW) model resembles an n = 3 Sérsic for reasonable concentrations. Finally, n = 4 Sérsic models are considered to resemble the Jaffe (1983) model (Ciotti et al. 2019). In Appendix A, we compare these models to the deprojected Sérsic.
Unfortunately, the deprojection of the Sérsic surface density profile to a 3D (mass or number)^{2} density profile, through Abel (1826) inversion
(e.g., Binney & Mamon 1982), as well as the corresponding 3D mass (or number) profile
both involve the complicated Meijer G special function (Mazure & Capelato 2002 for integer values of n, and Baes & Gentile 2011 for general values of n) or the other, complicated Fox H function (Baes & van Hese 2011), neither of which are available in popular computer languages.
Following the shape of the analytical approximation to the R^{1/4} law by Mellier & Mathez (1987), Prugniel & Simien (1997, hereafter, PS) proposed an analytical approximation for the 3D density of the Sérsic profile:
where p_{n} is a function depending only on n:
Equation (5) yields a simple analytical form for the 3D mass profile,
Lima Neto et al. (1999, hereafter LGM) later perfected the approximation of Eq. (6) with
According to LGM, Eq. (9) has 5% relative accuracy for 0.56 ≤ n ≤ 10 and −2 < log(r/R_{e}) < 3. However, the powerlaw approximation at small radii is unjustified for small n. Indeed, as shown by Baes & Gentile (2011), the central density profile converges to a finite value for n < 1 (and the inner density profile diverges only logarithmically for n = 1), as we illustrate in Sect. 3.
Simonneau & Prada (1999, 2004, hereafter SP) proposed the quasiGaussian expansion for the density profile
where
where x_{j} and w_{j} are ten fit parameters. The individual SP density profiles (the terms inside the sum of Eq. (10)) have a similar (but not the same) form to that of the PS/LGM one, hence the similar shape of the mass profile:
Trujillo et al. (2002) proposed an ellipsoid formula, which in the limit of spherical symmetry becomes
where K_{ν}(x) is the modified Bessel function of the second kind^{3} of index ν, while ν_{n}, p_{n}, a_{n, 0}, a_{n, 1}, and a_{n, 2} are empirical functions of index n. Trujillo et al. (2002) only provided their results for integer and halfinteger values of n for n ≤ 5 and only integer values of n beyond. Emsellem & van de Ven (2008, hereafter EV) repeated their analysis on a finer grid of n, with steps of 0.1 for 0.5 ≤ n ≤ 1.5 and with one more term, a_{n, 3}, in Eq. (14) involving 168 parameters. However, as shown in Fig. 1, these functions vary abruptly for n ≲ 1.2. Moreover, neither Trujillo et al. (2002) nor Emsellem & van de Ven (2008) provide analytical forms for the mass profile.
Fig. 1. Variation with Sérsic index of the different parameters of the analytical approximation of Emsellem & van de Ven (2008) for the deprojected Sérsic density profile (filled circles). The solid and dotted curves show the spline cubic and linear interpolations, respectively. At small n, the parameters vary abruptly and the interpolations (both linear and cubic) are therefore uncertain. 
In summary, all the previous approximations to the deprojected Sérsic model have drawbacks:

Both PS and LGM (and Márquez et al. 2000, which is the same as LGM, but with a slightly different last term for p_{n}, which was a typo) are inappropriate for low n (Baes & Gentile 2011) and less precise than claimed (Emsellem & van de Ven 2008).

SP is limited to n ≥ 1, and is generally less precise than EV.

Trujillo et al. (2002) is only given for halfinteger values of n and their parameters vary wildly with n for n ≤ 1.5. These authors do not provide a formula for the mass profile.

EV also suffers from discrete values of n, even though the grid is finer (Δn = 0.1 for n ≤ 1.5). These authors also did not provide a formula for the mass profile.
In this article, we provide polynomial fits to the log residuals of the LGM approximation, which allow high accuracy to be reached for both the 3D density and 3D mass profiles in a wide range of Sérsic indices. In Sect. 2, we present the mathematical formalism and briefly explain our numerical integration method. We then show in Sect. 3 how our polynomial plus LGM approximation is orders of magnitude more precise than the formulae of LGM, SP, and Trujillo et al. (2002), as well as that of EV for low n, and is only slightly less precise for n ≳ 3. We conclude and discuss our results in Sect. 4.
2. Method
2.1. Equations using dimensionless profiles
We express the general surface density, 3D density, and 3D mass (or number) profiles in terms of dimensionless functions:
Hereafter, we use x = r/R_{e} and X = R/R_{e}. For the Sérsic model (see Graham & Driver 2005 for a thorough review of the Sérsic profile), the dimensionless surface density profile is
while for the PS model, one can write the dimensionless 3D density and mass profiles as
It is easy to show that the deprojection Eq. (3) becomes
where
The dimensionless mass profile is
where Eq. (24) is obtained by inversion of the order of integration in the second equality of Eq. (23).
2.2. Numerical integration
We numerically evaluated the dimensionless 3D density (Eq. (21)) and mass (Eq. (24)) profiles by performing the numerical integrations in cells 50 × 100 of [log n, log(r/R_{e})], with log 0.5 ≤ log n ≤ 1 and −3 ≤ log(r/R_{e}) ≤ 3. Numerical calculations were done with Python’s SCIPY.INTEGRATE.QUAD. For both density and mass profiles, we split the numerical integration in two, that is,
where and a ≤ X_{crit} ≤ b. We used a relative tolerance of epsrel = 10^{−4} and limit = 1000 in both integrals. If X_{crit} ∉ [a, b], we also used epsrel = 10^{−4} and limit = 1000, but for a single integral from a to b.
We performed our analysis using either the highly accurate approximations for b_{n} of Ciotti & Bertin (1999, hereafter, CB) or the exact (numerical) solutions of Eq. (2). We noticed that the difference between these two approaches was negligible (see Sect. 3).
We then fit twodimensional polynomials to both and , for geometrically spaced x and n, writing
with polynomial orders 2 ≤ k ≤ 12. For this, we used Python’s package NUMPY.LINALG.LSTSQ. We found the smallest residuals for tenthorder polynomials when using both the b_{n} approximation of CB and b_{n} by numerically solving Eq. (2). The coefficients are provided in Tables B.1–B.4. In the rest of the paper, we present the results relative to the CB approximation, since it is a simpler and more used model, and also because our tenthorder polynomial fits the exact b_{n} case remarkably well.
2.3. Numerical precision: tests for known simple analytical deprojections (n = 0.5 and 1)
For Sérsic indices n = 0.5 and n = 1, there are analytical solutions for the 3D density profile:
where K_{0}(x) is the modified Bessel function of the second kind of index 0. We can therefore verify the numerical integration of Eq. (21) for these two Sérsic indices.
For the interval −3 ≤ log(r/R_{e}) ≤ 3, we compared the densities from numerical integration with the analytical formulae of Eq. (27) using the CB approximation for b_{n}. The match is very good, with rootmeansquare (rms) values of of 1.5 × 10^{−7} and 2 × 10^{−8} for n = 0.5 and n = 1, respectively. The same comparison using the exact b_{n} yields 7 × 10^{−5} and 2 × 10^{−8}, respectively (with one particular value of r causing the higher rms for n = 0.5).
3. Results
As seen in Fig. 2, the 3D density profiles depart from the power laws proposed by LGM at low n, especially for low radii, as expected by the asymptotic expansions of Baes & Gentile (2011) for n < 1. Interestingly, the LGM formula is also inadequate at low radii for n = 1.25 and 2.25, although the asymptotic expansion of Baes & Gentile (2011) indicate powerlaw behavior at small radii. This poor accuracy of the LGM approximation at low radii is a serious concern when performing kinematic modeling of systems with possible central massive black holes. For example, Gaia second data release (DR2) positions and proper motions for stars in nearby globular clusters extend inwards to 0.7 arcsec from the center, which translates to 0.002 R_{e}.
Fig. 2. Illustration of the accuracy of the PS formula with the LGM coefficients for p_{n}. The solid curves show the numerically estimated profiles, while the coloreddashed curves show the LGM approximation. 
We now compare the accuracy of the different analytical approximations for the 3D density and 3D mass profiles. Figure 3 displays the ratio , for and , for the main analytical approximations available in the literature, along with our new model
where f is either the 3D density or 3D mass profile. We see that our model presents a more continuous behavior over the full range of Sérsic indices and radii. Our approximation displays the smallest residuals among all models for n ≲ 3 (except that SP outperforms our model for mass estimates at r > 3 R_{e} for n > 1.3).
The variation of accuracy with Sérsic index can be seen in more detail in Fig. 4, which displays the rms of , over the radial domain where ρ(r) > 10^{−30} ρ(R_{e}), of the main analytical approximations using 1000 logspaced Sérsic indices. Figure 4 indicates that the SP and EV approximations for density are less accurate than 2.3% (0.01 for ) for n < 1.6 and n < 1.3, respectively. Our approximation (Eq. (28)) is more accurate than SP for n < 4.3 (density) and n < 3.1 (mass), and is more accurate than EV for n < 3.4 (density), except for their particular choices of n. Figure 4 shows that the EV approximation is much more accurate at specific values of n (we note that our grid does not contain all of these values precisely, and therefore the EV approximation is even more accurate at these specific values of n). However, these specific values of n represent a negligible measure compared to the full continuous range of 0.5 ≤ n ≤ 10. Therefore, the EV approximation at low n is not reliable for estimating the 3D density profile.
Fig. 3. Accuracy of deprojected density (top 6 panels) and mass (3 bottom panels) of the different analytical approximations (PS: Prugniel & Simien 1997; LGM: Lima Neto et al. 1999; SP: Simonneau & Prada 1999, 2004; Trujillo+02: Trujillo et al. 2002; EV: Emsellem & van de Ven 2008) and our new one (Eq. (28), with greencolored titles) as a function of both Sérsic index (abscissae) and radius (ordinates). The color scale given in the vertical color bars is linear for log ratios between −0.001 and 0.001 and logarithmic beyond. The gray region and green curves in the upper left of the density panels are for regions where the numerical integration reached the underflow limit or density 10^{−30} times ρ(R_{e}), respectively, because of the very rapid decline of density at large radii for low n, and also covers n < 1, which is not covered by the SP model. We note that the EV and Trujillo+02 models perform better at specific values of n that are often missed in our grid. 
Fig. 4. Accuracy of different approximations (LGM: Lima Neto et al. 1999; SP: Simonneau & Prada 1999, 2004; EV: Emsellem & van de Ven 2008) and our new one (Eq. (28)) as a function of Sérsic index. We note that the EV model performs better at specific values of n that are often missed in our logarithmic grid of 1000 values of n. 
We analyzed the results shown in Fig. 4 using b_{n} from CB or by numerically solving Eq. (2), and the results were very similar. In fact, the results are similar if we adopt one form of b_{n} in the numerical integration and the other in the analytical approximations. This can be explained by the fact that is practically the same for both estimates of b_{n}, yielding a very similar fit of Eq. (26).
Finally, Table 1 provides the rms accuracies computed over the full range of radii −3 ≤ log(r/R_{e}) ≤ 3 and 0.5 ≤ n ≤ 10, except for the SP formula, which does not allow n ≤ 1, and also avoiding the domain where ρ(r) < 10^{−30} ρ(R_{e}). We see that, averaging over all Sérsic indices, our approximation is much more accurate than all others (with over ten times lower rms).
Accuracy of approximations to 3D density and mass profiles.
4. Conclusions and discussion
The Sérsic model is usually considered to provide excellent fits to the surface density (or surface brightness) profiles of elliptical galaxies, spiral bulges, and even dwarf spheroidal galaxies and globular clusters. In the past, many authors have used simple analytical models to describe these systems, arguing that their models, once projected, resemble Sérsic models. It is more relevant to compare the physically meaningful threedimensional density profiles of these simple models to the deprojected Sérsic model.
This comparison is made in Appendix A for the Plummer, Jaffe, Hernquist, Einasto, and NFW models. As seen in Fig. A.1, most of the simple models do not provide good fits to the deprojected Sérsic model, even for narrow ranges of the Sérsic index. The Plummer model requires a low index at small radii, but a much higher index at large radii, and the normalized density profile fits poorly at most radii. The Hernquist model resembles the n = 2.8 deprojected Sérsic model at low radii and the n = 5.7 Sérsic at large radii. The NFW models resemble the n = 2.8 deprojected Sérsic at low radii (consistent with the similarity of the projected Sérsic with NFW discovered by Łokas & Mamon 2001), but have a shallower slope at large radii than even the shallowest (n = 8) deprojected Sérsic model. On the other hand, the Jaffe model resembles the n = 5.7 model at all radii. Moreover, as seen in Fig. A.2, the Einasto model provides a fair representation (rms difference of density profiles normalized to value at halfmass radius less than 0.1 dex) of the deprojected Sérsic model for n > 6.5.
We reconsidered the different analytical deprojections of the Sérsic surface brightness (or surface density) profile. We found that the analytical approximations present in the literature do not show satisfying results when the Sérsic index is in the range 0.5 ≤ n ≲ 1.5 (apart from the specific values of n given by Emsellem & van de Ven 2008). In particular, the powerlaw times exponential density profile of Prugniel & Simien (1997) and Lima Neto et al. (1999) fails to reproduce the inner density profiles for low n, even up to n = 2.25 despite the powerlaw behavior expected at small radii for n > 1 (Baes & Gentile 2011).
With a tenthorder twodimensional polynomial fit, we propose a new analytical approximation (Eq. (28)) that is precise over the range log0.5 ≤ log n ≤ 1, for −3 ≤ log(r/R_{e}) ≤ 3. Our approximation provides the highest precision when averaging over all values of Sérsic indices and radii (Table 1). While the approximations of Simonneau & Prada (1999, 2004) on one hand and of Emsellem & van de Ven (2008) on the other are more accurate than ours for n > 4.3 and 3.4, respectively, ours is more accurate at lower Sérsic indices.
This is important for the study of astronomical sources with low Sérsic indexes, such as galaxy bulges, nuclear star clusters, dwarf spheroidal galaxies, and globular clusters. Moreover, our approximation of Eq. (28) to the density profile is sufficiently accurate for most scientific analyses for n > 3. Nevertheless, the user could use a hybrid approximation, combining either the Simonneau & Prada (2004) or Emsellem & van de Ven (2008) approximations for n ≥ 3.4 and ours for n < 3.4 (as shown in the last rows of Table 1). Finally, our analysis has the advantage of also providing a precise approximation for the mass profile, whereas no analytical expression can be derived from the density profile of Emsellem & van de Ven (2008). Our Python 3 codes are available at^{4} along with coefficients of Tables B.1 and B.2.
These results will be useful in future massorbit modeling analyses of lowmass spherical systems, as we are preparing for globular clusters (Vitral & Mamon, in prep.).
Navarro et al. (2004) showed how the Einasto model accurately represents the density profiles of dark matter halos in dissipationless cosmological simulations, while Merritt et al. (2005) were the first to note its similar form to the Sérsic model, and Merritt et al. (2006) were first to realize that this model had been previously introduced by Einasto.
Trujillo et al. (2002) call this the modified Bessel function of the third kind, as some others do.
Acknowledgments
We thank Alister Graham for providing important references.
References
 Abel, N. H. 1826, Journal für die reine und angewandte Mathematik, 1, 153 [Google Scholar]
 Andredakis, Y. C., Peletier, R. F., & Balcells, M. 1995, MNRAS, 275, 874 [NASA ADS] [CrossRef] [Google Scholar]
 Baes, M., & Gentile, G. 2011, A&A, 525, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baes, M., & van Hese, E. 2011, A&A, 534, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barmby, P., McLaughlin, D. E., Harris, W. E., Harris, G. L. H., & Forbes, D. A. 2007, AJ, 133, 2764 [Google Scholar]
 Battaglia, G., Tolstoy, E., Helmi, A., et al. 2006, A&A, 459, 423 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Binney, J., & Mamon, G. A. 1982, MNRAS, 200, 361 [NASA ADS] [CrossRef] [Google Scholar]
 Caon, N., Capaccioli, M., & D’Onofrio, M. 1993, MNRAS, 265, 1013 [NASA ADS] [CrossRef] [Google Scholar]
 Ciotti, L., & Bertin, G. 1999, A&A, 352, 447 [NASA ADS] [Google Scholar]
 Ciotti, L., Mancino, A., & Pellegrini, S. 2019, MNRAS, 490, 2656 [NASA ADS] [CrossRef] [Google Scholar]
 de Vaucouleurs, G. 1948, Annales d’Astrophysique, 11, 247 [NASA ADS] [Google Scholar]
 Einasto, J. 1965, Trudy Inst. Astroz. AlmaAta, 51, 87 [Google Scholar]
 Emsellem, E., & van de Ven, G. 2008, ApJ, 674, 653 [CrossRef] [Google Scholar]
 Graham, A. W., & Driver, S. P. 2005, PASA, 22, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Graham, A. W., & Spitler, L. R. 2009, MNRAS, 397, 2148 [NASA ADS] [CrossRef] [Google Scholar]
 Hernquist, L. 1990, ApJ, 356, 359 [NASA ADS] [CrossRef] [Google Scholar]
 Jaffe, W. 1983, MNRAS, 202, 995 [NASA ADS] [Google Scholar]
 Lima Neto, G. B., Gerbal, D., & Márquez, I. 1999, MNRAS, 309, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Łokas, E. L., & Mamon, G. A. 2001, MNRAS, 321, 155 [NASA ADS] [CrossRef] [Google Scholar]
 Márquez, I., Lima Neto, G. B., Capelato, H., Durret, F., & Gerbal, D. 2000, A&A, 353, 873 [NASA ADS] [Google Scholar]
 Mazure, A., & Capelato, H. V. 2002, A&A, 383, 384 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mellier, Y., & Mathez, G. 1987, A&A, 175, 1 [NASA ADS] [Google Scholar]
 Merritt, D., Navarro, J. F., Ludlow, A., & Jenkins, A. 2005, ApJ, 624, L85 [NASA ADS] [CrossRef] [Google Scholar]
 Merritt, D., Graham, A. W., Moore, B., Diemand, J., & Terzić, B. 2006, AJ, 132, 2685 [NASA ADS] [CrossRef] [Google Scholar]
 Muñoz, R. R., Côté, P., Santana, F. A., et al. 2018, ApJ, 860, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Hayashi, E., Power, C., et al. 2004, MNRAS, 349, 1039 [NASA ADS] [CrossRef] [Google Scholar]
 Nusser, A. 2019, ArXiv eprints [arXiv:1907.08035] [Google Scholar]
 Plummer, H. C. 1911, MNRAS, 71, 460 [CrossRef] [Google Scholar]
 Prugniel, P., & Simien, F. 1997, A&A, 321, 111 [NASA ADS] [Google Scholar]
 Sérsic, J. L. 1963, Bull. Assoc. Argent. Astron., 6, 41 [NASA ADS] [Google Scholar]
 Sersic, J. L. 1968, Atlas de Galaxias Australes (Cordoba, Argentina: Observatorio Astronomico) [Google Scholar]
 Simonneau, E., & Prada, F. 1999, ArXiv eprints, unpublished [arXiv:astroph/9906151] [Google Scholar]
 Simonneau, E., & Prada, F. 2004, Rev. Mex. Astron. Astrofis., 40, 69 [NASA ADS] [Google Scholar]
 Trujillo, I., Asensio Ramos, A., RubiñoMartín, J. A., et al. 2002, MNRAS, 333, 510 [Google Scholar]
 Widrow, L. M., & Dubinski, J. 2005, ApJ, 631, 838 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Comparison of deprojected Sérsic to other popular models
Figure A.1 compares the density profiles, normalized to the halfmass radius r_{h}, for which we applied the following relations:
Fig. A.1. Comparison of other known 3D density profiles and the deprojected Sérsic density profile for certain values of the Sérsic index n. All density profiles are normalized to the value at the 3D halfmass radius, r_{h} (see text). The different NFW models can be distinguished at low radii, where the density increases with r_{max}/a. 
The ratio of halfmass radius to scale radius r_{h}/a is given by
In Eq. (A.2) for NFW, c = r_{max}/a, where r_{max} is the maximum allowed radius (because, contrary to all other models discussed here, the NFW model has logarithmically divergent mass). Also, for Einasto, P^{(−1)}(a, y) is the inverse regularized lower incomplete gamma function, i.e., x = P^{(−1)}(a, y) satisfies γ(a, x)/Γ(a) = y^{5}. For Sérsic, the conversion was done by fitting a thirdorder polynomial and recovering the relation , where {a_{0}, a_{1}, a_{2}, a_{3}}={1.32491, 0.0545396, − 0.0286632, 0.0035086}.
The Einasto model, which is the 3D analog of the Sérsic model, resembles the deprojected Sérsic model. Figure A.2 shows the bestfit values of the Einasto index, n_{Einasto} in terms of the Sérsic index. The relation (red curve) is almost onetoone (dashed line). The figure also shows the rms over all radii and bestfit indices (blue curve).
Fig. A.2. Comparison of Einasto and deprojected Sérsic density profiles (both normalized to halfmass radius). Red: bestfit Einasto index (dashed line is n_{Einasto} = n_{Sersic}). Blue: rms of best fit. 
Appendix B: Coefficients of polynomials for new deprojected Sérsic models
In this section, we present Tables B.1–B.4, which contain the coefficients a_{ij} in Eq. (26), for and , as well as for both Ciotti & Bertin (1999) approximation for b_{n} and the exact solution from Eq. (2). The numbers in parentheses are the exponents: e.g., “1.234 (−3)” corresponds to 1.234 × 10^{−3}. Coefficients that are not followed by a number in parentheses have exponent zero.
Coefficients of Eqs. (26) and (28), for and b_{n} calculated from Ciotti & Bertin (1999) accurate approximation.
Coefficients of Eqs. (26) and (28), for and b_{n} calculated from Ciotti & Bertin (1999) accurate approximation.
Coefficients of Eqs. (26) and (28), for and b_{n} calculated from Eq. (2).
Coefficients of Eqs. (26) and (28), for and b_{n} calculated from Eq. (2).
All Tables
Coefficients of Eqs. (26) and (28), for and b_{n} calculated from Ciotti & Bertin (1999) accurate approximation.
Coefficients of Eqs. (26) and (28), for and b_{n} calculated from Ciotti & Bertin (1999) accurate approximation.
All Figures
Fig. 1. Variation with Sérsic index of the different parameters of the analytical approximation of Emsellem & van de Ven (2008) for the deprojected Sérsic density profile (filled circles). The solid and dotted curves show the spline cubic and linear interpolations, respectively. At small n, the parameters vary abruptly and the interpolations (both linear and cubic) are therefore uncertain. 

In the text 
Fig. 2. Illustration of the accuracy of the PS formula with the LGM coefficients for p_{n}. The solid curves show the numerically estimated profiles, while the coloreddashed curves show the LGM approximation. 

In the text 
Fig. 3. Accuracy of deprojected density (top 6 panels) and mass (3 bottom panels) of the different analytical approximations (PS: Prugniel & Simien 1997; LGM: Lima Neto et al. 1999; SP: Simonneau & Prada 1999, 2004; Trujillo+02: Trujillo et al. 2002; EV: Emsellem & van de Ven 2008) and our new one (Eq. (28), with greencolored titles) as a function of both Sérsic index (abscissae) and radius (ordinates). The color scale given in the vertical color bars is linear for log ratios between −0.001 and 0.001 and logarithmic beyond. The gray region and green curves in the upper left of the density panels are for regions where the numerical integration reached the underflow limit or density 10^{−30} times ρ(R_{e}), respectively, because of the very rapid decline of density at large radii for low n, and also covers n < 1, which is not covered by the SP model. We note that the EV and Trujillo+02 models perform better at specific values of n that are often missed in our grid. 

In the text 
Fig. 4. Accuracy of different approximations (LGM: Lima Neto et al. 1999; SP: Simonneau & Prada 1999, 2004; EV: Emsellem & van de Ven 2008) and our new one (Eq. (28)) as a function of Sérsic index. We note that the EV model performs better at specific values of n that are often missed in our logarithmic grid of 1000 values of n. 

In the text 
Fig. A.1. Comparison of other known 3D density profiles and the deprojected Sérsic density profile for certain values of the Sérsic index n. All density profiles are normalized to the value at the 3D halfmass radius, r_{h} (see text). The different NFW models can be distinguished at low radii, where the density increases with r_{max}/a. 

In the text 
Fig. A.2. Comparison of Einasto and deprojected Sérsic density profiles (both normalized to halfmass radius). Red: bestfit Einasto index (dashed line is n_{Einasto} = n_{Sersic}). Blue: rms of best fit. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.