protoplanetary disks are a product of the formation of protostars from the collapse of molecular cloud cores .
these disks remain a source of ongoing accretion for the central stellar object for a period of up to @xmath9  myr ( manning & sargent , 1997 ; wyatt , dent & greaves , 2003 ) and are the birthplaces of planetary systems .
the evolution of the disk and of its planetary progeny is controlled by heating and irradiation from the central star  from ir to x - rays  but in the vicinity of ob stars can also be significantly affected by the external radiation environment ( e.g. johnstone et al .
1998 ; adams et al . , 2004 ; fatuzzo & adams , 2008 ) .
the discovery of the `` proplyd '' phenomenon and subsequent observations of protoplanetary disks represent an important and growing part of the legacies of the _ hubble _ and _ spitzer _ space telescopes ( e.g. odell 1993 , bally et al .
1998 , evans et al .
2003 , meyer et al .
new instruments such as the _ mid - infrared instrument _ on board _ james webb space telescope _ , which will provide both imaging and spectroscopy , hold the promise for higher spatial resolution . on the ground , millimeter and submillimeter studies have evolved from examining disk size and orientation ( e.g. keene & masson 1990 , lay et al .
1994 ) to detailed studies of structure , chemistry and thermal properties ( e.g. wilner et al .
2003 , qi et al .
2006 ) , while future high spatial resolution submillimeter observations with alma promise further major advances ( e.g. van dishoeck & jrgensen 2008 ) .    these and other observations should provide insights into issues central to understanding disks and the likelihood of forming planetary systems like our own , such as angular momentum transport , and the timescales for grain growth and gas dissipation . of particular interest for interaction of the disk with its parent star and for the study of planetary origins
are the inner regions of disks ( @xmath10  au ) where both terrestrial and giant planets are believed to form .    while the observational database continues to grow , all attempts at theoretically modeling the emission expected from such environments have lacked several of the ingredients needed for a fully self - consistent approach .
the multidimensional , far - ultra violet ( uv ) to x - ray radiative transfer ( rt ) problem must be solved simultaneously with the disk hydrodynamical evolution , to which it is strongly coupled .
furthermore , disks comprise adjacent regions of ionised , photon - dominated and molecular gas intermixed with a non - uniform dust component .
no computer code currently available includes all the important physics and chemistry needed to properly interpret the observational diagnostics of these different regimes .
furthermore , full theoretical understanding of a number of important contributors to the energetics of the system , such as accretion and wind - disk interactions , is still lacking ; these ingredients can currently only be included through approximate and phenomenological approaches ( e.g.see discussion in glassgold et al . , 2004 ) .
several studies have focused on the properties of the dust component of disk atmospheres ( e.g. dalessio et al .
2001 , dullemond et al .
2002 , akeson et al 2005 ) .
however , studies of the gaseous component of the inner disks have only recently been attempted , and only an exploratory picture of the thermochemical structure of this crucial component has emerged . pioneering work from igea & glassgold ( 1999 ) , glassgold et al . , ( 2004 , 2005 ) and meijerink et al . (
2008 ) examined physical conditions and expected emission lines in the gaseous component of non - evolving t  tauri disks irradiated by x - rays .
nomura & millar ( 2005 ) modeled molecular hydrogen emission from protoplanetary disks , taking into account uv radiation from a central star ; their models , while calculating the density and temperature structures self - consistently , did not account for x - radiation and employed an approximate treatment of the scattered light that often leads to an overestimation of the uv radiation field in the disk ( glassgold , 2006 ) .
semenov et al .
( 2004 ) employed perhaps the most complete chemical network to date to calculate the ionisation fractions at the disk midplane . while x - ray ionisation was included , radiative transfer through scattering and diffusion was not .
simple column density arguments were used instead to estimate the radiation field attenuation at discrete points in the disk where the chemical model was applied .
very recently gorti & hollenbach presented 1 + 1d models of gas in optically thick disks , calculating thermal , density and chemistry structures , with special attention to the effects of fuv irradiation .
more detailed 2d or 3d calculations are now needed to provide a quantitative picture of these environments in relation to their effects on the disk spectral energy distribution ( sed ) , and on other observables available now and in the near future ( e.g. @xmath11o  i@xmath12 $ ] , @xmath11ne  ii@xmath12 $ ] , co , oh , and uv transitions of h@xmath13 ; najita et al , 2007 ) .
in particular , observations of forbidden lines in t  tauri stars provide evidence for a tenuous , hot outer layer , or  corona ( kwan & tademaru 1988 , 1995 ; kwan ( 1997 ) . the earlier work by by glassgold et al .
( 2004 , 2007 ) , semenov et al . (
2004 ) and meijerink et al .
( 2008 ) emphasised the chemistry in the photo - dissociation region ( pdr ) and molecular zones , and stopped short of this coronal layer .
also the recent work by gorti & hollenbach does not include the fully ionized layer discussed here .
in this paper , which is the first of a series aimed at building a more realistic and self - consistent model of irradiated gaseous t  tauri disks , we employ a fully 3d photoionisation and dust radiative transfer code , mocassin ( ercolano et al .
, 2003 , 2005 , 2008 ) , to calculate the detailed ionisation and temperature structure of a typical t  tauri disk , with special emphasis on the photoionisation - dominated outer layers and corona .
we show that hot coronal temperatures are expected from x - ray irradiation from the central pre - main sequence star .
our calculations aim to understand the thermochemical structure of the irradiated region and to identify emission lines that can be used as gas - phase diagnostics for current and future observations .
furthermore , the detailed 2d temperature structure calculated by our models allows us to estimate the efficiency of x - radiation to drive a photoevaporative wind .
our model is described in section 2 .
the ionisation and temperature structures are presented and discussed in section 3 , while the predicted emission line spectra are shown in section 4 .
section 5 contains a discussion of our results in the context of disk dispersal via x - ray photoevaporation .
a brief summary is given in section 6 .
photoionisation and temperature structure calculations were performed with version 3.02.00 of the 3d monte carlo photoionisation and dust radiative transfer code , mocassin ( ercolano et al .
2003 , 2005 , 2008 ) .
the code uses a stochastic approach to the transfer of radiation allowing for the self - consistent transfer of both the primary and secondary components of the radiation field .
the code was modified to include viscous heating , approximated using the prescription for a thin disk by pringle ( 1981 ) , and gas cooling by collisions of grains with a mixture of atomic and molecular hydrogen ( hollenbach & mckee , 1979 ) ; we discuss this further below .
mocassin was also adapted to run in 2d , using the gas and dust density distribution given by a 2d protoplanetary disk calculation of dalessio et al .
( 1999 ) .
the dalessio model used in our calculation was chosen to be that which best fits the median sed of t  tauri stars in taurus ( dalessio , 2005 ) .
we refer the reader to dalessio et al .
( 1999 ) for a detailed description of the model ingredients and calculation . in brief , the input parameters for this model include a central star of 0.7  m@xmath6 , 2.5  r@xmath6 irradiating the disk with an effective temperature of 4000  k. additional disk parameters consist of a mass accretion rate of 1.e-8 m@xmath6 yr@xmath8 and a viscosity parameter @xmath14  =  0.01 .
the total mass of the disk is 0.027  m@xmath6 .
the present study is particularly focused on determining the physical properties of the gas in the outer layers and hot disk corona .
molecular species that are expected to form deeper in the disk are not included in our model ; this poses a limitation on the accuracy of our gas temperature calculations in the inner layers .
in particular , the calculations by glassgold et al . , ( 2004 , gni04 ) have shown that co rovibrational lines may become important coolants , and , e.g. , at a radial distance of 1  au and column densities of 10@xmath15 - 10@xmath16  @xmath17 co was found to dominate the disk cooling .
co rotational lines are also expected to form at greater depths , though in these regions the thermal balance is dominated by dust - gas collisions ( gni04 ) .
dust - gas collisions can heat the gas / dust if the gas is cooler / warmer than the dust ; photoelectric emission from dust grains heats the gas and it provides a minor cooling channel for the dust .
mocassin is able to self - consistently calculate the dust and gas thermal balance , taking into account the main microphysical processes that couple the two phases . here
, however , we choose to keep the dust temperatures fixed to the values calculated by dalessio et al.(1999 ) and focus our discussion on the gas phase . in the current work
we do not treat the transfer of the stellar or interstellar far - uv field , and therefore the dust temperatures we would obtain from our models would be incorrect .
the effects of gas - dust interactions on the temperature structure of the gas , however , are still taken into account in our gas thermal balance as well as in the radiative transfer , where the competition between dust and gas for the absorption of x - ray radiation is properly treated .
we illuminate the disk using synthetic x - ray spectra representative of the high energy emission from a typical t  tauri star .
there are currently at least three different mechanisms by which t  tauri stars are thought to produce observable x - rays : magnetospheric accretion shocks , shocks resulting from polar jets , and magnetically - confined hot coronae analogous to those found on all late - type main - sequence stars including the sun ( e.g. kastner et al 2002 , gdel et al 2005 , walter & kuhi , 1981 ) .
the first two of these have only been identified in a small handful of stars , and since the latter appears to dominate the observed x - ray spectra of t  tauri stars we currently ignore the effects of x - rays from jets and accretion .
coronal x - ray emission arises from an optically - thin plasma dominated by impact excitation and ionization by thermal electrons .
the spectra from such a plasma comprises significant contributions from bound - bound , bound - free and free - free radiation .
we computed synthetic coronal spectra incorporating these processes for the energy range 13.6  ev-12.5  kev for all elements with atomic number z=1 to 30 using line and continuum emissivities from the chianti compilation of atomic data ( landi et al .
, 2006 , and references therein ) , together with ion populations from mazzotta et al .
( 1998 ) , as implemented in the pintofale idl software suit ( kashyap & drake , 2000 ) .
we adopted the solar chemical composition of grevesse & sauval ( 1998 ) for the calculated spectra and an x - ray temperature , @xmath18 .
we assume an x - ray luminosity , @xmath19(0.1 - 10  kev)@xmath20  erg / sec .
details on atomic data and physical processes used for the photoionisation calculations are given by ercolano et al .
( 2003 , 2008 ) and ercolano & storey ( 2006 ) .
the chemical composition of grevesse & sauval ( 1998 ) was also used for the gas abundances of the most common elements used in our photoionisation model of the disk corona . in view of the ongoing debate on the `` solar oxygen crisis '' ( e.g. ayres et al .
2006 ; drake & testa , 2005 ; socas - navarro & norton , 2007 ; socas - navarro & centeno , 2008 ) , we have chosen not to use the more recent compilations of solar abundances based on 3-d non - lte hydrodynamic photospheric modeling ( e.g. asplund et al . 2005 ) , which recommends 25 - 35% lower c , n , o and ne abundances compared to earlier assessment .
we also ignore the possible effects of grain depletion .
although some depletion of refractory elements onto dust grains is expected , the amounts are very uncertain and probably quite small in the disk corona where the x - ray field is more intense . here , undepleted values of c , mg , si and fe abundances correspond to efficient destruction of grains in the upper layers of the disk
. a test for the validity of this assumption could be provided by the si  ii  34.8@xmath21 m line , which is a good diagnostic for grain depletion , and which our models predict to be observable ( see table  1 ) .
the dust absorption and scattering coefficients are calculated for high energies using the dielectric constants for graphite and silicates of laor & draine ( 1993 ) , which extend to the x - ray domain .
we assume spherical grains and use standard mie scattering series expansion for @xmath22 , where @xmath23 is the complex refractive index and @xmath24 is the scattering parameter ( see laor & draine , 1993 ) . for @xmath25 and @xmath26
we use rayleigh - gans theory ( bohren & huffman , 1983 ) , and for @xmath25 and @xmath27 we use the treatment specified by laor & draine ( 1993 ) , which is based on geometric optics .
our dust model consists of a typical ism mixture of graphite and silicates with mrn size distribution ( mathis , rumpl & nordsiek , 1977 ) described by @xmath28 m and @xmath29 m .
the chosen grain size distribution does not take into account grain growth , and therefore over - emphasizes the cooling of the gas by the dust in the disk interior .
we also assume a dust to gas mass ratio of 0.01 throughout the disk .
we note that the density distribution model of dalessio et al .
( 1999 ) uses different dust prescriptions at different heights , in order to account for the effects of grain settling and growth . for the sake of simplicity , and because we assume the dust temperatures calculated by dalessio et al .
( 1999 ) , we chose to use a more standard dust model throughout our disk .
since matching a particular sed is not one of our main objectives , and we are mostly interested in the physical properties of the disk corona , the choice of the size distribution and species is not crucial in our work as it will only have a very small effect on the ionisation and temperature structure of the coronal gas .
as noted earlier , in addition to all the standard heating and cooling channels included in mocassin , we also include viscous accretion heating calculated assuming the standard dissipation rate for a thin disk ( pringle , 1981 ) @xmath30 where @xmath31 is the local mass density of the gas , @xmath32 is the isothermal sound speed and @xmath33 is the angular rotation speed . the parameter @xmath14
= 0.01 relates viscous heating to the gas pressure .
the formulation above is purely phenomenological and carries large uncertainties . in particular ,
the appropriate value of @xmath14 is unclear when the formula is applied to the warm disk atmosphere ( gni04 ) , where it is unknown if angular momentum transport such as the magneto - rotational instability ( mri ; balbus & hawley , 1991 ) even operates at all .
gni04 show the effect that varying the value of @xmath14 has on the temperature structure of the gas , but finally adopt a model with @xmath14  =  0.01 for their line emission predictions ( glassgold et al . 2007 ;
meijerink et al .
2008 ) . a value of @xmath34 is also chosen by dalessio et al .
( 1998 , 1999 , 2001 ) for their disk structure calculations ( which includes the model used in this paper ) .
this choice is is based on the connection between @xmath14 and the disk accretion rate : @xmath35 , where @xmath36 is the integrated surface density of the disk . with @xmath34 , which is also the value implemented in our models ,
x - ray heating is dominant in the disk corona , while accretion heating dominates in the inner regions .
we finally note that the dependence of the heating rate on the square of the sound speed implies a _ linear _ dependence on the gas temperature , which poses potential problems for the hot disk corona . in the regions
where the gas is heated to high temperatures by x - ray irradiation the functional form of accretion heating , as formulated above , fails completely , preventing an equilibrium temperature from being found .
in fact , x - ray irradiation in the upper atmosphere of the inner regions of the disk heats the gas to temperatures above the local escape temperature of the gas ( see section  [ s : disp ] ) .
the gas in these regions is to be considered unbound and part of a photoevaporative flow , hence not subject to the standard viscous accretion law . for unbound gas
we remove the accretion heating contribution altogether by setting the local @xmath14 coefficient to zero in our models .
in this section , we present the temperature and ionisation structure of the gas in the inner disk , focusing our attention on the disk corona .
we limit our discussion to radial distances @xmath37  au , beyond which , in our model , sets a rough outer limit to the region where the effects of x - ray irradiation on the temperature and ionisation structure are more evident .
we note , however , that in our models a low levels of ionisation is maintained out to a radius of @xmath0190  au .
figure  [ f : f1 ] shows the temperature structure ( asterisks ) in the vertical column calculated at @xmath38  =  0.07  au , together with the the dust temperature profile calculated by dalessio et al .
( 1999 ) . also shown
is the volumetric gas density @xmath11@xmath39 $ ] as a function of depth , the corresponding geometrical height above the midplane is shown in the top x - axis .
the gaseous disk corona at these very small radii is fully ionised .
since all elements are stripped of their outer electrons , cooling by collisionally excited lines is inefficient .
the balance between heating by photoionisation and cooling by line emission and recombination processes corresponds to a gas equilibrium temperature of @xmath010@xmath1k .    at 0.07  au
, the jump in gas temperature occurs at a vertical column density of @xmath40  @xmath17 , where the gas density is @xmath41 and the ionisation parameter , defined as @xmath42  =  l / nr@xmath43 , is log(@xmath42 )  @xmath01.78 .
this transition would look even sharper in a hydrostatically balanced model , since a rise in the temperature would induce a drop in density in order to mantain pressure equilibrium . in this case ,
any equilibria found in the 10@xmath2 - 10@xmath1k range may in fact be unstable ( hatchett , 1976 ) .
this has already been noticed in accretion disk models around x - ray binaries ( e.g. raymond 1993 ) .
alexander , clarke & pringle ( 2004 , acp04 ) have also presented a model to illustrate the effects of x - ray heating on the structure of circumstellar disks . in that work
, the density profile was truncated at the point where the density fell to @xmath44@xmath45 , which is below the fully ionised region shown in our figure  [ f : f1 ] .
this should be kept in mind when comparing our temperature structure to theirs ( as shown in their figure  3 ) , where the lack of the @xmath010@xmath1k temperature ( fully ionised ) region is not surprising and is not discrepant with our results .
the electron density distribution , shown in figure  [ f : f2 ] , indicates that the a high level of ionisation at 0.07  au is maintained for a column density of 10@xmath46  @xmath17 , after which it decays to lower values .
a study of the thermal balance at this radius indicates x - ray heating dominates for a column depth of @xmath010@xmath16  @xmath17 , while heating by viscous accretion dominates further down in the disk where it is balanced by dust - gas collisional cooling . in our calculations , the contribution from accretion heating is set to zero in regions where x - ray irradiation heats up the gas to temperatures above the local escape temperature ( see section  [ s : model ] ) . at @xmath47
au this occurs at a column depth of @xmath010@xmath48  @xmath17 , at the onset of the temperature discontinuity.    figure  [ f : f4 ] shows the temperature and volumetric density @xmath11@xmath39 $ ] structures in the vertical column calculated at @xmath38  =  1  au , together with the dalessio et al .
( 1999 ) dust temperature profile .
the coronal gas here appears to be only partially ionised and the electron density , illustrated in figure  [ f : f5 ] , reaches a level of just over 10% of the hydrogen density , but decreases very rapidly with column depth .
our temperature structure at 1  au is comparable to that found by gni04 , allowing for differences in the model input parameters , and bearing in mind that their calculation started at a column depth of 10@xmath49  @xmath17 .
however , the cooling and heating contributions calculated by the two models disagree somewhat and deserve further explanation . the various heating and cooling channels in our model at 1  au
are illustrated as a function of vertical column depth in figure  [ f : f6 ] , which can be compared with an analogous illustration in figure  3 of gni04 .
one conspicuous difference between the two models is the accretion heating contribution , which is simply due to the different choice of @xmath14 ( 0.01 in this work and 1.0 in the particular model plotted in figure  3 of gni04 ) .
note that gni04 , also investigated models with lower values of @xmath14 ( see their figure  2 ) , and indeed used a model with @xmath14  =  0.01 , for their later work ( glassgold et al .
, 2007 ; meijerink et al . ;
2008 ) , where emission line predictions are given .
unfortunately a plot showing the heating / cooling contributions from their models is only available for the @xmath14  =  1 case .
another important difference lies in the ly@xmath14 cooling contribution , which is only of secondary importance in our models but dominates the cooling down to a vertical column of 10@xmath15  @xmath17 in the gni04 model .
the dominant cooling mechanism in the upper disk atmosphere in our model is metal line emission .
the ly@xmath14 cooling in gni04 was overestimated because of an assumption of a thermal equilibrium population for the n  = 2 level .
the `` ly@xmath14 '' channel in the gni04 model also includes cooling by the fine structure lines of o  i , while other atomic coolants are ignored . in spite of the aforementioned differences , the total cooling calculated by gni04
is very similar to that obtained in this work for the regions where the calculations overlap and when the same @xmath14 parameter is employed ( glassgold , 2008 ; private communication ) .
as we move to increasingly larger radii , up to approximately 25 - 30  au , our model disk structure remains qualitatively similar to that at 1  au : high gas temperatures are reached in the corona , followed by a steep temperature decrease , with the gas and dust becoming thermally coupled at a vertical column depth of @xmath05@xmath5010@xmath15@xmath17 . beyond @xmath035 - 40  au , due to geometrical dilution of the central radiation field , the effects of x - ray irradiation become less evident , however a low level of ionisation is maintained out to a radius of @xmath0190  au . for t - tauri stars born in clusters containing massive stars , like the orion nebula cluster , the dominant irradiation mechanism at large disk radii
is probably provided by the interstellar or intracluster fuv and x - ray fields , which are not included in our current models .
this is not true for stars born in smaller clusters like @xmath31 ophiuchi and taurus - auriga , where indeed large disks ( @xmath0500 - 1000  au ) are found ( andrews & williams , 2007 ) .    the maximum temperature reached by the disk corona at a given radial distance is an immediate consequence of the local ionisation parameter .
the fact that at 1  au we do not see the @xmath010@xmath1k gas layer is simply due to the fact that the density distribution of dalessio et al . ( 1999 ) used in our models is truncated at a volumetric density of @xmath010@xmath2 . at 0.07  au , the temperature jump to @xmath010@xmath1k occurs at a log(@xmath42 )  =  1.78 .
the same value at 1  au implies a density of @xmath0150  @xmath45 , which is below the cut - off .
the general results found in our models are similar to those found by gorti & hollenbach ( 2008 ) , in particular that the gas and dust are thermally de - coupled down to a column of @xmath51 .
as an aid to future observational campaigns , we have calculated the emission line spectrum predicted by our model .
line luminosities are integrated out to a radius of 190  au . in tables
[ t : t1 ] and  [ t : t2 ] we list a subset of the strongest collisionally excited metal line and hydrogen recombination lines , respectively , down to a luminosity of @xmath010@xmath5l@xmath6 , which corresponds to a flux of @xmath01.8@xmath5010@xmath52erg / s / cm@xmath43 at a distance of 140pc .
the flux at a distance @xmath54 , assuming a face - on orientation for the disk , is obtained by dividing the luminosities listed by the appropriate dilution factor ( 4@xmath55 ) .
he  i recombination lines are found to be rather weak , with the exception of he  i  10.8@xmath21 m which has a luminosity of 8.8@xmath5010@xmath56l@xmath6 , he  i  7067 with luminosity of 3.2@xmath5010@xmath5l@xmath6 and he  i  3890 , he  i  5877 , he  i  6688 and he  i  20.5@xmath21 , all with luminosities of @xmath02@xmath5010@xmath56l@xmath6 .
the full set of lines , including over 2000 transitions is available on request from be .
the line list will be further improved and updated in a future work , when the effects of simultaneous stellar and interstellar fuv irradiation , as well as x - rays , will be considered .
some changes to the emission line intensities are also expected when the x - irradiation will be studied in the context of a self - consistent hydrostatic equilibrium calculation .
tables  [ t : t1 ] and  [ t : t2 ] provide an indication of the lines that should be targeted for observations of the gaseous phase of disk atmospheres
. a number of potential atomic and ionised gas diagnostic lines should already be detectable with current instrumentation .
the @xmath11ne  ii@xmath12 $ ]  12.8@xmath21 m line , in particular , has already been detected in t  tauri disks with _
spitzer _ by pascucci et al .
( 2007 ) , lahuis et al .
( 2007 ) and espaillat et al .
( 2007 ) , and with michelle on gemini north by herczeg et al .
the detections to date range from a few 10@xmath56 to 10@xmath57 l@xmath6 , in qualitative agreement with values predicted by our model .
the @xmath11ne  iii@xmath12 $ ]  15.5@xmath21 m line has so far only been detected in sz102 ( lahuis et al . 2007 ) , implying a @xmath11ne  iii@xmath12$]/@xmath11ne  ii@xmath12 $ ] of @xmath00.06 for this object , which is some 2.5 times larger than our model . as noted by glassgold et al .
( 2007 ) the abundance of ne@xmath58 is lowered by charge exchange reactions with neutral h , implying low ( @xmath590.1 ) ratios of @xmath11ne  iii@xmath12$]/@xmath11ne  ii@xmath12 $ ] in the x - ray heated gas .
gorti & hollenbach ( 2008 ) also discuss the effects of euv irradiation which would enhance the predicted @xmath11ne  ii@xmath12 $ ] and @xmath11ne  iii@xmath12 $ ] line luminosities , and affect their ratio .
they predict that the @xmath11ne  iii@xmath12$]/@xmath11ne  ii@xmath12 $ ] also depends on the details of the euv spectrum , with hard euv spectra naturally producing the higher ratios .
table  [ t : t1 ] clearly shows that a number of optical lines have fluxes that are significantly larger than the @xmath11ne  ii@xmath12 $ ]  12.8@xmath21 m fluxes .
unfortunately many optical and uv lines are not observable with current technology , as they are invisible against the stellar spectrum . in particular
, we predict a luminosity of 3.3@xmath5010@xmath57 for the mg  i line at 4573@xmath60 , for undepleted mg abundances .
this corresponds to only a few % of the continuum flux level for a photosphere of 4000k and for an assumed resolving power of 50,000 . while this line might be detectable under favourable circumstances and where disk gas is not strongly depleted , we are unaware of any such detections in the literature .
the predicted absolute flux of the various lines depends , of course , on the details of the disk structure , the gas abundances and the ionising spectrum . as an example , meijerink et al .
( 2008 , mgn08 ) , show , for a fixed disk structure , the dependence of the predicted flux of a number of fine structure and forbidden lines on the assumed x - ray luminosity .
furthermore , lines with very low critical densities like the c  i and c  ii fine structure lines ( see table  [ t : t4 ] ) are also very sensitive to conditions in the outer regions of the disk , where densities are small .
these regions are preferentially illuminated by the interstellar ( or intracluster ) radiation fields which are not included in our current simulations .
unfortunately , table  [ t : t1 ] is not yet complete , as reliable collision strengths were not available for all species , and a number of less abundant elements were omitted from our calculations in order to limit computational costs . in particular , the following species were not included , due to lack of atomic data : ne  i , mg  iii , si  i , s  i , fe  i , fe  iv , fe  ix , fe  xi , fe  xiii . due to its complexity , the fe  ii ion was also omitted from the present calculations but will be included in future models .
a further observational complication lies in the difficulty in separating the emission from the disk , the wind , the accretion funnels and the photosphere .
high resolution spectroscopy offers one possibility , in which line profiles that are shifted and/or broadened to different extents in the different emission regimes can be resolved ( e.g. herczeg et al . , 2007 ) .
following the suggestion in mgn08 , in future work we will present predicted profiles for a number of potentially observable lines emitted from different regions of the disk .
an alternative way to distinguish lines formed in the upper layers of the disk and those formed in the accretion funnel would involve a comparison of transitions with different critical densities ( e.g. o  i6302 with n@xmath61  =  1.5@xmath5010@xmath62 , compared to o  i63@xmath21 m with n@xmath61  =  2@xmath5010@xmath63 ) .
we note that some of our predicted line intensities differ somewhat from those calculated by mgn08 . in most cases ,
the differences may be explained by the fact that our work used a different disk model for the underlying gas density distribution , and different gas abundances . by way of comparison
, we have also computed a model that uses the same input parameters and disk model as those of mgn08 and we compare the predicted line fluxes to those of mgn08 in table  [ t : t4 ] .
most of the lines listed are in reasonable agreement , particularly when one considers the two different approaches to the problem and the use of different atomic data .
mgn08 predict larger @xmath11ne  ii@xmath12 $ ]  12.8@xmath21 m line luminosities than in our model ; this is explained by the fact that mgn08 do not include cooling by @xmath11ne  ii@xmath12 $ ] line , which , as already noted by gorti & hollenbach ( 2008 ) is found to be important . as a consequence
the gas temperature in the @xmath11ne  ii@xmath12 $ ] emitting region will be higher in mgn08 s calculations , producing a larger @xmath11ne  ii@xmath12 $ ]  12.8@xmath21 m luminosity .
the c  ii fine - structure line at 158@xmath21 m shows a larger discrepancy and deserves further attention . because the critical density of this line is low ( see table  [ t : t4 ] ) , it will be collisionally quenched at even modest depths into the disk , and will be preferentially emitted in the less dense upper layers of the disk corona .
mgn08 calculations , however , were focused on the characterisation of lower layers in the disk and stopped short of the corona , with the result that some of the flux from the c  ii bright regions is probably missed .
a detailed comparison of the line luminosities predicted by gorti & hollenbach ( 2008 ) to those predicted by our models is hindered by the different parameters and approaches employed by the two models .
however we note that the @xmath11ne  ii@xmath12 $ ]  12.8@xmath21 m luminosity predicted by our models is a factor or 3.3 lower than that derived by their fiducial model ( model a , their table  4 ) .
this can be explained by the fact that our density distribution is determined by the dust temperature structure .
indeed gorti & hollenbach ( 2008 ) comment that when they assume a dust - determined temperature structure they obtain a factor of @xmath04 lower @xmath11ne  ii@xmath12 $ ] line luminosities .
we also note that our predicted c  i and c  ii fine - structure lines are higher by a factor of @xmath02 than those given by model  a in gorti & hollenbach ( 2008 ) .
as discussed above , these lines have low critical densities and are therefore extremely sensitive to the details of the density distribution in the emitting region , furthermore , we use a factor of @xmath03 higher elemental abundance for c. these two factors are probably responsible for the small discrepancy at hand .
gorti & hollenbach ( 2008 ) use highly depleted values for the si elemental abundance ; this is the likely cause of the lower ( factor of @xmath010 ) si  ii  35@xmath21 m line luminosities predicted by their models .
finally the factor of @xmath03 lower o abundance used in their model is probably one of the causes of the small ( factor of @xmath01.5 ) discrepancy with the o  i  63@xmath21 m line luminosity .
circumstellar disks are observed around the majority of young stars at @xmath010@xmath1  yr ( e.g. lada & lada 2003 ) . by ages of @xmath010@xmath64  yr , however , only low - mass debris disks are generally observed ( e.g. manning & sargent , 1997 ; wyatt , dent & greaves , 2003 ) , implying disk lifetimes of up to a few @xmath010@xmath1 yr . indeed , recent observations indicate that about 90% of low mass stars lose their primordial disks at 5 - 7 myr ( e.g. haisch et al . , 2001 ; hartmann , 2005 ; hernandz et al . , 2007a ) . however , very few observations exist of objects in transition between classical t  tauri stars , which still have disks of a few percent of a solar mass , and weak - lined t  tauri stars , which have dispersed nearly all of their circumstellar disks . for example , hernandez et al .
( 2007b ) reported only a 10% frequency of transitional disk candidates among the disk - bearing low mass stars in their _ spitzer _ observations of the orion ob1 association .
this suggests that the mechanism responsible for the final stages of disk dispersal must operate on relatively short timescales  perhaps as short as @xmath010@xmath65  yr ( e.g. duvert et al . , 2000 , andrews & williams , 2005 ) .
photoevaporative heating from both the central star and external sources has been posited as a viable mechanism able to satisfy the two - timescale nature of the problem .
a dramatic illustration of the latter is the apparent photoevaporation of orion `` proplyds '' in the vicinity of ob stars ( e.g. johnstone et al .
1998 , adams et al .
, 2004 ; fatuzzo & adams , 2008 ) .
evidence is also accumulating that disk loss is more rapid in more o - star rich clusters , as judged from observed disk frequencies in low - mass stars ( e.g. albacete - colombo et al .
2007 , guarcello et al .
2007 , mayne et al . 2007 ) .
we will return to photoevaporation by external sources of radiation in future work , and concentrate here on the effect on the disk of irradiation by the central t  tauri star .
photoevaporation by ultraviolet ( uv ) radiation from the central star has been investigated by a number of authors ( see alexander , 2007 ; dullemond et al , 2007 ; hollenbach , yorke & johnstone , 2000 , for recent reviews ) .
the original models did not couple photoevaporation with viscous evolution , resulting in disk dispersal timescales exceeding the observational estimates .
subsequently , the uv switch model of clarke et al .
( 2001 ) , which coupled photoevaporative mass loss to viscous evolution , was able to reproduce the two - timescale behaviour implied by the observations .
however a number of problems still remain , as pointed out in clarke et al .
in particular , current models require that the lyman continuum emission from the central star does not derive from accretion , an issue that is currently unclear ( acp04 ) .
x - ray ionisation and heating of the disk corona provide an attractive alternative engine to drive photoevaporative flows .
the fact that weak - line t  tauri stars , for which accretion has nearly ceased , are known to be brighter in the x - ray than classical t - tauri stars , which are still accreting , suggests that the production of x - ray is not linked to accretion .
acp04 estimated an upper limit to the mass - loss driven by x - rays and found it to be at best only comparable to that derived for ultraviolet evaporation .
they concluded , therefore , that x - ray heating is unlikely to be the dominant disk dispersal mechanism .
acp04 took into account the effects of x - ray photoionisation and heating on the hydrostatic structure of a typical t  tauri circumstellar disk , but were limited by the fact that the radiative transfer problem could not be solved self - consistently in 2d .
while our models do not yet iterate over the disk density structure , which would be undoubtedly modified by x - ray irradiation , we do calculate the irradiated thermal structure in 2d .
it is therefore worthwhile to provide an estimate of the mass - loss rates and dispersal timescales implied by our models .
we estimate mass - loss rates by locating , at each radial distance , the base of the photoevaporative flow , defined as the height where the local gas temperature exceeds the escape temperature of the gas , defined as @xmath66 , where @xmath67 is the stellar mass , @xmath38 is the radial distance .
this height is illustrated as a function of radial distance in figure  [ f : f9 ] .
the figure shows that , between 8 and 40  au , the height of the photoevaporating envelope follows a nearly linear increase with radial distance from the star ( slope @xmath00.9 ) .
to a first approximation , the mass loss rate per unit area can be estimated as @xmath68 , where @xmath31 and @xmath32 are , respectively , the gas density and the sound speed evaluated at the base of the flow .
figure  [ f : f8 ] ( left panel ) shows @xmath69 , the mass loss rate per unit area and ( right panel ) the mass loss rate from an infinitely thin annulus as a function of radial distance in the disk .
@xmath69 peaks at approximately 12  au , near the gravitational radius , @xmath70 for gas with sound speed of @xmath71 ( @xmath72 ) .
we use the definition of @xmath70 given by hollenbach et al .
( 1994 ) , where @xmath73 . the value at the peak is 3.1@xmath5010@xmath74 , this is about a factor of ten higher than the estimate of acp04 ( 2.6@xmath5010@xmath75 ) and the value obtained by hollenbach et al .
( 1994 ) for uv photoevaporation ( 3.9@xmath5010@xmath75 ) .
we note that both studies considered the case of a 1  m@xmath6 star ; however , since photoevaporation rates scale only weakly with mass ( @xmath76 for euv photoevaporation and almost linearly for x - ray photoevaporation ) , this is not sufficient to explain the order of magnitude difference with our rates .
as discussed in the next section , our model has a number of limitations which may result in uncertainties in the derived photoevaporation rates .
taking into account statistical uncertainties in the calculated temperatures alone , we find that our models are also consistent with peak values as low as @xmath08@xmath5010@xmath77 , only a factor of @xmath0 4 higher than previous estimates when corrected for the different stellar masses .
the effects of x - ray photoevaporation are only significant between 8 and 40  au .
there is also a small inner region between 0.03 and 0.16  au , where , in spite of the large escape temperatures , the intense x - ray field causes the disk to lose mass at a rate of 3 - 4@xmath5010@xmath75 .
given the small area at small radii , this does not contribute significantly to the mass loss .
the reason for the mass loss being mainly concentrated between 8 - 40  au is that this region lies far enough from the central star for the escape temperatures to become sufficiently low and not so far that the effects of x - ray irradiation become unimportant . in other words ,
the steep drop in @xmath69 inside 8  au corresponds to the location where the sound speed at @xmath010@xmath2  k drops below the keplerian speed . any mass loss at smaller radii comes from coronal gas .
in fact , for radial distances smaller than 8  au , there is also a thin layer of gas in the disk corona where x - ray irradiation causes temperatures to rise above the local escape temperatures .
however , the contribution of this region to the total mass loss budget is negligible .
we also note here that by analogy with the x - ray heated accretion disk models of begelman , mckee and shields ( 1983 ) the wind region might extend inward to about one tenth of the radius where the sound speed equals the escape speed . in that case
the wind in figure  [ f : f8 ] would begin at about 0.8  au and the total mass loss rate would be increased .
more recent works by liffman ( 2003 ) and font et al .
( 2004 ) also find that an euv - driven wind is launched from well inside the gravitational radius ( as defined by hollenbach et al . , 1994 ) .
the dispersal timescale implied by our model is of a few myr , for a 0.027  m@xmath6 disk , and the total mass loss rate is of the order of m@xmath78  m@xmath6yr@xmath8 .
the disk however , as well as being photoevaporated , is also accreting .
therefore , we can assume that the effects of x - ray photoevaporation would not be felt until its viscous accretion rate , m@xmath79 , has fallen below m@xmath80 , at which point the disk would empty rapidly within 10 - 20  au .
understanding the later evolution of the disk , after the region within 10 - 20  au has been emptied , requires us to determine whether the disk beyond @xmath020  au ( i ) has to accrete in to the 10 - 20  au region and then be photoevaporated or ( ii ) is more or less fixed in density profile and then overtaken by the outward propagating ring of mass loss .
behaviour of type i ( or ii ) results if the mass loss rate increases ( decreases ) as the hole size increases .
this question will be answered in a following study where we plan to solve for the density structure and the 2d x - ray transport self - consistently .
our current calculations give m@xmath78  m@xmath6yr@xmath8 , implying that disks with m@xmath81  m@xmath6yr@xmath8 should disperse rapidly and therefore be seldom observed .
however this is in disagreement with observations of a considerable number of systems in the 0.5 - 1  m@xmath6 range that are reported to have accretion rates lower than @xmath7  m@xmath6yr@xmath8 .
indeed accretion rates as low as @xmath82  m@xmath6yr@xmath8 are seen in the @xmath83 mass bin in figure  8 of gregory et al .
( 2006 ) , which , collecting data from a number of authors , shows the correlation between observed mass accretion rates and stellar mass .
this indicates that our predicted m@xmath84  m@xmath6yr@xmath8 is probably too high , and may point to a weakness of our model ( see discussion in the next section ) .
a number of important approximations were used in our models , which may lead to uncertainties in the estimates of the photoevaporation rates .
one major caveat of our model , already mentioned in this paper , is that the influence of x - ray photoionisation and heating on the disk density structure has been ignored .
the high gas temperatures reached in the inner regions of the disk would cause an expansion of the gas in the z - direction , as shown , for example , by acp04 .
indeed , their lower value of mass loss rate is perhaps caused by a larger attenuation of the x - ray field due to the puffed up inner disk material that is not accounted for by our models .
a hot and inflated inner disk could in fact cause a shadowing effect of the regions at larger radial distances , reducing the effects of x - ray photoevaporation there .
the magnitude of this effect can only be properly estimated via the simultaneous solution of the 2d radiative transfer problem and the 1d hydrostatic equilibrium in the z - direction , which is the aim of a follow - up paper .
further uncertainties are introduced by the use of the `` @xmath85 '' approximation in the calculation of photoevaporation rates . in section  [ s : disp ] we discuss how the base of the photoevaporative envelope is defined as the height where the local gas temperature exceeds the escape temperature of the gas . while this is generally a very good approximation in cases where the temperature and density vary sharply at the ionisation front , as for euv photoevaporation ( e.g. hollenbach et al . , 1994 ) , in the case of x - ray photoevaporation presented here , the distinction between the bound and unbound gas is less clear .
in other words , a small error in the determination of temperatures in our models , would lead to much larger uncertainties in the density , and hence in the estimated photoevaporation rates .
we conclude that while we have shown that x - rays have the potential to drive an outflow , our derived rates should be treated with caution and we defer firm conclusions to future work .
we have presented a 2d photoionisation and dust radiative transfer calculation of a typical t  tauri protoplanetary disk irradiated by 1.5@xmath5010@xmath86  k ( kt  =  1.36 kev ) x - ray radiation field from the central pre - main sequence star .
a detailed characterisation of the temperature structure of the upper disk corona was provided , which uncovered the existence of a layer of highly ionised , hot gas in the inner disk regions .
the maximum temperatures reached by the gas in the upper corona were mapped across the radial dimension of the disk , showing peaks exceeding 10@xmath1k in the very inner regions .
the hotter component of the gas in the inner disk and at radii between @xmath08 to @xmath040  au was shown to be unbound , and probably part of a photoevaporative flow .    to aid future efforts to observe the hot / warm gas component of t  tauri disks
, we provide a substantial list of predicted emission line luminosities , many of which should already be detectable with current instrumentation .
the temperature structure of the irradiated disk implies a total mass loss rate of @xmath87  m@xmath6yr@xmath8 and a disk dispersal timescale of a few myr for a disk mass of 0.027  m@xmath6 , which is the mass of the disk structure model we employed .
we discuss the limitations of our model and highlight the need for further calculations that solve the 2d radiative transfer problem simultaneously with the 1d hydrostatic equilibrium of disk structure in the polar direction .
we extend our warmest thanks to al glassgold , for insightful discussion and a thorough assessment of our results .
we also thank the anonymous referee and the editor eric feigelson for helpful comments that added to the clarity of the paper and the interpretation of the results .
we thank paola dalessio for providing us with the electronic data for the gas density distribution in the disk .
we also thank bruce draine for help with the x - ray dust opacity calculations .
we finally extend our thanks to roger wesson for helpful discussion .
jjd was supported by the chandra x - ray center nasa contract nas8 - 39073 during the course of this research .
be was partially supported by c__chandra _ _ grants go6 - 7008x and go6-7009x.the simulations were partially run on the cosmos ( sgi altix 4700 ) supercomputer at damtp in cambridge .
cosmos is a uk - ccc facility which is supported by hefce and stfc .