Free Access
Issue
A&A
Volume 642, October 2020
Article Number A127
Number of page(s) 12
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201935610
Published online 13 October 2020

© ESO 2020

1 Introduction

Circumbinary disks make up a significant fraction of protoplanetary disks (Duquennoy & Mayor 1991; Kraus & Hillenbrand 2009). The additional torque provided by the binary is expected to alter various physical processes in protoplanetary disks, such as planetesimal growth (Marzari et al. 2013) and planetary migration behavior (Kley & Haghighipour 2014). Inner regions are expected to show characteristic structures resulting from tidal interactions between the time-dependent non-axisymmetric gravitational potential generated by the binary and the surrounding disk. Typical structures are cavities in the disk center, accretion arms extending into these cavities, and density waves at the inner rim of the disk. All of these features are potentially most pronounced in the innermost regions where the gravitational potential of the stellar binary shows the strongest deviation to that of a single star. However, not all of these structures are unique. For example, inner cavities are also present in transitional disks (Garufi et al. 2013; Dong et al. 2017), are caused by photoevaporation (Clarke et al. 2001), or are due to the presence of a Jupiter-mass planet (van der Marel et al. 2016a,b). In contrast to photoevaporation, the tidal force generated by a binary is expected to affect the dust and gas phase equally. Moreover, due to the stronger tidal forces, the impact is more pronounced than in the case of a planet.

In the case of a disk seen in face-on orientation, the structures resulting from binary-circumbinary disk interactions are expected to be detectable using proper combinations of infrared and sub-millimeter wavelength continuum observations (Avramenko et al. 2017; Ruge et al. 2015). However, disks inclined edge-on obscure the innermost disk regions at optical to infrared wavelengths. Even at submm/mm wavelengths, the disk may still be optically thick in the continuum, as can be seen in the case of the Butterfly Star (Gräfe et al. 2013), where a flux reduction for submm wavelengths is observed in the disk’s central region. This reduction is attributed to the dust grain growth and vertical settling of larger dust grains. In contrast to continuum observations, molecular line emission has the advantage of the Doppler shift in frequency, which makes the optical thickness dependent on the velocity distribution of gas in the disk. By choosing proper molecular lines (i.e., isotopologs and transitions), it is therefore possible to directly observe the inner disk regions. Such observations were conducted by Dutrey et al. (2017), who detected a gap in the gas density of the “Flying Saucer” disk. A prominent example of an edge-on circumbinary disk cavity in the center is HH30. The cavity is supposed to be caused by a binary with 0.31 M and 0.14 M and a separation of 18 au (Lynch et al. 2019). Louvet et al. (2018) proposed an alternative explanation of the cavity based on the modeling of dust continuum and molecular transition line Atacama Large Millimeter/submillimeter Array (ALMA) observations. According to Louvet et al. (2018), the cavity could be the result of a bow shock propagating outward.

The goal of our study is to evaluate the feasibility of identifying characteristic binary-induced structures in edge-on circumbinary disks using spatially resolved or unresolved molecular line observations. To address this goal, we conducted 2D hydrodynamic simulations of circumbinary disks with the code Fosite (Illenseer & Duschl 2009). The resulting density profiles were used to calculate the temperature of the system with the 3D grid-based Monte Carlo code Mol3D (Ober et al. 2015). This code was further applied to calculate the level populations for individual molecular lines and subsequently combined with the velocity fields to generate molecular line emission maps. Based on these maps, we evaluated the feasibility of detecting the cavity, accretion arms, and density wave.

The paper is divided into two parts. In the first part (Sect. 2), we give a short overview of the applied density distributions and radiative transfer simulations. Next, we discuss the simulated observations of the resulting maps in spatially resolved and unresolved cases (Sect. 3). The prime focus lies on the possibility of reconstructing characteristic structures for the circumbinary disks.

Table 1

Parameters of the hydrodynamic simulation.

Table 2

Parameters of the radiative transfer simulation.

2 Density distributions and radiative transfer simulations

In this section, we give a brief overview of the software used to conduct the hydrodynamic and radiative transfer simulations. A complete discussion can be found in Avramenko et al. (2017).

2.1 Density distributions

The density and gas velocity distributions were calculated using the 2D grid-based hydrodynamic code Fosite (Illenseer & Duschl 2009). A full discussion of the simulation setup, as well as the parameter space, is available in Avramenko et al. (2017), where we investigated the feasibility of observing binary-induced disk structures in the continuum. For the current study, we used a single hydrodynamic simulation of this sample, namely a disk with a semimajor axis of a = 30 au and binary masses of Mprim = Msec = 1 M, with Mprim and Msec as the masses of the primary and secondary components, respectively. As was found by Avramenko et al. (2017), these parameters result in particularly large accretion arms and density oscillations. Both factors increase the chances of detecting these features.

Because Fosite only calculates the 2D surface density distribution, we reconstructed the vertical density profile of the disk using the geometrical thin disk model (Lynden-Bell 1969; Shakura & Sunyaev 1973). The scale height of the disk was calculated during the hydrodynamic simulation under the assumption of a vertically isothermal disk model. To get the appropriate scale height, we iterated between (a) calculating the scale height with Fosite using a given temperature distribution and (b) subsequently calculating the temperature with Mol3D (Ober et al. 2015) for the derived vertical disk structure. This iteration was repeated until the change in the resulting scale height between two subsequent iterations was below 1 percent. The initial temperature distribution was set to 50 K (isothermal).

The vertical component of the velocity cannot be deduced from the simulations, and we assume that it is small compared to the orbital velocity (i.e., primarily governed by Keplerian motion). An indirect contribution was made, however, by including a general line broadening due to turbulence in the disk. We will further discuss this effect in the next section. The full parameter set is compiled in Table 1.

In preparation of the analysis of the circumbinary disk model, we also considered an analytically described density distribution of a Keplerian circumstellar disk (see Eq. (1)), with the quantity ρ being the volume density, Σ the surface density, H the pressure scale height,and z the vertical disk coordinate. The central star in the Keplerian case has the combined mass of the binary (2 M) of the circumbinary case. First studying the Keplerian disk, as well as selected modified versions (e.g., the inner cavity and gap), will help us to identify distinct features that are associated with a specific morphology of the disk: ϱ=Σ2πHexp[12(zH)2].\begin{equation*}\varrho\,{=}\,\frac{\Sigma}{\sqrt{2 \pi}H} \, \exp\left[-\frac{1}{2}\left(\frac{z}{H}\right)^2\right] .\end{equation*}(1)

thumbnail Fig. 1

Midplane temperature distribution for a circumbinary disk.

thumbnail Fig. 2

CO J = 4− 3 Velocity-resolved surface brightness distribution of a Keplerian disk seen edge-on (see Sect. 3.1.1 for details). (a) No features. (b) Cavity at rin = 60 au. The inclined black dashed line indicates the lack of flux inside the cavity. (c) Gap between rin = 40 and 60 au. The inclined black dashed lines indicate the lack of flux inside the gap. (d) Sin-wave density profile. Black lines indicate the flux maxima (minima) corresponding to the location of maxima (minima) of the density distribution.

2.2 Radiative transfer simulation

The density distribution calculated with Fosite was further used to calculate the corresponding temperature distribution. This was done with the Monte Carlo-based 3D continuum and line radiative transfer code Mol3D (Ober et al. 2015). As above, we refer to Avramenko et al. (2017) for a full discussion of the radiative transfer simulations. The parameters used in the current study are summarized in Table 2.

The synthetic molecular line maps for the individual velocity channels were calculated using the large velocity gradient (LVG) method (Pavlyuchenkov et al. 2007). For the line broadening, we assumed a model based on earlier line observations by, for example, Piétu et al. (2007) and Chapillon et al. (2012) according to the formula V=2kTm+vturb2,\begin{equation*}\bigtriangleup V\,{=}\,\sqrt{ \frac{2kT}{m} + v_{\textrm{turb}}^2}, \end{equation*}(2)

with △V being the total line broadening, k the Boltzmann constant, T the temperature, m the molecular mass, and vturb the turbulent contribution. Piétu et al. (2007) assumed vturb to be a power law with regard to the radius. To reduce the number of free parameters, we instead used a constant value of 100 m s−1.

For our study, we used two different isotopologs of the carbon monoxide molecule: CO and C18O. For the abundances we assumed molecule-to-hydrogen ratios (N/N(H)) of 10−4 for CO and 2 × 10−7 for C18O, as reported in Pavlyuchenkov et al. (2007). These were calculated under the assumption of the gas-grain chemical model with surface reactions from Semenov et al. (2004, 2005). The temperature in the disk midplane is above 20 K (Fig. 1), which should prevent CO from freezing out.

Transition frequencies and transition rates for the considered isotopologs are taken from the Leiden Atomic and Molecular Database (LAMDA; Schöier et al. 2005). The first three rotational transition lines for each of these isotopologs, J = 2−1, J = 3−2, and J = 4−3, were considered. We used a maximum velocity of 6.5 km s−1 for the Keplerian circumstellar disk and 7.0 km s−1 for the circumbinary disk, with 67 and 99 velocity channels, respectively. This results in a spectral resolution that is comparable to the maximum resolution achievable with ALMA in the corresponding bands.

thumbnail Fig. 3

Velocity-resolved surface brightness distribution of a circumbinary disk (see Sect. 3.1.2 for details). (a) CO J = 2− 1. (b) CO J = 4− 3. (c) C18O J = 2− 1. (d) C18O J = 4− 3.

3 Results

This section is divided in two parts. First we analyze simulated “spatially resolved” molecular line maps. We begin with a standard Keplerian disk as a reference case of a featureless disk. After that we discuss the case of a Keplerian disk with a large inner cavity and another with a gap in the density distribution located between 40 and 50 au. We proceed with a disk showing a radial wave structure in its density. Finally, we discuss the case of a circumbinary disk.

The second part deals with “spatially unresolved” observations. There we follow the same sequence in the discussion of our results (i.e., starting with the circumstellar case). Afterward, we connect the distinctive elements in our data with the density gradient expected to occur at the inner edge of a circumbinary disk.

The midplane of the disk coincides with the x-y plane. We are discussing edge-on disks with the observer’s line of sight along the y axis. All models are presumed to be located at a distance of 140 pc.

3.1 Spatially resolved line observations

3.1.1 Keplerian disk

Our simulation results are stored in a 3D data cube with two spatial dimensions (along the x and z axes) for each of the velocity channels. In this study, we are interested in the spatial distribution of gas in the disk. To achieve better contrast and reduce the number of necessary plots, we integrated the flux density along the vertical extent of the edge-on viewed disk and divide it by the number of contributing pixels (see Fig. 2a for an example of an unmodified edge-on Keplerian disk). Due to Keplerian motion, each part of the disk moves with a specific projected velocity to the line of sight of the observer, resulting in a frequency shift. The projected velocity depends on the position along the x axis. The values of the flux along a straight line through the origin correspond to the radiation of the gas at a specific radius in the disk. To determine the radius, one notes the x value in a velocity channel at which the flux falls to zero (for the derivation, see Appendix A in Dutrey et al. 2017). The contributions from different radii are given along lines with different slopes in the diagram.

This interpretation of the velocity-position diagram is illustrated by the case of a Keplerian disk with a large cavity from which the disk’s inner radius of 60 au and the corresponding Keplerian velocity of ~ 3.6 km s −1 can be derived (Fig. 2b). As both sides of the edge-on disk are observed, the extent of the cavity in each direction can be determined separately (important for disks with non-circular cavities).

In the case of the Keplerian disk with a radial gap, the inner and outer radii of this gap can be identified (see Fig. 2c). Here, the location and extent of the gap can be directly deduced from the velocity-position diagram as well. A similar observation was already made by Dutrey et al. (2017), resulting in the identification of a gap in the Flying Saucer disk using CO J = 2−1 observations.

In the final example, the Keplerian disk density distribution is modified with a sin-function to mimic a density wave propagating radially outward. The resulting velocity-position diagram (Fig. 2d) is similar to that generated by a superposition of “gap-like” features with varying radii. In this case, it would be possible to determine the wavelength of the density oscillation by measuring the distance between the flux minima (or maxima) in a similar manner, as in the case of the analysis of a gap structure.

3.1.2 Circumbinary disk

Now we focus our attention on the case of a circumbinary disk. Figure 3 show the vertically integrated fluxes for transition lines J = 2−1 and J = 4−3 of CO and C18O. As expected, we observe low fluxes in the region close to the central star (i.e., small radii) and thus high velocities. This is a clear indication of the cavity associated with the circumbinary disk. However, as can be best seen in the case of the CO J = 4−3 line, the cavity isnot completely empty. One can clearly see the fiber-like structure that can be attributed to the presence of accretion arms in the disk. In fact, if we extract the inner radius of the disk from the fluxes of each CO isotopolog and each line, we derive the inner radii shown in Fig. 4. Here we can clearly see that the more abundant isotopolog CO leads to a smaller derived radius. We determine a radius of 60 and 50 au for J = 2−1 and J = 4−3, respectively. For C18O, we derive inner radii of 72 and 70 au for transition lines J = 2−1 and J = 4−3, respectively.

The isotopolog C18O also allows us to constrain the inner disk structure (see Figs. 3c and d). It shows a harp-like structure, which is caused by the density waves propagating through the disk, similar to the case discussed in Sect. 3.1.1 and Fig. 2d. If the spectral resolution is high enough, this allows one to detect the flux differences between the corresponding density maxima and minima. Moreover, since we observe multiple maxima and minima, it is possible to derive the wavelength of the oscillation. As was shown by Avramenko et al. (2017), the wavelength of the density wave is connected to the basic parameters of the binary system, such as binary mass and binary separation.

thumbnail Fig. 4

Derived inner radii plotted onto the initial surface density distribution. CO : J = 2− 1 (black, dash-dotted); CO: J = 4− 3 (black, dashed); C18O: J = 2− 1 (white, dotted); C18O: J = 4− 3 (white, dashed).

thumbnail Fig. 5

Keplerian disk line CO, J = 4− 3 emission. Top: spatially resolved images for specific channels numbered 1 through 5. Bottom: CO line profile. The channels from the top image are highlighted in blue and denoted by a number.

thumbnail Fig. 6

Keplerian disk line profile (see Sect. 3.2.1 for details). (a) CO, J = 4− 3. (b) C18O, J = 4− 3. (c) Keplerian disk with cavity. CO J = 4− 3 line emission. For details, see Sect. 3.2. (d) Keplerian disk with cavity. C18O J = 4− 3 line emission. For details, see Sect. 3.2.

3.2 Spatially unresolved line observations

In the previous section, we discussed the potential of spatially resolved line observations to constrain selected features of the disk morphology, such as the presence and shape of a central cavity. However, the angular resolution and/or brightness sensitivity of the observations may be insufficient to study this morphology. For this reason, we investigate the feasibility of deriving selected disk parameters using spatially unresolved observations. To achieve this, we calculated the ratios between two spatially integrated molecular lines. In one case, we combined data from the same isotopolog but for two different transition lines. In the second case, we used the same transition line J but considered different isotopologs. This allows one to trace regions where the ratio shows a high degree of variation, which is expected to occur in areas with high temperature and density gradients inside the cavity. We will first briefly discuss the case of a molecular line in an edge-on disk. After that we will address the line ratios for the simple case of a Keplerian disk. This will allow us to better identify the features associated with the presence of a gap in the disk, followed by the case of a circumbinary disk.

In order to understand the shape of a molecular line, one has to understand the shape of the area of the disk in each velocity channel. A more detailed discussion of line formation in protoplanetary disks is available in Sect. 2.2 of Beckwith & Sargent (1993). The exact shape of the covered area can be calculated with r(ϕ)=GMstarvobssin2θcos2ϕ,\begin{equation*}r(\phi)\,{=}\,\frac{GM_{\textrm{star}}}{v_{\textrm{obs}}} \sin^2{\theta}\cos^2{\phi} \,, \end{equation*}(3)

with G being the gravitational constant, Mstar the mass of the central object, vobs the average projected velocity toward the observer, and ϕ and θ the angles in polar coordinates. Figure 5 shows the generic CO J = 4−3 line of a Keplerian disk seen edge-on. Additionally, the spatially resolved images of the edge-on disk are displayed for the highlighted velocity channels. The line profile shows two distinct parts. The first part corresponds to the area covered by a singular velocity channel that moves outward in the disk (Numbers 1, 2, and 3 in Fig. 5). After reaching the peak flux, the area contributing to the flux starts moving toward the center of the disk (Numbers 4 and 5 in Fig. 5). This behavior can be explained by the projection of the Keplerian velocity on the line of sight of the observer. The process repeats itself on the other side of the disk. In Fig. 1 of Beckwith & Sargent (1993), one can see the exact shape of the area covered by a velocity channel in connection with the location at the line profile. For the Keplerian disk, both sides are symmetric. This results in the line profiles as well as the line ratios being symmetric as well. Later we will see that the two sides of the line spectra are not exactly symmetric in the case of a binary disk. For that reason, we restrict figures of Keplerian disks to one side of the spectrum. For binary cases, both sides are shown.

3.2.1 Keplerian disk

In Fig. 5, we highlight the link between the spatial shape and position of individual velocity channels and the lineprofiles. However, the shapes of channels can change depending on the optical depth of the line. The optical depth changes depending on the abundance of the isotopolog and the observed transition line. Figures 6a and b show two examples of line profiles for the CO (Fig. 6a) and C18O (Fig. 6b) transition line J = 4−3 for a Keplerian disk. Although the peaks are located at the same velocities, the line profiles at those peaks differ significantly from each other. This can be attributed to the difference in optical depth between both lines. In our models, we assumed the abundance ratio of N(CO)∕N(C18O) = 500. The peaks of the line correspond to the velocity channel reaching the outer edge of the disk (Number 3 in Fig. 5). From this we can infer that the shape of the area in which the optical depth τ is greater than unity differs significantly between isotopologs for the those channels. Figure 7 shows the optical depth of the line C18O J = 4−3. The black contour indicates the optical depth τ = 1 seen from the observer (y-axis). In Fig. 8, one can see the difference in the location of the τ = 1 surface between different isotopologs and different lines. At this point we want to note that although the line C18O is optically thinner compared to CO, it still reaches an optical depth of τ > 1 in some parts of the disk, as can be seen in Figs. 7 and 8.

Following the discussion of the individual line profiles, we now discuss the line flux ratios. Figure 9a shows an example of the flux ratio for CO and C18O for the line J = 4−3. For the Keplerian disk, the flux ratio has only two features: at 0.3 and 2.8 km s−1, respectively. Comparing the line profiles, one finds that the feature at 2.8 km s−1 is caused by the difference in radiating areas between the isotopologs at the maximum extent of the velocity channel (Number 3 in Fig. 5). For this reason, it can be used to determine the velocity at the disk’s outer edge. However, for a disk with a less sharp edge, we expect this effect to be less pronounced. The second feature at 0.3 km s−1 is caused by the velocity channel obscuring itself (Number 5 in Fig. 5). Once again, the difference in abundance results in a difference in size of the τ ≥ 1 area. This, in turn, results in the self-obscuring occurring at slightly different velocities.

For the C18O J = 4−3 to J = 3−2 line ratio (Fig. 9b), the features occur at the same velocity as in Fig. 9a. The effect is smaller in both cases due tosmaller differences in optical depth between the two lines (see Fig. 8). In the previous case, this difference was the result of the difference in abundance between the two isotopologs. In this case, the difference in the optical depth is caused by the temperature gradient in the disk and the associated difference in excitation of the individual levels. Those features are a geometrical effect of the velocity channels distributed in the disk and the way we observe the line emissions. They will appear in some form in all graphs presented from now on and will be subsequently ignored in the main text.

thumbnail Fig. 7

Optical depth for the C18O J = 3− 2 line in the Keplerian disk with cavity at ~4.0 km s −1 (i.e., a radius of ~60 AU). The optical depth is measured (i.e., integrated) in the y direction. The black contour line denotes the τ = 1 boundary.

thumbnail Fig. 8

Normalized density distribution. Plotted are optical depth τ = 1 boundaries for transition lines CO J = 4− 3 (green), C18O J = 4− 3 (blue), and C18O J = 3− 2 (black).

thumbnail Fig. 9

Keplerian disk: selected line emission ratio (see Sect. 3.2.1 for details). (a) CO /C18O, J = 4− 3. The highlighted velocity channels at 0.3 and 2.8 km s−1 indicate thefeatures caused by self obscuring and differences in radiating areas between the isotopologs, respectively. (b) C18O, J = 4− 3/J = 3−2. The highlighted velocity channels at 0.3 and 2.8 km s−1 indicate thefeatures caused by self-obscuring and differences in radiating areas between the isotopologs, respectively. (c) CO /C18O J = 4− 3. The blue highlighted channel at 4.2 km s−1 indicates thevelocity at the cavity outer edge. (d) C18O J = 4− 3/J = 3−2 line emissionratio. The blue highlighted channel at 4.2 km s−1 indicates the velocity at the cavity’s outer edge.

3.2.2 Keplerian disk with cavity

We now discuss the case of a Keplerian disk with an inner cavity (R = 60 au), which is comparable in size to the depleted region in our circumbinary disk. Similar to the previous subsection, we will address the individual line profiles first before continuing with the flux ratios. Figures 6c and d show line profiles for line J = 4−3 for isotopologs CO and C18O. The lack of flux at higher velocities in the center of the disk is indicative of the presence of the cavity. However, the velocity at which this occurs differs for both isotopologs. Once again, we can attribute this difference to the difference in optical depth. Due to a higher abundance of CO compared to C^18O, the flux level increase is stronger for CO. This results in the maximum flux being reached at a higher velocity, corresponding to a smaller radius. Comparing this to the flux ratio shown in Fig. 9c, one finds that this difference results in a clearly defined maximum in the flux ratio diagram.

In the case of the C18O J = 4−3 to J = 3−2 (Fig. 9d) ratio, the maximum is much less pronounced. As in the previous case, the feature is caused by the difference in optical depth between the two lines and the relative level population of both levels. The lines have different excitation energies. This results in the J = 4−3 line having a higher excitation ratio compared to the J = 3−2 line at the warmer, directly illuminated inner edge of the disk. The excitation ratio changes further inside the disk. Here the temperature is lower compared to the disk’s outer layer. This change manifests itself in the feature at ~ 4.2 km s−1 in Fig. 9d.

thumbnail Fig. 10

Circumbinary disk line CO, J = 4− 3 flux (see Sect. 3.2.3 for details).

3.2.3 Circumbinary disk

We begin our discussion by once again introducing an unaltered line profile. The emission for CO J = 2−1 for a circumbinary disk is shown in Fig. 10. Similar to the case of a Keplerian disk with a cavity (Fig. 6c), we detect a flux deficiency for high velocities compared to the Keplerian disk (Fig. 6a). However, since the cavity in the circumbinary disk is not completely gas-free, we still detect a significant flux for velocities above 4 km s−1. It is possible to determine the size of the inner cavity from a simple line profile, as was done by Dutrey et al. (2008). To do so, one has to perform an extensive fitting of multiple parameters simultaneously. With the flux ratios, one has an additional way of determining the cavity size that does not require fitting. This can be used to further constrain the parameter range, since the resulting cavity size does not depend on the other fitting parameters.

For the circumbinary disk, we calculated all six flux ratios that can be derived from our synthetic data. In particular we have three different ratios of CO /C18O for each line transition and three ratios for each isotopolog: J = 4−3∕3−2, J = 4−3∕2−1, and J = 3−2∕2−1.

Comparing the flux ratio for different isotopologs in Fig. 11a to the example of a Keplerian disk with a cavity (Fig. 9a), one notices a number of major differences. The first is the asymmetry between the right and the left side of the graph. This is caused by the asymmetry in the density of the hydrodynamical distribution. This allows for an independent investigation of both disk sides.

We now discuss the flux ratio between different transition lines for the same isotopolog. The feature corresponding to the inner disk rim appears at a higher velocity of vmax = 6.1 km s−1 in the case of CO (Fig. 11c) compared to vmax = 4.2 km s−1 for C18O (Fig. 11d); this is in agreement with the result we derived from the simulated spatially resolved observations (see Sect. 3.1.2 and Fig. 4 for comparison).

3.2.4 Observability

Spatially unresolved observations do not have the constraint of minimum-required spatial resolution, and they reduce the required sensitivity for the method presented. However, we still have to evaluate whether the sensitivity of the current instruments, ALMA in particular, is high enough to observe the studied features, in this case the cavity in the center of the disk. The peak fluxes in the synthetic images presented range between 0.435 Jy for C18O line J = 2−1 and 40.9 Jy for CO line J = 4−3. However, weare not interested in the regions that contribute to the line maxima, but rather the wings of the line. At these wings, the flux reaches values of, at most, an order of magnitude below the line maximum.

Since we are not interested in a high spatial resolution, we chose the ALMA antenna configuration with the smallest baseline possible (i.e., configuration Number 1 from Observational Cycle 7). This results in a resolution of 3.38 as, 1.47 as, and 0.98 as for lines J = 1−2, J = 2−3, and J = 3−4, respectively.The difference in transition frequency between both isotopologs is small enough for the resolution to be considered the same in both cases. At an assumed object distance of 140 pc and disk diameter of 400 AU, the disk has an angular extent of 2.85 as. However, at a spectral resolution of 141 m s−1, the area contributing to each individual channel is much smaller. This makes it possible to cover the contributing regions with a single beam for all three lines. As was the case with the synthetic lines in the previous subsection, we integrated the flux over the spatial dimensions for each velocity channel. We simulated an ALMA observation of 2 h using the ALMA simulation tool kit CASA (McMullin et al. 2007) for three different lines for both isotopologs. For this, we assumed the default value for the precipitable water vapor of 0.5 mm and an ambient temperature of 269.0 K. By comparing the integrated flux values before and after applying CASA, we can assure the flux conservation up to the noise level. Figure 12 shows an example of a line ratio observed with ALMA. In this image, one can discern the general form of a flux ratio before it is processed with CASA (compare to Fig. 11a). The ratios coincide well in channels with higher fluxes. Since it is not possible to determine the location of the ratio maxima from Fig. 12, it is not possible to determine the radius of the cavity. However, one can clearly see the increase in flux ratios around ±4.5 km s−1. With this we can can still give an upper boundary for the cavity radius.

While we assumed typical disk properties in our study, the presented method becomes much more powerful if applied to a disk with higher apparent brightness, such as the Flying Saucer in Dutrey et al. (2017). Their observations were performed under conditions similar to those described here. The two major differences were a significantly higher maximal surface brightness of the Flying Saucer of about a factor of ten and a larger CO disk extent of 5.5 (compared to the 2.85 used here).

3.2.5 Spatial distribution

So far we have only considered the correlation between the location of the disk’s inner edge and the velocity v. Beyond this it is possible to match each velocity channel to a specific volume inside the disk. For this we need the velocity field of the circumbinary disk. Assuming that it is mostly governed by the balance between the gravitational and centrifugal forces, it is sufficient to know the binary masses M1 and M2 and their positions (in the x-y plane according to our reference coordinate system). It is possible to use the reduced mass of the binary μ = (M1M2)∕(M1 + M2) and presume thevelocity field to be near-Keplerian. The reduced mass can be estimated from the radius of the outer edge of the disk and the corresponding gas velocity. The gas velocity can be derived with the associated feature in the line ratio discussed in Sect. 3.2 or by measuring the total line width for each line. However, as the binary and single star velocity profiles differ, this leads to errors in determining the spatial location of the ratio features. In Fig. 13, the absolute value of the difference between the gravitational potential of the binary and a single star with the combined mass divided by the single star potential is shown. A contour plot of the midplane gas density is superimposed over the relative difference plot. The white X symbols mark the binary positions. At the inner rim of the disk, the region we are most interested in, the difference in the potentials still amounts to 20 percent and decreases further in the disk.

Figure 14 shows a reconstruction of the inner disk structure using the line flux ratio. For this we used the binary mass and binary position data from our simulations. The binary positions are once again marked with white X symbols. The binary positions were used to calculate the velocity field of the system. The velocity field was binned in the same manner as the velocity channels. Subsequently, we assigned the area of the disk associated with a specific channel the line ratio derived in that channel. The midplane density distribution was overlaid as a contour plot. In Fig. 14, one finds that the outer edge of the feature seen in ratio maps, corresponding to higher orbital velocity, coincides with the inner rim of the disk. Furthermore, one finds that the inner edge of the feature coincides with the second maximum of the density oscillation. The inner edge of the feature is the point at which both lines become optically thick.

thumbnail Fig. 11

Circumbinary disk: selected line emission ratio (see Sect. 3.2.3 for details). (a) CO to C18O line J = 4−3. The blue highlighted channel at −5.1 and 5.1 km s−1 indicates thevelocity at the cavity’s outer edge. (b) CO to C18O line J = 2−1. The blue highlighted channel at −4.7 and 4.8 km s−1 indicates thevelocity at the cavity’s outer edge. (c) CO line J = 3−2∕2−1. The blue highlighted channel at −6.1 and 6.1 km s−1 indicates thevelocity at the cavity’s outer edge. (d) C18O line J = 3−2∕2−1. The blue highlighted channel at −4.2 and 4.2 km s−1 indicates the velocity at the cavity’s outer edge.

3.2.6 Model scaling

In this subsection, we discuss the influence of some of the more important model parameters on the studied features. We then, in turn, discuss the feasibility of deriving them from the observations.

The two main features discussed here are the inner cavity and the density waves. The size of the cavity is determined by the location of the last outer Lindblad resonance. In the case of a binary of equal mass, it is only dependent on the binary separation a (see Artymowicz & Lubow 1994; Eq. (2.2)). A larger binary semimajor axis a would result in a larger cavity. This, in turn, would make it easier to detect. The density waves propagate outward from the disk’s inner edge with the sound speed cs as the phase velocity. The wavelength of the oscillation is therefore dependent on the binary period P (see Sect. 3.1 of Avramenko et al. 2017). The binary period can be influenced by changing the mass as well as by changing the semimajor axis. Thus, decreasing the binary separation would result in smaller and therefore less detectable density waves and cavities. Additionally, a smaller cavity would lead to higher temperatures at the disk’s inner edge. For that reason, one would have to turn to molecular lines with higher excitation energies. Changing the binary mass, and therefore the temperature, would, for once, have a similar effect: here we would expect a modified temperature at the inner rim. On the other side, though it does not change the cavity size, the binary mass has a significant impact on the wavelength of the density waves. A lower binary mass results in an increased wavelength of the density wave, making it easier to detect this phenomenon.

Regarding the disk, the two parameters of interest are the disk mass and the abundance of the used isotopologs. We do not expect the disk mass to have a major impact on the applicability of the method. In the case of a spatially resolved observation, one detects the density gradient while the spatially unresolved case locates the optical depth τ = 1 surface for different lines. As long as the disk mass remains large enough for the disk to remain optically thick, it is still possible to locate the density gradient. The isotopolog abundances have the potential to alter the results, since the location of the τ = 1 surface would be changed for each molecule individually. This would change the shape of the line ratio graphs significantly.

For our study, we assumed the binary system to have a distance of 140 pc, motivated by the average distance to the star forming regions in Taurus. As was shown in the previous section, the sensitivity of ALMA is not sufficient for the synthetic observations considered here. Decreasing the distance would result in higher flux levels for all lines (geometrical thinning). This, in turn, would result in an improved observability of the considered features in spatially resolved and unresolved cases.

thumbnail Fig. 12

Synthetic observation of a circumbinary disk CO∕C18O line J = 4−3 emission ratio. Configuration 1 of the ALMA Observational Cycle 7 is used with the integration time of 2 h and assumed distance of 140 pc. The black error bars indicate a 1σ uncertainty.

thumbnail Fig. 13

Absolute value of the difference between the gravitational potential of the binary and a single star with the combined mass divided by the single star potential.

thumbnail Fig. 14

Line emission ratio for CO and C18O J = 3− 2. The contour plot shows the surface density profile of the circumbinary disk. The white X symbols denote the binary positions.

4 Conclusions

The goal of this study was to investigate the feasibility of constraining disk parameters with line observations of a circumbinary disk with an edge-on orientation. In order to do that, we conducted 2D hydrodynamic simulations. The resulting density distribution as well as the associated velocity field were used as inputs for subsequent radiative transfer simulations. In the first step, the 3D temperature distribution was calculated, which, in turn, was used to calculate the level populations and resulting line emission maps.

The study itself was split into two parts. In the first part we considered the case of spatially resolved observations. We find that, despite the edge-on orientation, the Doppler effect allows us to distinguish all characteristic disk regions expected of the considered circumbinary disk model. In the case of a Keplerian disk with a cavity, the cavity size was determined. For the Keplerian disk with a gap, the gap’s inner and outer radii were derived. Using the same method, the wavelength of the density oscillation present in a disk can be extracted. Those features were subsequently shown to be present in synthetic circumbinary disk observations. In particular, we find clear indicators for the presence of a central cavity, accretion arms, and the density wave structure. The derived cavity sizes vary with the considered isotopolog and transition line. We have attributed it to the density gradient at the disk’s inner edge.

Because of the high sensitivity requirements for a spatially resolved observation, we also investigated the potential of deriving constraints on the circumbinary disk structures from spatially unresolved observations. For this we derived the line flux ratios for different isotopolog-transition lines. In the reference case of a Keplerian disk with a central cavity, we find that it is still possible to derive the orbital velocity at the disk’s inner edge from the flux ratio maps. As was the case for spatially resolved images, this velocity varies with isotopologs and transition lines. Subsequently, we identified the same features in flux ratio maps for the circumbinary disk. With regard to the observability of those features with the current instruments, we find that the sensitivity of ALMA is only sufficient to give an upper boundary for the cavity size. By assuming binary masses and positions, we connected the derived velocity with a specific radius. Furthermore, this allowed us to infer information about the gas distribution around the disk’s inner edge.

Acknowledgements

We thank all the members of the Astrophysics Department Kiel for helpful discussions and remarks, as well as the language corrections. This study was funded by the German Science Foundation (DFG), grant: WO 857/12-1.

References

  1. Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 [NASA ADS] [CrossRef] [Google Scholar]
  2. Avramenko, R., Wolf, S., & Illenseer, T. F. 2017, A&A, 604, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Beckwith, S. V. W., & Sargent, A. I. 1993, ApJ, 402, 280 [NASA ADS] [CrossRef] [Google Scholar]
  4. Chapillon, E., Guilloteau, S., Dutrey, A., Piétu, V., & Guélin, M. 2012, A&A, 537, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Clarke, C. J., Gendrin, A., & Sotomayor, M. 2001, MNRAS, 328, 485 [NASA ADS] [CrossRef] [Google Scholar]
  6. Dong, R., van der Marel, N., Hashimoto, J., et al. 2017, ApJ, 836, 201 [NASA ADS] [CrossRef] [Google Scholar]
  7. Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 [NASA ADS] [Google Scholar]
  8. Dutrey, A., Guilloteau, S., Piétu, V., et al. 2008, A&A, 490, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Dutrey, A., Guilloteau, S., Piétu, V., et al. 2017, A&A, 607, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Garufi, A., Quanz, S. P., Avenhaus, H., et al. 2013, A&A, 560, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Gräfe, C., Wolf, S., Guilloteau, S., et al. 2013, A&A, 553, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Illenseer, T. F., & Duschl, W. J. 2009, Comput. Phys. Commun., 180, 2283 [NASA ADS] [CrossRef] [Google Scholar]
  13. Kley, W., & Haghighipour, N. 2014, A&A, 564, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Kraus, A. L., & Hillenbrand, L. A. 2009, ApJ, 704, 531 [Google Scholar]
  15. Louvet, F., Dougados, C., Cabrit, S., et al. 2018, A&A, 618, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Lynch, C. J. R., Smith, M. D., & Glover, S. C. O. 2019, MNRAS, 491, 3082 [Google Scholar]
  17. Lynden-Bell, D. 1969, Nature, 223, 690 [NASA ADS] [CrossRef] [Google Scholar]
  18. Marzari, F., Thebault, P., Scholl, H., Picogna, G., & Baruteau, C. 2013, A&A, 553, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127 [Google Scholar]
  20. Ober, F., Wolf, S., Uribe, A. L., & Klahr, H. H. 2015, A&A, 579, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Pavlyuchenkov, Y., Semenov, D., Henning, T., et al. 2007, ApJ, 669, 1262 [NASA ADS] [CrossRef] [Google Scholar]
  22. Piétu, V., Dutrey, A., & Guilloteau, S. 2007, A&A, 467, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Ruge, J. P., Wolf, S., Demidova, T., & Grinin, V. 2015, A&A, 579, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Semenov, D., Pavlyuchenkov, Y., Henning, T., Herbst, E., & van Dishoeck, E. 2004, Balt. Astron., 13, 454 [Google Scholar]
  26. Semenov, D., Pavlyuchenkov, Y., Schreyer, K., et al. 2005, ApJ, 621, 853 [Google Scholar]
  27. Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
  28. van der Marel, N., Cazzoletti, P., Pinilla, P., & Garufi, A. 2016a, ApJ, 832, 178 [Google Scholar]
  29. van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2016b, A&A, 585, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Parameters of the hydrodynamic simulation.

Table 2

Parameters of the radiative transfer simulation.

All Figures

thumbnail Fig. 1

Midplane temperature distribution for a circumbinary disk.

In the text
thumbnail Fig. 2

CO J = 4− 3 Velocity-resolved surface brightness distribution of a Keplerian disk seen edge-on (see Sect. 3.1.1 for details). (a) No features. (b) Cavity at rin = 60 au. The inclined black dashed line indicates the lack of flux inside the cavity. (c) Gap between rin = 40 and 60 au. The inclined black dashed lines indicate the lack of flux inside the gap. (d) Sin-wave density profile. Black lines indicate the flux maxima (minima) corresponding to the location of maxima (minima) of the density distribution.

In the text
thumbnail Fig. 3

Velocity-resolved surface brightness distribution of a circumbinary disk (see Sect. 3.1.2 for details). (a) CO J = 2− 1. (b) CO J = 4− 3. (c) C18O J = 2− 1. (d) C18O J = 4− 3.

In the text
thumbnail Fig. 4

Derived inner radii plotted onto the initial surface density distribution. CO : J = 2− 1 (black, dash-dotted); CO: J = 4− 3 (black, dashed); C18O: J = 2− 1 (white, dotted); C18O: J = 4− 3 (white, dashed).

In the text
thumbnail Fig. 5

Keplerian disk line CO, J = 4− 3 emission. Top: spatially resolved images for specific channels numbered 1 through 5. Bottom: CO line profile. The channels from the top image are highlighted in blue and denoted by a number.

In the text
thumbnail Fig. 6

Keplerian disk line profile (see Sect. 3.2.1 for details). (a) CO, J = 4− 3. (b) C18O, J = 4− 3. (c) Keplerian disk with cavity. CO J = 4− 3 line emission. For details, see Sect. 3.2. (d) Keplerian disk with cavity. C18O J = 4− 3 line emission. For details, see Sect. 3.2.

In the text
thumbnail Fig. 7

Optical depth for the C18O J = 3− 2 line in the Keplerian disk with cavity at ~4.0 km s −1 (i.e., a radius of ~60 AU). The optical depth is measured (i.e., integrated) in the y direction. The black contour line denotes the τ = 1 boundary.

In the text
thumbnail Fig. 8

Normalized density distribution. Plotted are optical depth τ = 1 boundaries for transition lines CO J = 4− 3 (green), C18O J = 4− 3 (blue), and C18O J = 3− 2 (black).

In the text
thumbnail Fig. 9

Keplerian disk: selected line emission ratio (see Sect. 3.2.1 for details). (a) CO /C18O, J = 4− 3. The highlighted velocity channels at 0.3 and 2.8 km s−1 indicate thefeatures caused by self obscuring and differences in radiating areas between the isotopologs, respectively. (b) C18O, J = 4− 3/J = 3−2. The highlighted velocity channels at 0.3 and 2.8 km s−1 indicate thefeatures caused by self-obscuring and differences in radiating areas between the isotopologs, respectively. (c) CO /C18O J = 4− 3. The blue highlighted channel at 4.2 km s−1 indicates thevelocity at the cavity outer edge. (d) C18O J = 4− 3/J = 3−2 line emissionratio. The blue highlighted channel at 4.2 km s−1 indicates the velocity at the cavity’s outer edge.

In the text
thumbnail Fig. 10

Circumbinary disk line CO, J = 4− 3 flux (see Sect. 3.2.3 for details).

In the text
thumbnail Fig. 11

Circumbinary disk: selected line emission ratio (see Sect. 3.2.3 for details). (a) CO to C18O line J = 4−3. The blue highlighted channel at −5.1 and 5.1 km s−1 indicates thevelocity at the cavity’s outer edge. (b) CO to C18O line J = 2−1. The blue highlighted channel at −4.7 and 4.8 km s−1 indicates thevelocity at the cavity’s outer edge. (c) CO line J = 3−2∕2−1. The blue highlighted channel at −6.1 and 6.1 km s−1 indicates thevelocity at the cavity’s outer edge. (d) C18O line J = 3−2∕2−1. The blue highlighted channel at −4.2 and 4.2 km s−1 indicates the velocity at the cavity’s outer edge.

In the text
thumbnail Fig. 12

Synthetic observation of a circumbinary disk CO∕C18O line J = 4−3 emission ratio. Configuration 1 of the ALMA Observational Cycle 7 is used with the integration time of 2 h and assumed distance of 140 pc. The black error bars indicate a 1σ uncertainty.

In the text
thumbnail Fig. 13

Absolute value of the difference between the gravitational potential of the binary and a single star with the combined mass divided by the single star potential.

In the text
thumbnail Fig. 14

Line emission ratio for CO and C18O J = 3− 2. The contour plot shows the surface density profile of the circumbinary disk. The white X symbols denote the binary positions.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.