The interaction between a YSO stellar magnetic field and its protostellar disc can result in stellar accretional flows and outflows from the inner disc rim. Gas flows with a velocity component perpendicular to disc midplane subject particles to centrifugal acceleration away from the protostar, resulting in particles being catapulted across the face of the disc. The ejected material can produce a "dust fan", which may be dense enough to mimic the appearance of a "puffed-up" inner disc rim. We derive analytic equations for the time dependent disc toroidal field, the disc magnetic twist, the size of the stable toroidal disc region, the jet speed and the disc region of maximal jet flow speed. We show how the observed infrared variability of the pre-transition disc system LRLL~31 can be modelled by a dust ejecta fan from the inner-most regions of the disc whose height is partially dependent on the jet flow speed. The greater the jet flow speed, the higher is the potential dust fan scale height. An increase in mass accretion onto the star tends to increase the height and optical depth of the dust ejection fan, increasing the amount of 1--8~$\mu$m radiation. The subsequent shadow reduces the amount of light falling on the outer disc and decreases the 8-- 40~$\mu$m radiation. A decrease in the accretion rate reverses this scenario, thereby producing the observed "see-saw" infrared variability.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

A small subset of young stellar objects (YSOs) exhibit "see-saw" temporal variations in their mid-infrared SED; as the flux short-ward of a fulcrum wavelength ($\lambda_{f}$) increases the flux long-wards of this wavelength decreases (and vice-versa) over timescales of weeks to years. While previous studies have shown that an opaque, axisymmetric occulter of variable height can cause this behaviour in the SED of these objects, the conditions under which a single $\lambda_{f}$ occurs have not previously been determined, nor the factors determining its value. Using radiative transfer modelling, we conduct a parametric study of the exemplar of this class, LRLL 31 to explore this phenomenon, and confirm that the cause of this flux variation is likely due to the change in height of the optically thick inner rim of the accretion disc at the dust sublimation radius, or some other phenomenon which results in a similar appearance. We also determine that a fulcrum wavelength only occurs for high inclinations, where the line of sight intersects the accretion disc. Accepting that the disc of LRLL 31 is highly inclined, the inner rim radius, radial and vertical density profiles are independently varied to gauge what effect this had on $\lambda_{f}$ and its position relative to the silicate feature near $10 \mu$m. While $\lambda_{f}$ is a function of each of these parameters, it is found to be most strongly dependent on the vertical density exponent $\beta$. All other factors being held constant, only for flatter discs ($\beta < 1.2$) did we find a $\lambda_{f}$ beyond the silicate feature.

This work presents a study of two Herbig Ae transitional discs, Oph IRS 48 and HD 169142; which both have reported rings in their dust density distributions. We use Keck-II/NIRC2 adaptive optics imaging observations in the L' filter (3.8 micron) to probe the regions of these discs inwards of ~20AU from the star. We introduce our method for investigating these transitional discs, which takes a forward modelling approach: making a model of the disc (using the Monte Carlo radiative transfer code RADMC), convolving it with point-spread functions of calibrator stars, and comparing the convolved models with the observational data. The disc surface density parameters are explored with a Monte Carlo Markov Chain technique. Our analysis recovers emission from both of the discs interior to the well known optically thick walls, modelled as a ring of emission at ~15AU in Oph IRS 48, and ~7AU for HD 169142, and identifies asymmetries in both discs. Given the brightness of the near-symmetric rings compared to the reported companion candidates, we suggest that the reported companion candidates can be interpreted as slightly asymmetric disc emission or illumination.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We present an approach that is able to both rapidly assess the dynamical stability of multiple planet systems, and determine whether an exoplanet system would be capable of hosting a dynamically stable Earth-mass companion in its habitable zone. We conduct a suite of numerical simulations using a swarm of massless test particles in the vicinity of the orbit of a massive planet, in order to develop a predictive tool which can be used to achieve these desired outcomes. In this work, we outline both the numerical methods we used to develop the tool, and demonstrate its use. We find that the test particles survive in systems either because they are unperturbed due to being so far removed from the massive planet, or due to being trapped in stable mean motion resonant orbits with the massive planet. The resulting unexcited test particle swarm produces a unique signature in (a,e) space that represents the stable regions within the system. We are able to scale and translate this stability signature, and combine several together in order to conservatively assess the dynamical stability of newly discovered multiple planet systems. We also assess the stability of a system's habitable zone and determine whether an Earth-mass companion could remain on a stable orbit, without the need for exhaustive numerical simulations.

In this work, we hunt for the best places to find exo-Earths in the currently known exoplanet population. While it is still unclear whether Jupiter had a beneficial or detrimental effect on the creation of the right environment for a habitable Earth to develop, we focus on the 51 multiple planet systems that have at least one Jupiter-like planet and aim to identify which would be good candidates to host an exo-Earth. We conduct a series of numerical simulations to identify dynamically stable regions of the habitable zone of the multiple exoplanet systems capable of hosting an Earth-mass planet. We produce a candidate list of 16 systems that could host such a stable exo-Earth in their habitable zone, and for which the induced radial velocity signal of a hypothetical one, two or four Earth-mass planet on the host star would be detectable with the ESPRESSO spectrograph. We find that whilst the gravitational interactions with the massive planet nearest the habitable zone are critical in determining stability, the secular resonant interactions between multiple planets can also have a dramatic influence on the overall stability of the habitable zone.

While the number of exoplanets discovered continues to increase at a rapid rate, we are still to discover any system that truly resembles the Solar system. Existing and near future surveys will likely continue this trend of rapid discovery. To see if these systems are Solar system analogues, we will need to efficiently allocate resources to carry out intensive follow-up observations. We seek to uncover the properties and trends across systems that indicate how much of the habitable zone is stable in each system to provide focus for planet hunters. We study the dynamics of all known single Jovian planetary systems, to assess the dynamical stability of the habitable zone around their host stars. We perform a suite of simulations of all systems where the Jovian planet will interact gravitationally with the habitable zone, and broadly classify these systems. Besides the system's mass ratio (Mpl/Mstar), and the Jovian planet's semi-major axis (apl) and eccentricity (epl), we find that there are no underlying system properties which are observable that indicate the potential for planets to survive within the system's habitable zone. We use Mpl/Mstar, apl and epl to generate a parameter space over which the unstable systems cluster, thus allowing us to predict which systems to exclude from future observational or numerical searches for habitable exoplanets. We also provide a candidate list of 20 systems that have completely stable habitable zones and Jovian planets orbiting beyond the habitable zone as potential first order Solar system analogues.

Planet formation is thought to occur in discs around young stars by the aggregation of small dust grains into much larger objects. The growth from grains to pebbles and from planetesimals to planets is now fairly well understood. The intermediate stage has however been found to be hindered by the radial-drift and fragmentation barriers. We identify a powerful mechanism in which dust overcomes both barriers. Its key ingredients are i) backreaction from the dust onto the gas, ii) grain growth and fragmentation, and iii) large-scale gradients. The pile-up of growing and fragmenting grains modifies the gas structure on large scales and triggers the formation of pressure maxima, in which particles are trapped. We show that these self-induced dust traps are robust: they develop for a wide range of disc structures, fragmentation thresholds and initial dust-to-gas ratios. They are favored locations for pebbles to grow into planetesimals, thus opening new paths towards the formation of planets.

With continued improvement in telescope sensitivity and observational techniques, the search for rocky planets in stellar habitable zones is entering an exciting era. With so many exoplanetary systems available for follow-up observations to find potentially habitable planets, one needs to prioritise the ever-growing list of candidates. We aim to determine which of the known planetary systems are dynamically capable of hosting rocky planets in their habitable zones, with the goal of helping to focus future planet search programs. We perform an extensive suite of numerical simulations to identify regions in the habitable zones of single Jovian planet systems where Earth mass planets could maintain stable orbits, specifically focusing on the systems in the Catalog of Earth-like Exoplanet Survey Targets (CELESTA). We find that small, Earth-mass planets can maintain stable orbits in cases where the habitable zone is largely, or partially, unperturbed by a nearby Jovian, and that mutual gravitational interactions and resonant mechanisms are capable of producing stable orbits even in habitable zones that are significantly or completely disrupted by a Jovian. Our results yield a list of 13 single Jovian planet systems in CELESTA that are not only capable of supporting an Earth-mass planet on stable orbits in their habitable zone, but for which we are also able to constrain the orbits of the Earth-mass planet such that the induced radial velocity signals would be detectable with next generation instruments.

The presence of dusty debris around main sequence stars denotes the existence of planetary systems. Such debris disks are often identified by the presence of excess continuum emission at infrared and (sub-)millimetre wavelengths, with measurements at longer wavelengths tracing larger and cooler dust grains. The exponent of the slope of the disk emission at sub-millimetre wavelengths, `q', defines the size distribution of dust grains in the disk. This size distribution is a function of the rigid strength of the dust producing parent planetesimals. As part of the survey `PLAnetesimals around TYpical Pre-main seqUence Stars' (PLATYPUS) we observed six debris disks at 9-mm using the Australian Telescope Compact Array. We obtain marginal (~3-\sigma) detections of three targets: HD 105, HD 61005, and HD 131835. Upper limits for the three remaining disks, HD20807, HD109573, and HD109085, provide further constraint of the (sub-)millimetre slope of their spectral energy distributions. The values of q (or their limits) derived from our observations are all smaller than the oft-assumed steady state collisional cascade model (q = 3.5), but lie well within the theoretically expected range for debris disks q ~ 3 to 4. The measured q values for our targets are all < 3.3, consistent with both collisional modelling results and theoretical predictions for parent planetesimal bodies being `rubble piles' held together loosely by their self-gravity.

Protoplanetary disc systems observed at radio wavelengths often show excess emission above that expected from a simple extrapolation of thermal dust emission observed at short millimetre wavelengths. Monitoring the emission at radio wavelengths can be used to help disentangle the physical mechanisms responsible for this excess, including free-free emission from a wind or jet, and chromospheric emission associated with stellar activity. We present new results from a radio monitoring survey conducted with Australia Telescope Compact Array over the course of several years with observation intervals spanning days, months and years, where the flux variability of 11 T Tauri stars in the Chamaeleon and Lupus star forming regions was measured at 7 and 15 mm and 3 and 6 cm. Results show that for most sources are variable to some degree at 7 mm, indicating the presence of emission mechanisms other than thermal dust in some sources. Additionally, evidence of grain growth to cm-sized pebbles was found for some sources that also have signs of variable flux at 7 mm. We conclude that multiple processes contributing to the emission are common in T Tauri stars at 7 mm and beyond, and that a detection at a single epoch at radio wavelengths should not be used to determine all processes contributing to the emission.

Radial and azimuthal features (such as disc offsets and eccentric rings) seen in high resolution images of debris discs, provide us with the unique opportunity of finding potential planetary companions which betray their presence by gravitationally sculpting such asymmetric features. The young debris disc around HD 115600, imaged recently by the Gemini Planet Imager, is such a disc with an eccentricity 0.1<e<0.2 and a projected offset from the star of 4 AU. Using our modified N-body code which incorporates radiation forces, we firstly aim to determine the orbit of a hidden planetary companion potentially responsible for shaping the disc. We run a suite of simulations covering a broad range of planetary parameters using a Monte Carlo Markov Chain sampling method and create synthetic images from which we extract the geometric disc parameters to be compared with the observed and model-derived quantities. We then repeat the study using a traditional grid to explore the planetary parameter space and aim secondly to compare the efficiency of both sampling methods. We find a planet of 7.8 Jupiter mass orbiting at 30 AU with an eccentricity of e = 0.2 to be the best fit to the observations of HD 115600. Technically, such planet has a contrast detectable by direct imaging, however the system's orientation does not favour such detection. In this study, at equal number of explored planetary configurations, the Monte Carlo Markov Chain not only converges faster but provides a better fit than a traditional grid.

Context: Gaps, cavities and rings in circumstellar disks are signposts of disk evolution and planet-disk interactions. We follow the recent suggestion that Herbig Ae/Be disks with a flared disk harbour a cavity, and investigate the disk around HD~97048. Aims: We aim to resolve the 34$\pm$ 4 au central cavity predicted by Maaskant et al. (2013) and to investigate the structure of the disk. Methods: We image the disk around HD~97048 using ALMA at 0.85~mm and 2.94~mm, and ATCA (multiple frequencies) observations. Our observations also include the 12CO J=1-0, 12CO J=3-2 and HCO+ J=4-3 emission lines. Results: A central cavity in the disk around HD~97048 is resolved with a 40-46 au radius. Additional radial structure present in the surface brightness profile can be accounted for either by an opacity gap at ~90 au or by an extra emitting ring at ~150 au. The continuum emission tracing the dust in the disk is detected out to 355 au. The 12CO J=3-2 disk is detected 2.4 times farther out. The 12CO emission can be traced down to $\approx$ 10 au scales. Non-Keplerian kinematics are detected inside the cavity via the HCO+ J=4-3 velocity map. The mm spectral index measured from ATCA observations suggests that grain growth has occurred in the HD~97048 disk. Finally, we resolve a highly inclined disk out to 150 au around the nearby 0.5~$M_{\odot}$ binary ISO-ChaI 126. Conclusions: The data presented here reveal a cavity in the disk of HD 97048, and prominent radial structure in the surface brightness. The cavity size varies for different continuum frequencies and gas tracers. The gas inside the cavity follows non-Keplerian kinematics seen in HCO+ emission. The variable cavity size along with the kinematical signature suggests the presence of a substellar companion or massive planet inside the cavity.

We model the behaviour of dust grains entrained by photoevaporation-driven winds from protoplanetary discs assuming a non-rotating, plane-parallel disc. We obtain an analytic expression for the maximum entrainable grain size in extreme-UV radiation-driven winds, which we demonstrate to be proportional to the mass loss rate of the disc. When compared with our hydrodynamic simulations, the model reproduces almost all of the wind properties for the gas and dust. In typical turbulent discs, the entrained grain sizes in the wind are smaller than the theoretical maximum everywhere but the inner disc due to dust settling.

We investigate dust entrainment by photoevaporative winds in protoplanetary discs using dusty smoothed particle hydrodrodynamics (SPH). We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in EUV driven winds. The maximum grain size is set by dust settling in the disc rather than aerodynamic drag in the wind. More generally, there is a linear relationship between the base flow density and the maximum entrainable grain size in the wind. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between $10$--$20\,$AU. The peak dust density for each grain size occurs at the inner edge of its own entrainment region.

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and $\sim5$ $\mu$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)\propto a^{-q}$, of $\langle q \rangle = 3.36\pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

Resolved debris disc images can exhibit a range of radial and azimuthal structures, including gaps and rings, which can result from planetary companions shaping the disc by their gravitational influence. Currently there are no tools available to determine the architecture of potential companions from disc observations. Recent work by Rodigas et al. (2014) presents how one can estimate the maximum mass and minimum semi major axis of a hidden planet empirically from the width of the disc in scattered light. In this work, we use the predictions of Rodigas et al. applied to two debris discs HD 202628 and HD 207129. We aim to test if the predicted orbits of the planets can explain the features of their debris disc, such as eccentricity and sharp inner edge. We first run dynamical simulations using the predicted planetary parameters of Rodigas et al., and then numerically search for better parameters. Using a modified N-body code including radiation forces, we perform simulations over a broad range of planet parameters and compare synthetics images from our simulations to the observations. We find that the observational features of HD 202628 can be reproduced with a planet five times smaller than expected, located 30 AU beyond the predicted value, while the best match for HD 207129 is for a planet located 5-10 AU beyond the predicted location with a smaller eccentricity. We conclude that the predictions of Rodigas et al. provide a good starting point but should be complemented by numerical simulations.

Infrared observations provide the dust composition in the protoplanetary discs surface layers, but can not probe the dust chemistry in the midplane, where planet formation occurs. Meteorites show that dynamics was important in determining the dust distribution in the Solar Nebula and needs to be considered if we are to understand the global chemistry in discs. 1D radial condensation sequences can only simulate one disc layer at a time and cannot describe the global chemistry or the complexity of meteorites. To address these limitations, we compute for the first time the two dimensional distribution of condensates in the inner Solar Nebula using a thermodynamic equilibrium model, and derive timescales for vertical settling and radial migration of dust. We find two enstatite-rich zones within 1 AU from the young Sun: a band ~0.1 AU thick in the upper optically-thin layer of the disc interior to 0.8 AU, and in the optically-thick disc midplane out to ~0.4 AU. The two enstatite-rich zones support recent evidence that Mercury and enstatite chondrites shared a bulk material with similar composition. Our results are also consistent with infrared observation of protoplanetary disc which show emission of enstatite-rich dust in the inner surface of discs. The resulting chemistry and dynamics suggests that the formation of the bulk material of enstatite chondrites occurred in the inner surface layer of the disc, within 0.4~AU. We also propose a simple alternative scenario in which gas fractionation and vertical settling of the condensates lead to an enstatite-chondritic bulk material.

We present ATCA continuum observations at a wavelength of 6.8 mm of five debris disks: $\beta$ Pictoris, q$^1$ Eridani, HD 107146, HD 181327, and HD 95086. These observations provide the detection at the longest wavelengths obtained to date for all these debris disks. By combining our 6.8 mm data with previous detections at shorter sub-millimeter/millimeter wavelengths we measure the long wavelength spectral index of these sources. We then use previous estimates for the temperature of the emitting dust to derive the spectral index of the dust emissivity. Under the assumption that all the detected flux comes from dust only, we constrain the slope of the solid size distribution, assumed to be a power-law. The values that we infer for the slope of the size distribution range between about 3.36 and 3.50. We compare our findings with the case of the Fomalhaut debris disk and use these results to test the predictions of collisional cascades of planetesimal belts.

The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the north-west of the central star, which matches the VLA clump position. The compactness of the Ka emission is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution with parameters adequately selected to reproduce the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarised emission. Continuum observations at a different wavelength are necessary to conclude if the VLA-ALMA difference is an opacity or a real dust segregation.

Numerical simulations are a crucial tool to understand the relationship between debris discs and planetary companions. However, simulations throughout the literature have been conducted with various initial conditions often with little or no justification. In this paper, we aim to study the dependence on the initial conditions of N-body simulations modelling the interaction between a massive and eccentric planet on an exterior debris disc. To achieve this, we first classify three broad approaches used in the literature and provide some physical context for when each category should be used. We then run a series of N-body simulations, that include radiation forces acting on small grains, with varying initial conditions across the three categories. We test the influence of the initial parent body belt width, eccentricity, and alignment with the planet on the resulting debris disc structure and compare the final peak emission location, disc width and offset of synthetic disc images produced with a radiative transfer code. We also track the evolution of the forced eccentricity of the dust grains induced by the planet, as well as resonance dust trapping. We find that an initially broad parent body belt always results in a broader debris disc than an initially narrow parent body belt. While simulations with a parent body belt with low initial eccentricity (e ~ 0) and high initial eccentricity (0 < e < 0.3) resulted in similar broad discs, we find that purely secular forced initial conditions, where the initial disc eccentricity is set to the forced value and the disc is aligned with the planet, always result in a narrower disc. We conclude that broad debris discs can be modelled by using either a dynamically cold or dynamically warm parent belt, while in contrast eccentric narrow debris rings are reproduced using a secularly forced parent body belt.

Protoplanetary discs are now routinely observed and exoplanets, after the numerous indirect discoveries, are starting to be directly imaged. To better understand the planet formation process, the next step is the detection of forming planets or of signposts of young planets still in their disc, such as gaps. A spectacular example is the ALMA science verification image of HL Tau showing numerous gaps and rings in its disc. To study the observability of planet gaps, we ran 3D hydrodynamical simulations of a gas and dust disc containing a 5 M J gap-opening planet and characterised the spatial distribution of migrating, growing and fragmenting dust grains. We then computed the corresponding synthetic images for ALMA. For a value of the dust fragmentation threshold of 15 m s --1 for the collisional velocity, we identify for the first time a self-induced dust pile up in simulations taking fragmentation into account. This feature, in addition to the easily detected planet gap, causes a second apparent gap that could be mistaken for the signature of a second planet. It is therefore essential to be cautious in the interpretation of gap detections.

We present observations of Epsilon Eridani from the Submillimeter Array (SMA) at 1.3 millimeters and from the Australia Telescope Compact Array (ATCA) at 7 millimeters that reach an angular resolution of ~4" (13 AU). These first millimeter interferometer observations of Epsilon Eridani, which hosts the closest debris disk to the Sun, reveal two distinct emission components: (1) the well-known outer dust belt, which, although patchy, is clearly resolved in the radial direction, and (2) an unresolved source coincident with the position of the star. We use direct model-fitting of the millimeter visibilities to constrain the basic properties of these two components. A simple Gaussian shape for the outer belt fit to the SMA data results in a radial location of $64.4^{+2.4}_{-3.0}$ AU and FWHM of $20.2^{+6.0}_{-8.2}$ AU (fractional width $\Delta R/R = 0.3$. Similar results are obtained taking a power law radial emission profile for the belt, though the power law index cannot be usefully constrained. Within the noise obtained (0.2 mJy/beam), these data are consistent with an axisymmetric belt model and show no significant azimuthal structure that might be introduced by unseen planets in the system. These data also limit any stellocentric offset of the belt to $<9$ AU, which disfavors the presence of giant planets on highly eccentric ($>0.1$) and wide (10's of AU) orbits. The flux density of the unresolved central component exceeds predictions for the stellar photosphere at these long wavelengths, by a marginally significant amount at 1.3 millimeters but by a factor of a few at 7 millimeters (with brightness temperature $13000 \pm 1600$ K for a source size of the optical stellar radius). We attribute this excess emission to ionized plasma from a stellar corona or chromosphere.

There is much evidence that planet formation is occurring in the disk around the Herbig Be star HD100546. To learn more about the processes occurring in this disk we conducted high resolution imaging at 43/45 GHz with the Australia Telescope Compact Array (ATCA). Multiple array configurations were used, providing a best spatial resolution of $\sim$ 0.15 arcsec, or 15 AU at HD100546's distance of $\sim$ 100 pc. Significant structure is revealed, but its precise form is dependent on the $u-v$ plane sampling used for the image reconstruction. At a resolution of $\leq$ 30 AU we detected an inner gap in the disk with a radius of $\sim$ 25 AU and a position angle approximately along the known disk major axis. With different weighting, and an achieved resolution of $\sim$ 15 AU, emission appears at the centre and the disk takes on the shape of an incomplete ring, much like a horseshoe, again with a gap radius of $\sim$ 25 AU. The position angle of the disk major axis and its inclination from face-on are determined to be $140^{\circ}\pm5^{\circ}$ and $40^{\circ}\pm5^{\circ}$ respectively. The $\sim$ 25 AU gap radius is confirmed by a null in the real part of the binned visibilities at 320$\pm$10 k$\lambda$, whilst the non-axisymmetric nature is also confirmed through significant structure in the imaginary component. The emission mechanism at the central peak is most likely to be free-free emission from a stellar or disk wind. Overall our data support the picture of at least one, but probably several, giant planets orbiting HD100546 within 25 AU.

We model the dust evolution in protoplanetary disks with full 3D, Smoothed Particle Hydrodynamics (SPH), two-phase (gas+dust) hydrodynamical simulations. The gas+dust dynamics, where aerodynamic drag leads to the vertical settling and radial migration of grains, is consistently treated. In a previous work, we characterized the spatial distribution of non-growing dust grains of different sizes in a disk containing a gap-opening planet and investigated the gap's detectability with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we take into account the effects of grain growth and fragmentation and study their impact on the distribution of solids in the disk. We show that rapid grain growth in the two accumulation zones around planet gaps is strongly affected by fragmentation. We discuss the consequences for ALMA observations.

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply confined and furthermore the trapped grains should follow a distribution that is markedly different from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD142527, from 34GHz to 700GHz, that reveal a compact concentration of ~cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of ~mm-sized particles. The emission peaks at wavelengths shorter than ~1mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reflect the gas distribution.

With improvements in exoplanet detection techniques, the number of multiple planet systems discovered is increasing, while the detection of potentially habitable Earth-mass planets remains complicated and thus requires new search strategies. Dynamical studies of known multiple planet systems are therefore a vital tool in the search for stable and habitable planet candidates. Here, we present a dynamical study of the three-planet system HD 204313 to determine whether it could harbour an Earth-like planet within its habitable zone for a sufficient time to develop life. We found two semi-stable regions in the system, but neither prove stable for long enough for a terrestrial planet to develop life. Our investigations suggest that overlapping weak and high order resonances may be responsible for these semi-stable regions. This study established a framework for a larger project that will study the dynamical stability of the habitable zone of multiple planet systems, providing a list of interesting targets for future habitable low-mass planet searches.

We aim to derive a simple analytic model to understand the essential properties of vertically settling growing dust grains in laminar protoplanetary discs. Separating the vertical dynamics from the motion in the disc midplane, we integrate the equations of motion for both a linear and an exponential grain growth rate. Numerical integrations are performed for more complex growth models. We find that the settling efficiency depends on the value of the dimensionless parameter gamma, which characterises the relative efficiency of grain growth with respect to the gas drag. Since gamma is expected to be of order as the initial dust-to-gas ratio in the disc (of order 10^-2), grain growth enhances the energy dissipation of the dust particles and improve the settling efficiency in protoplanetary discs. This behaviour is mostly independent of the growth model considered as well as of the radial drift of the particles.

In a series of papers, we present a comprehensive analytic study of the global motion of growing dust grains in protoplanetary discs, addressing both the radial drift and the vertical settling of the particles. Here we study how the radial drift of dust particles is affected by grain growth. In a first step, toy models in which grain growth can either be constant, accelerate or decelerate are introduced. The equations of motion are analytically integrable and therefore the grains dynamics is easy to understand. The radial motion of growing grains is governed by the relative efficiency of the growth and migration processes which is expressed by the dimensionless parameter Lambda, as well as the exponents for the gas surface density and temperature profiles, denoted p and q respectively. When Lambda is of order unity, growth and migration are strongly coupled, providing the most efficient radial drift. For the toy models considered, grains pile up when -p+q+1/2<0. Importantly, we show the existence of a second process which can help discs to retain their solid materials. For accelerating growth, grains end up their migration at a finite radius, thus avoiding being accreted onto the central star.

Star formation does not occur until the onset of gravitational collapse inside giant molecular clouds. However, the conditions that initiate cloud collapse and regulate the star formation process remain poorly understood. Local processes such as turbulence and magnetic fields can act to promote or prevent collapse. On larger scales, the galactic potential can also influence cloud stability and is traditionally assessed by the tidal and shear effects. In this paper, we examine the stability of giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) against shear and the galactic tide using CO data from the Magellanic Mopra Assessment (MAGMA) and rotation curve data from the literature. We calculate the tidal acceleration experienced by individual GMCs and determine the minimum cloud mass required for tidal stability. We also calculate the shear parameter, which is a measure of a clouds susceptibility to disruption via shearing forces in the galactic disk. We examine whether there are correlations between the properties and star forming activity of GMCs and their stability against shear and tidal disruption. We find that the GMCs are in approximate tidal balance in the LMC, and that shear is unlikely to affect their further evolution. GMCs with masses close to the minimal stable mass against tidal disruption are not unusual in terms of their mass, location, or CO brightness, but we note that GMCs with large velocity dispersion tend to be more sensitive to tidal instability. We also note that GMCs with smaller radii, which represent the majority of our sample, tend to more strongly resist tidal and shear disruption. Our results demonstrate that star formation in the LMC is not inhibited by to tidal or shear instability.

Among the numerous known extrasolar planets, only a handful have been imaged directly so far, at large orbital radii and in rather evolved systems. The Atacama Large Millimeter/submillimeter Array (ALMA) will have the capacity to observe these wide planetary systems at a younger age, thus bringing a better understanding of the planet formation process. Here we explore the ability of ALMA to detect the gaps carved by planets on wide orbits.

A number of objects in primitive meteorites have oxygen isotopic compositions that place them on a distinct, mass-independent fractionation line with a slope of one on a three-isotope plot. The most popular model for describing how this fractionation arose assumes that CO self-shielding produced 16O-rich CO and 16O-poor H2O, where the H2O subsequently combined with interstellar dust to form relatively 16O-poor solids within the Solar Nebula. Another model for creating the different reservoirs of 16O-rich gas and 16O-poor solids suggests that these reservoirs were produced by Galactic chemical evolution (GCE) if the Solar System dust component was somewhat younger than the gas component and both components were lying on the line of slope one in the O three-isotope plot. We argue that GCE is not the cause of mass-independent fractionation of the oxygen isotopes in the Solar System. The GCE scenario is in contradiction with observations of the 18O/17O ratios in nearby molecular clouds and young stellar objects. It is very unlikely for GCE to produce a line of slope one when considering the effect of incomplete mixing of stellar ejecta in the interstellar medium. Furthermore, the assumption that the Solar System dust was younger than the gas requires unusual timescales or the existence of an important stardust component that is not theoretically expected to occur nor has been identified to date.

The composition of the most primitive Solar System condensates, such as calcium-aluminum-rich inclusions (CAI) and micron-sized corundum grains, show that short-lived radionuclides (SLR), e.g., 26Al, were present in the early Solar System. Their abundances require a local origin, which however is far from being understood. We present for the first time the abundances of several SLR up to 60Fe predicted from stars with initial mass in the range roughly 7-11 Msun. These stars evolve through core H, He, and C burning. After core C burning they go through a "Super"-asymptotic giant branch (Super-AGB) phase, with the H and He shells activated alternately, episodic thermal pulses in the He shell, a very hot temperature at the base of the convective envelope (~ 10^8 K), and strong stellar winds driving the H-rich envelope into the surrounding interstellar medium. The final remnants of the evolution of Super-AGB stars are mostly O-Ne white dwarfs. Our Super-AGB models produce 26Al/27Al yield ratios ~ 0.02 - 0.26. These models can account for the canonical value of the 26Al/27Al ratio using dilutions with the Solar Nebula of the order of 1 part of Super-AGB mass per several 10^2 to several 10^3 of Solar Nebula mass, resulting in associated changes in the O composition in the range 3 permil to 20 permil. This is in agreement with observations of the O isotopic ratios in primitive Solar System condensated, which do not carry the signature of a stellar pollutor. The radionuclides 41Ca and 60Fe are produced by neutron captures in Super-AGB stars and their meteoritic abundances are also matched by some of our models, depending on the nuclear and stellar physics uncertainties as well as the meteoritic experimental data. We also expect and are currently investigating Super-AGB production of SLR heavier than iron, such as 107Pd.

[Abridged] Aims: We provide predictions for ALMA observations of planet gaps that account for the specific spatial distribution of dust that results from consistent gas+dust dynamics. Methods: In a previous work, we ran full 3D, two-fluid Smoothed Particle Hydrodynamics (SPH) simulations of a planet embedded in a gas+dust T Tauri disk for different planet masses and grain sizes. In this work, the resulting dust distributions are passed to the Monte Carlo radiative transfer code MCFOST to construct synthetic images in the ALMA wavebands. We then use the ALMA simulator to produce images that include thermal and phase noise for a range of angular resolutions, wavelengths, and integration times, as well as for different inclinations, declinations and distances. We also produce images which assume that gas and dust are well mixed with a gas-to-dust ratio of 100 to compare with previous ALMA predictions, all made under this hypothesis. Results: Our findings clearly demonstrate the importance of correctly incorporating the dust dynamics. We show that the gap carved by a 1 M_J planet orbiting at 40 AU is visible with a much higher contrast than the well-mixed assumption would predict. In the case of a 5 M_J planet, we clearly see a deficit in dust emission in the inner disk, and point out the risk of interpreting the resulting image as that of a transition disk with an inner hole if observed in unfavorable conditions. Planet signatures are fainter in more distant disks but declination or inclination to the line-of-sight have little effect on ALMA's ability to resolve the gaps. Conclusions: ALMA has the potential to see signposts of planets in disks of nearby star-forming regions. We present optimized observing parameters to detect them in the case of 1 and 5 M_J planets on 40 AU orbits.

As gas flowed from the solar accretion disk or solar nebula onto the proto-Sun, magnetic pressure gradients in the solar magnetosphere and the inner solar nebula provided an environment where some of this infalling flow was diverted to produce a low pressure, high temperature, gaseous, "infall" atmosphere around the inner solar nebula. The pressure in this inner disk atmosphere was mainly dependant on the accretion flow rate onto the star. High flow rates implied relatively high pressures, which decreased over time as the accretion rate decreased. In the first hundred thousand years after the formation of the solar nebula, accretional flow gas pressures were high enough to create submicron-sized Refractory Metal Nuggets (RMNs) - the precursors to Calcium Aluminum Inclusions (CAIs). Optimal temperatures and pressures for RMN formation may have occurred between 20,000 to 100,000 years after the formation of the solar nebula. It is possible that conditions were conducive to RMN/CAI formation over an eighty thousand year timescale. The "infall" atmosphere and the condensation of refractory particles within this atmosphere may be observable around the inner disks of other protostellar systems. The interaction of forces from magnetic fields with the radiation pressure from the proto-Sun and the inner solar accretion disk potentially produced an optical-magnetic trap above and below the inner solar nebula, which provided a relatively stable environment in which the RMNs/proto-CAIs could form and grow. These RMN formation sites only existed during accretion events from the proto-solar disk onto the proto-Sun. As such, the formation and growth time of a particular RMN was dependent on the timescale of its nascent accretion event.

We present ATCA results of a 3 and 7 mm continuum survey of 20 T Tauri stars in the Chamaeleon and Lupus star forming regions. This survey aims to identify protoplanetary discs with signs of grain growth. We detected 90% of the sources at 3 and 7 mm, and determined the spectral slopes, dust opacity indices and dust disc masses. We also present temporal monitoring results of a small sub-set of sources at 7, 15 mm and 3+6 cm to investigate grain growth to cm sizes and constrain emission mechanisms in these sources. Additionally, we investigated the potential correlation between grain growth signatures in the infrared (10 \mu m silicate feature) and millimetre (1-3 mm spectral slope, {\alpha}). Eleven sources at 3 and 7 mm have dominant thermal dust emission up to 7 mm, with 7 of these having a 1-3 mm dust opacity index less than unity, suggesting grain growth up to at least mm sizes. The Chamaeleon sources observed at 15 mm and beyond show the presence of excess emission from an ionised wind and/or chromo- spheric emission. Long-timescale monitoring at 7 mm indicated that cm-sized pebbles are present in at least four sources. Short-timescale monitoring at 15 mm suggests the excess emission is from thermal free-free emission. Finally, a weak correlation was found between the strength of the 10 \mum feature and {\alpha}, suggesting simultaneous dust evolution of the inner and outer parts of the disc. This survey shows that grain growth up to cm-sized pebbles and the presence of excess emission at 15 mm and beyond are common in these systems, and that temporal monitoring is required to disentangle these emission mechanisms.

We present a new catalogue of 5,106 infrared bubbles created through visual classification via the online citizen science website 'The Milky Way Project'. Bubbles in the new catalogue have been independently measured by at least 5 individuals, producing consensus parameters for their position, radius, thickness, eccentricity and position angle. Citizen scientists - volunteers recruited online and taking part in this research - have independently rediscovered the locations of at least 86% of three widely-used catalogues of bubbles and H ii regions whilst finding an order of magnitude more objects. 29% of the Milky Way Project catalogue bubbles lie on the rim of a larger bubble, or have smaller bubbles located within them, opening up the possibility of better statistical studies of triggered star formation. Also outlined is the creation of a 'heat map' of star-formation activity in the Galactic plane. This online resource provides a crowd-sourced map of bubbles and arcs in the Milky Way, and will enable better statistical analysis of Galactic star-formation sites.

We present new spatially resolved observations of the dust thermal emission at 7 mm from the Fomalhaut debris disk obtained with the Australia Telescope Compact Array. These observations provide the longest wavelength detection of the Fomalhaut debris disk to date. We combined the new data to literature sub-mm data to investigate the spectral index of the dust thermal emission in the sub-millimeter and constrained the $q$-slope of the power-law grain size distribution. We derived a value for $q = 3.48 \pm 0.14$ for grains with sizes around 1 mm. This is consistent with the classical prediction for a collisional cascade at the steady-state. The same value cannot be explained by more recent collisional models of planetesimals in which either the velocity distribution of the large bodies or their tensile strength is a strong function of the body size.

HD 181327 is a young F5/F6V star belonging to the Beta Pictoris moving group (12 Myr). It harbors an optically thin belt of circumstellar material at 90 AU. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. We obtained far-IR observations with the Herschel/PACS instrument, and 3.2 mm observations with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS images that break the degeneracy between the disk geometry and the dust properties. We use the radiative transfer code GRaTer to compute a large grid of models, and we identify the grain models that best reproduce the Spectral Energy Distribution through a Bayesian analysis. We attempt to detect the [OI] and [CII] lines with PACS spectroscopy, providing observables to our photochemical code ProDiMo. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected in the far-IR with PACS and the disk is resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of 0.05 M\oplus (up to 1 mm). We discuss evidences that the grains are fluffy aggregates. The upper limits on the gas lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than 17 M\oplus. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

Context. To form metre-sized pre-planetesimals in protoplanetary discs, growing grains have to decouple from the gas before they are accreted onto the central star during their phase of fast radial migration and thus overcome the so-called "radial-drift barrier" (often inaccurately referred to as the "metre-size barrier"). Aims. To predict the outcome of the radial motion of dust grains in protoplanetary discs whose surface density and temperature follow power-law profiles, with exponent p and q respectively. We investigate both the Epstein and the Stokes drag regimes which govern the motion of the dust. Methods. We analytically integrate the equations of motion obtained from perturbation analysis. We compare these results with those from direct numerical integration of the equations of motion. Then, using data from observed discs, we predict the fate of dust grains in real discs. Results. When a dust grain reaches the inner regions of the disc, the acceleration due to the increase of the pressure gradient is counterbalanced by the increase of the gas drag. We find that most grains in the Epstein (resp. the Stokes) regime survive their radial migration if-p+q+1/2 \leq0 (resp. if q\leq 2/3). The majority of observed discs satisfies both-p+q+ 1/2 \leq0 and q\leq 2/3: a large fraction of both their small and large grains remain in the disc, for them the radial drift barrier does not exist.

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) mapped at 11 pc resolution in the CO(1-0) line. We identify clouds as regions of connected CO emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decomposition parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL \propto L^{-2}, suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, while apparent for the largest emission structures, breaks down when those structures are decomposed into smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of comparing two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more luminous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between giant molecular clouds and massive star formation.

Here we report on the Women in Astronomy Workshop (http://asawomeninastronomy.org/meetings/wia2011/), which was held on 13 May 2011 in Sydney, Australia. The workshop was organised by the Astronomical Society of Australia's Chapter on Women in Astronomy, to discuss some of the issues that face women in astronomy and make recommendations to help support the success of women in Australian astronomy but came to broader conclusions that have value for the whole astronomical community. The workshop consisted of four sessions, with presentations by invited speakers on demographics, leadership, varied career paths, and how institutions & individuals can help. The workshop ended with a discussion panel that summarised the day's debate and presented a list of recommendations for the Australian astronomical community (both individuals and institutions) that are provided in this report.

We are undertaking a multi-frequency Expanded Very Large Array (EVLA) survey of edge-on protoplanetary disks to probe the growth of solids in each disk, sedimentation of such material into the disk midplane, and the connection of these phenomena to the planet formation process. The projection of edge-on disk systems along our line of sight enables a study of the vertical stratification of large grains with fewer model dependencies than would be required for disks that are more face-on. Robust studies of the spatial distribution of grains up to ~1 cm in size are possible with the wavelength range and sensitivity of the EVLA. In this contribution we describe target selection and observational strategies. First results concerning the Class 0 source IRAS04368+2557 (L1527 IRS) are presented, including a study of this source's 8.46 GHz continuum variability over short and long time baselines and an indication that its protoplanetary disk may have a dearth of pebble-sized grains.

Here we describe the Compact Array Broadband Backend (CABB) and present first results obtained with the upgraded Australia Telescope Compact Array (ATCA). The 16-fold increase in observing bandwidth, from 2 x 128 MHz to 2 x 2048 MHz, high bit sampling, and addition of 16 zoom windows (each divided into a further 2048 channels) provide major improvements for all ATCA observations. The benefits of the new system are: (1) hugely increased radio continuum and polarization sensitivity as well as image fidelity, (2) substantially improved capability to search for and map emission and absorption lines over large velocity ranges, (3) simultaneous multi-line and continuum observations, (4) increased sensitivity, survey speed and dynamic range due to high-bit sampling, and (5) high velocity resolution, while maintaining full polarization output. The new CABB system encourages all observers to make use of both spectral line and continuum data to achieve their full potential. Given the dramatic increase of the ATCA capabilities in all bands (ranging from 1.1 to 105 GHz) CABB enables scientific projects that were not feasible before the upgrade, such as simultaneous observations of multiple spectral lines, on-the-fly mapping, fast follow-up of radio transients (e.g., the radio afterglow of new supernovae) and maser observations at high velocity resolution and full polarization. The first science results presented here include wide-band spectra, high dynamic-range images, and polarization measurements, highlighting the increased capability and discovery potential of the ATCA.

We utilize a chemical equilibrium code in order to study the condensation process which occurs in protoplanetary discs during the formation of the first solids. The model specifically focuses on the thermodynamic behaviour on the solid species assuming the regular solution model. For each solution, we establish the relationship between the activity of the species, the composition and the temperature using experimental data from the literature. We then apply the Gibbs free energy minimization method and study the resulting condensation sequence for a range of temperatures and pressures within a protoplanetary disc. Our results using the regular solution model show that grains condense over a large temperature range and therefore throughout a large portion of the disc. In the high temperature region (T > 1400 K) Ca-Al compounds dominate and the formation of corundum is sensitive to the pressure. The mid-temperature region is dominated by Fe(s) and silicates such as Mg2SiO4 and MgSiO3 . The chemistry of forsterite and enstatite are strictly related, and our simulations show a sequence of forsterite-enstatite-forsterite with decreasing temperature. In the low temperature regions (T < 600 K) a range of iron compounds and sulfides form. We also run simulations using the ideal solution model and see clear differences in the resulting condensation sequences with changing solution model In particular, we find that the turning point in which forsterite replaces enstatite in the low temperature region is sensitive to the solution model. Our results show that the ideal solution model is often a poor approximation to experimental data at most temperatures important in protoplanetary discs. We find some important differences in the resulting condensation sequences when using the regular solution model, and suggest that this model should provide a more realistic condensation sequence.

Context: While sub-micron- and micron-sized dust grains are generally well mixed with the gas phase in protoplanetary disks, larger grains will be partially decoupled and as a consequence have a different distribution from that of the gas. This has ramifications for predictions of the observability of protoplanetary disks, for which gas-only studies will provide an inaccurate picture. Specifically, criteria for gap opening in the presence of a planet have generally been studied for the gas phase, whereas the situation can be quite different in the dust layer once grains reach mm sizes, which is what will be observed by ALMA. Aims: We aim to investigate the formation and structure of a planetary gap in the dust layer of a protoplanetary disk with an embedded planet. Methods: We perform 3D, gas+dust SPH simulations of a protoplanetary disk with a planet on a fixed circular orbit at 40 AU to study the evolution of both the gas and dust distributions and densities in the disk. We run a series of simulations in which the planet mass and the dust grain size varies. Results: We show that the gap in the dust layer is more striking than in the gas phase and that it is deeper and wider for more massive planets as well as for larger grains. For a massive enough planet, we note that cm-sized grains remain inside the gap in corotation and that their population in the outer disk shows an asymmetric structure, a signature of disk-planet interactions even for a circular planetary orbit, which should be observable with ALMA.

Context: While sub-micron- and micron-sized dust grains are generally well mixed with the gas phase in protoplanetary disks, larger grains will be partially decoupled and as a consequence have a different distribution from that of the gas. This has ramifications for predictions of the observability of protoplanetary disks, for which gas-only studies will provide an inaccurate picture. Specifically, criteria for gap opening in the presence of a planet have generally been studied for the gas phase, whereas the situation can be quite different in the dust layer once grains reach mm sizes, which is what will be observed by ALMA. Aims: We aim to investigate the formation and structure of a planetary gap in the dust layer of a protoplanetary disk with an embedded planet. Methods: We perform 3D, gas+dust SPH simulations of a protoplanetary disk with a planet on a fixed circular orbit at 40 AU to study the evolution of both the gas and dust distributions and densities in the disk. We run a series of simulations in which the planet mass and the dust grain size varies. Results: We show that the gap in the dust layer is more striking than in the gas phase and that it is deeper and wider for more massive planets as well as for larger grains. For a massive enough planet, we note that cm-sized grains remain inside the gap in corotation and that their population in the outer disk shows an asymmetric structure, a signature of disk-planet interactions even for a circular planetary orbit, which should be observable with ALMA.

The Magellanic Mopra Assessment (MAGMA) is a high angular resolution CO mapping survey of giant molecular clouds (GMCs) in the Large and Small Magellanic Clouds using the Mopra Telescope. Here we report on the basic physical properties of 125 GMCs in the LMC that have been surveyed to date. The observed clouds exhibit scaling relations that are similar to those determined for Galactic GMCs, although LMC clouds have narrower linewidths and lower CO luminosities than Galactic clouds of a similar size. The average mass surface density of the LMC clouds is 50 Msol/pc2, approximately half that of GMCs in the inner Milky Way. We compare the properties of GMCs with and without signs of massive star formation, finding that non-star-forming GMCs have lower peak CO brightness than star-forming GMCs. We compare the properties of GMCs with estimates for local interstellar conditions: specifically, we investigate the HI column density, radiation field, stellar mass surface density and the external pressure. Very few cloud properties demonstrate a clear dependence on the environment; the exceptions are significant positive correlations between i) the HI column density and the GMC velocity dispersion, ii) the stellar mass surface density and the average peak CO brightness, and iii) the stellar mass surface density and the CO surface brightness. The molecular mass surface density of GMCs without signs of massive star formation shows no dependence on the local radiation field, which is inconsistent with the photoionization-regulated star formation theory proposed by McKee (1989). We find some evidence that the mass surface density of the MAGMA clouds increases with the interstellar pressure, as proposed by Elmegreen (1989), but the detailed predictions of this model are not fulfilled once estimates for the local radiation field, metallicity and GMC envelope mass are taken into account.

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.

Grains in disks around young stars grow from interstellar submicron sizes to planetesimals over the course of several Myr. Thermal emission of large grains or pebbles can be best observed at cm wavelengths. However, other emission mechanisms can contribute. We aim to determine the mechanisms of cm emission for 3 T Tauri stars. WW Cha and RU Lup were recently found to have grain growth at least up to mm sizes in their circumstellar disks. CS Cha has similar indications for grain growth in its circumbinary disk. The T Tauri stars WW Cha and RU Lup were monitored over several years at mm and cm wavelengths, using ATCA. The new ATCA 7 mm system was also used to observe CS Cha. WW Cha was detected on several occasions at 7 and 16 mm. We obtained one detection of WW Cha at 3 cm and upper limits only at 6 cm. The emission at 16 mm was stable over days, months and years, but the emission at 3 cm is found to be variable. RU Lup was detected at 7 mm. It was observed at 16 mm 3 times and at 3 and 6 cm 4 times and found to be variable in all 3 wavebands. CS Cha was detected at 7 mm, but the S/N was too low to resolve the gap in the circumbinary disk. The emission at 3, 7 and 16 mm for WW Cha is well explained by thermal emission from mm and cm-sized pebbles. The cm spectral index is consistent with the emission from an optically-thick ionised wind, but the high variability of the cm emission points to a non-thermal contribution. The SEDs of RU Lup and CS Cha from 1 to 7 mm are consistent with thermal emission from mm-sized grains. The variability of the longer-wavelength emission for RU Lup and the negative spectral index suggest non-thermal emission.

Aims: In order to understand the first stages of planet formation, when tiny grains aggregate to form planetesimals, one needs to simultaneously model grain growth, vertical settling and radial migration of dust in protoplanetary disks. In this study, we implement an analytical prescription for grain growth into a 3D two-phase hydrodynamics code to understand its effects on the dust distribution in disks. Methods: Following the analytic derivation of Stepinski & Valageas (1997), which assumes that grains stick perfectly upon collision, we implement a convenient and fast method of following grain growth in our 3D, two-phase (gas+dust) SPH code. We then follow the evolution of the size and spatial distribution of a dust population in a classical T Tauri star disk. Results: We find that the grains go through various stages of growth due to the complex interplay between gas drag, dust dynamics, and growth. Grains initially grow rapidly as they settle to the mid-plane, then experience a fast radial migration with little growth through the bulk of the disk, and finally pile-up in the inner disk where they grow more efficiently. This results in a bimodal distribution of grain sizes. Using this simple prescription of grain growth, we find that grains reach decimetric sizes in 10^5 years in the inner disk and survive the fast migration phase.

We investigate the behaviour of dust in protoplanetary disks under the action of gas drag using our 3D, two-fluid (gas+dust) SPH code. We present the evolution of the dust spatial distribution in global simulations of planetless disks as well as of disks containing an already formed planet. The resulting dust structures vary strongly with particle size and planetary gaps are much sharper than in the gas phase, making them easier to detect with ALMA than anticipated. We also find that there is a range of masses where a planet can open a gap in the dust layer whereas it doesn't in the gas disk. Our dust distributions are fed to the radiative transfer code MCFOST to compute synthetic images, in order to derive constraints on the settling and growth of dust grains in observed disks.

We present the first results of the treatment of grain growth in our 3D, two-fluid (gas+dust) SPH code describing protoplanetary disks. We implement a scheme able to reproduce the variation of grain sizes caused by a variety of physical processes and test it with the analytical expression of grain growth given by Stepinski & Valageas (1997) in simulations of a typical T Tauri disk around a one solar mass star. The results are in agreement with a turbulent growing process and validate the method. We are now able to simulate the grain growth process in a protoplanetary disk given by a more realistic physical description, currently under development. We discuss the implications of the combined effect of grain growth and dust vertical settling and radial migration on subsequent planetesimal formation.

We present the first results of the treatment of grain growth in our 3D, two-fluid (gas+dust) SPH code describing protoplanetary disks. We implement a scheme able to reproduce the variation of grain sizes caused by a variety of physical processes and test it with the analytical expression of grain growth given by Stepinski & Valageas (1997) in simulations of a typical T Tauri disk around a one solar mass star. The results are in agreement with a turbulent growing process and validate the method. We are now able to simulate the grain growth process in a protoplanetary disk given by a more realistic physical description, currently under development. We discuss the implications of the combined effect of grain growth and dust vertical settling and radial migration on subsequent planetesimal formation.

Aims: We investigate the behaviour of dust in protoplanetary disks under the action of gas drag in the presence of a planet. Our goal is twofold: to determine the spatial distribution of dust depending on grain size and planet mass, and therefore to provide a framework for interpretation of coming observations and future studies of planetesimal growth. Method: We numerically model the evolution of dust in a protoplanetary disk using a two-fluid (gas + dust) Smoothed Particle Hydrodynamics (SPH) code, which is non-self-gravitating and locally isothermal. The code follows the three dimensional distribution of dust in a protoplanetary disk as it interacts with the gas via aerodynamic drag. In this work, we present the evolution of a minimum mass solar nebula (MMSN) disk comprising 1% dust by mass in the presence of an embedded planet. We run a series of simulations which vary the grain size and planetary mass to see how they affect the resulting disk structure. Results: We find that gap formation is much more rapid and striking in the dust layer than in the gaseous disk and that a system with a given stellar, disk and planetary mass will have a completely different appearance depending on the grain size. For low mass planets in our MMSN disk, a gap can open in the dust disk while not in the gas disk. We also note that dust accumulates at the external edge of the planetary gap and speculate that the presence of a planet in the disk may enhance the formation of a second planet by facilitating the growth of planetesimals in this high density region.

We numerically model the evolution of dust in a protoplanetary disk using a two-phase (gas+dust) Smoothed Particle Hydrodynamics (SPH) code, which is non-self-gravitating and locally isothermal. The code follows the three dimensional distribution of dust in a protoplanetary disk as it interacts with the gas via aerodynamic drag. In this work, we present the evolution of a disk comprising 1% dust by mass in the presence of an embedded planet for two different disk configurations: a small, minimum mass solar nebular (MMSN) disk and a larger, more massive Classical T Tauri star (CTTS) disk. We then vary the grain size and planetary mass to see how they effect the resulting disk structure. We find that gap formation is much more rapid and striking in the dust layer than in the gaseous disk and that a system with a given stellar, disk and planetary mass will have a different appearance depending on the grain size and that such differences will be detectable in the millimetre domain with ALMA. For low mass planets in our MMSN models, a gap can open in the dust disk while not in the gas disk. We also note that dust accumulates at the external edge of the planetary gap and speculate that the presence of a planet in the disk may facilitate the growth of planetesimals in this high density region.

CONTEXT - Low-mass stars form with disks in which the coagulation of grains may eventually lead to the formation of planets. It is not known when and where grain growth occurs, as models that explain the observations are often degenerate. A way to break this degeneracy is to resolve the sources under study. AIMS - To find evidence for the existence of grains of millimetre sizes in disks around in T Tauri stars, implying grain growth. METHODS - The Australia Telescope Compact Array (ATCA) was used to observe 15 southern T Tauri stars, five in the constellation Lupus and ten in Chamaeleon, at 3.3 millimetre. The five Lupus sources were also observed with the Submillimeter Array (SMA) at 1.4 millimetre. Our new data are complemented with data from the literature to determine the slopes of the spectral energy distributions in the millimetre regime. RESULTS - Ten sources were detected at better than 3sigma with the ATCA, with sigma ~1-2 mJy, and all sources that were observed with the SMA were detected at better than 15sigma, with sigma ~4 mJy. Six of the sources in our sample are resolved to physical radii of ~100 AU. Assuming that the emission from such large disks is predominantly optically thin, the millimetre slope can be related directly to the opacity index. For the other sources, the opacity indices are lower limits. Four out of six resolved sources have opacity indices <~1, indicating grain growth to millimetre sizes and larger. The masses of the disks range from < 0.01 to 0.08 MSun, which is comparable to the minimum mass solar nebula. A tentative correlation is found between the millimetre slope and the strength and shape of the 10-micron silicate feature, indicating that grain growth occurs on similar (short) timescales in both the inner and outer disk.

We investigate the local correlation betwen the 1.4 GHz radio continuum and 60 micron far-infrared (FIR) emission within the Large Magellanic Cloud (LMC) on spatial scales between 0.05 and 1.5 kpc. On scales below ~1 kpc, the radio-FIR correlation is clearly better than the correlation of the cold gas tracers with either the radio or the FIR emission. For the LMC as a whole, there is a tight correlation between the radio and FIR emission on spatial scales above ~50 pc. By decomposing the radio emission into thermal and non-thermal components, however, we show that the scale on which the radio-FIR correlation breaks down is inversely proportional to the thermal fraction of the radio emission: regions that show a strong correlation to very small scales are the same regions where the thermal fraction of the radio emission is high. Contrary to previous studies of the local radio-FIR correlation in the LMC, we show that the slope of the relation between the radio and FIR emission is non-linear. In bright star-forming regions, the radio emission increases faster than linearly with respect to the FIR emission (power-law slope of ~1.2), whereas a flatter slope of ~0.6-0.9 applies more generally across the LMC. Our results are consistent with a scenario in which the UV photons and cosmic rays in the LMC have a common origin in massive star formation, but the cosmic rays are able to diffuse away from their production sites. Our results do not provide direct evidence for coupling between the magnetic field and the local gas density, but we note that synchrotron emission may not be a good tracer of the magnetic field if cosmic rays can readily escape the LMC.

We present the results of a three dimensional, locally isothermal, non-self-gravitating SPH code which models protoplanetary disks with two fluids: gas and dust. We ran simulations of a 1 Msun star surrounded by a 0.01 Msun disk comprising 99% gas and 1% dust in mass and extending from 0.5 to ~300 AU. The grain size ranges from 0.001 mm to 10 m for the low resolution (~25 000 SPH particles) simulations and from 0.1 mm to 10 cm for the high resolution (~160 000 SPH particles) simulations. Dust grains are slowed down by the sub-Keplerian gas and lose angular momentum, forcing them to migrate towards the central star and settle to the midplane. The gas drag efficiency varies according to the grain size, with the larger bodies being weakly influenced and following marginally perturbed Keplerian orbits, while smaller grains are strongly coupled to the gas. For intermediate sized grains, the drag force decouples the dust and gas, allowing the dust to preferentially migrate radially and efficiently settle to the midplane. The resulting dust distributions for each grain size will indicate, when grain growth is added, the regions when planets are likely to form.

The cold, dry and stable air above the summits of the Antarctic plateau provides the best ground-based observing conditions from optical to sub-mm wavelengths to be found on the Earth. PILOT is a proposed 2 m telescope, to be built at Dome C in Antarctica, able to exploit these conditions for conducting astronomy at optical and infrared wavelengths. While PILOT is intended as a pathfinder towards the construction of future grand-design facilities, it will also be able to undertake a range of fundamental science investigations in its own right. This paper provides the performance specifications for PILOT, including its instrumentation. It then describes the kinds of science projects that it could best conduct. These range from planetary science to the search for other solar systems, from star formation within the Galaxy to the star formation history of the Universe, and from gravitational lensing caused by exo-planets to that produced by the cosmic web of dark matter. PILOT would be particularly powerful for wide-field imaging at infrared wavelengths, achieving near-diffraction limited performance with simple tip-tilt wavefront correction. PILOT would also be capable of near-diffraction limited performance in the optical wavebands, as well be able to open new wavebands for regular ground based observation; in the mid-IR from 17 to 40 microns and in the sub-mm at 200 microns.

Infrared and submillimetre observations of nearby Vega-like stars have revealed a number of clumpy, asymmetric dust debris disks. Previous studies using semi-analytical and numerical methods have suggested planetary companions of various mass as the likely cause of most examples of disk asymmetry. In this paper, we modify an N-body symplectic gravitational integrator to include radiation terms and conduct medium-resolution parameter searches to identify likely planetary candidates in observed Vega-like systems. We also present high resolution models of Vega and epsilon Eridani, comparing our results to those of previous authors, and a new model for Fomalhaut.

We examine the historical development of astrophysical science in Antarctica from the early 20th century until today. We find three temporally overlapping eras, each having a rather distinct beginning. These are the astrogeological era of meteorite discovery, the high energy era of particle detectors, and the photon astronomy era of microwave, sub--mm and infrared telescopes, sidelined by a few niche experiments at optical wavelengths. The favourable atmospheric and geophysical conditions are briefly examined, followed by an account of the major experiments and a summary of their results.

Hierarchical clustering represents the favoured paradigm for galaxy formation throughout the Universe; due to its proximity, the Magellanic system offers one of the few opportunities for astrophysicists to decompose the full six-dimensional phase-space history of a satellite in the midst of being cannibalised by its host galaxy. The availability of improved observational data for the Magellanic Stream and parallel advances in computational power has led us to revisit the canonical tidal model describing the disruption of the Small Magellanic Cloud and the consequent formation of the Stream. We suggest improvements to the tidal model in light of these recent advances.

We present the results of three-dimensional numerical simulations that include the effects of hydrodynamical forces and gas drag upon an evolving dusty gas disk. We briefly describe a new parallel, two phase numerical code based upon the smoothed particle hydrodynamics (SPH) technique in which the gas and dust phases are represented by two distinct types of particles. We use the code to follow the dynamical evolution of a population of grains in a gaseous protoplanetary disk in order to understand the distribution of grains of different sizes within the disk. Our ``grains'' range from metre to submillimetre in size.

On examining the historical development of astrophysical science at the bottom of the world from the early 20th century until today we find three temporally overlapping eras of which each has a rather distinct beginning. These are the eras of Astrogeology, High Energy Astrophysics and Photon Astronomy. A full version of this paper will be submitted to PASA.

We have developed a new numerical technique for simulating dusty-gas flows. Our unique code incorporates gas hydrodynamics, self-gravity and dust drag to follow the dynamical evolution of a dusty-gas medium. We have incorporated several descriptions for the drag between gas and dust phases and can model flows with submillimetre, centimetre and metre size "dust". We present calculations run on the APAC* supercomputer following the evolution of the dust distribution in the pre-solar nebula. *Australian Partnership for Advanced Computing

We are in a privileged location in the Universe which allows us to observe galactic interactions from close range -- the merger of our two nearest dwarf satellite galaxies, the LMC and SMC. It is important to understand the local merger process before we can have confidence in understanding mergers at high redshift. We present high resolution Nbody+SPH simulations of the disruption of the LMC and SMC and the formation of the Magellanic Stream, and discuss the implications for galaxy formation and evolution.

Motivated by the apparent order-of-magnitude discrepancy between the observed number of Local Group satellite galaxies, and that predicted by Lambda-CDM hierarchical clustering cosmologies, we explore an alternate suggestion - perhaps the missing satellites are not actually ``missing'', but are instead ``in disguise''. The disguise we consider here is that of the classical HI High-Velocity Clouds. Is it possible that what have been thought of traditionally as a ``Galactic'' phenomenon, are actually the building blocks of the Local Group? We discuss the strengths and weaknesses of this hypothesis, and highlight avenues of future research which may provide an unequivocal resolution to this contentious issue.

We model a massless viscous disk using Smoothed Particle Hydrodynamics (SPH) and note that it evolves according to the Lynden-Bell \& Pringle theory (1974) until a non-axisymmetric instability develops at the inner edge of the disk. This instability may have the same origin as the instability of initially axisymmetric viscous disks discussed by Lyubarskij et al. (1994). To clarify the evolution we evolved single and double rings of particles. It is actually inconsistent with the SPH scheme to set up a single ring as an initial condition because SPH assumes a smoothed initial state. As would be expected from an SPH simulation, the ring rapidly breaks up into a band. We analyse the stability of the ring and show that the predictions are confirmed by the simulation.