Issue 
A&A
Volume 557, September 2013



Article Number  A132  
Number of page(s)  10  
Section  Interstellar and circumstellar matter  
DOI  https://doi.org/10.1051/00046361/201321600  
Published online  20 September 2013 
Online material
Appendix A: Modeling details
Fig. A.1
Spectral energy distribution (black line) of our best fit model, compared to dereddened photometry (A_{V} = 0.3, Q11) from the literature (Q11, Tilling et al. 2012; Cutri et al. 2011; Ishihara et al. 2010; Egan et al. 2003; Helou & Walker 1988). The stellar photosphere is shown as the thin dashed line. 

Open with DEXTER 
Comparison of observed and model fluxes.
Adopted model.
Appendix A.1: Parametric dust and gas structure
We adopt the overall disk properties derived by Q11 from fitting the broadband spectral energy distribution and spatial extent of their mm data: a disk mass M_{disk} = 0.089 M_{⊙}, with an inner edge R_{in} = 0.6 AU and a critical radius R_{C} = 150 AU. We define the surface density as (e.g. Andrews et al. 2009): (A.1)where we set the surface density power law index, γ, equal to 1, the value for a selfsimilar accretion disk in steadystate (Hartmann et al. 1998). We adopt the gastodust ratio of 154 from Q11 as well as their two component dust grain population (small and large) used to simulate dust settling. Both populations are made of a 60%/40% mix of astronomical silicates and graphite, and follow a power law grain size distribution n(a) ∝ a^{3.5} (Draine 2006) with a minimum grain size a_{min} = 0.005 μm for both populations. The small grain population has a maximum grain size a_{max,small} = 0.25 μm, and the large population has a_{max,large} = 1 mm. The small grain component is set to include 1.6% of the dust mass (f_{small} = 0.016).
The vertical distribution is modeled as a Gaussian with angular scale height (A.2)where h_{C} is the angular scale height at the characteristic radius R_{C} and ψ describes the powerlaw disk flaring. We set ψ = 0.066, for consistency with the modeling of Tilling et al. (2012). The two dust component populations have independent scale heights (h_{C}) of h_{small} and h_{large}, respectively, which we adjust to match the spectral energy distribution (discussed below).
The gas density distribution produced by the physical model of Q11 is approximately Gaussian at low heights and has a long “tail” to large heights (cf. their Fig. 7). In order to approximate this, we have constructed the gas density distribution using a twocomponent model with independently varying scale heights for the main component (h_{main}) and the lower mass, larger scale height tail (h_{tail}). In addition, we vary the distribution of gas mass between these two components, with a term f_{tail} describing the fraction of the total gas mass in the tail.
Appendix A.2: Determining the disk parameters
With this framework in place, we adjust parameters to match the SED and CO line fluxes in the literature for HD 163296. We carried out continuum radiative transfer modeling using the 2D code RADMC (Dullemond & Dominik 2004), which receives as input stellar properties and a dust density structure and outputs the resulting temperature structure. It includes a raytracing code, Raytran, for producing the model spectral energy distribution. We adjusted h_{small} and h_{large} in order to approximate the temperature structure found in Q11 (cf. their Fig. 9), which came from a selfconsistent vertical structure calculation. We made further small adjustments to these parameters in order to match the SED. Our adopted disk model has small and largegrain scale heights of 0.08 and 0.06, respectively, and we show the resulting SED in Fig. A.1. In our further modeling, we assume the gas temperature matches the mean dust temperature calculated in this modeling, which is a reasonable approximation in the dense regions upon which we focus.
After setting the dust parameters, we adjust h_{main}, h_{tail}, and f_{tail} to approximate the gas density structure found in Q11 (their Fig. 9). We finetune these parameters by comparing our modeled fluxes of CO and its isotopologues with the observed fluxes of Q11, focusing on the lowJ rotational lines that are expected to come largely from the cool outer disk regions studied here. We assume a relative CO abundance [CO]/[H] ≈ 10^{4} (Frerking et al. 1982), which drops to zero at hydrogen column
densities less than 2 × 10^{21} cm^{2} due to photodissociation (Visser et al. 2009) and at temperatures less than 19 K due to freezeout (Qi et al. 2011). We also assume a ratio of ^{12}C to ^{13}C of 75, and ratios of ^{16}O to ^{18}O and ^{17}O of 500 and 1750, respectively (Frerking et al. 1982). To generate our model line emission, we input our adopted dust density and temperature structure, along with the model gas density structure, to the line radiative transfer code LIME (Brinch & Hogerheijde 2010), and we achieve a match to the observed line fluxes within ~50%.
The gas “tail” has a scale height h_{tail} = 0.2 and includes 5% of the gas mass (f_{tail}), while the “main” gas component has a scale height h_{main} = 0.1. In Table A.1, we list the integrated fluxes of the CO lines which we use to determine the disk gas structure. Table A.2 lists and categorizes all disk modeling parameters, as well as the range of explored values or the reference from which the fixed value was adopted. In both tables, we repeat the HCO^{+}, H^{13}CO^{+}, and DCO^{+} properties discussed in the main body of the text.
© ESO, 2013
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.