We have surveyed a sample of massive star-forming regions located over a range of distances from the Galactic centre for methyl formate, HCOOCH$_{3}$, and its isotopologues H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$. The observations were carried out with the APEX telescope in the frequency range 283.4--287.4~GHz. Based on the APEX observations, we report tentative detections of the $^{13}$C-methyl formate isotopologue HCOO$^{13}$CH$_{3}$ towards the following four massive star-forming regions: Sgr~B2(N-LMH), NGC~6334~IRS~1, W51 e2 and G19.61-0.23. In addition, we have used the 1~mm ALMA science verification observations of Orion-KL and confirm the detection of the $^{13}$C-methyl formate species in Orion-KL and image its spatial distribution. Our analysis shows that the $^{12}$C/$^{13}$C isotope ratio in methyl formate toward Orion-KL Compact Ridge and Hot Core-SW components (68.4$\pm$10.1 and 71.4$\pm$7.8, respectively) are, for both the $^{13}$C-methyl formate isotopologues, commensurate with the average $^{12}$C/$^{13}$C ratio of CO derived toward Orion-KL. Likewise, regarding the other sources, our results are consistent with the $^{12}$C/$^{13}$C in CO. We also report the spectroscopic characterization, which includes a complete partition function, of the complex H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$ species. New spectroscopic data for both isotopomers H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$, presented in this study, has made it possible to measure this fundamentally important isotope ratio in a large organic molecule for the first time.

A fundamental question about the early evolution of low-mass protostars is when circumstellar disks may form. High angular resolution observations of molecular transitions in the (sub)millimeter wavelength windows make it possible to investigate the kinematics of the gas around newly-formed stars, for example to identify the presence of rotation and infall. IRAS 16293-2422 was observed with the extended Submillimeter Array (eSMA) resulting in subarcsecond resolution (0.46" x 0.29", i.e. $\sim$ 55 $\times$ 35~AU) images of compact emission from the C$^{17}$O (3-2) and C$^{34}$S (7-6) transitions at 337~GHz (0.89~mm). To recover the more extended emission we have combined the eSMA data with SMA observations of the same molecules. The emission of C$^{17}$O (3-2) and C$^{34}$S (7-6) both show a velocity gradient oriented along a northeast-southwest direction with respect to the continuum marking the location of one of the components of the binary, IRAS16293A. Our combined eSMA and SMA observations show that the velocity field on the 50--400~AU scales is consistent with a rotating structure. It cannot be explained by simple Keplerian rotation around a single point mass but rather needs to take into account the enclosed envelope mass at the radii where the observed lines are excited. We suggest that IRAS 16293-2422 could be among the best candidates to observe a pseudo-disk with future high angular resolution observations.

The fast outflow emerging from a region associated with massive star formation in the Orion Molecular Cloud 1 (OMC-1), located behind the Orion Nebula, appears to have been set in motion by an explosive event. Here we study the structure and dynamics of outflows in OMC-1. We combine radial velocity and proper motion data for near-IR emission of molecular hydrogen to obtain the first 3-dimensional (3D) structure of the OMC-1 outflow. Our work illustrates a new diagnostic tool for studies of star formation that will be exploited in the near future with the advent of high spatial resolution spectro-imaging in particular with data from the Atacama Large Millimeter Array (ALMA). We use published radial and proper motion velocities obtained from the shock-excited vibrational emission in the H2 v=1-0 S(1) line at 2.122 $\mu$m obtained with the GriF instrument on the Canada-France-Hawaii Telescope, the Apache Point Observatory, the Anglo-Australian Observatory and the Subaru Telescope. These data give the 3D velocity of ejecta yielding a 3D reconstruction of the outflows. This allows one to view the material from different vantage points in space giving considerable insight into the geometry. Our analysis indicates that the ejection occurred <720 years ago from a distorted ring-like structure of ~15" (6000 AU) in diameter centered on the proposed point of close encounter of the stars BN, source I and maybe also source n. We propose a simple model involving curvature of shock trajectories in magnetic fields through which the origin of the explosion and the centre defined by extrapolated proper motions of BN, I and n may be brought into spatial coincidence.

We present near-infrared imaging and spectroscopic high spatial resolution observations of the SMC region N88 containing the bright, excited, extincted and compact H II region N88A of size ~ 1 pc. To investigate its stellar content and reddening, N88 was observed using spectroscopy and imagery in the JHKs- and L'-band at a spatial resolution of ~ 0.1-0.3", using the VLT UT4 equipped with the NAOS adaptive optics system. In order to attempt to establish if the origin of the infra-red (IR) excess is due to bright nebulosity, circumstellar material and/or local dust, we used Ks vs J-K colour-magnitude (CM) and JHK colour-colour (CC) diagrams, as well as L' imagery.Our IR-data reveal in the N88 area an IR-excess fraction of geq 30 per cent of the detected stars,as well as an unprecedently detailed morphology of N88A. It consists of an embedded cluster of ~3.5" (~ 1 pc) in diameter, of at least thirteen resolved stars superposed with an unusual bright continuum centered on a very bright star. The four brightest stars in this cluster lie red-ward of H-K geq 0.45 mag, and could be classified as young stellar object (YSO) candidates. Four other probable YSO candidates are also detected in N88 along a south-north bow-shaped thin H2 filament at ~ 7" east of the young central bright star. At 0.2" east of this star, a heavily embedded core is detected that could be a massive class I protostar candidate. The 2.12 mu H2 image of N88A resembles a shell of diameter ~ 3" ~ 0.9 pc) centered on the bright star. The line ratios of H2 2-1 S(1) and 1-0 S(0) relative to 1-0 S(1), as well as the presence of high v lines, are indicative of photodissociation regions, rather than shocks.

Context. In order to understand the evolution of molecular clouds it is important to identify the departures from self-similarity associated with the scales of self-gravity and the driving of turbulence. Aims. A method is described based on structure functions for determining whether a region of gas, such as a molecular cloud, is fractal or contains structure with characteristic scale sizes. Methods. Using artificial data containing structure it is shown that derivatives of higher order structure functions provide a powerful way to detect the presence of characteristic scales should any be present and to estimate the size of such structures. The method is applied to observations of hot H2 in the Kleinman-Low nebula, north of the Trapezium stars in the Orion Molecular Cloud, including both brightness and velocity data. The method is compared with other techniques such as Fourier transform and histogram techniques. Results. It is found that the density structure, represented by H2 emission brightness in the K-band (2-2.5micron), exhibits mean characteristic sizes of 110, 550, 1700 and 2700AU. The velocity data show the presence of structure at 140, 1500 and 3500AU. Compared with other techniques such as Fourier transform or histogram, the method appears both more sensitive to characteristic scales and easier to interpret.

Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions. Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. These are compared with the predictions of two-dimensional projections of simulations of supersonic hydrodynamic turbulence. Results: The structure functions of OMC1 are not well represented by power laws, but show clear deviations below 2000AU. However, using the technique of extended self-similarity, power laws are recovered at scales down to 160AU. The scaling of the higher order structure functions with order deviates from the standard scaling for supersonic turbulence. This is explained as a selection effect of preferentially observing the shocked part of the gas and the scaling can be reproduced using line-of-sight integrated velocity data from subsets of supersonic turbulence simulations. These subsets select regions of strong flow convergence and high density associated with shock structure. Deviations of the structure functions in OMC1 from power laws cannot however be reproduced in simulations and remains an outstanding issue.

Observations are reported of H2 IR emission in the S(1) v=1-0 line at 2.121 microns in the Orion Molecular Cloud, OMC1, using the GriF instrument on the Canada-France-Hawaii Telescope. GriF is a combination of adaptive optics and Fabry-Perot interferometry, yielding a spatial resolution of 0.15" to 0.18" and a velocity discrimination as high as 1 km/s. Thanks to the high spatial and velocity resolution of the GriF data, 193 bright H2 emission regions can be identified in OMC1. The general characteristics of these features are described in terms of radial velocities, brightness and spatial displacement of maxima of velocity and brightness, the latter to yield the orientation of flows in the plane of the sky. Strong spatial correlation between velocity and bright H2 emission is found and serves to identify many features as shocks. Important results are: (i) velocities of the excited gas illustrate the presence of a zone to the south of BN-IRc2 and Peak 1, and the west of Peak 2, where there is a powerful blue-shifted outflow with an average velocity of -18 km/s. This is shown to be the NIR counterpart of an outflow identified in the radio from source I, a very young O-star. (ii) There is a band of weak velocity features (<5 km/s) in Peak 1 which may share a common origin through an explosive event, in the BN-IRc2 region, with the fast-moving fingers (or bullets) to the NW of OMC1. (iii) A proportion of the flows are likely to represent sites of low mass star formation and several regions show multiple outflows, probably indicative of multiple star formation within OMC1. The high spatial and velocity resolution of the GriF data show these and other features in more detail than has previously been possible.

A statistical analysis is presented of the turbulent velocity structure in the Orion Molecular Cloud at scales ranging from 70 AU to 30000 AU. Results are based on IR Fabry-Perot interferometric observations of shock and photon-excited H2 in the K-band S(1) v=1-0 line at 2.121micron and refer to the dynamical characteristics of warm perturbed gas. Observations establish that the Larson size-linewidth relation is obeyed to the smallest scales studied here extending the range of validity of this relationship by nearly 2 orders of magnitude. The velocity probability distribution function (PDF) is constructed showing extended exponential wings, providing evidence of intermittency, further supported by the skewness and kurtosis of the velocity distribution. Variance and kurtosis of the PDF of velocity differences are constructed as a function of lag. The variance shows an approximate power law dependence on lag, with exponent significantly lower than the Kolmogorov value, and with deviations below 2000AU which are attributed to outflows and possibly disk structures associated with low mass star formation within OMC1. The kurtosis shows strong deviation from a gaussian velocity field, providing evidence of velocity correlations at small lags. Results agree accurately with semi-empirical simulations in Eggers & Wang (1998). In addition, 170 individual H2 emitting clumps have been analysed with sizes between 500 and 2200 AU. These show considerable diversity with regard to PDFs and variance functions. Our analysis constitutes the first characterization of the turbulent velocity field at the scale of star formation and provide a dataset which models of star-forming regions should aim to reproduce.

We present MERLIN observations of OH maser and radio continuum emission from the Ultra Luminous IR Galaxy Markarian 231. The 1665- and 1667-MHz transitions have a combined velocity extent of 720 km/s and show a similar position-velocity structure including a gradient of 1.7 km/s/pc from NW to SE along the 420-pc major axis, steeper in the inner few tens of pc. The maser distribution is modelled as a torus rotating about an axis inclined at ~45deg. We estimate the enclosed mass density to be 320(90) Msun in a flattened distribution, including a central unresolved mass of </=8E+06 Msun. All the maser emission is projected against a region with a radio continuum brightness temperature >/=1E+05 K, giving a maser gain of </=2.2. The 1667:1665-MHz line ratio is close to the LTE ratio of 1.8 consistent with radiatively pumped, unsaturated masers. The size of individual masing regions is in the range 0.25-4 pc with a covering factor close to unity. There are no very bright compact masers, in contrast to galaxies such as the Seyfert 2 Markarian 273 where the masing torus is viewed nearer edge-on. The comparatively modest maser amplification seen from Markarian 231 is consistent with its classification as a Seyfert 1. Most of the radio continuum emission on 50-500 pc scales is probably of starburst origin but the compact peak is 0.4 per cent polarized by a magnetic field running north-south, similar to the jet direction on these scales. There is no close correlation between maser and continuum intensity. Comparisons with other data show that the jet changes direction close the nucleus and suggest that the sub-kpc disc hosting the masers and starburst activity is severely warped.

Near-infrared observations of line emission from excited H2 and in the continuum are reported in the direction of the Orion molecular cloud OMC1, using the European Southern Observatory Very Large Telescope UT4, equipped with the NAOS adaptive optics system and the CONICA infrared array camera. Spatial resolution has been achieved at close to the diffraction limit of the telescope (0.08" - 0.12") and images show a wealth of morphological detail. Structure is not fractal but shows two preferred scale sizes of 2.4" (1100 AU) and 1.2" (540 AU), where the larger scale may be associated with star formation.

A stellar hydrodynamic pulsation model has been combined with a SiO maser model in an attempt to calculate the temporal variability of SiO maser emission in the circumstellar envelope (CE) of a model AGB star. This study investigates whether the variations in local physical conditions brought about by shocks are the predominant contributing factor to SiO maser variability because, in this work, the radiative part of the pump is constant. We find that some aspects of the variability are not consistent with a pump provided by shock-enhanced collisions alone. In these simulations, gas parcels of relatively enhanced SiO abundance are distributed in a model CE by a Monte Carlo method, at a single epoch of the stellar cycle. From this epoch on, Lagrangian motions of individual parcels are calculated according to the velocity fields encountered in the model CE during the stellar pulsation cycle. The potentially masing gas parcels therefore experience different densities and temperatures, and have varying line-of-sight velocity gradients throughout the stellar cycle, which may or may not be suitable to produce maser emission. At each epoch (separated by 16.6 days), emission lines from the parcels are combined to produce synthetic spectra and VLBI-type images. We report here the results for v=1, J=1-0 (43-GHz) and J=2-1 (86-GHz) masers.