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

Thomas J. Haworth

Thomas J. Haworth 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.GA

Catalog footprint

What is connected

27works
3topics
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

27 published item(s)

preprint2026arXiv

Formation of multi-planetary systems via pebble accretion in externally photoevaporating discs in stellar clusters

In this paper, we investigate how external photo-evaporation influences the formation, dynamical evolution and the resultant planetary architecture of multi-planet systems born in stellar clusters. We use a model of N-body simulations of multiple planet formation via pebble accretion coupled with a 1-D viscous disc subject to external photo-evaporation. We found that external photo-evaporation reduces the planet growth by reducing the pebble mass reservoir in discs containing multiple planetary embryos across a wide range of disc masses, and is particularly effective in suppressing planet growth in less initially massive discs (< 0.1 M$_{\odot}$). However, in more initially massive ($\geq$ 0.1 M$_{\oplus}$) discs planets lost due to planet-planet interactions dominate the outcome for final resultant total planet mass, masking the effects of external photo-evaporation in curbing the planet mass growth. In terms of the final resulting planetary architectures, the signature of external photo-evaporation is visible in less massive (< 0.1 M$_{\odot}$) discs, with fewer numbers and lower masses of planets surviving in discs irradiated with stronger external FUV radiation. External photo-evaporation also leaves a signature for the wide orbit (> 10 au) terrestrial planets (0.1 - 1 M$_{\oplus}$), with fewer planets populating this region for stronger FUV field. Finally, the 1st-order resonant pairs fraction decreases with stronger FUV radiation, although the resonant pairs occur rarely regardless of the FUV radiation environment, due to the small number of planets that survive gravitational encounters.

0 citations0 reviewsNo code linkedOpen access
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

Determining dispersal mechanisms of protoplanetary disks using accretion and wind mass loss rates

Understanding the origin of accretion and dispersal of protoplanetary disks is fundamental for investigating planet formation. Recent numerical simulations show that launching winds are unavoidable when disks undergo magnetically driven accretion and/or are exposed to external UV radiation. Observations also hint that disk winds are common. We explore how the resulting wind mass loss rate can be used as a probe of both disk accretion and dispersal. As a proof-of-concept study, we focus on magnetocentrifugal winds, MRI (magnetorotational instability) turbulence, and external photoevapotaion. By developing a simple, yet physically motivated disk model and coupling it with simulation results available in the literature, we compute the mass loss rate as a function of external UV flux for each mechanism. We find that different mechanisms lead to different levels of mass loss rate, indicating that the origin of disk accretion and dispersal can be determined, by observing the wind mass loss rate resulting from each mechanism. This determination provides important implications for planet formation. This work thus shows that the ongoing and future observations of the wind mass loss rate for protoplanetary disks are paramount to reliably constrain how protoplanetary disks evolve with time and how planet formation takes place in the disks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dispersal of protoplanetary discs: How stellar properties and the local environment determine the pathway of evolution

We study the evolution and final dispersal of protoplanetary discs that evolve under the action of internal and external photoevaporation, and different degrees of viscous transport. We identify five distinct dispersal pathways, which are i) very long lived discs ($>20\,$Myr), ii) inside-out dispersal where internal photoevaporation dominates and opens inner holes, iii) outside-in dispersal where external photoevaporation dominates through disc truncation and two intermediate regimes characterised by lingering material in the inner disc with the outer disc dispersed predominantly by either internal or external photoevaporation. We determine how the lifetime, relative impact of internal and external winds and clearing pathway varies over a wide, plausible, parameter space of stellar/disc/radiation properties. There are a number of implications, for example in high UV environments because the outer disc lifetime is shorter than the time-scale for clearing the inner disc we do not expect transition discs to be common, which appears to be reflected in the location of transition disc populations towards the Orion Nebular Cluster. Irrespective of environment, we find that ongoing star formation is required to reproduce observed disc fractions as a function of stellar cluster age. This work demonstrates the importance of including both internal and external winds for understanding protoplanetary disc evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The evolution of protoplanetary discs in star formation and feedback simulations

We couple star cluster formation and feedback simulations of a Carina-like star forming region with 1D disc evolutionary models to study the impact of external photoevaporation on disc populations in massive star forming regions. To investigate the effect of shielding of young stellar objects by star forming material, we track the FUV field history at each star in the cluster with two methods: i) Monte Carlo radiative transfer accounting for the shielding of stars from the FUV by the star forming cloud ii) Geometric dilution of the radiation from other stars which ignores shielding effects. We found that significant shielding only occurs for a small fraction of discs and offers protection from external photoevaporation for < 0.5 Myr. However, this initial protection can prevent significant early gas/dust mass loss and disc radius reduction due to external photoevaporation. Particularly, shielding for 0.5 Myr is sufficient for much of the solid reservoir to evolve to larger sizes where it will not be entrained in an external wind. Shielding is therefore potentially significant for terrestrial planet formation in retaining the solid mass budget, but the continued stripping of gas when shielding ends could still impact migration and the gas reservoir for giant planet atmospheres. Our models highlight issues with treating all discs in a cluster with a single characteristic age, since shielded objects are typically only the youngest. Our model predicts that the majority of discs in a 2 Myr Carina-like environment are subject to strong external photoevaporation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The growth and migration of massive planets under the influence of external photoevaporation

The formation of gas giant planets must occur during the first few Myr of a star's lifetime, when the protoplanetary disc still contains sufficient gas to be accreted onto the planetary core. The majority of protoplanetary discs are exposed to strong ultraviolet irradiation from nearby massive stars, which drives winds and depletes the mass budget for planet formation. It remains unclear to what degree external photoevaporation affects the formation of massive planets. In this work, we present a simple one dimensional model for the growth and migration of a massive planet under the influence of external FUV fields. We find that even moderate FUV fluxes $F_\mathrm{FUV}\gtrsim 100 \, G_0$ have a strong influence on planet mass and migration. By decreasing the local surface density and shutting off accretion onto the planet, external irradiation suppresses planet masses and halts migration early. The distribution of typical stellar birth environments can therefore produce an anti-correlation between semi-major axis and planet mass, which may explain the apparent decrease in planet occurrence rates at orbital periods $P_\mathrm{orb}\gtrsim 10^3$ days. Even moderate fluxes $F_\mathrm{FUV}$ strongly suppress giant planet formation and inward migration for any initial semi-major axis if the stellar host mass $M_*\lesssim 0.5\, M_\odot$, consistent with findings that massive planet occurrence is much lower around such stars. The outcomes of our prescription for external disc depletion show significant differences to the current approximation adopted in state-of-the-art population synthesis models, motivating future careful treatment of this important process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Three Dimensional View of Gomez's Hamburger

Unraveling the 3D physical structure, the temperature and density distribution, of protoplanetary discs is an essential step if we are to confront simulations of embedded planets or dynamical instabilities. In this paper we focus on Submillimeter Array observations of the edge-on source, Gomez's Hamburger, believed to host an over-density hypothesised to be a product of gravitational instability in the disc, GoHam b. We demonstrate that, by leveraging the well characterised rotation of a Keplerian disc to deproject observations of molecular lines in position-position-velocity space into disc-centric coordinates, we are able to map out the emission distribution in the (r, z) plane and (x, y, z) space. We show that 12CO traces an elevated layer of $z\,/\,r \sim 0.3$, while 13CO traces deeper in the disc at $z\,/\,r \lesssim 0.2$. We localize emission associated with GoHam b, finding it at deprojected radius of approximately 500 au and at polar angle of +\- 30 degrees from the disc major axis. At the spatial resolution of $\sim 1.5^{\prime\prime}$, GoHam b is spatially unresolved, with an upper limit to its radius of $<190$~au.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fragmentation favoured in discs around higher mass stars

We investigate how a protoplanetary disc's susceptibility to gravitational instabilities and fragmentation depends on the mass of its host star. We use 1D disc models in conjunction with 3D SPH simulations to determine the critical disc-to-star mass ratios at which discs become unstable against fragmentation, finding that discs become increasingly prone to the effects of self-gravity as we increase the host star mass. The actual limit for stability is sensitive to the disc temperature, so if the disc is optically thin stellar irradiation can dramatically stabilise discs against gravitational instability. However, even when this is the case we find that discs around $2$M$_{\odot}$ stars are prone to fragmentation, which will act to produce wide-orbit giant planets and brown dwarfs. The consequences of this work are two-fold: that low mass stars could in principle support high disc-to-star mass ratios, and that higher mass stars have discs that are more prone to fragmentation, which is qualitatively consistent with observations that favour high-mass wide-orbit planets around higher mass stars. We also find that the initial masses of these planets depends on the temperature in the disc at large radii, which itself depends on the level of stellar irradiation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Illuminating a tadpole's metamorphosis II: observing the on-going transformation with ALMA

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the tadpole, a small globule in the Carina Nebula that hosts the HH 900 jet+outflow system. Our data include $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1, $^{13}$CO, C$^{18}$O J=3-2, and serendipitous detections of DCN J=3-2 and CS J=7-6. With angular resolution comparable to the Hubble Space Telescope (HST), our data reveal for the first time the bipolar molecular outflow in CO, seen only inside the globule, that is launched from the previously unseen jet-driving protostar (the HH 900 YSO). The biconical morphology joins smoothly with the externally irradiated outflow seen in ionized gas tracers outside the globule, tracing the overall morphology of a jet-driven molecular outflow. Continuum emission at the location of the HH 900 YSO appears to be slightly flattened perpendicular to outflow axis. Model fits to the continuum have a best-fit spectral index of $\sim 2$, suggesting cold dust and the onset of grain growth. In position-velocity space, $^{13}$CO and C$^{18}$O gas kinematics trace a C-shaped morphology, similar to infall profiles seen in other sources, although the global dynamical behaviour of the gas remains unclear. Line profiles of the CO isotopologues display features consistent with externally heated gas. We estimate a globule mass of $\sim 1.9$ M$_{\odot}$, indicating a remaining lifetime of $\sim 4$ Myr, assuming a constant photoevaporation rate. This long globule lifetime will shield the disk from external irradiation perhaps prolonging its life and enabling planet formation in regions where disks are typically rapidly destroyed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Illuminating a tadpole's metamorphosis III: quantifying past and present environmental impact

We combine MUSE and ALMA observations with theoretical models to evaluate how a tadpole-shaped globule located in the Carina Nebula has been influenced by its environment. This globule is now relatively small (radius ~2500 au), hosts a protostellar jet+outflow (HH 900) and, with a blue-shifted velocity of ~10 km/s, is travelling faster than it should be if its kinematics were set by the turbulent velocity dispersion of the precursor cloud. Its outer layers are currently still subject to heating, but comparing the internal and external pressures implies that the globule is in a post-collapse phase. Intriguingly the outflow is bent, implying that the YSO responsible for launching it is comoving with the globule, which requires that the star formed after the globule was up to speed since otherwise it would have been left behind. We conclude that the most likely scenario is one in which the cloud was much larger before being subject to radiatively-driven implosion, which accelerated the globule to the high observed speeds under the photoevaporative rocket effect and triggered the formation of the star responsible for the outflow. The globule may now be in a quasi-steady state following collapse. Finally, the HH 900 YSO is likely $\gtrsim$1 M$_{\odot}$ and may be the only star forming in the globule. It may be that this process of triggered star formation has prevented the globule from fragmenting to form multiple stars (e.g., due to heating) and has produced a single higher mass star.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massive discs around low-mass stars

We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Peter Pan Discs: finding Neverland's parameters

Peter Pan discs are a recently discovered class of long-lived discs around low-mass stars that survive for an order of magnitude longer than typical discs. In this paper we use disc evolutionary models to determine the required balance between initial conditions and the magnitude of dispersal processes for Peter Pan discs to be primordial. We find that we require low transport ($α\sim10^{-4}$), extremely low external photoevaporation ($\leq10^{-9}{\rm M_{\odot}/yr}$), and relatively high disc masses ($>0.25M_*$) to produce discs with ages and accretion rates consistent with Peter Pan discs. Higher transport ($α= 10^{-3}$) results in disc lifetimes that are too short and even lower transport ($α= 10^{-5}$) leads to accretion rates smaller than those observed. The required external photoevaporation rates are so low that primordial Peter Pan discs will have formed in rare environments on the periphery of low-mass star-forming regions, or deeply embedded, and as such have never subsequently been exposed to higher amounts of UV radiation. Given that such an external photoevaporation scenario is rare, the required disc parameters and accretion properties may reflect the initial conditions and accretion rates of a much larger fraction of the discs around low-mass stars.

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

Testing viscous disc theory using the balance between stellar accretion and external photoevaporation of protoplanetary discs

The nature and rate of (viscous) angular momentum transport in protoplanetary discs (PPDs) has important consequences for the formation process of planetary systems. While accretion rates onto the central star yield constraints on such transport in the inner regions of a PPD, empirical constraints on viscous spreading in the outer regions remain challenging to obtain. Here we demonstrate a novel method to probe the angular momentum transport at the outer edge of the disc. This method applies to PPDs that have lost a significant fraction of their mass due to thermal winds driven by UV irradiation from a neighbouring OB star. We demonstrate that this external photoevaporation can explain the observed depletion of discs in the 3-5 Myr old $σ$ Orionis region, and use our model to make predictions motivating future empirical investigations of disc winds. For populations of intermediate-age PPDs, in viscous models we show that the mass flux outwards due to angular momentum redistribution is balanced by the mass-loss in the photoevaporative wind. A comparison between wind mass-loss and stellar accretion rates therefore offers an independent constraint on viscous models in the outer regions of PPDs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Observational Anatomy of Externally Photoevaporating Planet-Forming Discs I: Atomic Carbon

We demonstrate the utility of CI as a tracer of photoevaporative winds that are being driven from discs by their ambient UV environment. Commonly observed CO lines only trace these winds in relatively weak UV environments and are otherwise dissociated in the wind at the intermediate to high UV fields that most young stars experience. However, CI traces unsubtle kinematic signatures of a wind in intermediate UV environments ($\sim1000$G$_0$) and can be used to place constraints on the kinematics and temperature of the wind. In CI position-velocity diagrams external photoevaporation results in velocities that are faster than those from Keplerian rotation alone, as well as emission from quadrants of position-velocity space in which there would be no Keplerian emission. This is independent of viewing angle because the wind has components that are perpendicular to the azimuthal rotation of the disc. At intermediate viewing angles ($\sim 30-60^\circ$) moment 1 maps also exhibit a twisted morphology over large scales (unlike other processes that result in twists, which are typically towards the inner disc). CI is readily observable with ALMA, which means that it is now possible to identify and characterise the effect of external photoevaporation on planet-forming discs in intermediate UV environments.

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

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

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

Detecting stellar-wind bubbles through infrared arcs in HII regions

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the relative importance of different microphysics on the D-type expansion of galactic HII regions

Radiation hydrodynamics (RHD) simulations are used to study many astrophysical phenomena, however they require the use of simplified radiation transport and thermal prescriptions to reduce computational cost. In this paper we present a systematic study of the importance of microphysical processes in RHD simulations using the example of D-type HII region expansion. We compare the simplest hydrogen-only models with those that include: ionisation of H, He, C, N, O, S and Ne, different gas metallicity, non-LTE metal line blanketed stellar spectral models of varying metallicity, radiation pressure, dust and treatment of photodissociation regions. Each of these processes are explicitly treated using modern numerical methods rather than parameterisation. In line with expectations, changes due to microphysics in either the effective number of ionising photons or the thermal structure of the gas lead to differences in D-type expansion. In general we find that more realistic calculations lead to the onset of D-type expansion at smaller radii and a slower subsequent expansion. Simulations of star forming regions using simplified microphysics are therefore likely overestimating the strength of radiative feedback. We find that both variations in gas metallicity and the inclusion of dust can affect the ionisation front evolution at the 10-20 per cent level over 500kyr, which could substantially modify the results of simplified 3D models including feedback. Stellar metallicity, radiation pressure and the inclusion of photodissociation regions are all less significant effects at the 1 per cent level or less, rendering them of minor importance in the modelling the dynamical evolution of HII regions.

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

Assessing molecular line diagnostics of triggered star formation using synthetic observations

We investigate observational signatures of triggered star formation in bright rimmed clouds (BRCs) by using molecular line transfer calculations based on radiation-hydrodynamic radiatively-driven-implosion models. We find that for BRCs the separation in velocity between the line profile peak of an optically thick and an optically thin line is determined by both the observer viewing angle and the density of the shell driving into the cloud. In agreement with observations, we find that most BRC line profiles are symmetric and that asymmetries can be either red or blue, in contrast to the blue-dominance expected for a collapsing cloud. Asymmetries in the line profiles arise when an optically thick line is dominated by the shell and an optically thin line is dominated by the cloud interior to the shell. The asymmetries are red or blue depending on whether the shell is moving towards or away from the observer respectively. Using the known motions of the molecular gas in our models we rule out the envelope expansion with core collapse mechanism as the cause of the lack of blue-asymmetry in our simulated observations. We show that the absence of a strong photon dominated region (PDR) around a BRC may not rule out the presence of triggered star formation: if the BRC line profile has a strong blue component then the shell is expected to be driving towards the observer, suggesting that the cloud is being viewed from behind and the PDR is obstructed. This could explain why BRCs such as SFO 80, 81 and 86 have a blue secondary peak and only a weak PDR inferred at 8 microns. Finally we also test the use of 12CO, 13CO and C18O as diagnostics of cloud mass, temperature and column density. We find that the inferred conditions are in reasonable agreement with those from the models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Testing diagnostics of triggered star formation

We produce synthetic images and SEDs from radiation hydrodynamical simulations of radiatively driven implosion. The synthetically imaged bright rimmed clouds (BRCs) are morphologically similar to those observed in star forming regions. Using nebular diagnostic line-ratios, simulated Very Large Array (VLA) radio images, Hα imaging and SED fitting we compute the neutral cloud and ionized boundary layer gas densities and temperatures and perform a virial stability analysis for each model cloud. We determine that the neutral cloud temperatures derived by SED fitting are hotter than the dominant neutral cloud temperature by 1 - 2 K due to emission from warm dust. This translates into a change in the calculated cloud mass by 8-35 %. Using a constant mass conversion factor (Cν) for BRCs of different class is found to give rise to errors in the cloud mass of up to a factor of 3.6. The ionized boundary layer (IBL) electron temperature calculated using diagnostic line ratios is more accurate than assuming the canonical value adopted for radio diagnostics of 10^4 K. Both radio diagnostics and diagnostic line ratios are found to underestimate the electron density in the IBL. Each system is qualitatively correctly found to be in a state in which the pressure in the ionized boundary layer is greater than the supporting cloud pressure, implying that the objects are being compressed. We find that observationally derived mass loss estimates agree with those on the simulation grid and introduce the concept of using the mass loss flux to give an indication of the relative strength of photo-evaporative flow between clouds. The effect of beam size on these diagnostics in radio observations is found to be a mixing of the bright rim and ambient cloud and HII region fluxes, which leads to an underestimate of the cloud properties relative to a control diagnostic.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Testing models of triggered star formation: theory and observation

One of the main reasons that triggered star formation is contentious is the failure to accurately link the observations with models in a detailed, quantitative, way. It is therefore critical to continuously test and improve the model details and methods with which comparisons to observations are made. We use a Monte Carlo radiation transport and hydrodynamics code TORUS to show that the diffuse radiation field has a significant impact on the outcome of radiatively driven implosion (RDI) models. We also calculate SEDs and synthetic images from the models to test observational diagnostics that are used to determine bright rimmed cloud conditions and search for signs of RDI.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation hydrodynamics of triggered star formation: the effect of the diffuse radiation field

We investigate the effect of including diffuse field radiation when modelling the radiatively driven implosion of a Bonnor-Ebert sphere (BES). Radiation-hydrodynamical calculations are performed by using operator splitting to combine Monte Carlo photoionization with grid-based Eulerian hydrodynamics that includes self-gravity. It is found that the diffuse field has a significant effect on the nature of radiatively driven collapse which is strongly coupled to the strength of the driving shock that is established before impacting the BES. This can result in either slower or more rapid star formation than expected using the on-the-spot approximation depending on the distance of the BES from the source object. As well as directly compressing the BES, stronger shocks increase the thickness and density in the shell of accumulated material, which leads to short, strong, photo-evaporative ejections that reinforce the compression whenever it slows. This happens particularly effectively when the diffuse field is included as rocket motion is induced over a larger area of the shell surface. The formation and evolution of 'elephant trunks' via instability is also found to vary significantly when the diffuse field is included. Since the perturbations that seed instabilities are smeared out elephant trunks form less readily and, once formed, are exposed to enhanced thermal compression.

0 citations0 reviewsNo code linkedOpen access