Abstract:Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to E-u/k = 4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms.Methods. The model consists of a spherical envelope with a power-law density structure and a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, ultraviolet irradiation of the outflow cavity walls, and small-scale C-type shocks along the cavity walls. Most of the model parameters are constrained from independent observations; the two remaining free parameters considered here are the protostellar UV luminosity and the shock velocity. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC 1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases (early Stage 0 to late Stage I) and physical characteristics such as luminosity and envelope mass.Results. The bulk of the gas in the envelope, heated by the protostellar luminosity, accounts for 3-10% of the CO luminosity summed over all rotational lines up to J = 40-39; it is best probed by low-J CO isotopologue lines such as (CO)-O-18 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by (CO)-C-12 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (> 99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level and within the context of our model, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.

Modelling H2 Infrared Emission of the Helix Nebula Cometary Knots

Kinematics of Molecular Hydrogen Emission from Planetary and Pre-planetary Nebulae

The Orion fingers: H$_2$ temperatures and excitation in an explosive outflow

Molecular Hydrogen in the Ionized Region of Planetary Nebulae

Temperature stratification in a molecular shock: Analysis of the emission of H2 pure rotational lines in IC443G

A classification of HII regions based on Oxygen and Helium lines: The cases of TOL2146-391 and TOL0357-3915

Modelling Herschel Observations of Hot Molecular Gas Emission from Embedded Low-Mass Protostars

The kinematics of the large western knot in the halo of the young planetary nebula NGC 6543

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

Pushing the limits: detecting H2 emission from faint bipolar planetary nebulae in the IPHAS sample

The Infrared Knots Around Ss 433

Modelling the abundance structure of isocyanic acid (HNCO) toward the low-mass solar type protostar IRAS 16293-2422

The excitation mechanism of H2 in bipolar planetary nebulae

H$_2$ molecular gas in the old planetary nebula NGC 3587

CO in the C1 globule of the Helix nebula with ALMA

Molecular shock tracers in NGC 1068: SiO and HNCO

Feedback of molecular outflows from protostars in NGC 1333, revealed by Herschel and Spitzer spectro-imaging observations

The excitation of near-infrared H2 emission in NGC 253

Structure and Physical Conditions in the Huygens Region of the Orion Nebula

Discovery of atomic and molecular mid infra-red emission lines in off-nuclear regions of NGC 1275 and NGC 4696 with the Spitzer Space Telescope

An analysis of the isomers HCN and HNC in the evolution of high-mass star-forming regions