fit_modified_blackbody
Fits a single-temperature modified blackbody to SED fluxes in GHz and mJy, returning dust mass, temperature, beta, and covariance errors for SED modeling.
Instructions
Fit a single-temperature modified blackbody to an SED (observed frequencies in GHz, fluxes in mJy).
form='thin': S = (1+z) M kappa(nu_rest) B(nu_rest,T) / D_L^2. form='general': S = (1+z) A (1 - exp(-tau)) B / D_L^2 with tau = (nu_rest/nu_0)^beta, nu_0 = c/lambda0_um, and A = M kappa(nu_0) (tau = 1 at lambda0 defines the area). The mass in this form depends on that area assumption. Free parameters: log10 M_dust (Msun), T_d (K), beta (unless fixed). Errors in quadrature with flux_cal_frac. L_IR is integrated over rest 8-1000 um (and L_FIR over 42.5-122.5 um) from the intrinsic B(T) model. cmb_contrast subtracts B(T_CMB(z)) (da Cunha 2013 contrast term only; T_d is then the heated temperature). Returns least-squares values with covariance errors; mcmc=True adds an emcee posterior (needs [fit] extra). Limits: one temperature misses cold dust mass; T and beta are degenerate; upper limits are not handled.
Tool call template
Replace each placeholder with a value matching the input schema. Only required parameters are included; add optional parameters as needed.
{
"id": 1,
"jsonrpc": "2.0",
"method": "tools/call",
"params": {
"arguments": {
"freq_ghz": "<freq_ghz>",
"flux_mjy": "<flux_mjy>",
"flux_err_mjy": "<flux_err_mjy>"
},
"name": "fit_modified_blackbody"
}
}Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| form | No | thin | |
| mcmc | No | ||
| seed | No | ||
| nu0_ghz | No | ||
| fix_beta | No | ||
| flux_mjy | Yes | ||
| freq_ghz | Yes | ||
| redshift | No | ||
| lambda0_um | No | ||
| mcmc_steps | No | ||
| distance_pc | No | ||
| cmb_contrast | No | ||
| fix_t_dust_k | No | ||
| flux_err_mjy | Yes | ||
| kappa0_cm2_g | No | ||
| kappa_preset | No | galaxy_850um | |
| mcmc_walkers | No | ||
| flux_cal_frac | No |