Researcher profile

Thomas J. Haworth

Thomas J. Haworth contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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.

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.

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&#39;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.

preprint2020arXiv

A Three Dimensional View of Gomez&#39;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&#39;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.

preprint2020arXiv

Fragmentation favoured in discs around higher mass stars

We investigate how a protoplanetary disc&#39;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.

preprint2020arXiv

Illuminating a tadpole&#39;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.

preprint2020arXiv

Illuminating a tadpole&#39;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.

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.

preprint2020arXiv

Peter Pan Discs: finding Neverland&#39;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.

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.

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.

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.

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.