BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Cathie J. Clarke

Cathie J. Clarke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.GAastro-ph.IMastro-ph.HE

Catalog footprint

What is connected

44works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

44 published item(s)

preprint2022arXiv

An APEX search for carbon emission from NGC 1977 proplyds

We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3$σ$ detections of the CI line were made despite reaching sensitivities deeper than the anticipated requirement for detection from prior APEX CI observations of nearby discs and models of external photoevaporation of quite massive discs. By comparing both the proplyd mass loss rates and CI flux constraints with a large grid of external photoevaporation simulations, we determine that the non-detections are in fact fully consistent with the models if the proplyd discs are very low mass. Deeper observations in CI and probes of the disc mass with other tracers (e.g. in the continuum and CO) can test this. If such a test finds higher masses, this would imply carbon depletion in the outer disc, as has been proposed for other discs with surprisingly low CI fluxes, though more massive discs would also be incompatible with models that can explain the observed mass loss rates and CI non-detections. The expected remaining lifetimes of the proplyds are estimated to be similar to those of proplyds in the ONC at 0.1Myr. Rapid destruction of discs is therefore also a feature of common, intermediate UV environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Forming short period sub-stellar companions in 47 Tucanae -- II. Analytic expressions for the orbital evolution of planets in dense environments

Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularising tides in the planet. We present an analytic framework for quantifying the role of dynamical encounters in high density environments during HEM. We show that encounters can enhance or suppress HEM, depending on the local stellar density and the initial semi-major axis $a_0$. For moderate densities, external perturbations can excite large eccentricities that allow a planet to circularise over the stellar lifetime. At extremely high densities, these perturbations can instead result in tidal disruption of the planet, thus yielding no HJ. This may explain the apparent excess of HJs in M67 compared with their local field star abundance versus their apparent deficit in 47 Tuc. Applying our analytic framework, we demonstrate that for an initial massive planet population similar to the field, the expected HJ occurrence rate in 47 Tuc is $f_\mathrm{HJ}=2.2\times 10^{-3}$, which remains consistent with present constraints. Future large (sample sizes $\gtrsim 10^5$) or sensitive transit surveys of stars in globular clusters are required to refute the hypothesis that the initial planet population is similar to the solar neighbourhood average. Non-detection in such surveys would have broad consequences for planet formation theory, implying planet formation rates in globular clusters must be suppressed across a wide range of $a_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modelling the secular evolution of proto-planetary disc dust sizes -- A comparison between the viscous and magnetic wind case

For many years proto-planetary discs have been thought to evolve viscously: angular momentum redistribution leads to accretion and outward disc spreading. Recently, the hypothesis that accretion is due, instead, to angular momentum removal by magnetic winds gained new popularity: no disc spreading is expected in this case. In this paper, we run several one-dimensional gas and dust simulations to make predictions on the time-evolution of disc sizes \textit{in the dust} and to assess whether they can be used to understand how discs evolve. We show that viscous and magnetic wind models have very different dust disc radii. In particular, MHD wind models are compact and their sizes either remain constant or decrease with time. On the contrary, discs become larger with time in the viscous case (when $α\gtrsim10^{-3}$). Although current observations lack enough sensitivity to discriminate between these two scenarios, higher-sensitivity surveys could be fruitful to this goal on a $1\,{\rm to}\,10\,{\rm Myr}$ age range. When compared with the available ALMA Band~7 data, both viscous and magnetic wind models are compatible with the observationally-inferred dust radii in Lupus, Chamaeleon~I and Upper Sco. Furthermore, in the drift-dominated regime, the size-luminosity correlation is reproduced in Lupus, both in Band~7 and 3, while in Upper Sco a different slope than in the data is predicted. Sub-structures (potentially undetected) can explain several outliers with large observed sizes. Higher-angular-resolution observations will be helpful to test our predictions in the case of more compact discs, expected in both frameworks, particularly at the age of Upper Sco.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stellar multiplicity affects the correlation between proto-planetary disc masses and accretion rates: binaries explain high-accretors in Upper Sco

In recent years a correlation between mass accretion rates onto new-born stars and their proto-planetary disc masses was detected in nearby young star-forming regions. Although such a correlation can be interpreted as due to viscous-diffusion processes in the disc, highly-accreting sources with low disc masses in more evolved regions remain puzzling. In this paper, we hypothesise that the presence of a stellar companion truncating the disc can explain these outliers. Firstly, we searched the literature for information on stellar multiplicity in Lupus, Chamaeleon~I and Upper Sco, finding that roughly 20 per cent of the discs involved in the correlation are in binaries or higher-order multiple stellar systems. We prove with high statistical significance that at any disc mass these sources have systematically higher accretion rates than those in single-stars, with the bulk of the binary population being clustered around $M_\mathrm{disc}/\dot{M}_\mathrm{acc}\approx0.1\,\mathrm{Myr}$. We then run coupled gas and dust one-dimensional evolutionary models of tidally truncated discs to be compared with the data. We find that these models are able to reproduce well most of the population of observed discs in Lupus and Upper Sco, even though the unknown eccentricity of each binary prevents an object by object comparison. In the latter region, the agreement improves if the grain coagulation efficiency is reduced, as may be expected in discs around close binaries. Finally, we mention that thermal winds and sub-structures can be important in explaining few outlying sources.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Super-resolution trends in the ALMA Taurus survey: Structured inner discs and compact discs

The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank') to identify new substructure in three compact and seven extended sources. Among the new features we identify three trends: a higher occurrence rate of substructure in the survey's compact discs than previously seen, underresolved (potentially azimuthally asymmetric) substructure in the innermost disc of extended sources, and a 'shoulder' on the trailing edge of a ring in discs with strong depletion at small radii. Noting the shoulder morphology is present in multiple discs observed at higher resolution, we postulate it is tracing a common physical mechanism. We further demonstrate how a super-resolution frank brightness profile is useful in motivating an accurate parametric model, using the highly structured source DL Tau in which frank finds two new rings. Finally we show that sparse (u, v) plane sampling may be masking the presence of substructure in several additional compact survey sources.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The importance of X-ray frequency in driving photoevaporative winds

Photoevaporative winds are a promising mechanism for dispersing protoplanetary discs, but so far theoretical models have been unable to agree on the relative roles that the X-ray, Extreme Ultraviolet or Far Ultraviolet play in driving the winds. This has been attributed to a variety of methodological differences between studies, including their approach to radiative transfer and thermal balance, the choice of irradiating spectrum employed, and the processes available to cool the gas. We use the \textsc{mocassin} radiative transfer code to simulate wind heating for a variety of spectra on a static density grid taken from simulations of an EUV-driven wind. We explore the impact of choosing a single representative X-ray frequency on their ability to drive a wind by measuring the maximum heated column as a function of photon energy. We demonstrate that for reasonable luminosities and spectra, the most effective energies are at a few $100~\mathrm{eV}$, firmly in the softer regions of the X-ray spectrum, while X-rays with energies $\sim1000~\mathrm{eV}$ interact too weakly with disc gas to provide sufficient heating to drive a wind. We develop a simple model to explain these findings. We argue that further increases in the cooling above our models - for example due to molecular rovibrational lines - may further restrict the heating to the softer energies but are unlikely to prevent X-ray heated winds from launching entirely; increasing the X-ray luminosity has the opposite effect. The various results of photoevaporative wind models should therefore be understood in terms of the choice of irradiating spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modeling the delivery of dust from discs to ionized winds

A necessary first step for dust removal in protoplanetary disc winds is the delivery of dust from the disc to the wind. In the case of ionized winds, the disc and wind are sharply delineated by a narrow ionization front where the gas density and temperature vary by more than an order of magnitude. Using a novel method that is able to model the transport of dust across the ionization front in the presence of disc turbulence, we revisit the problem of dust delivery. Our results show that the delivery of dust to the wind is determined by the vertical gas flow through the disc induced by the mass loss, rather than turbulent diffusion (unless the turbulence is strong, i.e. $α\gtrsim 0.01$). Using these results we provide a simple relation between the maximum size of particle that can be delivered to the wind and the local mass-loss rate per unit area from the wind. This relation is independent of the physical origin of the wind and predicts typical sizes in the 0.01 -- $1\,μm$ range for EUV or X-ray driven winds. These values are a factor $\sim 10$ smaller than those obtained when considering only whether the wind is able to carry away the grains.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Frankenstein: Protoplanetary disc brightness profile reconstruction at sub-beam resolution with a rapid Gaussian process

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planet gap opening across stellar masses

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Proplyds in the Flame Nebula NGC 2024

A recent survey of the inner $0.35\times0.35$pc of the NGC 2024 star forming region revealed two distinct millimetre continuum disc populations that appear to be spatially segregated by the boundary of a dense cloud. The eastern (and more embedded) population is $\sim0.2-0.5$Myr old, with an ALMA mm continuum disc detection rate of about $45\,$per cent. However this drops to only $\sim15$per cent in the 1Myr western population. When presenting this result, van Terwisga et al. (2020) suggested that the two main UV sources, IRS 1 (a B0.5V star in the western region) and IRS 2b (an O8V star in the eastern region, but embedded) have both been evaporating the discs in the depleted western population. In this paper we report the firm discovery in archival HST data of 4 proplyds and 4 further candidate proplyds in NGC 2024, confirming that external photoevaporation of discs is occurring. However, the locations of these proplyds changes the picture. Only three of them are in the depleted western population and their evaporation is dominated by IRS 1, with no obvious impact from IRS 2b. The other 5 proplyds are in the younger eastern region and being evaporated by IRS 2b. We propose that both populations are subject to significant external photoevaporation, which happens throughout the region wherever discs are not sufficiently shielded by the interstellar medium. The external photoevaporation and severe depletion of mm grains in the 0.2-0.5Myr eastern part of NGC 2024 would be in competition even with very early planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unresolved stellar companions with Gaia DR2 astrometry

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Type II migration strikes back -- An old paradigm for planet migration in discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collision velocity of dust grains in self-gravitating protoplanetary discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Determining the mid-plane conditions of circumstellar discs using gas and dust modelling: a study of HD 163296

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

External photoevaporation of protoplanetary discs in sparse stellar groups: the impact of dust growth

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Grand challenges in protoplanetary disc modelling

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photochemical-dynamical models of externally FUV irradiated protoplanetary discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Proplyds around a B1 star - 42 Orionis in NGC 1977

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rapid radiative clearing of protoplanetary discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The minimum mass of detectable planets in protoplanetary discs and the derivation of planetary masses from high resolution observations

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The origin of the eccentricity of the hot Jupiter in CI Tau

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A tidal encounter caught in the act: modelling a star-disc fly-by in the young RW Aurigae system

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Binary accretion rates: dependence on temperature and mass-ratio

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dependence of fragmentation in self-gravitating accretion discs on small scale structure

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantification of stochastic fragmentation of self-gravitating discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Smoothed particle hydrodynamics simulations of gas and dust mixtures

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiplicity in Early Stellar Evolution

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The evolution of the mass ratio of accreting binaries: the role of gas temperature

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The theory of globulettes: candidate precursors of brown dwarfs and free floating planets in H II regions

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterising thermal sweeping: a rapid disc dispersal mechanism

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Do all stars in the solar neighbourhood form in clusters? A cautionary note on the use of the distribution of surface densities

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative Transfer in Star Formation: Testing FLD and Hybrid Methods

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reliable estimation of the column density in Smoothed Particle Hydrodynamic simulations

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Periodic Spectroscopic Variability of FU Orionis

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The rapid dispersal of low-mass virialised clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Debris disk size distributions: steady state collisional evolution with P-R drag and other loss processes

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the theory of disc photoevaporation

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing the accuracy of radiative cooling approximations in SPH simulations

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The dynamical state of stellar structure in star-forming regions

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The formation of permanent soft binaries in dispersing clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Collisional formation of very massive stars in dense clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The clearing of discs around late type T Tauri stars: constraints from the infrared two colour plane

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The imprint of photoevaporation on edge-on discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

X-ray irradiated protoplanetary disk atmospheres II: Predictions from models in hydrostatic equilibrium

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

0 citations0 reviewsNo code linkedOpen access