Package to calculate various profiles for deprojected, flattened (or elongated) Sérsic mass distributions, including: enclosed mass, circular velocity, density, log density slope, surface density, and projected enclosed mass.
These calculations follow and extend the derivation of rotation curves for flattened Sérsic bulges presented by Noordermeer, 2008, MNRAS, 385, 1359. Further details about the calculations included in this package are described in Price et al., 2022, A&A 665 A159.
Please see the documentation for this package for detailed information about installation, usage, and to download the set of pre-computed Sérsic profile tables.
import os
import numpy as np
import deprojected_sersic_models as deproj_sersic
table_dir = os.getenv('DEPROJECTED_SERSIC_MODELS_DATADIR')
# Sersic profile properties & radius array
total_mass = 1.e11
Reff = 5.0
n = 1.0
invq = 5. # Oblate, q = c/a = 0.2
r = np.arange(0., 30.1, 0.1)
# Load & interpolate all profiles in saved table:
table_interp = deproj_sersic.interpolate_entire_table(r=r, total_mass=total_mass,
Reff=Reff, n=n, invq=invq,
path=table_dir)- numpy
- scipy
- matplotlib
- astropy
- dill
If you use this package or the precomputed profile tables in a publication, please cite Price et al., 2022, A&A 665 A159 (ADS | arXiv).
This project is Copyright (c) Sedona Price / MPE IR/Submm Group and licensed under the terms of the BSD 3-Clause license. See the LICENSE.rst for more information.