Quickstart¶
Introduction¶
DESPOTIC is a tool to Derive the Energetics and SPectra of Optically Thick Interstellar Clouds. It can perform a variety of calculations regarding the chemical and thermal state of interstellar clouds, and predict their observable line emission. DESPOTIC treats clouds in a simple one-zone model, and is intended to allow rapid, interactive exploration of parameter space.
In this Quickstart, we will walk through a basic interactive python session using DESPOTIC. This will work equally well from an ipython shell or in an ipython notebook.
The cloud
Class¶
The basic object in DESPOTIC, which provides an interface to most of
its functionality, is the class cloud
. This class stores the basic
properties of an interstellar cloud, and provides methods to perform
calculations on those properties. The first step in most DESPOTIC
sessions is to import this class:
from despotic import cloud
The next step is generally to input data describing a cloud upon which
computations are to be performed. The input data describe the cloud’s
physical properties (density, temperature, etc.), the bulk composition
of the cloud, what emitting species it contains, and the radiation
field around it. While it is possible to enter the data manually, it
is usually easier to read the data from a file, using the
Cloud Files format. For this Quickstart, we’ll use one of
the configuration files that ship with DESPOTIC and that are
included in the cloudfiles
subdirectory of the DESPOTIC
distribution. To create a cloud whose properties are as given in
a particular cloud file, we simply invoke the constructor with
the fileName
optional argument set equal to a string containing
the name of the file to be read:
gmc = cloud(fileName="cloudfiles/MilkyWayGMC.desp", verbose=True)
Note that, if you’re not running DESPOTIC from the directory where you
installed it, you’ll need to include the full path to the cloudfiles
subdirectory in this command. Also note the optional argument
verbose
, which we have set to True
. Most DESPOTIC methods
accept the verbose
argument, which causes them to produce printed
output containing a variety of information. By default DESPOTIC
operations are silent.
Computing Temperatures¶
At this point most of the calculations one could want to do on a cloud
are provided as methods of the cloud
class. One of the most basic is
to set the cloud to its equilibrium dust and gas temperatures. This is
accomplished via the setTempEq
method:
gmc.setTempEq(verbose=True)
With verbose
set to True
, this command will produce variety of
output as it iterates to calculate the equilibrium gas and dust
temperatures, before finally printing True
. This illustrates
another feature of DESPOTIC commands: those that iterate return a
value of True
if they converge, and False
if they do not.
To see the gas and dust temperatures to which the cloud has been set, we can simply print them:
print gmc.Tg
print gmc.Td
This shows that DESPOTIC has calculated an equilibrium gas temperature of 10.2 K, and an equilibrium dust temperature of 14.4 K.
Line Emission¶
Next we might wish to compute the CO line emission emerging from the
cloud. We do this with the cloud
method lineLum
:
lines = gmc.lineLum("co")
The argument co
specifies that we are interested in the emission
from the CO molecule. This method returns a list
of dict
, each
of which gives information about one of the CO lines. The dict
contains a variety of fields, but one of them is the
velocity-integrated brightness temperature of the line. Again, we can
just print the values we want. The first element in the list is the
\(J = 1 \rightarrow 0\) line, and the velocity-integrated
brightness temperature is listed as intTB
in the dict
. Thus to
get the velocity-integrated brightness temperature of the first line,
we just do:
print lines[0][’intTB’]
This shows that the velocity-integrated brightness temperature of the CO \(J = 1 \rightarrow 0\) line is 79 K km/s.
Heating and Cooling Rates¶
Finally, we might wish to know the heating and cooling rates produced
by various processes, which lets us determined what sets the thermal
balance in the cloud. This may be computed using the method dEdt
,
as follows:
rates = gmc.dEdt()
This method returns a dict
that contains all the heating and
cooling terms for gas and dust. For example, we can print the rates of
cosmic ray heating and CO line cooling via:
print rates["GammaCR"]
print rates["LambdaLine"]["co"]