Package enrichpy :: Module initialmassfunction

Module initialmassfunction

source code

Implements stellar initial mass functions from Romano et al. 2005.

Full Citation: Romano D., Chiappini C., Matteucci F., Tosi M., 2005, A&A, 430, 491 (2005A&A...430..491R)

http://adsabs.harvard.edu/abs/2005A&A...430..491R

Functions
 
by_number(imf_mass_function)
Convert an IMF by mass to an IMF by number.
source code
 
broken_powerlaw_function(limits, coefficients, powers)
Return a function consisting of seperate powerlaws.
source code
 
imf_mass_Scalo86(mass)
Scalo 1986 stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
imf_mass_Scalo98(mass)
Scalo 1998 stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
imf_mass_Tinsley(mass)
Tinsley stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
imf_mass_Kroupa(mass)
Kroupa stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
imf_mass_Salpeter(mass)
Salpeter stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
imf_mass_Chabrier(mass)
Chabrier stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
imf_number_Chabrier(mass)
Chabrier stellar Initial Mass function by mass via Romano et al. (2005).
source code
 
test_plot_norm_imf(imf_function)
Plot and test the normalization of imf_function.
source code
 
test_plot_norm_all() source code
Variables
  coeff_Scalo86 = array([ 0.19, 0.24])
  pow_Scalo86 = array([-1.35, -1.7 ])
  lims_Scalo86 = array([ 0., 2., Inf])
  coeff_Scalo98 = array([ 0.39, 0.39, 0.16])
  pow_Scalo98 = array([-0.2, -1.7, -1.3])
  lims_Scalo98 = array([ 0., 1., 10., Inf])
  coeff_Tinsley = array([ 0.21, 0.26, 2.6 ])
  pow_Tinsley = array([-1. , -1.3, -2.3])
  lims_Tinsley = array([ 0., 2., 10., Inf])
  coeff_Kroupa = array([ 0.58, 0.31, 0.31])
  pow_Kroupa = array([-0.3, -1.2, -1.7])
  lims_Kroupa = array([ 0. , 0.5, 1. , Inf])
  mc = 0.079
  logmc = -1.10237290871
  sigma = 0.69
  twosigsq = 0.9522
  fa = 0.85
  fb = 0.24
  __package__ = 'enrichpy'
Function Details

broken_powerlaw_function(limits, coefficients, powers)

source code 

Return a function consisting of seperate powerlaws.

F(x) = a_i x^{p_i} for l_i < x < l_{i+1}

Parameters

limits: array (length n+1)
boundaries of the specified powerlaws. Must be one greater in length than coefficents and powers. Specify -numpy.infty for the first limit or numpy.infty for the last limit for unbounded powerlaws.
coefficients: array (length n)
values of the coefficient a_i
powers: array (length n)
values of the powerlaw indices p_i

The returned function takes a single, one-dimensional array of values on which to operate.

imf_mass_Scalo86(mass)

source code 

Scalo 1986 stellar Initial Mass function by mass via Romano et al. (2005).

Parameters

mass: 1d array
mass in Solar masses.

imf_mass_Scalo98(mass)

source code 

Scalo 1998 stellar Initial Mass function by mass via Romano et al. (2005).

Parameters

mass: 1d array
mass in Solar masses.

imf_mass_Tinsley(mass)

source code 

Tinsley stellar Initial Mass function by mass via Romano et al. (2005).

Parameters

mass: 1d array
mass in Solar masses.

imf_mass_Kroupa(mass)

source code 

Kroupa stellar Initial Mass function by mass via Romano et al. (2005).

Parameters

mass: 1d array
mass in Solar masses.

imf_mass_Chabrier(mass)

source code 

Chabrier stellar Initial Mass function by mass via Romano et al. (2005).

With slope -1.3 above 1 MSun.

imf_number_Chabrier(mass)

source code 

Chabrier stellar Initial Mass function by mass via Romano et al. (2005).

With slope -1.3 above 1 MSun.