We present a review of the TORUS radiation transfer and hydrodynamics code. TORUS uses a 1-D, 2-D or3-D adaptive mesh refinement scheme to store and manipulate the state variables, and solves the equation of radiative transfer using Monte Carlo techniques. A framework of microphysics modules is described, including atomic and molecular line transport in moving media, dust radiative equilibrium, photoionisation equilibrium, and time-dependent radiative transfer. These modules provide a flexible scheme for producing synthetic observations, either from analytical models or as post-processing of hydrodynamical simulations (both grid-based and Lagrangian). A hydrodynamics module is also presented, which maybe used in combination with the radiation-transport modules to perform radiation-hydrodynamics simulations. Benchmarking and validation tests of each major mode of operation are detailed, along with descriptions and performance/scaling tests of the various parallelisation schemes. We give examples of the uses of the code in the literature, including applications to low- and high-mass star formation, cluster feedback, and stellar winds, along with an Appendix listing the refereed papers that have used TORUS.

We simulate a self-gravitating, turbulent cloud of 1000 Msol with photoionization and radiation pressure feedback from a 34 Msol star. We use a detailed Monte Carlo radiative transfer scheme alongside the hydrodynamics to compute photoionization and thermal equilibrium with dust grains and multiple atomic species. Using these gas temperatures, dust temperatures, and ionization fractions, we produce self-consistent synthetic observations of line and continuum emission. We find that all material is dispersed from the (15.5 pc)$^3$ grid within 1.6 Myr or 0.74 free-fall times. Mass exits with a peak flux of $2 \times 10^{-3}$ Msol/yr, showing efficient gas dispersal. The model without radiation pressure has a slight delay in the breakthrough of ionization, but overall its effects are negligible. 85 per cent of the volume, and 40 per cent of the mass, become ionized -- dense filaments resist ionization and are swept up into spherical cores with pillars that point radially away from the ionizing star. We use free-free emission at 20 cm to estimate the production rate of ionizing photons. This is almost always underestimated: by a factor of a few at early stages, then by orders of magnitude as mass leaves the volume. We also test the ratio of dust continuum surface brightnesses at 450 and 850 micron to probe dust temperatures. This underestimates the actual temperature by more than a factor of 2 in areas of low column density or high line-of-sight temperature dispersion; the HII region cavity is particularly prone to this discrepancy. However, the probe is accurate in dense locations such as filaments.

We present a numerical simulation of the formation of a massive star using Monte- Carlo-based radiation hydrodynamics (RHD). The star forms via stochastic disc accretion and produces fast, radiation-driven bipolar cavities. We find that the evolution of the infall rate (considered to be the mass flux across a 1500 au spherical boundary), and the accretion rate onto the protostar, are broadly consistent with observational constraints. After 35kyr the star has a mass of 25 solar masses and is surrounded by a disc of mass 7 solar masses and 1500 au radius, and we find that the velocity field of the disc is close to Keplerian. Once again these results are consistent with those from recent high-resolution studies of discs around forming massive stars. Synthetic imaging of the RHD model shows good agreement with observations in the near- and far-IR, but may be in conflict with observations that suggests that MYSOs are typically circularly symmetric on the sky at 24.5 microns. Molecular line simulations of a CH3CN transition compare well with observations in terms of surface brightness and line width, and indicate that it should be possible to reliably extract the protostellar mass from such observations.

We present an extensive search for RR Lyrae stars in and around the ultra-faint Milky Way companions Segue 2 and Segue 3. The former (M_V = -2.5, Belokurov et al. 2009) appears to be an extremely faint dwarf galaxy companion of the Milky Way. The latter (M_V = 0.0, Fadely et al. 2011) is among the faintest star clusters known. We use B and V band time-series imaging obtained at the WIYN 0.9 meter telescope at Kitt Peak National Observatory to search for RR Lyrae in these objects. In our Segue 2 observations, we present a previously unknown fundamental mode (RRab) RR Lyrae star with a period of P_ab = 0.748 days. With this measurement, we revisit the inverse correlation between <P_ab> and <[Fe/H]> established in the literature for Milky Way dwarf galaxies and their RR Lyrae. In this context, the long period of Segue 2's RRab star as well as the known significant spread in metallicity in this dwarf galaxy are consistent with the observed trend in <P_ab> and <[Fe/H]>. We derive the first robust distance to Segue 2, using both its RRab star and spectroscopically confirmed blue horizontal branch stars. Using [Fe/H] = -2.16 and -2.44 dex, we find d_RRL = 36.6 +2.5/-2.4 and 37.7 +/- 2.7 kpc; assuming [Fe/H] = -2.257 dex, we find d_BHB = 34.4 +/- 2.6 kpc. Although no RR Lyrae were present in the Segue 3 field, we found a candidate eclipsing binary star system.

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.