smaug: optically-thin, spherically-symmetric thermal plasma.

This model performs an analytical deprojection of an extended, optically-thin and spherically-symmetric source. A thorough description of the model is given in Pizzolato et al. (2003). In this model the 3D distributions of hydrogen, metals and temperature throughout the source are given specific functional forms dependent on a number of parameters, whose values are determined by the fitting procedure. The user has to extract the spectra in annular sectors, concentric about the emission peak. The inner boundary (in arcmin), the outer by the fitting procedure. The user has to extract the spectra in annular sectors, concentric about the emission peak. Three additional XFLTnnnn keywords must be added (e.g. with the ftools fkeypar). These should take the values “inner: x”, “outer: y”, “width: z” where x, y, z are the inner boundary (in arcmin), the outer boundary (also in arcmin), and the width (in degrees), respectively, of each annular sector. Some parameters of smaug define the redshift and other options (see below). The other, “relevant” ones define the 3D distributions of hydrogen density, temperature and metal abundance, determined by a simultaneous fit of the spectra. The cosmological parameters can be set using the cosmo command.

par1	central temperature [keV]
par2	max difference of temperature [keV]
par3	exponent of the inner temperature
par4	radius of the inner temperature [Mpc]
par5	exponent of the middle temperature
par6	radius of the middle temperature [Mpc]
par7	exponent of the outer temperature
par8	radius of the outer temperature [Mpc]
par9	central hydrogen density [cm$^{-3}$]
par10	fraction of nH.cc relative to the 1st beta component
par11	exponent of the first beta component
par12	radius of the 1st beta component [Mpc]
par13	exponent of the 2nd beta component
par14	radius of the 2nd beta component [Mpc]
par15	central metallicity [solar units]
par16	exponent of the metal distribution
par17	radius of the metal distribution [Mpc]
par18	redshift of the source
par19	number of mesh-points of the dem summation grid
par20	cutoff radius for the calculation [Mpc]
par21	mode of spectral evaluation: 0 = calculate, 1 = interpolate, 2 = APEC interpolate
par22	type of plasma emission code, 1 = Raymond-Smith, 2 = Mekal, 3 = Meka, 4 = APEC
K	model normalisation (nH.cc squared [cm$^{-6}$] )

Note that if the interactive chattiness level in XSPEC is set to a value >10, smaug also prints on screen the following quantities:

H$_0$	Hubble constant [km/s/Mpc]
q$_0$	deceleration parameter
L$_0$	cosmological constant
DA	source angular distance [Mpc]
DSET	dataset no. to which the quantities listed below are
IN	inner rim of the projected annular sector [Mpc]
OUT	outer rim of the projected annular sector [Mpc]
WID	width of the projected annular sector [deg]
EVOL	emitting volume within the integration radius cutoff [Mpc$^3$]
EINT	emission integral within the integration radius cutoff [ Mpc$^3$/cm$^6$]. If nH.cc is frozen to 1, the actual EI is obtained by multiplying this figure by the square root of the model normalisation

HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public

Last modified: Friday, 23-Aug-2024 13:20:40 EDT

HEASARC Staff Scientist Position - Applications are now being accepted for a Staff Scientist with significant experience and interest in the technical aspects of astrophysics research, to work in the High Energy Astrophysics Science Archive Research Center (HEASARC) at NASA Goddard Space Flight Center (GSFC) in Greenbelt, MD. Refer to the AAS Job register for full details.