---
author:
- |
    People$^{1,2}$[^1]\
    $^{1}$place1\
    $^{2}$place2
date: 'Accepted XXX. Received YYY; in original form ZZZ'
title: 'X-ray photoevaporation gap in the closest planet-forming disc'
---

\[firstpage\]

protoplanetary discs

Introduction
============

Planets form from the reservoir of dust and gas residing in the protoplanetary discs which surround young stars. The physical and chemical characteristics of this material influence the process of planet formation and are in turn influenced by the energetic radiation from the young stellar object. Indeed, young stars are strong sources of X-ray radiation, with approximately one thousandth of their bolometric luminosity being emitted at energies higher than $100\ \mbox{eV}$. X-rays play an important role in the disc evolution, by providing ionisation over a large range of columns (e.g. Igea & Glassgold 1999; Ercolano & Glassgold 2013), which is important for the chemical evolution of the material as well as for angular momentum transport mechanisms. In particular, ionisation of the disc atmosphere by the “soft X-rays” ($100\ \mbox{eV}<\mbox{eV}<1\ \mbox{keV}$), heats the gas and launches a thermal wind (e.g. Owen et al. 2010), which is able to disperse the disc, via (a) the formation of a gap, followed by (b) a hole which quickly grows in radius, with the outer disc being corroded from the inside-out.

The X-ray photoevaporation model (Ercolano et al. 2008, 2009; Owen et al. 2010, 2011, 2012) is successful in reproducing the observed two-timescale dispersal observed (e.g. Luhman et al. ????; Koepferl et al. ????), and it can reproduce the intensities and profiles of emission lines produced in the wind (Ercolano & Owen, 2010, 2016). However, the direct observation of a disc undergoing gap formation via photoevaporation is still lacking. In this letter we propose that the gap at $1\ \mbox{au}$ in the nearest protoplanetary disc, TW Hya, observed in the $870\ \mu\mbox{m}$ continuum by Andrews et al. (2016) using the long-baseline Atacama Large Millimeter/submillimeter array (ALMA), may indeed be the first directly imaged photoevaporative gap in a solar-type star to date. Surprisingly, this possibility has not yet been properly explored in the recent literature.

This Letter is organised as follows: In Section 2 we summarise the relevant characteristics of the TW Hya system. In Section 3 we explore the possible origin of the unresolved 870$\mu$m emission from the inner 0.5 au of the disc. In Section 4 we show dust and gas density distributions from a numerical model of a 1D viscously evolving disc, subject to X-ray photoevaporation. Our conclusions are summarised in Section 5.

The TW Hya system: what is known about our closest neighbour 
=============================================================

TW Hya is considered a young solar analogue with a central stellar mass of approximately $0.7\ M_\odot$. Being the closest gas-rich protoplanetary disc ($54\ \mbox{pc}$ van Leeuwen 2007), this object has been subject of a number of observational campaigns dating back from the early nineties (e.g. Blondel et al. 1993). Its nearly pole-on inclination ($i \sim 7^o$, Qi et al. 2004) further contributes to making TW Hya one of the reference objects for the study of protoplanetary discs.

Disc structure
--------------

Calvet et al. (2002), first reported evidence for a dust-depleted inner hole, which they placed at a radius of $4\ \mbox{au}$ via detailed radiative transfer modelling of its spectral energy distribution (SED). On the basis of the $10\ \mu\mbox{m}$ silicon emission, Calvet et al (2002) concluded that the inner region cannot be completely devoid of dust, and they estimate that $\sim 0.5$ lunar mass of $\sim 1\ \mu\mbox{m}$ particles must be still present inside the cavity, and interpret this as evidence of a [*developing*]{} gap in the disc. The exact size of the gap of TW Hya has been object of dispute in the literature. Near-infrared (Eisner et al. 2006) and mm (Hughes et al. 2007) interferometry confirmed a size of $4\ \mbox{au}$ for the central hole, while a gap radius of only $0.7\ \mbox{au}$ was obtained by Ratzka et al. (2007) via mid-infrared intereferometric observations. The observations can however be reconciled if instead of an inner hole completely devoid of dust, an additional optically thick ring with a $0.5\ \mbox{au}$ radius is considered (e.g. Akeson et al. 2011). More recently Menu et al (2014) presented a new model, which was able to fit all interferometric observations from near-infrared to wavelengths available at the time. This model required an inner cavity with radius $0.3-0.5\ \mbox{au}$ with a maximum in the surface density occurring at $\sim 2.5\ \mbox{au}$.

While, the presence of a gap or at least of a dust depleted region in the inner disc has been confirmed by interferometric observation, the detection of a non-negligible accretion signature (4$\times 10^{-10}\ M_{\odot}/\mbox{yr}$ to $\sim 10^{-9}\ M_{\odot}/\mbox{yr}$, Muzerolle et al., 2000; Eisner et al. 2010), implies that a significant amount of gas must still be present inside the dust cavity.

Gaps and rings: the ALMA view
-----------------------------

Andrews et al. (2016) presented long-baseline ALMA observations of the $870\ \mu\mbox{m}$ continuum from TW Hya, which are sensitive to millimeter-sized dust grains down to a spatial scale of $1\ \mbox{au}$. The data show a series of concentric dark and bright shallow ring-like structures in the outer disc, while the inner disk is dominated by a dark annulus centred at $1\ \mbox{au}$, separating an bright ring that peaks at $2.4\ \mbox{au}$ and an unresolved continuum source less than $0.5\ \mbox{au}$ in radius and emitting $\sim 1\ \mbox{mJy}$ in flux. The modest brightness contrast of some of the concentric features may be associated with condensation fronts of abundant volatiles (Cuzzi & Zahnle, 2004; Okuzumi et al. 2015), or may be evidence of zonal flows (Johansen et al. 2009).

The dark annulus at $\sim 1\ \mbox{au}$ must have, however, a different origin. The ALMA observations together with the dip in the SED at mid-infrared wavelengths (Calvet et al. 2002) suggest that the gap must be depleted of grains from $\mu\mbox{m}$ to $\mbox{cm}$ size. On that basis, Andrews et al. (2016) argue against the dip being a consequence of changing solid properties at around the H$_2$O snow-line. As an alternative, they suggest that a more detailed modeling of the ALMA data may show whether a young super-Earth may be responsible for dynamically carving the observed gap. Surprisingly, the photoevaporation scenario is not considered at all in this work.

Evidence of X-ray photoevaporation in TW Hya 
---------------------------------------------

The TW Hya central object is known to emit substantial X-ray radiation, with a total X-ray flux ($0.1\ \mbox{keV}<E<10\ \mbox{keV}$) of approximately $2.5\times10^{30}\ \mbox{erg/sec}$, based on an extrapolation to lower and higher energies of the $0.3-5\ \mbox{keV}$ XMM-Newton spectrum (Robrade & Schmitt, 2006). Following Equation 9 of Owen, Clarke & Ercolano (2012), the corresponding mass loss rates due to X-ray photoevaporation is $\sim 2\times10^{-8}\ M_{\odot}/\mbox{yr}$.

Direct evidence of currently ongoing photoevaporation is provided by the detection of a few $\mbox{km/s}$ blueshift in the profile of the \[NeII\] $12.8\ \mu\mbox{m}$ line (Herczeg et al. 2007; Pascucci & Sterzik 2009), which is well fit by both EUV and X-ray photoevapoaration models (Alexander 2009, Ercolano & Owen, 2010). The \[OI\]$6300$ has also been shown to be produced in a slow-moving quasi-neutral X-ray photoevaporative wind (Ercolano & Owen, 2016), and is often observed to be blue-shifted with a low-velocity component of a few $\mbox{km/s}$. In the case of TW Hya the \[OI\]$6300$ was found to be centred on the stello-centric velocity (Pascucci et al. 2011). This is not surprising given TW Hya near face-on inclination and the fact that the \[OI\]$6300$ line is emitted from the very inner regions of the disc wind (Ercolano & Owen 2010, 2016), where the underlying disc is devoid of dust, allowing us to see the redshifted outflowing gas. The $10\ \mbox{km/s}$ FWHM and $\sim10^{-5}\ \mbox{L}_\odot$ line luminosity measured by Pascucci et al. (2011) for \[OI\]$6300$ are also well reproduced by the X-ray photoevaporation models presented by Ercolano & Owen (2016) for accretion and X-ray luminosities and the inclination of TW Hya.

The photoevaporative mass loss rates estimated from the X-ray luminosity of TW Hya are roughly a factor of $10$ lower than the measured accretion rate for this object. This is exactly the time at which the X-ray photoevaporation model of Owen et al. (2010) predicts that a gap will open at $1\ \mbox{au}$ for a $0.7\ \mbox{M}_\odot$ star and X-ray luminosity of $2\times10^{30}\ \mbox{erg/sec}$ (e.g. see their Figure 11). Hence, the dark annulus detected by Andrews et al. (2016), may indeed be the first stage of the dispersal phase of this object.

Origin of the unresolved 870$\mu$m emission from the inner 0.5 AU
=================================================================

Andrews et al. (2016) report $\sim 1\ \mbox{mJy}$ flux at $870\ \mu\mbox{m}$ as an unresolved continuum source from the inner $\sim 0.5\ \mbox{au}$ of the TW Hya disc. Free-free gas emission or dust continuum emission could both contribute to the flux at this wavelength.

Free-free emission
------------------

Pascucci et al. (2012) analysed the long wavelength portion of the SED from TW Hya covering from 870$\mu$m out to 6 cm (data from Wilner et al. 2005 and Andrews et al. 2012) and came to the conclusion that about half of the flux detected at 3.5cm could not be accounted for by dust emission. They attributed the extra 140$\pm$40$\mu$Jy at 3.5$\mu$m to a photoevaporative disc wind. However, emission from a disc wind is expected to originate from a region far more extended than 0.5 au, and thus it cannot explain the unresolved 870$\mu$m emission detected by Andrews et al. (2016). Dense ionised gas, however, exists also in the form of a stellar wind or an ionised bound inner disc atmosphere. One can obtain an upper limit to the contribution from gas free-free emission to the sub-mm flux by assuming a spherical optically thick wind. In this case the SED emission in the radio can be described by a power law with index 0.6 (Panagia & Felli, 1975)[^2]. A 0.93mJy flux at 870$\mu$m would then correspond to $\sim$100$\mu$Ju at 3.5cm, explaining thus most of the 140$\pm$40$\mu$Jy excess flux at 3.5cm reported by Pascucci et al. (2012). This would lead to the conclusion that the compact, unresolved source in Andrews et al. (2016), can be caused by a stellar wind and there is no need to invoke the presence of dust in the inner disk. Conversely, this is an upper limit to the emission of the wind at mm-wavelengths: if the wind is optically thin, the spectral index is also expected to be lower, meaning that its contribution in the (sub-)mm will be lower. It is thus not excluded that a fraction of the compact unresolved emission at 870$\mu$m may be caused by dust.

We have calculated photoionisation models of the draining inner disc of TW Hya with [mocassin]{} (Ercolano et al. 2003, 2005, 2008b) using use the gas and dust density distributions from Ercolano & Owen (2016), but limiting the calculations to the inner 0.5 au. We refer to Ercolano & Owen (2010, 2016) for specific details about the calculation methods. Assuming an X-ray luminosity of 2.5$\times$10$^{30}$erg/sec yields a flux of 1.3mJy at 870$\mu$m for a distance of 51pc, more than sufficient to explain the observed values. This is in agreement with the calculations of Owen, Scaife & Ercolano (???), who find a flux of 100$\mu$Jy at 150GHz for a distance of 140pm, assuming a similar X-ray luminosity ( L$_X$ = 2.5$\times$10$^{30}$erg/sec). Assuming the flux scales as $\nu^6$, this corresponds to 1.2 mJy at a distance of 51pc. The bound atmosphere of the draining inner disc of TW Hya could thus alone explain the $\sim$1mJy flux detection from Andrews et al. (2016).

Dust continuum emission
-----------------------

If dust is responsible for all of the emission, then following equation 5 of Andrews & Williams (2005) for a factor five higher temperature, appropriate for the inner $1\ \mbox{au}$ of the disc (e.g. Moulders & Dominik, 2011), we obtain that about $3$ lunar masses of mm-size grains must still be present in the inner disc. This is in contrast to the general assumption that once photoevaporation has disconnected the outer regions of the disc from the inner ones, the dust in the inner disc will disappear almost immediately due to the fast drift timescales. This general assumption originates from calculations by Alexander & Armitage (2007), who find significant radial drift, with larger grains disappearing from the inner disc on timescales of $10^4$. These estimates, however, are tailored to a significantly different disc and photoevaporation model: The mass accretion rate (and thus surface density) assumed in the Alexander & Armitage (2007) models are about an order of magnitude lower than TW Hya, and the $\alpha$ viscosity parameter about a factor ten higher than what would be needed to explain the longevity and high mass of the TW Hya disc ($\sim0.05\ \mbox{M}_{\odot}$ at $10\ \mbox{Myr}$ Bergin et al., 2013; Calvet et al. 2002 suggest $\alpha \sim 10^{-3}$ for a similar mass estimate of $\sim0.06\ \mbox{M}_{\odot}$ at $10\ \mbox{Myr}$)[^3]. These two facts alone result in two orders of magnitude difference in the drift timescales, as will be shown below.

An order of magnitude estimate for the disappearance of mm grains from the inner disc of TW Hya after the photoevaporative gap has opened at 1 au can be obtained as follows. Ignoring factors of order unity and the contribution due to the gas velocity, for small Stokes’ numbers ($St<<1$ ), the radial drift velocity can be expressed as,

$$v_\mathrm{r} = - St \left(\frac{h}{r} \right)^2 v_\mathrm{k},$$

where $h/r$ is the scale hight and $v_k$ is the Keplerian velocity. The radial drift timescale is thus $$t_\mathrm{drift}=\frac{R}{v_\mathrm{r}} = St^{-1} \left(\frac{h}{r}\right)^{-2} \Omega_\mathrm{k}^{-1}$$

With a tipical $h/r=0.033$ at $1\ \mbox{au}$, for $St=1$ particles the radial drift timescale is thus $10^3\ \mbox{yr}$. The Stokes number of $1\ \mbox{mm}$ particles is (e.g. Zhu 2014): $$St = 1.55 \times 10^{-3} \left(\frac{\Sigma}{100\ \mbox{g}/\mbox{cm}^2}\right)^{-1}$$

where $\Sigma$ is the surface density of the gas. If the accretion rate is $\sim 10^{-9}\ \mbox{M}_\odot/\mbox{yr}$ at the time of gap opening, then the surface density at $1\ \mbox{au}$ can be obtained from the following, $$\dot{M}=3 \pi \nu \Sigma = 3 \pi \alpha (h/r)^2 \Omega R^2 \Sigma$$ and therefore $$\Sigma = \frac{\dot{M}}{3 \pi \alpha (h/r)^2 \Omega R^2} =
  135\ \mbox{g}/\mbox{cm}^2 \alpha^{-1},$$ with $\alpha=10^{-3}$. This yields a timescale of $10^6\ \mbox{yr}$ for the drift of the dust from $1\ \mbox{au}$ onto the star. Note however that this is an *upper* limit because as the inner disc drains its surface density decreases, the Stokes number increases and the radial drift velocity increases.

In brief, dust will remain longer in the inner disc if the viscosity is low – which is probably the case in the long-lived TW Hya disc – and/or if the accretion rate at the time of gap opening is high – again a probable scenario in TW Hya given the strong X-ray flux from the central object which results in high photoevaporation rates. This lends credence to the fact that we might indeed be able to observe the gap opening phase in the TW Hya disc.

Gas and Dust surface density Evolution under viscosity and X-ray photoevaporation
=================================================================================

We have developed a simple 1D viscous evolution model which includes X-ray photoevaporation for a TW Hya-like system.

Our 1D gas surface density evolution code, [Spock]{}, is described in detail in Ercolano & Rosotti (2015). At $t=0$ the disc has a mass of $0.15\ \mbox{M}_\odot$, a radius of $18\ \mbox{au}$ and an accretion rate of $10^{-8}\ \mbox{M}_\odot/\mbox{yr}$. We solve the viscous evolution equation, taking into account the photo-evaporative mass-loss, on a grid of $2500$ points uniformly spaced in $R^{1/2}$ extending from $0.025\ \mbox{au}$ to $2500\ \mbox{au}$.

![Gas surface density as a function of disc radius in steps of $10^6$ yr up to 7.5 Myr, and then every $10^5$ yr.](fig1.pdf){width="47.00000%"}

Figure 1 shows the gas surface density as a function of disc radius at in steps of $10^6$ yr up to 7.5 Myr, and then every $10^5$ yr. In our model the gap opens at $1\ \mbox{au}$ at an age of ??$\mbox{Myr}$. At this point the mass accretion rate onto the star is $\sim 10 ^{-9}\ M_{\odot}/\mbox{yr}$. ???? years later (???% of the global disc lifetime and ????% of the transition phase lifetime) the mass accretion rate has already dropped to ????$\sim 10^{-11}\ \mbox{M}_{\odot}/\mbox{yr}$ and the gap outer radius is out to $???\ \mbox{au}$. This implies that one expects $???\%$ of the global protoplanetary disc population and $???\%$ of transition discs to be observed in this early clearing phase[^4]. From the above it appears thus that the initial dispersal phase of gap opening, while brief, should still be observed in a minority of objects given the current disc statistics.

In order to verify the timescales for dispersal of dust particles of various sizes from the inner disc once the gap is opened, we have used the gas surface density distributions obtained with Spock as an input to the <span style="font-variant:small-caps;">fargo</span> hydro code (Masset 2000) modified by Picogna & Kley (2015). We model the dust as individual particles, sampling the evolution of a larger ensamble of single dust particles and their interaction with the surrounding environment.

In Figure 2 we show the drift timescales for particles of sizes between 100$\mu$m and 1 mm just before and just after the opening of the gap, when the gas mass accretion rates are ????? and ?????, respectively. At about 0.5 au particles of 880$mu$m have drift timescales of $\sim$ ?????? yr. The 2D 880$\mu$m particle distribution at the same evolutionary times are shown in Figure 3.

The drift timescales vary steeply with particle size in the vicinity of the gap, however all particles are quickly removed from the edge of the gap and pushed into the inner 0.5 au, where the drift timescales then start slowing down significantly. It is thus likely that some of the 880 $\mu$ size particles still linger in the inner gas disk of TW Hya and contribute to the emission of the unresolved 1mJy flux detected by Andrews et al. (2016).

Conclusions
===========

Acknowledgements
================

We thank people

[99]{}

Alexander, R. D., & Pascucci, I. 2012, , 422, L82

Andr[é]{}, P., Men’shchikov, A., Bontemps, S., et al. 2010, , 518, L102

Andrews, S. M., Wilner, D. J., Espaillat, C., et al. 2011, , 732, 42

Armitage, P. J., & Hansen, B. M. S. 1999, , 402, 633

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, , 558, A33

Bjorkman, J. E. 1997, Stellar Atmospheres: Theory and Observations, 497, 239

Bustamante, I., Mer[í]{}n, B., Ribas, [Á]{}., et al. 2015, arXiv:1501.05204

Calvet, N., D’Alessio, P., Watson, D. M., et al. 2005, , 630, L185

Christiaens, V., Casassus, S., Perez, S., van der Plas, G., & M[é]{}nard, F. 2014, , 785, L12

Cohen, M., Wheaton, W. A., & Megeath, S. T. 2003, , 126, 1090

Comer[ó]{}n, F. 2008, Handbook of Star Forming Regions, Volume II, 295

Ercolano, B., Drake, J. J., Raymond, J. C., & Clarke, C. C. 2008, , 688, 398

Ercolano, B., Clarke, C. J., & Drake, J. J. 2009, , 699, 1639

Ercolano, B., Clarke, C. J., & Hall, A. C. 2011, , 410, 671

Ercolano, B., & Rosotti, G. 2015, arXiv:1504.03433

Espaillat, C., Calvet, N., D’Alessio, P., et al. 2007, , 664, L111

Espaillat, C., D’Alessio, P., Hern[á]{}ndez, J., et al. 2010, , 717, 441

Espaillat, C., Furlan, E., D’Alessio, P., et al. 2011, , 728, 49

Espaillat, C., Muzerolle, J., Najita, J., et al. 2014, Protostars and Planets VI, 497

Gauvin, L. S., & Strom, K. M. 1992, , 385, 217

Gorti, U., & Hollenbach, D. 2009, , 690, 1539

Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, , 495, 385

Hunte, J.D. 2007, Computing In Science & Engineering, 9, 90

Keane, J. T., Pascucci, I., Espaillat, C., et al. 2014, , 787, 153

Kenyon, S. J., & Hartmann, L. 1995, , 101, 117

Kim, K. H., Watson, D. M., Manoj, P., et al. 2009, , 700, 1017

Koepferl 2012, Master Thesis, University Observatory at the Faculty of Physics of the Ludwig-Maximilians University Munich

Koepferl, C. M., Ercolano, B., Dale, J., et al. 2013, , 428, 3327

Koepferl, Robitaille, Morales, 2015, ApJ, 799, 53

Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, , 495, 385

Luhman, K. L. 2007, , 173, 104

Luhman, K. L., Allen, L. E., Allen, P. R., et al. 2008, , 675, 1375

Luhman, K. L., Allen, P. R., Espaillat, C., Hartmann, L., & Calvet, N. 2010, , 186, 111

Lynden-Bell, D., & Pringle, J. E. 1974, , 168, 603

Owen, J. E., & Clarke, C. J. 2012, , 426, L96

Owen, J. E., Ercolano, B., Clarke, C. J., & Alexander, R. D. 2010, , 401, 1415

Owen, J. E., Ercolano, B., & Clarke, C. J. 2011, , 412, 13

Owen, J. E., Clarke, C. J., & Ercolano, B. 2012, , 422, 1880

Owen J. E., Wu Y., 2016, ApJ, 817, 107

P[é]{}rez, L. M., Isella, A., Carpenter, J. M., & Chandler, C. J. 2014, , 783, L13

Picogna, G., & Kley, W. 2015, , 584, A110

Pringle, J. E. 1981, , 19, 137

Ribas, [Á]{}., Mer[í]{}n, B., Bouy, H., et al. 2013, , 552, A115

Rice, W. K. M., Armitage, P. J., Bonnell, I. A., et al. 2003, , 346, L36

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, , 167, 256

Rosotti, G. P., Ercolano, B., Owen, J. E., & Armitage, P. J. 2013, , 430, 1392

Rydgren, A. E. 1980, , 85, 444

Shakura, N. I., & Sunyaev, R. A. 1973, , 24, 337

Skrutskie, M. F., Dutkevitch, D., Strom, S. E., et al. 1990, , 99, 1187

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, , 97, 1451

Whitney, B. A., Indebetouw, R., Bjorkman, J. E., & Wood, K. 2004, , 617, 1177

Whitney, B. A., Wood, K., Bjorkman, J. E., & Cohen, M. 2003, , 598, 1079

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, , 140, 1868

\[lastpage\]

[^1]: E-mail: ercolano@usm.lmu.de (BE)

[^2]: Steeper indices are obtained by e.g. Franco et al. (2000), but the latter are unlikely in this case and more appropriate to hypercompact HII regions.

[^3]: Note that CO isotopologue mass measurements are much lower, due to the probable depletion of volatile C in the disc. e.g. Schwarz et al. 2016; Miotello et al. 2016; Kama et al. 2016

[^4]: It is important to note at this point, when talking about transition discs demographics, that transition discs are a very diverse class of objects. We refer here to objects that are on the verge on dispersal and have developed a hole as a consequence of photoevaporation, which do not include the class of mm-bright “transition” discs, which accrete vigorously and simultaneously have large holes. It is unclear if the latter are objects undergoing a fast [*transition*]{} phase at all; they may be more long-lived but rare systems, whose large gaps have been carved by dynamical interaction with a companion (see also discussion by Owen & Clarke, 2012).
