Continuum and line modelling of discs around young stars - I. 300 000 disc models for HERSCHEL/GASPS

We have combined the thermo-chemical disc code ProDiMo with the Monte Carlo radiative transfer code MCFOST to calculate a grid of similar to 300 000 circumstellar disc models, systematically varying 11 stellar, disc and dust parameters including the total disc mass, several disc shape parameters and the dust-to-gas ratio. For each model, dust continuum and line radiative transfer calculations are carried out for 29 far-infrared, sub-mm and mm lines of [O i], [C ii], 12CO and o/p-H(2)O under five inclinations. The grid allows us to study the influence of the input parameters on the observables, to make statistical predictions for different types of circumstellar discs and to find systematic trends and correlations between the parameters, the continuum fluxes and the line fluxes. The model grid, comprising the calculated disc temperature and chemical structures, the computed spectral energy distributions, line fluxes and profiles, will be used in particular for the data interpretation of the HERSCHEL open time-key program GASPS. The calculated line fluxes show a strong dependence on the assumed ultraviolet excess of the central star and on the disc flaring. The fraction of models predicting [O i] and [C ii] fine-structure lines fluxes above HERSCHEL/PACS and SPICA/SAFARI detection limits is calculated as a function of disc mass. The possibility of deriving the disc gas mass from line observations is discussed.