Source author record

Sarah T. Maddison

Sarah T. Maddison appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Constraints on Planetesimal Collision Models in Debris Disks

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.

preprint2016arXiv

On dust entrainment in photoevaporative winds

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.

preprint2016arXiv

On the maximum grain size entrained by photoevaporative winds

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

Growing dust grains in protoplanetary discs - III. Vertical settling

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.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Grain growth signatures in the protoplanetary discs of Chamaeleon and Lupus

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

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

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.

preprint2010arXiv

Grain growth across protoplanetary discs: 10-micron silicate feature versus millimetre slope

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

preprint2010arXiv

Planet gaps in the dust layer of 3D protoplanetary disks: I. Hydrodynamical simulations of T Tauri disks

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

preprint2005arXiv

Dust distribution in protoplanetary disks - Vertical settling and radial migration

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