Modeling options#
There are several modeling options that can be used for the ModifiedBlackbody fitting.
All of these follow the form of the Casey (2012) and/or Drew & Casey (2022) models,
with optional mid-infrared powerlaw (slope defaults to alpha=2.0) and the choice between assuming optically thin dust or the general opacity model (turnover wavelength defaults to l0 = 200, in microns).
The model includes the effects of the CMB as described in Da Cunha et al. (2013).
To choose a given model, use the pl (powerlaw), opthin (optically thin), and/or pl_piecewise (join power law piecewise?) flags when intializing a ModifiedBlackbody instance.
The documentation for the overall model is below:
- mbb.mbb_funcs.mbb_func(l, N=12, beta=1.8, T=35, z=0, alpha=2.0, l0=200, opthin=True, pl=False, pl_piecewise=False)#
MBB function with optional powerlaw and variable opacity assumptions
- Parameters:
l (float) – the rest-frame wavelengths, in microns, at which to evaluate the model.
N (float) – log10 normalization factor of blackbody
beta (float) – emissivity index
T (float) – effective dust temperature
z (float) – the redshift of the model
alpha (float) – power-law slope
l0 (float) – turnover wavelength at which dust is optically thin
opthin (bool) – Whether or not the model should assume optically thin dust emission.
pl (bool) – Whether or not the model should include a MIR power law
pl_piecewise (bool) – if the powerlaw should be attached piecewise (as in Drew&Casey 2022) or fall off exponentially (as in Casey+ 2012)
- Returns:
the value(s) of the model in Jy at wavelengths
l, in microns- Return type:
float