Synthetic observations are playing an increasingly important role across astrophysics, both for interpreting real observations and also for making meaningful predictions from models. In this review, we provide an overview of methods and tools used for generating, manipulating and analysing synthetic observations and their application to problems involving star formation and the interstellar medium. We also discuss some possible directions for future research using synthetic observations.

When modeling astronomical objects throughout the universe, it is important to correctly treat the limitations of the data, for instance finite resolution and sensitivity. In order to simulate these effects, and to make radiative transfer models directly comparable to real observations, we have developed an open-source Python package called the FluxCompensator that enables the post-processing of the output of 3-d Monte-Carlo radiative transfer codes, such as HYPERION. With the FluxCompensator, realistic synthetic observations can be generated by modelling the effects of convolution with arbitrary point-spread functions (PSFs), transmission curves, finite pixel resolution, noise and reddening. Pipelines can be applied to compute synthetic observations that simulate observatories, such as the Spitzer Space Telescope or the Herschel Space Observatory. Additionally, this tool can read in existing observations (e.g. FITS format) and use the same settings for the synthetic observations. In this paper, we describe the package as well as present examples of such synthetic observations.

Through an extensive set of realistic synthetic observations (produced in Paper I), we assess in this part of the paper series (Paper III) how the choice of observational techniques affects the measurement of star-formation rates (SFRs) in star-forming regions. We test the accuracy of commonly used techniques and construct new methods to extract the SFR, so that these findings can be applied to measure the SFR in real regions throughout the Milky Way. We investigate diffuse infrared SFR tracers such as those using 24 {\mu}m, 70 {\mu}m and total infrared emission, which have been previously calibrated for global galaxy scales. We set up a toy model of a galaxy and show that the infrared emission is consistent with the intrinsic SFR using extra-galactic calibrated laws (although the consistency does not prove their reliability). For local scales, we show that these techniques produce completely unreliable results for single star-forming regions, which are governed by different characteristic timescales. We show how calibration of these techniques can be improved for single star-forming regions by adjusting the characteristic timescale and the scaling factor and give suggestions of new calibrations of the diffuse star-formation tracers. We show that star-forming regions that are dominated by high-mass stellar feedback experience a rapid drop in infrared emission once high-mass stellar feedback is turned on, which implies different characteristic timescales. Moreover, we explore the measured SFRs calculated directly from the observed young stellar population. We find that the measured point sources follow the evolutionary pace of star formation more directly than diffuse star-formation tracers.

We use a large data-set of realistic synthetic observations (PaperI) to assess how observational techniques affect the measurement of physical properties of star-forming regions. In this paper (PaperII), we explore the reliability of the measured total gas mass, dust surface density and dust temperature maps derived from modified blackbody fitting of synthetic Herschel observations. We found from our pixel-by-pixel analysis of the measured dust surface density and dust temperature a worrisome error spread especially close to star-formation sites and low-density regions, where for those "contaminated" pixels the surface densities can be under/overestimated by up to three orders of magnitude. In light of this, we recommend to treat the pixel-based results from this technique with caution in regions with active star formation. In regions of high background typical in the inner Galactic plane, we are not able to recover reliable surface density maps of individual synthetic regions, since low-mass regions are lost in the FIR background. When measuring the total gas mass of regions in moderate background, we find that modified blackbody fitting works well (absolute error:+9%;-13%) up to 10kpc distance (errors increase with distance). Commonly, the initial images are convolved to the largest common beam-size, which smears contaminated pixels over large areas. The resulting information loss makes this commonly-used technique less verifiable as now chi^2-values cannot be used as a quality indicator of a fitted pixel. Our control measurements of the total gas mass (without the step of convolution to the largest common beam size) produce similar results (absolute error:+20%;-7%) while having much lower median errors especially for the high-mass stellar feedback phase. In upcoming papers (III&IV) we test the reliability of measured star-formation rate with direct and indirect techniques.

Through synthetic observations of a hydrodynamical simulation of an evolving star-forming region, we assess how the choice of observational techniques affects the measurements of properties which trace star formation. Testing and calibrating observational measurements requires synthetic observations which are as realistic as possible. In this part of the paper series (Paper I), we explore different techniques for how to map the distributions of densities and temperatures from the particle-based simulations onto a Voronoi mesh suitable for radiative transfer and consequently explore their accuracy. We further test different ways to set up the radiative transfer in order to produce realistic synthetic observations. We give a detailed description of all methods and ultimately recommend techniques. We have found that the flux around 20 microns is strongly overestimated when blindly coupling the dust radiative transfer temperature with the hydrodynamical gas temperature. We find that when instead assuming a constant background dust temperature in addition to the radiative transfer heating, the recovered flux is consistent with actual observations. We present around 5800 realistic synthetic observations for Spitzer and Herschel bands, at different evolutionary time-steps, distances and orientations. In the upcoming papers of this series (Paper II, Paper III and Paper IV), we will test and calibrate measurements of the star-formation rate (SFR), gas mass and the star-formation efficiency (SFE) using our realistic synthetic observations.

Filaments represent a key structure during the early stages of the star formation process. Simulations show filamentary structure commonly formed before and during the formation of cores. Aims. The Serpens Core represents an ideal laboratory to test the state-of-the-art of simulations of turbulent Giant Molecular Clouds. We use Herschel observations of the Serpens Core to compute temperature and column density maps of the region. Among the simulations of Dale et al. (2012), we select the early stages of their Run I, before stellar feedback is initiated, with similar total mass and physical size as the Serpens Core. We derive temperature and column density maps also from the simulations. The observed distribution of column densities of the filaments has been analysed first including and then masking the cores. The same analysis has been performed on the simulations as well. A radial network of filaments has been detected in the Serpens Core. The analysed simulation shows a striking morphological resemblance to the observed structures. The column density distribution of simulated filaments without cores shows only a log-normal distribution, while the observed filaments show a power-law tail. The power-law tail becomes evident in the simulation if one focuses just on the column density distribution of the cores. In contrast, the observed cores show a flat distribution. Even though the simulated and observed filaments are subjectively similar-looking, we find that they behave in very different ways. The simulated filaments are turbulence-dominated regions, the observed filaments are instead self-gravitating structures that will probably fragment into cores.

By means of radiative transfer simulation we study the evolution of the far-infrared colours of protoplanetary discs undergoing inside-out dispersal, often referred to as transition discs. We show that a brightening of the mid and far-infrared emission from these objects is a natural consequence of the removal of the inner disc. Our results can fully explain recent observations of transition discs in the Chamaleon and Lupus star forming regions from the Herschel Gould Belt Survey, which show a higher median for the 70?um (Herschel PACS 1) band of known transition objects compared with primordial discs. Our theoretical results hence support the suggestion that the 70?um band may be a powerful diagnostic for the identification of transition discs from photometry data, provided that the inner hole is larger than tens of AU, depending on spectral type. Furthermore we show that a comparison of photometry in the K , 12?um and 7u0?m bands to model tracks can provide a rough, but quick estimate of the inner hole size of these objects, provided their inclination is below ?85 degrees and the inner hole size is again larger than tens of AU.

In contrast to most other galaxies, star-formation rates in the Milky Way can be estimated directly from Young Stellar Objects (YSOs). In the Central Molecular Zone (CMZ) the star-formation rate calculated from the number of YSOs with 24 microns emission is up to order of magnitude higher than the value estimated from methods based on diffuse emission (such as free-free emission). Whether this effect is real or whether it indicates problems with either or both star formation rate measures is not currently known. In this paper, we investigate whether estimates based on YSOs could be heavily contaminated by more evolved objects such as main-sequence stars. We present radiative transfer models of YSOs and of main-sequence stars in a constant ambient medium which show that the main-sequence objects can indeed mimic YSOs at 24 microns. However, we show that in some cases the main-sequence models can be marginally resolved at 24 microns, whereas the YSO models are always unresolved. Based on the fraction of resolved MIPS 24 microns sources in the sample of YSOs previously used to compute the star formation rate, we estimate the fraction of misclassified "YSOs" to be at least 63%, which suggests that the star-formation rate previously determined from YSOs is likely to be at least a factor of three too high.

The time-scale over which and the modality by which young stellar objects (YSOs) disperse their circumstellar discs dramatically influences the eventual formation and evolution of planetary systems. By means of extensive radiative transfer (RT) modelling, we have developed a new set of diagnostic diagrams in the infrared colour-colour plane (K-[24] vs. K-[8]), to aid with the classification of the evolutionary stage of YSOs from photometric observations. Our diagrams allow the differentiation of sources with unevolved (primordial) discs from those evolving according to different clearing scenarios (e.g. homologous depletion vs. inside-out dispersal), as well as from sources that have already lost their disc. Classification of over 1500 sources in 15 nearby star-forming regions reveals that approximately 39 % of the sources lie in the primordial disc region, whereas between 31 % and 32 % disperse from the inside-out and up to 22 % of the sources have already lost their disc. Less than 2 % of the objects in our sample lie in the homogeneous draining regime. Time-scales for the transition phase are estimated to be typically a few 10^5 years independent of stellar mass. Therefore, regardless of spectral type, we conclude that currently available infrared photometric surveys point to fast (of order 10 % of the global disc lifetime) inside-out clearing as the preferred mode of disc dispersal.

The time-scale over which and modality by which young stellar objects (YSOs) disperse their circumstellar discs dramatically influences the eventual formation and evolution of planetary systems. By means of extensive radiative transfer (RT) modelling, we have developed a new set of diagnostic diagrams in the infrared colour-colour plane (K-[24] vs. K-[8]), to aid with the classiffication of the evolutionary stage of YSOs from photometric observations. Our diagrams allow the differentiation of sources with un-evolved (primordial) discs from those evolving according to different clearing scenarios (e.g. homologous depletion vs. inside-out dispersal), as well as from sources that have already lost their disc. Classification of over 1500 sources in 15 nearby star-forming regions reveals that approximately 39% of the sources lie in the primordial disc region, whereas between 31% and 32% disperse from the inside-out and up to 22% of the sources have already lost their disc. Less than 2% of the objects in our sample lie in the homogeneous draining regime. Time-scales for the transition phase are estimated to be typically a few 10^5 years independent of stellar mass. Therefore, regardless of spectral type, we conclude that currently available infrared photometric surveys point to fast (of order 10% of the global disc lifetime) inside-out clearing as the preferred mode of disc dispersal.