Showing 1–50 of 73 results for author: Maddison, S

The Origin of the Doppler-flip in HD 100546: a large scale spiral arm generated by an inner binary companion

Authors: Brodie J. Norfolk, Christophe Pinte, Josh Calcino, Iain Hammond, Nienke van der Marel, Daniel J. Price, Sarah T. Maddison, Valentin Christiaens, Jean-Francois Gonzalez, Dori Blakely, Giovanni Rosotti, Christian Ginski

Abstract: Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner sp… ▽ More Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner spiral and the previously reported CO emission Doppler-flip, which has been interpreted as the signature of an embedded protoplanet. Comparisons with hydrodynamical models indicate that this Doppler-flip is instead the kinematic counterpart of the spiral, which is likely generated by an inner companion inside the disk cavity. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL, 14 pages, 10 figures

Four new PLanetesimals Around TYpical and Pre-main seqUence Stars (PLATYPUS) Debris Discs at 8.8mm

Authors: Brodie J. Norfolk, Sarah T. Maddison, Jonathan P. Marshall, Grant M. Kennedy, Gaspard Duchêne, David J. Wilner, Christophe Pinte, Attila Moór, Brenda Matthews, Péter Ábrahám, Ágnes Kóspál, Nienke van der Marel

Abstract: Millimetre continuum observations of debris discs can provide insights into the physical and dynamical properties of the unseen planetesimals that these discs host. The material properties and collisional models of planetesimals leave their signature on the grain size distribution, which can be traced through the millimetre spectral index. We present 8.8 mm observations of the debris discs HD 4837… ▽ More Millimetre continuum observations of debris discs can provide insights into the physical and dynamical properties of the unseen planetesimals that these discs host. The material properties and collisional models of planetesimals leave their signature on the grain size distribution, which can be traced through the millimetre spectral index. We present 8.8 mm observations of the debris discs HD 48370, CPD 72 2713, HD 131488, and HD 32297 using the Australian Telescope Compact Array (ATCA) as part of the PLanetesimals Around TYpicalPre-main seqUence Stars (PLATYPUS) survey. We detect all four targets with a characteristic beam size of 5 arcseconds and derive a grain size distribution parameter that is consistent with collisional cascade models and theoretical predictions for parent planetesimal bodies where binding is dominated by self-gravity. We combine our sample with 19 other millimetre-wavelength detected debris discs from the literature and calculate a weighted mean grain size power law index which is close to analytical predictions for a classical steady state collisional cascade model. We suggest the possibility of two distributions of q in our debris disc sample; a broad distribution (where q is approximately 3.2 to 3.7) for "typical" debris discs (gas-poor/non-detection), and a narrow distribution (where q is less than 3.2) for bright gas-rich discs. Or alternatively, we suggest that there exists an observational bias between the grain size distribution parameter and absolute flux which may be attributed to the detection rates of faint debris discs at cm wavelengths. △ Less

Submitted 1 July, 2021; v1 submitted 21 June, 2021; originally announced June 2021.

Comments: 10 pages, 4 figures, accepted

Dust Traps and the Formation of Cavities in Transition Discs: A millimetre to sub-millimetre comparison survey

Authors: Brodie J. Norfolk, Sarah T. Maddison, Christophe Pinte, Nienke van der Marel, Richard A. Booth, Logan Francis, Jean-François Gonzalez, François Ménard, Chris M. Wright, Gerrit van der Plas, Himanshi Garg

Abstract: The origin of the inner dust cavities observed in transition discs remains unknown. The segregation of dust and size of the cavity is expected to vary depending on which clearing mechanism dominates grain evolution. We present the results from the Discs Down Under program, an 8.8 mm continuum Australia Telescope Compact Array (ATCA) survey targeting 15 transition discs with large (> 20 au) cavitie… ▽ More The origin of the inner dust cavities observed in transition discs remains unknown. The segregation of dust and size of the cavity is expected to vary depending on which clearing mechanism dominates grain evolution. We present the results from the Discs Down Under program, an 8.8 mm continuum Australia Telescope Compact Array (ATCA) survey targeting 15 transition discs with large (> 20 au) cavities, and compare the resulting dust emission to Atacama Large millimetre/sub-millimetre Array (ALMA) observations. Our ATCA observations resolve the inner cavity for 8 of the 14 detected discs. We fit the visibilities and reconstruct 1D radial brightness models for 10 sources with a S/N > 5sigma. We find that, for sources with a resolved cavity in both wavebands, the 8.8 mm and sub-mm brightness distributions peak at the same radius from the star. We suggest that a similar cavity size for 8.8 mm and sub-mm dust grains is due to a dust trap induced by the presence of a companion. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 18 pages, 11 figures

High-resolution survey for planetary companions to young stars in the Taurus Molecular Cloud

Authors: A. L. Wallace, J. Kammerer, M. J. Ireland, C. Federrath, A. L. Kraus, S. T. Maddison, A. C. Rizzuto, E. K. Birchall, F. Martinache

Abstract: Direct imaging in the infrared at the diffraction limit of large telescopes is a unique probe of the properties of young planetary systems. We survey 55 single class I and class II stars in Taurus in the L' filter using natural and laser guide star adaptive optics and the near-infrared camera (NIRC2) of the Keck II telescope, in order to search for planetary mass companions. We use both reference… ▽ More Direct imaging in the infrared at the diffraction limit of large telescopes is a unique probe of the properties of young planetary systems. We survey 55 single class I and class II stars in Taurus in the L' filter using natural and laser guide star adaptive optics and the near-infrared camera (NIRC2) of the Keck II telescope, in order to search for planetary mass companions. We use both reference star differential imaging and kernel phase techniques, achieving typical 5-sigma contrasts of ~6 magnitudes at separations of 0.2" and ~8 magnitudes beyond 0.5". Although we do not detect any new faint companions, we constrain the frequency of wide separation massive planets, such as HR 8799 analogues. We find that, assuming hot-start models and a planet distribution with power-law mass and semi-major axis indices of -0.5 and -1, respectively, less than 20% of our target stars host planets with masses >2 MJ at separations >10 AUえーゆー. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

Infrared Variability due to Magnetic Pressure Driven Jets, Dust Ejection and Quasi-Puffed-Up Inner Rims

Authors: Kurt Liffman, Geoffrey Bryan, Mark Hutchison, Sarah T. Maddison

Abstract: 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 "… ▽ More 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~$μみゅー$m radiation. The subsequent shadow reduces the amount of light falling on the outer disc and decreases the 8-- 40~$μみゅー$m radiation. A decrease in the accretion rate reverses this scenario, thereby producing the observed "see-saw" infrared variability. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 18 pages, 21 figures

VLA cm-wave survey of young stellar objects in the Oph A cluster: constraining extreme UV- and X-ray-driven disk photo-evaporation -- A pathfinder for Square Kilometre Array studies

Authors: A. Coutens, H. B. Liu, I. Jiménez-Serra, T. L. Bourke, J. Forbrich, M. Hoare, L. Loinard, L. Testi, M. Audard, P. Caselli, A. Chacón-Tanarro, C. Codella, J. Di Francesco, F. Fontani, M. Hogerheijde, A. Johansen, D. Johnstone, S. Maddison, O. Panić, L. M. Pérez, L. Podio, A. Punanova, J. M. C. Rawlings, D. Semenov, M. Tazzari , et al. (6 additional authors not shown)

Abstract: 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 Larg… ▽ More 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 $θしーた_{maj}$ $\times$ $θしーた_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $μみゅー$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. △ Less

Submitted 8 September, 2019; originally announced September 2019.

Comments: Accepted in A&A

Journal ref: A&A 631, A58 (2019)

The fulcrum wavelength of young stellar objects -- the case of LRLL 31

Authors: Geoffrey R. Bryan, Sarah T. Maddison, Kurt Liffman

Abstract: 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 ($λらむだ_{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 t… ▽ More 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 ($λらむだ_{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 $λらむだ_{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 $λらむだ_{f}$ and its position relative to the silicate feature near $10 μみゅー$m. While $λらむだ_{f}$ is a function of each of these parameters, it is found to be most strongly dependent on the vertical density exponent $βべーた$. All other factors being held constant, only for flatter discs ($βべーた< 1.2$) did we find a $λらむだ_{f}$ beyond the silicate feature. △ Less

Submitted 23 August, 2019; originally announced August 2019.

Comments: 12 pages, 11 figures

Tiny Grains Shining Bright in the Gaps of Herbig Ae Transitional Discs

Authors: Eloise K. Birchall, Michael J. Ireland, Christoph Federrath, John D. Monnier, Stefan Kraus, Matthew Willson, Adam L. Kraus, Aaron Rizzuto, Matthew T. Agnew, Sarah T. Maddison

Abstract: 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 forw… ▽ More 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. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: Accepted at MNRAS. 22 pages, 12 figures

Planet formation: The case for large efforts on the computational side

Authors: Wladimir Lyra, Thomas Haworth, Bertram Bitsch, Simon Casassus, Nicolás Cuello, Thayne Currie, Andras Gáspár, Hannah Jang-Condell, Hubert Klahr, Nathan Leigh, Giuseppe Lodato, Mordecai-Mark Mac Low, Sarah Maddison, George Mamatsashvili, Colin McNally, Andrea Isella, Sebastián Pérez, Luca Ricci, Debanjan Sengupta, Dimitris Stamatellos, Judit Szulágyi, Richard Teague, Neal Turner, Orkan Umurhan, Jacob White , et al. (32 additional authors not shown)

Abstract: 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 a… ▽ More 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. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: White paper submitted to the Astro2020 decadal survey

Predicting multiple planet stability and habitable zone companions in the TESS era

Authors: Matthew T. Agnew, Sarah T. Maddison, Jonathan Horner, Stephen R. Kane

Abstract: 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… ▽ More 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. △ Less

Submitted 31 January, 2019; originally announced January 2019.

Comments: 24 pages, 16 figures

Prospecting for exo-Earths in multiple planet systems with a gas giant

Authors: Matthew T. Agnew, Sarah T. Maddison, Jonathan Horner

Abstract: 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 t… ▽ More 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. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 20 pages, 12 figures, 6 tables

Properties of the single Jovian planet population and the pursuit of Solar system analogues

Authors: Matthew T. Agnew, Sarah T. Maddison, Jonathan Horner

Abstract: 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 un… ▽ More 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. △ Less

Submitted 18 April, 2018; originally announced April 2018.

Comments: 14 pages, 10 figures

Stable habitable zones of single Jovian planet systems

Authors: Matthew T. Agnew, Sarah T. Maddison, Elodie Thilliez, Jonathan Horner

Abstract: 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 d… ▽ More 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. △ Less

Submitted 19 June, 2017; originally announced June 2017.

Comments: 15 pages, 12 figures, Accepted for publication by MNRAS

New constraints on the millimetre emission of six debris disks

Authors: Jonathan P. Marshall, S. T. Maddison, E. Thilliez, B. C. Matthews, D. J. Wilner, J. S. Greaves, W. S. Holland

Abstract: 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… ▽ More 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-σしぐま) 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. △ Less

Submitted 27 March, 2017; originally announced March 2017.

Comments: 7 pages, 5 tables, 2 figures, accepted for publication in MNRAS

Radio Monitoring of Protoplanetary Discs

Authors: C. Ubach, S. T. Maddison, C. M. Wright, D. J. Wilner, D. J. P. Lommen, B. Koribalski

Abstract: 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 emi… ▽ More 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. △ Less

Submitted 8 January, 2017; originally announced January 2017.

Comments: Accepted 2017 January 4, Monthly Notices of the Royal Astronomical Society

Self-induced dust traps: overcoming planet formation barriers

Authors: Jean-François Gonzalez, Guillaume Laibe, Sarah T. Maddison

Abstract: 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… ▽ More 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. △ Less

Submitted 31 August, 2017; v1 submitted 4 January, 2017; originally announced January 2017.

Comments: 15 pages, 16 figures. Added a 1-page erratum correcting an inconsequential error in Appendix C. A higher resolution version of the paper, as well as the online videos, are available at https://cral.univ-lyon1.fr/labo/perso/jean-francois.gonzalez/Papers/SIDT/

Journal ref: MNRAS 467, 1984 (2017)

Numerical search for a potential planet sculpting the young disc of HD 115600

Authors: E. Thilliez, S. T. Maddison

Abstract: 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 eccentr… ▽ More 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. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

A cavity and further radial substructures in the disk around HD~97048

Authors: G. van der Plas, C. M. Wright, F. Ménard, S. Casassus, H. Canovas, C. Pinte, S. T. Maddison, K. Maaskant, H. Avenhaus, L. Cieza, S. Perez, C. Ubach

Abstract: 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. Meth… ▽ More 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. △ Less

Submitted 8 September, 2016; originally announced September 2016.

Comments: 18 pages, 16 figures. Accepted to Astronomy & Astrophysics September 3rd 2016

Journal ref: A&A 597, A32 (2017)

On the maximum grain size entrained by photoevaporative winds

Authors: Mark A. Hutchison, Guillaume Laibe, Sarah T. Maddison

Abstract: 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 rep… ▽ More 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. △ Less

Submitted 29 August, 2016; originally announced August 2016.

Comments: 11 pages, 9 figures, 1 table, 1 animation. Accepted for publication in MNRAS

Constraints on Planetesimal Collision Models in Debris Disks

Authors: Meredith A. MacGregor, David J. Wilner, Claire Chandler, Luca Ricci, Sarah T. Maddison, Steven R. Cranmer, Sean M. Andrews, A. Meredith Hughes, Amy Steele

Abstract: 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… ▽ More 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$ $μみゅー$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. △ Less

Submitted 17 March, 2016; originally announced March 2016.

Comments: 34 pages, 5 figures

On dust entrainment in photoevaporative winds

Authors: Mark A. Hutchison, Daniel J. Price, Guillaume Laibe, Sarah T. Maddison

Abstract: 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… ▽ More 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. △ Less

Submitted 10 May, 2016; v1 submitted 9 March, 2016; originally announced March 2016.

Comments: 20 pages, 13 figures, Accepted for publication in MNRAS

Journal ref: MNRAS 2016 461: 742-759

Numerical predictions for planets in the debris discs of HD 202628 and HD 207129

Authors: E. Thilliez, S. T. Maddison

Abstract: 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… ▽ More 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. △ Less

Submitted 10 January, 2016; originally announced January 2016.

Comments: 15 Pages, 8 Figures, Accepted in MNRAS

2D condensation model for the inner Solar Nebula: an enstatite-rich environment

Authors: Francesco C. Pignatale, Kurt Liffman, Sarah T. Maddison, Geoffrey Brooks

Abstract: 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… ▽ More 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. △ Less

Submitted 7 January, 2016; originally announced January 2016.

Comments: Astroph version. Accepted in MNRAS on 2015 December 23. Received 2015 December 16; in original form 2015 May 04

An ATCA survey of debris disks at 7 millimeters

Authors: L. Ricci, S. T. Maddison, D. Wilner, M. A. MacGregor, C. Ubach, J. M. Carpenter, L. Testi

Abstract: We present ATCA continuum observations at a wavelength of 6.8 mm of five debris disks: $βべーた$ 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 s… ▽ More We present ATCA continuum observations at a wavelength of 6.8 mm of five debris disks: $βべーた$ 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. △ Less

Submitted 12 October, 2015; originally announced October 2015.

Comments: 7 pages, 3 figures, accepted for publication on ApJ

The impact of initial conditions in N-body simulations of debris discs

Authors: E. Thilliez, S. T. Maddison

Abstract: 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 e… ▽ More 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. △ Less

Submitted 29 September, 2015; originally announced September 2015.

Comments: 17 Pages, 13 Figures, Accepted in PASA

ALMA images of discs: are all gaps carved by planets?

Authors: Jean-François Gonzalez, Guillaume Laibe, Sarah T. Maddison, Christophe Pinte, François Ménard

Abstract: 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 ga… ▽ More 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. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Comments: Accepted for publication in MNRAS Letters. 5 pages, 6 figures. MNRAS, Oxford University press, 2015, in press

The Epsilon Eridani System Resolved by Millimeter Interferometry

Authors: Meredith A. MacGregor, David J. Wilner, Sean M. Andrews, Jean-Francois Lestrade, Sarah Maddison

Abstract: 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 oute… ▽ More 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 $Δでるた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. △ Less

Submitted 6 July, 2015; originally announced July 2015.

Comments: 27 pages, 6 figures

Resolving structure of the disk around HD100546 at 7 mm with ATCA

Authors: Christopher Wright, Sarah Maddison, David Wilner, Michael Burton, Dave Lommen, Ewine van Dishoeck, Paola Pinilla, Tyler Bourke, Francois Menard, Catherine Walsh

Abstract: 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… ▽ More 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$λらむだ$, 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えーゆー. △ Less

Submitted 3 June, 2015; originally announced June 2015.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

The accumulation and trapping of grains at planet gaps: effects of grain growth and fragmentation

Authors: J. -F Gonzalez, G Laibe, S. T. Maddison, C Pinte, F Ménard

Abstract: 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… ▽ More 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. △ Less

Submitted 1 June, 2015; originally announced June 2015.

Comments: Accepted for publication in Planetary and Space Science. 13 pages, 4 figures

A compact concentration of large grains in the HD142527 protoplanetary dust trap

Authors: Simon Casassus, Chris Wright, Sebastian Marino, Sarah T. Maddison, Al Wootten, Pablo Roman, Sebastian Perez, Paola Pinilla, Mark Wyatt, Victor Moral, Francois Menard, Valentin Christiaens, Lucas Cieza, Gerrit van der Plas

Abstract: 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 seg… ▽ More 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. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: Submitted to ApJ

Journal ref: 2015, ApJ, 812, 126

Compact dust concentration in the MWC 758 protoplanetary disk

Authors: S. Marino, S. Casassus, S. Perez, W. Lyra, P. E. Roman, H. Avenhaus, C. M. Wright, S. T. Maddison

Abstract: 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 K… ▽ More 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. △ Less

Submitted 2 October, 2015; v1 submitted 25 May, 2015; originally announced May 2015.

Comments: Accepted for publication in ApJ

A dynamical test for terrestrial planets in the habitable zone of HD 204313

Authors: E. Thilliez, L. Jouvin, S. T. Maddison, J. Horner

Abstract: 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 dyn… ▽ More 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. △ Less

Submitted 11 February, 2014; originally announced February 2014.

Comments: 13 pages, 6 figures. Accepted for publication in the peer-reviewed proceedings of the 13th annual Australian Space Science Conference

Tidal stability of giant molecular clouds in the Large Magellanic Cloud

Authors: E. Thilliez, S. T. Maddison, A. Hughes, T. Wong

Abstract: 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… ▽ More 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. △ Less

Submitted 24 October, 2013; originally announced October 2013.

Comments: 18 pages, 10 Figures, Accepted in PASA

Growing dust grains in protoplanetary discs - III. Vertical settling

Authors: Guillaume Laibe, Jean-François Gonzalez, Sarah T. Maddison, Elisabeth Crespe

Abstract: 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 t… ▽ More 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. △ Less

Submitted 6 November, 2013; v1 submitted 9 October, 2013; originally announced October 2013.

Comments: 8 pages, 9 figures. Accepted for publication in MNRAS. v2: typos corrected and minor changes

Growing dust grains in protoplanetary discs - I. Radial drift with toy growth models

Authors: Guillaume Laibe, Jean-François Gonzalez, Sarah T. Maddison

Abstract: 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 intr… ▽ More 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. △ Less

Submitted 6 November, 2013; v1 submitted 9 October, 2013; originally announced October 2013.

Comments: 12 pages, 9 figures. Accepted for publication in MNRAS. v2: typos corrected

Planet gaps in the dust layer of 3D proto-planetary disks: Observability with ALMA

Authors: Jean-François Gonzalez, Christophe Pinte, Sarah T. Maddison, François Ménard

Abstract: 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… ▽ More 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. △ Less

Submitted 5 July, 2013; originally announced July 2013.

Comments: 2 pages, 2 figures, to appear in the Proceedings of IAU Symp. 299: Exploring the Formation and Evolution of Planetary Systems (Victoria, Canada)

Can Galactic chemical evolution explain the oxygen isotopic variations in the Solar System?

Authors: Maria Lugaro, Kurt Liffman, Trevor R. Ireland, Sarah T. Maddison

Abstract: 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 relati… ▽ More 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. △ Less

Submitted 11 September, 2012; originally announced September 2012.

Comments: 28 pages, 3 figures, accepted for publication on The Astrophysical Journal

Short-lived radioactivity in the early Solar System: the Super-AGB star hypothesis

Authors: Maria Lugaro, Carolyn L. Doherty, Amanda. I. Karakas, Sarah T. Maddison, Kurt Liffman, D. A. García-Hernández, Lionel Siess, John C. Lattanzio

Abstract: 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 pred… ▽ More 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. △ Less

Submitted 28 August, 2012; originally announced August 2012.

Comments: 23 pages, 5 figures, 6 tables, accepted for publication on Meteoritics and Planetary Science

Planet gaps in the dust layer of 3D protoplanetary disks. II. Observability with ALMA

Authors: J. -F. Gonzalez, C. Pinte, S. T. Maddison, F. Ménard, L. Fouchet

Abstract: [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… ▽ More [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. △ Less

Submitted 27 August, 2012; originally announced August 2012.

Comments: 15 pages, 21 figures, accepted by Astronomy & Astrophysics, a higher resolution version of the paper is available at http://www-obs.univ-lyon1.fr/labo/perso/jean-francois.gonzalez/Papers/Gaps_ALMA.pdf

Refractory Metal Nuggets -- Formation of the First Condensates in the Solar Nebula

Authors: Kurt Liffman, Francesco C. Pignatale, Sarah Maddison, Geoffrey Brooks

Abstract: 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 dependan… ▽ More 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. △ Less

Submitted 23 July, 2012; originally announced July 2012.

Comments: 51 pages, 22 figures

Grain growth signatures in the protoplanetary discs of Chamaeleon and Lupus

Authors: Catarina Ubach, Sarah T. Maddison, Christopher M. Wright, David J. Wilner, Dave J. P. Lommen, Baerbel Koribalski

Abstract: 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 sour… ▽ More 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 μみゅーm silicate feature) and millimetre (1-3 mm spectral slope, αあるふぁ). 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 αあるふぁ, 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. △ Less

Submitted 1 July, 2012; originally announced July 2012.

Comments: 31 pages, 13 figures, 11 tables, accepted for publication in MNRAS

Journal ref: 2012, MNRAS, 425, 3137-3161

The Milky Way Project First Data Release: A Bubblier Galactic Disk

Authors: R. J. Simpson, M. S. Povich, S. Kendrew, C. J. Lintott, E. Bressert, K. Arvidsson, C. Cyganowski, S. Maddison, K. Schawinski, R. Sherman, A. M. Smith, G. Wolf-Chase

Abstract: 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… ▽ More 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. △ Less

Submitted 30 January, 2012; originally announced January 2012.

Comments: 19 pages, 20 figures, submitted to MNRAS

Fomalhaut debris disk emission at 7 millimeters: constraints on the collisional models of planetesimals

Authors: L. Ricci, L. Testi, S. T. Maddison, D. J. Wilner

Abstract: 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-millimete… ▽ More 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. △ Less

Submitted 16 January, 2012; originally announced January 2012.

Comments: 4 pages, 2 figures; A\&A Letters, in press

An icy Kuiper-Belt around the young solar-type star HD 181327

Authors: J. Lebreton, J. -C. Augereau, W. -F. Thi, A. Roberge, J. Donaldson, G. Schneider, S. T. Maddison, F. Ménard, P. Riviere-Marichalar, G. S. Mathews, I. Kamp, C. Pinte, W. R. F. Dent, D. Barrado, G. Duchêne, J. -F. Gonzalez, C. A. Grady, G. Meeus, E. Pantin, J. P. Williams, P. Woitke

Abstract: 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 con… ▽ More 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. △ Less

Submitted 14 December, 2011; originally announced December 2011.

Comments: 17 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics (2011 december 7)

Revisiting the "radial-drift barrier" of planet formation and its relevance in observed protoplanetary discs

Authors: Guillaume Laibe, Jean-François Gonzalez, Sarah T. Maddison

Abstract: 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 protopl… ▽ More 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. △ Less

Submitted 13 November, 2011; originally announced November 2011.

Comments: Accepted for publication in A&A

The Magellanic Mopra Assessment (MAGMA). I. The Molecular Cloud Population of the Large Magellanic Cloud

Authors: Tony Wong, Annie Hughes, Jürgen Ott, Erik Muller, Jorge L. Pineda, Jean-Philippe Bernard, You-Hua Chu, Yasuo Fukui, Robert A. Gruendl, Christian Henkel, Akiko Kawamura, Ulrich Klein, Leslie W. Looney, Sarah Maddison, Yoji Mizuno, Deborah Paradis, Jonathan Seale, Daniel E. Welty

Abstract: 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 steep… ▽ More 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. △ Less

Submitted 29 August, 2011; originally announced August 2011.

Comments: Accepted by ApJS; 22 pages in emulateapj format; full-resolution version and data tables available at http://mmwave.astro.illinois.edu/magma/

Women in Astronomy Workshop Report

Authors: Sarah Brough, Amanda E. Bauer, Kate Brooks, Andrew Hopkins, Sarah Maddison

Abstract: 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… ▽ More 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. △ Less

Submitted 29 June, 2011; originally announced June 2011.

Comments: 9 pages

Microwave Observations of Edge-on Protoplanetary Disks: Program Overview and First Results

Authors: Carl Melis, G. Duchene, Laura Chomiuk, Patrick Palmer, M. D. Perrin, S. T. Maddison, F. Menard, K. Stapelfeldt, C. Pinte, G. Duvert

Abstract: 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 gr… ▽ More 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. △ Less

Submitted 24 June, 2011; v1 submitted 22 June, 2011; originally announced June 2011.

Comments: 17 pages, 2 tables, 3 figures. Accepted for ApJ Letters EVLA special issue. Updated affiliations

The Australia Telescope Compact Array Broadband Backend (CABB)

Authors: Warwick E. Wilson, R. H. Ferris, P. Axtens, A. Brown, E. Davis, G. Hampson, M. Leach, P. Roberts, S. Saunders, B. S. Koribalski, J. L. Caswell, E. Lenc, J. Stevens, M. A. Voronkov, M. H. Wieringa, K. Brooks, P. G. Edwards, R. D. Ekers, B. Emonts, L. Hindson, S. Johnston, S. T. Maddison, E. K. Mahony, S. S. Malu, M. Massardi , et al. (8 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 18 May, 2011; originally announced May 2011.

Comments: 27 pages, 34 figures, MNRAS, accepted

The effect of the regular solution model in the condensation of protoplanetary dust

Authors: Francesco C. Pignatale, Sarah T. Maddison, Vianney Taquet, Geoffrey Brooks, Kurt Liffman

Abstract: 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 temperatu… ▽ More 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. △ Less

Submitted 28 February, 2011; originally announced February 2011.

Comments: MNRAS: Accepted 2011 February 16. Received 2011 February 14; in original form 2010 July 22