We present models of the inner region of the circumstellar disk of RY Tau which aim to explain our near-infrared ($K$-band: $2.1\,\mu$m) interferometric observations while remaining consistent with the optical to near-infrared portions of the spectral energy distribution. Our sub-milliarcsecond resolution CHARA Array observations are supplemented with shorter baseline, archival data from PTI, KI and VLTI/GRAVITY and modeled using an axisymmetric Monte Carlo radiative transfer code. The $K$-band visibilities are well-fit by models incorporating a central star illuminating a disk with an inner edge shaped by dust sublimation at $0.210\pm0.005\,$au, assuming a viewing geometry adopted from millimeter interferometry ($65^{\circ}$ inclined with a disk major axis position angle of $23^{\circ}$). This sublimation radius is consistent with that expected of Silicate grains with a maximum size of $0.36-0.40\,\mu$m contributing to the opacity and is an order of magnitude further from the star than the theoretical magnetospheric truncation radius. The visibilities on the longest baselines probed by CHARA indicate that we lack a clear line-of-sight to the stellar photosphere. Instead, our analysis shows that the central star is occulted by the disk surface layers close to the sublimation rim. While we do not see direct evidence of temporal variability in our multi-epoch CHARA observations, we suggest the aperiodic photometric variability of RY~Tau is likely related temporal and/or azimuthal variations in the structure of the disk surface layers.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

We investigate radiative feedback from a 34 M$_\odot$ star in a $10^4$ M$_\odot$ turbulent cloud using three-dimensional radiation-hydrodynamics (RHD) models. We use Monte Carlo radiative transfer to accurately compute photoionization equilibrium and radiation pressure, with multiple atomic species and silicate dust grains. We include the diffuse radiation field, dust absorption/re-emission, and scattering. The cloud is efficiently dispersed, with 75 per cent of the mass leaving the (32.3 pc)$^3$ grid within 4.3 Myr (1.1 $t_{ff}$). This compares to all mass exiting within 1.6 Myr (0.74 $t_{ff}$) in our previously published $10^3$ M$_\odot$ cloud. At most 20 per cent of the mass is ionized, compared to 40 per cent in the lower mass model, despite the ionized volume fraction being 80 per cent in both, implying the higher mass cloud is more resilient to feedback. The total Jeans-unstable mass increases linearly up to 1500 M$_\odot$ before plateauing after 2 Myr, corresponding to a core formation efficiency of 15 per cent. We also measure the time-variation of the far-ultraviolet (FUV) radiation field, $G_0$, impinging on other cluster members, taking into account for the first time how this changes in a dynamic cluster environment with intervening opacity sources and stellar motions. Many objects remain shielded in the first 0.5 Myr whilst the massive star is embedded, after which $G_0$ increases by orders of magnitude. Gas motions later on cause comparable drops which happen instantaneously and last for $\sim$ 1 Myr before being restored. This highly variable UV field will influence the photoevaporation of protoplanetary discs near massive stars.

We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_\mathrm{disc}/M_*\gtrsim 0.25$ and $R_\mathrm{outer}\lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $\sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales

The first stable object to develop in the low-mass star formation process has long been predicted to be the first hydrostatic core (FHSC). Despite much effort, it has still yet to be definitively observed in nature. More specific observational signatures are required to enable observers to distinguish the FHSC from young, faint, but more evolved protostars. Here we present synthetic spectral line observations for CO, SO, CS and HCO$^+$ that were calculated from radiation (magneto)hydrodynamical models, chemical modelling and Monte Carlo radiative transfer. HCO$^+$ $(1-0)$ and SO $(8_7 - 7_6)$ spectra of the FHSC show variations for observations at a low inclination which may allow a candidate FHSC to be distinguished from a more evolved object. We find that the FHSC outflow is unlikely to be detectable with ALMA, which would discount the observed sources with slow outflows that are currently identified as candidate FHSCs. We compare the results of simulated ALMA observations with observed candidate FHSCs and recommend Oph A SM1N and N6-mm as the most promising candidates to follow up.

T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer 3 x higher resolution and better baseline position angle coverage over previous observations. We investigate the characteristics of circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The CHARA array offers opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB 3-telescope combiner in the K-band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and geometric modeling. Additionally, the fitting of radiative transfer models constrains the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Image reconstruction reveals an inclined disk with slight asymmetry along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission/scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17+/-0.02 au with an inclination of 50.9+/-1.0 degrees and minor axis position angle 60.8+/-1.2 degrees. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 um and a scale height of 15.0 au at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K-band excess by introducing dust above the mid-plane.

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 present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.

V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to identify the origin of the previously detected asymmetries. We have observed V1247 Ori at three epochs spanning $\sim678$ days using sparse aperture masking interferometry with Keck/NIRC2 and VLT/NACO. In addition, we search for signs of accretion through VLT/SPHERE-ZIMPOL spectral differential imaging in H$\alpha$ and R-band continuum. Our SMA sub-millimetre interferometry in 880 $\mu$m continuum and in the CO(3-2) line allows us to constrain the orientation and direction of rotation of the outer disc. We find the L'-band emission to be dominated by static features which trace forward-scattered dust emission from the inner edge of the outer disc located to the north-east. In H- and K-band, we see evidence for a companion candidate that moved systematically by 45$^{\circ}$ within the first $\sim$345 days. The separation of the companion candidate is not well constrained, but the observed position angle change is consistent with Keplerian motion of a body located on a 6 au orbit. From the SMA CO moment map, the location of the disc rim, and the detected orbital motion, we deduced the three-dimensional orientation of the disc. We see no indication of accretion in H$\alpha$ and set upper limits for an accreting companion. The measured contrast of the companion candidate in H and K is consistent with an actively accreting protoplanet. Hence, we identify V1247 Ori as a unique laboratory for studying companion-disc interactions and disc clearing.

We present high resolution ($\sim$300 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment of MM 1 in 1.3 mm continuum emission and CH$_{3}$CN emission for the first time. The object divides into two main sources - MM 1a, which is the source of a bipolar molecular outflow, and MM 1b, located 0.57'' (1920 au) to the South-East. The main component of MM 1a is an elongated continuum structure, perpendicular to the bipolar outflow, with a size of $0.141'' \times 0.050''$ ($480\times170$ au). The gas kinematics toward MM 1a probed via CH$_{3}$CN trace a variety of scales. The lower energy $J=12-11$ $K=3$ line traces extended, rotating gas within the outflow cavity, while the $v$8=1 line shows a clearly-resolved Keplerian rotation signature. Analysis of the gas kinematics and dust emission shows that the total enclosed mass in MM 1a is $40\pm5$ M$_{\odot}$ (where between 2.2-5.8 M$_{\odot}$ is attributed to the disk), while MM 1b is $<0.6$ M$_{\odot}$. The extreme mass ratio and orbital properties of MM 1a and MM 1b suggest that MM 1b is one of the first observed examples of the formation of a binary star via disk fragmentation around a massive young (proto)star.

We present comprehensive models of Herbig Ae star, HD 142666, which aim to simultaneously explain its spectral energy distribution (SED) and near-infrared (NIR) interferometry. Our new sub-milliarcsecond resolution CHARA (CLASSIC and CLIMB) interferometric observations, supplemented with archival shorter baseline data from VLTI/PIONIER and the Keck Interferometer, are modeled using centro-symmetric geometric models and an axisymmetric radiative transfer code. CHARA's 330 m baselines enable us to place strong constraints on the viewing geometry, revealing a disk inclined at 58 degrees from face-on with a 160 degree major axis position angle. Disk models imposing vertical hydrostatic equilibrium provide poor fits to the SED. Models accounting for disk scale height inflation, possibly induced by turbulence associated with magneto-rotational instabilities, and invoking grain growth to >1 micron size in the disk rim are required to simultaneously reproduce the SED and measured visibility profile. However, visibility residuals for our best model fits to the SED indicate the presence of unexplained NIR emission, particularly along the apparent disk minor axis, while closure phase residuals indicate a more centro-symmetric emitting region. In addition, our inferred 58 degree disk inclination is inconsistent with a disk-based origin for the UX Ori-type variability exhibited by HD 142666. Additional complexity, unaccounted for in our models, is clearly present in the NIR-emitting region. We propose the disk is likely inclined toward a more edge-on orientation and/or an optically thick outflow component also contributes to the NIR circumstellar flux.

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.

High angular resolution observations of young stellar objects are required to study the inner astronomical units of protoplanetary disks in which the majority of planets form. As they evolve, gaps open up in the inner disk regions and the disks are fully dispersed within ~10 Myrs. MWC 614 is a pre-transitional object with a ~10au radius gap. We present a set of high angular resolution observations of this object including SPHERE/ZIMPOL polarimetric and coronagraphic images in the visible, KECK/NIRC2 near-infrared aperture masking observations and VLTI (AMBER, MIDI, and PIONIER) and CHARA (CLASSIC and CLIMB) long-baseline interferometry at infrared wavelengths. We find that all the observations are compatible with an inclined disk (i ~55deg at a position angle of ~20-30deg). The mid-infrared dataset confirms the disk inner rim to be at 12.3+/-0.4 au from the central star. We determined an upper mass limit of 0.34 Msun for a companion inside the cavity. Within the cavity, the near-infrared emission, usually associated with the dust sublimation region, is unusually extended (~10 au, 30 times larger than the theoretical sublimation radius) and indicates a high dust temperature (T~1800 K). As a possible result of companion-induced dust segregation, quantum heated dust grains could explain the extended near-infrared emission with this high temperature. Our observations confirm the peculiar state of this object where the inner disk has already been accreted onto the star exposing small particles inside the cavity to direct stellar radiation.

A recent ALMA observation of the Elias 2-27 system revealed a two-armed structure extending out to ~300 au in radius. The protostellar disc surrounding the central star is unusually massive, raising the possibility that the system is gravitationally unstable. Recent work has shown that the observed morphology of the system can be explained by disc self-gravity, so we examine the physical properties of the disc necessary to detect self-gravitating spiral waves. Using three-dimensional Smoothed Particle Hydrodynamics, coupled with radiative transfer and synthetic ALMA imaging, we find that observable spiral structure can only be explained by self-gravity if the disc has a low opacity (and therefore efficient cooling), and is minimally supported by external irradiation. This corresponds to a very narrow region of parameter space, suggesting that, although it is possible for the spiral structure to be due to disc self-gravity, other explanations, such as an external perturbation, may be preferred.

Recent observational work has indicated that mechanisms for accretion and outflow in Herbig Ae/Be star-disk systems may differ from magnetospheric accretion (MA) as it is thought to occur in T Tauri star-disk systems. In this work, we assess the temporal evolution of spectral lines probing accretion and mass loss in Herbig Ae/Be systems and test for consistency with the MA paradigm. For two Herbig Ae/Be stars, HD 98922 (B9e) and V1295 Aql (A2e), we have gathered multi-epoch (~years) and high-cadence (~minutes) high-resolution optical spectra to probe a wide range of kinematic processes. Employing a line equivalent width evolution correlation metric introduced here, we identify species co-evolving (indicative of common line origin) via novel visualization. We interferometrically constrain often problematically degenerate parameters, inclination and inner disk radius, allowing us to focus on the structure of the wind, magnetosphere, and inner gaseous disk in radiative transfer models. Over all timescales sampled, the strongest variability occurs within the blueshifted absorption components of the Balmer series lines; the strength of variability increases with the cadence of the observations. Finally, high-resolution spectra allow us to probe substructure within the Balmer series' blueshifted absorption components: we observe static, low-velocity features and time-evolving features at higher velocities. Overall, we find the observed line morphologies and variability are inconsistent with a scaled-up T Tauri MA scenario. We suggest that as magnetic field structure and strength change dramatically with increasing stellar mass from T Tauri to Herbig Ae/Be stars, so too may accretion and outflow processes.

Herbig Ae/Be stars are young, pre-main-sequence stars that sample the transition in structure and evolution between low- and high-mass stars, providing a key test of accretion processes in higher-mass stars. Few Herbig Ae/Be stars have detected magnetic fields, calling into question whether the magnetospheric accretion paradigm developed for low-mass stars can be scaled to higher masses. We present He I 10830 \AA\ line profiles for 64 Herbig Ae/Be stars with a magnetic field measurement in order to test magnetospheric accretion in the physical regime where its efficacy remains uncertain. Of the 5 stars with a magnetic field detection, 1 shows redshifted absorption, indicative of infall, and 2 show blueshifted absorption, tracing mass outflow. The fraction of redshifted and blueshifted absorption profiles in the non-magnetic Herbig Ae/Be stars is remarkably similar, suggesting that the stellar magnetic field does not affect gas kinematics traced by He I 10830 \AA. Line profile morphology does not correlate with the luminosity, rotation rate, mass accretion rate, or disk inclination. Only the detection of a magnetic field and a nearly face-on disk inclination show a correlation (albeit for few sources). This provides further evidence for weaker dipoles and more complex field topologies as stars develop a radiative envelope. The small number of magnetic Herbig Ae/Be stars has already called into question whether magnetospheric accretion can be scaled to higher masses; accretion signatures are not substantially different in magnetic Herbig Ae/Be stars, casting further doubt that they accrete in the same manner as classical T Tauri stars.

The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier. Here, we present ALMA 0.04"-resolution imaging of the pre-transitional disk of V1247 Orionis that reveals an asymmetric ring as well as a sharply-confined crescent structure, resembling morphologies seen in theoretical models of vortex formation. The asymmetric ring (at 0.17"=54 au separation from the star) and the crescent (at 0.38"=120 au) seem smoothly connected through a one-armed spiral arm structure that has been found previously in scattered light. We propose a physical scenario with a planet orbiting at $\sim0.3$"$\approx$100 au, where the one-armed spiral arm detected in polarised light traces the accretion stream feeding the protoplanet. The dynamical influence of the planet clears the gap between the ring and the crescent and triggers two vortices that trap mm-sized particles, namely the crescent and the bright asymmetry seen in the ring. We conducted dedicated hydrodynamics simulations of a disk with an embedded planet, which results in similar spiral-arm morphologies as seen in our scattered light images. At the position of the spiral wake and the crescent we also observe $^{12}$CO (3-2) and H$^{12}$CO$^{+}$ (4-3) excess line emission, likely tracing the increased scale-height in these disk regions.

The first hydrostatic core (FHSC) is the first stable object to form in simulations of star formation. This stage has yet to be observed definitively, although several candidate FHSCs have been reported. We have produced synthetic spectral energy distributions (SEDs) from 3D hydrodynamical simulations of pre-stellar cores undergoing gravitational collapse for a variety of initial conditions. Variations in the initial rotation rate, radius and mass lead to differences in the location of the SED peak and far-infrared flux. Secondly, we attempt to fit the SEDs of five FHSC candidates from the literature and five newly identified FHSC candidates located in the Serpens South molecular cloud with simulated SEDs. The most promising FHSC candidates are fitted by a limited number of model SEDs with consistent properties, which suggests the SED can be useful for placing constraints on the age and rotation rate of the source. The sources we consider most likely to be in FHSC phase are B1-bN, CB17-MMS, Aqu-MM1 and Serpens South candidate K242. We were unable to fit SerpS-MM22, Per-Bolo 58 and Chamaeleon-MMS1 with reasonable parameters, which indicates that they are likely to be more evolved.

We present near-infrared AMBER (R = 12, 000) and CRIRES (R = 100, 000) observations of the Herbig B[e] star MWC297 in the hydrogen Br-gamma-line. Using the VLTI unit telescopes, we obtained a uv-coverage suitable for aperture synthesis imaging. We interpret our velocity-resolved images as well as the derived two-dimensional photocenter displacement vectors, and fit kinematic models to our visibility and phase data in order to constrain the gas velocity field on sub-AU scales. The measured continuum visibilities constrain the orientation of the near-infrared-emitting dust disk, where we determine that the disk major axis is oriented along a position angle of 99.6 +/- 4.8 degrees. The near-infrared continuum emission is 3.6 times more compact than the expected dust-sublimation radius, possibly indicating the presence of highly refractory dust grains or optically thick gas emission in the inner disk. Our velocity-resolved channel maps and moment maps reveal the motion of the Br-gamma-emitting gas in six velocity channels, marking the first time that kinematic effects in the sub-AU inner regions of a protoplanetary disk could be directly imaged. We find a rotation-dominated velocity field, where the blue- and red-shifted emissions are displaced along a position angle of 24 +/- 3 degrees and the approaching part of the disk is offset west of the star. The visibility drop in the line as well as the strong non-zero phase signals can be modeled reasonably well assuming a Keplerian velocity field, although this model is not able to explain the 3 sigma difference that we measure between the position angle of the line photocenters and the position angle of the dust disk. We find that the fit can be improved by adding an outflowing component to the velocity field, as inspired by a magneto-centrifugal disk-wind scenario.

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.

The spin of Wolf-Rayet (WR) stars at low metallicity (Z) is most relevant for our understanding of gravitational wave sources such as GW 150914, as well as the incidence of long-duration gamma-ray bursts (GRBs). Two scenarios have been suggested for both phenomena: one of them involves rapid rotation and quasi-chemical homogeneous evolution (CHE), the other invokes classical evolution through mass loss in single and binary systems. WR spin rates might enable us to test these two scenarios. In order to obtain empirical constraints on black hole progenitor spin, we infer wind asymmetries in all 12 known WR stars in the Small Magellanic Cloud (SMC) at Z = 1/5 Zsun, as well as within a significantly enlarged sample of single and binary WR stars in the Large Magellanic Cloud (LMC at Z = 1/2 Zsun), tripling the sample of Vink (2007). This brings the total LMC sample to 39, making it appropriate for comparison to the Galactic sample. We measure WR wind asymmetries with VLT-FORS linear spectropolarimetry. We report the detection of new line effects in the LMC WN star BAT99-43 and the WC star BAT99-70, as well as the famous WR/LBV HD 5980 in the SMC, which might be evolving chemically homogeneously. With the previous reported line effects in the late-type WNL (Ofpe/WN9) objects BAT99-22 and BAT99-33, this brings the total LMC WR sample to 4, i.e. a frequency of ~10%. Perhaps surprisingly, the incidence of line effects amongst low-Z WR stars is not found to be any higher than amongst the Galactic WR sample, challenging the rotationally-induced CHE model. As WR mass loss is likely Z-dependent, our Magellanic Cloud line-effect WR stars may maintain their surface rotation and fulfill the basic conditions for producing long GRBs, both via the classical post-red supergiant (RSG) or luminous blue variable (LBV) channel, as well as resulting from CHE due to physics specific to very massive stars (VMS).

In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.

High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9\pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/GRAVITY+AMBER near-infrared interferometry allows us to image the circumstellar disks around the individual components with 3 milliarcsecond resolution. We estimate the component masses to $\sim20$ and $\sim18 M_{\odot}$ and find that the radial intensity profiles can be reproduced with an irradiated disk model, where the inner regions are excavated of dust, likely tracing the dust sublimation region in these disks. The circumstellar disks are strongly misaligned with respect to the binary separation vector, which indicates that the tidal forces did not have time to realign the disks, pointing towards a young dynamical age of the system. We constrain the distribution of the Br$\gamma$ and CO-emitting gas using VLTI/GRAVITY spectro-interferometry and VLT/CRIRES spectro-astrometry and find that the secondary is accreting at a higher rate than the primary. VLT/NACO imaging shows $L'$-band emission on 3-4 times larger scales than the binary separation, matching the expected dynamical truncation radius for the circumbinary disk. The IRAS17216-3801 system is $\sim3\times$ more massive and $\sim5\times$ more compact than other high-mass multiplies imaged at infrared wavelengths and the first high-mass protobinary system where circumstellar and circumbinary dust disks could be spatially resolved. This opens exciting new opportunities for studying star-disk interactions and the role of multiplicity in high-mass star formation.

TYC 8241 2652 1 is a young star that showed a strong mid-infrared (mid-IR, 8-25 mu) excess in all observations before 2008 consistent with a dusty disk. Between 2008 and 2010 the mid-IR luminosity of this system dropped dramatically by at least a factor of 30 suggesting a loss of dust mass of an order of magnitude or more. We aim to constrain possible models including removal of disk material by stellar activity processes, the presence of a binary companion, or other explanations suggested in the literature. We present new X-ray observations, optical spectroscopy, near-IR interferometry, and mid-IR photometry of this system to constrain its parameters and further explore the cause of the dust mass loss. In X-rays TYC 8241 2652 1 has all properties expected from a young star: Its luminosity is in the saturation regime and the abundance pattern shows enhancement of O/Fe. The photospheric Ha line is filled with a weak emission feature, indicating chromospheric activity consistent with the observed level of coronal emission. Interferometry does not detect a companion and sets upper limits on the companion mass of 0.2, 0.35, 0.1 and 0.05 M_sun at projected physical separations of 0.1-4 AU,4-5 AU, 5-10 AU, and 10-30 AU, respectively (assuming a distance of 120.9 pc). Our mid-IR measurements, the first of the system since 2012, are consistent with the depleted dust level seen after 2009. The new data confirms that stellar activity is unlikely to destroy the dust in the disk and shows that scenarios where either TYC 8241 2652 1 heats the disk of a binary companion or a potential companion heats the disk of TYC 8241 2652 1 are unlikely.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

We present near and mid-infrared interferometric observations made with the Keck Interferometer Nuller and near-contemporaneous spectro-photometry from the IRTF of 11 well known young stellar objects, several observed for the first time in these spectral and spatial resolution regimes. With AU-level spatial resolution, we first establish characteristic sizes of the infrared emission using a simple geometrical model consisting of a hot inner rim and mid-infrared disk emission. We find a high degree of correlation between the stellar luminosity and the mid-infrared disk sizes after using near-infrared data to remove the contribution from the inner rim. We then use a semi-analytical physical model to also find that the very widely used "star + inner dust rim + flared disk" class of models strongly fails to reproduce the SED and spatially-resolved mid-infrared data simultaneously; specifically a more compact source of mid-infrared emission is required than results from the standard flared disk model. We explore the viability of a modification to the model whereby a second dust rim containing smaller dust grains is added, and find that the two-rim model leads to significantly improved fits in most cases. This complexity is largely missed when carrying out SED modelling alone, although detailed silicate feature fitting by McClure et al. 2013 recently came to a similar conclusion. As has been suggested recently by Menu et al. 2015, the difficulty in predicting mid-infrared sizes from the SED alone might hint at "transition disk"-like gaps in the inner AU; however, the relatively high correlation found in our mid-infrared disk size vs. stellar luminosity relation favors layered disk morphologies and points to missing disk model ingredients instead.

We use a simple, self-consistent, self-gravitating semi-analytic disc model to conduct an examination of the parameter space in which self-gravitating discs may exist. We then use Monte-Carlo radiative transfer to generate synthetic ALMA images of these self-gravitating discs to determine the subset of this parameter space in which they generate non-axisymmetric structure that is potentially detectable by ALMA. Recently, several transition discs have been observed to have non-axisymmetric structure that extends out to large radii. It has been suggested that one possible origin of these asymmetries could be spiral density waves induced by disc self-gravity. We use our simple model to see if these discs exist in the region of parameter space where self-gravity could feasibly explain these spiral features. We find that for self-gravity to play a role in these systems typically requires a disc mass around an order of magnitude higher than the observed disc masses for the systems. The spiral amplitudes produced by self-gravity in the local approximation are relatively weak when compared to amplitudes produced by tidal interactions, or spirals launched at Lindblad resonances due to embedded planets in the disc. As such, we ultimately caution against diagnosing spiral features as being due to self-gravity, unless the disc exists in the very narrow region of parameter space where the spiral wave amplitudes are large enough to produce detectable features, but not so large as to cause the disc to fragment.

Deuterium-to-hydrogen (D/H) enrichments in molecular species provide clues about their original formation environment. The organic materials in primitive solar system bodies have generally higher D/H ratios and show greater D/H variation when compared to D/H in solar system water. We propose this difference arises at least in part due to 1) the availability of additional chemical fractionation pathways for organics beyond that for water, and 2) the higher volatility of key carbon reservoirs compared to oxygen. We test this hypothesis using detailed disk models, including a sophisticated, new disk ionization treatment with a low cosmic ray ionization rate, and find that disk chemistry leads to higher deuterium enrichment in organics compared to water, helped especially by fractionation via the precursors CH$_2$D$^+$/CH$_3^+$. We also find that the D/H ratio in individual species varies significantly depending on their particular formation pathways. For example, from $\sim20-40$ AU, CH$_4$ can reach $\rm{D/H\sim2\times10^{-3}}$, while D/H in CH$_3$OH remains locally unaltered. Finally, while the global organic D/H in our models can reproduce intermediately elevated D/H in the bulk hydrocarbon reservoir, our models are unable to reproduce the most deuterium-enriched organic materials in the solar system, and thus our model requires some inheritance from the cold interstellar medium from which the Sun formed.

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

We assess the multi-wavelength observable properties of the bow shock around a runaway early type star using a combination of hydrodynamical modelling, radiative transfer calculations and synthetic imaging. Instabilities associated with the forward shock produce dense knots of material which are warm, ionised and contain dust. These knots of material are responsible for the majority of emission at far infra-red, H alpha and radio wavelengths. The large scale bow shock morphology is very similar and differences are primarily due to variations in the assumed spatial resolution. However infra-red intensity slices (at 22 microns and 12 microns) show that the effects of a temperature gradient can be resolved at a realistic spatial resolution for an object at a distance of 1 kpc.

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.

Radiation hydrodynamics (RHD) simulations are used to study many astrophysical phenomena, however they require the use of simplified radiation transport and thermal prescriptions to reduce computational cost. In this paper we present a systematic study of the importance of microphysical processes in RHD simulations using the example of D-type HII region expansion. We compare the simplest hydrogen-only models with those that include: ionisation of H, He, C, N, O, S and Ne, different gas metallicity, non-LTE metal line blanketed stellar spectral models of varying metallicity, radiation pressure, dust and treatment of photodissociation regions. Each of these processes are explicitly treated using modern numerical methods rather than parameterisation. In line with expectations, changes due to microphysics in either the effective number of ionising photons or the thermal structure of the gas lead to differences in D-type expansion. In general we find that more realistic calculations lead to the onset of D-type expansion at smaller radii and a slower subsequent expansion. Simulations of star forming regions using simplified microphysics are therefore likely overestimating the strength of radiative feedback. We find that both variations in gas metallicity and the inclusion of dust can affect the ionisation front evolution at the 10-20 per cent level over 500kyr, which could substantially modify the results of simplified 3D models including feedback. Stellar metallicity, radiation pressure and the inclusion of photodissociation regions are all less significant effects at the 1 per cent level or less, rendering them of minor importance in the modelling the dynamical evolution of HII regions.

We examine changes in the molecular abundances resulting from increased heating due to a self-luminous planetary companion embedded within a narrow circumstellar disk gap. Using 3D models that include stellar and planetary irradiation, we find that luminous young planets locally heat up the parent circumstellar disk by many tens of Kelvin, resulting in efficient thermal desorption of molecular species that are otherwise locally frozen out. Furthermore, the heating is deposited over large regions of the disk, $\pm5$ AU radially and spanning $\lesssim60^\circ$ azimuthally. From the 3D chemical models, we compute rotational line emission models and full ALMA simulations, and find that the chemical signatures of the young planet are detectable as chemical asymmetries in $\sim10h$ observations. HCN and its isotopologues are particularly clear tracers of planetary heating for the models considered here, and emission from multiple transitions of the same species is detectable, which encodes temperature information in addition to possible velocity information from the spectra itself. We find submillimeter molecular emission will be a useful tool to study gas giant planet formation in situ, especially beyond $R\gtrsim10$ AU.

We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo estimator and enables a complete treatment of the detailed microphysics, including polychromatic radiation and anisotropic scattering, in both the free-streaming and optically-thick limits. Since the new method is computationally demanding we have developed two new methods that speed up the algorithm. The first is a photon packet splitting algorithm that enables efficient treatment of the Monte Carlo process in very optically thick regions. The second is a parallelisation method that distributes the Monte Carlo workload over many instances of the hydrodynamic domain, resulting in excellent scaling of the radiation step. We also describe the implementation of a sink particle method that enables us to follow the accretion onto, and the growth of, the protostars. We detail the results of extensive testing and benchmarking of the new algorithms.

We present photometric and spectroscopic models of the Classical T Tauri star AA Tau. Photometric and spectroscopic variability present in observations of AA Tau is attributed to a magnetically induced warp in the accretion disc, periodically occulting the photosphere on an 8.2--day timescale. Emission line profiles show signatures of both infall, attributed to magnetospherically accreting material, and outflow. Using the radiative transfer code TORUS, we have investigated the geometry and kinematics of AA Tau's circumstellar disc and outflow, which is modelled here as a disc wind. Photometric models have been used to constrain the aspect ratio of the disc, the offset angle of the magnetosphere dipole with respect to the stellar rotation axis, and the inner radius of the circumstellar disc. Spectroscopic models have been used to constrain the wind and magnetosphere temperatures, wind acceleration parameter, and mass loss rate. We find observations are best fitted by models with a mass accretion rate of $5\times10^{-9}$ M$_\odot$ yr$^{-1}$, a dipole offset of between $10^\circ$ and $20^\circ$, a magnetosphere that truncates the disc from 5.2 to 8.8 R$_\star$, a mass-loss-rate to accretion-rate ratio of ~ 0.1, a magnetosphere temperature of 8500 -- 9000 K, and a disc wind temperature of 8000 K.

Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk. Utilizing a comprehensive treatment of disk ionization, we find that ion-driven deuterium pathways are inefficient, curtailing the disk's deuterated water formation and its viability as the sole source for the solar system's water. This finding implies that if the solar system's formation was typical, abundant interstellar ices are available to all nascent planetary systems.

Among the most fascinating and hotly-debated areas in contemporary astrophysics are the means by which planetary systems are assembled from the large rotating disks of gas and dust which attend a stellar birth. Although important work has already been, and is still being done both in theory and observation, a full understanding of the physics of planet formation can only be achieved by opening observational windows able to directly witness the process in action. The key requirement is then to probe planet-forming systems at the natural spatial scales over which material is being assembled. By definition, this is the so-called Hill Sphere which delineates the region of influence of a gravitating body within its surrounding environment. The Planet Formation Imager project (PFI) has crystallized around this challenging goal: to deliver resolved images of Hill-Sphere-sized structures within candidate planet-hosting disks in the nearest star-forming regions. In this contribution we outline the primary science case of PFI. For this purpose, we briefly review our knowledge about the planet-formation process and discuss recent observational results that have been obtained on the class of transition disks. Spectro-photometric and multi-wavelength interferometric studies of these systems revealed the presence of extended gaps and complex density inhomogeneities that might be triggered by orbiting planets. We present detailed 3-D radiation-hydrodynamic simulations of disks with single and multiple embedded planets, from which we compute synthetic images at near-infrared, mid-infrared, far-infrared, and sub-millimeter wavelengths, enabling a direct comparison of the signatures that are detectable with PFI and complementary facilities such as ALMA. From these simulations, we derive some preliminary specifications that will guide the array design and technology roadmap of the facility.

We study the evaporation and condensation of CO and CO_2 during the embedded stages of low-mass star formation by using numerical simulations. We focus on the effect of luminosity bursts, similar in magnitude to FUors and EXors, on the gas-phase abundance of CO and CO_2 in the protostellar disk and infalling envelope. The evolution of a young protostar and its environment is followed based on hydrodynamical models using the thin-disk approximation, coupled with a stellar evolution code and phase transformations of CO and CO_2. The accretion and associated luminosity bursts in our model are caused by disk gravitational fragmentation followed by quick migration of the fragments onto the forming protostar. We found that bursts with luminosity on the order of 100-200 L_sun can evaporate CO ices in part of the envelope. The typical freeze-out time of the gas-phase CO onto dust grains in the envelope (a few kyr) is much longer than the burst duration (100-200 yr). This results in an increased abundance of the gas-phase CO in the envelope long after the system has returned into a quiescent stage. In contrast, luminosity bursts can evaporate CO_2 ices only in the disk, where the freeze-out time of the gas-phase CO_2 is comparable to the burst duration. We thus confirm that luminosity bursts can leave long-lasting traces in the abundance of gas-phase CO in the infalling envelope, enabling the detection of recent bursts as suggested by previous semi-analytical studies.

Pre-transitional disks are protoplanetary disks with a gapped disk structure, potentially indicating the presence of young planets in these systems. In order to explore the structure of these objects and their gap-opening mechanism, we observed the pre-transitional disk V1247 Orionis using the Very Large Telescope Interferometer, the Keck Interferometer, Keck-II, Gemini South, and IRTF. This allows us spatially resolve the AU-scale disk structure from near- to mid-infrared wavelengths (1.5 to 13 {\mu}m), tracing material at different temperatures and over a wide range of stellocentric radii. Our observations reveal a narrow, optically-thick inner-disk component (located at 0.18 AU from the star) that is separated from the optically thick outer disk (radii >46 AU), providing unambiguous evidence for the existence of a gap in this pre-transitional disk. Surprisingly, we find that the gap region is filled with significant amounts of optically thin material with a carbon-dominated dust mineralogy. The presence of this optically thin gap material cannot be deduced solely from the spectral energy distribution, yet it is the dominant contributor at mid-infrared wavelengths. Furthermore, using Keck/NIRC2 aperture masking observations in the H, K', and L' band, we detect asymmetries in the brightness distribution on scales of about 15-40 AU, i.e. within the gap region. The detected asymmetries are highly significant, yet their amplitude and direction changes with wavelength, which is not consistent with a companion interpretation but indicates an inhomogeneous distribution of the gap material. We interpret this as strong evidence for the presence of complex density structures, possibly reflecting the dynamical interaction of the disk material with sub-stellar mass bodies that are responsible for the gap clearing.

We investigate observational signatures of triggered star formation in bright rimmed clouds (BRCs) by using molecular line transfer calculations based on radiation-hydrodynamic radiatively-driven-implosion models. We find that for BRCs the separation in velocity between the line profile peak of an optically thick and an optically thin line is determined by both the observer viewing angle and the density of the shell driving into the cloud. In agreement with observations, we find that most BRC line profiles are symmetric and that asymmetries can be either red or blue, in contrast to the blue-dominance expected for a collapsing cloud. Asymmetries in the line profiles arise when an optically thick line is dominated by the shell and an optically thin line is dominated by the cloud interior to the shell. The asymmetries are red or blue depending on whether the shell is moving towards or away from the observer respectively. Using the known motions of the molecular gas in our models we rule out the envelope expansion with core collapse mechanism as the cause of the lack of blue-asymmetry in our simulated observations. We show that the absence of a strong photon dominated region (PDR) around a BRC may not rule out the presence of triggered star formation: if the BRC line profile has a strong blue component then the shell is expected to be driving towards the observer, suggesting that the cloud is being viewed from behind and the PDR is obstructed. This could explain why BRCs such as SFO 80, 81 and 86 have a blue secondary peak and only a weak PDR inferred at 8 microns. Finally we also test the use of 12CO, 13CO and C18O as diagnostics of cloud mass, temperature and column density. We find that the inferred conditions are in reasonable agreement with those from the models.

(shortened) We investigate the properties of Galactic WR stars and their environment to identify evolutionary channels that may lead to the formation of LGRBs. To this purpose we compile available information on the spectropolarimetric properties, the presence of CS ejecta, and the CS velocities in the environment of Galactic WR stars. We use linear line-depolarization as an indicator of rotation, nebular morphology as an indicator of stellar ejecta, and velocity patterns in UV absorption features as an indicator of increased velocities in the CS environment. We find that the ~23% WR stars with "possible ejecta nebulae" dominate the population of WR stars with spectropolarimetric signatures of rotation, while WR stars without such nebulae only rarely show indications of rotation. The corresponding objects are most likely in an early stage after a preceding RSG or LBV phase, and have not yet lost their angular momenta due to the strong mass-loss in the WR phase. From their photometric periods we estimate rotation parameters in the range omega = 0.04...0.25, corresponding to moderate rotation speeds of 36...120 km/s. These values are very uncertain, but comply with the specific surface angular momentum requirement for LGRB progenitors. Our results indicate that, in the Galaxy, mainly "young" WR stars shortly after a RSG/LBV phase show spectropolarimetric signatures of rotation. Their rotation rates are thus likely enhanced with respect to the majority of Galactic WR stars. According to their estimated specific surface angular momenta, a subgroup of stars exploding in this phase may represent an evolutionary channel towards LGRBs at high metallicities, comparable to the Galaxy.

We present a radiation hydrodynamics simulation of the formation of a massive star using a Monte Carlo treatment for the radiation field. We find that strong, high speed bipolar cavities are driven by the radiation from the protostar, and that accretion occurs stochastically from a circumstellar disc. We have computed spectral energy distributions and images at each timestep, which may in future be used to compare our models with photometric, spectroscopic, and interferometric observations of young massive stellar objects.

One of the main reasons that triggered star formation is contentious is the failure to accurately link the observations with models in a detailed, quantitative, way. It is therefore critical to continuously test and improve the model details and methods with which comparisons to observations are made. We use a Monte Carlo radiation transport and hydrodynamics code TORUS to show that the diffuse radiation field has a significant impact on the outcome of radiatively driven implosion (RDI) models. We also calculate SEDs and synthetic images from the models to test observational diagnostics that are used to determine bright rimmed cloud conditions and search for signs of RDI.

We produce synthetic images and SEDs from radiation hydrodynamical simulations of radiatively driven implosion. The synthetically imaged bright rimmed clouds (BRCs) are morphologically similar to those observed in star forming regions. Using nebular diagnostic line-ratios, simulated Very Large Array (VLA) radio images, H{\alpha} imaging and SED fitting we compute the neutral cloud and ionized boundary layer gas densities and temperatures and perform a virial stability analysis for each model cloud. We determine that the neutral cloud temperatures derived by SED fitting are hotter than the dominant neutral cloud temperature by 1 - 2 K due to emission from warm dust. This translates into a change in the calculated cloud mass by 8-35 %. Using a constant mass conversion factor (C{\nu}) for BRCs of different class is found to give rise to errors in the cloud mass of up to a factor of 3.6. The ionized boundary layer (IBL) electron temperature calculated using diagnostic line ratios is more accurate than assuming the canonical value adopted for radio diagnostics of 10^4 K. Both radio diagnostics and diagnostic line ratios are found to underestimate the electron density in the IBL. Each system is qualitatively correctly found to be in a state in which the pressure in the ionized boundary layer is greater than the supporting cloud pressure, implying that the objects are being compressed. We find that observationally derived mass loss estimates agree with those on the simulation grid and introduce the concept of using the mass loss flux to give an indication of the relative strength of photo-evaporative flow between clouds. The effect of beam size on these diagnostics in radio observations is found to be a mixing of the bright rim and ambient cloud and HII region fluxes, which leads to an underestimate of the cloud properties relative to a control diagnostic.

We present derived stellar and disc parameters for a sample of Taurus brown dwarfs both with and without evidence of an associated disc. These parameters have been derived using an online fitting tool (http://bd-server.astro.ex.ac.uk/), which includes a statistically robust derivation of uncertainties, an indication of pa- rameter degeneracies, and a complete treatment of the input photometric and spectroscopic observations. The observations of the Taurus members with indications of disc presence have been fitted using a grid of theoretical models including detailed treatments of physical processes accepted for higher mass stars, such as dust sublimation, and a simple treatment of the accretion flux. This grid of models has been designed to test the validity of the adopted physical mechanisms, but we have also constructed models using parameterisation, for example semi-empirical dust sublimation radii, for users solely interested in parameter derivation and the quality of the fit. The parameters derived for the naked and disc brown dwarf systems are largely consistent with literature observations. However, our inner disc edge locations are consistently closer to the star than previous results and we also derive elevated accretion rates over non-SED based accretion rate derivations. For inner edge locations we attribute these differences to the detailed modelling we have performed of the disc structure, particularly at the crucial inner edge where departures in geometry from the often adopted vertical wall due to dust sublimation (and therefore accretion flux) can compensate for temperature (and therefore distance) changes to the inner edge of the dust disc. In the case of the elevated derived accretion rates, in some cases, this may be caused by the intrinsic stellar luminosities of the targets exceeding that predicted by the isochrones we have adopted.

Long gamma-ray bursts involve the most powerful cosmic explosions since the Big Bang. Whilst it has been established that GRBs are related to the death throes of massive stars, the identification of their progenitors has proved challenging. Theory suggests that rotating Wolf-Rayet stars are the best candidates, but their strong stellar winds shroud their surfaces, preventing a direct measurement of their rotation. Fortunately, linear spectropolarimetry may be used to probe the flattening of their winds due to stellar spin. Spectropolarimetry surveys show that an 80% majority of WR stars have spherically symmetric winds and are thus rotating slowly, yet a small 20% minority display a spectropolarimetric signature indicative of rotation. Here we find a highly significant correlation between WR objects that carry the signature of stellar rotation and the subset of WR stars with ejecta nebulae that have only recently transitioned from a red sugergiant or luminous blue variable phase. As these youthful WR stars have yet to spin-down due to mass loss, they are the best candidate GRB progenitors identified to date. When we take recently published WR ejecta nebula numbers we find that five out of the six line-effect WR stars are surrounded by ejecta nebulae. The statistics imply that the null hypothesis of no correlation between line-effect WR stars and ejecta nebulae can be rejected at the 0.0004% level. Given that four line-effect and WR ejecta nebula have spectroscopically been confirmed to contain nucleosynthetic products, we argue that the correlation is both statistically significant and physically convincing. The implication is that we have identified a WR sub-population that fulfills the necessary criteria for making GRBs. Finally, we discuss the potential of identifying GRB progenitors via spectropolarimetry with extremely large telescopes.

We investigate the effect of including diffuse field radiation when modelling the radiatively driven implosion of a Bonnor-Ebert sphere (BES). Radiation-hydrodynamical calculations are performed by using operator splitting to combine Monte Carlo photoionization with grid-based Eulerian hydrodynamics that includes self-gravity. It is found that the diffuse field has a significant effect on the nature of radiatively driven collapse which is strongly coupled to the strength of the driving shock that is established before impacting the BES. This can result in either slower or more rapid star formation than expected using the on-the-spot approximation depending on the distance of the BES from the source object. As well as directly compressing the BES, stronger shocks increase the thickness and density in the shell of accumulated material, which leads to short, strong, photo-evaporative ejections that reinforce the compression whenever it slows. This happens particularly effectively when the diffuse field is included as rocket motion is induced over a larger area of the shell surface. The formation and evolution of 'elephant trunks' via instability is also found to vary significantly when the diffuse field is included. Since the perturbations that seed instabilities are smeared out elephant trunks form less readily and, once formed, are exposed to enhanced thermal compression.

We present multidimensional non-LTE radiative transfer models of hydrogen and helium line profiles formed in the accretion flows and the outflows near the star-disk interaction regions of classical T Tauri stars (CTTSs). The statistical equilibrium calculations, performed under the assumption of the Sobolev approximation using the radiative transfer code TORUS, has been improved to include He I and He II energy levels. This allows us to probe the physical conditions of the inner wind of CTTSs by simultaneously modelling the robust wind diagnostic line He I (10830) and the accretion diagnostic lines such as Pa-beta, Br-gamma and He I (5876). The code has been tested in 1 and 2-D problems, and we have shown that the results are in agreement with established codes. We apply the model to the complex flow geometries of CTTSs. Example model profiles are computed using the combinations of (1) magnetospheric accretion and disc wind, and (2) magnetospheric accretion and the stellar wind. In both cases, the model produces line profiles which are qualitatively similar to those found in observations. Our models are consistent with the scenario in which the narrow blueshifted absorption component of He I (10830) seen in observations is caused by a disc wind, and the wider blueshifted absorption component (the P-Cygni profile) is caused by a bipolar stellar wind. However, we do not have a strong constraint on the relative importance of the wind and the magnetosphere for the `emission' component. Our preliminary calculations suggest that the temperatures of the disc wind, stellar wind and the magnetosphere cannot be much higher than ~10,000 K, on the basis of the strengths of hydrogen lines. With these low temperatures, we find that the photoionzation by high energy photons (e.g. X-rays) is necessary to produce He I (10830) in emission and to produce the blueshifted absorption components.

A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick and free-streaming regimes simultaneously. The algorithm is tested on a variety of 1-d and 2-d problems, and good agreement with benchmark solutions is found. The method is used to calculate the time-varying spectral energy distribution from a circumstellar disc illuminated by a protostar whose accretion luminosity is varying. It is shown that the time lag between the optical variability and the infrared variability results from a combination of the photon travel time and the thermal response in the disc, and that the lag is an approximately linear function of wavelength.

We present a grid of models of accreting brown dwarf systems with circumstellar discs. The calculations involve a self-consistent solution of both vertical hydrostatic and radiative equilibrium along with a sophisticated treatment of dust sublimation. We have simulated observations of the spectral energy distributions and several broadband photometric systems. Analysis of the disc structures and simulated observations reveal a natural dichotomy in accretion rates, with \logmdot $>-$9 and $\leq -$9 classed as extreme and typical accretors respectively. Derivation of ages and masses from our simulated photometry using isochrones is demonstrated to be unreliable even for typical accretors. Although current brown dwarf disc candidate selection criteria have been shown to be largely reliable when applied to our model grid we suggest improved selection criteria in several colour indices. We show that as accretion rates increase brown dwarf disc systems are less likely to be correctly identified. This suggests that, within our grid, systems with higher accretion rates would be preferentially lost during brown dwarf target selection. We suggest that observations used to assert a $\dot{M}\propto M_*^2$ relationship may contain an intrinsic selection bias.

We present the first non-LTE, co-moving frame molecular line calculations of a star-forming cluster simulated using smoothed particle hydrodynamics (SPH), from which we derive high-resolution synthetic observations. We have resampled a particle representation onto an adaptive mesh and self-consistently solved the equations of statistical equilibrium in the co-moving frame, using TORUS, a three-dimensional adaptive mesh refined (AMR) radiative transfer (RT) code. We verified the applicability of the code to the conditions of the SPH simulation by testing its output against other codes. We find that the level populations obtained for optically thick and thin scenarios closely match the ensemble average of the other codes. We have used the code to obtain non-LTE level populations of multiple molecular species throughout the cluster and have created three-dimensional velocity-resolved spatial maps of the emergent intensity. Line profiles of cores traced by N2H+ (1-0) are compared to probes of low density gas, 13CO (1-0) and C18O (1-0), surrounding the cores along the line of sight. The relative differences of the line-centre velocities are shown to be small compared to the velocity dispersion, matching recent observations. We conclude that one cannot reject competitive accretion as a viable theory of star formation based on observed velocity profiles.

We present results from combining a grid-based radiative transfer code with a Smoothed Particle Hydrodynamics code to produce a flexible system for modelling radiation hydrodynamics. We use a benchmark model of a circumstellar disc to determine a robust method for constructing a gridded density distribution from SPH particles. The benchmark disc is then used to determine the accuracy of the radiative transfer results. We find that the SED and the temperature distribution within the disc are sensitive to the representation of the disc inner edge, which depends critically on both the grid and SPH resolution. The code is then used to model a circumstellar disc around a T-Tauri star. As the disc adjusts towards equilibrium vertical motions in the disc are induced resulting in scale height enhancements which intercept radiation from the central star. Vertical transport of radiation enables these perturbations to influence the mid-plane temperature of the disc. The vertical motions decay over time and the disc ultimately reaches a state of simultaneous hydrostatic and radiative equilibrium.

(abridged) As the favoured progenitors of long-duration gamma-ray bursts, massive stars may represent our best signposts of individual objects in the early Universe, but special conditions seem required to make these bursters, which might originate from the progenitor's rapid rotation and associated asymmetry. To obtain empirical constraints on the interplay between stellar rotation and wind asymmetry, we perform linear Halpha spectropolarimetry on a sample of 18 spectroscopically peculiar massive O stars, including OVz, Of?p, Oe, and Onfp stars, supplemented by an earlier sample of 20 O supergiants. Despite their rapid rotation (with vsin(i) up to 400 km/s) most O-type stars are found to be spherically symmetric, but with notable exceptions amongst specific object classes. We divide the peculiar O stars into four distinct categories: Groups III and IV include the Oe stars and Onfp stars, which are on the high-end tail of the O star rotation distribution and have in the past been claimed to be embedded in disks. Here we report the detection of a classical depolarization ``line effect'' in the Oe star HD 45314, but the overall incidence of line effects amongst Oe stars is significantly lower (1 out of 6) than amongst Be stars. The chance that the Oe and Be datasets are drawn from the same parent population is negligible (with 95% confidence). This implies there is as yet no evidence for a disk hypothesis in Oe stars, providing relevant constraints on the physical mechanism that is responsible for the Be phenomenon. Finally, we find that 3 out of 4 of the group IV Onfp stars show evidence for complex polarization effects, which are likely related to rapid rotation, and we speculate on the evolutionary links to B[e] stars.

We have produced brightness and magnetic field maps of the surfaces of CV Cha and CR Cha: two actively accreting G and K-type T Tauri stars in the Chamaeleon I star-forming cloud with ages of 3-5 Myr. Our magnetic field maps show evidence for strong, complex multi-polar fields similar to those obtained for young rapidly rotating main sequence stars. Brightness maps indicate the presence of dark polar caps and low latitude spots -- these brightness maps are very similar to those obtained for other pre-main sequence and rapidly rotating main sequence stars. Only two other classical T Tauri stars have been studied using similar techniques so far: V2129 Oph and BP Tau. CV Cha and CR Cha show magnetic field patterns that are significantly more complex than those recovered for BP Tau, a fully convective T Tauri star. We discuss possible reasons for this difference and suggest that the complexity of the stellar magnetic field is related to the convection zone; with more complex fields being found in T Tauri stars with radiative cores (V2129 Oph, CV Cha and CR Cha). However, it is clearly necessary to conduct magnetic field studies of T Tauri star systems, exploring a wide range of stellar parameters in order to establish how they affect magnetic field generation, and thus how these magnetic fields are likely to affect the evolution of T Tauri star systems as they approach the main sequence.

In response to ESO's call for proposals for second generation instruments for the Very Large Telescope Interferometer (VLTI), a consortium is currently developing the VLTI-Spectro-Imager (VSI). In the context of the Phase A study, a science group has prepared a science case taking advantage of the expected performances of VSI. Among several science topics, the case of dust making Wolf-Rayet binaries producing the so-called pinwheel nebulae has been considered. Here, we review the main specifications of VSI, and we provide preliminary results expected to illustrate the imaging capabilities of VSI, and the interest for the study of pinwheel nebulae similar to those formed close to well-known systems such as WR98a and WR104.

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.

Classically, optical and near-infrared interferometry have relied on closure phase techniques to produce images. Such techniques allow us to achieve modest dynamic ranges. In order to test the feasibility of next generation optical interferometers in the context of the VLTI-spectro-imager (VSI), we have embarked on a study of image reconstruction and analysis. Our main aim was to test the influence of the number of telescopes, observing nights and distribution of the visibility points on the quality of the reconstructed images. Our results show that observations using six Auxiliary Telescopes (ATs) during one complete night yield the best results in general and is critical in most science cases; the number of telescopes is the determining factor in the image reconstruction outcome. In terms of imaging capabilities, an optical, six telescope VLTI-type configuration and ~200 meter baseline will achieve 4 mas spatial resolution, which is comparable to ALMA and almost 50 times better than JWST will achieve at 2.2 microns. Our results show that such an instrument will be capable of imaging, with unprecedented detail, a plethora of sources, ranging from complex stellar surfaces to microlensing events.

One of the aims of next generation optical interferometric instrumentation is to be able to make use of information contained in the visibility phase to construct high dynamic range images. Radio and optical interferometry are at the two extremes of phase corruption by the atmosphere. While in radio it is possible to obtain calibrated phases for the science objects, in the optical this is currently not possible. Instead, optical interferometry has relied on closure phase techniques to produce images. Such techniques allow only to achieve modest dynamic ranges. However, with high contrast objects, for faint targets or when structure detail is needed, phase referencing techniques as used in radio interferometry, should theoretically achieve higher dynamic ranges for the same number of telescopes. Our approach is not to provide evidence either for or against the hypothesis that phase referenced imaging gives better dynamic range than closure phase imaging. Instead we wish to explore the potential of this technique for future optical interferometry and also because image reconstruction in the optical using phase referencing techniques has only been performed with limited success. We have generated simulated, noisy, complex visibility data, analogous to the signal produced in radio interferometers, using the VLTI as a template. We proceeded with image reconstruction using the radio image reconstruction algorithms contained in AIPS IMAGR (CLEAN algorithm). Our results show that image reconstruction is successful in most of our science cases, yielding images with a 4 milliarcsecond resolution in K band. (abridged)

We present comprehensive models for the Herbig Ae stars MWC275 and AB Aur that aim to explain their spectral energy distribution (from UV to millimeter) and long baseline interferometry (from near-infrared to millimeter) simultaneously. Data from the literature, combined with new mid-infrared (MIR) interferometry from the Keck Segment Tilting Experiment, are modeled using an axisymmetric Monte Carlo radiative transfer code. Models in which most of the near-infrared (NIR) emission arises from a dust rim fail to fit the NIR spectral energy distribution (SED) and sub-milli-arcsecond NIR CHARA interferometry. Following recent work, we include an additional gas emission component with similar size scale to the dust rim, inside the sublimation radius, to fit the NIR SED and long-baseline NIR interferometry on MWC275 and AB Aur. In the absence of shielding of star light by gas, we show that the gas-dust transition region in these YSOs will have to contain highly refractory dust, sublimating at ~1850K. Despite having nearly identical structure in the thermal NIR, the outer disks of MWC275 and AB Aur differ substantially. In contrast to the AB Aur disk, MWC275 lacks small grains in the disk atmosphere capable of producing significant 10-20micron emission beyond ~7AU, forcing the outer regions into the "shadow" of the inner disk

Using the longest optical-interferometeric baselines currently available, we have detected strong near-infrared (NIR) emission from inside the dust-destruction radius of Herbig Ae stars MWC275 and AB Aur. Our sub-milli-arcsecond resolution observations unambiguously place the emission between the dust-destruction radius and the magnetospheric co-rotation radius. We argue that this new component corresponds to hot gas inside the dust-sublimation radius, confirming recent claims based on spectrally-resolved interferometry and dust evaporation front modeling.

From observations collected with the ESPaDOnS and NARVAL spectropolarimeters, we report the detection of Zeeman signatures on the classical T Tauri star BP Tau. Circular polarisation signatures in photospheric lines and in narrow emission lines tracing magnetospheric accretion are monitored throughout most of the rotation cycle of BP Tau at two different epochs in 2006. We observe that rotational modulation dominates the temporal variations of both unpolarised and circularly polarised spectral proxies tracing the photosphere and the footpoints of accretion funnels. From the complete data sets at each epoch, we reconstruct the large-scale magnetic topology and the location of accretion spots at the surface of BP Tau using tomographic imaging. We find that the field of BP Tau involves a 1.2 kG dipole and 1.6 kG octupole, both slightly tilted with respect to the rotation axis. Accretion spots coincide with the two main magnetic poles at high latitudes and overlap with dark photospheric spots; they cover about 2% of the stellar surface. The strong mainly-axisymmetric poloidal field of BP Tau is very reminiscent of magnetic topologies of fully-convective dwarfs. It suggests that magnetic fields of fully-convective cTTSs such as BP Tau are likely not fossil remants, but rather result from vigorous dynamo action operating within the bulk of their convective zones. Preliminary modelling suggests that the magnetosphere of BP Tau extends to distances of at least 4 R* to ensure that accretion spots are located at high latitudes, and is not blown open close to the surface by a putative stellar wind. It apparently succeeds in coupling to the accretion disc as far out as the corotation radius, and could possibly explain the slow rotation of BP Tau.

We present three-dimensional (3-D) simulations of rotationally induced line variability arising from complex circumstellar environment of classical T Tauri stars (CTTS) using the results of the 3-D magnetohydrodynamic (MHD) simulations of Romanova et al., who considered accretion onto a CTTS with a misaligned dipole magnetic axis with respect to the rotational axis. The density, velocity and temperature structures of the MHD simulations are mapped on to the radiative transfer grid, and corresponding line source function and the observed profiles of neutral hydrogen lines (H-beta, Pa-beta and Br-gamma) are computed using the Sobolev escape probability method. We study the dependency of line variability on inclination angles (i) and magnetic axis misalignment angles (Theta). By comparing our models with the Pa-beta profiles of 42 CTTS observed by Folha & Emerson, we find that models with a smaller misaligngment angle (Theta<~15 deg.) are more consistent with the observations which show that majority of Pa-beta are rather symmetric around the line centre. For a high inclination system with a small dipole misalignment angle (Theta ~ 15 deg.), only one accretion funnel (on the upper hemisphere) is visible to an observer at any given rotational phase. This can cause an anti-correlation of the line equivalent width in the blue wing (v<0) and that in the red wing (v>0) over a half of a rotational period, and a positive correlation over other half. We find a good overall agreement of the line variability behaviour predicted by our model and those from observations. (Abridged)

From observations collected with the ESPaDOnS spectropolarimeter, we report the discovery of magnetic fields at the surface of the mildly accreting classical T Tauri star V2129 Oph. Zeeman signatures are detected, both in photospheric lines and in the emission lines formed at the base of the accretion funnels linking the disc to the protostar, and monitored over the whole rotation cycle of V2129 Oph. We observe that rotational modulation dominates the temporal variations of both unpolarized and circularly polarized line profiles. We reconstruct the large-scale magnetic topology at the surface of V2129 Oph from both sets of Zeeman signatures simultaneously. We find it to be rather complex, with a dominant octupolar component and a weak dipole of strengths 1.2 and 0.35 kG, respectively, both slightly tilted with respect to the rotation axis. The large-scale field is anchored in a pair of 2-kG unipolar radial field spots located at high latitudes and coinciding with cool dark polar spots at photospheric level. This large-scale field geometry is unusually complex compared to those of non-accreting cool active subgiants with moderate rotation rates. As an illustration, we provide a first attempt at modelling the magnetospheric topology and accretion funnels of V2129 Oph using field extrapolation. We find that the magnetosphere of V2129 Oph must extend to about 7R* to ensure that the footpoints of accretion funnels coincide with the high-latitude accretion spots on the stellar surface. It suggests that the stellar magnetic field succeeds in coupling to the accretion disc as far out as the corotation radius, and could possibly explain the slow rotation of V2129 Oph. The magnetospheric geometry we derive produces X-ray coronal fluxes typical of those observed in cTTSs.

Stellar magnetic fields can be investigated using several, very complementary approaches. While conventional spectroscopy is capable of estimating the average magnetic strength of potentially complex field configurations thanks to its low sensitivity to the vector properties of the field, spectropolarimetry can be used to map the medium- and large-scale structure of magnetic topologies. In particular, the latter approach allows one to retrieve information about the poloidal and toroidal components of the large-scale dynamo fields in low-mass stars, and thus to investigate the physical processes that produce them. Similarly, this technique can be used to explore how magnetic fields couple young stars to their massive accretion disc and thus to estimate how much mass and angular momentum are transfered to the newly-born low-mass star. We present here the latest results in this field obtained with spectropolarimetry, with special emphasis on the surprising discoveries obtained on very-low mass fully-convective stars and classical T Tauri stars thanks to the ESPaDOnS spectropolarimeter recently installed on the 3.6m Canada-France-Hawaii Telescope.

Dust-grain growth and settling are the first steps towards planet formation. An understanding of dust physics is therefore integral to a complete theory of the planet formation process. In this paper, we explore the possibility of using the dust evaporation front in YSO disks (`the inner rim') as a probe of the dust physics operating in circumstellar disks. The geometry of the rim depends sensitively on the composition and spatial distribution of dust. Using radiative transfer and hydrostatic equilibrium calculations we demonstrate that dust growth and settling can curve the evaporation front dramatically (from a cylindrical radius of about 0.5 AU in the disk mid-plane to 1.2 AU in the disk upper layers for an A0 star). We compute synthetic images and interferometric visibilities for our representative rim models and show that the current generation of near-IR long-baseline interferometers (VLTI, CHARA) can strongly constrain the dust properties of circumstellar disks, shedding light on the relatively poorly understood processes of grain growth, settling and turbulent mixing.

Working with the submillimetre continuum map of the Perseus molecular cloud (Hatchell et al. 2005), we aimed to determine the evolutionary stage of each submm core in Perseus, and investigate the lifetimes of these phases. We compile spectral energy distributions (SEDs) from 2MASS, Spitzer IRAC, Michelle, IRAS, SCUBA and Bolocam data. Sources are classified starless/protostellar on the basis of infrared and/or outflow detections and Class I/Class 0 on the basis of Tbol, Lbol/Lsmm and F_{3.6}/F_{850}. In order to investigate the dependence of these evolutionary indicators on mass, we construct radiative transfer models of Class 0 sources. Of the submm cores, 56/103 (54%) are confirmed protostars on the basis of infrared emission or molecular outflows. Of these, 22 are classified Class 1 on the basis of three evolutionary indicators, 34 are Class 0, and the remaining 47 are assumed starless. Perseus contains a much greater fraction of Class 0 sources than either Taurus or Rho Oph. Comparing the protostellar with the T Tauri population, the lifetime of the protostellar phase in Perseus is 0.25-0.67 Myr (95% confidence limits). The relative lifetime of the Class 0 and Class 1 phases are similar. We find that for the same source geometry but different masses, evolutionary indicators such as Tbol vary their value. It is therefore not always appropriate to use a fixed threshold to separate Class 0 and Class I sources. More modelling is required to determine the observational characteristics of the Class 0/Class I boundary over a range of masses.

We report the results of a high angular resolution near-infrared survey of dusty Wolf-Rayet stars using the Keck-1 Telescope, including new multi-wavelength images of the pinwheel nebulae WR 98a, WR 104, and WR 112. Angular sizes were measured for an additional 8 dusty WR stars using aperture masking interferometry, allowing us to probe characteristics sizes down to ~20 milliarcseconds (~40 AU for typical sources). With angular sizes and specific fluxes, we can directly measure the wavelength-dependent surface brightness and size relations for our sample. We discovered tight correlations of these properties within our sample which could not be explained by simple spherically-symmetric dust shells or even the more realistic ``pinwheel nebula'' (3-D) radiative transfer model, when using optical constants of Zubko. While the tightly-correlated surface brightness relations we uncovered offer compelling indirect evidence of a shared and distinctive dust shell geometry amongst our sample, long-baseline interferometers should target the marginally-resolved objects in our sample in order to conclusively establish the presence or absence of the putative underyling colliding wind binaries thought to produce the dust shells around WC Wolf-Rayets.

We present the results of a radial velocity (RV) survey of 14 brown dwarfs (BDs) and very low-mass (VLM) stars in the Upper Scorpius OB association (UScoOB) and 3 BD candidates in the rho Ophiuchi dark cloud core. We obtained high-resolution echelle spectra at the Very Large Telescope using Ultraviolet and Visual Echelle Spectrograph (UVES) at two different epochs for each object, and measured the shifts in their RVs to identify candidates for binary/multiple systems in the sample. The average time separation of the RV measurements is 21.6d, and our survey is sensitive to the binaries with separation < 0.1 au. We found that 4 out of 17 objects (or 24^{+16}_{-13} per cent by fraction) show a significant RV change in 4-33d time scale, and are considered as binary/multiple `candidates.' We found no double-lined spectroscopic binaries in our sample, based on the shape of cross-correlation curves. The RV dispersion of the objects in UScoOB is found to be very similar to that of the BD and VLM stars in Chamaeleon I (Cha I). We also found the distribution of the mean rotational velocities (v sin i) of the UScoOB objects is similar to that of the Cha I, but the dispersion of v sin i is much larger than that of the Cha I objects.

We present radiative transfer models of the circumstellar environment of classical T Tauri stars, concentrating on the formation of the H-alpha emission. The wide variety of line profiles seen in observations are indicative of both inflow and outflow, and we therefore employ a circumstellar structure that includes both magnetospheric accretion and a disc wind. We perform systematic investigations of the model parameters for the wind and the magnetosphere to search for possible geometrical and physical conditions which lead to the types of profiles seen in observations. We find that the hybrid models can reproduce the wide range of profile types seen in observations, and that the most common profile types observed occupy a large volume of parameter space. Conversely, the most infrequently observed profile morphologies require a very specific set of models parameters. We find our model profiles are consistent with the canonical value of the mass-loss rate to mass-accretion rate ratio (mu=0.1) found in earlier magneto-hydrodynamic calculations and observations, but the models with 0.05<mu<0.2 are still in accord with observed H-alpha profiles. We investigate the wind contribution to the line profile as a function of model parameters, and examine the reliability of H-alpha as a mass accretion diagnostic. Finally, we examine the H-alpha spectroscopic classification used by Reipurth et. al, and discuss the basic physical conditions that are required to reproduce the profiles in each classified type.

The inner 0.1 AU around accreting T Tauri stars hold clues to many physical processes that characterize the early evolution of solar-type stars. The accretion-ejection connection takes place at least in part in this compact magnetized region around the central star, with the inner disk edge interacting with the star's magnetosphere thus leading simultaneously to magnetically channeled accretion flows and to high velocity winds and outflows. The magnetic star-disk interaction is thought to have strong implications for the angular momentum evolution of the central system, the inner structure of the disk, and possibly for halting the migration of young planets close to the stellar surface. We review here the current status of magnetic field measurements in T Tauri stars, the recent modeling efforts of the magnetospheric accretion process, including both radiative transfer and multi-D numerical simulations, and summarize current evidence supporting the concept of magnetically-channeled accretion in young stars. We also discuss the limits of the models and highlight observational results which suggest that the star-disk interaction is a highly dynamical and time variable process in young stars.

VITRUV is a second generation spectro-imager for the PRIMA enabled Very Large Telescope Interferometer. By combining simultaneously up to 8 telescopes VITRUV makes the VLTI up to 6 times more efficient. This operational gain allows two novel scientific methodologies: 1) massive surveys of sizes; 2) routine interferometric imaging. The science cases presented concentrate on the qualitatively new routine interferometric imaging methodology. The science cases are not exhaustive but complementary to the PRIMA reference mission. The focus is on: a) the close environment of young stars probing for the initial conditions of planet formation and disk evolution; b) the surfaces of stars tackling dynamos, activity, pulsation, mass-loss and evolution; c) revealing the origin of the extraordinary morphologies of Planetary Nebulae and related stars; d) studying the accretion-ejection structures of stellar black-holes (microquasars) in our galaxy; e) unveiling the different interacting components (torus, jets, BLRs) of Active Galactic Nuclei; and f) probing the environment of nearby supermassive black-holes and relativistic effects in the Galactic Center black-hole.

We present the results of an optical spectroscopic study of the massive Wolf-Rayet binary WR137. These data cover the dust-formation maximum in 1997. Combining all available measurements of radial velocities, we derive, for the first time, a spectroscopic orbit with period 13.05 +/- 0.18 years. The resulting masses, adopting i=67degrees, are M(O)= 20 +/- 2 Mo and M(WR)= 4.4+/- 1.5 Mo. These appear, respectively, around normal and on the low side for the given spectral types. Analysis of the intense multi-site spectroscopic monitoring in 1999 shows that the CIII5696 and CIV5802/12 lines have the highest intrinsic variability levels. The periodogram analysis yields a small-amplitude modulation in the absorption troughs of the CIV5802/12 and HeI5876 lines with a period of 0.83 days, which could be related either to pulsations or large-scale rotating structures as seen in the WN4 star EZ Canis Majoris (WR6).Wavelet analysis of the strong emission lines of CIII5696 and CIV5802/12 enabled us to isolate and follow for several hours small structures associated with density enhancements within the wind of the Wolf-Rayet star. Cross-correlating the variability patterns seen in different lines, we find a weak but significant correlation between the varability in emission lines with different ionization potential, i.e. in lines formed at different distances from the WR stellar core. Adopting a beta wind-velocity law, from the motion of individual subpeaks we find beta around 5, which is significantly larger than the canonical value of 1 usually found in O-star winds.

We present measurements of the longitudinal magnetic field in the circumstellar environment of seven classical T Tauri stars. The measurements are based on high-resolution circular spectropolarimetry of the He I 5876 emission line, which is thought to form in accretion streams controlled by a stellar magnetosphere. We detect magnetic fields in BP Tau, DF Tau and DN Tau, and detect statistically significant fields in GM Aur and RW Aur A at one epoch but not at others. We detect no field for DG Tau and GG Tau, with the caveat that these objects were observed at one epoch only. Our measurements for BP Tau and DF Tau are consistent, both in terms of sign and magnitude, with previous studies, suggesting that the characteristics of T Tauri magnetospheres are persistent over several years. We observed the magnetic field of BP Tau to decline monotonically over three nights, and have detected a peak field of 4kG in this object, the highest magnetic field yet observed in a T Tauri star. We combine our observations with results from the literature in order to perform a statistical analysis of the magnetospheric fields in BP Tau and DF Tau. Assuming a dipolar field, we determine a polar field of ~3kG and a dipole offset of 40deg for BP Tau, while DF Tau's field is consistent with a polar field of ~-4.5kG and a dipole offset of 10deg. We conclude that many classical T Tauri stars have circumstellar magnetic fields that are both strong enough and sufficiently globally-ordered to sustain large-scale magnetospheric accretion flows.

We present Halpha spectropolarimetry observations of a sample of 10 bright T Tauri stars, supplemented with new Herbig Ae/Be star data. A change in the linear polarisation across Halpha is detected in most of the T Tauri (9/10) and Herbig Ae (9/11) objects, which we interpret in terms of a compact source of line photons that is scattered off a rotating accretion disk. We find consistency between the position angle (PA) of the polarisation and those of imaged disk PAs from infrared and millimetre imaging and interferometry studies, probing much larger scales. For the Herbig Ae stars AB Aur, MWC 480 and CQ Tau, we find the polarisation PA to be perpendicular to the imaged disk, which is expected for single scattering. On the other hand, the polarisation PA aligns with the outer disk PA for the T Tauri stars DR Tau and SU Aur and FU Ori, conforming to the case of multiple scattering. This difference can be explained if the inner disks of Herbig Ae stars are optically thin, whilst those around our T Tauri stars and FU Ori are optically thick. Furthermore, we develop a novel technique that combines known inclination angles and our recent Monte Carlo models to constrain the inner rim sizes of SU Aur, GW Ori, AB Aur, and CQ Tau. Finally, we consider the connection of the inner disk structure with the orientation of the magnetic field in the foreground interstellar medium: for FU Ori and DR Tau, we infer an alignment of the stellar axis and the larger magnetic field direction.

We present time-series echelle spectra of the Pa-beta line of the T Tauri star SU Aur, observed over three consecutive nights. The line shows strong variability (~10 per cent) over the velocity range (100 km/s, 420 km/s) in the red broad absorption component, and weaker variability (~2 per cent) over the velocity range (-200 km/s, 0 km/s) in the blue wing. The variability in the velocity range (-200 km/s, 0 km/s) is correlated with that in (200 km/s, 400 km/s), and the variability in these velocity ranges anti-correlates with that in (0 km/s, 100 km/s). The mean spectrum from the second night shows a suggestion of a blue-shifted absorption component at about -150 km/s, similar to that found in the H-alpha and H-beta lines. We find the position of the subpeak in the red absorption component changes steadily with time, and its motion modulates on half the rotational period. We also find that the modulation of the line equivalent width is associated with a half and a third of the rotational period, which is consistent with the surface Doppler images of SU Aur. Radiative transfer models of a rotationally modulated Pa-beta line, produced in the shock-heated magnetospheric accretion flow, are also presented. Models with a magnetic dipole offset reproduce the overall characteristics of the observed line variability, including the line equivalent width and the motion of the subpeak in the red absorption trough.

We present hydrogen emission line profile models of magnetospheric accretion onto Classical T Tauri stars. The models are computed under the Sobolev approximation using the three-dimensional Monte Carlo radiative-transfer code TORUS. We have calculated four illustrative models in which the accretion flows are confined to azimuthal curtains - a geometry predicted by magneto-hydrodynamical simulations. Properties of the line profile variability of our models are discussed, with reference to dynamic spectra and cross-correlation images. We find that some gross characteristics of observed line profile variability are reproduced by our models, although in general the level of variability predicted is larger than that observed. We conclude that this excessive variability probably excludes dynamical simulations that predict accretion flows with low degrees of axisymmetry.

We have conducted a programme to determine the fundamental parameters of a substantial number of eclipsing binaries of spectral types O and B in the Small Magellanic Cloud. New spectroscopic data, obtained with the two-degree-field multi-object spectrograph on the 3.9-m Anglo-Australian Telescope, have been used in conjunction with photometry from the Optical Gravitational Lens Experiment (OGLE-II) database of SMC eclipsing binaries. Previously we reported results for 10 systems; in this second and concluding paper we present spectral types, masses, radii, temperatures, surface gravities and luminosities for the components of a further 40 binaries. The full sample of 50 OB-type eclipsing systems is the largest single set of fundamental parameters determined for high-mass binaries in any galaxy. We find that 21 of the systems studied are in detached configurations, 28 are in semi-detached post-mass-transfer states, and one is a contact binary. Each system provides a primary distance indicator. We find a mean distance modulus to the SMC of 18.91+/-0.03+/-0.1 (internal and external uncertainties; D=60.6+/-1.0 kpc). This value represents one of the most precise available determinations of the distance to the SMC.

Spectropolarimetry is a photon-hungry technique that will reach fruition in the 8-m telescope age. Here I summarize some of the stellar spectropolarimetric research that my collaborators and I have undertaken, with particular emphasis on the circumstellar environment of massive stars, symbiotic binaries, and star formation.

We predict polarimetric line profiles for scattering off rotating disks using a Monte Carlo technique. We have discovered that there is a marked difference between scattering of line emission by a disk that reaches the stellar surface, and a disk with an inner hole. For the case with an inner hole, we find single position-angle rotations, similar to those predicted by analytic models. For the case of an undisrupted disk, we find double rotations in the position angle -- an effect not reported before. We show that this new effect is due to the finite-sized star interacting with the disk's rotational velocity field. Since a gradual increase of the hole size transforms the double rotations smoothly back into single ones -- as the line emission object approaches that of a point source -- our models demonstrate the diagnostic potential of line polarimetry in determining not only the disk inclination, but also the size of the disk inner hole. Our models are generic, and relevant to any line emitting object that is embedded in a rotating disk. Examples are: Herbig Ae stars, T Tauri stars, other young stars, early-type stars with disks, post-AGB stars, cataclysmic variables and other binary systems, as well as extra-galactic objects, such as the disks around super-massive black holes.

We present a catalogue of new spectral types for hot, luminous stars in the Small Magellanic Cloud. The catalogue contains 4161 objects, giving an order of magnitude increase in the number of SMC stars with published spectroscopic classifications. The targets are primarily B- and A-type stars (2862 and 853 objects respectively), with 1 Wolf-Rayet, 139 O-type, and 306 FG stars, sampling the main sequence to ~mid-B. The selection and classification criteria are described, and objects of particular interest are discussed, including UV-selected targets from the Ultraviolet Imaging Telescope (UIT) experiment, Be and B[e] stars, `anomalous A supergiants', and composite-spectrum systems. We examine the incidence of Balmer-line emission, and the relationship between H-gamma equivalent width and absolute magnitude for BA stars.

We present emission line polarimetry data and modelling relevant to the circumstellar geometry and kinematics around pre-main sequence stars. For a sample of both Herbig Ae/Be stars and T Tauri stars, we find that most show polarization changes across Halpha, implying that flattened structures are common on the smallest scales -- and over a range of stellar masses. We also present Monte Carlo calculations of spectral line profiles scattered in rotating accretion disks. We consider both the case of a central star that emits line photons uniformly, as well as via hot spots. Intriguingly, the switch between a uniform point source and a finite-sized star results in a marked difference in the position angle variation across the line. Our models demonstrate the diagnostic potential of line polarimetry in determining the disk inclination and the size of the inner hole -- a spatial scale no other technique currently accesses.

The aim of this article is to demonstrate the useful role that can be played by spectropolarimetric observations of young and evolved emission line stars that analyse the linearly polarized component in their spectra. At the time of writing, this demonstration has to be made on the basis of optical data since there is no common-user infrared facility, in operation, that offers the desired combination of spectral resolution and sensitivity. Here we focus on what can be learned from linear spectropolarimetry alone at reasonably high spectral resolution and at $10^3 < $S/N$ < 10^4$. And we remind that the near infrared (1--2 micron) has the potential to out-perform the optical as a domain to work in because of the greatly reduced interstellar obscuration at these wavelengths. This point has been reached at a time when theory, exploiting flexible Monte Carlo methods, is fast becoming a powerful tool. In short we have the complex phenomena, and the rise of the modelling capability to match -- good data are the missing link.

We present three-dimensional Monte Carlo radiative transfer models of a very young (<10^5 years old) low mass (50 Msun) stellar cluster containing 23 stars and 27 brown dwarfs. The models use the density and the stellar mass distributions from the large-scale smoothed particle hydrodynamics (SPH) simulation of the formation of a low-mass stellar cluster by Bate, Bonnell and Bromm. We compute the spectral energy distributions of the whole cluster and the individual objects, simulated far-infrared Spitzer Space Telescope (SST) images and construct colour-colour diagrams (near-infrared HKL and SST mid-infrared bands). The presence of accretion discs around the light sources influences the morphology of the dust temperature structure on a large scale (up to a several 10^4 au). A considerable fraction of the interstellar dust is underheated compared to a model without the accretion discs because the radiation from the light sources is blocked/shadowed by the discs. The spectral energy distribution (SED) of the model cluster with accretion discs shows excess emission at lambda = 3-30 micron and lambda > 500 micron compared to that without accretion discs. Our model with accretion discs around each object shows a similar distribution of spectral index (2.2-20 micron) values (i.e. Class 0--III sources) as seen in the Rho Ophiuchus cloud. We confirm that the best diagnostics for identifying objects with accretion discs are mid-infrared (lambda = 3-10 micron) colours (e.g. SST IRAC bands) rather than HKL colours.

We present radiative-transfer modelling of the dusty spiral Pinwheel Nebula observed around the Wolf-Rayet/OB-star binary WR104. The models are based on the three-dimensional radiative-transfer code TORUS, modified to include an adaptive mesh that allows us to adequately resolve both the inner spiral turns (sub-AU scales) and the outer regions of the nebula (distances of 10^4 AU from the central source). The spiral model provides a good fit to both the spectral energy distribution and Keck aperture masking interferometry, reproducing both the maximum entropy recovered images and the visibility curves. We deduce a dust creation rate of 8+-1 x 10^{-7} solar masses per year, corresponding to approximately 2% by mass of the carbon produced by the Wolf-Rayet star. Simultaneous modelling of the imaging and spectral data enables us to constrain both the opening-angle of the wind-wind collision interface and the dust grain size. We conclude that the dust grains in the inner part of the Pinwheel nebula are small (~100A), in agreement with theoretical predictions, although we cannot rule out the presence of larger grains (~1 micron) further from the central binary. The opening angle of the wind-wind collision interface appears to be about 40 degrees, in broad agreement with the wind parameters estimated for the central binary. We discuss the success and deficiencies of the model, and the likely benefits of applying similar techniques to the more the more complex nebulae observed around other WR/O star binaries.

We have discovered T Tauri stars which show startling spectral variability between observations seperated by 20 years. In spectra published by Bouvier & Appenzeller (1992) these objects showed very weak H-alpha emission, broad CaII absorption and so called ``composite spectra'', where the spectral type inferred from the blue region is earlier than that inferred from the red. We present here new spectroscopy which shows that all four stars now exhibit strong H-alpha emission, narrow CaII emission and a spectral type which is consistent at all wavelengths. We propose a scheme to understand these changes whereby the composite spectra of these stars can be explained by a period of active accretion onto the central, young star. In this scheme the composite spectrum consists of a contribution from the stellar photosphere and a contribution from a hot, optically thick, accretion component. The optically thick nature of the accretion flow explains the weakness of the H-alpha emission during this phase. Within this scheme, the change to a single spectral type at all wavelengths and emergence of strong H-alpha emission are consistent with the accretion columns becoming optically thin, as the accretion rate drops. There is a strong analogy here with the dwarf novae class of interacting binaries, which show similar behaviour during the decline from outbursts of high mass-transfer rate. The most important consequence of this interpretation is that these objects bring into question the association of Weak-Line T Tauri stars (WTTs) with non-accreting or discless objects. In light of this result we consider the justification for this paradigm.

We present Halpha spectropolarimetry of a large sample of pre-main sequence (PMS) stars of low and intermediate mass, and argue that the technique is a powerful tool in studying the circumstellar geometry around these objects. For the intermediate mass (2 -- 15 Msun) Herbig Ae/Be stars we find that 16 out of 23 show a line effect, which immediately implies that flattening is common among these objects. Furthermore, we find a significant difference in Halpha spectropolarimetry behaviour between the Herbig Be and Ae groups. For the Herbig Be stars, the concept of an electron scattering disc is shown to be a useful concept to explain the depolarizations seen in this spectral range. At lower masses, more complex Halpha polarimetry behaviour starts to appear. The concept of a compact source of Halpha emission that is formed close to the stellar surface, for instance by hot spots due to magnetospheric accretion, is postulated as a working hypothesis to qualitatively explain the Halpha spectropolarimetry behaviour around Herbig Ae and lower mass (M < 2 Msun) T Tauri stars. The striking resemblance in spectropolarimetric behaviour between the T Tauri star RY Tau and the Herbig Ae stars suggests a common origin of the polarized line photons, and hints that low and higher mass pre-main sequence stars may have more in common than had hitherto been suspected.

We present linear Halpha spectropolarimetry of the classical T Tauri star RY Tau. A change in the polarization percentage and position angle across Halpha is detected, which suggests that line photons are scattered in a rotating disc. We derive the position angle from the slope of the loop in the (Q,U) diagram and find it to be 146 +/- 3 degrees. This is perpendicular to the position angle of the disc of 48 +/- 5 degrees as deduced from submillimeter imaging by Koerner & Sargent (1995). This finding is consistent, as scattering off the imaged millimeter disc is expected to yield a polarization signature in a direction that is rotated by 90 degrees from this disc. The observed spectropolarimetric behaviour of RY Tau is reminiscent of that seen in a large group of Herbig Ae stars, suggesting a common circumstellar origin of the polarized photons.

We present the first results of an observational programme to measure the fundamental parameters of over 100 eclipsing binaries in the Small Magellanic Cloud (SMC). The spectroscopic data were obtained by using the 2dF multi-object spectrograph on the 3.9-m Anglo-Australian Telescope, and are used in conjunction with photometry from the OGLE database of SMC eclipsing binaries. Ten systems are discussed in this first paper. Three are detached early-B binaries, six are in a semi-detached configuration, and one is in a marginal contact state. We conclude that the semi-detached systems are undergoing the slow mass-transfer phase of case-A binary evolution, in which the mass donor has reached its Roche lobe while still on the main sequence. Each system provides a primary distance indicator. By constructing a new calibration between spectral type and temperature for O and early B stars, we find a mean distance modulus to the SMC of 18.89 +- 0.04 (statistical) +- 0.10 (systematic). This value represents one of the most precise determinations to date of the distance to the SMC.

We present Halpha spectropolarimetry observations of a sample of 23 Herbig Ae/Be stars. A change in the linear polarisation across Halpha is detected in a large fraction of the objects, which indicates that the regions around Herbig stars are flattened (disc-like) on small scales. A second outcome of our study is that the spectropolarimetric signatures for the Ae stars differ from those of the Herbig Be stars, with characteristics changing from depolarisation across Halpha in the Herbig Be stars, to line polarisations in the Ae group. The frequency of depolarisations detected in the Herbig Be stars (7/12) is particularly interesting as, by analogy to classical Be stars, it may be the best evidence to date that the higher mass Herbig stars are surrounded by flattened structures. For the Herbig Ae stars, 9 out of 11 show a line polarisation effect that can be understood in terms of a compact Halpha emission that is itself polarised by a rotating disc-like circumstellar medium. The spectropolarimetric difference between the Herbig Be and Ae stars may be the first indication that there is a transition in the Hertzsprung-Russell Diagram from magnetic accretion at spectral type A to disc accretion at spectral type B. Alternatively, the interior polarised line emission apparent in the Ae stars may be masked in the Herbig Be stars due to their higher levels of Halpha emission.

H_alpha spectropolarimetry on Herbig Ae/Be stars shows that the innermost regions of intermediate mass (2 -- 15 M_sun) Pre-Main Sequence stars are flattened. This may be the best evidence to date that the higher mass Herbig Be stars are embedded in circumstellar discs. A second outcome of our study is that the spectropolarimetric signatures for the lower mass Herbig Ae stars differ from those of the higher mass Herbig Be stars. Depolarisations across H_alpha are observed in the Herbig Be group, whereas line polarisations are common amongst the Herbig Ae stars in our sample. These line polarisation effects can be understood in terms of a compact H_alpha source that is polarised by a rotating disc-like configuration. The difference we detect between the Herbig Be and Ae stars may be the first indication that there is a transition in the Hertzsprung-Russell Diagram from magnetic accretion at spectral type A to disc accretion at spectral type B. However, it is also possible that the compact polarised line component, present in the Herbig Ae stars, is masked in the Herbig Be stars due to their higher levels of H_alpha emission.

We now understand how low mass stars evolve through cloud collapse and disc accretion, but whether higher mass stars are also disc accretors is as yet unknown. Spectropolarimetry observations can help in answering this basic question, as they probe the first few stellar radii around young stars.