Photoevaporation driven by high energy radiation from the central star plays an important role in the evolution of protoplanetary discs. Photoevaporative winds have been unambiguously detected through blue-shifted emission lines, but their detailed properties remain uncertain. Here we present a new empirical approach to make observational predictions of these thermal winds, seeking to fill the gap between theory and observations. We use a self-similar model of an isothermal wind to compute line profiles of several characteristic emission lines (the [Ne${\rm{\scriptsize II}}$] line at 12.81 $\mu$m, and optical forbidden lines such as [O${\rm{\scriptsize I}}$] 6300 $\mathring{A}$ and [S${\rm{\scriptsize II}}$] 4068/4076 $\mathring{A}$), studying how the lines are affected by parameters such as the gas temperature, disc inclinations, and density profile. Our model successfully reproduces blue-shifted lines with $v_{\rm peak} \lesssim 10$ km/s, which decrease with increasing disc inclination. The line widths increase with increasing disc inclinations and range from $\Delta v \sim 15-30$ km/s. The predicted blue-shifts are mostly sensitive to the gas sound speed. The observed [Ne${\rm{\scriptsize II}}$] line profiles are consistent with a thermal wind and point towards a relatively high sound speed, as expected for EUV photoevaporation. However, the observed [O${\rm{\scriptsize I}}$] line profiles require lower temperatures, as expected in X-ray photoevaporation, and show a wider scatter that is difficult to reconcile with a single wind model; it seems likely that these lines trace different components of a multi-phase wind. We also note that the spectral resolution of current observations remains an important limiting factor in these studies, and that higher resolution spectra are required if emission lines are to further our understanding of protoplanetary disc winds.

Interferometric observations of the mm dust distribution in protoplanetary discs are now showing a ubiquity of annular gap and ring substructures. Their identification and accurate characterization is critical to probing the physical processes responsible. We present Frankenstein (frank), an open source code that recovers axisymmetric disc structures at sub-beam resolution. By fitting the visibilities directly, the model reconstructs a disc's 1D radial brightness profile nonparametrically using a fast (<~1 min) Gaussian process. The code avoids limitations of current methods that obtain the radial brightness profile by either extracting it from the disc image via nonlinear deconvolution at the cost of reduced fit resolution, or by assumptions placed on the functional forms of disc structures to fit the visibilities parametrically. We use mock ALMA observations to quantify the method's intrinsic capability and its performance as a function of baseline-dependent signal-to-noise. Comparing the technique to profile extraction from a CLEAN image, we motivate how our fits accurately recover disc structures at a sub-beam resolution. Demonstrating the model's utility in fitting real high and moderate resolution observations, we conclude by proposing applications to address open questions on protoplanetary disc structure and processes.

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $\alpha/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($\alpha/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $\alpha$ turbulent parameter ($\alpha \sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

Annular structures in proto-planetary discs, such as gaps and rings, are now ubiquitously found by high-resolution ALMA observations. Under the hypothesis that they are opened by planets, in this paper we investigate how the minimum planet mass needed to open a gap varies across different stellar host masses and distances from the star. The dependence on the stellar host mass is particularly interesting because, at least in principle, gap opening around low mass stars should be possible for lower mass planets, giving us a look into the young, low mass planet population. Using dusty hydrodynamical simulations, we find however the opposite behaviour, as a result of the fact that discs around low mass stars are geometrically thicker: gap opening around low mass stars can require more massive planets. Depending on the theoretical isochrone employed to predict the relationship between stellar mass and luminosity, the gap opening planet mass could also be independent of stellar mass, but in no case we find that gap opening becomes easier around low mass stars. This would lead to the expectation of a {\it lower} incidence of such structures in lower mass stars, since exoplanet surveys show that low mass stars have a lower fraction of giant planets. More generally, our study enables future imaging observations as a function of stellar mass to be interpreted using information on the mass vs. luminosity relations of the observed samples.

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time.

In this paper we analyse giant gap-opening planet migration in protoplanetary discs, focusing on the type II migration regime. According to standard type II theory, planets migrate at the same rate as the gas in the disc, as they are coupled to the disc viscous evolution; however, recent studies questioned this paradigm, suggesting that planets migrate faster than the disc material. We study the problem through 2D long-time simulations of systems consistent with type II regime, using the hydrodynamical grid code FARGO3D. Even though our simulations confirm the presence of an initial phase characterised by fast migration, they also reveal that the migration velocity slows down and eventually reaches the theoretical prediction if we allow the system to evolve for enough time. We find the same tendency to evolve towards the theoretical predictions at later times when we analyse the mass flow through the gap and the torques acting on the planet. This transient is related to the initial conditions of our (and previous) simulations, and is due to the fact that the shape of the gap has to adjust to a new profile, once the planet is set into motion. Secondly, we test whether the type II theory expectation that giant planet migration is driven by viscosity is consistent with our simulation by comparing simulations with the same viscosity and different disc mass (or viceversa). We find a good agreement with the theory, since when the discs are characterised by the same viscosity, the migration properties are the same.

Protoplanetary discs form and evolve in a wide variety of stellar environments and are accordingly exposed to a wide range of ambient far ultraviolet (FUV) field strengths. Strong FUV fields are known to drive vigorous gaseous flows from the outer disc. In this paper we conduct the first systematic exploration of the evolution of the solid component of discs subject to external photoevaporation. We find that the main effect of photoevaporation is to reduce the reservoir of dust at large radii and this leads to more efficient subsequent depletion of the disc dust due to radial drift. Efficient radial drift means that photoevaporation causes no significant increase of the dust to gas ratio in the disc. We show that the disc lifetime in both dust and gas is strongly dependent on the level of the FUV background and that the relationship between these two lifetimes just depends on the Shakura-Sunyaev $\alpha$ parameter, with the similar lifetimes observed for gas and dust in discs pointing to higher $\alpha$ values ($\sim 10^{-2}$). On the other hand the distribution of observed discs in the plane of disc size versus flux at $850~\mu$m is better reproduced by lower $\alpha$ ($\sim 10^{-3}$). We find that photoevaporation does not assist rocky planet formation but need not inhibit mechanisms (such as pebble accretion at the water snow line) which can be effective sufficiently early in the disc's lifetime (i.e. well within a Myr).

Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.

We examine the distribution of on-sky relative velocities for wide binaries previously assembled from GAIA DR2 data and focus on the origin of the high velocity tail of apparently unbound systems which may be interpreted as evidence for non-Newtonian gravity in the weak field limit. We argue that this tail is instead explicable in terms of a population of hidden triples, i.e. cases where one of the components of the wide binary is itself a close binary unresolved in the GAIA data. In this case the motion of the photocentre of the inner pair relative to its barycentre affects the apparent relative proper motion of the wide pair and can make pairs that are in fact bound appear to be unbound. We show that the general shape of the observed distributions can be reproduced using simple observationally motivated assumptions about the population of hidden triples.

Protoplanetary discs (PPDs) in the Orion Nebula Cluster (ONC) are irradiated by UV fields from the massive star $\theta^1$C. This drives thermal winds, inducing mass loss rates of up to $\dot{M}_\mathrm{wind}\sim 10^{-7}\,M_\odot$/yr in the `proplyds' (ionised PPDs) close to the centre. For the mean age of the ONC and reasonable initial PPD masses, such mass loss rates imply that discs should have been dispersed. However, ~80% of stars still exhibit a NIR excess, suggesting that significant circumstellar mass remains. This `proplyd lifetime problem' has persisted since the discovery of photoevaporating discs in the core of the ONC by O'Dell & Wen (1994). In this work, we demonstrate how an extended period of star formation can solve this problem. Coupling N-body calculations and a viscous disc evolution model, we obtain high disc fractions at the present day. This is partly due to the migration of older stars outwards, and younger stars inwards such that the most strongly irradiated PPDs are also the youngest. We show how the disc mass distribution can be used to test the recent claims in the literature for multiple stellar populations in the ONC. Our model also explains the recent finding that host mass and PPD mass are only weakly correlated, in contrast with other regions of similar age. We conclude that the status of the ONC as the archetype for understanding the influence of environment on planet formation is undeserved; the complex star formation history (involving star formation episodes within ~0.8 Myr of the present day) results in confusing signatures in the PPD population.

The presence of planets in binary systems poses interesting problems for planet formation theories, both in cases where planets must have formed in very compact discs around the individual stars and where they are located near the edge of the stable circumbinary region, where in situ formation is challenging. Dust dynamics is expected to play an important role in such systems, since dust trapping at the inner edge of circumbinary discs could aid in situ formation, but would simultaneously starve the circumstellar discs of the solid material needed to form planets. Here we investigate the dynamics of dust in binary systems using Smooth Particle Hydrodynamics. We find that all our simulations tend towards dust trapping in the circumbinary disc, but the timescale on which trapping begins depends on binary mass ratio ($q$) and eccentricity as well as the angular momentum of the infalling material. For $q \gtrsim 0.1$, we find that dust can initially accrete onto the circumstellar discs, but as the circumbinary cavity grows in radius, dust eventually becomes trapped in the circumbinary disc. For $q = 0.01$, we find that increasing the binary eccentricity increases the time required for dust trapping to begin. However, even this longer timescale is likely to be shorter than the planet formation timescale in the inner disc and is insufficient to explain the observed pre-transitional discs. This indicates that increase in companion eccentricity alone is not enough to allow significant transfer of solids from the outer to the inner disc.

The large majority of protoplanetary disks have very compact ($\lesssim15\,$AU) continuum emission at mm wavelengths. However, high angular resolution observations that resolve these small disks are still lacking, due to their intrinsically fainter emission compared with large bright disks. In this letter, we present $1.3\,$mm ALMA data of the faint ($\sim10\,$mJy) disk orbiting the TTauri star CX Tau at a resolution of $\sim40\,$mas, $\sim5\,$AU in diameter. The mm-dust disk is compact, with a 68$\%$ enclosing flux radius of 14$\,$AU, and the intensity profile exhibits a sharp drop between 10-20$\,$AU, and a shallow tail between 20-40$\,$AU. No clear signatures of substructure in the dust continuum are observed, down to the same sensitivity level of the DSHARP large program. However, the angular resolution does not allow to detect substructures at a scale of the disk aspect ratio in the inner regions. The radial intensity profile resembles well the inner regions of more extended disks imaged at the same resolution in DSHARP, but with no rings present in the outer disk. No inner cavity is detected, even though the disk has been classified as a transition disk from the spectral energy distribution in the near infrared. The emission of $^{12}$CO is much more extended, with a 68$\%$ enclosing flux radius of 75$\,$AU. The large difference of the mm dust and gas extents ($>5$) strongly points to the occurrence of radial drift, and matches well the predictions of theoretical models.

A correlation between proto-planetary disc radii and sub-mm fluxes has been recently reported. In this Letter we show that the correlation is a sensitive probe of grain growth processes. Using models of grain growth and drift, we have shown in a companion paper that the observed disc radii trace where the dust grains are large enough to have a significant sub-mm opacity. We show that the observed correlation emerges naturally if the maximum grain size is set by radial drift, implying relatively low values of the viscous $\alpha$ parameter $ \lesssim 0.001$. In this case the relation has an almost universal normalisation, while if the grain size is set by fragmentation the flux at a given radius depends on the dust-to-gas ratio. We highlight two observational consequences of the fact that radial drift limits the grain size. The first is that the dust masses measured from the sub-mm could be overestimated by a factor of a few. The second is that the correlation should be present also at longer wavelengths (e.g. 3mm), with a normalisation factor that scales as the square of the observing frequency as in the optically thick case.

Proto-planetary disc surveys conducted with ALMA are measuring disc radii in multiple star forming regions. The disc radius is a fundamental quantity to diagnose whether discs undergo viscous spreading, discriminating between viscosity or angular momentum removal by winds as drivers of disc evolution. Observationally, however, the sub-mm continuum emission is dominated by the dust, which also drifts inwards, complicating the picture. In this paper we investigate, using theoretical models of dust grain growth and radial drift, how the radii of dusty viscous proto-planetary discs evolve with time. Despite the existence of a sharp outer edge in the dust distribution, we find that the radius enclosing most of the dust $\textit{mass}$ increases with time, closely following the evolution of the gas radius. This behaviour arises because, although dust initially grows and drifts rapidly onto the star, the residual dust retained on Myr timescales is relatively well coupled to the gas. Observing the expansion of the dust disc requires using definitions based on high fractions of the disc $\textit{flux}$ (e.g. 95 per cent) and very long integrations with ALMA, because the dust grains in the outer part of the disc are small and have a low sub-mm opacity. We show that existing surveys lack the sensitivity to detect viscous spreading. The disc radii they measure do not trace the mass radius or the sharp outer edge in the dust distribution, but the outer limit of where the grains have significant sub-mm opacity. We predict that these observed radii should shrink with time.

Recent observations show that rings and gaps are ubiquitous in protoplanetary discs. These features are often interpreted as being due to the presence of planets; however, the effect of planetary migration on the observed morphology has not been investigated hitherto. In this work we investigate whether multiwavelength mm/submm observations can detect signatures of planet migration, using 2D dusty hydrodynamic simulations to model the structures generated by migrating planets and synthesising ALMA continuum observations at 0.85 and 3 mm. We identify three possible morphologies for a migrating planet: a slowly migrating planet is associated with a single ring outside the planet's orbit, a rapidly migrating planet is associated with a single ring inside the planet's orbit while a planet migrating at intermediate speed generates one ring on each side of the planet's orbit. We argue that multiwavelength data can distinguish multiple rings produced by a migrating planet from other scenarios for creating multiple rings, such as multiple planets or discs with low viscosity. The signature of migration is that the outer ring has a lower spectral index, due to larger dust grains being trapped there. Of the recent ALMA observations revealing protoplanetary discs with multiple rings and gaps, we suggest that Elias 24 is the best candidate for a planet migrating in the intermediate speed regime.

Computing the flow from externally FUV irradiated protoplanetary discs requires solving complicated and expensive photodissociation physics iteratively in conjunction with hydrodynamics. Previous studies have therefore been limited to 1D models of this process. In this paper we compare 2D-axisymmetric models of externally photoevaporating discs with their 1D analogues, finding that mass loss rates are consistent to within a factor four. The mass loss rates in 2D are higher, in part because half of the mass loss comes from the disc surface (which 1D models neglect). 1D mass loss rates used as the basis for disc viscous evolutionary calculations are hence expected to be conservative. We study the anatomy of externally driven winds including the streamline morphology, kinematic, thermal and chemical structure. A key difference between the 1D and 2D models is in the chemical abundances. For instance in the 2D models CO can be dissociated at smaller radial distances from the disc outer edge than in 1D calculations because gas is photodissociated by radiation along trajectories that are assumed infinitely optically thick in 1D models. Multidimensional models will hence be critical for predicting observable signatures of environmentally photoevaporating protoplanetary discs.

In this paper we present new ALMA observations towards the proto-planet hosting transitional disc of Herbig Ae/Be star HD 100546. This includes resolved 1.3 mm continuum, $^{13}$CO and the first detection of C$^{18}$O in this disc, which displays azimuthal asymmetry in regions spatially coincident with structures previously identified in HST images related to spiral arms. The lower limit on the mass of the dust disc is calculated to be 9.6x10$^{-4}$M$_\odot$. A firm lower-limit on the total gas mass calculated from optically thin, mid-plane tracing C$^{18}$O (2-1) emission is 0.018M$_\odot$ assuming ISM abundances. These mass estimates provide an estimate of gas-to-dust ratio in the disc of 19, the ratio will increase if C$^{18}$O is relatively under-abundant in the disc compared to CO and H2. Through deprojection and azimuthal averaging of the image plane we detect 1.3 mm continuum emission out to 290+/-10 au,$^{13}$CO to 390+/-10 au and C$^{18}$O to 300+/-10au. We measure a radially increasing millimetre spectral index between wavelengths of 867$\mu$m and 1.3 mm, which shows that grain sizes increase towards the star, with solid particles growing to cm scales in the inner disc.

Many stars form in regions of enhanced stellar density, wherein the influence of stellar neighbours can have a strong influence on a protoplanetary disc (PPD) population. In particular, far ultraviolet (FUV) flux from massive stars drives thermal winds from the outer edge of PPDs, accelerating disc destruction. In this work, we present a novel technique for constraining the dynamical history of a star forming environment using PPD properties in a strongly FUV irradiated environment. Applying recent models for FUV induced mass loss rates to the PPD population of Cygnus OB2, we constrain how long ago primordial gas was expelled from the region; $ 0.5$ Myr ago if the Shakura & Sunyaev $\alpha$-viscosity parameter is $\alpha = 10^{-2}$ (corresponding to a viscous timescale of $\tau_\mathrm{visc} \approx 0.5$ Myr for a disc of scale radius $40$ au around a $1\, M_\odot$ star). This value of $\alpha$ is effectively an upper limit, since it assumes efficient extinction of FUV photons throughout the embedded phase. With this gas expulsion timescale we are able to produce a full dynamical model that fits kinematic and morphological data as well as disc fractions. We suggest Cygnus OB2 was originally composed of distinct massive clumps or filaments, each with a stellar mass $\sim 10^4 \, M_\odot$. Finally we predict that in regions of efficient FUV induced mass loss, disc mass $M_\mathrm{disc}$ as a function of stellar host mass $m_\mathrm{star}$ follows a power law with $M_\mathrm{disc} \propto m_\mathrm{star}^\beta$, where $\beta \gtrsim 2.7$ (depending on disc initial conditions and FUV exposure). This is steeper than observed correlations in regions of moderate FUV flux ($1 < \beta <1.9$), and offers a promising diagnostic to establish the influence of external photoevaporation in a given region.

We have investigated the properties of gravito-turbulent discs in 3D using high-resolution shearing-box simulations. For large enough domain sizes, $L_y \gtrsim 60H$, the disc settles down into a quasi-steady state, showing no long term trends in properties or variation with box size. For smaller boxes, we find that the azimuthal wavelength of the dominant spiral modes are limited to the domain size. This is associated with a bursty behaviour that differs from the quasi-steady dynamics at larger sizes. We point out that a similar transition may be expected in global simulations at the point where the range of azimuthal wavelengths is limited by the finite disc size. This condition (i.e. when $60 H \sim 2 \pi R$, i.e. $H/R \sim 0.1$) correctly predicts the transition to bursty behaviour previously found in global simulations for disc-to-star mass ratios in excess of 0.25. We recover a transition in the dynamics from two- to three-dimensional behaviour, characterized by a turbulence that becomes more isotropic on small scales. This turbulence likely plays an important role in the evolution of dust in self-gravitating discs, potentially dominating the collision velocity for particles with Stokes number $< 1$. Finally, we consider the stability of gravito-turbulence against fragmentation, finding that discs which cool faster than a few dynamical times fragment immediately, supporting previous results. We also find hints of stochastic fragmentation at marginally longer cooling times, in which a fragment forms from a quasi-equilibrium state. However, this makes little practical difference to region where planet formation via gravitational instability may occur.

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.

Planet migration in protoplanetary discs plays an important role in the longer term evolution of planetary systems, yet we currently have no direct observational test to determine if a planet is migrating in its gaseous disc. We explore the formation and evolution of dust rings - now commonly observed in protoplanetary discs by ALMA - in the presence of relatively low mass (12-60 Earth masses) migrating planets. Through two dimensional hydrodynamical simulations using gas and dust we find that the importance of perturbations in the pressure profile interior and exterior to the planet varies for different particle sizes. For small sizes a dust enhancement occurs interior to the planet, whereas it is exterior to it for large particles. The transition between these two behaviours happens when the dust drift velocity is comparable to the planet migration velocity. We predict that an observational signature of a migrating planet consists of a significant outwards shift of an observed midplane dust ring as the wavelength is increased.

We derive a self-similar description for the 2D streamline topology and flow structure of an axi-symmetric, thermally driven wind originating from a disc in which the density is a power law function of radius. Our scale-free solution is strictly only valid in the absence of gravity or centrifugal support; comparison with 2D hydrodynamic simulations of winds from Keplerian discs however demonstrates that the scale-free solution is a good approximation also in the outer regions of such discs, and can provide a reasonable description even for launch radii well within the gravitational radius of the flow. Although other authors have considered the flow properties along streamlines whose geometry has been specified in advance, this is the first isothermal calculation in which the flow geometry and variation of flow variables along streamlines is determined self-consistently. It is found that the flow trajectory is very sensitive to the power-law index of radial density variation in the disc: the steeper the density gradient, the stronger is the curvature of streamlines close to the flow base that is required in order to maintain momentum balance perpendicular to the flow. Steeper disc density profiles are also associated with more rapid acceleration, and a faster fall-off of density, with height above the disc plane. The derivation of a set of simple governing equations for the flow structure of thermal winds from the outer regions of power law discs offers the possibility of deriving flow observables without having to resort to hydrodynamical simulation.

We present high resolution millimeter continuum imaging of the disc surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disc system. The system has extended mm emission on which are superposed three prominent annular gaps at radii ~ 13, 39 and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disc; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exo-planets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius.Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage.

We present an open access grid of 3930 calculations of externally evaporating protoplanetary discs. This spans a range of disc sizes (1-400AU), disc masses, UV field strengths (10-10$^4$G$_0$) and stellar masses (0.05-1.9M$_\odot$). The grid is publicly available for download, and offers a means of cheaply including external photoevaporation in disc evolutionary calculations. It can also be queried using an online tool for quick estimates of instantaneous mass loss rates (e.g for convenient evaluation of real observed systems). The `FRIED' grid itself illustrates that for discs around stars $\leq0.3$M$_\odot$ external photoevaporation is effective down to small radii ($<50$AU) down to UV fields at least as weak as 10G$_0$. At the other end of the scale, in a $10^4$G$_0$ environment photoevaporation is effective down to 1AU even for stellar masses at least as high as 1.9M$_\odot$. We also illustrate in which regimes CO survives in the photoevaporative outflow for significant mass loss rates; marking a system a good candidate to detect external photoevaporation in weak-intermediate UV environments through sub-Keplerian rotation. Finally we make illustrative mass loss rate estimates for discs in Taurus based on the Guilloteau et al. (2011) star-disc parameters, finding that around half are expected to have both significant mass loss and retain CO in the photoevaporative outflow.

Theory and observations suggest that star formation occurs hierarchically due to the fragmentation of giant molecular clouds. In this case we would expect substructure and enhanced stellar multiplicity in the primordial cluster. This substructure is expected to decay quickly in most environments, however historic stellar encounters might leave imprints in a protoplanetary disc (PPD) population. In a low density environment such as Taurus, tidal tails from violent star-disc or disc-disc encounters might be preserved over time-scales sufficient to be observed. In this work, we investigate the possibility that just such an event occured between HV Tau C (itself a component of a triple system) and DO Tau $\sim 0.1$ Myr ago, as evidenced by an apparent `bridge' structure evident in the $160$ $\mu$m emission. By modelling the encounter using smoothed particle hydrodynamics (SPH) we reproduce the main features of the observed extended structure (`V'-shaped emission pointing west of HV Tau and a tail-like structure extending east of DO Tau). We suggest that HV Tau and DO Tau formed together in a quadruple system on a scale of $\sim 5000$ au ($0.025$ pc).

We present high-resolution ALMA Band 6 and 7 observations of the tidally disrupted protoplanetary disks of the RW Aurigae binary. Our observations reveal the presence of additional tidal streams to the previously observed tidal arm around RW Aur A. The observed configuration of tidal streams surrounding RW Aur A and B is incompatible with a single star--disk tidal encounter, suggesting that the RW Aurigae system has undergone multiple fly-by interactions. We also resolve the circumstellar disks around RW Aur A and B, with CO radii of 58 au and 38 au consistent with tidal truncation, and 2.5 times smaller dust emission radii. The disks appear misaligned by 12$^{\circ}$ or 57$^{\circ}$. Using new photometric observations from the American Association of Variable Star Observers (AAVSO) and All Sky Automated Survey for SuperNovae (ASAS-SN) archives, we have also identified an additional dimming event of the primary that began in late 2017 and is currently ongoing. With over a century of photometric observations, we are beginning to explore the same spatial scales as ALMA.

Most stars form and spend their early life in regions of enhanced stellar density. Therefore the evolution of protoplanetary discs (PPDs) hosted by such stars are subject to the influence of other members of the cluster. Physically, PPDs might be truncated either by photoevaporation due to ultraviolet flux from massive stars, or tidal truncation due to close stellar encounters. Here we aim to compare the two effects in real cluster environments. In this vein we first review the properties of well studied stellar clusters with a focus on stellar number density, which largely dictates the degree of tidal truncation, and far ultraviolet (FUV) flux, which is indicative of the rate of external photoevaporation. We then review the theoretical PPD truncation radius due to an arbitrary encounter, additionally taking into account the role of eccentric encounters that play a role in hot clusters with a 1D velocity dispersion $\sigma_v > 2$ km/s. Our treatment is then applied statistically to varying local environments to establish a canonical threshold for the local stellar density ($n_{c} > 10^4$ pc$^{-3}$) for which encounters can play a significant role in shaping the distribution of PPD radii over a timescale $\sim 3$ Myr. By combining theoretical mass loss rates due to FUV flux with viscous spreading in a PPD we establish a similar threshold for which a massive disc is completely destroyed by external photoevaporation. Comparing these thresholds in local clusters we find that if either mechanism has a significant impact on the PPD population then photoevaporation is always the dominating influence.

Protoplanetary discs around brown dwarfs and very low mass stars offer some of the best prospects for forming Earth-sized planets in their habitable zones. To this end, we study the nature of the disc around the very low mass star V410 X-ray 1, whose SED is indicative of an optically thick and very truncated dust disc, with our modelling suggesting an outer radius of only 0.6 au. We investigate two scenarios that could lead to such a truncation, and find that the observed SED is compatible with both. The first scenario involves the truncation of both the dust and gas in the disc, perhaps due to a previous dynamical interaction or the presence of an undetected companion. The second scenario involves the fact that a radial location of 0.6 au is close to the expected location of the H$_2$O snowline in the disc. As such, a combination of efficient dust growth, radial migration, and subsequent fragmentation within the snowline leads to an optically thick inner dust disc and larger, optically thin outer dust disc. We find that a firm measurement of the CO $J=2$--1 line flux would enable us to distinguish between these two scenarios, by enabling a measurement of the radial extent of gas in the disc. Many models we consider contain at least several Earth-masses of dust interior to 0.6 au, suggesting that V410 X-ray 1 could be a precursor to a system with tightly-packed inner planets, such as TRAPPIST-1.

We study the evolution of protoplanetary discs that would have been precursors of a Trappist-1 like system under the action of accretion and external photoevaporation in different radiation environments. Dust grains swiftly grow above the critical size below which they are entrained in the photoevaporative wind, so although gas is continually depleted, dust is resilient to photoevaporation after only a short time. This means that the ratio of the mass in solids (dust plus planetary) to the mass in gas rises steadily over time. Dust is still stripped early on, and the initial disc mass required to produce the observed $4\,M_{\oplus}$ of Trappist-1 planets is high. For example, assuming a Fatuzzo & Adams (2008) distribution of UV fields, typical initial disc masses have to be $>30\,$per cent the stellar (which are still Toomre $Q$ stable) for the majority of similar mass M dwarfs to be viable hosts of the Trappist-1 planets. Even in the case of the lowest UV environments observed, there is a strong loss of dust due to photoevaporation at early times from the weakly bound outer regions of the disc. This minimum level of dust loss is a factor two higher than that which would be lost by accretion onto the star during 10 Myr of evolution. Consequently even in these least irradiated environments, discs that are viable Trappist-1 precursors need to be initially massive ($>10\,$per cent of the stellar mass).

The majority of stars form in a clustered environment. This has an impact on the evolution of surrounding protoplanetary discs (PPDs) due to either photoevaporation or tidal truncation. Consequently, the development of planets depends on formation environment. Here we present the first thorough investigation of tidally induced angular momentum loss in PPDs in the distant regime, partly motivated by claims in the literature for the importance of distant encounters in disc evolution. We employ both theoretical predictions and dynamical/hydrodynamical simulations in 2D and 3D. Our theoretical analysis is based on that of Ostriker (1994) and leads us to conclude that in the limit that the closest approach distance $x_{min} \gg r$, the radius of a particle ring, the fractional change in angular momentum scales as $(x_{min}/r)^{-5}$. This asymptotic limit ensures that the cumulative effect of distant encounters is minor in terms of its influence on disc evolution. The angular momentum transfer is dominated by the $m=2$ Lindblad resonance for closer encounters and by the $m=1$, $\omega = 0$ Lindblad resonance at large $x_{min}/r$. We contextualise these results by comparing expected angular momentum loss for the outer edge of a PPD due to distant and close encounters. Contrary to the suggestions of previous works we do not find that distant encounters contribute significantly to angular momentum loss in PPDs. We define an upper limit for closest approach distance where interactions are significant as a function of arbitrary host to perturber mass ratio $M_2/M_1$.

For grain growth to proceed effectively and lead to planet formation a number of barriers to growth must be overcome. One such barrier, relevant for compact grains in the inner regions of the disc, is the `bouncing barrier' in which large grains ($\sim$ mm size) tend to bounce off each other rather than sticking. However, by maintaining a population of small grains it has been suggested that cm-size particles may grow rapidly by sweeping up these small grains. We present the first numerically resolved investigation into the conditions under which grains may be lucky enough to grow beyond the bouncing barrier by a series of rare collisions leading to growth (so-called `breakthrough'). Our models support previous results, and show that in simple models breakthrough requires the mass ratio at which high velocity collisions transition to growth instead of causing fragmentation to be low, $\phi \lesssim 50$. However, in models that take into account the dependence of the fragmentation threshold on mass-ratio, we find breakthrough occurs more readily, even if mass transfer is relatively inefficient. This suggests that bouncing may only slow down growth, rather than preventing growth beyond a threshold barrier. However, even when growth beyond the bouncing barrier is possible, radial drift will usually prevent growth to arbitrarily large sizes.

During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.

We present high resolution ($0.11"\times0.06"$) 3mm ALMA observations of the highly inclined transition disk around the star T Cha. Our continuum image reveals multiple dust structures: an inner disk, a spatially resolved dust gap, and an outer ring. When fitting sky-brightness models to the real component of the 3mm visibilities, we infer that the inner emission is compact ($\le1$au in radius), the gap width is between 18-28 au, and the emission from the outer ring peaks at $\sim36$ au. We compare our ALMA image with previously published 1.6$\mu$m VLT/SPHERE imagery. This comparison reveals that the location of the outer ring is wavelength dependent. More specifically, the peak emission of the 3mm ring is at a larger radial distance than that of the 1.6$\mu$m ring, suggesting that millimeter-sized grains in the outer disk are located further away from the central star than micron-sized grains. We discuss different scenarios to explain our findings, including dead zones, star-driven photoevaporation, and planet-disk interactions. We find that the most likely origin of the dust gap is from an embedded planet, and estimate --- for a single planet scenario --- that T Cha's gap is carved by a $1.2M_{jup}$ planet.

The paper presents new high angular resolution ALMA 1.3 mm dust continuum observations of the protoplanetary system AS 209 in the Ophiuchus star forming region. The dust continuum emission is characterized by a main central core and two prominent rings at $r = 75\,$au and $r = 130\,$au intervaled by two gaps at at $r = 62\,$au and $r = 103\,$au. The two gaps have different widths and depths, with the inner one being narrower and shallower. We determined the surface density of the millimeter dust grains using the 3D radiative transfer disk code \textsc{dali}. According to our fiducial model the inner gap is partially filled with millimeter grains while the outer gap is largely devoid of dust. The inferred surface density is compared to 3D hydrodynamical simulations (FARGO-3D) of planet-disk interaction. The outer dust gap is consistent with the presence of a giant planet ($M_{\rm planet} \sim 0.8\,M_{\rm Staturn}$); the planet is responsible for the gap opening and for the pile-up of dust at the outer edge of the planet orbit. The simulations also show that the same planet can give origin to the inner gap at $r = 62\,$au. The relative position of the two dust gaps is close to the 2:1 resonance and we have investigated the possibility of a second planet inside the inner gap. The resulting surface density (including location, width and depth of the two dust gaps) are in agreement with the observations. The properties of the inner gap pose a strong constraint to the mass of the inner planet ($M_{\rm planet} < 0.1\,M_{\rm J}$). In both scenarios (single or pair of planets), the hydrodynamical simulations suggest a very low disk viscosity ($\alpha < 10^{-4}$). Given the young age of the system (0.5 - 1 Myr), this result implies that the formation of giant planets occurs on a timescale of $\lesssim$ 1\,Myr.

A large fraction of stars are in binary systems, yet the evolution of proto-planetary discs in binaries has been little explored from the theoretical side. In this paper we investigate the evolution of the discs surrounding the primary and secondary components of binary systems on the assumption that this is driven by photoevaporation induced by X-rays from the respective star. We show how for close enough separations (20-30 AU for average X-ray luminosities) the tidal torque of the companion changes the qualitative behaviour of disc dispersal from inside out to outside in. Fewer transition discs created by photoevaporation are thus expected in binaries. We also demonstrate that in close binaries the reduction in viscous time leads to accelerated disc clearing around both components, consistent with $\textit{unresolved}$ observations. When looking at the $\textit{differential}$ disc evolution around the two components, in close binaries discs around the secondary clear first due to the shorter viscous timescale associated with the smaller outer radius. In wide binaries instead the difference in photo-evaporation rate makes the secondaries longer lived, though this is somewhat dependent on the assumed scaling of viscosity with stellar mass. We find that our models are broadly compatible with the growing sample of $\textit{resolved}$ observations of discs in binaries. We also predict that binaries have higher accretion rates than single stars for the same disc mass. Thus binaries probably contribute to the observed scatter in the relationship between disc mass and accretion rate in young stars.

We study the sub-Keplerian rotation and dust content of the circumstellar material around the asymptotic giant branch (AGB) star L$_2$ Puppis. We find that the thermal pressure gradient alone cannot explain the observed rotation profile. We find that there is a family of possible dust populations for which radiation pressure can drive the observed sub-Keplerian rotation. This set of solutions is further constrained by the spectral energy distribution (SED) of the system, and we find that a dust-to-gas mass ratio of $\sim10^{-3}$ and a maximum grain size that decreases radially outwards can satisfy both the rotation curve and SED. These dust populations are dynamically tightly coupled to the gas azimuthally. However grains larger than $\sim0.5\mu$m are driven outward radially by radiation pressure at velocities $\sim5$kms$^{-1}$, which implies a dust replenishment rate of $\sim3\times10^{-9}$M$_\odot$yr$^{-1}$. This replenishment rate is consistent with observational estimates to within uncertainties. Coupling between the radial motion of the dust and gas is weak and hence the gas does not share in this rapid outward motion. Overall we conclude that radiation pressure is a capable and necessary mechanism to explain the observed rotation profile of L$_2$ Pup, and offers other additional constraints on the dust properties.

The chemical conditions in the planet forming regions of protoplanetary discs remain difficult to observe directly. Gas accreting from the disc on to the star provides a way to measure the elemental abundances because even refractory species are in an atomic gaseous form. Here we compare the abundance ratios derived from UV lines probing T Tauri accretion streams to simple models of disc evolution. Although the interpretation of line ratios in terms of abundances is highly uncertain, discs with large cavities in mm images tend to have lower Si emission. Since this can naturally be explained by the suppressed accretion of dust, this suggests that abundance variations are at least partially responsible for the variations seen in the line ratios. Our models of disc evolution due to grain growth, radial drift and the flux of volatile species carried as ices on grain surfaces, give rise to a partial sorting of the atomic species based on the volatility of their dominant molecular carriers. This arises because volatiles are left behind at their snow lines while the grains continue to drift. We show that this reproduces the main features seen in the accretion line ratio data, such as C/N ratios that are a few times solar and the correlation between the Si to volatile ratio with mm-flux. We highlight the fact that developing a more robust linkage between line ratios and abundance ratios and acquiring data for larger samples has the potential to cast considerable light on the chemical history of protoplanetary discs.

(Abridged) We explore the role of X-ray photoevaporation in the evolution and dispersal of viscously evolving T-Tauri discs. We show that the X-ray photoevaporation wind rates scale linearly with X-ray luminosity, such that the observed range of X-ray luminosities for solar-type T-Tauri stars (10e28-10e31 erg\s) gives rise to vigorous disc winds with rates of order 10e-10-10e-7 M_sun/yr. We use the wind solutions from radiation-hydrodynamic models, coupled to a viscous evolution model to construct a population synthesis model so that we may study the physical properties of evolving discs and so-called `transition discs'. Current observations of disc lifetimes and accretion rates can be matched by our model assuming a viscosity parameter alpha = 2.5e-3. Our models confirm that X-rays play a dominant role in the evolution and dispersal of protoplanetary discs giving rise to the observed diverse population of inner hole `transition' sources which include those with massive outer discs, those with gas in their inner holes and those with detectable accretion signatures. To help understand the nature of observed transition discs we present a diagnostic diagram based on accretion rates versus inner hole sizes that demonstrate that, contrary to recent claims, many of the observed accreting and non accreting transition discs can easily be explained by X-ray photoevaporation. Finally, we confirm the conjecture of Drake et al. (2009), that accretion is suppressed by the X-rays through `photoevaporation starved accretion' and predict this effect can give rise to a negative correlation between X-ray luminosity and accretion rate, as reported in the Orion data.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

Chemical compositions of giant planets provide a means to constrain how and where they form. Traditionally, super-stellar elemental abundances in giant planets were thought to be possible due to accretion of metal-rich solids. Such enrichments are accompanied by oxygen-rich compositions (i.e. C/O below the disc's value, assumed to be solar, C/O=0.54). Without solid accretion the planets are expected to have sub-solar metallicity, but high C/O ratios. This arises because the solids are dominated by oxygen-rich species, e.g. H$_2$O and CO$_2$, which freeze out in the disk earlier than CO, leaving the gas metal poor but carbon-rich. Here we demonstrate that super-solar metallicities can be achieved by gas accretion alone when growth and radial drift of pebbles are considered in protoplanetary discs. Through this mechanism planets may simultaneously acquire super-solar metallicities and super-solar C/O ratios. This happens because the pebbles transport volatile species inward as they migrate through the disc, enriching the gas at snow lines where the volatiles sublimate. Furthermore, the planet's composition can be used to constrain where it formed. Since high C/H and C/O ratios cannot be created by accreting solids, it may be possible to distinguish between formation via pebble accretion and planetesimal accretion by the level of solid enrichment. Finally, we expect that Jupiter's C/O ratio should be near or above solar if its enhanced carbon abundance came through accreting metal rich gas. Thus Juno's measurement of Jupiter's C/O ratio should determine whether Jupiter accreted its metals from carbon rich gas or oxygen rich solids.

This prospective chapter gives our view on the evolution of the study of circumstellar discs within the next 20 years from both observational and theoretical sides. We first present the expected improvements in our knowledge of protoplanetary discs as for their masses, sizes, chemistry, the presence of planets as well as the evolutionary processes shaping these discs. We then explore the older debris disc stage and explain what will be learnt concerning their birth, the intrinsic links between these discs and planets, the hot dust and the gas detected around main sequence stars as well as discs around white dwarfs.

The young star Elias 2-27 has recently been observed to posses a massive circumstellar disc with two prominent large-scale spiral arms. In this Letter we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modelling, synthetic ALMA imaging and an unsharped masking technique to explore three possibilities for the origin of the observed structures -- an undetected companion either internal or external to the spirals, and a self-gravitating disc. We find that a gravitationally unstable disc and a disc with an external companion can produce morphology that is consistent with the observations. In addition, for the latter, we find that the companion could be a relatively massive planetary mass companion (less than approximately 10 - 13 MJup) and located at large radial distances (between approximately 300 - 700 au). We therefore suggest that Elias 2-27 may be one of the first detections of a disc undergoing gravitational instabilities, or a disc that has recently undergone fragmentation to produce a massive companion.

We model the radiatively driven flow from IM Lup -- a large protoplanetary disc expected to be irradiated by only a weak external radiation field (at least 10$^4$ times lower than the UV field irradiating the Orion Nebula Cluster proplyds). We find that material at large radii ($>400$AU) in this disc is sufficiently weakly gravitationally bound that significant mass loss can be induced. Given the estimated values of the disc mass and accretion rate, the viscous timescale is long ($\sim10$Myr) so the main evolutionary behaviour for the first Myr of the disc's lifetime is truncation of the disc by photoevaporation, with only modest changes effected by viscosity. We also produce approximate synthetic observations of our models, finding substantial emission from the flow which can explain the CO halo observed about IM Lup out to $\geq1000$AU. Solutions that are consistent with the extent of the observed CO emission generally imply that IM Lup is still in the process of having its disc outer radius truncated. We conclude that IM Lup is subject to substantial external photoevaporation, which raises the more general possibility that external irradiation of the largest discs can be of significant importance even in low mass star forming regions.

We explore how measurements of protoplanetary disc masses and accretion rates provided by surveys of star forming regions can be analysed via the dimensionless accretion parameter, which we define as the product of the accretion rate and stellar age divided by the disc mass. By extending and generalising the study of Jones et al (2012), we demonstrate that this parameter should be less than or of order unity for a wide range of evolutionary scenarios, rising above unity only during the final stages of outside in clearing by external photoevaporation. We use this result to assess the reliability of disc mass estimates derived from CO isotopologues and submm continuum emission by examining the distribution of accretion efficiencies in regions which are not subject to external photoevaporation. We find that while dust based mass estimates produce results compatible with theoretical expectations assuming canonical dust/gas ratio, the systematically lower CO based estimates yield accretion efficiencies significantly above unity in contrast with the theory. This finding provides additional evidence that CO based disc masses are an under-estimate, in line with arguments that have been made on the basis of chemical modelling of relatively small samples. On the other hand, we demonstrate that dust based mass estimates are sufficiently accurate to reveal distinctly higher accretion efficiencies in the Trapezium cluster, where this result is expected given the evident importance of external photoevaporation. We therefore propose the dimensionless accretion parameter as a new diagnostic of external photoevaporation in other star forming regions.

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($\Delta$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

Following the recent discovery of the first radial velocity planet in a star still possessing a protoplanetary disc (CI Tau), we examine the origin of the planet's eccentricity (e $\sim 0.3$). We show through long timescale ($10^5$ orbits) simulations that the planetary eccentricity can be pumped by the disc, even when its local surface density is well below the threshold previously derived from short timescale integrations. We show that the disc may be able to excite the planet's orbital eccentricity in $<$ a Myr for the system parameters of CI Tau. We also perform two planet scattering experiments and show that alternatively the observed planet may plausibly have acquired its eccentricity through dynamical scattering of a migrating lower mass planet, which has either been ejected from the system or swallowed by the central star. In the latter case the present location and eccentricity of the observed planet can be recovered if it was previously stalled within the disc's magnetospheric cavity.

There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their "critical radius". This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes.

Observations of exoplanetary spectra are leading to unprecedented constraints on their atmospheric elemental abundances, particularly O/H, C/H, and C/O ratios. Recent studies suggest that elemental ratios could provide important constraints on formation and migration mechanisms of giant exoplanets. A fundamental assumption in such studies is that the chemical composition of the planetary envelope represents the sum-total of compositions of the accreted gas and solids during the formation history of the planet. We investigate the efficiency with which accreted planetesimals ablate in a giant planetary envelope thereby contributing to its composition rather than sinking to the core. From considerations of aerodynamic drag causing `frictional ablation' and the envelope temperature structure causing `thermal ablation', we compute mass ablations for impacting planetesimals of radii 30 m to 1 km for different compositions (ice to iron) and a wide range of velocities and impact angles, assuming spherical symmetry. Icy impactors are fully ablated in the outer envelope for a wide range of parameters. Even for Fe impactors substantial ablation occurs in the envelope for a wide range of sizes and velocities. For example, iron impactors of sizes below ~0.5 km and velocities above ~30 km/s are found to ablate by ~60-80% within the outer envelope at pressures below 10^3 bar due to frictional ablation alone. For deeper pressures (~10^7 bar), substantial ablation happens over a wider range of parameters. Therefore, our exploratory study suggests that atmospheric abundances of volatile elements in giant planets reflect their accretion history during formation.

The mass of gas in protoplanetary discs is a quantity of great interest for assessing their planet formation potential. Disc gas masses are, however, traditionally inferred from measured dust masses by applying an assumed standard gas-to-dust ratio of $g/d=100$. Furthermore, measuring gas masses based on CO observations has been hindered by the effects of CO freeze-out. Here we present a novel approach to study the mid-plane gas by combining C$^{18}$O line modelling, CO snowline observations and the spectral energy distribution (SED) and selectively study the inner tens of au where freeze-out is not relevant. We apply the modelling technique to the disc around the Herbig Ae star HD 163296 with particular focus on the regions within the CO snowline radius, measured to be at 90 au in this disc. Our models yield the mass of C$^{18}$O in this inner disc region of $M_{\text{C}^{18}\text{O}}(<90\,\text{au})\sim 2\times10^{-8}$ M$_\odot$. We find that most of our models yield a notably low $g/d<20$, especially in the disc mid-plane ($g/d<1$). Our only models with a more interstellar medium (ISM)-like $g/d$ require C$^{18}$O to be underabundant with respect to the ISM abundances and a significant depletion of sub-micron grains, which is not supported by scattered light observations. Our technique can be applied to a range of discs and opens up a possibility of measuring gas and dust masses in discs within the CO snowline location without making assumptions about the gas-to-dust ratio.

We explore the relationship between active galactic nuclei and star formation in a sample of 513 optically luminous type 1 quasars up to redshifts of $\sim$4 hosting extremely high star formation rates (SFRs). The quasars are selected to be individually detected by the \textit{Herschel} SPIRE instrument at $> $3$\sigma$ at 250 $\mu$m, leading to typical SFRs of order of 1000 M$_{\odot}$yr$^{-1}$. We find the average SFRs to increase by almost a factor 10 from $z\sim0.5$ to $z\sim3$, mirroring the rise in the comoving SFR density over the same epoch. However, we find that the SFRs remain approximately constant with increasing accretion luminosity for accretion luminosities above 10$^{12}$ L$_{\odot}$. We also find that the SFRs do not correlate with black hole mass. Both of these results are most plausibly explained by the existence of a self-regulation process by the starburst at high SFRs, which controls SFRs on time-scales comparable to or shorter than the AGN or starburst duty cycles. We additionally find that SFRs do not depend on Eddington ratio at any redshift, consistent with no relation between SFR and black hole growth rate per unit black hole mass. Finally, we find that high-ionisation broad absorption line (HiBAL) quasars have indistinguishable far-infrared properties to those of classical quasars, consistent with HiBAL quasars being normal quasars observed along a particular line of sight, with the outflows in HiBAL quasars not having any measurable effect on the star formation in their hosts.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

We present the discovery of seven new proplyds (i.e. sources surrounded by cometary Halpha emission characteristic of offset ionization fronts) in NGC 1977, located about 30' north of the Orion Nebula Cluster at a distance of ~400 pc. Each of these proplyds are situated at projected distances 0.04-0.27 pc from the B1V star 42 Orionis (c Ori), which is the main source of UV photons in the region. In all cases the ionization fronts of the proplyds are clearly pointing toward the common ionizing source, 42 Ori, and 6 of the 7 proplyds clearly show tails pointing away from it. These are the first proplyds to be found around a B star, with previously known examples instead being located around O stars, including those in the Orion Nebula Cluster around $\theta^1$ Ori C. The radii of the offset ionization fronts in our proplyds are between ~200 and 550 AU; two objects also contain clearly resolved central sources that we associate with disks of radii 50-70 AU. The estimated strength of the FUV radiation field impinging on the proplyds is around 10-30 times less than that incident on the classic proplyds in the Orion Nebula Cluster. We show that the observed proplyd sizes are however consistent with recent models for FUV photoevaporation in relatively weak FUV radiation fields.

We have developed a new prior-based source extraction tool, XID+, to carry out photometry in the Herschel SPIRE maps at the positions of known sources. XID+ is developed using a probabilistic Bayesian framework which provides a natural framework in which to include prior information, and uses the Bayesian inference tool, Stan, to obtain the full posterior probability distribution on flux estimates. In this paper, we discuss the details of XID+ and demonstrate the basic capabilities and performance by running it on simulated SPIRE maps resembling the COSMOS field, and comparing to the current prior-based source extraction tool DESPHOT. We show that not only does XID+ perform better on metrics such as flux accuracy and flux uncertainty accuracy, we illustrate how obtaining the posterior probability distribution can help overcome some of the issues inherent with maximum likelihood based source extraction routines. We run XID+ on the COSMOS SPIRE maps from HerMES, using a 24 $\mathrm{\mu m}$ catalogue as a prior and show the marginalised SPIRE colour-colour plot and marginalised contribution to the cosmic infrared background at the SPIRE wavelengths. XID+ is a core tool arising from the Herschel Extragalactic Legacy Project (HELP) and we discuss how additional work within HELP providing prior information on fluxes can and will be utilised. The software is available at https://github.com/H-E-L-P/XID_plus. We also provide the data product for COSMOS. We believe this is the first time that the full posterior probability of galaxy photometry has been provided as a data product.

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

We investigate the minimum planet mass that produces observable signatures in infrared scattered light and submm continuum images and demonstrate how these images can be used to measure planet masses to within a factor of about two. To this end we perform multi-fluid gas and dust simulations of discs containing low mass planets, generating simulated observations at $1.65 \mu$m, $10 \mu$m and $850 \mu$m. We show that the minimum planet mass that produces a detectable signature is $\sim 15 M_\oplus$: this value is strongly dependent on disc temperature and changes slightly with wavelength (favouring the submm). We also confirm previous results that there is a minimum planet mass of $\sim 20 M_\oplus$ that produces a pressure maximum in the disc: only planets above this threshold mass generate a dust trap that can eventually create a hole in the submm dust. Below this mass, planets produce annular enhancements in dust outward of the planet and a reduction in the vicinity of the planet. These features are in steady state and can be understood in terms of variations in the dust radial velocity, imposed by the perturbed gas pressure radial profile, analogous to a traffic jam. We also show how planet masses can be derived from structure in scattered light and sub-mm images. We emphasise that simulations with dust need to be run over thousands of planetary orbits so as to allow the gas profile to achieve a steady state and caution against the estimation of planet masses using gas only simulations.

We have conducted the first comprehensive numerical investigation of the relative velocity distribution of dust particles in self-gravitating protoplanetary discs with a view to assessing the viability of planetesimal formation via direct collapse in such environments. The viability depends crucially on the large sizes that are preferentially collected in pressure maxima produced by transient spiral features (Stokes numbers, $St \sim 1$); growth to these size scales requires that collision velocities remain low enough that grain growth is not reversed by fragmentation. We show that, for a single sized dust population, velocity driving by the disc's gravitational perturbations is only effective for $St > 3$, while coupling to the gas velocity dominates otherwise. We develop a criterion for understanding this result in terms of the stopping distance being of order the disc scale height. Nevertheless, the relative velocities induced by differential radial drift in multi-sized dust populations are too high to allow the growth of silicate dust particles beyond $St \sim 10^{-2}$ or $10^{-1}$ ($10\,\mathrm{cm}$ to $\mathrm{m}$ sizes at $30\,\mathrm{au}$), such Stokes numbers being insufficient to allow concentration of solids in spiral features. However, for icy solids (which may survive collisions up to several $10\,\mathrm{m\,s}^{-1}$), growth to $St \sim 1$ ($10\,\mathrm{m}$ size) may be possible beyond $30\,\mathrm{au}$ from the star. Such objects would be concentrated in spiral features and could potentially produce larger icy planetesimals/comets by gravitational collapse. These planetesimals would acquire moderate eccentricities and remain unmodified over the remaining lifetime of the disc.

We estimate the mass loss rates of photoevaporative winds launched from the outer edge of protoplanetary discs impinged by an ambient radiation field. We focus on mild/moderate environments (the number of stars in the group/cluster is N ~ 50), and explore disc sizes ranging between 20 and 250 AU. We evaluate the steady-state structures of the photoevaporative winds by coupling temperature estimates obtained with a PDR code with 1D radial hydrodynamical equations. We also consider the impact of dust dragging and grain growth on the final mass loss rates. We find that these winds are much more significant than have been appreciated hitherto when grain growth is included in the modelling: in particular, mass loss rates > 1e-8 M_sun/yr are predicted even for modest background field strengths ( ~ 30 G_0) in the case of discs that extend to R > 150 AU. Grain growth significantly affects the final mass loss rates by reducing the average cross section at FUV wavelengths, and thus allowing a much more vigorous flow. The radial profiles of observable quantities (in particular surface density, temperature and velocity patterns) indicate that these winds have characteristic features that are now potentially observable with ALMA. In particular, such discs should have extended gaseous emission that is dust depleted in the outer regions, characterised by a non-Keplerian rotation curve, and with a radially increasing temperature gradient.

The lack of observed transition discs with inner gas holes of radii greater than ~50AU implies that protoplanetary discs dispersed from the inside out must remove gas from the outer regions rapidly. We investigate the role of photoevaporation in the final clearing of gas from low mass discs with inner holes. In particular, we study the so-called "thermal sweeping" mechanism which results in rapid clearing of the disc. Thermal sweeping was originally thought to arise when the radial and vertical pressure scale lengths at the X-ray heated inner edge of the disc match. We demonstrate that this criterion is not fundamental. Rather, thermal sweeping occurs when the pressure maximum at the inner edge of the dust heated disc falls below the maximum possible pressure of X-ray heated gas (which depends on the local X-ray flux). We derive new critical peak volume and surface density estimates for rapid radiative clearing which, in general, result in rapid dispersal happening less readily than in previous estimates. This less efficient clearing of discs by X-ray driven thermal sweeping leaves open the issue of what mechanism can clear gas from the outer disc sufficiently quickly to explain the non-detection of cold gas around weak line T Tauri stars.

Protoplanetary disks dissipate rapidly after the central star forms, on time-scales comparable to those inferred for planet formation. In order to allow the formation of planets, disks must survive the dispersive effects of UV and X-ray photoevaporation for at least a few Myr. Viscous accretion depletes significant amounts of the mass in gas and solids, while photoevaporative flows driven by internal and external irradiation remove most of the gas. A reasonably large fraction of the mass in solids and some gas get incorporated into planets. Here, we review our current understanding of disk evolution and dispersal, and discuss how these might affect planet formation. We also discuss existing observational constraints on dispersal mechanisms and future directions.

Using 2D smoothed particle hydrodynamics, we investigate the distribution of wait times between strong shocks in a turbulent, self-gravitating accretion disc. We show the resulting distributions do not depend strongly on the cooling time or resolution of the disc and that they are consistent with the predictions of earlier work (Young & Clarke 2015; Cossins et al. 2009, 2010). We use the distribution of wait times between shocks to estimate the likelihood of stochastic fragmentation by gradual contraction of shear-resistant clumps on the cooling time scale. We conclude that the stochastic fragmentation mechanism (Paardekooper 2012) cannot change the radius at which fragmentation is possible by more than ~20%, restricting direct gravitational collapse as a mechanism for giant planet formation to the outer regions of protoplanetary discs.

We perform a series of 2D smoothed particle hydrodynamics (SPH) simulations of gas accretion onto binaries via a circumbinary disc, for a range of gas temperatures and binary mass ratios ($q$). We show that increasing the gas temperature increases the accretion rate onto the primary for all values of the binary mass ratio: for example, for $q=0.1$ and a fixed binary separation, an increase of normalised sound speed by a factor of $5$ (from our "cold" to "hot" simulations) changes the fraction of the accreted gas that flows on to the primary from $ 10\%$ to $\sim40\%$. We present a simple parametrisation for the average accretion rate of each binary component accurate to within a few percent and argue that this parametrisation (rather than those in the literature based on warmer simulations) is relevant to supermassive black hole accretion and all but the widest stellar binaries. We present trajectories for the growth of $q$ during circumbinary disc accretion and argue that the period distribution of stellar "twin" binaries is strong evidence for the importance of circumbinary accretion. We also show that our parametrisation of binary accretion increases the minimum mass ratio needed for spin alignment of supermassive black holes to $q \sim 0.4$, with potentially important implications for the magnitude of velocity kicks acquired during black hole mergers.

We present a 'two-fluid' implementation of dust in smoothed particle hydrodynamics (SPH) in the test particle limit. The scheme is able to handle both short and long stopping times and reproduces the short friction time limit, which is not properly handled in other implementations. We apply novel tests to verify its accuracy and limitations, including multi-dimensional tests that have not been previously applied to the drag-coupled dust problem and which are particularly relevant to self-gravitating protoplanetary discs. Our tests demonstrate several key requirements for accurate simulations of gas-dust mixtures. Firstly, in standard SPH particle jitter can degrade the dust solution, even when the gas density is well reproduced. The use of integral gradients, a Wendland kernel and a large number of neighbours can control this, albeit at a greater computational cost. Secondly, when it is necessary to limit the artificial viscosity we recommend using the Cullen & Dehnen (2010) switch, since the alternative, using ${\alpha} \sim 0.1$, can generate a large velocity noise up to ${{\sigma}_v} \lesssim 0.3 c_s$ in the dust particles. Thirdly, we find that an accurate dust density estimate requires $>400$ neighbours, since, unlike the gas, the dust particles do not feel regularization forces. This density noise applies to all particle-based two-fluid implementations of dust, irrespective of the hydro solver and could lead to numerically induced fragmentation. Although our tests show accurate dusty gas simulations are possible, care must be taken to minimize the contribution from numerical noise.

We propose a framework for understanding the fragmentation criterion for self-gravitating discs which, in contrast to studies that emphasise the `gravoturbulent' nature of such discs, instead focuses on the properties of their quasi-regular spiral structures. Within this framework there are two evolutionary paths to fragmentation: i) collapse on the free-fall time, which requires that the ratio of cooling time to dynamical time ($\beta$) $< 3$ and ii) quasistatic collapse on the cooling time at a rate that is sufficiently fast that fragments are compact enough to withstand disruption when they encounter spiral features in the disc. We perform 2D grid simulations which demonstrate numerically converged fragmentation at $\beta < 3$ (in good agreement with Paardekooper et al. (2011) and others) and argue that this is a consequence of the fact that such simulations smooth the gravitational force on the scale $H$, the scale height of the disc. Such simulations thus only allow fragmentation via route i) above since they suppress the quasistatic contraction of fragments on scales $< H$; the inability of fragments to contract to significantly smaller scales then renders them susceptible to disruption at the next spiral arm encounter. On the other hand, 3D simulations indeed show fragmentation at higher $\beta$ via route ii). We derive an analytic prediction of fragmentation by route ii) when $\beta \lesssim 12$, based on the requirement that fragments must contract sufficiently to withstand disruption by spiral arms. We also discuss the necessary numerical requirements on both grid based and SPH codes if they are to model fragmentation via route ii).

RW Aurigae (RW Aur) is a binary star system with a long molecular arm trailing the primary star. Cabrit et al. (2006) noted the resemblance between this extended structure and the tidal arm stripped from the primary star in the simulations of star-disc encounters by Clarke & Pringle (1993). In this paper we use new hydrodynamical models and synthetic observations to fit many of the parameters of RW Aur. Using hydrodynamic models we find that the morphological appearance of RW Aur can be indeed explained by a tidal encounter with the secondary star. We reproduce all the major morphological and kinematic features of the system. Using radiative transfer calculations, we find that synthetic CO and dust continuum observations of our hydrodynamic models agree well with observations. We reproduce all the main features of the line profiles, from emission fluxes to the optical depth of the different components of the system. The agreement between observations and simulations thus lends strong support to the hypothesis of a tidal encounter scenario. Finally, we propose a possible solution for the origin of the dimming of the primary star observed in 2010/2011 by Rodriguez et al. (2013).

We explore an unresolved controversy in the literature about the accuracy of Smoothed Particle Hydrodynamics (SPH) in modeling the accretion of gas onto a binary system, a problem with important applications to the evolution of proto-binaries as well as accreting binary super massive black holes. It has previously been suggested that SPH fails to model the flow of loosely bound material from the secondary to primary Roche lobe and that its general prediction that accretion drives mass ratios upwards is numerically flawed. Here we show with 2D SPH that this flow from secondary to primary Roche lobe is a sensitive function of gas temperature and that this largely explains the conflicting claims in the literature which have hitherto been based on either 'cold' SPH simulations or 'hot' grid based calculations. We present simulations of a specimen 'cold' and 'hot' accretion scenario which are numerically converged and evolved into a steady state. Our analysis of the conservation of the Jacobi integral of accreting particles indicates that our results are not strongly compromised by numerical dissipation. We also explore the low resolution limit and find that simulations where the ratio of SPH smoothing length to disc scale height at the edge of the circumsecondary is less than 1 accurately capture binary accretion rates.

This paper presents a detailed analysis of two-armed spiral structure in a sample of galax- ies from the Spitzer Infrared Nearby Galaxies Survey (SINGS), with particular focus on the relationships between the properties of the spiral pattern in the stellar disc and the global struc- ture and environment of the parent galaxies. Following Paper I we have used a combination of Spitzer Space Telescope mid-infrared imaging and visible multi-colour imaging to isolate the spiral pattern in the underlying stellar discs, and we examine the systematic behaviours of the observed amplitudes and shapes (pitch angles) of these spirals. In general, spiral morphology is found to correlate only weakly at best with morphological parameters such as stellar mass, gas fraction, disc/bulge ratio, and vflat. In contrast to weak correlations with galaxy structure a strong link is found between the strength of the spiral arms and tidal forcing from nearby companion galaxies. This appears to support the longstanding suggestion that either a tidal interaction or strong bar is a necessary condition for driving grand-design spiral structure. The pitch angles of the stellar arms are only loosely correlated with the pitch angles of the corresponding arms traced in gas and young stars. We find that the strength of the shock in the gas and the contrast in the star formation rate are strongly correlated with the stellar spiral amplitude.

We show that if young low mass stars undergo vigorous X-ray driven disc winds, these may be detected in clusters through their interaction with ionising radiation from massive stars. We argue that in the ONC (Orion Nebula Cluster) one should see $\sim$ 10s of `X-ray proplyds' ( objects with optically imaged offset ionisation fronts) in the range $0.3-0.6$pc from $\theta_1$C Ori (the dominant O star in the ONC). Such objects lie outside the central `FUV zone' in the ONC where proplyds are instead well explained by neutral winds driven by external Far Ultraviolet (FUV) emission from $\theta_1$C. We show that the predicted numbers and sizes of X-ray proplyds are compatible with those observed and that this may also explain at least some of the far flung proplyds seen in the Carina nebula. We compare the sizes of observed proplyds outside the FUV region of the ONC with model predictions based on the current observed X-ray luminosities of these sources ( bearing in mind that the current size is actually set by the X-ray luminosity a few hundred years ago, $\sim$ the flow time to the ionisation front). We discuss whether variability on this timescale can plausibly explain proplyd size data on a case by case basis. We also show that the predicted radio free-free emission signature of X-ray proplyds is readily detectable. Monitoring is however required to distinguish such emission from non-thermal emission from active coronae. We also predict that it is only at $> 1$ pc from $\theta_1$C that the free-free emission from such offset ionised structures would be clearly distinguishable from an externally driven ionised disc wind. We argue that the fortuitous proximity of massive stars in the ONC can be used as a beacon to light up internally driven X-ray winds and that this is a promising avenue for observationally testing X-ray photoevaporation.

Large numbers of small opaque dust clouds - termed 'globulettes' by Gahm et al - have been observed in the H II regions surrounding young stellar clusters. With masses typically in the planetary (or low mass brown dwarf) regime, these objects are so numerous in some regions (e.g. the Rosette) that, if only a small fraction of them could ultimately collapse, then they would be a very significant source of free floating planets. Here we review the properties of globulettes and present a theoretical framework for their structure and evolution. We demonstrate that their interior structure is well described by a pressure confined isothermal Bonnor-Ebert sphere and that the observed mass-radius relation (mass approximately proportional to the radius squared) is a systematic consequence of a column density threshold below which components of the globulette are not identified. We also find that globulettes with this interior structure are very stable against collapse within H II regions. We follow Gahm et al in assuming that globulettes are detached from the tips of pillars protruding in from the swept up shell that borders the expanding H II region and produce a model for their dynamics, finding that globulettes will eventually impact the shell. We derive an expression for the time it takes to do so and show that dissipation of energy via dust cooling allows all globulettes to survive this encounter and escape into the wider ISM. Once there the ambient pressure drops and they disperse on timescales around 30-300 kyr and should be observable using ALMA out to distances of order a parsec. Since we find that globulettes are stable, the only route via which they might still form brown dwarfs or planets is during their collision with the shell or some other violent perturbative event.

We have studied in detail a sample of 967 SPIRE sources with 5-sigma detections at 350 and 500 micron and associations with Spitzer-SWIRE 24 micron galaxies in the HerMES-Lockman survey area, fitting their mid- and far-infrared, and submillimetre, SEDs in an automatic search with a set of six infrared templates. For almost 300 galaxies we have modelled their SEDs individually to ensure the physicality of the fits. We confirm the need for the new cool and cold cirrus templates, and also of the young starburst template, introduced in earlier work. We also identify 109 lensing candidates via their anomalous SEDs and provide a set of colour-redshift constraints which allow lensing candidates to be identified from combined Herschel and Spitzer data. The picture that emerges of the submillimetre galaxy population is complex, comprising ultraluminous and hyperluminous starbursts, lower luminosity galaxies dominated by interstellar dust emission, lensed galaxies and galaxies with surprisingly cold (10-13K) dust. 11 % of 500$\mu$m selected sources are lensing candidates. 70 % of the unlensed sources are ultraluminous infrared galaxies and 26 % are hyperluminous. 34 % are dominated by optically thin interstellar dust ('cirrus') emission, but most of these are due to cooler dust than is characteristic of our Galaxy. At the highest infrared luminosities we see SEDs dominated by M82, Arp220 and young starburst types, in roughly equal proportions.

Fomalhaut is a triple system, with all components widely separated (~1E5 au). Such widely separated binaries are thought to form during cluster dissolution, but that process is unlikely to form such a triple system. We explore an alternative scenario, where A and C form as a tighter binary from a single molecular cloud core (with semimajor axis ~1E4 au), and B is captured during cluster dispersal. We use N-body simulations augmented with the Galactic tidal forces to show that such a system naturally evolves into a Fomalhaut-like system in about half of cases, on a timescale compatible with the age of Fomalhaut. From initial non-interacting orbits, Galactic tides drive cycles in B's eccentricity that lead to a close encounter with C. After several close encounters, typically lasting tens of millions of years, one of the stars is ejected. The Fomalhaut-like case with both components at large separations is almost invariably a precursor to the ejection of one component, most commonly Fomalhaut C. By including circumstellar debris in a subset of the simulations, we also show that such an evolution usually does not disrupt the coherently eccentric debris disk around Fomalhaut A, and in some cases can even produce such a disk. We also find that the final eccentricity of the disk around A and the disk around C are correlated, which may indicate that the dynamics of the three stars stirred C's disk, explaining its unusual brightness.

Observations from optical to centimeter wavelengths have demonstrated that multiple systems of two or more bodies is the norm at all stellar evolutionary stages. Multiple systems are widely agreed to result from the collapse and fragmentation of cloud cores, despite the inhibiting influence of magnetic fields. Surveys of Class 0 protostars with mm interferometers have revealed a very high multiplicity frequency of about 2/3, even though there are observational difficulties in resolving close protobinaries, thus supporting the possibility that all stars could be born in multiple systems. Near-infrared adaptive optics observations of Class I protostars show a lower binary frequency relative to the Class 0 phase, a declining trend that continues through the Class II/III stages to the field population. This loss of companions is a natural consequence of dynamical interplay in small multiple systems, leading to ejection of members. We discuss observational consequences of this dynamical evolution, and its influence on circumstellar disks, and we review the evolution of circumbinary disks and their role in defining binary mass ratios. Special attention is paid to eclipsing PMS binaries, which allow for observational tests of evolutionary models of early stellar evolution. Many stars are born in clusters and small groups, and we discuss how interactions in dense stellar environments can significantly alter the distribution of binary separations through dissolution of wider binaries. The binaries and multiples we find in the field are the survivors of these internal and external destructive processes, and we provide a detailed overview of the multiplicity statistics of the field, which form a boundary condition for all models of binary evolution. Finally we discuss various formation mechanisms for massive binaries, and the properties of massive trapezia.

We analyse a hydrodynamical simulation of star formation. Sink particles in the simulations which represent stars show episodic growth, which is presumably accretion from a core that can be regularly replenished in response to the fluctuating conditions in the local environment. The accretion rates follow $\dot{m} \propto m^{2/3}$, as expected from accretion in a gas-dominated potential, but with substantial variations over-laid on this. The growth times follow an exponential distribution which is tapered at long times due to the finite length of the simulation. The initial collapse masses have an approximately lognormal distribution with already an onset of a power-law at large masses. The sink particle mass function can be reproduced with a non-linear stochastic process, with fluctuating accretion rates $\propto m^{2/3}$, a distribution of seed masses and a distribution of growth times. All three factors contribute equally to the form of the final sink mass function. We find that the upper power law tail of the IMF is unrelated to Bondi-Hoyle accretion.

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

(Abridged) We consider the properties of protoplanetary discs that are undergoing inside-out clearing by photoevaporation. In particular, we aim to characterise the conditions under which a protoplanetary disc may undergo `thermal sweeping', a rapid (< 1e4 years) disc destruction mechanism proposed to occur when a clearing disc reaches sufficiently low surface density at its inner edge and where the disc is unstable to runaway penetration by the X-rays. We use a large suite of 1D radiation-hydrodynamic simulations to probe the observable parameter space, which is unfeasible in higher dimensions. These models allow us to determine the surface density at which thermal sweeping will take over the disc's evolution and to evaluate this critical surface density as a function of X-ray luminosity, stellar mass and inner hole radius. We find that this critical surface density scales linearly with X-ray luminosity, increases with inner hole radius and decreases with stellar mass and we develop an analytic model that reproduces these results. This surface density criterion is then used to determine the evolutionary state of protoplanetary discs at the point that they become unstable to destruction by thermal sweeping. We find that transition discs created by photoevaporation will undergo thermal sweeping when their inner holes reach 20-40 AU, implying that transition discs with large holes and no accretion (which were previously a predicted outcome of the later stages of all flavours of photoevaporation model) will not form. We emphasise that the surface density criteria that we have developed apply to all situations where the disc develops an inner hole that is optically thin to X-rays. It thus applies not only to the case of holes originally created by photoevaporation but also to holes formed, for example, by the tidal influence of planets.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

We assume a scenario in which transition discs (i.e. discs around young stars that have signatures of cool dust but lack significant near infra-red emission from warm dust) are associated with the presence of planets (or brown dwarfs). These are assumed to filter the dust content of any gas flow within the planetary orbit and produce an inner `opacity hole'. In order to match the properties of transition discs with the largest (~50 A.U. scale) holes, we place such `planets' at large radii in massive discs and then follow the evolution of the tidally coupled disc-planet system, comparing the system's evolution in the plane of mm flux against hole radius with the properties of observed transition discs. We find that, on account of the high disc masses in these systems, all but the most massive `planets' (100 Jupiter masses) are conveyed to small radii by Type II migration without significant fading at millimetre wavelengths. Such behaviour would contradict the observed lack of mm bright transition discs with small (<10 A.U.) holes. On the other hand, imaging surveys clearly rule out the presence of such massive companions in transition discs. We conclude that this is a serious problem for models that seek to explain transition discs in terms of planetary companions unless some mechanism can be found to halt inward migration and/or suppress mm flux production. We suggest that the dynamical effects of substantial accretion on to the planet/through the gap may offer the best prospect for halting such migration but that further long term simulations are required to clarify this issue.

We investigate the kinematic evolution of planetesimals in self-gravitating discs, combining Smoothed Particle Hydrodynamical (SPH) simulations of the disc gas with a gravitationally coupled population of test particle planetesimals. We find that at radii of 10s of au (which is where we expect planetesimals to be possibly formed in such discs) the planetesimals' eccentricities are rapidly pumped to values $>$ 0.1 within the timescales for which the disc is in the self-gravitating regime. The high resulting velocity dispersion and the lack of planetesimal concentration in the spiral arms means that the collision timescale is very long and that the effect of those collisions that do occur is destructive rather than leading to further planetesimal growth. We also use the SPH simulations to calibrate Monte Carlo dynamical experiments: these can be used to evolve the system over long timescales and can be compared with analytical solutions of the diffusion equation in particle angular momentum space. We find that if planetesimals are only formed in a belt at large radius then there is significant scattering of objects to small radii; nevertheless the majority of planetesimals remain at large radii. If planetesimals indeed form at early evolutionary stages, when the disc is strongly self-gravitating, then the results of this study constrain their spatial and kinematic distribution at the end of the self-gravitating phase.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

Bressert et al. recently showed that the surface density distribution of low-mass, young stellar objects (YSOs) in the solar neighbourhood is approximately log-normal. The authors conclude that the star formation process is hierarchical and that only a small fraction of stars form in dense star clusters. Here, we show that the peak and the width of the density distribution are also what follow if all stars form in bound clusters which are not significantly affected by the presence of gas and expand by two-body relaxation. The peak of the surface density distribution is simply obtained from the typical ages (few Myr) and cluster membership number (few hundred) typifying nearby star-forming regions. This result depends weakly on initial cluster sizes, provided that they are sufficiently dense (initial half mass radius of <0.3 pc) for dynamical evolution to be important at an age of a few Myr. We conclude that the degeneracy of the YSO surface density distribution complicates its use as a diagnostic of the stellar formation environment.

FU Orionis systems are young stars undergoing outbursts of disc accretion and where the optical spectrum contains lines associated with both the disc photosphere and a wind component. Previous observations of the prototype FU Orionis have suggested that the wind lines and the photospheric lines are modulated with periods of 14.54 and 3.54 days respectively (Herbig et al. 2003). We have re-observed the system at higher spectral resolution, by monitoring variations of optical line profiles over 21 nights in 2007 and have found periods of 13.48 and 3.6 days in the wind and disc components consistent with the above: this implies variability mechanisms that are stable over at least a decade. In addition we have found: i) that the variations in the photospheric absorption lines are confined to the blue wing of the line (around -9km/s): we tentatively ascribe this to an orbiting hotspot in the disc which is obscured by a disc warp during its receding phase. ii) The wind period is manifested not only in blue-shifted Halpha absorption, but also in red-shifted emission of Halpha and Hbeta, as well as in blue-shifted absorption of Na I D, Li I and Fe II. iii) We find that the periodic modulation of blue-shifted Halpha absorption at around -100km/s, is phase lagged with respect to variations in the other lines by ~1.8days. This is consistent with a picture in which variations at the wind base first affect chromospheric emission and then low velocity blue-shifted absorption, followed - after a lag equal to the propagation time of disturbances across the wind's acceleration region - by a response in high velocity blue-shifted absorption. Such arguments constrain the size of the acceleration region to ~10^12cm. We discuss possible mechanisms for periodic variations within the innermost 0.1AU of the disc, including the possibility that these variations indicate the presence of an embedded hot Jupiter.

We perform a comparison between two radiative transfer algorithms commonly employed in hydrodynamical calculations of star formation: grey flux limited diffusion and the hybrid scheme, in addition we compare these algorithms to results from the Monte-Carlo radiative transfer code MOCASSIN. In disc like density structures the hybrid scheme performs significantly better than the FLD method in the optically thin regions, with comparable results in optically thick regions. In the case of a forming high mass star we find the FLD method significantly underestimates the radiation pressure by a factor of ~100.

We examine the distribution of transition discs as a function of mm flux. We confirm that as expected in any model in which most primordial discs turn into transition discs and in which mm flux declines with time, transition discs have lower mm fluxes on average than primordial discs. However, we find that the incidence of transition discs does not, as expected, fall monotonically towards large mm fluxes and we investigate the hypothesis that these mm bright transition discs may have a distinct physical origin. We find that mm bright transition discs occupy a separate region of parameter space. Transition discs in the bright mm sub-sample have systematically higher accretion rates and inner hole radii than those in the faint mm sub-sample, along with being systematically weighted to earlier spectral types.

We describe a simple method for estimating the vertical column density in Smoothed Particle Hydrodynamics (SPH) simulations of discs. As in the method of Stamatellos et al. (2007), the column density is estimated using pre-computed local quantities and is then used to estimate the radiative cooling rate. The cooling rate is a quantity of considerable importance, for example, in assessing the probability of disc fragmentation. Our method has three steps: (i) the column density from the particle to the mid plane is estimated using the vertical component of the gravitational acceleration, (ii) the "total surface density" from the mid plane to the surface of the disc is calculated, (iii) the column density from each particle to the surface is calculated from the difference between (i) and (ii). This method is shown to greatly improve the accuracy of column density estimates in disc geometry compared with the method of Stamatellos. On the other hand, although the accuracy of our method is still acceptable in the case of high density fragments formed within discs, we find that the Stamatellos method performs better than our method in this regime. Thus, a hybrid method (where the method is switched in regions of large over-density) may be optimal.

Infant mortality brought about by the expulsion of a star cluster's natal gas is widely invoked to explain cluster statistics at different ages. While a well studied problem, most recent studies of gas expulsion's effect on a cluster have focused on massive clusters, with stellar counts of order $10^4$. Here we argue that the evolutionary timescales associated with the compact low-mass clusters typical of the median cluster in the Solar neighborhood are short enough that significant dynamical evolution can take place over the ages usually associated with gas expulsion. To test this we perform {\it N}-body simulations of the dynamics of a very young star forming region, with initial conditions drawn from a large-scale hydrodynamic simulation of gravitational collapse and fragmentation. The subclusters we analyse, with populations of a few hundred stars, have high local star formation efficiencies and are roughly virialised even after the gas is removed. Over 10 Myr they expand to a similar degree as would be expected from gas expulsion if they were initially gas-rich, but the expansion is purely due to the internal stellar dynamics of the young clusters. The expansion is such that the stellar densities at 2 Myr match those of YSOs in the Solar neighborhood. We argue that at the low-mass end of the cluster mass spectrum, a deficit of clusters at 10s of Myr does not necessarily imply gas expulsion as a disruption mechanism.

We discuss a hydrodynamical model for the dispersal of protoplanetary discs around young, low mass (<1.5 M_sun) stars by photoevaporation from the central object's energetic radiation, which considers the far-ultraviolet as well as the X-ray component of the radiation field. We present analytical scaling relations and derive estimates for the total mass-loss rates, as well as discussing the existence of similarity solutions for flows from primordial discs and discs with inner holes. Furthermore, we perform numerical calculations, which span a wide range of parameter space and allow us to provide accurate scalings of the mass-loss rates with the physical parameters of the systems (X-ray and FUV luminosity, stellar mass, disc mass, disc temperature and inner hole radius). The model suggest that the X-ray component dominates the photoevaporative mass-loss rates from the inner disc. The mass-loss rates have values in the range from 10e-11 to 10e-7 M_sun/yr and scale linearly with X-ray luminosity, with only a weak dependence on the other parameters explored. However, in the case of high FUV to X-ray (L_FUV/L_X>100) luminosity ratios, the FUV constricts the X-ray flow and may dominate the mass-loss. Simulations of low mass discs with inner holes demonstrate a further stage of disc clearing, which we call `thermal sweeping'. This process occurs when the mid-plane pressure drops to sufficiently low values. At this stage a bound, warm, X-ray heated region becomes sufficiently large and unstable, such that the remaining disc material is cleared on approximately dynamical time-scales. This process significantly reduces the time taken to clear the outer regions of the disc, resulting in an expected transition disc population that will be dominated by accreting objects, as indicated by recent observations.

Hydrodynamical simulations of star formation have stimulated a need to develop fast and robust algorithms for evaluating radiative cooling. Here we undertake a critical evaluation of what is currently a popular method for prescribing cooling in SPH simulations, i.e. the polytropic cooling due originally to Stamatellos et al. This method uses the local density and potential to estimate the column density and optical depth to each particle and then uses these quantities to evaluate an approximate expression for the net radiative cooling. We evaluate the algorithm by considering both spherical and disc-like systems with analytic density and temperature structures. In spherical systems, the total cooling rate computed by the method is within around 20 for the astrophysically relevant case of opacity dominated by ice grains and is correct to within a factor of order unity for a range of opacity laws. In disc geometry, however, the method systematically under-estimates the cooling by a large factor at all heights in the disc. For the self-gravitating disc studied, we find that the method under-estimates the total cooling rate by a factor of 200. This discrepancy may be readily traced to the method's systematic over-estimate of the disc column density and optical depth, since (being based only on the local density and potential) it does not take into account the low column density route for photon escape normal to the disc plane. These results raise an obvious caution about the method's use in disc geometry whenever an accurate cooling rate is required, although we note that there are situations where the discrepancies highlighted above may not significantly affect the global outcome of simulations. Finally, we draw attention to our introduction of an analytic self-gravitating disc structure that may be of use in the calibration of future cooling algorithms.

The fraction of star formation that results in bound star clusters is influenced by the density spectrum in which stars are formed and by the response of the stellar structure to gas expulsion. We analyse hydrodynamical simulations of turbulent fragmentation in star-forming regions to assess the dynamical properties of the resulting population of stars and (sub)clusters. Stellar subclusters are identified using a minimum spanning tree algorithm. When considering only the gravitational potential of the stars and ignoring the gas, we find that the identified subclusters are close to virial equilibrium (the typical virial ratio Q_vir~0.59, where virial equilibrium would be Q_vir~0.5). This virial state is a consequence of the low gas fractions within the subclusters, caused by the accretion of gas onto the stars and the accretion-induced shrinkage of the subclusters. Because the subclusters are gas-poor, up to a length scale of 0.1-0.2 pc at the end of the simulation, they are only weakly affected by gas expulsion. The fraction of subclusters that reaches the high density required to evolve to a gas-poor state increases with the density of the star-forming region. We extend this argument to star cluster scales, and suggest that the absence of gas indicates that the early disruption of star clusters due to gas expulsion (infant mortality) plays a smaller role than anticipated, and is potentially restricted to star-forming regions with low ambient gas densities. We propose that in dense star-forming regions, the tidal shocking of young star clusters by the surrounding gas clouds could be responsible for the early disruption. This `cruel cradle effect' would work in addition to disruption by gas expulsion. We suggest possible methods to quantify the relative contributions of both mechanisms.

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

Recent analyses of mass segregation diagnostics in star forming regions invite a comparison with the output of hydrodynamic simulations of star formation. In this work we investigate the state of mass segregation of 'stars' (i.e. sink particles in the simulations) in the case of hydrodynamical simulations which omit feedback. We first discuss methods to quantify mass segregation in substructured regions, either based on the minimum spanning tree (Allison's Lambda), or through analysis of correlations between stellar mass and local stellar surface number densities. We find that the presence of even a single 'outlier' (i.e. a massive object far from other stars) can cause the Allison Lambda method to describe the system as inversely mass segregated, even where in reality the most massive sink particles are overwhelmingly in the centres of the subclusters. We demonstrate that a variant of the Lambda method is less susceptible to this tendency but also argue for an alternative representation of the data in the plane of stellar mass versus local surface number density. The hydrodynamical simulations show global mass segregation from very early times which continues throughout the simulation, being only mildly influenced during sub-cluster merging. We find that up to approx. 2-3% of the "massive" sink particles (m > 2.5 Msun) are in relative isolation because they have formed there, although other sink particles can form later in their vicinity. Ejections of massive sinks from subclusters do not contribute to the number of isolated massive sink particles, as the gravitational softening in the calculation suppresses this process.

We present a new scheme for determining the shape of the size distribution, and its evolution, for collisional cascades of planetesimals undergoing destructive collisions and loss processes like Poynting-Robertson drag. The scheme treats the steady state portion of the cascade by equating mass loss and gain in each size bin; the smallest particles are expected to reach steady state on their collision timescale, while larger particles retain their primordial distribution. For collision-dominated disks, steady state means that mass loss rates in logarithmic size bins are independent of size. This prescription reproduces the expected two phase size distribution, with ripples above the blow-out size, and above the transition to gravity-dominated planetesimal strength. The scheme also reproduces the expected evolution of disk mass, and of dust mass, but is computationally much faster than evolving distributions forward in time. For low-mass disks, P-R drag causes a turnover at small sizes to a size distribution that is set by the redistribution function (the mass distribution of fragments produced in collisions). Thus information about the redistribution function may be recovered by measuring the size distribution of particles undergoing loss by P-R drag, such as that traced by particles accreted onto Earth. Although cross-sectional area drops with 1/age^2 in the PR-dominated regime, dust mass falls as 1/age^2.8, underlining the importance of understanding which particle sizes contribute to an observation when considering how disk detectability evolves. Other loss processes are readily incorporated; we also discuss generalised power law loss rates, dynamical depletion, realistic radiation forces and stellar wind drag.

Wide, fragile binary stellar systems are found in the galactic field, and have recently been noted in the outskirts of expanding star clusters in numerical simulations. Energetically soft, with semi-major axes exceeding the initial size of their birth cluster, it is puzzling how these binaries are created and preserved. We provide an interpretation of the formation of these binaries that explains the total number formed and their distribution of energies. A population of weakly bound binaries can always be found in the cluster, in accordance with statistical detailed balance, limited at the soft end only by the current size of the cluster and whatever observational criteria are imposed. At any given time, the observed soft binary distribution is predominantly a snapshot of a transient population. However, there is a constantly growing population of long-lived soft binaries that are removed from the detailed balance cycle due to the changing density and velocity dispersion of an expanding cluster. The total number of wide binaries that form, and their energy distribution, are insensitive to the cluster population; the number is approximately one per cluster. This suggests that a population composed of many dissolved small-N clusters will more efficiently populate the field with wide binaries than that composed of dissolved large-N clusters. Locally such binaries are present at approximately the 2% level; thus the production rate is consistent with the field being populated by clusters with a median of a few hundred stars rather than a few thousand.

This paper, the first of two, introduces an observational study of spiral structure in galaxies chosen from the SINGS survey. Near infrared (NIR) and optical data are used to produce mass surface density maps, and from these the morphology of the disc is examined. The aim of this work is to characterise the prevalence of spiral structure in this sample and, in the cases where a clear spiral pattern is found, include the findings in a comparative study (reported in paper II). A two-armed (`grand design') spiral pattern is found in approximately half the galaxies studied, including all those that are designated as grand design in the optical, but also including some, but not all, optically flocculent galaxies. It is found that the level of non-axisymmetric structure in the galaxies' mass distributions is only modestly higher in those galaxies that are classified as `grand design' compared with those that are not, implying that non-grand design galaxies possess significant power in higher order modes. There is no evidence that bars preferentially trigger the spirals, but they do appear to stir up non-axisymmetric structure in the disc. In contrast, there is evidence that strong/close tidal interactions with companion galaxies are associated with strong two-armed spiral structure in the infrared, though there are a number of galaxies with relatively weak infrared spiral structure that do not possess such companions.

Stimulated by recent results by Meru and Bate (2010a,b), we revisit the issue of resolution requirements for simulating self-gravitating accretion discs with Smoothed Particle Hydrodynamics (SPH). We show that the results by Meru and Bate (2010a) are consistent with those of Meru and Bate (2010b) if they are both interpreted as driven by resolution effects, therefore implying that the resolution criterion for cooling gaseous discs is a function of the imposed cooling rate. We discuss two possible numerical origins of such dependence, which are both consistent with the limited number of available data. Our results tentatively indicate that convergence for current simulations is being reached for a number of SPH particles approaching 10 millions (for a disc mass of order 10 per cent of the central object mass), which would set the critical cooling time for fragmentation at about $15\Omega^{-1}$, roughly a factor two larger than previously thought. More in general, we discuss the extent to which the large number of recent numerical results are reliable or not. We argue that those results that pertain to the dynamics associated with gravitational instabilities (such as the locality of angular momentum transport, and the relationship between density perturbation and induced stress) are robust, while those pertaining to the thermodynamics of the system (such as the determination of the critical cooling time for fragmentation) can be affected by poor resolution.

The highest-mass stars have the lowest frequency in the stellar IMF, and they are also the most easily observed stars. Thus, the counting statistics for OB stars provide important tests for the fundamental nature and quantitative parameters of the IMF. We first examine some local statistics for the stellar upper-mass limit itself. Then, we examine the parameter space and statistics for extremely sparse clusters that contain OB stars, in the SMC. We find that thus far, these locally observed counting statistics are consistent with a constant stellar upper-mass limit. The sparse OB star clusters easily fall within the parameter space of Monte Carlo simulations of cluster populations. If the observed objects are representative of their cluster birth masses, their existence implies that the maximum stellar mass is largely independent of the parent cluster mass.

We have performed hydrodynamic and radiative transfer calculations of a photoevaporating disc around a Herbig Ae/Be star to determine the evolution and observational impact of dust entrained in the wind. We find that the wind selectively entrains grains of different sizes at different radii resulting in a dust population that varies spatially and increases with height above the disc at radii > 10 AU. This variable grain population results in a 'wingnut' morphology to the dust density distribution. We calculate images of this dust distribution at NIR wavelengths that also show a wingnut morphology at all wavelengths considered. We have also considered the contribution that crystalline dust grains will have in the wind and show that a photoevaporative wind can result in a significant crystallinity fraction at all radii, when the disc is edge-on. However, when the disc's photosphere is unobscured, a photoevaporative wind makes no contribution to the observable crystallinity fraction in the disc. Finally, we conclude that the analysis of extended emission around edge-on discs could provide a new and independent method of testing photoevaporation models.

We investigate the contraction of accreting protoclusters using an extension of n-body techniques that incorporates the accretional growth of stars from the gaseous reservoir in which they are embedded. Following on from Monte Carlo studies by Davis et al., we target our experiments toward populous clusters likely to experience collisions as a result of accretion-driven contraction. We verify that in less extreme star forming environments, similar to Orion, the stellar density is low enough that collisions are unimportant, but that conditions suitable for stellar collisions are much more easily satisfied in large-n clusters, i.e. n ~ 30,000 (we argue, however, that the density of the Arches cluster is insufficient for us to expect stellar collisions to have occurred in the cluster's prior evolution). We find that the character of the collision process is not such that it is a route toward smoothly filling the top end of the mass spectrum. Instead, runaway growth of one or two extreme objects can occur within less than 1 Myr after accretion is shut off, resulting in a few objects with masses several times the maximum reached by accretion. The rapid formation of these objects is due to not just the post-formation dynamical evolution of the clusters, but an interplay of dynamics and the accretional growth of the stars. We find that accretion-driven cluster shrinkage results in a distribution of gas and stars that offsets the disruptive effect of gas expulsion, and we propose that the process can lead to massive binaries and early mass segregation in star clusters.

We attempt to evaluate whether the integrated galactic IMF (IGIMF) is expected to be steeper than the IMF within individual clusters through direct evaluation of whether there is a systematic dependence of maximum stellar mass on cluster mass. We show that the result is sensitive to observational selection biases and requires an accurate knowledge of cluster ages, particularly in more populous clusters. At face value there is no compelling evidence for non-random selection of stellar masses in low mass clusters but there is arguably some evidence that the maximum stellar mass is anomalously low (compared with the expectations of random mass selection) in clusters containing more than several thousand stars. Whether or not this effect is then imprinted on the IGIMF then depends on the slope of the cluster mass function. We argue that a more economical approach to the problem would instead involve direct analysis of the upper IMF in clusters using statistical tests for truncation of the mass function. When such an approach is applied to data from hydrodynamic simulations we find evidence for truncated mass functions even in the case of simulations without feedback.

We have undertaken SED modeling of discs around low mass T Tauri stars in order to delineate regions of the infrared two colour plane (K - [8] versus K - [24]) that correspond to discs in different evolutionary stages. This provides a ready tool for classifying the nature of star-disc systems based on infrared photometry. In particular we demonstrate the distinct loci followed by discs that undergo `uniform draining' (reduction in surface density by a spatially uniform factor) from those that clear from the inside out. We draw attention to the absence of objects on this `draining locus' in those star forming regions where the 24um sensitivity would permit their detection, as compared with the ~20 objects in these regions with colours suggestive of inner holes. We thus conclude that discs predominantly clear from the inside out. We also apply our classification of the infrared two colour plane to stars of spectral type M3-M5 in the IC 348 cluster and show that some of its members (dubbed `weak excess' sources by Muzerolle et al 2010) that had previously been assumed to be in a state of clearing are instead likely to be optically thick discs in which the dust is well settled towards the mid-plane. Nevertheless, there are many discs in a state of partial clearing in IC 348, with their abundance (relative to the total population of disc bearing stars) being around four times higher than for comparable stars in Taurus. However, the number of partially cleared discs relative to the total number of late type stars is similarly low in both regions (10 and 20 % respectively). We therefore conclude that IC 348 represents a more evolved version of the Taurus population (with more of its discs being highly settled or partially cleared) but that the timescale for clearing is similarly short (a few times 10^5 years) in both cases.

We explore the behaviour of accreting protoclusters with a Monte Carlo dynamical code in order to evaluate the relative roles of accretion, two body relaxation and stellar collisions in the cluster evolution. We corroborate the suggestion of Clarke & Bonnell that the number of stellar collisions should scale as $N^{5/3} \dot M^{2/3}$ (independent of other cluster parameters, where N is the number of stars in the cluster and $\dot M$ the rate of mass accretion) and thus strengthen the argument that stellar collisions are more likely in populous (large N) clusters. We however find that the estimates of Clarke & Bonnell were pessimistic in the sense that we find that more than 99 % of the stellar collisions occur within the post-adiabatic regime as the cluster evolves towards core collapse, driven by a combination of accretion and two-body relaxation. We discuss how the inclusion of binaries may reduce the number of collisions through the reversal of core collapse but also note that it opens up another collisional channel involving the merger of stars within hard binaries; future Nbody simulations are however required in order to explore this issue.

We undertake a systematic analysis of the early (< 0.5 Myr) evolution of clustering and the stellar initial mass function in turbulent fragmentation simulations. These large scale simulations for the first time offer the opportunity for a statistical analysis of IMF variations and correlations between stellar properties and cluster richness. The typical evolutionary scenario involves star formation in small-n clusters which then progressively merge; the first stars to form are seeds of massive stars and achieve a headstart in mass acquisition. These massive seeds end up in the cores of clusters and a large fraction of new stars of lower mass is formed in the outer parts of the clusters. The resulting clusters are therefore mass segregated at an age of 0.5 Myr, although the signature of mass segregation is weakened during mergers. We find that the resulting IMF has a smaller exponent (alpha=1.8-2.2) than the Salpeter value (alpha=2.35). The IMFs in subclusters are truncated at masses only somewhat larger than the most massive stars (which depends on the richness of the cluster) and an universal upper mass limit of 150 Msun is ruled out. We also find that the simulations show signs of the IGIMF effect proposed by Weidner & Kroupa, where the frequency of massive stars is suppressed in the integrated IMF compared to the IMF in individual clusters. We identify clusters through the use of a minimum spanning tree algorithm which allows easy comparison between observational survey data and the predictions of turbulent fragmentation models. In particular we present quantitative predictions regarding properties such as cluster morphology, degree of mass segregation, upper slope of the IMF and the relation between cluster richness and maximum stellar mass. [abridged]

We review progress in numerical simulations of star cluster formation. These simulations involve the bottom-up assembly of clusters through hierarchical mergers, which produces a fractal stellar distribution at young (~0.5 Myr) ages. The resulting clusters are predicted to be mildly aspherical and highly mass-segregated, except in the immediate aftermath of mergers. The upper initial mass function within individual clusters is generally somewhat flatter than for the aggregate population. Recent work has begun to clarify the factors that control the mean stellar mass in a star-forming cloud and also the efficiency of star formation. The former is sensitive to the thermal properties of the gas while the latter depends both on the magnetic field and the initial degree of gravitational boundedness of the natal cloud. Unmagnetized clouds that are initially bound undergo rapid collapse, which is difficult to reverse by ionization feedback or stellar winds.

The formation of planets within a disc must operate within the time frame of disc dispersal, it is thus crucial to establish what is the dominant process that disperses the gaseous component of discs around young stars. Planet formation itself as well as photoevaporation by energetic radiation from the central young stellar object have been proposed as plausible dispersal mechanisms. [abridged]. In this paper we use the different metallicity dependance of X-ray photoevaporation and planet formation to discriminate between these two processes. We study the effects of metallicity, Z, on the dispersal timescale, t_phot, in the context of a photoevaporation model, by means of detailed thermal calculations of a disc in hydrostatic equilibrium irradiated by EUV and X-ray radiation from the central source. Our models show t_phot \propto Z^0.52 for a pure photoevaporation model. By means of analytical estimates we derive instead a much stronger negative power dependance on metallicity of the disc lifetime for a dispersal model based on planet formation. A census of disc fractions in lower metallicity regions should therefore be able to distinguish between the two models. A recent study by Yasui et al. in low metallicity clusters of the extreme outer Galaxy ([O/H] ~- 0.7dex and dust to gas ratio of ~0.001) provides preliminary observational evidence for shorter disc lifetimes at lower metallicities, in agreement with the predictions of a pure photoevaporation model. [abridged] We finally develop an analytical framework to study the effects of metallicity dependent photoevaporation on the formation of gas giants in the core accretion scenario. We show that accounting for this effect strengthens the conclusion that planet formation is favoured at higher metallicity. [abridged]

In this paper we consider the effects of opacity regimes on the stability of self-gravitating protoplanetary discs to fragmentation into bound objects. Using a self-consistent 1-D viscous disc model, we show that the ratio of local cooling to dynamical timescales Omega*tcool has a strong dependence on the local temperature. We investigate the effects of temperature-dependent cooling functions on the disc's gravitational stability through controlled numerical experiments using an SPH code. We find that such cooling functions raise the susceptibility of discs to fragmentation through the influence of temperature perturbations - the average value of Omega*tcool has to increase to prevent local variability leading to collapse. We find the effects of temperature dependence to be most significant in the "opacity gap" associated with dust sublimation, where the average value of Omega*tcool at fragmentation is increased by over an order of magnitude. We then use this result to predict where protoplanetary discs will fragment into bound objects, in terms of radius and accretion rate. We find that without temperature dependence, for radii < ~10AU a very large accretion rate ~10^-3 Msun/yr is required for fragmentation, but that this is reduced to 10^-4 Msun/yr with temperature-dependent cooling. We also find that the stability of discs with accretion rates < ~10^-7 Msun/yr at radii > ~50AU is enhanced by a lower background temperature if the disc becomes optically thin.

(Abridged) We present the first radiation-hydrodynamic model of a protoplanetary disc irradiated with an X-EUV spectrum. In a model where the total ionizing luminosity is divided equally between X-ray and EUV luminosity, we find a photoevaporation rate of 1.4e-8 M_sun/yr, which is two orders of magnitude greater than the case of EUV photoevaporation alone. Thus it is clear that the X-rays are the dominant driving mechanism for photoevaporation. This can be understood inasmuch as X-rays are capable of penetrating much larger columns (~1e22 cm^-2) and can thus effect heating in denser regions and at larger radius than the EUV can. The radial extent of the launching region of the X-ray heated wind is 1-70AU compared with the pure EUV case where the launch region is concentrated around a few AU. When we couple our wind mass-loss rates with models for the disc's viscous evolution, we find that, as in the pure EUV case, there is a photoevaporative switch, such that an inner hole develops at ~ 1 AU at the point that the accretion rate in the disc drops below the wind mass loss rate. At this point, the remaining disc material is quickly removed in the final 15-20% of the disc's lifetime. This is consistent with the 1e5 yr transitional timescale estimated from observations of T-Tauri stars. We also caution that although our mass-loss rates are high compared to some accretion rates observed in young stars, our model has a rather large X-ray luminosity of 2e30 erg/s; further modeling is required in order to investigate the evolutionary implications of the large observed spread of X-ray luminosities in T-Tauri stars.

In this Letter we show that if planetesimals form in spiral features in self-gravitating discs, as previously suggested by the idealised simulations of Rice et al, then in realistic protostellar discs, this process will be restricted to the outer regions of the disc (i.e. at radii in excess of several tens of A.U.). This restriction relates to the requirement that dust has to be concentrated in spiral features on a timescale that is less than the (roughly dynamical) lifetime of such features, and that such rapid accumulation requires spiral features whose fractional amplitude is not much less than unity. This in turn requires that the cooling timescale of the gas is relatively short, which restricts the process to the outer disc. We point out that the efficient conversion of a large fraction of the primordial dust in the disc into planetesimals could rescue this material from the well known problem of rapid inward migration at a $\sim$ metre size scale and that in principle the collisional evolution of these objects could help to re-supply small dust to the protostellar disc. We also point out the possible implications of this scenario for the location of planetesimal belts inferred in debris discs around main sequence stars, but stress that further dynamical studies are required in order to establish whether the disc retains a memory of the initial site of planetesimal creation.

We present new models for the X-ray photoevaporation of circumstellar discs which suggest that the resulting mass loss (occurring mainly over the radial range 10-40 AU) may be the dominant dispersal mechanism for gas around low mass pre-main sequence stars, contrary to the conclusions of previous workers. Our models combine use of the MOCASSIN Monte Carlo radiative transfer code and a self-consistent solution of the hydrostatic structure of the irradiated disc. We estimate the resulting photoevaporation rates assuming sonic outflow at the surface where the gas temperature equals the local escape temperature and derive mass loss rates of ~10^{-9} M_sun/yr, typically a factor 2-10 times lower than the corresponding rates in our previous work (Ercolano et al., 2008) where we did not adjust the density structure of the irradiated disc. The somewhat lower rates, and the fact that mass loss is concentrated towards slightly smaller radii, result from the puffing up of the heated disc at a few AU which partially screens the disc at tens of AU. (.....) We highlight the fact that X-ray photoevaporation has two generic advantages for disc dispersal compared with photoevaporation by extreme ultraviolet (EUV) photons that are only modestly beyond the Lyman limit: the demonstrably large X-ray fluxes of young stars even after they have lost their discs and the fact that X-rays are effective at penetrating much larger columns of material close to the star (abridged).

We present analytic models for the local structure of self-regulated self-gravit ating accretion discs that are subject to realistic cooling. Such an approach can be used to predict the secular evolution of self-gravitating discs (which can usefully be compared with future radiation hydrodynamical simulations) and to define various physical regimes as a function of radius and equivalent steady state accretion rate. We show that fragmentation is inevitable, given realistic rates of infall into the disc, once the disc extends to radii $> 70$ A.U. (in the case of a solar mass central object). Owing to the outward redistribution of disc material by gravitational torques, we also predict fragmentation at $> 70$ A.U. even in the case of low angular momentum cores which initially collapse to a much smaller radius. We point out that 70 A.U. is close to the median binary separation and propose that such delayed fragmentation, at the point that the disc expands to $> 70$ A.U., ensures the creation of low mass ratio companions that can avoid substantial further growth and consequent evolution towards unit mass ratio. We thus propose this as a promising mechanism for producing low mass ratio binaries, which, while abundant observationally, are severely underproduced in hydrodynamical models.

We re-examine the recent suggestion of a high fraction of transition discs (i.e. those with a cleared inner hole) in M stars, motivated by the fact that we expect that, for M stars, even discs without inner holes should exhibit very weak excess shortward of around 10um. Our analysis of spectral energy distribution models suggest that this indeed means that M stars where a detectable excess begins at around 6um may be mis-classified as transition discs when in fact they have optically thick dust extending in to the dust sublimation radius. Consequently, we estimate that the transition disc fraction among M stars in the Coronet cluster is ~15 +/-10 % (rather than the recently claimed value of 50%). This revised figure would imply that the transition disc fraction is not after all markedly higher in later type stars. We suggest that for M stars, transition discs can only be readily identified if they have emission that is close to photospheric out to > 10um.

We investigate the potential of using ratios of fine structure and near-infrared forbidden line transitions of atomic carbon to diagnose protoplanetary disk extension. Using results from 2D photoionisation and radiative transfer modeling of a realistic protoplanetary disk structure irradiated by X-rays from a T Tauri star, we obtain theoretical emission maps from which we construct radial distributions of the strongest emission lines produced in the disk. We show that ratios of fine structure to near-infrared forbidden line emission of atomic carbon are especially promising to constrain the minimum size of gaseous protoplanetary disks. While theoretically viable, the method presents a number of observational difficulties that are also discussed here.

We perform numerical analyses of the structure induced by gravitational instabilities in cooling gaseous accretion discs. For low enough cooling rates a quasi-steady configuration is reached, with the instability saturating at a finite amplitude in a marginally stable disc. We find that the saturation amplitude scales with the inverse square root of the cooling parameter beta = t_cool / t_dyn, which indicates that the heating rate induced by the instability is proportional to the energy density of the induced density waves. We find that at saturation the energy dissipated per dynamical time by weak shocks due is of the order of 20 per cent of the wave energy. From Fourier analysis of the disc structure we find that while the azimuthal wavenumber is roughly constant with radius, the mean radial wavenumber increases with radius, with the dominant mode corresponding to the locally most unstable wavelength. We demonstrate that the density waves excited in relatively low mass discs are always close to co-rotation, deviating from it by approximately 10 per cent. This can be understood in terms of the flow Doppler-shifted phase Mach number -- the pattern speed self-adjusts so that the flow into spiral arms is always sonic. This has profound effects on the degree to which transport through self-gravity can be modelled as a viscous process. Our results thus provide (a) a detailed description of how the self-regulation mechanism is established for low cooling rates, (b) a clarification of the conditions required for describing the transport induced by self-gravity through an effective viscosity, (c) an estimate of the maximum amplitude of the density perturbation before fragmentation occurs, and (d) a simple recipe to estimate the density perturbation in different thermal regimes.

We have made a new compilation of observations of maximum stellar mass versus cluster membership number from the literature, which we analyse for consistency with the predictions of a simple random drawing hypothesis for stellar mass selection in clusters. Previously, Weidner and Kroupa have suggested that the maximum stellar mass is lower, in low mass clusters, than would be expected on the basis of random drawing, and have pointed out that this could have important implications for steepening the integrated initial mass function of the Galaxy (the IGIMF) at high masses. Our compilation demonstrates how the observed distribution in the plane of maximum stellar mass versus membership number is affected by the method of target selection; in particular, rather low n clusters with large maximum stellar masses are abundant in observational datasets that specifically seek clusters in the environs of high mass stars. Although we do not consider our compilation to be either complete or unbiased, we discuss the method by which such data should be statistically analysed. Our very provisional conclusion is that the data is not indicating any striking deviation from the expectations of random drawing.

We present MOCASSIN 2D photoionisation and dust radiative transfer models of a prototypical T Tauri disk irradiated by X-rays from the young pre-main sequence star. The calculations demonstrate a layer of hot gas reaching temperatures of ~10^6 K at small radii and ~10^4 K at a distance of 1 AU. The gas temperatures decrease sharply with depth, but appear to be completely decoupled from dust temperatures down to a column depth of ~5*10^21 cm^-2. We predict that several fine-structure and forbidden lines of heavy elements, as well as recombination lines of hydrogen and helium, should be observable with current and future instrumentation, although optical lines may be smothered by the stellar spectrum. Predicted line luminosities are given for the the brightest collisionally excited lines (down to ~10^-8L_sun, and for recombination transitions from several levels of HI and HeI. The mass loss rate due to X-ray photoevaporation estimated from our models is of the order of 10^-8 M_sun yr^-1, implying a dispersal timescale of a few Myr for a disk of mass 0.027 M_sun, which is the mass of the disk structure model we employed. We discuss the limitations of our model and highlight the need for further calculations that should include the simultaneous solution of the 2D radiative transfer problem and the 1D hydrostatic equilibrium in the polar direction.

We use SPITZER IRAC 3.6 and 4.5micron near infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), optical B, V and I and 2MASS Ks band data to produce mass surface density maps of M81. The IRAC 3.6 and 4.5micron data, whilst dominated by emission from old stellar populations, is corrected for small-scale contamination by young stars and PAH emission. The I band data are used to produce a mass surface density map by a B-V colour-correction, following the method of Bell and de Jong. We fit a bulge and exponential disc to each mass map, and subtract these components to reveal the non-axisymmetric mass surface density. From the residual mass maps we are able to extract the amplitude and phase of the density wave, using azimuthal profiles. The response of the gas is observed via dust emission in the 8micron IRAC band, allowing a comparison between the phase of the stellar density wave and gas shock. The relationship between this angular offset and radius suggests that the spiral structure is reasonably long lived and allows the position of corotation to be determined.

We investigate the evolution of accretion luminosity $L_{\rm acc}$ and stellar luminosity ${L_\ast}$ in pre-mainsequence stars. We make the assumption that when the star appears as a Class II object, the major phase of accretion is long past, and the accretion disc has entered its asymptotic phase. We use an approximate stellar evolution scheme for accreting pre-mainsequence stars based on Hartmann, Cassen & Kenyon, 1997. We show that the observed range of values $k = L_{\rm acc}/L_\ast$ between 0.01 and 1 can be reproduced if the values of the disc mass fraction $M_{\rm disc}/M_*$ at the start of the T Tauri phase lie in the range 0.01 -- 0.2, independent of stellar mass. We also show that the observed upper bound of $L_{\rm acc} \sim L_\ast$ is a generic feature of such disc accretion. We conclude that as long as the data uniformly fills the region between this upper bound and observational detection thresholds, then the degeneracies between age, mass and accretion history severely limit the use of this data for constraining possible scalings between disc properties and stellar mass.

We investigate the orbital evolution of planetesimals in a self-gravitating circumstellar disc in the size regime ($\sim 1-5000$ km) where the planetesimals behave approximately as test particles in the disc's non-axisymmetric potential. We find that the particles respond to the stochastic, regenerative spiral features in the disc by executing large random excursions (up to a factor of two in radius in $\sim 1000$ years), although typical random orbital velocities are of order one tenth of the Keplerian speed. The limited time frame and small number of planetesimals modeled does not permit us to discern any {\it net} direction of planetesimal migration. Our chief conclusion is that the high eccentricities ($\sim 0.1$) induced by interaction with spiral features in the disc is likely to be highly unfavourable to the collisional growth of planetesimals in this size range while the disc is in the self-gravitating regi me. Thus {\it if}, as recently argued by Rice et al 2004, 2006, the production of planetesimals gets under way when the disc is in the self-gravitating regime (either at smaller planetesimal size scales, where gas drag is important, or via gravitational fragmentation of the solid component), then the planetesimals thus produced would not be able to grow collisionally until the disc ceased to be self-gravitating. It is unclear, however,given the large amplitude excursions undergone by planetesimals in the self-gravitating disc, whether they would be retained in the disc throughout this period, or whether they would instead be lost to the central star.

A number of previous studies of the fragmentation of self-gravitating protostellar discs have modeled radiative cooling with a cooling timescale (t_{cool}) parameterised as a simple multiple (beta_{cool}) of the local dynamical timescale. Such studies have delineated the `fragmentation boundary' in terms of a critical value of beta_{cool} (beta_{crit}), where the disc fragments if beta_{cool} < beta_{crit}. Such an approach however begs the question of how in reality a disc could ever be assembled with beta_{cool} < beta_{crit}. Here we adopt the more realistic approach of gradually reducing beta_{cool}, as might correspond to changes in thermal regime due to secular changes in the disc density profile. We find that when beta_{cool} is gradually reduced (on a timescale longer than t_{cool}), the disc is stabilised against fragmentation, compared with models in which beta_{cool} is reduced rapidly. We therefore conclude that a disc's ability to remain in a self-regulated, self-gravitating state (without fragmentation) is partly dependent on its thermal history, as well as its current cooling rate. Nevertheless, a slow reduction in t_{cool} appears only to lower the fragmentation boundary by about a factor two in t_{cool} and thus only permits maximum alpha values (parameterising the efficiency of angular momentum transfer in the disc) that are about a factor two higher than determined hitherto. Our results therefore do not undermine the notion of a fundamental upper limit to the heating rate that can be delivered by gravitational instabilities before the disc is subject to fragmentation. An important implication of this work, therefore, is that self-gravitating discs can enter into the regime of fragmentation via secular evolution and it is not necessary to invoke rapid (impulsive) events to trigger fragmentation.

The probability of a star hosting a planet that is detectable in radial velocity surveys increases Ppl(Z) oc 10^2Z, where Z is metallicity. Core accretion models reproduce this trend, since the protoplanetary disk of a high metallicity star has a high density of solids and so forms cores which accrete gas before the primordial gas disk dissipates. This paper considers the origin of the form of Ppl(Z). We introduce a simple model in which detectable planets form when the mass of solids in the protoplanetary disk, Ms, exceeds a critical value. In this model the form of Ppl(Z) is a direct reflection of the distribution of protoplanetary disk masses, Mg, and the observed Ppl(Z) is reproduced if P(Mg>Mg') oc 1/Mg'^2. We argue that a protoplanetary disk's sub-mm dust mass is a pristine indicator of the mass available for planet-building and find the observed sub-mm disk mass distribution is consistent with the observed Ppl(Z) if Ms>0.5M_J is required to form detectable planets. Any planet formation model which imposes a critical solid mass for planet formation would reproduce the observed Ppl(Z), and core accretion models are empirically consistent with a threshold criterion. We identify 7 protoplanetary disks which, by rigid application of this criterion, would be expected to form detectable planets. A testable prediction is that Ppl(Z) should flatten both for Z>0.5dex and as more distant and lower mass planets are discovered. Further, combining this model with one in which the evolution of a star's debris disk is also influenced by the solid mass in its protoplanetary disk, results in the prediction that debris disks detected around stars with planets should be more infrared luminous than those around stars without planets in tentative agreement with recent observations.

Eta Chamaeleontis is a unique young (~9 Myr) association with 18 systems concentrated in a radius of ~35 arcmin, i.e. 1pc at the cluster distance of 97pc. No other members have been found up to 1.5 degrees from the cluster centre. The cluster mass function is consistent with the IMF of other rich young open clusters in the higher mass range but shows a clear deficit of low mass stars and brown dwarfs with no objects below 0.1Msun. The aim of this paper is to test whether this peculiar mass function could result from dynamical evolution despite the young age of the cluster. We performed N-body numerical calculations starting with a log-normal IMF and different initial conditions in terms of number of systems and cluster radius using the code NBODY3. We simulated the cluster dynamical evolution over 10 Myr and compared the results to the observations. We found that it is possible to reproduce eta Cha when starting with a very compact configuration (with Ninit=40 and R0=0.005pc) which suggests that the IMF of the association might not be abnormal. The high initial density might also explain the deficit of wide binaries that is observed in the cluster.

Eduardo Delgado was due to have presented a poster at this meeting on his latest results on the formation of extreme mass ratio binaries. Tragically, Eduardo was among those killed in a hiking accident in Tenerife earlier this year. As his PhD supervisor, and as a longstanding collaborator, the organisers of this meeting kindly invited me to incorporate a report on his most recent work into a more general tribute to his life and work. I will reflect on Eduardo's scientific career, the problems that motivated him and his achievements, focusing particularly on a problem which had intrigued us both for several years and on which Eduardo was making important progress at the time of his death. Finally, I will mention the personal qualities that Eduardo brought to his work and the acute sense of loss that is shared by all those - friends and collaborators - who were privileged to know him.

We present a new fast algorithm which allows the simulation of ionising radiation emitted from point sources to be included in high-resolution three-dimensional smoothed particle hydrodynamics simulations of star cluster formation. We employ a Str\"omgren volume technique in which we use the densities of particles near the line-of-sight between the source and a given target particle to locate the ionisation front in the direction of the target. Along with one--dimensional tests, we present fully three--dimensional comparisons of our code with the three--dimensional Monte-Carlo radiative transfer code, MOCASSIN, and show that we achieve good agreement, even in the case of highly complex density fields.

We examine self-consistent parameterizations of the high-mass stellar population and resulting feedback, including mechanical, radiative, and chemical feedback, as we understand them locally. To date, it appears that the massive star population follows simple power-law clustering that extends down to individual field OB stars, and the robust stellar IMF seems to have a constant upper-mass limit. These properties result in specific patterns in the HII region LF and ionization of the diffuse, warm ionized medium. The resulting SNe generate superbubbles whose size distribution is also described by a simple power law, and from which a galaxy's porosity parameter is easily derived. A critical star-formation threshold can then be estimated, above which the escape of Lyman continuum photons, hot gas, and nucleosynthetic products is predicted. A first comparison with a large H-alpha sample of galaxies is broadly consistent with this prediction, and suggests that ionizing photons likely escapes from starburst galaxies. The superbubble size distribution also offers a basis for a Simple Inhomogeneous Model for galactic chemical evolution, which is especially relevant to metal-poor systems and instantaneous metallicity distributions. This model offers an alternative interpretation of the Galactic halo metallicity distribution and emphasizes the relative importance of star-formation intensity, in addition to age, in a system's evolution. The fraction of zero-metallicity, Population III stars is easily predicted for any such model. We emphasize that all these phenomena can be modeled in a simple, analytic framework over an extreme range in scale, offering powerful tools for understanding the role of massive stars in the cosmos. (Abridged)

We report on the results of an i-band time-series photometric survey of NGC 2516 using the CTIO 4m Blanco telescope and 8k Mosaic-II detector, achieving better than 1% photometric precision per data point over 15<~i<~19. Candidate cluster members were selected from a V vs V-I colour magnitude diagram over 16<V<26 (covering masses from 0.7 M_solar down to below the brown dwarf limit), finding 1685 candidates, of which we expect ~1000 to be real cluster members, taking into account contamination from the field (which is most severe at the extremes of our mass range). Searching for periodic variations in these gave 362 detections over the mass range 0.15<~M/M_solar<~0.7. The rotation period distributions were found to show a remarkable morphology as a function of mass, with the fastest rotators bounded by P>0.25 days, and the slowest rotators for M<~0.5 M_solar bounded by a line of P \propto M^3, with those for M>~0.5 M_solar following a flatter relation closer to P ~ constant. Models of the rotational evolution were investigated, finding that the evolution of the fastest rotators was well-reproduced by a conventional solid body model with a mass-dependent saturation velocity, whereas core-envelope decoupling was needed to reproduce the evolution of the slowest rotators. None of our models were able to simultaneously reproduce the behaviour of both populations.

We present models in which the photoevaporation of discs around young stars by an external ultraviolet source (as computed by Adams et al 2004) is coupled with the internal viscous evolution of the discs. These models are applied to the case of the Orion Nebula Cluster, where the presence of a strong ultraviolet field from the central OB stars, together with a detailed census of circumstellar discs and photoevaporative flows, is well established. In particular we investigate the constraints that are placed on the initial disc properties in the ONC by the twin requirement that most stars possess a disc on a scale of a few A.U., but that only a minority ($< 20 %$) are resolved by HST at a scale of 50 A.U.. We find that these requirements place very weak constraints on the initial radius distribution of circumstellar discs: the resulting size distribution readily forgets the initial radius distribution, owing to the strong positive dependence of the photoevaporation rate on disc radius. Instead, the scarcity of large discs reflects the relative scarcity of initially massive discs (with mass $> 0.1 M_\odot$). The ubiquity of discs on a small scale, on the other hand, mainly constrains the timespan over which the discs have been exposed to the ultraviolet field ($< 2 $Myr). We argue that the discs that are resolved by HST represent a population of discs in which self-gravity was important at the time that the dominant central OB star switched on, but that, according to our models, self-gravity is unlikely to be important in these discs at the present time. We discuss the implications of our results for the so-called proplyd lifetime problem.

In this paper, we use high-resolution smoothed particle hydrodynamics (SPH) simulations to investigate the response of a marginally stable self-gravitating protostellar disc to a close parabolic encounter with a companion discless star. Our main aim is to test whether close brown dwarfs or massive planets can form out of the fragmentation of such discs. We follow the thermal evolution of the disc by including the effects of heating due to compression and shocks and a simple prescription for cooling and find results that contrast with previous isothermal simulations. In the present case we find that fragmentation is inhibited by the interaction, due to the strong effect of tidal heating, which results in a strong stabilization of the disc. A similar behaviour was also previously observed in other simulations involving discs in binary systems. As in the case of isolated discs, it appears that the condition for fragmentation ultimately depends on the cooling rate.

We study the response of the gaseous component of a galactic disc to the time dependent potential generated by N-body simulations of a spiral galaxy. The results show significant variation of the spiral structure of the gas which might be expected to result in significant fluctuations in the Star Formation Rate (SFR). Pronounced local variations of the SFR are anticipated in all cases. Bursty histories for the global SFR, however, require that the mean surface density is much less (around an order of magnitude less) than the putative threshold for star formation. We thus suggest that bursty star formation histories, normally attributed to mergers and/or tidal interactions, may be a normal pattern for gas poor isolated spiral galaxies.

We report on the results of a V and i-band time-series photometric survey of M34 (NGC 1039) using the Wide Field Camera (WFC) on the Isaac Newton Telescope (INT), achieving better than 1% precision per data point for 13 <~ i <~ 17. Candidate cluster members were selected from a V vs V-I colour-magnitude diagram over 14 < V < 24 (0.12 <~ M/Msun <~ 1.0), finding 714 candidates, of which we expect ~ 400 to be real cluster members (taking into account contamination from the field). The mass function was computed, and found to be consistent with a log-normal distribution in dN/dlogM. Searching for periodic variable objects in the candidate members gave 105 detections over the mass range 0.25 < M/Msun < 1.0. The distribution of rotation periods for 0.4 < M/Msun < 1.0 was found to peak at ~ 7 days, with a tail of fast rotators down to periods of ~ 0.8 days. For 0.25 < M/Msun < 0.4 we found a peak at short periods, with a lack of slow rotators (eg. P >~ 5 days, consistent with the work of other authors (eg. Scholz & Eisloffel 2004) at very low masses. Our results are interpreted in the context of previous work, finding that we reproduce the same general features in the rotational period distributions. A number of rapid rotators were found with velocities ~ a factor of two lower than in the Pleiades, consistent with models of angular momentum evolution assuming solid body rotation without needing to invoke core-envelope decoupling.

We examine the recent data and analysis of Natta et al. concerning the accretion rate on to young stars as a function of stellar mass, and conclude that the apparently steep dependence of accretion rate on mass is strongly driven by selection/detection thresholds. We argue that a convincing demonstration of a physical relationship between accretion and stellar mass requires further studies which, as is the case for Natta et al., include information on upper limits, and which quantify the possible incompleteness of the sample, at both low and high accretion rates. We point out that the distribution of detections in the ($M, \dot{M}$)-plane can in principle be used to test conventional accretion disc evolutionary models, and that higher sensitivity observations might be able to test the hypothesis of accelerated disc clearing at late times.

The formation of a binary system surrounded by disks is the most common outcome of stellar formation. Hence studying and understanding the formation and the evolution of binary systems and associated disks is a cornerstone of star formation science. Moreover, since the components within binary systems are coeval and the sizes of their disks are fixed by the tidal truncation of their companion, binary systems provide an ideal "laboratory" in which to study disk evolution under well defined boundary conditions. In this paper, we review observations of several inner disk diagnostics in multiple systems, including hydrogen emission lines (indicative of ongoing accretion), $K-L$ and $K-N$ color excesses (evidence of warm inner disks), and polarization (indicative of the relative orientations of the disks around each component). We examine to what degree these properties are correlated within binary systems and how this degree of correlation depends on parameters such as separation and binary mass ratio. These findings will be interpreted both in terms of models that treat each disk as an isolated reservoir and those in which the disks are subject to re-supply from some form of circumbinary reservoir, the observational evidence for which we will also critically review. The planet forming potential of multiple star systems is discussed in terms of the relative lifetimes of disks around single stars, binary primaries and binary secondaries. Finally, we summarize several potentially revealing observational problems and future projects that could provide further insight into disk evolution in the coming decade

We use numerical simulations of the fragmentation of a 1000 solar mass molecular cloud and the formation of a stellar cluster to study how the initial conditions for star formation affect the resulting initial mass function (IMF). In particular, we are interested in the relation between the thermal Jeans mass in a cloud and the knee of the initial mass function, i.e. the mass separating the region with a flat IMF slope from that typified by a steeper, Salpeter-like, slope. In three isothermal simulations with M_J=1 solar mass, M_J=2 solar masses and M_J=5 solar masses, the number of stars formed, at comparable dynamical times, scales roughly with the number of initial Jeans masses in the cloud. The mean stellar mass also increases (though less than linearly) with the initial Jeans mass in the cloud. It is found that the IMF in each case displays a prominent knee, located roughly at the mass scale of the initial Jeans mass. Thus clouds with higher initial Jeans masses produce IMFs which are shallow to higher masses. This implies that a universal IMF requires a physical mechanism that sets the Jeans mass to be near 1 solar mass. Simulations including a barotropic equation of state as suggested by Larson, with cooling at low densities followed by gentle heating at higher densities, are able to produce realistic IMFs with the knee located at approximately 1 solar mass, even with an initial M_J=5 solar masses. We therefore suggest that the observed universality of the IMF in the local Universe does not require any fine tuning of the initial conditions in star forming clouds but is instead imprinted by details of the cooling physics of the collapsing gas.

In this paper we consider the effect of the direct ionizing stellar radiation field on the evolution of protoplanetary discs subject to photoevaporative winds. We suggest that models which combine viscous evolution with photoevaporation of the disc (e.g. Clarke, Gendrin & Sotomayor 2001) incorrectly neglect the direct field after the inner disc has drained, at late times in the evolution. We construct models of the photoevaporative wind produced by the direct field, first using simple analytic arguments and later using detailed numerical hydrodynamics. We find that the wind produced by the direct field at late times is much larger than has previously been assumed, and we show that the mass-loss rate scales as $R_{in}^{1/2}$ (where $R_{in}$ is the radius of the instantaneous inner disc edge). We suggest that this result has important consequences for theories of disc evolution, and go on to consider the effects of this result on disc evolution in detail in a companion paper (Alexander, Clarke & Pringle 2006b).

We present a new model for protoplanetary disc evolution. This model combines viscous evolution with photoevaporation of the disc, in a manner similar to Clarke, Gendrin & Sotomayor (2001). However in a companion paper (Alexander, Clarke & Pringle 2006a) we have shown that at late times such models must consider the effect of stellar radiation directly incident on the inner disc edge, and here we model the observational implications of this process. We find that the entire disc is dispersed on a time-scale of order $10^5$yr after a disc lifetime of a few Myr, consistent with observations of T Tauri (TT) stars. We use a simple prescription to model the spectral energy distribution of the evolving disc, and demonstrate that the model is consistent with observational data across a wide range of wavelengths. We note also that the model predicts a short ``inner hole'' phase in the evolution of all TT discs, and make predictions for future observations at mid-infrared and millimetre wavelengths.

In this paper we discuss the nature and the possible formation scenarios of the companion of the brown dwarf 2MASS 1207334-393254. We initially discuss the basic physical properties of this object and conclude that, although from its absolute mass ($5M_{\rm Jup}$), it is a planetary object, in terms of its mass ratio $q$ and of its separation $a$ with respect to the primary brown dwarf, it is consistent with the statistical properties of binaries with higher primary mass. We then explore the possible formation mechanism for this object. We show that the standard planet formation mechanism of core accretion is far too slow to form this object within 10 Myr, the observed age of the system. On the other hand, the alternative mechanism of gravitational instability (proposed both in the context of planet and of binary formation) may, in principle, work and form a system with the observed properties.

We present the results of a photometric monitoring campaign of three well studied FU Orionis systems (FU Orionis, V1057 Cygni and V1515 Cygni) undertaken at Maidanak Observatory between 1981 and 2003. When combined with photometric data in the literature, this database provides a valuable resource for searching for short timescale variability - both periodic and aperiodic - as well as for studying the secular evolution of these systems. In the case of V1057 Cyg (which is the system exhibiting the largest changes in brightness since it went into outburst) we compare the photometric data with time dependent models. We show that prior to the end of the `plateau' stage in 1996, the evolution of V1057 Cyg in the $V-(B-V)$ colour-magnitude diagram is well represented by disc instability models in which the outburst is triggered by some agent - such as an orbiting planet - in the inner disc. Following the end of the plateau phase in 1996, the dimming and irregular variations are consistent with occultation of the source by a variable dust screen, which has previously been interpreted in terms of dust condensation events in the observed disc wind. Here we instead suggest that this effect results from the interaction between the wind and an infalling dusty envelope, the existence of this envelope having been previously invoked in order to explain the mid infrared emissio of FU Orionis systems. We discuss how this model may explain some of the photometric and spectroscopic characteristics of FU Orionis systems in general.

We construct a protostellar disk model that takes into account the combined effect of viscous evolution, photoevaporation and the differential radial motion of dust grains and gas. For T Tauri disks, the lifetimes of dust disks that are mainly composed of millimeter sized grains are always shorter than the gas disks' lifetimes, and become similar only when the grains are fluffy (density < 0.1 g cm^{-3}). If grain growth during the classical T Tauri phase produces plenty of millimeter sized grains, such grains completely accrete onto the star in 10^7 yr, before photoevaporation begins to drain the inner gas disk and the star evolves to the weak line T Tauri phase. In the weak line phase, only dust-poor gas disks remain at large radii (> 10 AU), without strong signs of gas accretion nor of millimeter thermal emission from the dust. For Herbig AeBe stars, the strong photoevaporation clears the inner disks in 10^6 yr, before the dust grains in the outer disk migrate to the inner region. In this case, the grains left behind in the outer gas disk accumulate at the disk inner edge (at 10-100 AU from the star). The dust grains remain there even after the entire gas disk has been photoevaporated, and form a gas-poor dust ring similar to that observed around HR 4796A. Hence, depending on the strength of the stellar ionizing flux, our model predicts opposite types of products around young stars. For low mass stars with a low photoevaporation rate, dust-poor gas disks with an inner hole would form, whereas for high mass stars with a high photoevaporation rate, gas-poor dust rings would form. This prediction should be examined by observations of gas and dust around weak line T Tauri stars and evolved Herbig AeBe stars.

We present the first ever hydrodynamic calculations of star cluster formation that incorporate the effect of feedback from ionising radiation. In our simulations, the ionising source forms in the cluster core at the intersection of several dense filaments of inflowing gas. We show that these filaments collimate ionised outflows and suggest such an environmental origin for at least some observed outflows in regions of massive star formation. Our simulations show both positive feedback (i.e. promotion of star formation in neutral gas compressed by expanding HII regions) and negative feedback (i.e. suppression of the accretion flow in to the central regions). We show that the volume filling factor of ionised gas is very different in our simulations than would result from the case where the central source interacted with an azimuthally smoothed gas density distribution. As expected, gas density is the key parameter in determining whether clusters are unbound by photoionising radiation. Nevertheless, we find - on account of the acceleration of a small fraction of the gas to high velocities in the outflows - that the deposition in the gas of an energy that exceeds the binding energy of the cluster is not a sufficient criterion for unbinding the bulk of the cluster mass.

We derive the expectation value for the maximum stellar mass (m_max) in an ensemble of N stars, as a function of the IMF upper-mass cutoff (m_up) and N. We statistically demonstrate that the upper IMF of the local massive star census observed thus far in the Milky Way and Magellanic Clouds clearly exhibits a universal upper mass cutoff around 120 - 200 M_sun for a Salpeter IMF, although the result is more ambiguous for a steeper IMF.

We present the results of an analysis of ultraviolet observations of T Tauri Stars (TTS). By analysing emission measures taken from the literature we derive rates of ionizing photons from the chromospheres of 5 classical TTS in the range ~10^41-10^44 photons/s, although these values are subject to large uncertainties. We propose that the HeII/CIV line ratio can be used as a reddening-independent indicator of the hardness of the ultraviolet spectrum emitted by TTS. By studying this line ratio in a much larger sample of objects we find evidence for an ionizing flux which does not decrease, and may even increase, as TTS evolve. This implies that a significant fraction of the ionizing flux from TTS is not powered by the accretion of disc material onto the central object, and we discuss the significance of this result and its implications for models of disc evolution. The presence of a significant ionizing flux in the later stages of circumstellar disc evolution provides an important new constraint on disc photoevaporation models.

We have undertaken a series of hydrodynamic + N-body simulations in order to explore the properties of young stars. Our results suggest that the IMF may be sensitive to environment in its substellar region, with more brown dwarfs being formed where clusters are denser or more compact. We find that multiple stars are a natural outcome of collapsing turbulent flows, with a high incidence of N > 2 multiples. We find a positive correlation of multiplicity with primary mass but a companion frequency that decreases with age. Binary brown dwarfs are rarely formed, in conflict with observations. Brown dwarfs as companions are predominantly found orbiting binaries or triples at large separations.

We present the results of a theoretical study investigating the effects of photoionization and heating by X-rays on discs around low-mass stars. In particular we address the question of whether or not X-rays can drive a disc wind. First, we construct a 1-dimensional ``quasi-hydrostatic'' model, which solves for the vertical structure introduced by X-ray heating. We consider uniform X-ray illumination of the disc, but the X-ray fluxes required to heat the disc significantly are much greater than those seen by recent observations. When the model is extended to consider heating from a central X-ray source we find that the 1-dimensional model is only valid very close to the star. We extend our analysis to consider a simple 2-dimensional model, treating the disc as a two-layered structure and solving for its density profile self-consistently. For T Tauri stars we are able to set a crude upper limit on the mass-loss rate that can be driven by X-ray photoevaporation, with a value of ~10^-13 g/cm^2/s. Our model is designed to maximise this value, and most likely over-estimates it significantly. However we still find a mass-loss rate which is less than that found in studies of ultraviolet photoevaporation. We conclude that in the presence of a significant UV field, X-ray driven disc winds are unlikely to play a significant role in the evolution of discs around low-mass stars.

FU Orionis are young stellar objects undergoing episodes of enhanced luminosity, which are generally ascribed to a sudden increase of mass accretion rate in the surrounding protostellar disc. Models invoking a thermal instability in the disc are able to reproduce many features of the outburst, but cannot explain the rapid rise time-scale observed in many cases. Here we explore the possibility (originally suggested by Clarke & Syer 1996) that the thermal instability is triggered away from the disc inner edge (at a distance of $\approx 10R_{\odot}$ from the central protostar) due to the presence of a massive planet embedded in the disc. We have constructed simple, one-dimensional, time-dependent models of the disc evolution, taking into account both the interaction between the disc and the planet, and the thermal evolution of the disc. We are indeed able to reproduce rapid rise outbursts (with rise time-scale $\approx 1$ yr), with a planet mass $M_{\mathrm{s}}=10-15M_{\mathrm{Jupiter}}$. We show that the luminosity and the duration of the outbursts are increasing functions of planet mass. We also show that the inward migration of the planet is significantly slowed once it reaches the radius where it is able to trigger FU Orionis outbursts, thus suggesting that a single planet may be involved in triggering several outbursts.

We show how the viscous evolution of Keplerian accretion discs can be understood in terms of simple kinetic theory. Although standard physics texts give a simple derivation of momentum transfer in a linear shear flow using kinetic theory, many authors, as detailed by Hayashi & Matsuda 2001, have had difficulties applying the same considerations to a circular shear flow. We show here how this may be done, and note that the essential ingredients are to take proper account of, first, isotropy locally in the frame of the fluid and, second, the geometry of the mean flow.

Abridged/ We present numerical results on the properties of young binary and multiple stellar systems. Our analysis is based on a series of SPH + Nbody simulations of the fragmentation of small molecular clouds, that fully resolve the opacity limit for fragmentation. We have produced a statistically significant number of stable multiple systems, with components separations in the range 1-1000 AU. At the end of the hydrodynamical evolution (0.5 Myr) we find that ~60% of stars and brown dwarfs are members of multiples systems, with about a third of these being low mass, weakly bound outliers in wide eccentric orbits. Our results imply that in the stellar regime most stars are in multiples (~80%) and that this fraction is an increasing function of primary mass. After Nbody integration to 10.5 Myr, the percentage of bound objects has dropped to ~40%, as most very low mass stars and brown dwarfs have been released to the field. Brown dwarfs are never found to be very close companions to stars (brown dwarf desert at very small separations), but one case exists of a brown dwarf companion at intermediate separations (10 AU). Our simulations can accommodate the existence of brown dwarf companions at large separations, but only if the primaries of these systems are themselves multiples. We have compared the outcome of our simulations with the properties of real stellar systems as deduced from the IR CM diagram of the Praesepe cluster and from spectroscopic and high-resolution imaging surveys of young clusters and the field.

We present a new method for estimating the corotation radius in tightly wound spiral galaxies, through analysis of the radial variation of the offset between arms traced by the potential (P-arms) and those traced by dust (D-arms). We have verified the predictions of semi-analytical theory through hydrodynamical simulations and have examined the uniqueness of the galactic parameters that can be deduced by this method. We find that if the range of angular offsets measured at different radii in a galaxy is greater than around pi/4, it is possible to locate the radius of corotation to within ~ 25%. We argue that the relative location of the P- and D-arms provides more robust constraints on the galactic parameters than can be inferred from regions of enhanced star formation (SF-arms), since interpretation of the latter involves uncertainties due to reddening and the assumed star formation law. We thus stress the importance of K-band studies of spiral galaxies.

We apply an algorithm for radiative feedback on a dusty flow (detailed in Edgar and Clarke (2003)) to the problem of Bondi--Hoyle accretion. This calculation is potentially relevant to the formation of massive stars in ultradense cores of stellar clusters. We find that radiative feedback is \emph{more effective} than in the case of previous calculations in spherical symmetry. The Bondi-Hoyle geometry implies that material is flowing nearly tangentially when it experiences the sharp radiative impulse at the dust destruction radius, and consequently it is readily perturbed into outflowing orbits. We find that it is difficult for stellar masses to grow beyond around 10 M_sol (for standard interstellar dust abundances). We discuss the possible implications of this result for the formation mechanism of OB stars in cluster cores. We end by proposing a series of conditions which must be fulfilled if Bondi--Hoyle accretion is to continue.

We address the issue of the production of Lyman continuum photons by T Tauri stars, in an attempt to provide constraints on theoretical models of disc photoionisation. By treating the accretion shock as a hotspot on the stellar surface we show that Lyman continuum photons are produced at a rate approximately three orders of magnitude lower than that produced by a corresponding black body, and that a strong Lyman continuum is only emitted for high mass accretion rates. When our models are extended to include a column of material accreting on to the hotspot we find that the accretion column is extremely optically thick to Lyman continuum photons. Further, we find that radiative recombination of hydrogen atoms within the column is not an efficient means of producing photons with energies greater than 13.6eV, and find that an accretion column of any conceivable height suppresses the emission of Lyman continuum photons to a level below or comparable to that expected from the stellar photosphere. The photospheric Lyman continuum is itself much too weak to affect disc evolution significantly, and we find that the Lyman continuum emitted by an accretion shock is similarly unable to influence disc evolution significantly. This result has important consequences for models which use photoionisation as a mechanism to drive the dispersal of circumstellar discs, essentially proving that an additional source of Lyman continuum photons must exist if disc photoionisation is to be significant.

Abridged: We have undertaken a series of hydrodynamical simulations of multiple star formation in small turbulent molecular clouds. Our goal is to determine the sensitivity of the properties of the resulting stars and brown dwarfs to variations in the initial conditions imposed. In this paper we report on the results obtained by applying two different initial turbulent velocity fields. The slope of the turbulent power-law spectrum alpha is set to -3 in half of the calculations and to -5 in the other half. We find that, whereas the stellar mass function seems to only be weakly dependent on the value of alpha, the sub-stellar mass function turns out to be more sensitive to the initial slope of the velocity field. We argue that, since the role of turbulence is to create substructure from which gravitational instabilities may grow, variations in other initial conditions that also determine the fragmentation process are likely to affect the shape of the sub-stellar mass function as well. The absence of many planetary mass `free-floaters' in our simulations, especially in the mass range 1-10 MJ, suggests that, if these objects are abundant, they are likely to form by similar mechanisms to those thought to operate in quiescent accretion discs, instead of via instabilities in gravitationally unstable discs. Finally, we find that multiple and single stars share comparable kinematical properties, both populations being able to attain velocities in the range 1-10 km/s. From these values we draw the conclusion that only low-mass star-forming regions such as Taurus-Auriga or Ophiuchus, where the escape speed is low, might have suffered some depletion of its single and binary stellar population.

We show that if FU Orionis objects harbour hot Jupiters embedded in their discs, the resulting non-axisymmetric dissipation profile in the disc would be manifest as time-dependent distortions in the absorption line profiles of these objects. In order to affect the infrared line profiles, planets must lie within approximately 0.5 au of the central star, whereas only planets within 0.1 au would influence the optical line profiles. The timescale for modulation of the line profiles is relatively short (months) in each case, so that the effect could not have been discovered from published spectra (which combine data taken in different observing seasons). The detection of hot Jupiters in FU Orionis objects would be in line with the expectations of tidal migration theories (which predict a high incidence of close planets around young stars) and would also lend support to models which link the triggering of rapid rise FU Orionis events to the existence of a close massive planet.

We present an analytical method for constraining the substellar initial mass function in globular clusters, based on the observed frequency of transit events. Globular clusters typically have very high stellar densities where close encounters are relatively common, and thus tidal capture can occur to form close binary systems. Encounters between main sequence stars and lower-mass objects can result in tidal capture if the mass ratio is > 0.01. If brown dwarfs exist in significant numbers, they too will be found in close binaries, and some fraction of their number should be revealed as they transit their stellar companions. We calculate the rate of tidal capture of brown dwarfs in both segregated and unsegregated clusters, and find that the tidal capture is more likely to occur over an initial relaxation time before equipartition occurs. The lack of any such transits in recent HST monitoring of 47 Tuc implies an upper limit on the frequency of brown dwarfs (< 15 % relative to stars) which is significantly below that measured in the galactic field and young clusters.

We present a model for the star formation process during the initial collapse of dark matter haloes at redshifts z=0-30. We derive a simple expression for the characteristic stellar mass scale during this initial burst of star formation. In our picture, this characteristic scale reflects both the minimum temperature to which the gas can cool (determined by the metallicity and the temperature of the cosmic microwave background) and the pressure of overlying baryons in the collapsing halo. This prescription reproduces both the large mass scales found in simulations of Population III star formation and the near solar values observed for star formation at low redshift.

We present results from high-resolution hydrodynamical simulations which explore the effects of small scale clustering in star-forming regions. A large ensemble of small-N clusters with 5 stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data. Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses which is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed IMF, both in the stellar and sub-stellar regime. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities > 2 km/s, which might deprive low mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long period (uncircularized) binaries. We predict that many close binary stars should have wide brown dwarfs companions. Brown dwarfs and, in general, very low mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or sub-stellar primaries being scarce. Where such binaries exist, they are either expected to be close enough (semi-major axis ~ 10 AU) to survive strong interactions with more massive binaries or else born in very small molecular cloud cores.

We compare evolutionary models for protoplanetary discs that include disc winds with observational determinations of the disc lifetime and accretion rate in Taurus. Using updated estimates for stellar ages in Taurus, together with published classifications, we show that the evolution of the disc fraction with stellar age is similar to that derived for ensembles of stars within young clusters. Around 30 percent of stars lose their discs within 1 Myr, while the remainder have disc lifetimes that are typically in the 1-10 Myr range. We show that the latter range of ages is consistent with theoretical models for disc evolution, provided that there is a dispersion of around 0.5 in the log of the initial disc mass. The same range of initial conditions brackets the observed variation in the accretion rate of Classical T Tauri stars at a given age. We discuss the expected lifetime of discs in close binary systems, and show that our models predict that the disc lifetime is almost constant for separations exceeding 10 au. This implies a low predicted fraction of binaries that pair a Classical T Tauri star with a Weak-lined T Tauri star, and is in better agreement with observations of the disc lifetime in binaries than disc models that do not include disc mass loss in a wind.

We present a full hydrodynamical investigation of a system of self-contained, interacting gas `cloudlets'. Calculations are performed using SPH, and cloudlets are allowed to collide, dissipate kinetic energy, merge and undergo local gravitational collapse. The numerical feasibility of such a technique is explored, and the resolution requirements examined. We examine the effect of the collision time and the velocity field on the timescale for evolution of the system, and show that almost all star formation will be confined to dense regions much smaller than the entire cloud. We also investigate the possibility, discussed by various authors, that such a system could lead to a power-law IMF by a process of repeated cloudlet coagulation, finding that in fact the inter-cloudlet collisions occur at too high a Mach number for merging to play an important part.

We have developed a simplified method of treating the radiative acceleration of dusty flows. This method retains the sharp impulse at the dust destruction radius that is a feature of frequency dependent radiative transfer, whilst placing minimal demands on computing resources. As such, it is suitable for inclusion in hydrodynamic codes. We have applied this method to the formation of massive stars in spherical geometry, and find that the fraction of a cloud which can accrete on to the central star is a strong function of the Jeans' Number and density profile of the cloud. Massive star formation is favoured by cold homogeneous conditions, as might result in regions where gas is swept up by some external triggering agent. We find (in contrast to previous authors) that massive star formation does not require a depleted dust mixture, although the use of dust at typical interstellar abundances does reduce the final stellar mass compared to cores formed from a depleted mixture.

We discuss a scenario in which brown dwarfs are formed like stars, except that their full collapse phases are interrupted through dynamical interactions in small multiple systems, leading to the ejection of the lightest member. This disintegration is a stochastic process, often resulting in the expulsion of newborn low mass stars, but when it occurs early enough the ejected stellar embryo will be a substellar object. This process may be so common at early ages that a large fraction of the ubiquitous brown dwarfs could have formed in this manner. Detailed gas dynamical simulations are required in order to better understand the details of the decay of small newborn multiple systems. We discuss the observational consequences of the ejection hypothesis, noting especially the importance of binaries with brown dwarf components as an observational test. Finally, we note that brown dwarfs that have recently been ejected may be so disturbed, by infall from the collapsing core and also by heavy accretion from perturbed circumstellar disks, that traditional spectral and luminosity criteria may fail to identify their substellar nature.

The spatial distribution of star formation within galaxies strongly affects the resulting feedback processes. Previous work has considered the case of a single, concentrated nuclear starburst, and also that of distributed single supernovae (SNe). Here, we consider ISM structuring by SNe originating in spatially distributed clusters having a cluster membership spectrum given by the observed HII region luminosity function. We show that in this case, the volume of HI cleared per SN is considerably greater than in either of the two cases considered hitherto. We derive a simple relationship between the ``porosity'' of the ISM and the star formation rate (SFR), and deduce a critical SFR_crit, at which the ISM porosity is unity. This critical value describes the case in which the SN mechanical energy output over a timescale t_e is comparable with the ISM ``thermal'' energy contained in random motions; t_e is the duration of SN mechanical input per superbubble. This condition also defines a critical gas consumption timescale t_exh, which for a Salpeter IMF and random velocities of \simeq 10 km s-1 is roughly 10e10 years. We draw a link between porosity and the escape of ionising radiation from galaxies, arguing that high escape fractions are expected if SFR >~ SFR_crit. The Lyman Break Galaxies, which are presumably subject to infall on a timescale < t_exh, meet this criterion, as is consistent with the significant leakage of ionising photons inferred in these systems. We suggest the utility of this simple parameterisation of escape fraction in terms of SFR for semi-empirical models of galaxy formation and evolution and for modeling mechanical and chemical feedback effects.

We use HST broad band images, VLA and VLBI continuum data to study the three dimensional orientation of jets relative to nuclear dust disks in 20 radio galaxies. The comparison between the position angles of the jets with those of the dust disks major axes shows a wide distribution, suggesting that they are not aligned preferentially perpendicular to each other. We use a statistical technique to determine the 3 dimensional distribution of angles between jets and dust disks rotation axes. This analysis shows that the observations are consistent with jets homogeneously distributed over a large region, extending over polar caps of 55 degrees to 77 degrees, but seem to avoid lying close to the plane of the dust disks. We argue that the lack of close alignment between jets and dust disks axes is not likely to be caused by feeding the nucleus with gas from mergers originated from random directions. We suggest that the misalignment can be due by a warping mechanism in the accretion disk, like self-irradiation instability or the Bardeen-Petterson effect, or that the gravitational potential in the inner regions of the galaxy is misaligned with respect to that of the dust disk.

We explore a mechanism for the formation of the first globular clusters, operating during the assembly of dwarf galaxies at high redshifts, z > 10. The substructure in the dark matter and the corresponding potential wells are responsible for setting the cluster scale of ~10^5 M_sun. The second mass scale in the formation problem, the stellar scale of ~1 M_sun, is determined in turn by the processes that cool the gas. We address the origin of the first, cluster scale by means of three-dimensional numerical simulations of the collapsing dark matter and gaseous components. We find that the gas falls into the deepest dark subhalos, resulting in a system of ~ 5 proto-globular clouds. The incipient globular clusters lose their individual dark halos in the process of violent relaxation, leading to the build-up of the general dark halo around the dwarf galaxy.

We present a series of statistical tests using homogeneous data and measurements for a sample of Seyfert galaxies. These galaxies were selected from mostly isotropic properties, their far infrared fluxes and warm infrared colors, which provide a considerable advantage over the criteria used by most investigators in the past, like ultraviolet excess. Our results provide strong support for a Unified Model in which Seyferts 2's contain a torus seen more edge-on than in Seyferts 1's and show that previous results showing the opposite were most likely due to selection effects.

Brown dwarfs may have such low masses because they are prematurely ejected from small unstable multiple systems, while the members are still actively building up their masses. We demonstrate that this scenario is consistent with all currently existing observations of brown dwarfs, and propose further observational tests. We review the status of the latest realistic numerical simulations of disintegrating small N clusters, which show that many of the ejected members end up with masses that are substellar, drifting away from their birth region with velocities rarely exceeding 2 km/s.

The current understanding of mechanical feedback is reviewed by evaluating the standard, adiabatic model for shell formation and evolution. This model is relevant to phenomena ranging from individual stellar-wind bubbles to galactic superwinds, forming the basis for our understanding of the multiphase ISM, IGM, and galactic evolutionary processes. Although significant discrepancies between the model and observation have been identified, to date there are none that require a fundamental revision. A variety of evidence, ranging over three orders of magnitude in spatial scale, is broadly consistent with the standard model. This includes kinematics of individual objects, observations of hot gas, the size distribution of HI shells, and outflow rates from starburst galaxies. However, some of the most pressing issues relating to shell evolution are still outstanding and obstruct efforts to resolve key questions like the fate of the hot gas.

We present a series of statistical tests done to a sample of 29 Seyfert 1 and 59 Seyfert 2 galaxies selected from mostly isotropic properties, their far infrared fluxes and warm infrared colors. Such selection criteria provide a profound advantage over the criteria used by most investigators in the past, such as ultraviolet excess. These tests were done using ground based high resolution VLA A-configuration 3.6 cm radio and optical B and I imaging data. From the relative number of Seyfert 1's and Seyfert 2's we calculate that the torus half opening angle is 48deg. We show that, as seen in previous papers, there is a lack of edge-on Seyfert 1 galaxies, suggesting dust and gas along the host galaxy disk probably play an important role in hiding some nuclei from direct view. We find that there is no statistically significant difference in the distribution of host galaxy morphological types and radio luminosities of Seyfert 1's and Seyfert 2's, suggesting that previous results showing the opposite may have been due to selection effects. The average extension of the radio emission of Seyfert 1's is smaller than that of Seyfert 2's by a factor of ~2-3, as predicted by the Unified Model. A search for galaxies around our Seyferts allows us to put a lower and an upper limit on the possible number of companions around these galaxies of 19% and 28%, respectively, with no significant difference in the number of companion galaxies between Seyfert 1's and Seyfert 2's. We also show that there is no preference for the radio jets to be aligned closer to the host galaxy disk axis in late type Seyferts, unlike results claimed by previous papers. These results, taken together, provide strong support for a Unified Model in which type 2 Seyferts contain a torus seen more edge-on than the torus in type 1 Seyferts.

Models of disc galaxies which invoke viscosity-driven radial flows have long been known to provide a natural explanation for the origin of stellar exponential discs, under the assumption that the star formation and viscous timescales are comparable. We present models which invoke simultaneous star formation, viscous redistribution of gas and cosmologically-motivated gaseous infall and explore the predictions such models make for the scale length evolution and radial star formation history of galactic stellar discs. While the inclusion of viscous flows is essential for ensuring that the stellar disc is always exponential over a significant range in radius, we find that such flows play essentially no role in determining the evolution of the disc scale length. In models in which the main infall phase precedes the onset of star formation and viscous evolution, we find the exponential scale length to be rather invariant with time. On the other hand, models in which star formation/viscous evolution and infall occur concurrently result in a smoothly increasing scale length with time, reflecting the mean angular momentum of material which has fallen in at any given epoch. The disc stellar populations in these models are predominantly young (ie. ages < 5 Gyr) beyond a few scale lengths. In both cases, viscous flows are entirely responsible for transporting material to very large radii. We discuss existing observational constraints on these models from studies of both local and moderate redshift disc galaxies. In particular, a good agreement is found between the solar neighbourhood star formation history predicted by our infall model and the recent observational determination of this quantity by Rocha-Pinto et al (2000).

We conjecture that brown dwarfs are substellar objects because they have been ejected from small newborn multiple systems which have decayed in dynamical interactions. In this view, brown dwarfs are stellar embryos for which the star formation process was aborted before the hydrostatic cores could build up enough mass to eventually start hydrogen burning. The disintegration of a small multiple system is a stochastic process, which can be described only in terms of the half-life of the decay. A stellar embryo competes with its siblings in order to accrete infalling matter, and the one that grows slowest is most likely to be ejected. With better luck, a brown dwarf would therefore have become a normal star. This interpretation of brown dwarfs readily explains the rarity of brown dwarfs as companions to normal stars (aka the ``brown dwarf desert''), the absence of wide brown dwarf binaries, and the flattening of the low mass end of the initial mass function. Possible observational tests of this scenario include statistics of brown dwarfs near Class 0 sources, and the kinematics of brown dwarfs in star forming regions while they still retain a kinematic signature of their expulsion. Because the ejection process limits the amount of gas brought along in a disk, it is predicted that substellar equivalents to the classical T Tauri stars should be very rare.

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses $\ms \simgreat 0.5 \solm$. The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a $\rho\propto R^{-2}$ distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of $\gamma=-3/2$ (where Salpeter is $\gamma=-2.35$). Once the stars dominate the potential and are virialised, which occurs first in the cluster core, the accretion radius is the Bondi-Hoyle radius. The resultant mass spectrum has an asymptotic limit of $\gamma=-2$ with slightly steeper slopes ($\gamma\approx-2.5$) if the stars are already mass-segr egated. Simulations of accretion onto clusters containing 1000 stars show that as expected, the low-mass stars accumulate the majority of their masses during the gas dominated phase whereas the high-mass stars accumulate the majority of their massed during the stellar dominated phase. This results in a mass spectrum with a relatively shallow $\gamma\approx 3/2$ power-law for low-mass stars and a steeper, power-law for high-mass stars $ -2.5\simless\gamma\le -2$. This competitive accretion model also results in a mass segregated cluster.

We investigate the physics of gas accretion in young stellar clusters. Accretion in clusters is a dynamic phenomenon as both the stars and the gas respond to the same gravitational potential. Accretion rates are highly non-uniform with stars nearer the centre of the cluster, where gas densities are higher, accreting more than others. This competitive accretion naturally results in both initial mass segregation and a spectrum of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal-lobe radius instead of the commonly used Bondi-Hoyle accretion radius. This works as both the stellar and gas velocities are under the influence of the same gravitational potential and are thus comparable. The low relative velocity that results means that the tidal radius is smaller than the Bondi-Hoyle radius in these systems. In contrast, when the stars dominate the potential and are virialised, the Bondi-Hoyle radius is smaller than the tidal radius and thus Bondi-Hoyle accretion is a better fit to the accretion rates.

We use numerical N-body simulations of the Orion Nebula Cluster (ONC) to investigate the destruction of protoplanetary disks by close stellar encounters and UV radiation from massive stars. The simulations model a cluster of 4000 stars, and we consider separately cases in which the disks have fixed radii of 100 AU and 10 AU. In the former case, depending on a star's position and orbit in the cluster over 10^7 years, UV photoevaporation removes at least 0.01 Msol from its disk, and can remove up to 1 Msol. We find no dynamical models of the ONC consistent with the suggestion of Storzer and Hollenbach that the observed distribution and abundance of proplyds could be explained by a population of stars on radial orbits which spend relatively little time near Theta 1C Ori (the most massive star in the ONC). Instead the observations require either massive disks (e.g. a typical initial disk mass of 0.4 Msol) or a very recent birth for Theta 1C Ori. When we consider the photoevaporation of the inner 10 AU of disks in the ONC, we find that planet formation would be hardly affected. Outside that region, planets would be prevented from forming in about half the systems, unless either the initial disk masses were very high or they formed in less than ~ 2 Myr and Theta 1C Ori has only very recently appeared. We also present statistics on the distribution of minimum stellar encounter separations. This peaks at 1000 AU, with less than 10% of stars having had an encounter closer than 100 AU after 10^7 years. We conclude that stellar encounters are unlikely to play a significant role in destroying protoplanetary disks. In the absence of any disruption mechanism other than those considered here, we would thus predict planetary systems like our own to be common amongst stars forming in ONC-like environments.

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.

We present the study of the relative angle between the accretion disk (or radio jet) and the galaxy disk for a sample of Seyfert galaxies selected from a mostly isotropic property, the 60$\mu$m flux. For each galaxy we have a pair of points ($i$,$\delta$), which are the inclination of the galaxy relative to the line of sight and the angle between the jet projected into the plane of the sky and the host galaxy major axis, respectively. This data is combined with a statistical technique to determine the distribution of $\beta$ angles {\it in 3 dimensions}, the angle between the jet and the host galaxy plane axis. The analysis of the data, not differentiating between Seyfert 1's and 2's, showed that the observed distribution of $i$ and $\delta$ values can be well represented by a homogeneous $\beta-$distribution in the range $0^{\circ}\leq\beta\leq90^{\circ}$. However, when we distinguished between Seyfert 1's and 2's, the models could not represent Seyfert 1's adequately. It was necessary to introduce viewing angle restrictions, that a galaxy can only be recognized as a Seyfert 1 if the angle between the jet and the line of sight ($|\phi|$) is smaller than a given angle $\phi_c$ and that the galaxy inclination $i$ is smaller than an angle $i_c$, in order to have statistically acceptable models. This is a direct and independent confirmation of the underlying concepts of the Unified Model. We discuss several ways to explain the misalignment between the accretion disk axis and the host galaxy disk axis.

We demonstrate that, by combining optical, radio and X-ray observations of a Seyfert, it is possible to provide a direct measurement of the angle $\beta$ between the direction of the radio jet and the normal to the plane of the spiral host galaxy. To do so, we make the assumptions that the inner radio jet is perpendicular to the X-ray observed inner accretion disk, and that the observed jet (or the stronger component, if the jet is two-sided) is physically closer to Earth than the plane of the galaxy. We draw attention to the possibility of measurement producing a result which is not self-consistent, in which case for that galaxy, one of the assumptions must fail.

We review recent work that investigates the formation of stellar clusters, ranging in scale from globular clusters through open clusters to the small scale aggregates of stars observed in T associations. In all cases, recent advances in understanding have been achieved through the use of state of the art stellar dynamical and gas dynamical calculations, combined with the possibility of intercomparison with an increasingly large dataset on young clusters. Among the subjects that are highlighted are the frequency of cluster-mode star formation, the possible relationship between cluster density and the highest stellar mass, subclustering and the dynamical interactions that occur in compact aggregates, such as binary star formation. We also consider how the spectrum of stellar masses may be shaped by the process of competitive accretion in dense environments and how cluster properties, such as mass segregation and cluster morphology, can be used in conjunction with numerical simulations to investigate the initial conditions for cluster formation. Lastly, we contrast bottom-up and top-down scenarios for cluster formation and discuss their applicability to the formation of clusters on a range of scales.

Using proper motion data for 894 stars in the Orion Nebula Cluster (ONC) compiled by Jones & Walker in 1988, we search for binaries with apparent separations in the range 1000-5000 AU, and find an upper limit of three. Using a Monte Carlo method, we test the consistency of this result with two hypotheses: i) that the cluster contains a binary population identical to that found in the solar neighbourhood, and ii) that the cluster contains no binaries at all in this separation range. We obtain results strongly favouring the latter hypothesis. Star formation in the Galaxy is seen to occur in a variety of different environments, but it has been proposed that most stars may be formed in dense regions similar to the ONC, rather than in less dense groupings like that found in Taurus-Auriga. Since roughly 15 per cent of galactic field stars are known to be in binaries with separations greater than 1000 AU, the apparent absence of such binaries in the ONC places an upper limit on the contribution that dense clusters can make to galactic star formation.

We present theoretical models for the evolution of T Tauri stars surrounded by circumstellar discs. The models include the effects of pre-main-sequence stellar and time dependent disc evolution, and incorporate the effects of stellar magnetic fields acting on the inner disc. For single stars, consistency with observations in Taurus-Auriga demands that disc dispersal occurs rapidly, on much less than the viscous timescale of the disc, at roughly the epoch when heating by stellar radiation first dominates over internal viscous dissipation. Applying the models to close binaries, we find that because the initial conditions for discs in binaries are uncertain, studies of extreme mass ratio systems are required to provide a stringent test of theoretical disc evolution models. We also note that no correlation of the infra-red colours of T Tauri stars with their rotation rate is observed, in apparent contradiction to the predictions of simple magnetospheric accretion models.

We study the interpretation of the mean surface density of stellar companions as a function of separation (or the two point correlation function of stars) in star-forming regions. First, we consider the form of the functions for various simple stellar distributions (binaries, global density profiles, clusters, and fractals) and the effects of survey boundaries. Following this, we study the dependencies of the separation at which a transition from the binary to the large-scale clustering regime occurs. Larson (1995) found that the mean surface density of companions follows different power-law functions of separation in the two regimes. He identified the transition separation with the typical Jeans length in the molecular cloud. However, we show that this is valid only for special cases. In general, the transition separation depends on the volume density of stars, the depth of the star-forming region, the volume-filling nature of the stellar distribution, and on the parameters of the binaries. Furthermore, the transition separation evolves with time. We then apply these results to the Taurus-Auriga, Ophiuchus, and Orion Trapezium star-forming regions. We find that while the transition separation in the Taurus-Auriga star-forming region may indicate a typical Jeans length, this is not true of the Orion Trapezium Cluster. We also show that there is no evidence for sub-clustering or fractal structure in the stars of the Orion Trapezium Cluster. This is consistent with the fact that, if such structure were present when the stars formed, it would have been erased by the current age of the cluster due to the stellar velocity dispersion. [Abstract abridged]

We use the standard, adiabatic shell evolution to predict the size distribution N(R) for populations of SN-driven superbubbles in a uniform ISM. We derive N(R) for simple cases of superbubble creation rate and mechanical luminosity function. We then compare our predictions for N(R) with the largely complete HI hole catalogue for the SMC, with a view toward the global structure of the ISM in that galaxy. We also present a preliminary derivation for N(v), the distribution of shell expansion velocities.

Observed variations in the HII region luminosity function (HII LF) seen in spiral arm vs. interarm regions, and different galactic Hubble type, can be explained simply by evolutionary effects and maximum number of ionizing stars per cluster. We present Monte Carlo simulations of the HII LF, drawing the number of ionizing stars N_* from a power-law distribution of constant slope, and the stellar masses from a Salpeter IMF with an upper-mass limit of 100 M_sol. We investigate the evolution of the HII LF, as determined by stellar main-sequence lifetimes and ionizing luminosities, for a single burst case and continuous creation of the nebular population. Shallower HII LF slopes measured for the arms of spiral galaxies can be explained as a composite slope, expected for a zero-age burst population, whereas the interarm regions tend to be dominated by evolved rich clusters described by a single, steeper slope. Steeper slopes in earlier-type galaxies can be explained simply by a lower maximum N_* cutoff found for the parent OB associations. The form of the HII LF can reveal features of the most recent (~< 10 Myr) star formation history in nearby galaxies.

We present a new method for the statistical investigation into the distributions of the angle beta between the radio axis and the normal to the galactic disk for a sample of Seyfert galaxies. We discuss how further observations of the sample galaxies can strengthen the conclusions. Our data are consistent with the hypothesis that AGN jets are oriented randomly in space, independent of the position of the plane of the galaxy. By making the simple assumption that the Standard Model of AGN holds, with a universal opening angle of the thick torus of phi_c, we demonstrate a statistical method to obtain an estimate of phi_c. Our data are not consistent with the simple-minded idea that Seyfert 1s and Seyfert 2s are differentiated solely by whether or not our line of sight lies within some fixed angle of the jet axis. Our result is significant on the 2 sigma level and can thus be considered only suggestive, not conclusive. A complete sample of Seyfert galaxies selected on an isotropic property is required to obtain a conclusive result.

We use the standard, adiabatic shell evolution to predict the size distribution N(R) for populations of OB superbubbles in a uniform ISM. We derive N(R) for simple cases of superbubble creation rate and mechanical luminosity function (MLF). For R < the characteristic radius R_e, N(R) is dominated by stalled objects, while for R>R_e it is dominated by growing objects. We also briefly investigate N(R) resulting from momentum-conserving shell evolution. We predict a peak in N(R) corresponding to individual SNRs. To estimate the MLF, we also examine evolutionary effects on the HII region luminosity function (HII LF), finding that for nebular luminosity fading as a power law in time, there is a minimum observed slope for the HII LFs. Comparison with the largely complete HI hole catalog for the SMC shows surprising agreement in the predicted and observed slope of N(R), suggesting that no other fundamental process is needed to explain the size distribution of shells in the SMC. Further comparison with largely incomplete HI data for M31, M33, and Holmberg II is also encouraging. We present expressions for the ISM porosity parameters, and estimate that they are substantially <1 for all of the galaxies except Holmberg II. Most of these galaxies therefore may not be strongly dominated by a hot interstellar component. However, porosity results for the Galaxy remain inconclusive.

We investigate the evolution of circumstellar discs around T Tauri stars that are ejected from small stellar clusters within molecular clouds. In particular, we study how the interaction that leads to ejection may hasten the transition between Classical and Weak-lined T Tauri status. In our models, ejections of T Tauri stars at velocities of 3-10 km/s truncate the accretion disc at radii between 1 and 10 a.u., reducing the viscous evolution time of the disc so that accretion rapidly ceases. The observational appearance of the resulting systems is then dependent on the presence or absence of a stellar magnetic field. For non-magnetic stars we find that a near-infra red excess should persist due to reprocessing of stellar radiation, but that this is greatly diminished for magnetic T Tauri stars by the presence of a magnetosphere extending to corotation. In either case, there is a period when ejected stars should appear as non-accreting systems with detectable circumstellar material at wavelengths of 5 microns and beyond. We discuss the implications of these results for models in which ejected stars contribute to the halo of pre-main-sequence objects discovered from ROSAT observations of star forming regions and the All-Sky Survey.

We construct models for the rotation rates of T Tauri stars whose spin is regulated by magnetic linkage between the star and a surrounding accretion disc. Our models utilise a time-dependent disc code to follow the accretion process and include the effects of pre-main-sequence stellar evolution. We find that the initial disc mass controls the evolution of the star-disc system. For sufficiently massive discs, a stellar field of $\sim$ 1 kG is able to regulate the spin rate to the observed values during the classical T Tauri phase. The field then acts to expel the disc and the star spins up at constant angular momentum as a weak-line system. Lower mass discs are ejected at an early epoch and fail to brake the star significantly. We extend the model to close binary systems, and find that the removal of angular momentum from the disc by the secondary significantly prolongs the inner disc lifetime. Such systems should therefore be relatively slow rotators. We also discuss the implications of our model for the spectral energy distributions and variability of T Tauri stars.

We show that outside-in disc instabilities, which can produce rapid-\null rise FU Orionis outbursts in T~Tauri systems, are a natural consequence of the existence of protoplanetary protostellar companions. The \pcomp{} would be formed through gravitational instability on a dynamical timescale ($\sim 10^5\y$), and provided it was formed at less than a few \AU{}, it would be swept in to $\sim15\Rsun$ in less than $10^5\y$. A companion of mass $10^{-2}$\Msun at this radius would act as a flood gate which stores material upstream and periodically releases it following the ignition of the thermal ionisation instability.

We investigate the response of a circumstellar accretion disc to the fly-by of a perturbing mass on a parabolic orbit. The energy and angular momentum transferred during the encounter are calculated using a reduced three-body method. In almost all close encounters the energy and angular momentum transfer is dominated by disc material becoming unbound from the system, with the contributions from close disc particle -- star encounters being significant. For more distant encounters with some prograde element to the motion the disc material loses energy and angular momentum to the perturber's orbit through a resonance feature. The magnitude of the energy transfer calculated in our simulations is greater than that of the binding energy of material exterior to periastron by a factor of two in the prograde case, and up to a factor of five in the case of the retrograde encounter. The destructive nature of the encounters indicates that a non-linear treatment is essential in all but the most distant encounters.

We demonstrate, using a simple analytic model, that the presence of a massive satellite can globally modify the structure and emission properties of an accretion disc to which it is tidally coupled. We show, using two levels of numerical approximation, that the analytic model gives reasonable results. The results are applicable to two astrophysical situations. In the case of an active galactic nucleus, we consider the case of a $\sim 10^3\Msun$ compact companion to the central black-hole and show that it could modulate the emitted spectrum on a timescale of $\sim10^5$ years. In the case of a T Tauri accretion disc, a satellite such as a sub-dwarf or giant planet could modify the disc spectral energy distribution over a substantial fraction of the T Tauri star lifetime.

We examine how accretion on to T Tauri stars may be modulated by a time-dependent `magnetic gate' where the inner edge of the accretion disc is disrupted by a varying stellar field. We show that magnetic field variations on time-scales shorter than 10^5 yr can modulate the accretion flow, thus providing a possible mechanism both for the marked photometric variability of T Tauri stars and for the possible conversion of T Tauri stars between classical and weak line status. We thus suggest that archival data relating to the spectrophotometric variability of T Tauri stars may provide an indirect record of magnetic activity cycles in low-mass pre-main-sequence stars.

The presence of protostellar disks can greatly increase the dissipation during close stellar encounters, leading to the formation of a significant population of binaries during the initial collapse and virialization of a cluster. We have used N-body simulations of collapsing globular clusters to find the major factors that determine the efficiency of binary formation through star-disk captures. This work serves the dual purpose of verifying the results of earlier analytic work as well as examining parameters not testable by that work. As in the earlier work, typical binary fractions of a few percent are found. For the parameters studied, the results are found to depend remarkably little upon disk evolution, the mass distribution of the stars, or their spatial distribution, though distributions in which the stars are highly clumped yield binary fractions larger by a factor of a few. The direct N-body integrations limit the models to relatively small values of N. Semiempirical relations are derived, however, which allow the results to be extrapolated to values of N appropriate to globular clusters.