rotational_diagram
Fit LTE line intensities to estimate rotational temperature and total column density, optionally correcting for opacity and using partition tables.
Instructions
LTE rotational (population) diagram: fit ln(N_u/g_u) = ln(N_tot/Q(T_rot)) - E_u/T_rot.
Each line dict: freq_ghz, a_ul (s^-1), e_u_k (K), g_u, w_k_kms (integrated intensity, K km/s) and w_err_k_kms. Instead of w_k_kms you can give flux_jy_kms with bmaj_arcsec, bmin_arcsec (converted with the RJ 1.222e6 factor). flux_jy_kms must be the per-beam integrated intensity in Jy/beam km/s (moment 0 at the position), not a spatially integrated flux in Jy km/s, unless the source is smaller than the beam. For tau_correction also give fwhm_kms per line. N_u = 8 pi k nu^2 W / (h c^3 A_ul) / beam_filling (thin, RJ, background neglected). Errors on W get flux_cal_frac in quadrature. Q(T) from partition_table (log-log interpolation; include vibrational states) or a linear rotor approximation Q = kT/(hB) + 1/3. tau_correction fits ln N_tot and T_rot jointly with each line's thin-assumption column modelled as N_u (1 - exp(-tau))/tau, tau the Gaussian peak opacity from N_u, T_rot and fwhm_kms (Goldsmith and Langer 1999 C_tau). Limits: needs >= 3 lines over a range of E_u; T_rot is a lower limit to T_kin for subthermal gas; beam filling < 1 and opacity make N a lower limit; blends bias it.
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": {
"lines": "<lines>"
},
"name": "rotational_diagram"
}
}Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| lines | Yes | ||
| beam_filling | No | ||
| flux_cal_frac | No | ||
| tau_correction | No | ||
| partition_table | No | ||
| linear_rotor_b_ghz | No |