Source author record

Cathie Clarke

Cathie Clarke 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

15works
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

15 published item(s)

preprint2022arXiv

Accretion rates in hierarchical triple systems with discs

Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical configurations. Despite its astrophysical relevance, differential accretion in hierarchical systems remains to be understood. In this work, we investigate whether the accretion trends expected in binaries are valid for higher order multiples. We performed a set of 3D Smoothed Particle Hydrodynamics simulations of binaries and of hierarchical triples (HTs) embedded in an accretion disc, with the code Phantom. We identify for the first time accretion trends in HTs and their deviations compared to binaries. These deviations, due to the interaction of the small binary with the infalling material from the circum-triple disc, can be described with a semi-analytical prescription. Generally, the smaller binary of a HT accretes more mass than a single star of the same mass as the smaller binary. We found that in a HT, if the small binary is heavier than the third body, the standard differential accretion scenario (whereby the secondary accretes more of the mass) is hampered. Reciprocally, if the small binary is lighter than the third body, the standard differential accretion scenario is enhanced. The peculiar differential accretion mechanism we find in HTs is expected to affect their mass ratio distribution.

preprint2020arXiv

The efficiency of dust trapping in ringed proto-planetary discs

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

preprint2019arXiv

The distribution of relative proper motions of wide binaries in GAIA DR2: MOND or multiplicity?

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

preprint2016arXiv

Efficiency of Planetesimal Ablation in Giant Planetary Envelopes

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

preprint2014arXiv

Dancing with the Stars: Formation of the Fomalhaut triple system and its effect on the debris disks

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

preprint2013arXiv

The evolution of planetesimal swarms in self-gravitating protoplanetary discs

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

preprint2012arXiv

Two populations of transition discs?

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

preprint2010arXiv

Stellar collisions in accreting protoclusters: a Monte Carlo dynamical study

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

preprint2009arXiv

Metallicity, planet formation, and disc lifetimes

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

preprint2009arXiv

The effects of opacity on gravitational stability in protoplanetary discs

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

preprint2009arXiv

The physics and modes of star cluster formation: simulations

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

preprint2007arXiv

A new algorithm for modelling photoionising radiation in smoothed particle hydrodynamics

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

preprint2000arXiv

Destruction of protoplanetary disks in the Orion Nebula Cluster

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

preprint1995arXiv

Low Mass Companions to T~Tauri Stars: a Mechanism for Rapid-rise FU~Orionis Outbursts

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

preprint1995arXiv

Satellites in Discs: Regulating the Accretion Luminosity

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