Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

The Herschel Space Observatory was the fourth cornerstone mission in the European Space Agency (ESA) science programme. It had excellent broad band imaging capabilities in the far-infrared (FIR) and sub-millimetre part of the electromagnetic spectrum. Although the spacecraft finished observing in 2013, it left a large legacy dataset that is far from having been fully explored and still has a great potential for new scientific discoveries. The PACS and SPIRE photometric cameras observed about 8% of the sky in six different wavebands. This document describes the Herschel/PACS Point Source Catalogue (HPPSC), a FIR catalogue based on the broad-band photometric observations of the PACS instrument with filters centred at 70, 100 and 160 microns. We analysed 14842 combined, Level 2.5/Level 3 Herschel/PACS photometric observations. The PACS photometer maps were generated by the JScanam task of the Herschel Interactive Processing Environment (HIPE) v13.0.0. Sources were identified with the HIPE implementation of SUSSEXtractor, and the flux densities obtained by aperture photometry. We found a total of 108 319 point sources that are considered to be reliable in the 70 micron maps, 131 322 at 100 micron and 251 392 point sources in the 160 micron maps. In addition, our quality control algorithm identified 546 587 candidate sources that were found to be extended and 7 185 160 features which did not pass the signal-to-noise and other criteria to be considered reliable sources. These sources were included in the Extended Source List and Rejected Source List of the HPPSC, respectively. The calculated completeness and photometric accuracy values are based on simulations, where artificial sources were injected into the observational timeline with well controlled flux density values. The actual completeness is a complex function of the source flux, photometric band and the background complexity.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

The Herschel Space Observatory was the fourth cornerstone mission in the European Space Agency (ESA) science programme with excellent broad band imaging capabilities in the sub-mm and far-infrared part of the spectrum. Although the spacecraft finished its observations in 2013, it left a large legacy dataset that is far from having been fully scrutinised and still has a large potential for new scientific discoveries. This is specifically true for the photometric observations of the PACS and SPIRE instruments. Some source catalogues have already been produced by individual observing programs, but there are many observations that risk to remain unexplored. To maximise the science return of the SPIRE and PACS data sets, we are in the process of building the Herschel Point Source Catalogue (HPSC) from all primary and parallel mode observations. Our homogeneous source extraction enables a systematic and unbiased comparison of sensitivity across the different Herschel fields that single programs will generally not be able to provide. The catalogue will be made available online through archives like the Herschel Science Archive (HSA), the Infrared Science Archive (IRSA), and the Strasbourg Astronomical Data Center (CDS).

We observed the giant HII region around the NGC 3603 YC with the 5 broad bands (70, 160, 250, 350, 500 micron) of the SPIRE and PACS instruments, on-board the Herschel Space Observatory. Together with what is currently known of the stellar, atomic, molecular and warm dust components, this additional and crucial information should allow us to better understand the details of the star formation history in this region. The main objective of the investigation is to study, at high spatial resolution, the distribution and main physical characteristics of the cold dust. By reconstructing the temperature and density maps, we found respectively a mean value of 36 K and Log(NH)=22.0+/-0.1 cm^(-2). We carried out a photometric analysis detecting 107 point-like sources, mostly confined to the North and South of the cluster. By comparing our data with SED models we found that 35 sources are well represented by YSOs in early evolutionary phases, from Classes0 to ClassI. The Herschel detections also provided far-IR counterparts for 4 H2O masers and 11 objects previously known from mid-IR observations. The existence of so many embedded sources confirms the hypothesis of an intense and ongoing star formation activity in the region around NGC3603YC.

We present spectroscopic observations of FIR fine-structure lines of 26 Seyfert galaxies obtained with the Herschel-PACS spectrometer. These observations are complemented by spectroscopy with Spitzer-IRS and Herschel-SPIRE. The ratios of the OIII, NII, SIII and NeV lines have been used to determine electron densities in the ionised gas regions. The CI lines, observed with SPIRE, have been used to measure the densities in the neutral gas, while the OI lines provide a measure of the gas temperature, at densities below 10000 cm-3. Using the OI145/63um and SIII33/18um line ratios we find an anti-correlation of the temperature with the gas density. Using various fine-structure line ratios, we find that density stratification is common in these active galaxies. On average, the electron densities increase with the ionisation potential of the ions producing the NII, SIII and NeV emission. The infrared emission lines arise partly in the Narrow Line Region (NLR) photoionised by the AGN central engine, partly in HII regions photo ionised by hot stars and partly in neutral gas in photo-dissociated regions (PDRs). We attempt to separate the contributions to the line emission produced in these different regions by comparing our emission line ratios to empirical and theoretical values. In particular, we tried to separate the contribution of AGN and star formation by using a combination of Spitzer and Herschel lines, and we found that, besides the well known mid-IR line ratios, the mixed mid-IR/far-IR line ratio of OIII88um/OIV26um can reliably discriminate the two emission regimes, while the far-IR line ratio of CII157um/OI63um is only able to mildly separate the two regimes. By comparing the observed CII157um/NII205um ratio with photoionisation models, we also found that most of the CII emission in the galaxies we examined is due to PDRs.

Context: Herschel observations suggest that the H$_2$O distribution in outflows from low-mass stars resembles the H$_2$ emission. It is still unclear which of the different excitation components that characterise the mid- and near-IR H$_2$ distribution is associated with H$_2$O. Aim: The aim is to spectrally resolve the different excitation components observed in the H$_2$ emission. This will allow us to identify the H$_2$ counterpart associated with H$_2$O and finally derive directly an H$_2$O abundance estimate with respect to H$_2$. Methods: We present new high spectral resolution observations of H$_2$ 0-0 S(4), 0-0 S(9), and 1-0 S(1) towards HH 54, a bright nearby shock region in the southern sky. In addition, new Herschel-HIFI H$_2$O (2$_{12}$$-$1$_{01}$) observations at 1670~GHz are presented. Results: Our observations show for the first time a clear separation in velocity of the different H$_2$ lines: the 0-0 S(4) line at the lowest excitation peaks at $-$7~km~s$^{-1}$, while the more excited 0-0 S(9) and 1-0 S(1) lines peak at $-$15~km~s$^{-1}$. H$_2$O and high-$J$ CO appear to be associated with the H$_2$ 0-0 S(4) emission, which traces a gas component with a temperature of 700$-$1000 K. The H$_2$O abundance with respect to H$_2$ 0-0 S(4) is estimated to be $X$(H$_2$O)$<$1.4$\times$10$^{-5}$ in the shocked gas over an area of 13$^{\prime\prime}$. Conclusions: We resolve two distinct gas components associated with the HH 54 shock region at different velocities and excitations. This allows us to constrain the temperature of the H$_2$O emitting gas ($\leq$1000 K) and to derive correct estimates of H$_2$O abundance in the shocked gas, which is lower than what is expected from shock model predictions.

We present a study of the structure of the Galactic interstellar medium through the Delta-variance technique, related to the power spectrum and the fractal properties of infrared/sub-mm maps. Through this method, it is possible to provide quantitative parameters which are useful to characterize different morphological and physical conditions, and to better constrain the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey carried out at five photometric bands from 70 to 500 \mu m constitutes an unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2{\deg}x2{\deg} Hi-GAL tiles located in the third Galactic quadrant (217{\deg} < l < 225{\deg}, -2{\deg} < b < 0{\deg}). The low level of confusion along the line of sight testified by CO observations makes this region an ideal case. We find very different values of the power spectrum slope from tile to tile but also from wavelength to wavelength (2 < \beta < 3), with similarities between fields attributable to components located at the same distance. Thanks to the comparison with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed.

We present the data reduction pipeline for the Hi-GAL survey. Hi-GAL is a key project of the Herschel satellite which is mapping the inner part of the Galactic plane (|l| <= 70\cdot and |b| <= 1\cdot), using 2 PACS and 3 SPIRE frequency bands, from 70{\mu}m to 500{\mu}m. Our pipeline relies only partially on the Herschel Interactive Standard Environment (HIPE) and features several newly developed routines to perform data reduction, including accurate data culling, noise estimation and minimum variance map-making, the latter performed with the ROMAGAL algorithm, a deep modification of the ROMA code already tested on cosmological surveys. We discuss in depth the properties of the Hi-GAL Science Demonstration Phase (SDP) data.

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

Variations in the dust emissivity are critical for gas mass determinations derived from far-infrared observations, but also for separating dust foreground emission from the Cosmic Microwave Background (CMB). Hi-GAL observations allow us for the first time to study the dust emissivity variations in the inner regions of the Galactic plane at resolution below 1 degree. We present maps of the emissivity spectral index derived from the combined Herschel PACS 160 \mu m, SPIRE 250 \mu m, 350 \mu m, and 500 \mu m data, and the IRIS 100 \mu m data, and we analyze the spatial variations of the spectral index as a function of dust temperature and wavelength in the two Science Demonstration Phase Hi-GAL fields, centered at l=30{\deg} and l=59{\deg}. Applying two different methods, we determine both dust temperature and emissivity spectral index between 100 and 500 \mu m, at an angular resolution of 4'. Combining both fields, the results show variations of the emissivity spectral index in the range 1.8-2.6 for temperatures between 14 and 23 K. The median values of the spectral index are similar in both fields, i.e. 2.3 in the range 100-500 \mu m, while the median dust temperatures are equal to 19.1 K and 16.0 K in the l=30{\deg} and l=59{\deg} field, respectively. Statistically, we do not see any significant deviations in the spectra from a power law emissivity between 100 and 500 \mu m. We confirm the existence of an inverse correlation between the emissivity spectral index and dust temperature, found in previous analyses.

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.

We present a first study of the star-forming compact dust condensations revealed by Herschel in the two 2 \times 2 \degr Galactic Plane fields centered at [l;b] = [30\degr; 0 \degr] and [l;b] = [59\degr; 0 \degr], respectively, and observed during the Science Demonstration Phase for the Herschel infrared Galactic Plane survey (Hi-GAL) Key-Project. Compact source catalogs extracted for the two fields in the five Hi-GAL bands (70, 160, 250, 350 and 500 $\mu$m) were merged based on simple criteria of positional association and spectral energy distribution (SED) consistency into a final catalog which contains only coherent SEDs with counterparts in at least three adjacent Herschel bands. These final source lists contain 528 entries for the l = 30\degr field, and 444 entries for the l = 59\degr field. The SED coverage has been augmented with ancillary data at 24 $\mu$m and 1.1 mm. SED modeling for the subset of 318 and 101 sources (in the two fields, respectively) for which the distance is known was carried out using both a structured star/disk/envelope radiative transfer model and a simple isothermal grey-body. Global parameters like mass, luminosity, temperature and dust properties have been estimated. The Lbol/Menv ratio spans four orders of magnitudes from values compatible with the pre-protostellar phase to embedded massive zero-age main sequence stars. Sources in the l = 59\degr field have on average lower L/M, possibly outlining an overall earlier evolutionary stage with respect to the sources in the l = 30\degr field. Many of these cores are actively forming high-mass stars, although the estimated core surface densities appear to be an order of magnitude below the 1 g cm$^{-2}$ critical threshold for high-mass star formation.

Infrared dark clouds (IRDCs) are cold and dense reservoirs of gas potentially available to form stars. Many of these clouds are likely to be pristine structures representing the initial conditions for star formation. The study presented here aims to construct and analyze accurate column density and dust temperature maps of IRDCs by using the first Herschel data from the Hi-GAL galactic plane survey. These fundamental quantities, are essential for understanding processes such as fragmentation in the early stages of the formation of stars in molecular clouds. We have developed a simple pixel-by-pixel SED fitting method, which accounts for the background emission. By fitting a grey-body function at each position, we recover the spatial variations in both the dust column density and temperature within the IRDCs. This method is applied to a sample of 22 IRDCs exhibiting a range of angular sizes and peak column densities. Our analysis shows that the dust temperature decreases significantly within IRDCs, from background temperatures of 20-30 K to minimum temperatures of 8-15 K within the clouds, showing that dense molecular clouds are not isothermal. Temperature gradients have most likely an important impact on the fragmentation of IRDCs. Local temperature minima are strongly correlated with column density peaks, which in a few cases reach NH2 = 1 x 10^{23} cm^{-2}, identifying these clouds as candidate massive prestellar cores. Applying this technique to the full Hi-GAL data set will provide important constraints on the fragmentation and thermal properties of IRDCs, and help identify hundreds of massive prestellar core candidates.

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

We present a multi-epoch (20 years baseline) kinematical investigation of HH52, 53, and 54 at optical and near-IR wavelengths, along with medium and high- resolution spectroscopic analyses, probing the kinematical and physical time variability conditions of the gas along the flows. By means of multi-epoch and multi-wavelength narrow-band images, we derived proper motions, tangential velocities, velocity and flux variability of the knots. Radial velocities and physical parameters of the gas were derived from spectroscopy. Finally, spatial velocities and inclination of the flows were obtained by combining both imaging and spectroscopy. The P.M. analysis reveals three distinct, partially overlapping outflows. In 20 years, about 60% of the knots show some degree of flux variability. Our set of observations apparently indicates acceleration and deceleration in a variety of knots along the jets. For about 20% of the knots, mostly coincident with working surfaces or interacting knots along the flows, a relevant variability in both flux and velocity is observed. We argue that both variabilities are related and that all or part of the kinetic energy lost by the interacting knots is successively radiated. The analysis indicates the presence of very light, ionised, and hot flows, impacting a denser medium. Several knots are deflected. At least for a couple of them (HH54 G and G0), the deflection originates from the collision of the two. For the more massive parts of the flow, the deflection is likely the result of the flow collision with a dense cloud or with clumps.

We aim to characterise the morphology and the physical parameters governing the shock physics of the Herbig-Haro object HH99B. We have obtained SINFONI-SPIFFI IFU spectroscopy between 1.10 and 2.45 um detecting more than 170 emission lines, Most of them come from ro-vibrational transitions of H_2 and [FeII]. All the brightest lines appear resolved in velocity. Intensity ratios of ionic lines have been compared with predictions of NLTE models to derive bi-dimensional maps of extinction and electron density, along with estimates of temperature, fractional ionisation and atomic hydrogen post-shock density. H_2 line intensities have been interpreted in the framework of Boltzmann diagrams, from which we have derived maps of extinction and temperature of the molecular gas. From the intensity maps of bright lines the kinematical properties of the shock(s) at work in the region have been delineated. Finally, from selected [FeII] lines, constraints on the spontaneous emission coefficients of the 1.257, 1.321 and 1.644 um lines are provided. The kinematical properties derived for the molecular gas substantially confirm those published in Davis et al.(1999), while new information (e.g. v_shock ~115 km s^-1 is provided for the shock component responsible for the ionic emission. We also provide an indirect measure of the H_2 breakdown speed (between 70 and 90 km s^-1) and compute the inclination angle with respect to the line of sight. The map parameters, along with images of the observed line intensities, will be used to put stringent constraints on up-to-date shock models.

Aims. We present a detailed study of the infrared line emission (1-200 micron) in the Herbig-Haro object HH54. Our database comprises: high- (R~9000) and low- (R~600) resolution spectroscopic data in the near-infrared band (1-2.5 micron); mid-infrared spectrophotometric images (5-12 micron); and, far-IR (45-200 micron, R~200) spectra acquired with the ISO satellite. As a result, we provide the detection of and the absolute fluxes for more than 60 molecular features (mainly from H2 in the near- and mid-infrared and from H2O, CO and OH in the far-infrared) and 23 ionic lines. Methods. The H2 lines, coming from levels from v=0 to v=4 have been interpreted in the context of a state-of-art shock code, whose output parameters are adopted as input to a Large Velocity Gradient computation in order to interpret the FIR emission of CO, H2O and OH. Results. The H2 emission can be interpreted as originating in either steady-state J-type shocks or in quasi-steady J-type shocks with magnetic precursor. However, our multi-species analysis shows that only a model of a J-type shock with magnetic precursor (v_shock=18 km/s, n_H=10^4 cm^-3, B=100 microG, age=400 yr) can account for both the observed H2 emission and the CO and H2O lines. Such a model predicts a H2O abundance of ~ 7 10^-5, in agreement with estimations from other shock models for outflows associated with low mass protostars. We can exclude the possibility that the observed atomic lines arise in the same shock as the molecular lines, and give arguments in favour of the presence of a further high-velocity, fully dissociative shock component in the region. Finally, in view of the forthcoming spectroscopic facilities on board of the Herschel satellite, we provide predictions for H2O lines considered to be the most suitable for diagnostic purposes.

We present a near-IR study of the EXor variable V1118 Ori, performed by following a slightly declining phase after a recent outburst. In particular, the near-IR (0.8 - 2.3 micron) spectrum, obtained for the first time, shows a large variety of emission features of the HI and HeI recombination and CO overtone. By comparing the observed spectrum with a wind model, a mass loss rate value is derived along with other parameters whose values are typical of an accreting T Tauri star. In addition, we have used X-ray data from the XMM archive, taken in two different epochs during the declining phase monitored in IR. X-ray emission (in the range 0.5 - 10 keV) permits to derive several parameters which confirm the T Tauri nature of the source. In the near-IR the object maintains a low visual extinction during all the activity phases, confirming that variable extinction does not contribute to brightness variations. The lack of both a significant amount of circumstellar material and any evidence of IR cooling from collimated jet/outflow driven by the source, indicates that, at least this member of the EXor class, is in a late stage of the Pre-Main Sequence evolution. In the X-ray regime, an evident fading is present, detected in the post-outburst phase, that cannot be reconciled with the presence of any absorbing material. This circumstance, combined with the persistence (in the pre- and post-outburst phases) of a temperature component at about 10 MK, suggests that accretion has some influence in regulating the coronal activity.