Sunday, April 23, 2017

a little bit more about fractional molecular abundances






















2.4 Practical Issues for Deriving Molecular Abundances
2.4.1 Source Structure
The above formulation contains various assumptions. First of all, the cloud is assumed to be homogeneous in physical and chemical conditions along a line of sight. In reality, physical conditions and molecular abundances vary within a cloud. For instance, H2 density will be lower at the periphery of the cloud. Furthermore, molecular abundances will also vary within a cloud because of different chemical processes in different cloud regions. For example, some molecules are abundant in dense regions; other molecules are abundant in the periphery. Observations trace the emission integrated over these effects along a line of sight.
This situation is further complicated by the limited spatial resolution of tele- scopes. If the source size is smaller than the telescope resolution (beam size), the observed emission is significantly diluted. For instance, we consider a circular source with a diameter of θs with a uniform surface brightness on the celestial sphere. If this source is observed with a Gaussian beam with a full-width at half maximum (FWHM) beam size of θb, the observed intensity becomes weaker by a factor of:
!
1 exp 1:44θ2s : ð2:52Þ θ2b
If the surface brightness distribution is represented by a Gaussian function with an FWHM width of θs, the intensity becomes weaker by a factor of:
θ2s : ð2:53Þ θ 2s þ θ 2b
The decrease in the intensity is called the beam dilution effect. The limited resolution causes an additional problem when we employ the intensity data observed with different beam sizes. The beam size of a telescope depends on the observed frequency (Chap. 11). Hence, we have to be careful when the intensities of two or more lines observed at different frequencies are used to derive the column densities and excitation temperatures. If the source size is comparable to or smaller than the beam size, the correction for the beam dilution effect is very important so as not to introduce systematic error into the interpretation of the data.
Moreover, the source structure can be different for different lines of sight. If the cloud has an internal structure whose size is smaller than the telescope resolution (beam size), the observing emission is “beam averaged.” If the molecular abun- dances vary among small-scale structures, the derived column density is also “beam averaged.”
Fractional abundances are typically derived relative to the H2 molecule for comparison with chemical models. In this case, the observed column density has to be divided by the H2 column density observed with the same beam size. Because we cannot directly observe the H2 molecule in interstellar clouds as it lacks electric dipole-allowed transitions in the radio and infrared regimes, the H2 column density is derived from observations of rare isotopic species lines of CO such as C18O by assuming a relative abundance. However, it is known that CO is depleted onto dust grains in cold and dense parts of clouds, as described in Chap. 6. Hence, the H2 column density derived from the rare isotopic species lines of CO would be underestimated. Furthermore, the critical densities of the CO lines are lower than those of other molecular lines, and hence CO emission tends to trace longer distances along the line of sight than the other molecular lines.
Instead of CO, the millimeter-wave and submillimeter-wave continuum emis- sion of dust grains is often used to derive the H2 column density. In this case, we need to assume the emissivity of dust grains and the gas-to-dust ratio (about 100). In any case, we have to consider that the H2 column density and temperature are again beam averaged. Hence, the derived fractional abundances are also beam-averaged values. If the source is expected to have internal structure unresolved with the telescope beam, we have to be careful in interpreting the abundances.
 


What is "fractional abundance" (X) explicitly when talking about molecules?

It's X/H2.




see Jørgensen et al. (2002).


Thursday, April 20, 2017

Summary of meeting with Ted on 4/20

Today I had Ted read my "Section 3" of my paper draft. Some thoughts and action items that fell out of that:


  • The paper will really have to deal with the ammonia issue head-on. Probably with a dedicated section.
  • We'll have to relate *why* hot cores tell you anything about ices.
    • References: the early HDO detection by Jacq, T
      & Nathan's HDO/H2O paper in Orion. (tsr – actually it appears that this is Justin's paper. )
      D/H in Orion is enhanced, inconsistent with the observed temperature, indicating that it was created in an "earlier phase". That's also a big reason why we think the rich organics are found.
    • There *may* be a gas-phase reshuffle of the evaporated ices, but we expect this won't ruin the overall nitrogen accounting.
  • We tried to see if the DCN/HCN ratio told us anything. It's inconclusive. We couldn't find strong DCN vibrational lines in the hexos spectra. (Lines @ 507ish GHz.)
  • Revision of the central point of this work: Ammonia and HCN are the carriers to follow. HCN is the organic carrier. Ammonia needs to be chemically processed before it could be the key progenitor of N in rocks.
    • Actually – amendment to the above point – Scott Messenger mentions "Organic globules" from the planetary science / meteoritics literature that contain possibly ammonia-formaldehyde polymers, so maybe ammonia can be a direct progenitor after all.
  • So to bring it back – Tom, can we use the following statement as an argument for HCN over NH3?
    • "Nh3 doesn't have the right 15N enrichment to be a progenitor of Earth's nitrogen"
    • Tom: check the Hily-Blant paper for anything pointing to this.
  • NH3 *is* present in the ices, but it's not clear whether it's with the water or the organic material... separate or together. (Tom doesn't fully get this last point.)

Thursday, April 6, 2017

First approaches towards getting silicates out of TSC model parameters

So. I'm interested in getting silicon measurements out of my SED models.

I wrote down two approaches in my notebook, for transcription here. For some context: The SED-generating models I'm looking into are consisting of some parameters like this:



  • Version "zero" is to not even assume a gas-to-dust ratio. Version 0 is "take the total mass of the disk/envelope model inside of a relevant radius, assume that *all* the mass is apportioned among elements along cosmic abundance ratios, and take the Si directly from that."
  • Version 1:
    • Assume a gas/dust ratio* (by mass?) within the envelope/disk model.
    • For the dust component of that mass, assume a composition distribution following some vaguely credible reference (e.g. Poteet 2015), and figure out what fraction of the dust mass is in silicon atoms.