Showing 1–50 of 80 results for author: Harries, T J

Imaging the warped dusty disk wind environment of SU Aurigae with MIRC-X

Authors: Aaron Labdon, Stefan Kraus, Claire L. Davies, Alexander Kreplin, Sebastian Zarrilli, John D. Monnier, Jean-Baptiste le Bouquin, Narsireddy Anugu, Benjamin Setterholm, Tyler Gardner, Jacob Ennis, Cyprien Lanthermann, Theo ten Brummelaar, Gail Schaefer, Tim J. Harries

Abstract: SU Aurigae is a widely studied T Tauri star and here we present original state-of-the-art interferometric observations with better uv and baseline coverage than previous studies. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The MIRC-X instrument at CHARA is a 6 telescope optical beam co… ▽ More SU Aurigae is a widely studied T Tauri star and here we present original state-of-the-art interferometric observations with better uv and baseline coverage than previous studies. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The MIRC-X instrument at CHARA is a 6 telescope optical beam combiner offering baselines up to 331 m. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrains the physical parameters of the disk. Image reconstruction reveals a highly inclined disk with a slight asymmetry consistent with inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission/scattering of the far-side. Geometric models find that the underlying brightness distribution is best modelled as a Gaussian with a FWHM of $1.53\pm0.01 \mathrm{mas}$ at an inclination of $56.9\pm0.4^\circ$ and minor axis position angle of $55.9\pm0.5^\circ$. Radiative transfer modelling shows a flared disk with an inner radius at 0.16 au which implies a grain size of $0.14 \mathrm{μm}$ assuming astronomical silicates and a scale height of 9.0 au at 100 au. In agreement with literature, only the dusty disk wind successfully accounts for the NIR excess by introducing dust above the mid-plane. Our results confirm and provide better constraints than previous inner disk studies of SU Aurigae. We confirm the presence of a dusty disk wind in the cicumstellar environment, the strength of which is enhanced by a late infall event which also causes very strong misalignments between the inner and outer disks. △ Less

Submitted 9 June, 2023; originally announced June 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2111.06205, arXiv:1905.11907

Journal ref: A&A 678, A6 (2023)

Hydrogen Emission from Accretion and Outflow in T Tauri Stars

Authors: Tom J. G. Wilson, S. Matt, T. J. Harries, G. J. Herczeg

Abstract: Radiative transfer modelling offers a powerful tool for understanding the enigmatic hydrogen emission lines from T Tauri stars. This work compares optical and near-IR spectroscopy of 29 T Tauri stars with our grid of synthetic line profiles. The archival spectra, obtained with VLT's X-Shooter, provide simultaneous coverage of many optical and infrared hydrogen lines. The observations exhibit simil… ▽ More Radiative transfer modelling offers a powerful tool for understanding the enigmatic hydrogen emission lines from T Tauri stars. This work compares optical and near-IR spectroscopy of 29 T Tauri stars with our grid of synthetic line profiles. The archival spectra, obtained with VLT's X-Shooter, provide simultaneous coverage of many optical and infrared hydrogen lines. The observations exhibit similar morphologies of line profiles seen in other studies. We used the radiative transfer code TORUS to create synthetic H$α$, Pa$β$, Pa$γ$, and Br$γ$ emission lines for a fiducial T Tauri model that included axisymmetric magnetospheric accretion and a polar stellar wind. The distribution of Reipurth types and line widths for the synthetic H$α$ lines is similar to the observed results. However, the modelled infrared lines are narrower than the observations by $\approx 80{~\rm kms}^{-1}$, and our models predict a significantly higher proportion ($\approx 90$ per cent) of inverse P-Cygni profiles. Furthermore, our radiative transfer models suggest that the frequency of P-Cygni profiles depends on the ratio of the mass loss to mass accretion rates and blue-shifted sub-continuum absorption was predicted for mass loss rates as low as $10^{-12}~M_{\odot}{\rm~ yr}^{-1}$. We explore the effect of rotation, turbulence, and the contributions from red-shifted absorption in an attempt to explain the discrepancy in widths. Our findings show that, singularly, none of these effects is sufficient to explain the observed disparity. However, a combination of rotation, turbulence, and non-axisymmetric accretion may improve the fit of the models to the observed data. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

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

Authors: C. J. Cyganowski, J. D. Ilee, C. L. Brogan, T. R. Hunter, S. Zhang, T. J. Harries, T. J. Haworth

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 25 April, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: 11 pages, 3 figures, 2 tables; accepted for publication in ApJL, no changes to manuscript

Scattering and sublimation: a multi-scale view of $μ$m-sized dust in the inclined disc of HD 145718

Authors: Claire L. Davies, Evan A. Rich, Tim J. Harries, John D. Monnier, Anna S. E. Laws, Sean M. Andrews, Jaehan Bae, David J. Wilner, Narsireddy Anugu, Jacob Ennis, Tyler Gardner, Stefan Kraus, Aaron Labdon, Jean-Baptiste le Bouquin, Cyprien Lanthermann, Gail H. Schaefer, Benjamin R. Setterholm, Theo ten Brummelaar

Abstract: We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angl… ▽ More We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angle, PA\,$=-1.0-0.6^{\circ}$, consistent with previous estimates. The NIR scattering surface extends out to $\sim75\,$au and we infer an aspect ratio, $h_{\rm{scat}}(r)/r\sim0.24$ in $J$-band; $\sim0.22$ in $H$-band. Our GPI images and VLTI+CHARA NIR interferometry suggest that the disc surface layers are populated by grains $\gtrsim λ/2π$ in size, indicating these grains are aerodynamically supported against settling and/or the density of smaller grains is relatively low. We demonstrate that our geometric analysis provides a reasonable assessment of the height of the NIR scattering surface at the outer edge of the disc and, if the inclination can be independently constrained, has the potential to probe the flaring exponent of the scattering surface in similarly inclined ($i\gtrsim70^{\circ}$) discs. In re-evaluating HD~145718's stellar properties, we found that the object's dimming events - previously characterised as UX~Or and dipper variability - are consistent with dust occultation by grains larger, on average, than found in the ISM. This occulting dust likely originates close to the inferred dust sublimation radius at $0.17\,$au. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 18 pages, 16 figures; Accepted for publication in MNRAS

Investigating the Relative Gas and Small Dust Grain Surface Heights in Protoplanetary Disks

Authors: Evan A. Rich, Richard Teague, John D. Monnier, Claire L. Davies, Arthur Bosman, Tim J. Harries, Nuria Calvet, Fred C. Adams, Zhaohuan Zhu

Abstract: Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $μm$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $β$ of the small dust grain sc… ▽ More Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $μm$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $β$ of the small dust grain scattering surface and $^{12}$CO gas emission surface for three protoplanetary disks IM Lup, HD 163296, and HD 97048 (CU Cha). For two systems, IM Lup and HD 163296, the $^{12}$CO gas and small dust grains at small radii from the star have similar heights but at larger radii ($>$100 au) the dust grain scattering surface height is lower than the $^{12}$CO gas emission surface height. In the case of HD 97048, the small dust grain scattering surface has similar heights to the $^{12}$CO gas emission surface at all radii. We ran a protoplanetary disk radiative transfer model of a generic protoplanetary disk with TORUS and showed that there is no difference between the observed scattering surface and $^{12}$CO emission surface. We also performed analytical modeling of the system and found that gas-to-dust ratios larger than 100 could explain the observed difference in IM Lup and HD 163296. This is the first direct comparison of observations of gas and small dust grain heights distribution in protoplanetary disks. Future observations of gas emission and near-infrared scattered light instruments are needed to look for similar trends in other protoplanetary disks. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: Accepted ApJ, 13 pages, 3 figures, 4 tables

A Dust Trap in the Young Multiple System HD 34700

Authors: Peyton Benac, Luca Matra, David J. Wilner, Maria J. Jimenez-Donaire, John D. Monnier, Evan A. Rich, Tim J. Harries, Anna Laws, Qizhou Zhang

Abstract: Millimeter observations of disks around young stars reveal substructures indicative of gas pressure traps that may aid grain growth and planet formation. We present Submillimeter Array observations of HD 34700- two Herbig Ae stars in a close binary system (Aa/Ab, $\sim$0.25 AU), surrounded by a disk presenting a large cavity and spiral arms seen in scattered light, and two distant, lower mass comp… ▽ More Millimeter observations of disks around young stars reveal substructures indicative of gas pressure traps that may aid grain growth and planet formation. We present Submillimeter Array observations of HD 34700- two Herbig Ae stars in a close binary system (Aa/Ab, $\sim$0.25 AU), surrounded by a disk presenting a large cavity and spiral arms seen in scattered light, and two distant, lower mass companions. These observations include 1.3 mm continuum emission and the $^{12}$CO 2-1 line at $\sim0.5$" (178 AU) resolution. They resolve a prominent azimuthal asymmetry in the continuum, and Keplerian rotation of a circumbinary disk in the $^{12}$CO line. The asymmetry is located at a radius of $155^{+11}_{-7}$ AU, consistent with the edge of the scattered light cavity, being resolved in both radius ($72 ^{+14}_{-15}$ AU) and azimuth (FWHM = $64 ^{\circ +8}_{-7}$). The strong asymmetry in millimeter continuum emission could be evidence for a dust trap, together with the more symmetric morphology of $^{12}$CO emission and small grains. We hypothesize an unseen circumbinary companion, responsible for the cavity in scattered light and creating a vortex at the cavity edge that manifests in dust trapping. The disk mass has limitations imposed by the detection of $^{12}$CO and non-detection of $^{13}$CO. We discuss its consequences for the potential past gravitational instability of this system, likely accounting for the rapid formation of a circumbinary companion. We also report the discovery of resolved continuum emission associated with HD 34700B (projected separation $\sim1850$AU), which we explain through a circumstellar disk. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: 12 pages, 5 figures, accepted to ApJ

The observational impact of dust trapping in self-gravitating discs

Authors: James Cadman, Cassandra Hall, Ken Rice, Tim J. Harries, Pamela D. Klaassen

Abstract: We present a 3D semi-analytic model of self-gravitating discs, and include a prescription for dust trapping in the disc spiral arms. Using Monte-Carlo radiative transfer we produce synthetic ALMA observations of these discs. In doing so we demonstrate that our model is capable of producing observational predictions, and able to model real image data of potentially self-gravitating discs. For a dis… ▽ More We present a 3D semi-analytic model of self-gravitating discs, and include a prescription for dust trapping in the disc spiral arms. Using Monte-Carlo radiative transfer we produce synthetic ALMA observations of these discs. In doing so we demonstrate that our model is capable of producing observational predictions, and able to model real image data of potentially self-gravitating discs. For a disc to generate spiral structure that would be observable with ALMA requires that the disc's dust mass budget is dominated by millimetre and centimetre-sized grains. Discs in which grains have grown to the grain fragmentation threshold may satisfy this criterion, thus we predict that signatures of gravitational instability may be detectable in discs of lower mass than has previously been suggested. For example, we find that discs with disc-to-star mass ratios as low as $0.10$ are capable of driving observable spiral arms. Substructure becomes challenging to detect in discs where no grain growth has occurred or in which grain growth has proceeded well beyond the grain fragmentation threshold. We demonstrate how we can use our model to retrieve information about dust trapping and grain growth through multi-wavelength observations of discs, and using estimates of the opacity spectral index. Applying our disc model to the Elias 27, WaOph 6 and IM Lup systems we find gravitational instability to be a plausible explanation for the observed substructure in all 3 discs, if sufficient grain growth has indeed occurred. △ Less

Submitted 22 August, 2020; originally announced August 2020.

Comments: 19 pages, 21 figures, accepted for publication in MNRAS

The inner disk of RY Tau: evidence of stellar occultation by the disk atmosphere at the sublimation rim from K-band continuum interferometry

Authors: Claire L. Davies, Stefan Kraus, Tim J. Harries, John D. Monnier, Brian Kloppenborg, Alicia Aarnio, Fabien Baron, Rebeca Garcia Lopez, Rafael Millan-Gabet, Robert Parks, Ettore Pedretti, Karine Perraut, Judit Sturmann, Laszlo Sturmann, Theo A. ten Brummelaar, Yamina Touhami

Abstract: We present models of the inner region of the circumstellar disk of RY Tau which aim to explain our near-infrared ($K$-band: $2.1\,μ$m) interferometric observations while remaining consistent with the optical to near-infrared portions of the spectral energy distribution. Our sub-milliarcsecond resolution CHARA Array observations are supplemented with shorter baseline, archival data from PTI, KI and… ▽ More We present models of the inner region of the circumstellar disk of RY Tau which aim to explain our near-infrared ($K$-band: $2.1\,μ$m) interferometric observations while remaining consistent with the optical to near-infrared portions of the spectral energy distribution. Our sub-milliarcsecond resolution CHARA Array observations are supplemented with shorter baseline, archival data from PTI, KI and VLTI/GRAVITY and modeled using an axisymmetric Monte Carlo radiative transfer code. The $K$-band visibilities are well-fit by models incorporating a central star illuminating a disk with an inner edge shaped by dust sublimation at $0.210\pm0.005\,$au, assuming a viewing geometry adopted from millimeter interferometry ($65^{\circ}$ inclined with a disk major axis position angle of $23^{\circ}$). This sublimation radius is consistent with that expected of Silicate grains with a maximum size of $0.36-0.40\,μ$m contributing to the opacity and is an order of magnitude further from the star than the theoretical magnetospheric truncation radius. The visibilities on the longest baselines probed by CHARA indicate that we lack a clear line-of-sight to the stellar photosphere. Instead, our analysis shows that the central star is occulted by the disk surface layers close to the sublimation rim. While we do not see direct evidence of temporal variability in our multi-epoch CHARA observations, we suggest the aperiodic photometric variability of RY~Tau is likely related temporal and/or azimuthal variations in the structure of the disk surface layers. △ Less

Submitted 14 May, 2020; originally announced May 2020.

Comments: Accepted for publication in The Astrophysical Journal

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Authors: Stefan Kraus, Alexander Kreplin, Alison K. Young, Matthew R. Bate, John D. Monnier, Tim J. Harries, Henning Avenhaus, Jacques Kluska, Anna S. E. Laws, Evan A. Rich, Matthew Willson, Alicia N. Aarnio, Fred C. Adams, Sean M. Andrews, Narsireddy Anugu, Jaehan Bae, Theo ten Brummelaar, Nuria Calvet, Michel Curé, Claire L. Davies, Jacob Ennis, Catherine Espaillat, Tyler Gardner, Lee Hartmann, Sasha Hinkley , et al. (7 additional authors not shown)

Abstract: Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple… ▽ More Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits. △ Less

Submitted 4 September, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 63 pages, 4+13 Figures, 6 Tables, published at Science

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

Authors: Anna S. E. Laws, Tim J. Harries, Benjamin R. Setterholm, John D. Monnier, Evan A. Rich, Alicia N. Aarnio, Fred C. Adams, Sean Andrews, Jaehan Bae, Nuria Calvet, Catherine Espaillat, Lee Hartmann, Sasha Hinkley, Andrea Isella, Stefan Kraus, David Wilner, Zhaohuan Zhu

Abstract: We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagr… ▽ More We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Massive star feedback in clusters: variation of the FUV interstellar radiation field in time and space

Authors: Ahmad A. Ali, Tim J. Harries

Abstract: We investigate radiative feedback from a 34 M$_\odot$ star in a $10^4$ M$_\odot$ turbulent cloud using three-dimensional radiation-hydrodynamics (RHD) models. We use Monte Carlo radiative transfer to accurately compute photoionization equilibrium and radiation pressure, with multiple atomic species and silicate dust grains. We include the diffuse radiation field, dust absorption/re-emission, and s… ▽ More We investigate radiative feedback from a 34 M$_\odot$ star in a $10^4$ M$_\odot$ turbulent cloud using three-dimensional radiation-hydrodynamics (RHD) models. We use Monte Carlo radiative transfer to accurately compute photoionization equilibrium and radiation pressure, with multiple atomic species and silicate dust grains. We include the diffuse radiation field, dust absorption/re-emission, and scattering. The cloud is efficiently dispersed, with 75 per cent of the mass leaving the (32.3 pc)$^3$ grid within 4.3 Myr (1.1 $t_{ff}$). This compares to all mass exiting within 1.6 Myr (0.74 $t_{ff}$) in our previously published $10^3$ M$_\odot$ cloud. At most 20 per cent of the mass is ionized, compared to 40 per cent in the lower mass model, despite the ionized volume fraction being 80 per cent in both, implying the higher mass cloud is more resilient to feedback. The total Jeans-unstable mass increases linearly up to 1500 M$_\odot$ before plateauing after 2 Myr, corresponding to a core formation efficiency of 15 per cent. We also measure the time-variation of the far-ultraviolet (FUV) radiation field, $G_0$, impinging on other cluster members, taking into account for the first time how this changes in a dynamic cluster environment with intervening opacity sources and stellar motions. Many objects remain shielded in the first 0.5 Myr whilst the massive star is embedded, after which $G_0$ increases by orders of magnitude. Gas motions later on cause comparable drops which happen instantaneously and last for $\sim$ 1 Myr before being restored. This highly variable UV field will influence the photoevaporation of protoplanetary discs near massive stars. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Comments: 12 pages, 12 figures, accepted for publication in MNRAS

Synthetic molecular line observations of the first hydrostatic core from chemical calculations

Authors: Alison K. Young, Matthew R. Bate, Tim J. Harries, David M. Acreman

Abstract: The first stable object to develop in the low-mass star formation process has long been predicted to be the first hydrostatic core (FHSC). Despite much effort, it has still yet to be definitively observed in nature. More specific observational signatures are required to enable observers to distinguish the FHSC from young, faint, but more evolved protostars. Here we present synthetic spectral line… ▽ More The first stable object to develop in the low-mass star formation process has long been predicted to be the first hydrostatic core (FHSC). Despite much effort, it has still yet to be definitively observed in nature. More specific observational signatures are required to enable observers to distinguish the FHSC from young, faint, but more evolved protostars. Here we present synthetic spectral line observations for CO, SO, CS and HCO$^+$ that were calculated from radiation (magneto)hydrodynamical models, chemical modelling and Monte Carlo radiative transfer. HCO$^+$ $(1-0)$ and SO $(8_7 - 7_6)$ spectra of the FHSC show variations for observations at a low inclination which may allow a candidate FHSC to be distinguished from a more evolved object. We find that the FHSC outflow is unlikely to be detectable with ALMA, which would discount the observed sources with slow outflows that are currently identified as candidate FHSCs. We compare the results of simulated ALMA observations with observed candidate FHSCs and recommend Oph A SM1N and N6-mm as the most promising candidates to follow up. △ Less

Submitted 4 June, 2019; originally announced June 2019.

Comments: Accepted to MNRAS. 23 pages, 25 figures

Dusty disk winds at the sublimation rim of the highly inclined, low mass YSO SU Aurigae

Authors: Aaron Labdon, Stefan Kraus, Claire L Davies, Alexander Kreplin, Jacques Kluska, Tim J Harries, John D Monnier, Theo ten Brummelaar, Fabian Baron, Rafael Millan-Gabet, Brian Kloppenborg, Joshua Eisner, Judit Sturmann, Laszlo Sturmann

Abstract: T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer 3 x higher resolution and better baseline position angle coverage over previous observations. We investigate the characteristics of circumstellar m… ▽ More T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer 3 x higher resolution and better baseline position angle coverage over previous observations. We investigate the characteristics of circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. The CHARA array offers opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB 3-telescope combiner in the K-band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and geometric modeling. Additionally, the fitting of radiative transfer models constrains the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Image reconstruction reveals an inclined disk with slight asymmetry along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission/scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17+/-0.02 au with an inclination of 50.9+/-1.0 degrees and minor axis position angle 60.8+/-1.2 degrees. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 um and a scale height of 15.0 au at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K-band excess by introducing dust above the mid-plane. △ Less

Submitted 28 May, 2019; originally announced May 2019.

Comments: Accepted for publication in Astronomy \& Astrophysics

Journal ref: A&A 627, A36 (2019)

Multiple spiral arms in the disk around intermediate-mass binary HD 34700A

Authors: John D. Monnier, Tim J. Harries, Jaehan Bae, Benjamin R. Setterholm, Anna Laws, Alicia Aarnio, Fred C. Adams, Sean Andrews, Nuria Calvet, Catherine Espaillat, Lee Hartmann, Stefan Kraus, Melissa McClure, Chris Miller, Rebecca Oppenheimer, David Wilner, Zhaohuan Zhu

Abstract: We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearl… ▽ More We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: Accepted by the Astrophysical Journal on Dec 15, 2018

The temporal requirements of directly observing self-gravitating spiral waves in protoplanetary discs with ALMA

Authors: Cassandra Hall, Ruobing Dong, Ken Rice, Tim J. Harries, Joan Najita, Richard Alexander, Sean Brittain

Abstract: We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We… ▽ More We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_\mathrm{disc}/M_*\gtrsim 0.25$ and $R_\mathrm{outer}\lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $\sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales △ Less

Submitted 8 June, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: 13 pages, 6 figures

Journal ref: 2019ApJ...871..228H

Imaging the disc rim and a moving close-in companion candidate in the pre-transitional disc of V1247 Orionis

Authors: Matthew Willson, Stefan Kraus, Jacques Kluska, John D. Monnier, Michel Cure, Mike Sitko, Alicia Aarnio, Michael J. Ireland, Aaron Rizzuto, Edward Hone, Alexander Kreplin, Sean Andrews, Nuria Calvet, Catherine Espaillat, Misato Fukagawa, Tim J. Harries, Sasha Hinkley, Samer Kanaan, Takayuki Muto, David J. Wilner

Abstract: V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to… ▽ More V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to identify the origin of the previously detected asymmetries. We have observed V1247 Ori at three epochs spanning $\sim678$ days using sparse aperture masking interferometry with Keck/NIRC2 and VLT/NACO. In addition, we search for signs of accretion through VLT/SPHERE-ZIMPOL spectral differential imaging in H$α$ and R-band continuum. Our SMA sub-millimetre interferometry in 880 $μ$m continuum and in the CO(3-2) line allows us to constrain the orientation and direction of rotation of the outer disc. We find the L'-band emission to be dominated by static features which trace forward-scattered dust emission from the inner edge of the outer disc located to the north-east. In H- and K-band, we see evidence for a companion candidate that moved systematically by 45$^{\circ}$ within the first $\sim$345 days. The separation of the companion candidate is not well constrained, but the observed position angle change is consistent with Keplerian motion of a body located on a 6 au orbit. From the SMA CO moment map, the location of the disc rim, and the detected orbital motion, we deduced the three-dimensional orientation of the disc. We see no indication of accretion in H$α$ and set upper limits for an accreting companion. The measured contrast of the companion candidate in H and K is consistent with an actively accreting protoplanet. Hence, we identify V1247 Ori as a unique laboratory for studying companion-disc interactions and disc clearing. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 17 pages, 14 figures

Journal ref: A&A 621, A7 (2019)

G11.92-0.61 MM 1: A fragmented Keplerian disk surrounding a proto-O star

Authors: J. D. Ilee, C. J. Cyganowski, C. L. Brogan, T. R. Hunter, D. H. Forgan, T. J. Haworth, C. J. Clarke, T. J. Harries

Abstract: We present high resolution ($\sim$300 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment of MM 1 in 1.3 mm continuum emission and CH$_{3}$CN emission for the first time. The object divides into two main sources - MM 1a, which is the source of a bipolar molecular outflow, and M… ▽ More We present high resolution ($\sim$300 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment of MM 1 in 1.3 mm continuum emission and CH$_{3}$CN emission for the first time. The object divides into two main sources - MM 1a, which is the source of a bipolar molecular outflow, and MM 1b, located 0.57'' (1920 au) to the South-East. The main component of MM 1a is an elongated continuum structure, perpendicular to the bipolar outflow, with a size of $0.141'' \times 0.050''$ ($480\times170$ au). The gas kinematics toward MM 1a probed via CH$_{3}$CN trace a variety of scales. The lower energy $J=12-11$ $K=3$ line traces extended, rotating gas within the outflow cavity, while the $v$8=1 line shows a clearly-resolved Keplerian rotation signature. Analysis of the gas kinematics and dust emission shows that the total enclosed mass in MM 1a is $40\pm5$ M$_{\odot}$ (where between 2.2-5.8 M$_{\odot}$ is attributed to the disk), while MM 1b is $<0.6$ M$_{\odot}$. The extreme mass ratio and orbital properties of MM 1a and MM 1b suggest that MM 1b is one of the first observed examples of the formation of a binary star via disk fragmentation around a massive young (proto)star. △ Less

Submitted 13 November, 2018; originally announced November 2018.

Comments: 8 pages, 3 figures, accepted for publication in ApJ Letters

Simultaneous spectral energy distribution and near-infrared interferometry modeling of HD 142666

Authors: Claire L Davies, Stefan Kraus, Tim J Harries, Alexander Kreplin, John D Monnier, Aaron Labdon, Brian Kloppenborg, David M Acreman, Fabien Baron, Rafael Millan-Gabet, Judit Sturmann, Laszlo Sturmann, Theo A Ten Brummelaar

Abstract: We present comprehensive models of Herbig Ae star, HD 142666, which aim to simultaneously explain its spectral energy distribution (SED) and near-infrared (NIR) interferometry. Our new sub-milliarcsecond resolution CHARA (CLASSIC and CLIMB) interferometric observations, supplemented with archival shorter baseline data from VLTI/PIONIER and the Keck Interferometer, are modeled using centro-symmetri… ▽ More We present comprehensive models of Herbig Ae star, HD 142666, which aim to simultaneously explain its spectral energy distribution (SED) and near-infrared (NIR) interferometry. Our new sub-milliarcsecond resolution CHARA (CLASSIC and CLIMB) interferometric observations, supplemented with archival shorter baseline data from VLTI/PIONIER and the Keck Interferometer, are modeled using centro-symmetric geometric models and an axisymmetric radiative transfer code. CHARA's 330 m baselines enable us to place strong constraints on the viewing geometry, revealing a disk inclined at 58 degrees from face-on with a 160 degree major axis position angle. Disk models imposing vertical hydrostatic equilibrium provide poor fits to the SED. Models accounting for disk scale height inflation, possibly induced by turbulence associated with magneto-rotational instabilities, and invoking grain growth to >1 micron size in the disk rim are required to simultaneously reproduce the SED and measured visibility profile. However, visibility residuals for our best model fits to the SED indicate the presence of unexplained NIR emission, particularly along the apparent disk minor axis, while closure phase residuals indicate a more centro-symmetric emitting region. In addition, our inferred 58 degree disk inclination is inconsistent with a disk-based origin for the UX Ori-type variability exhibited by HD 142666. Additional complexity, unaccounted for in our models, is clearly present in the NIR-emitting region. We propose the disk is likely inclined toward a more edge-on orientation and/or an optically thick outflow component also contributes to the NIR circumstellar flux. △ Less

Submitted 31 August, 2018; originally announced August 2018.

Comments: Accepted for publication in The Astrophysical Journal

Modelling massive-star feedback with Monte Carlo radiation hydrodynamics: photoionization and radiation pressure in a turbulent cloud

Authors: Ahmad Ali, Tim J. Harries, Thomas A. Douglas

Abstract: We simulate a self-gravitating, turbulent cloud of 1000 Msol with photoionization and radiation pressure feedback from a 34 Msol star. We use a detailed Monte Carlo radiative transfer scheme alongside the hydrodynamics to compute photoionization and thermal equilibrium with dust grains and multiple atomic species. Using these gas temperatures, dust temperatures, and ionization fractions, we produc… ▽ More We simulate a self-gravitating, turbulent cloud of 1000 Msol with photoionization and radiation pressure feedback from a 34 Msol star. We use a detailed Monte Carlo radiative transfer scheme alongside the hydrodynamics to compute photoionization and thermal equilibrium with dust grains and multiple atomic species. Using these gas temperatures, dust temperatures, and ionization fractions, we produce self-consistent synthetic observations of line and continuum emission. We find that all material is dispersed from the (15.5 pc)$^3$ grid within 1.6 Myr or 0.74 free-fall times. Mass exits with a peak flux of $2 \times 10^{-3}$ Msol/yr, showing efficient gas dispersal. The model without radiation pressure has a slight delay in the breakthrough of ionization, but overall its effects are negligible. 85 per cent of the volume, and 40 per cent of the mass, become ionized -- dense filaments resist ionization and are swept up into spherical cores with pillars that point radially away from the ionizing star. We use free-free emission at 20 cm to estimate the production rate of ionizing photons. This is almost always underestimated: by a factor of a few at early stages, then by orders of magnitude as mass leaves the volume. We also test the ratio of dust continuum surface brightnesses at 450 and 850 micron to probe dust temperatures. This underestimates the actual temperature by more than a factor of 2 in areas of low column density or high line-of-sight temperature dispersion; the HII region cavity is particularly prone to this discrepancy. However, the probe is accurate in dense locations such as filaments. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: 16 pages, 17 figures, accepted for publication in MNRAS

Linking Signatures of Accretion with Magnetic Field Measurements - Line Profiles are not Significantly Different in Magnetic and Non-Magnetic Herbig Ae/Be Stars

Authors: Megan Reiter, Nuria Calvet, Thanawuth Thanathibodee, Stefan Kraus, P. Wilson Cauley, John Monnier, Adam Rubinstein, Alicia Aarnio, Tim J. Harries

Abstract: Herbig Ae/Be stars are young, pre-main-sequence stars that sample the transition in structure and evolution between low- and high-mass stars, providing a key test of accretion processes in higher-mass stars. Few Herbig Ae/Be stars have detected magnetic fields, calling into question whether the magnetospheric accretion paradigm developed for low-mass stars can be scaled to higher masses. We presen… ▽ More Herbig Ae/Be stars are young, pre-main-sequence stars that sample the transition in structure and evolution between low- and high-mass stars, providing a key test of accretion processes in higher-mass stars. Few Herbig Ae/Be stars have detected magnetic fields, calling into question whether the magnetospheric accretion paradigm developed for low-mass stars can be scaled to higher masses. We present He I 10830 Å line profiles for 64 Herbig Ae/Be stars with a magnetic field measurement in order to test magnetospheric accretion in the physical regime where its efficacy remains uncertain. Of the 5 stars with a magnetic field detection, 1 shows redshifted absorption, indicative of infall, and 2 show blueshifted absorption, tracing mass outflow. The fraction of redshifted and blueshifted absorption profiles in the non-magnetic Herbig Ae/Be stars is remarkably similar, suggesting that the stellar magnetic field does not affect gas kinematics traced by He I 10830 Å. Line profile morphology does not correlate with the luminosity, rotation rate, mass accretion rate, or disk inclination. Only the detection of a magnetic field and a nearly face-on disk inclination show a correlation (albeit for few sources). This provides further evidence for weaker dipoles and more complex field topologies as stars develop a radiative envelope. The small number of magnetic Herbig Ae/Be stars has already called into question whether magnetospheric accretion can be scaled to higher masses; accretion signatures are not substantially different in magnetic Herbig Ae/Be stars, casting further doubt that they accrete in the same manner as classical T Tauri stars. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Comments: accepted to ApJ; 17 pages, 4 figures, 3 tables

What can the SEDs of first hydrostatic core candidates reveal about their nature?

Authors: Alison K. Young, Matthew R. Bate, Chris F. Mowat, Jennifer Hatchell, Tim J. Harries

Abstract: The first hydrostatic core (FHSC) is the first stable object to form in simulations of star formation. This stage has yet to be observed definitively, although several candidate FHSCs have been reported. We have produced synthetic spectral energy distributions (SEDs) from 3D hydrodynamical simulations of pre-stellar cores undergoing gravitational collapse for a variety of initial conditions. Varia… ▽ More The first hydrostatic core (FHSC) is the first stable object to form in simulations of star formation. This stage has yet to be observed definitively, although several candidate FHSCs have been reported. We have produced synthetic spectral energy distributions (SEDs) from 3D hydrodynamical simulations of pre-stellar cores undergoing gravitational collapse for a variety of initial conditions. Variations in the initial rotation rate, radius and mass lead to differences in the location of the SED peak and far-infrared flux. Secondly, we attempt to fit the SEDs of five FHSC candidates from the literature and five newly identified FHSC candidates located in the Serpens South molecular cloud with simulated SEDs. The most promising FHSC candidates are fitted by a limited number of model SEDs with consistent properties, which suggests the SED can be useful for placing constraints on the age and rotation rate of the source. The sources we consider most likely to be in FHSC phase are B1-bN, CB17-MMS, Aqu-MM1 and Serpens South candidate K242. We were unable to fit SerpS-MM22, Per-Bolo 58 and Chamaeleon-MMS1 with reasonable parameters, which indicates that they are likely to be more evolved. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: 26 pages, 28 figures. Accepted for publication in MNRAS

Radiation-hydrodynamical simulations of massive star formation using Monte Carlo radiative transfer II. The formation of a 25 solar-mass star

Authors: Tim J Harries, Tom A Douglas, Ahmad Ali

Abstract: We present a numerical simulation of the formation of a massive star using Monte- Carlo-based radiation hydrodynamics (RHD). The star forms via stochastic disc accretion and produces fast, radiation-driven bipolar cavities. We find that the evolution of the infall rate (considered to be the mass flux across a 1500 au spherical boundary), and the accretion rate onto the protostar, are broadly consi… ▽ More We present a numerical simulation of the formation of a massive star using Monte- Carlo-based radiation hydrodynamics (RHD). The star forms via stochastic disc accretion and produces fast, radiation-driven bipolar cavities. We find that the evolution of the infall rate (considered to be the mass flux across a 1500 au spherical boundary), and the accretion rate onto the protostar, are broadly consistent with observational constraints. After 35kyr the star has a mass of 25 solar masses and is surrounded by a disc of mass 7 solar masses and 1500 au radius, and we find that the velocity field of the disc is close to Keplerian. Once again these results are consistent with those from recent high-resolution studies of discs around forming massive stars. Synthetic imaging of the RHD model shows good agreement with observations in the near- and far-IR, but may be in conflict with observations that suggests that MYSOs are typically circularly symmetric on the sky at 24.5 microns. Molecular line simulations of a CH3CN transition compare well with observations in terms of surface brightness and line width, and indicate that it should be possible to reliably extract the protostellar mass from such observations. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 11 pages, 10 figures, accepted for publication in MNRAS

Wolf-Rayet spin at low metallicity and its implication for Black Hole formation channels

Authors: Jorick S. Vink, Tim J. Harries

Abstract: The spin of Wolf-Rayet (WR) stars at low metallicity (Z) is most relevant for our understanding of gravitational wave sources such as GW 150914, as well as the incidence of long-duration gamma-ray bursts (GRBs). Two scenarios have been suggested for both phenomena: one of them involves rapid rotation and quasi-chemical homogeneous evolution (CHE), the other invokes classical evolution through mass… ▽ More The spin of Wolf-Rayet (WR) stars at low metallicity (Z) is most relevant for our understanding of gravitational wave sources such as GW 150914, as well as the incidence of long-duration gamma-ray bursts (GRBs). Two scenarios have been suggested for both phenomena: one of them involves rapid rotation and quasi-chemical homogeneous evolution (CHE), the other invokes classical evolution through mass loss in single and binary systems. WR spin rates might enable us to test these two scenarios. In order to obtain empirical constraints on black hole progenitor spin, we infer wind asymmetries in all 12 known WR stars in the Small Magellanic Cloud (SMC) at Z = 1/5 Zsun, as well as within a significantly enlarged sample of single and binary WR stars in the Large Magellanic Cloud (LMC at Z = 1/2 Zsun), tripling the sample of Vink (2007). This brings the total LMC sample to 39, making it appropriate for comparison to the Galactic sample. We measure WR wind asymmetries with VLT-FORS linear spectropolarimetry. We report the detection of new line effects in the LMC WN star BAT99-43 and the WC star BAT99-70, as well as the famous WR/LBV HD 5980 in the SMC, which might be evolving chemically homogeneously. With the previous reported line effects in the late-type WNL (Ofpe/WN9) objects BAT99-22 and BAT99-33, this brings the total LMC WR sample to 4, i.e. a frequency of ~10%. Perhaps surprisingly, the incidence of line effects amongst low-Z WR stars is not found to be any higher than amongst the Galactic WR sample, challenging the rotationally-induced CHE model. As WR mass loss is likely Z-dependent, our Magellanic Cloud line-effect WR stars may maintain their surface rotation and fulfill the basic conditions for producing long GRBs, both via the classical post-red supergiant (RSG) or luminous blue variable (LBV) channel, as well as resulting from CHE due to physics specific to very massive stars (VMS). △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 12 pages, Accepted by Astronomy & Astrophysics

Journal ref: A&A 603, A120 (2017)

Polarized disk emission from Herbig Ae/Be stars observed using Gemini Planet Imager: HD 144432, HD 150193, HD 163296, and HD 169142

Authors: John D. Monnier, T. J. Harries, A. Aarnio, F. Adams, S. Andrews, N. Calvet, C. Espaillat, L. Hartmann, S. Hinkley, S. Kraus, M. McClure, R. Oppenheimer, M. Perrin, D. Wilner

Abstract: In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the… ▽ More In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted Astrophysical Journal

A high-mass protobinary system with spatially resolved circumstellar accretion disks and circumbinary disk

Authors: Stefan Kraus, Jacques Kluska, Alexander Kreplin, Matthew Bate, Tim J. Harries, Karl-Heinz Hofmann, Edward Hone, John D. Monnier, Gerd Weigelt, Narsireddy Anugu, Willem-Jan de Wit, Markus Wittkowski

Abstract: High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9\pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/G… ▽ More High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9\pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/GRAVITY+AMBER near-infrared interferometry allows us to image the circumstellar disks around the individual components with 3 milliarcsecond resolution. We estimate the component masses to $\sim20$ and $\sim18 M_{\odot}$ and find that the radial intensity profiles can be reproduced with an irradiated disk model, where the inner regions are excavated of dust, likely tracing the dust sublimation region in these disks. The circumstellar disks are strongly misaligned with respect to the binary separation vector, which indicates that the tidal forces did not have time to realign the disks, pointing towards a young dynamical age of the system. We constrain the distribution of the Br$γ$ and CO-emitting gas using VLTI/GRAVITY spectro-interferometry and VLT/CRIRES spectro-astrometry and find that the secondary is accreting at a higher rate than the primary. VLT/NACO imaging shows $L'$-band emission on 3-4 times larger scales than the binary separation, matching the expected dynamical truncation radius for the circumbinary disk. The IRAS17216-3801 system is $\sim3\times$ more massive and $\sim5\times$ more compact than other high-mass multiplies imaged at infrared wavelengths and the first high-mass protobinary system where circumstellar and circumbinary dust disks could be spatially resolved. This opens exciting new opportunities for studying star-disk interactions and the role of multiplicity in high-mass star formation. △ Less

Submitted 22 December, 2016; originally announced December 2016.

Comments: 9 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters

Planet Formation Imager (PFI): science vision and key requirements

Authors: Stefan Kraus, John D. Monnier, Michael J. Ireland, Gaspard Duchene, Catherine Espaillat, Sebastian Hoenig, Attila Juhasz, Chris Mordasini, Johan Olofsson, Claudia Paladini, Keivan Stassun, Neal Turner, Gautam Vasisht, Tim J. Harries, Matthew R. Bate, Jean-Francois Gonzalez, Alexis Matter, Zhaohuan Zhu, Olja Panic, Zsolt Regaly, Alessandro Morbidelli, Farzana Meru, Sebastian Wolf, John Ilee, Jean-Philippe Berger , et al. (53 additional authors not shown)

Abstract: The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the prot… ▽ More The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs. △ Less

Submitted 16 August, 2016; v1 submitted 1 August, 2016; originally announced August 2016.

Comments: 13 pages, 5 figures, Proceedings of SPIE 2016

Confronting Standard Models of Proto--Planetary Disks With New Mid--Infrared Sizes from the Keck Interferometer

Authors: Rafael Millan-Gabet, Xiao Che, John D. Monnier, Michael L. Sitko, Ray W. Russell, Carol A. Grady, Amanda N. Day, R. B. Perry, Tim J. Harries, Alicia N. Aarnio, Mark M. Colavita, Peter L. Wizinowich, Sam Ragland, Julien Woillez

Abstract: We present near and mid-infrared interferometric observations made with the Keck Interferometer Nuller and near-contemporaneous spectro-photometry from the IRTF of 11 well known young stellar objects, several observed for the first time in these spectral and spatial resolution regimes. With AU-level spatial resolution, we first establish characteristic sizes of the infrared emission using a simple… ▽ More We present near and mid-infrared interferometric observations made with the Keck Interferometer Nuller and near-contemporaneous spectro-photometry from the IRTF of 11 well known young stellar objects, several observed for the first time in these spectral and spatial resolution regimes. With AU-level spatial resolution, we first establish characteristic sizes of the infrared emission using a simple geometrical model consisting of a hot inner rim and mid-infrared disk emission. We find a high degree of correlation between the stellar luminosity and the mid-infrared disk sizes after using near-infrared data to remove the contribution from the inner rim. We then use a semi-analytical physical model to also find that the very widely used "star + inner dust rim + flared disk" class of models strongly fails to reproduce the SED and spatially-resolved mid-infrared data simultaneously; specifically a more compact source of mid-infrared emission is required than results from the standard flared disk model. We explore the viability of a modification to the model whereby a second dust rim containing smaller dust grains is added, and find that the two-rim model leads to significantly improved fits in most cases. This complexity is largely missed when carrying out SED modelling alone, although detailed silicate feature fitting by McClure et al. 2013 recently came to a similar conclusion. As has been suggested recently by Menu et al. 2015, the difficulty in predicting mid-infrared sizes from the SED alone might hint at "transition disk"-like gaps in the inner AU; however, the relatively high correlation found in our mid-infrared disk size vs. stellar luminosity relation favors layered disk morphologies and points to missing disk model ingredients instead. △ Less

Submitted 22 April, 2016; originally announced April 2016.

Directly observing continuum emission from self-gravitating spiral waves

Authors: Cassandra Hall, Duncan Forgan, Ken Rice, Tim J. Harries, Pamela D. Klaassen, Beth Biller

Abstract: We use a simple, self-consistent, self-gravitating semi-analytic disc model to conduct an examination of the parameter space in which self-gravitating discs may exist. We then use Monte-Carlo radiative transfer to generate synthetic ALMA images of these self-gravitating discs to determine the subset of this parameter space in which they generate non-axisymmetric structure that is potentially detec… ▽ More We use a simple, self-consistent, self-gravitating semi-analytic disc model to conduct an examination of the parameter space in which self-gravitating discs may exist. We then use Monte-Carlo radiative transfer to generate synthetic ALMA images of these self-gravitating discs to determine the subset of this parameter space in which they generate non-axisymmetric structure that is potentially detectable by ALMA. Recently, several transition discs have been observed to have non-axisymmetric structure that extends out to large radii. It has been suggested that one possible origin of these asymmetries could be spiral density waves induced by disc self-gravity. We use our simple model to see if these discs exist in the region of parameter space where self-gravity could feasibly explain these spiral features. We find that for self-gravity to play a role in these systems typically requires a disc mass around an order of magnitude higher than the observed disc masses for the systems. The spiral amplitudes produced by self-gravity in the local approximation are relatively weak when compared to amplitudes produced by tidal interactions, or spirals launched at Lindblad resonances due to embedded planets in the disc. As such, we ultimately caution against diagnosing spiral features as being due to self-gravity, unless the disc exists in the very narrow region of parameter space where the spiral wave amplitudes are large enough to produce detectable features, but not so large as to cause the disc to fragment. △ Less

Submitted 3 February, 2016; originally announced February 2016.

Comments: 13 pages, 9 figures

Exploring the Origins of Deuterium Enrichments in Solar Nebular Organics

Authors: L. Ilsedore Cleeves, Edwin A. Bergin, Conel M. O'D. Alexander, Fujun Du, Dawn Graninger, Karin I. Öberg, Tim J. Harries

Abstract: Deuterium-to-hydrogen (D/H) enrichments in molecular species provide clues about their original formation environment. The organic materials in primitive solar system bodies have generally higher D/H ratios and show greater D/H variation when compared to D/H in solar system water. We propose this difference arises at least in part due to 1) the availability of additional chemical fractionation pat… ▽ More Deuterium-to-hydrogen (D/H) enrichments in molecular species provide clues about their original formation environment. The organic materials in primitive solar system bodies have generally higher D/H ratios and show greater D/H variation when compared to D/H in solar system water. We propose this difference arises at least in part due to 1) the availability of additional chemical fractionation pathways for organics beyond that for water, and 2) the higher volatility of key carbon reservoirs compared to oxygen. We test this hypothesis using detailed disk models, including a sophisticated, new disk ionization treatment with a low cosmic ray ionization rate, and find that disk chemistry leads to higher deuterium enrichment in organics compared to water, helped especially by fractionation via the precursors CH$_2$D$^+$/CH$_3^+$. We also find that the D/H ratio in individual species varies significantly depending on their particular formation pathways. For example, from $\sim20-40$ AU, CH$_4$ can reach $\rm{D/H\sim2\times10^{-3}}$, while D/H in CH$_3$OH remains locally unaltered. Finally, while the global organic D/H in our models can reproduce intermediately elevated D/H in the bulk hydrocarbon reservoir, our models are unable to reproduce the most deuterium-enriched organic materials in the solar system, and thus our model requires some inheritance from the cold interstellar medium from which the Sun formed. △ Less

Submitted 21 January, 2016; originally announced January 2016.

Comments: 11 pages, 7 figures, accepted for publication in ApJ

Detecting stellar-wind bubbles through infrared arcs in HII regions

Authors: Jonathan Mackey, Thomas J. Haworth, Vasilii V. Gvaramadze, Shazrene Mohamed, Norbert Langer, Tim J. Harries

Abstract: 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 regio… ▽ More 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. △ Less

Submitted 21 December, 2015; originally announced December 2015.

Comments: 16 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 586, A114 (2016)

Modelling multi-wavelength observational characteristics of bow shocks from runaway early type stars

Authors: David M. Acreman, Ian R. Stevens, Tim J. Harries

Abstract: We assess the multi-wavelength observable properties of the bow shock around a runaway early type star using a combination of hydrodynamical modelling, radiative transfer calculations and synthetic imaging. Instabilities associated with the forward shock produce dense knots of material which are warm, ionised and contain dust. These knots of material are responsible for the majority of emission at… ▽ More We assess the multi-wavelength observable properties of the bow shock around a runaway early type star using a combination of hydrodynamical modelling, radiative transfer calculations and synthetic imaging. Instabilities associated with the forward shock produce dense knots of material which are warm, ionised and contain dust. These knots of material are responsible for the majority of emission at far infra-red, H alpha and radio wavelengths. The large scale bow shock morphology is very similar and differences are primarily due to variations in the assumed spatial resolution. However infra-red intensity slices (at 22 microns and 12 microns) show that the effects of a temperature gradient can be resolved at a realistic spatial resolution for an object at a distance of 1 kpc. △ Less

Submitted 10 November, 2015; originally announced November 2015.

Comments: 11 pages, 6 figures. Accepted for publication in MNRAS

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

Authors: T. G. Bisbas, T. J. Haworth, M. J. Barlow, S. Viti, T. J. Harries, T. Bell, J. A. Yates

Abstract: 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 eith… ▽ More 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. △ Less

Submitted 16 September, 2015; originally announced September 2015.

Comments: 17 pages, 10 Figures, accepted for publication in MNRAS. Comments welcome

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

Authors: Thomas J. Haworth, Tim J. Harries, David M. Acreman, Thomas G. Bisbas

Abstract: 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-… ▽ More 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. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: 14 pages, 9 figures, accepted for publication in the MNRAS main journal

Indirect Detection of Forming Protoplanets via Chemical Asymmetries in Disks

Authors: L. Ilsedore Cleeves, Edwin A. Bergin, Tim J. Harries

Abstract: We examine changes in the molecular abundances resulting from increased heating due to a self-luminous planetary companion embedded within a narrow circumstellar disk gap. Using 3D models that include stellar and planetary irradiation, we find that luminous young planets locally heat up the parent circumstellar disk by many tens of Kelvin, resulting in efficient thermal desorption of molecular spe… ▽ More We examine changes in the molecular abundances resulting from increased heating due to a self-luminous planetary companion embedded within a narrow circumstellar disk gap. Using 3D models that include stellar and planetary irradiation, we find that luminous young planets locally heat up the parent circumstellar disk by many tens of Kelvin, resulting in efficient thermal desorption of molecular species that are otherwise locally frozen out. Furthermore, the heating is deposited over large regions of the disk, $\pm5$ AU radially and spanning $\lesssim60^\circ$ azimuthally. From the 3D chemical models, we compute rotational line emission models and full ALMA simulations, and find that the chemical signatures of the young planet are detectable as chemical asymmetries in $\sim10h$ observations. HCN and its isotopologues are particularly clear tracers of planetary heating for the models considered here, and emission from multiple transitions of the same species is detectable, which encodes temperature information in addition to possible velocity information from the spectra itself. We find submillimeter molecular emission will be a useful tool to study gas giant planet formation in situ, especially beyond $R\gtrsim10$ AU. △ Less

Submitted 27 May, 2015; originally announced May 2015.

Comments: 14 pages, 14 figures, accepted for publication in ApJ

Radiation-hydrodynamical simulations of massive star formation using Monte Carlo radiative transfer: I. Algorithms and numerical methods

Authors: Tim J. Harries

Abstract: We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo estimator and enables a complete treatment of the detailed microphysics, including polychromatic radiation and anisotropic scattering, in both the free-streaming and… ▽ More We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo estimator and enables a complete treatment of the detailed microphysics, including polychromatic radiation and anisotropic scattering, in both the free-streaming and optically-thick limits. Since the new method is computationally demanding we have developed two new methods that speed up the algorithm. The first is a photon packet splitting algorithm that enables efficient treatment of the Monte Carlo process in very optically thick regions. The second is a parallelisation method that distributes the Monte Carlo workload over many instances of the hydrodynamic domain, resulting in excellent scaling of the radiation step. We also describe the implementation of a sink particle method that enables us to follow the accretion onto, and the growth of, the protostars. We detail the results of extensive testing and benchmarking of the new algorithms. △ Less

Submitted 23 January, 2015; originally announced January 2015.

Comments: 12 pages, 15 figures. Accepted by MNRAS

The ancient heritage of water ice in the solar system

Authors: L. Ilsedore Cleeves, Edwin A. Bergin, Conel M. O'D. Alexander, Fujun Du, Dawn Graninger, Karin I. Öberg, Tim J. Harries

Abstract: Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk.… ▽ More Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk. Utilizing a comprehensive treatment of disk ionization, we find that ion-driven deuterium pathways are inefficient, curtailing the disk's deuterated water formation and its viability as the sole source for the solar system's water. This finding implies that if the solar system's formation was typical, abundant interstellar ices are available to all nascent planetary systems. △ Less

Submitted 25 September, 2014; originally announced September 2014.

Comments: 33 pages, 7 figures including main text and supplementary materials. Published in Science

Journal ref: Science, Vol. 345 no. 6204 pp. 1590-1593 (2014)

Line and continuum radiative transfer modelling of AA Tau

Authors: C. F. Esau, T. J. Harries, J. Bouvier

Abstract: We present photometric and spectroscopic models of the Classical T Tauri star AA Tau. Photometric and spectroscopic variability present in observations of AA Tau is attributed to a magnetically induced warp in the accretion disc, periodically occulting the photosphere on an 8.2--day timescale. Emission line profiles show signatures of both infall, attributed to magnetospherically accreting materia… ▽ More We present photometric and spectroscopic models of the Classical T Tauri star AA Tau. Photometric and spectroscopic variability present in observations of AA Tau is attributed to a magnetically induced warp in the accretion disc, periodically occulting the photosphere on an 8.2--day timescale. Emission line profiles show signatures of both infall, attributed to magnetospherically accreting material, and outflow. Using the radiative transfer code TORUS, we have investigated the geometry and kinematics of AA Tau's circumstellar disc and outflow, which is modelled here as a disc wind. Photometric models have been used to constrain the aspect ratio of the disc, the offset angle of the magnetosphere dipole with respect to the stellar rotation axis, and the inner radius of the circumstellar disc. Spectroscopic models have been used to constrain the wind and magnetosphere temperatures, wind acceleration parameter, and mass loss rate. We find observations are best fitted by models with a mass accretion rate of $5\times10^{-9}$ M$_\odot$ yr$^{-1}$, a dipole offset of between $10^\circ$ and $20^\circ$, a magnetosphere that truncates the disc from 5.2 to 8.8 R$_\star$, a mass-loss-rate to accretion-rate ratio of ~ 0.1, a magnetosphere temperature of 8500 -- 9000 K, and a disc wind temperature of 8000 K. △ Less

Submitted 14 November, 2014; v1 submitted 16 June, 2014; originally announced June 2014.

Comments: 22 pages, 32 figures, 4 tables. Accepted by MNRAS. V3: Corrected typo

Resolving the gap and AU-scale asymmetries in the pre-transitional disk of V1247 Orionis

Authors: Stefan Kraus, Michael J. Ireland, Michael L. Sitko, John D. Monnier, Nuria Calvet, Catherine Espaillat, Carol A. Grady, Tim J. Harries, Sebastian F. Hoenig, Ray W. Russell, Jeremy R. Swearingen, Chelsea Werren, David J. Wilner

Abstract: Pre-transitional disks are protoplanetary disks with a gapped disk structure, potentially indicating the presence of young planets in these systems. In order to explore the structure of these objects and their gap-opening mechanism, we observed the pre-transitional disk V1247 Orionis using the Very Large Telescope Interferometer, the Keck Interferometer, Keck-II, Gemini South, and IRTF. This allow… ▽ More Pre-transitional disks are protoplanetary disks with a gapped disk structure, potentially indicating the presence of young planets in these systems. In order to explore the structure of these objects and their gap-opening mechanism, we observed the pre-transitional disk V1247 Orionis using the Very Large Telescope Interferometer, the Keck Interferometer, Keck-II, Gemini South, and IRTF. This allows us spatially resolve the AU-scale disk structure from near- to mid-infrared wavelengths (1.5 to 13 μm), tracing material at different temperatures and over a wide range of stellocentric radii. Our observations reveal a narrow, optically-thick inner-disk component (located at 0.18 AU from the star) that is separated from the optically thick outer disk (radii >46 AU), providing unambiguous evidence for the existence of a gap in this pre-transitional disk. Surprisingly, we find that the gap region is filled with significant amounts of optically thin material with a carbon-dominated dust mineralogy. The presence of this optically thin gap material cannot be deduced solely from the spectral energy distribution, yet it is the dominant contributor at mid-infrared wavelengths. Furthermore, using Keck/NIRC2 aperture masking observations in the H, K', and L' band, we detect asymmetries in the brightness distribution on scales of about 15-40 AU, i.e. within the gap region. The detected asymmetries are highly significant, yet their amplitude and direction changes with wavelength, which is not consistent with a companion interpretation but indicates an inhomogeneous distribution of the gap material. We interpret this as strong evidence for the presence of complex density structures, possibly reflecting the dynamical interaction of the disk material with sub-stellar mass bodies that are responsible for the gap clearing. △ Less

Submitted 9 April, 2013; originally announced April 2013.

Comments: 16 pages, 17 Figures, accepted by Astrophysical Journal

Assessing molecular line diagnostics of triggered star formation using synthetic observations

Authors: Thomas J. Haworth, Tim J. Harries, David M. Acreman, David A. Rundle

Abstract: 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 d… ▽ More 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. △ Less

Submitted 21 March, 2013; v1 submitted 6 March, 2013; originally announced March 2013.

Comments: 17 pages, Accepted for publication in MNRAS, revision following editor proof

Rotating Wolf-Rayet stars in a post RSG/LBV phase. An evolutionary channel towards long-duration GRBs?

Authors: G. Gräfener, J. S. Vink, T. J. Harries, N. Langer

Abstract: (shortened) We investigate the properties of Galactic WR stars and their environment to identify evolutionary channels that may lead to the formation of LGRBs. To this purpose we compile available information on the spectropolarimetric properties, the presence of CS ejecta, and the CS velocities in the environment of Galactic WR stars. We use linear line-depolarization as an indicator of rotation,… ▽ More (shortened) We investigate the properties of Galactic WR stars and their environment to identify evolutionary channels that may lead to the formation of LGRBs. To this purpose we compile available information on the spectropolarimetric properties, the presence of CS ejecta, and the CS velocities in the environment of Galactic WR stars. We use linear line-depolarization as an indicator of rotation, nebular morphology as an indicator of stellar ejecta, and velocity patterns in UV absorption features as an indicator of increased velocities in the CS environment. We find that the ~23% WR stars with "possible ejecta nebulae" dominate the population of WR stars with spectropolarimetric signatures of rotation, while WR stars without such nebulae only rarely show indications of rotation. The corresponding objects are most likely in an early stage after a preceding RSG or LBV phase, and have not yet lost their angular momenta due to the strong mass-loss in the WR phase. From their photometric periods we estimate rotation parameters in the range omega = 0.04...0.25, corresponding to moderate rotation speeds of 36...120 km/s. These values are very uncertain, but comply with the specific surface angular momentum requirement for LGRB progenitors. Our results indicate that, in the Galaxy, mainly "young" WR stars shortly after a RSG/LBV phase show spectropolarimetric signatures of rotation. Their rotation rates are thus likely enhanced with respect to the majority of Galactic WR stars. According to their estimated specific surface angular momenta, a subgroup of stars exploding in this phase may represent an evolutionary channel towards LGRBs at high metallicities, comparable to the Galaxy. △ Less

Submitted 12 October, 2012; v1 submitted 3 October, 2012; originally announced October 2012.

Comments: 13 pages, 2 figures, A&A accepted

Radiation hydrodynamics simulations of massive star formation using Monte Carlo radiation transfer

Authors: Tim J. Harries, Tom J. Haworth, David Acreman

Abstract: We present a radiation hydrodynamics simulation of the formation of a massive star using a Monte Carlo treatment for the radiation field. We find that strong, high speed bipolar cavities are driven by the radiation from the protostar, and that accretion occurs stochastically from a circumstellar disc. We have computed spectral energy distributions and images at each timestep, which may in future b… ▽ More We present a radiation hydrodynamics simulation of the formation of a massive star using a Monte Carlo treatment for the radiation field. We find that strong, high speed bipolar cavities are driven by the radiation from the protostar, and that accretion occurs stochastically from a circumstellar disc. We have computed spectral energy distributions and images at each timestep, which may in future be used to compare our models with photometric, spectroscopic, and interferometric observations of young massive stellar objects. △ Less

Submitted 7 September, 2012; originally announced September 2012.

Comments: 4 pages, 2 Figures; to appear in The Labyrinth of Star Formation, (eds.) D. Stamatellos, S. Goodwin, and D. Ward-Thompson, Springer, in press

Testing models of triggered star formation: theory and observation

Authors: Thomas J. Haworth, Tim J. Harries, David M. Acreman

Abstract: 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 ra… ▽ More 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. △ Less

Submitted 16 August, 2012; v1 submitted 8 August, 2012; originally announced August 2012.

Comments: 2 pages, 1 figure; to appear in The Labyrinth of Star Formation, (eds.) D. Stamatellos, S. Goodwin, and D. Ward-Thompson, Springer, in press. Citation issues resolved

Testing diagnostics of triggered star formation

Authors: Thomas J. Haworth, Tim J. Harries, David M. Acreman

Abstract: 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 la… ▽ More 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. △ Less

Submitted 10 July, 2012; v1 submitted 31 May, 2012; originally announced May 2012.

Comments: 17 pages. Post reviewer comments

Bayesian fitting of Taurus brown dwarf spectral energy distributions

Authors: N. J. Mayne, Tim J. Harries, John Rowe, David M. Acreman

Abstract: We present derived stellar and disc parameters for a sample of Taurus brown dwarfs both with and without evidence of an associated disc. These parameters have been derived using an online fitting tool (http://bd-server.astro.ex.ac.uk/), which includes a statistically robust derivation of uncertainties, an indication of pa- rameter degeneracies, and a complete treatment of the input photometric and… ▽ More We present derived stellar and disc parameters for a sample of Taurus brown dwarfs both with and without evidence of an associated disc. These parameters have been derived using an online fitting tool (http://bd-server.astro.ex.ac.uk/), which includes a statistically robust derivation of uncertainties, an indication of pa- rameter degeneracies, and a complete treatment of the input photometric and spectroscopic observations. The observations of the Taurus members with indications of disc presence have been fitted using a grid of theoretical models including detailed treatments of physical processes accepted for higher mass stars, such as dust sublimation, and a simple treatment of the accretion flux. This grid of models has been designed to test the validity of the adopted physical mechanisms, but we have also constructed models using parameterisation, for example semi-empirical dust sublimation radii, for users solely interested in parameter derivation and the quality of the fit. The parameters derived for the naked and disc brown dwarf systems are largely consistent with literature observations. However, our inner disc edge locations are consistently closer to the star than previous results and we also derive elevated accretion rates over non-SED based accretion rate derivations. For inner edge locations we attribute these differences to the detailed modelling we have performed of the disc structure, particularly at the crucial inner edge where departures in geometry from the often adopted vertical wall due to dust sublimation (and therefore accretion flux) can compensate for temperature (and therefore distance) changes to the inner edge of the dust disc. In the case of the elevated derived accretion rates, in some cases, this may be caused by the intrinsic stellar luminosities of the targets exceeding that predicted by the isochrones we have adopted. △ Less

Submitted 29 March, 2012; originally announced March 2012.

Comments: The paper contains 35 pages with 15 figures and 17 tables. Accepted for publication in MNRAS

In pursuit of gamma-ray burst progenitors: the identification of a sub-population of rotating Wolf-Rayet stars

Authors: Jorick S. Vink, G. Graefener, T. J. Harries

Abstract: Long gamma-ray bursts involve the most powerful cosmic explosions since the Big Bang. Whilst it has been established that GRBs are related to the death throes of massive stars, the identification of their progenitors has proved challenging. Theory suggests that rotating Wolf-Rayet stars are the best candidates, but their strong stellar winds shroud their surfaces, preventing a direct measurement o… ▽ More Long gamma-ray bursts involve the most powerful cosmic explosions since the Big Bang. Whilst it has been established that GRBs are related to the death throes of massive stars, the identification of their progenitors has proved challenging. Theory suggests that rotating Wolf-Rayet stars are the best candidates, but their strong stellar winds shroud their surfaces, preventing a direct measurement of their rotation. Fortunately, linear spectropolarimetry may be used to probe the flattening of their winds due to stellar spin. Spectropolarimetry surveys show that an 80% majority of WR stars have spherically symmetric winds and are thus rotating slowly, yet a small 20% minority display a spectropolarimetric signature indicative of rotation. Here we find a highly significant correlation between WR objects that carry the signature of stellar rotation and the subset of WR stars with ejecta nebulae that have only recently transitioned from a red sugergiant or luminous blue variable phase. As these youthful WR stars have yet to spin-down due to mass loss, they are the best candidate GRB progenitors identified to date. When we take recently published WR ejecta nebula numbers we find that five out of the six line-effect WR stars are surrounded by ejecta nebulae. The statistics imply that the null hypothesis of no correlation between line-effect WR stars and ejecta nebulae can be rejected at the 0.0004% level. Given that four line-effect and WR ejecta nebula have spectroscopically been confirmed to contain nucleosynthetic products, we argue that the correlation is both statistically significant and physically convincing. The implication is that we have identified a WR sub-population that fulfills the necessary criteria for making GRBs. Finally, we discuss the potential of identifying GRB progenitors via spectropolarimetry with extremely large telescopes. △ Less

Submitted 5 December, 2011; v1 submitted 24 November, 2011; originally announced November 2011.

Comments: 5 pages, accepted for publication in Astronomy & Astrophysics Letters (small textual changes)

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

Authors: Thomas J. Haworth, Tim J. Harries

Abstract: 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 natu… ▽ More 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. △ Less

Submitted 23 November, 2011; v1 submitted 24 October, 2011; originally announced October 2011.

Comments: Accepted for publication in MNRAS. 19 pages, 14 figures, 8 tables

An algorithm for Monte-Carlo time-dependent radiation transfer

Authors: Tim J. Harries

Abstract: A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick and free-streaming regimes simultaneously. The algorithm is tested on a variety of 1-d and 2-d problems, and good agreement with benchmark solutions is found.… ▽ More A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick and free-streaming regimes simultaneously. The algorithm is tested on a variety of 1-d and 2-d problems, and good agreement with benchmark solutions is found. The method is used to calculate the time-varying spectral energy distribution from a circumstellar disc illuminated by a protostar whose accretion luminosity is varying. It is shown that the time lag between the optical variability and the infrared variability results from a combination of the photon travel time and the thermal response in the disc, and that the lag is an approximately linear function of wavelength. △ Less

Submitted 27 May, 2011; originally announced May 2011.

Comments: 10 pages, 10 figures, accepted for publication by MNRAS

Multidimensional models of hydrogen and helium emission line profiles for classical T Tauri Stars: method, tests and examples

Authors: Ryuichi Kurosawa, Marina M. Romanova, Tim J. Harries

Abstract: We present multidimensional non-LTE radiative transfer models of hydrogen and helium line profiles formed in the accretion flows and the outflows near the star-disk interaction regions of classical T Tauri stars (CTTSs). The statistical equilibrium calculations, performed under the assumption of the Sobolev approximation using the radiative transfer code TORUS, has been improved to include He I an… ▽ More We present multidimensional non-LTE radiative transfer models of hydrogen and helium line profiles formed in the accretion flows and the outflows near the star-disk interaction regions of classical T Tauri stars (CTTSs). The statistical equilibrium calculations, performed under the assumption of the Sobolev approximation using the radiative transfer code TORUS, has been improved to include He I and He II energy levels. This allows us to probe the physical conditions of the inner wind of CTTSs by simultaneously modelling the robust wind diagnostic line He I (10830) and the accretion diagnostic lines such as Pa-beta, Br-gamma and He I (5876). The code has been tested in 1 and 2-D problems, and we have shown that the results are in agreement with established codes. We apply the model to the complex flow geometries of CTTSs. Example model profiles are computed using the combinations of (1) magnetospheric accretion and disc wind, and (2) magnetospheric accretion and the stellar wind. In both cases, the model produces line profiles which are qualitatively similar to those found in observations. Our models are consistent with the scenario in which the narrow blueshifted absorption component of He I (10830) seen in observations is caused by a disc wind, and the wider blueshifted absorption component (the P-Cygni profile) is caused by a bipolar stellar wind. However, we do not have a strong constraint on the relative importance of the wind and the magnetosphere for the `emission' component. Our preliminary calculations suggest that the temperatures of the disc wind, stellar wind and the magnetosphere cannot be much higher than ~10,000 K, on the basis of the strengths of hydrogen lines. With these low temperatures, we find that the photoionzation by high energy photons (e.g. X-rays) is necessary to produce He I (10830) in emission and to produce the blueshifted absorption components. △ Less

Submitted 9 June, 2011; v1 submitted 3 February, 2011; originally announced February 2011.

Comments: 18 pages, 12 figures, accepted for publication in MNRAS

Three-dimensional molecular line transfer: A simulated star-forming region

Authors: David Rundle, Tim J. Harries, David M. Acreman, Matthew R. Bate

Abstract: We present the first non-LTE, co-moving frame molecular line calculations of a star-forming cluster simulated using smoothed particle hydrodynamics (SPH), from which we derive high-resolution synthetic observations. We have resampled a particle representation onto an adaptive mesh and self-consistently solved the equations of statistical equilibrium in the co-moving frame, using TORUS, a three-dim… ▽ More We present the first non-LTE, co-moving frame molecular line calculations of a star-forming cluster simulated using smoothed particle hydrodynamics (SPH), from which we derive high-resolution synthetic observations. We have resampled a particle representation onto an adaptive mesh and self-consistently solved the equations of statistical equilibrium in the co-moving frame, using TORUS, a three-dimensional adaptive mesh refined (AMR) radiative transfer (RT) code. We verified the applicability of the code to the conditions of the SPH simulation by testing its output against other codes. We find that the level populations obtained for optically thick and thin scenarios closely match the ensemble average of the other codes. We have used the code to obtain non-LTE level populations of multiple molecular species throughout the cluster and have created three-dimensional velocity-resolved spatial maps of the emergent intensity. Line profiles of cores traced by N2H+ (1-0) are compared to probes of low density gas, 13CO (1-0) and C18O (1-0), surrounding the cores along the line of sight. The relative differences of the line-centre velocities are shown to be small compared to the velocity dispersion, matching recent observations. We conclude that one cannot reject competitive accretion as a viable theory of star formation based on observed velocity profiles. △ Less

Submitted 10 May, 2010; originally announced May 2010.

Comments: 19 pages, 17 figures, 3 tables

Modelling circumstellar discs with 3D radiation hydrodynamics

Authors: David M. Acreman, Tim J. Harries, David A. Rundle

Abstract: We present results from combining a grid-based radiative transfer code with a Smoothed Particle Hydrodynamics code to produce a flexible system for modelling radiation hydrodynamics. We use a benchmark model of a circumstellar disc to determine a robust method for constructing a gridded density distribution from SPH particles. The benchmark disc is then used to determine the accuracy of the radi… ▽ More We present results from combining a grid-based radiative transfer code with a Smoothed Particle Hydrodynamics code to produce a flexible system for modelling radiation hydrodynamics. We use a benchmark model of a circumstellar disc to determine a robust method for constructing a gridded density distribution from SPH particles. The benchmark disc is then used to determine the accuracy of the radiative transfer results. We find that the SED and the temperature distribution within the disc are sensitive to the representation of the disc inner edge, which depends critically on both the grid and SPH resolution. The code is then used to model a circumstellar disc around a T-Tauri star. As the disc adjusts towards equilibrium vertical motions in the disc are induced resulting in scale height enhancements which intercept radiation from the central star. Vertical transport of radiation enables these perturbations to influence the mid-plane temperature of the disc. The vertical motions decay over time and the disc ultimately reaches a state of simultaneous hydrostatic and radiative equilibrium. △ Less

Submitted 10 December, 2009; originally announced December 2009.

Comments: MNRAS accepted; 15 pages; 17 figures, 4 in colour