Source author record

T. J. Haworth

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

9 published item(s)

preprint2022arXiv

Discovery of a 500 au Protobinary in the Massive Prestellar Core G11.92-0.61 MM2

We present high-resolution (< ~160 au) Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm observations of the high-mass prestellar core candidate G11.92-0.61 MM2, which reveal that this source is in fact a protobinary system with a projected separation of 505 au. The binary components, MM2E and MM2W, are compact (radii < 140 au) sources within the partially optically thick dust emission with alpha$_{0.9 cm-1.3 mm}$=2.47-2.94. The 1.3 mm brightness temperatures, T$_b$=68.4/64.6 K for MM2E/MM2W, imply internal heating and minimum luminosities L$_*$ > 24.7 Lsun for MM2E and L$_*$ > 12.6 Lsun for MM2W. The compact sources are connected by a "bridge" of lower-surface-brightness dust emission and lie within more extended emission that may correspond to a circumbinary disk. The circumprotostellar gas mass, estimated from ~0.2" resolution VLA 0.9 cm observations assuming optically thin emission, is 6.8+/-0.9 Msun. No line emission is detected towards MM2E and MM2W in our high-resolution 1.3 mm ALMA observations. The only line detected is 13CO J=2-1, in absorption against the 1.3 mm continuum, which likely traces a layer of cooler molecular material surrounding the protostars. We also report the discovery of a highly asymmetric bipolar molecular outflow that appears to be driven by MM2E and/or MM2W in new deep, ~0.5" resolution (1680 au) ALMA 0.82 mm observations. This outflow, traced by low-excitation CH3OH emission, indicates ongoing accretion onto the protobinary system. Overall, the super-Alfvenic models of Mignon-Risse et al. (2021) agree well with the observed properties of the MM2E/MM2W protobinary, suggesting that this system may be forming in an environment with a weak magnetic field.

preprint2022arXiv

Thermal emission from bow shocks II: 3D magnetohydrodynamic models of Zeta Ophiuchi

The nearby, massive, runaway star Zeta Ophiuchi has a large bow shock detected in optical and infrared, and, uniquely among runaway O stars, diffuse X-ray emission is detected from the shocked stellar wind. Here we make the first detailed computational investigation of the bow shock of Zeta Ophiuchi, to test whether a simple model of the bow shock can explain the observed nebula, and to compare the detected X-ray emission with simulated emission maps. We re-analysed archival {\it Chandra} observations of the thermal diffuse X-ray emission from the shocked wind region of the bow shock, finding total unabsorbed X-ray flux (0.3-2 keV band) corresponding to a diffuse luminosity of $L_\mathrm{X}=2.33~(0.79-3.45)\times10^{29}$ergs$^{-1}$. 3D MHD simulations were used to model the interaction of the star's wind with a uniform ISM using a range of stellar and ISM parameters motivated by observational constraints. Synthetic infrared, Ha, soft X-ray, emission measure, and radio 6\,GHz emission maps were generated from three simulations, for comparison with relevant observations. Simulations where the space velocity of Zeta Ophiuchi has a significant radial velocity produce infrared emission maps with opening angle of the bow shock in better agreement with observations than for the case where motion is fully in the plane of the sky. The simulation with the highest pressure has the closest match, with flux level within a factor of 2 of the observational lower limit, and emission weighted temperature of $\log_{10}(T_\mathrm{A}/\mathrm{K})=6.4$, although the morphology of the diffuse emission appears somewhat different. Observed X-ray emission is a filled bubble brightest near the star whereas simulations predict brightening towards the contact discontinuity as density increases.

preprint2020arXiv

Searching for wide-orbit gravitational instability protoplanets with ALMA in the dust continuum

Searches for young gas giant planets at wide separations have so far focused on techniques appropriate for compact (Jupiter sized) planets. Here we point out that protoplanets born through Gravitational Instability (GI) may remain in an initial pre-collapse phase for as long as the first $ 10^5-10^7$ years after formation. These objects are hundreds of times larger than Jupiter and their atmospheres are too cold ($T\sim$ tens of K) to emit in the NIR or H$α$ via accretion shocks. However, it is possible that their dust emission can be detected with ALMA, even around Class I and II protoplanetary discs. In this paper we produce synthetic observations of these protoplanets. We find that making a detection in a disc at 140 parsecs would require a few hundred minutes of ALMA band 6 observation time. Protoplanets with masses of 3-5 $M_J$ have the highest chance of being detected; less massive objects require unreasonably long observation times (1000 minutes) while more massive ones collapse into giant planets before $10^5$ years. We propose that high resolution surveys of young ($10^5-10^6$ years), massive and face on discs offer the best chance for observing protoplanets. Such a detection would help to place constraints on the protoplanet mass spectrum, explain the turnover in the occurrence frequency of gas giants with system metallicity and constrain the prevalence of GI as a planet formation mechanism. Consistent lack of detection would be evidence against GI as a common planet formation mechanism.

preprint2016arXiv

Violent environment of the inner disk of RW Aur A probed by the 2010 and 2015 dimming events

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($Δ$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

preprint2015arXiv

Isolating signatures of major cloud-cloud collisions II: The lifetimes of broad bridge features

We investigate the longevity of broad bridge features in position-velocity diagrams that appear as a result of cloud-cloud collisions. Broad bridges will have a finite lifetime due to the action of feedback, conversion of gas into stars and the timescale of the collision. We make a series of analytic arguments with which to estimate these lifetimes. Our simple analytic arguments suggest that for collisions between clouds larger than R~10 pc the lifetime of the broad bridge is more likely to be determined by the lifetime of the collision rather than the radiative or wind feedback disruption timescale. However for smaller clouds feedback becomes much more effective. This is because the radiative feedback timescale scales with the ionising flux Nly as R^{7/4}Nly^{-1/4} so a reduction in cloud size requires a relatively large decrease in ionising photons to maintain a given timescale. We find that our analytic arguments are consistent with new synthetic observations of numerical simulations of cloud-cloud collisions (including star formation and radiative feedback). We also argue that if the number of observable broad bridges remains ~ constant, then the disruption timescale must be roughly equivalent to the collision rate. If this is the case our analytic arguments also provide collision rate estimates, which we find are readily consistent with previous theoretical models at the scales they consider (clouds larger than about 10 pc) but are much higher for smaller clouds.

preprint2015arXiv

Isolating signatures of major cloud-cloud collisions using position-velocity diagrams

Collisions between giant molecular clouds are a potential mechanism for triggering the formation of massive stars, or even super star clusters. The trouble is identifying this process observationally and distinguishing it from other mechanisms. We produce synthetic position-velocity diagrams from models of: cloud-cloud collisions, non-interacting clouds along the line of sight, clouds with internal radiative feedback and a more complex cloud evolving in a galactic disc, to try and identify unique signatures of collision. We find that a broad bridge feature connecting two intensity peaks, spatially correlated but separated in velocity, is a signature of a high velocity cloud-cloud collision. We show that the broad bridge feature is resilient to the effects of radiative feedback, at least to around 2.5Myr after the formation of the first massive (ionising) star. However for a head on 10km/s collision we find that this will only be observable from 20-30 per cent of viewing angles. Such broad-bridge features have been identified towards M20, a very young region of massive star formation that was concluded to be a site of cloud-cloud collision by Torii et al (2011), and also towards star formation in the outer Milky Way by Izumi et al (2014).

preprint2015arXiv

StarBench: The D-type expansion of an HII region

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

preprint2015arXiv

The dangers of being trigger--happy

We examine the evidence offered for triggered star formation against the backdrop provided by recent numerical simulations of feedback from massive stars at or below giant molecular cloud sizescales. We compile a catalogue of sixty--seven observational papers, mostly published over the last decade, and examine the signposts most commonly used to infer the presence of triggered star formation. We then determine how well these signposts perform in a recent suite of hydrodynamic simulations of star formation including feedback from O--type stars performed by Dale et al (2012a, b, 2013a, b, 2014). We find that none of the observational markers improve the chances of correctly identifying a given star as triggered by more than factors of two at most. This limits the fidelity of these techniques in interpreting star formation histories. We therefore urge caution in interpreting observations of star formation near feedback--driven structures in terms of triggering.

preprint2015arXiv

TORUS-3DPDR: A self-consistent code treating three-dimensional photoionization and photodissociation regions

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.