sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Alexander, Richard

Normalized to: Alexander, R.

61 article(s) in total. 210 co-authors, from 1 to 11 common article(s). Median position in authors list is 3,0.

[1]  oai:arXiv.org:2006.09811  [pdf] - 2120087
Forbidden line diagnostics of photoevaporative disc winds
Comments: 14 pages, 10 figures. Accepted for publication in MNRAS
Submitted: 2020-06-17
Photoevaporation driven by high energy radiation from the central star plays an important role in the evolution of protoplanetary discs. Photoevaporative winds have been unambiguously detected through blue-shifted emission lines, but their detailed properties remain uncertain. Here we present a new empirical approach to make observational predictions of these thermal winds, seeking to fill the gap between theory and observations. We use a self-similar model of an isothermal wind to compute line profiles of several characteristic emission lines (the [Ne${\rm{\scriptsize II}}$] line at 12.81 $\mu$m, and optical forbidden lines such as [O${\rm{\scriptsize I}}$] 6300 $\mathring{A}$ and [S${\rm{\scriptsize II}}$] 4068/4076 $\mathring{A}$), studying how the lines are affected by parameters such as the gas temperature, disc inclinations, and density profile. Our model successfully reproduces blue-shifted lines with $v_{\rm peak} \lesssim 10$ km/s, which decrease with increasing disc inclination. The line widths increase with increasing disc inclinations and range from $\Delta v \sim 15-30$ km/s. The predicted blue-shifts are mostly sensitive to the gas sound speed. The observed [Ne${\rm{\scriptsize II}}$] line profiles are consistent with a thermal wind and point towards a relatively high sound speed, as expected for EUV photoevaporation. However, the observed [O${\rm{\scriptsize I}}$] line profiles require lower temperatures, as expected in X-ray photoevaporation, and show a wider scatter that is difficult to reconcile with a single wind model; it seems likely that these lines trace different components of a multi-phase wind. We also note that the spectral resolution of current observations remains an important limiting factor in these studies, and that higher resolution spectra are required if emission lines are to further our understanding of protoplanetary disc winds.
[2]  oai:arXiv.org:2003.14314  [pdf] - 2093268
NGTS J214358.5-380102 -- NGTS discovery of the most eccentric known eclipsing M-Dwarf binary system
Comments: Accepted for Publication in MNRAS
Submitted: 2020-03-31, last modified: 2020-05-11
We present the discovery of NGTS J214358.5-380102, an eccentric M-dwarf binary discovered by the Next Generation Transit Survey. The system period of 7.618 days is greater than many known eclipsing M-dwarf binary systems. Its orbital eccentricity of $0.323^{+0.0014}_{-0.0037}$, is large relative to the period and semi-major axis of the binary. Global modelling of photometry and radial velocities indicate stellar masses of $M_A$=$0.426 ^{+0.0056}_{-0.0049}$, $M_B$=$0.455 ^{+0.0058}_{-0.0052}$ and stellar radii $R_A$=$0.461 ^{+0.038}_{-0.025}$ $R_B$=$0.411 ^{+0.027}_{-0.039}$, respectively. Comparisons with stellar models for low mass stars show that one star is consistent with model predictions whereas the other is substantially oversized. Spectral analysis of the system suggests a primary of spectral type M3V, consistent with both modelled masses and radii, and with SED fitting of NGTS photometry. As the most eccentric eclipsing M-dwarf binary known, NGTS J214358.5-380102 provides an interesting insight into the strength of tidal effects in the circularisation of stellar orbits.
[3]  oai:arXiv.org:1911.05760  [pdf] - 2026342
Flyby-induced misalignments in planet-hosting discs
Comments: Accepted for publication in MNRAS. 9 pages with 6 figures
Submitted: 2019-11-13
We now have several observational examples of misaligned broken protoplanetary discs, where the disc inner regions are strongly misaligned with respect to the outer disc. Current models suggest that this disc structure can be generated with an internal misaligned companion (stellar or planetary), but the occurrence rate of these currently unobserved companions remains unknown. Here we explore whether a strong misalignment between the inner and outer disc can be formed without such a companion. We consider a disc that has an existing gap --- essentially separating the disc into two regions --- and use a flyby to disturb the discs, leading to a misalignment. Despite considering the most optimistic parameters for this scenario, we find maximum misalignments between the inner and outer disc of $\sim$45$^{\circ}$ and that these misalignments are short-lived. We thus conclude that the currently observed misaligned discs must harbour internal, misaligned companions.
[4]  oai:arXiv.org:1911.00297  [pdf] - 2026223
A Swift view of X-ray and UV radiation in the planet-forming T-Tauri system PDS 70
Comments: 5 pages, 3 figures, accepted for publication in MNRAS Letters
Submitted: 2019-11-01
PDS 70 is a $\sim$5 Myr old star with a gas and dust disc in which several proto-planets have been discovered. We present the first UV detection of the system along with X-ray observations taken with the \textit{Neil Gehrels Swift Observatory} satellite. PDS 70 has an X-ray flux of 3.4$\times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 0.3-10.0 keV range, and UV flux (U band) of 3.5$\times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ . At the distance of 113.4 pc determined from Gaia DR2 this gives luminosities of 5.2$\times 10^{29}$ erg s$^{-1}$ and 5.4$\times 10^{29}$ erg s$^{-1}$ respectively. The X-ray luminosity is consistent with coronal emission from a rapidly rotating star close to the log $\frac{L_{\mathrm{X}}}{L_{\mathrm{bol}}} \sim -3$ saturation limit. We find the UV luminosity is much lower than would be expected if the star were still accreting disc material and suggest that the observed UV emission is coronal in origin.
[5]  oai:arXiv.org:1910.06822  [pdf] - 1994321
Flybys in protoplanetary discs: II. Observational signatures
Comments: 11 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2019-10-15
Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa).
[6]  oai:arXiv.org:1909.04093  [pdf] - 1967076
A solution to the proplyd lifetime problem
Comments: 16 pages, 11 figures, accepted for publication in MNRAS
Submitted: 2019-09-09
Protoplanetary discs (PPDs) in the Orion Nebula Cluster (ONC) are irradiated by UV fields from the massive star $\theta^1$C. This drives thermal winds, inducing mass loss rates of up to $\dot{M}_\mathrm{wind}\sim 10^{-7}\,M_\odot$/yr in the `proplyds' (ionised PPDs) close to the centre. For the mean age of the ONC and reasonable initial PPD masses, such mass loss rates imply that discs should have been dispersed. However, ~80% of stars still exhibit a NIR excess, suggesting that significant circumstellar mass remains. This `proplyd lifetime problem' has persisted since the discovery of photoevaporating discs in the core of the ONC by O'Dell & Wen (1994). In this work, we demonstrate how an extended period of star formation can solve this problem. Coupling N-body calculations and a viscous disc evolution model, we obtain high disc fractions at the present day. This is partly due to the migration of older stars outwards, and younger stars inwards such that the most strongly irradiated PPDs are also the youngest. We show how the disc mass distribution can be used to test the recent claims in the literature for multiple stellar populations in the ONC. Our model also explains the recent finding that host mass and PPD mass are only weakly correlated, in contrast with other regions of similar age. We conclude that the status of the ONC as the archetype for understanding the influence of environment on planet formation is undeserved; the complex star formation history (involving star formation episodes within ~0.8 Myr of the present day) results in confusing signatures in the PPD population.
[7]  oai:arXiv.org:1901.02407  [pdf] - 1897015
The temporal requirements of directly observing self-gravitating spiral waves in protoplanetary discs with ALMA
Comments: 13 pages, 6 figures
Submitted: 2019-01-08, last modified: 2019-06-08
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
[8]  oai:arXiv.org:1902.00036  [pdf] - 1830688
Scattered light shadows in warped protoplanetary discs
Comments: Accepted for publication in MNRAS, 13 pages and 9 figures. Movies available at https://youtu.be/vuzX0KbhgGA and https://youtu.be/q4XeK60NGOE
Submitted: 2019-01-31
Three-dimensional hydrodynamic numerical simulations have demonstrated that the structure of a protoplanetary disc may be strongly affected by a planet orbiting in a plane that is misaligned to the disc. When the planet is able to open a gap, the disc is separated into an inner, precessing disc and an outer disc with a warp. In this work, we compute infrared scattered light images to investigate the observational consequences of such an arrangement. We find that an inner disc misaligned by a less than a degree to the outer disc is indeed able to cast a shadow at larger radii. In our simulations a planet of around 6 Jupiter masses inclined by around 2 degrees is enough to warp the disc and cast a shadow with a depth of more than 10% of the average flux at that radius. We also demonstrate that warp in the outer disc can cause a variation in the azimuthal brightness profile at large radii. Importantly, this latter effect is a function of the distance from the star and is most prominent in the outer disc. We apply our model to the TW Hya system, where a misaligned, precessing inner disc has been invoked to explain an recently observed shadow in the outer disc. Consideration of the observational constraints suggest that an inner disc precessing due to a misaligned planet is an unlikely explanation for the features found in TW Hya.
[9]  oai:arXiv.org:1806.10148  [pdf] - 1790793
Gas and multi-species dust dynamics in viscous protoplanetary discs: the importance of the dust back-reaction
Comments: Accepted for publication in MNRAS, 20 pages, 9 figures
Submitted: 2018-06-26, last modified: 2018-11-21
We study the dynamics of a viscous protoplanetary disc hosting a population of dust grains with a range of sizes. We compute steady-state solutions, and show that the radial motion of both the gas and the dust can deviate substantially from those for a single-size dust population. Although the aerodynamic drag from the dust on the gas is weaker than in the case where all grains are optimally coupled to the gas, the cumulative "back-reaction" of the dust particles can still alter the gas dynamics significantly. In typical protoplanetary discs, the net effect of the dust back-reaction decreases the gas accretion flow compared to the dust-free (viscous) case, even for dust-to-gas ratios of order $1\%$. In the outer disc, where dust grains are typically less strongly coupled to the gas and settle towards the midplane, the dust back-reaction can even drive outward gas flow. Moreover, the radial inward drift of large grains is reduced below the gas motion in the inner disc regions, while small dust grains follow the gas dynamics over all the disc extent. The resulting dust and gas dynamics can give rise to observable structures, such as gas and dust cavities. Our results show that the dust back-reaction can play a major role in both the dynamics and observational appearance of protoplanetary discs, and cannot be ignored in models of protoplanetary disc evolution.
[10]  oai:arXiv.org:1605.04089  [pdf] - 1765956
A self-similar solution for thermal disc winds
Comments: Typographical errors in eqns 7,8,17 and 18 corrected. Revised first para. of 3.1; results unchanged
Submitted: 2016-05-13, last modified: 2018-10-12
We derive a self-similar description for the 2D streamline topology and flow structure of an axi-symmetric, thermally driven wind originating from a disc in which the density is a power law function of radius. Our scale-free solution is strictly only valid in the absence of gravity or centrifugal support; comparison with 2D hydrodynamic simulations of winds from Keplerian discs however demonstrates that the scale-free solution is a good approximation also in the outer regions of such discs, and can provide a reasonable description even for launch radii well within the gravitational radius of the flow. Although other authors have considered the flow properties along streamlines whose geometry has been specified in advance, this is the first isothermal calculation in which the flow geometry and variation of flow variables along streamlines is determined self-consistently. It is found that the flow trajectory is very sensitive to the power-law index of radial density variation in the disc: the steeper the density gradient, the stronger is the curvature of streamlines close to the flow base that is required in order to maintain momentum balance perpendicular to the flow. Steeper disc density profiles are also associated with more rapid acceleration, and a faster fall-off of density, with height above the disc plane. The derivation of a set of simple governing equations for the flow structure of thermal winds from the outer regions of power law discs offers the possibility of deriving flow observables without having to resort to hydrodynamical simulation.
[11]  oai:arXiv.org:1805.10449  [pdf] - 1753727
NGTS-2b: An inflated hot-Jupiter transiting a bright F-dwarf
Comments: 12 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2018-05-26, last modified: 2018-09-20
We report the discovery of NGTS-2b, an inflated hot-Jupiter transiting a bright F5V star (2MASS J14202949-3112074; $T_{\rm eff}$=$6478^{+94}_{-89}$ K), discovered as part of the Next Generation Transit Survey (NGTS). The planet is in a P=4.51 day orbit with mass $0.74^{+0.13}_{-0.12}$ M$_{J}$, radius $1.595^{+0.047}_{-0.045}$ R$_{J}$ and density $0.226^{+0.040}_{-0.038}$ g cm$^{-3}$; therefore one of the lowest density exoplanets currently known. With a relatively deep 1.0% transit around a bright V=10.96 host star, NGTS-2b is a prime target for probing giant planet composition via atmospheric transmission spectroscopy. The rapid rotation ($vsin$i=$15.2\pm0.8$ km s$^{-1}$) also makes this system an excellent candidate for Rossiter-McLaughlin follow-up observations, to measure the sky-projected stellar obliquity. NGTS-2b was confirmed without the need for follow-up photometry, due to the high precision of the NGTS photometry.
[12]  oai:arXiv.org:1808.05693  [pdf] - 1740132
Warping a protoplanetary disc with a planet on an inclined orbit
Comments: Accepted for publication in MNRAS. 17 pages, 12 Figures
Submitted: 2018-08-16
Recent observations of several protoplanetary discs have found evidence of departures from flat, circular motion in the inner regions of the disc. One possible explanation for these observations is a disc warp, which could be induced by a planet on a misaligned orbit. We present three-dimensional numerical simulations of the tidal interaction between a protoplanetary disc and a misaligned planet. For low planet masses we show that our simulations accurately model the evolution of inclined planet orbit (up to moderate inclinations). For a planet massive enough to carve a gap, the disc is separated into two components and the gas interior and exterior to the planet orbit evolve separately, forming an inner and outer disc. Due to the inclination of the planet, a warp develops across the planet orbit such that there is a relative tilt and twist between these discs. We show that when other parameters are held constant, the relative inclination that develops between the inner and outer disc depends on the outer radius of the total disc modelled. For a given disc mass, our results suggest that the observational relevance of the warp depends more strongly on the mass of the planet rather than the inclination of the orbit.
[13]  oai:arXiv.org:1802.09451  [pdf] - 1648939
Is the spiral morphology of the Elias 2-27 circumstellar disc due to gravitational instability?
Comments: 12 pages, 5 figures
Submitted: 2018-02-26
A recent ALMA observation of the Elias 2-27 system revealed a two-armed structure extending out to ~300 au in radius. The protostellar disc surrounding the central star is unusually massive, raising the possibility that the system is gravitationally unstable. Recent work has shown that the observed morphology of the system can be explained by disc self-gravity, so we examine the physical properties of the disc necessary to detect self-gravitating spiral waves. Using three-dimensional Smoothed Particle Hydrodynamics, coupled with radiative transfer and synthetic ALMA imaging, we find that observable spiral structure can only be explained by self-gravity if the disc has a low opacity (and therefore efficient cooling), and is minimally supported by external irradiation. This corresponds to a very narrow region of parameter space, suggesting that, although it is possible for the spiral structure to be due to disc self-gravity, other explanations, such as an external perturbation, may be preferred.
[14]  oai:arXiv.org:1801.05812  [pdf] - 1630364
Rings and gaps in the disc around Elias 24 revealed by ALMA
Comments: 17 pages, 11 figures. Accepted for publication in MNRAS
Submitted: 2018-01-17
We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.
[15]  oai:arXiv.org:1710.11099  [pdf] - 1648671
NGTS-1b: A hot Jupiter transiting an M-dwarf
Comments: 10 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2017-10-30
We present the discovery of NGTS-1b, a hot-Jupiter transiting an early M-dwarf host ($T_{eff}=3916^{+71}_{-63}~K$) in a P=2.674d orbit discovered as part of the Next Generation Transit Survey (NGTS). The planet has a mass of $0.812^{+0.066}_{-0.075}~M_{J}$, making it the most massive planet ever discovered transiting an M-dwarf. The radius of the planet is $1.33^{+0.61}_{-0.33}~R_{J}$. Since the transit is grazing, we determine this radius by modelling the data and placing a prior on the density from the population of known gas giant planets. NGTS-1b is the third transiting giant planet found around an M-dwarf, reinforcing the notion that close-in gas giants can form and migrate similar to the known population of hot Jupiters around solar type stars. The host star shows no signs of activity, and the kinematics hint at the star being from the thick disk population. With a deep (2.5%) transit around a $K=11.9$ host, NGTS-1b will be a strong candidate to probe giant planet composition around M-dwarfs via JWST transmission spectroscopy.
[16]  oai:arXiv.org:1710.06435  [pdf] - 1608519
Breaking mean-motion resonances during Type I planet migration
Comments: MNRAS; Accepted 16/10/2017, 12 pages, 9 figures, 1 table
Submitted: 2017-10-17
We present two-dimensional hydrodynamical simulations of pairs of planets migrating simultaneously in the Type I regime in a protoplanetary disc. Convergent migration naturally leads to the trapping of these planets in mean-motion resonances. Once in resonance the planets' eccentricity grows rapidly, and disc-planet torques cause the planets to escape resonance on a time-scale of a few hundred orbits. The effect is more pronounced in highly viscous discs, but operates efficiently even in inviscid discs. We attribute this resonance-breaking to overstable librations driven by moderate eccentricity damping, but find that this mechanism operates differently in hydrodynamic simulations than in previous analytic calculations. Planets escaping resonance in this manner can potentially explain the observed paucity of resonances in Kepler multi-transiting systems, and we suggest that simultaneous disc-driven migration remains the most plausible means of assembling tightly-packed planetary systems.
[17]  oai:arXiv.org:1708.03360  [pdf] - 1586983
Disk evolution and the fate of water
Comments: 26 pages, 8 figures: in the proceedings of the ISSI workshop on "The delivery of water to proto-planets, planets and satellites", to appear in Space Science Reviews
Submitted: 2017-08-10
We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolution, including dust growth, settling, and radial drift, and the level of turbulence and viscous dissipation within disks. Interferometric observations are now providing constraints on the positions of CO snow lines, but extrapolation to the unresolved regions where water ice sublimates will require much better theoretical understanding of mass and angular momentum transport in disks as well as more refined comparison of observations with sophisticated disk models.
[18]  oai:arXiv.org:1512.02649  [pdf] - 1342984
There might be giants: unseen Jupiter-mass planets as sculptors of tightly-packed planetary systems
Comments: 7 pages, 3 figures, 3 tables, accepted for publication in MNRAS
Submitted: 2015-12-08
The limited completeness of the Kepler sample for planets with orbital periods $\gtrsim$ 1 yr leaves open the possibility that exoplanetary systems may host undetected giant planets. Should such planets exist, their dynamical interactions with the inner planets may prove vital in sculpting the final orbital configurations of these systems. Using an $N$-body code with additional forces to emulate the effects of a protoplanetary disc, we perform simulations of the assembly of compact systems of super-Earth-mass planets with unseen giant companions. The simulated systems are analogous to Kepler-11 or Kepler-32 in that they contain 4 or 5 inner super-Earths, but our systems also contain longer-period giant companions which are unlikely to have been detected by Kepler. We find that giant companions tend to break widely-spaced, first-order mean-motion resonances, allowing the inner planets to migrate into tighter resonances. This leads to more compact architectures and increases the occurrence rate of Laplace resonant chains.
[19]  oai:arXiv.org:1512.01555  [pdf] - 1347754
Magnetospheres of hot Jupiters: hydrodynamic models & ultraviolet absorption
Comments: 13 pages, 12 figures. Accepted for publication in MNRAS
Submitted: 2015-12-04
We present hydrodynamic simulations of stellar wind-magnetosphere interactions in hot Jupiters such as WASP-12b. For fiducial stellar wind rates we find that a planetary magnetic field of a few G produces a large magnetospheric cavity, which is typically 6-9 planetary radii in size. A bow shock invariably forms ahead of the magnetosphere, but the pre-shock gas is only mildly supersonic (with typical Mach numbers of $\simeq$1.6-1.8) so the shock is weak. This results in a characteristic signature in the ultraviolet light curve: a broad absorption feature that leads the optical transit by 10-20% in orbital phase. The shapes of our synthetic light-curves are consistent with existing observations of WASP-12b, but the required near-UV optical depth ($\tau \sim 0.1$) can only be achieved if the shocked gas cools rapidly. We further show that radiative cooling is inefficient, so we deem it unlikely that a magnetospheric bow shock is responsible for the observed near-UV absorption. Finally, we apply our model to two other well-studied hot Jupiters (WASP-18b and HD209458b), and suggest that UV observations of more massive short-period planets (such as WASP-18b) will provide a straightforward test to distinguish between different models of circumplanetary absorption.
[20]  oai:arXiv.org:1509.02172  [pdf] - 1292338
The Evolution of Inner Disk Gas in Transition Disks
Comments: 18 pages, 14 figures, Accepted for publication in ApJ
Submitted: 2015-09-07
Investigating the molecular gas in the inner regions of protoplanetary disks provides insight into how the molecular disk environment changes during the transition from primordial to debris disk systems. We conduct a small survey of molecular hydrogen (H$_2$) fluorescent emission, using 14 well-studied Classical T Tauri stars at two distinct dust disk evolutionary stages, to explore how the structure of the inner molecular disk changes as the optically thick warm dust dissipates. We simulate the observed HI-Lyman $\alpha$-pumped H$_2$ disk fluorescence by creating a 2D radiative transfer model that describes the radial distributions of H$_{2}$ emission in the disk atmosphere and compare these to observations from the Hubble Space Telescope. We find the radial distributions that best describe the observed H$_2$ FUV emission arising in primordial disk targets (full dust disk) are demonstrably different than those of transition disks (little-to-no warm dust observed). For each best-fit model, we estimate inner and outer disk emission boundaries (r$_{in}$ and r$_{out}$), describing where the bulk of the observed H$_2$ emission arises in each disk, and we examine correlations between these and several observational disk evolution indicators, such as n$_{13-31}$, r$_{in,CO}$, and the mass accretion rate. We find strong, positive correlations between the H$_2$ radial distributions and the slope of the dust SED, implying the behavior of the molecular disk atmosphere changes as the inner dust clears in evolving protoplanetary disks. Overall, we find that H$_2$ inner radii are $\sim$4 times larger in transition systems, while the bulk of the H$_2$ emission originates inside the dust gap radius for all transitional sources.
[21]  oai:arXiv.org:1504.05562  [pdf] - 1300169
Probing the Final Stages of Protoplanetary Disk Evolution with ALMA
Comments: submitted to A&A
Submitted: 2015-04-21, last modified: 2015-08-17
The evolution of a circumstellar disk from its gas-rich protoplanetary stage to its gas-poor debris stage is not understood well. It is apparent that disk clearing progresses from the inside-out on a short time scale and models of photoevaporation are frequently used to explain this. However, the photoevaporation rates predicted by models differ by up to two orders of magnitude, resulting in uncertain time scales for the final stages of disk clearing. Photoevaporation theories predict that the final stages of disk-clearing progress in objects that have ceased accretion but still posses considerable material at radii far from the star. Weak-line T Tauri stars (WTTS) with infrared excess are likely in this configuration. We aim to provide observational constraints on theories of disk-clearing by measuring the dust masses and CO content of a sample of young (1.8-26.3 Myr) WTTS. We used ALMA Band 6 to obtain continuum and 12 CO(2-1) line fluxes for a sample of 24 WTTS stars with known infrared excess. We detect continuum emission in only four of 24 WTTS, and no 12 CO(2-1) emission in any. For those WTTS where no continuum was detected, their ages and derived upper limits suggest they are debris disks, which makes them some of the youngest debris disks known. Of those where continuum was detected, three are possible photoevaporating disks, although the lack of CO detection suggests a severely reduced gas-to-dust ratio. The low fraction of continuum detections implies that, once accretion onto the star stops, the clearing of the majority of dust progresses very rapidly. Most WTTS with infrared excess are likely not in transition but are instead young debris disks, whose dust is either primordial and has survived disk-clearing, or is of second-generation origin.
[22]  oai:arXiv.org:1502.07489  [pdf] - 1231882
Hubble Space Telescope observations of the NUV transit of WASP-12b
Comments: Accepted into the Astrophysical Journal
Submitted: 2015-02-26
We present new observations of four closely-spaced NUV transits of the hot Jupiter-like exoplanet WASP-12b using HST/COS, significantly increasing the phase resolution of the observed NUV light curve relative to previous observations, while minimising the temporal variation of the system. We observe significant excess NUV absorption during the transit, with mean normalised in-transit fluxes of $F_\mathrm{norm}\simeq0.97$, i.e. $\simeq$2-5 $\sigma$ deeper than the optical transit level of $\simeq0.986$ for a uniform stellar disk (the exact confidence level depending on the normalisation method used). We further observe an asymmetric transit shape, such that the post-conjunction fluxes are overall $\simeq$2-3 $\sigma$ higher than pre-conjunction values, and characterised by rapid variations in count rate between the pre-conjunction and out of transit levels. We do not find evidence for an early ingress to the NUV transit as suggested by earlier HST observations. However, we show that the NUV count rate observed prior to the optical transit is highly variable, but overall $\simeq$2.2-3.0 $\sigma$ below the post-transit values and comparable in depth to the optical transit, possibly forming a variable region of NUV absorption from at least phase $\phi\simeq$0.83, limited by the data coverage.
[23]  oai:arXiv.org:1409.3842  [pdf] - 1216921
Misaligned accretion on to supermassive black hole binaries
Comments: 13 pages, 10 figures, 3 tables. Accepted for publication in MNRAS. Simulation movies can be found at: http://www.acdunhill.com/binary-smbhs
Submitted: 2014-09-12
We present the results of high-resolution numerical simulations of gas clouds falling onto binary supermassive black holes to form circumbinary accretion discs, with both prograde and retrograde cloud orbits. We explore a range of clouds masses and cooling rates. We find that for low mass discs that cool fast enough to fragment, prograde discs are significantly shorter-lived than similar discs orbiting retrograde with respect to the binary. For fragmenting discs of all masses, we also find that prograde discs fragment across a narrower radial region. If the cooling is slow enough that the disc enters a self-regulating gravitoturbulent regime, we find that alignment between the disc and binary planes occurs on a timescale primarily dictated by the disc thickness. We estimate realistic cooling times for such discs, and find that in the majority of cases we expect fragmentation to occur. The longer lifetime of low-mass fragmenting retrograde discs allows them to drive significant binary evolution, and may provide a mechanism for solving the 'last parsec problem'.
[24]  oai:arXiv.org:1409.0532  [pdf] - 1216616
Understanding the assembly of Kepler's compact planetary systems
Comments: 12 pages, 4 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2014-09-01
The Kepler mission has recently discovered a number of exoplanetary systems, such as Kepler-11 and Kepler-32, in which ensembles of several planets are found in very closely packed orbits (often within a few percent of an AU of one another). These compact configurations present a challenge for traditional planet formation and migration scenarios. We present a dynamical study of the assembly of these systems, using an N-body method which incorporates a parametrized model of planet migration in a turbulent protoplanetary disc. We explore a wide parameter space, and find that under suitable conditions it is possible to form compact, close-packed planetary systems via traditional disc-driven migration. We find that simultaneous migration of multiple planets is a viable mechanism for the assembly of tightly-packed planetary systems, as long as the disc provides significant eccentricity damping and the level of turbulence in the disc is modest. We discuss the implications of our preferred parameters for the protoplanetary discs in which these systems formed, and comment on the occurrence and significance of mean-motion resonances in our simulations.
[25]  oai:arXiv.org:1311.1819  [pdf] - 1180576
The Dispersal of Protoplanetary Disks
Comments: 23 pages, 6 figures. Refereed review chapter, accepted for publication in Protostars & Planets VI, University of Arizona Press (2014), eds. H.Beuther, C.Dullemond, Th.Henning, R.Klessen
Submitted: 2013-11-07
Protoplanetary disks are the sites of planet formation, and the evolution and eventual dispersal of these disks strongly influences the formation of planetary systems. Disk evolution during the planet-forming epoch is driven by accretion and mass-loss due to winds, and in typical environments photoevaporation by high-energy radiation from the central star is likely to dominate final gas disk dispersal. We present a critical review of current theoretical models, and discuss the observations that are used to test these models and inform our understanding of the underlying physics. We also discuss the role disk dispersal plays in shaping planetary systems, considering its influence on both the process(es) of planet formation and the architectures of planetary systems. We conclude by presenting a schematic picture of protoplanetary disk evolution and dispersal, and discussing prospects for future work.
[26]  oai:arXiv.org:1308.1791  [pdf] - 1173278
Planet formation in evolving protoplanetary discs
Comments: 10 pages, 4 figures. Invited review to appear in Proceedings of IAU Symposium 299, "Exploring the formation and evolution of planetary systems"
Submitted: 2013-08-08
I attempt to summarize our knowledge of planet formation in evolving protoplanetary discs. I first review the physics of disc evolution and dispersal. For most of the disc lifetime evolution is driven by accretion and photoevaporation, and I discuss how the interplay between these processes shapes protoplanetary discs. I also discuss the observations that we use to test these models, and the major uncertainties that remain. I will then move on to consider planet formation and migration in evolving discs, and discuss how observations of both discs and planets can be used to inform our understanding of protoplanetary disc evolution.
[27]  oai:arXiv.org:1308.0596  [pdf] - 1173176
The curiously circular orbit of Kepler-16b
Comments: 7 pages, 6 figures. Accepted for publication in MNRAS. Movies of the simulations can be found at: http://www.astro.le.ac.uk/~acd23/K16.html
Submitted: 2013-08-02
The recent discovery of a number of circumbinary planets lends a new tool to astrophysicists seeking to understand how and where planet formation takes place. Of the increasingly numerous circumbinary systems, Kepler-16 is arguably the most dynamically interesting: it consists of a planet on an almost perfectly circular orbit (e = 0.0069) around a moderately eccentric binary (e = 0.16). We present high-resolution 3D smoothed-particle hydrodynamics simulations of a Kepler-16 analogue embedded in a circumbinary disc, and show that the planet's eccentricity is damped by its interaction with the protoplanetary disc. We use this to place a lower limit on the gas surface density in the real disc through which Kepler-16b migrated of \Sigma_min ~ 10 g cm^-2. This suggests that Kepler-16b, and other circumbinary planets, formed and migrated in relatively massive discs. We argue that secular evolution of circumbinary discs requires that these planets likely formed early on in the lifetime of the disc and migrated inwards before the disc lost a significant amount of its original mass.
[28]  oai:arXiv.org:1306.5241  [pdf] - 1172211
Formation of circumbinary planets in a dead zone
Comments: Accepted for publication in ApJ
Submitted: 2013-06-21
Circumbinary planets have been observed at orbital radii where binary perturbations may have significant effects on the gas disk structure, on planetesimal velocity dispersion, and on the coupling between turbulence and planetesimals. Here, we note that the impact of all of these effects on planet formation is qualitatively altered if the circumbinary disk structure is layered, with a non-turbulent midplane layer (dead zone) and strongly turbulent surface layers. For close binaries, we find that the dead zone typically extends from a radius close to the inner disk edge up to a radius of around 10-20 au from the centre of mass of the binary. The peak in the surface density occurs within the dead zone, far from the inner disk edge, close to the snow line, and may act as a trap for aerodynamically coupled solids. We suggest that circumbinary planet formation may be easier near this preferential location than for disks around single stars. However, dead zones around wide binaries are less likely and hence planet formation may be more difficult there.
[29]  oai:arXiv.org:1304.3746  [pdf] - 1165943
Hot Gas Lines in T Tauri Stars
Comments: Accepted in the ApJSS. 52 pages, 9 tables, 20 figures
Submitted: 2013-04-12
For Classical T Tauri Stars (CTTSs), the resonance lines of N V, Si IV, and C IV, as well as the He II 1640 A line, act as diagnostics of the accretion process. Here we assemble a large high-resolution dataset of these lines in CTTSs and Weak T Tauri Stars (WTTSs). We present data for 35 stars: one Herbig Ae star, 28 CTTSs, and 6 WTTSs. We decompose the C IV and He II lines into broad and narrow Gaussian components (BC & NC). The most common (50 %) C IV line morphology in CTTSs is that of a low-velocity NC together with a redshifted BC. The velocity centroids of the BCs and NCs are such that V_BC > 4 * V_NC, consistent with the predictions of the accretion shock model, in at most 12 out of 22 CTTSs. We do not find evidence of the post-shock becoming buried in the stellar photosphere due to the pressure of the accretion flow. The He II CTTSs lines are generally symmetric and narrow, less redshifted than the CTTSs C IV lines, by ~10 km/sec. The flux in the BC of the He II line is small compared to that of the C IV line, consistent with models of the pre-shock column emission. The observations are consistent with the presence of multiple accretion columns with different densities or with accretion models that predict a slow-moving, low-density region in the periphery of the accretion column. For HN Tau A and RW Aur A, most of the C IV line is blueshifted suggesting that the C IV emission is produced by shocks within outflow jets. In our sample, the Herbig Ae star DX Cha is the only object for which we find a P-Cygni profile in the C IV line, which argues for the presence of a hot (10^5 K) wind. For the overall sample, the Si IV and N V line luminosities are correlated with the C IV line luminosities, although the relationship between Si IV and C IV shows large scatter about a linear relationship and suggests that TW Hya, V4046 Sgr, AA Tau, DF Tau, GM Aur, and V1190 Sco are silicon-poor.
[30]  oai:arXiv.org:1303.0769  [pdf] - 1164985
Accretion Rates for T Tauri Stars Using Nearly Simultaneous Ultraviolet and Optical Spectra
Comments: 47 pages, 16 figures, accepted for publication in ApJ
Submitted: 2013-03-04
We analyze the accretion properties of 21 low mass T Tauri stars using a dataset of contemporaneous near ultraviolet (NUV) through optical observations obtained with the Hubble Space Telescope Imaging Spectrograph (STIS) and the ground based Small and Medium Aperture Research Telescope System (SMARTS), a unique dataset because of the nearly simultaneous broad wavelength coverage. Our dataset includes accreting T Tauri stars (CTTS) in Taurus, Chamaeleon I, $\eta$ Chamaeleon and the TW Hydra Association. For each source we calculate the accretion rate by fitting the NUV and optical excesses above the photosphere, produced in the accretion shock, introducing multiple accretion components characterized by a range in energy flux (or density) for the first time. This treatment is motivated by models of the magnetospheric geometry and accretion footprints, which predict that high density, low filling factor accretion spots co-exist with low density, high filling factor spots. By fitting the UV and optical spectra with multiple accretion components, we can explain excesses which have been observed in the near infrared. Comparing our estimates of the accretion rate to previous estimates, we find some discrepancies; however, they may be accounted for when considering assumptions for the amount of extinction and variability in optical spectra. Therefore, we confirm many previous estimates of the accretion rate. Finally, we measure emission line luminosities from the same spectra used for the accretion rate estimates, to produce correlations between accretion indicators (H$\beta$, Ca II K, C II] and Mg II) and accretion properties obtained simultaneously.
[31]  oai:arXiv.org:1210.6035  [pdf] - 1152349
A limit on eccentricity growth from global 3-D simulations of disc-planet interactions
Comments: 12 pages, 9 figures, 2 tables. Accepted for publication in MNRAS. Movies of some of the simulations can be found at: http://www.astro.le.ac.uk/~acd23/eccentricity.html
Submitted: 2012-10-22
We present high resolution 3-D simulations of the planet-disc interaction using smoothed particle hydrodynamics, to investigate the possibility of driving eccentricity growth by this mechanism. For models with a given disc viscosity (\alpha = 0.01), we find that for small planet masses (a few Jupiter masses) and canonical surface densities, no significant eccentricity growth is seen over the duration of our simulations. This contrasts with the limiting case of large planet mass (over twenty Jupiter masses) and extremely high surface densities, where we find eccentricity growth in agreement with previously published results. We identify the cause of this as being a threshold surface density for a given planet mass below which eccentricity growth cannot be excited by this method. Further, the radial profile of the disc surface density is found to have a stronger effect on eccentricity growth than previously acknowledged, implying that care must be taken when contrasting results from different disc models. We discuss the implication of this result for real planets embedded in gaseous discs, and suggest that the disc-planet interaction does not contribute significantly to observed exoplanet eccentricities.
[32]  oai:arXiv.org:1210.0446  [pdf] - 1151762
WD0837+185:the formation and evolution of an extreme mass ratio white dwarf-brown dwarf binary in Praesepe
Comments: Accepted for publication in ApJL
Submitted: 2012-10-01, last modified: 2012-10-02
There is a striking and unexplained dearth of brown dwarf companions in close orbits (< 3AU) around stars more massive than the Sun, in stark contrast to the frequency of stellar and planetary companions. Although rare and relatively short-lived, these systems leave detectable evolutionary end points in the form of white dwarf - brown dwarf binaries and these remnants can offer unique insights into the births and deaths of their parent systems. We present the discovery of a close (orbital separation ~ 0.006 AU) substellar companion to a massive white dwarf member of the Praesepe star cluster. Using the cluster age and the mass of the white dwarf we constrain the mass of the white dwarf progenitor star to lie in the range 3.5 - 3.7 Msun (B9). The high mass of the white dwarf means the substellar companion must have been engulfed by the B star's envelope while it was on the late asymptotic giant branch (AGB). Hence, the initial separation of the system was ~2 AU, with common envelope evolution reducing the separation to its current value. The initial and final orbital separations allow us to constrain the combination of the common envelope efficiency (alpha) and binding energy parameters (lambda) for the AGB star to alpha lambda ~3. We examine the various formation scenarios and conclude that the substellar object was most likely to have been captured by the white dwarf progenitor early in the life of the cluster, rather than forming in situ.
[33]  oai:arXiv.org:1209.0779  [pdf] - 634182
The dispersal of protoplanetary disks around binary stars
Comments: 5 pages, 2 figures. Accepted for publication in ApJ Letters
Submitted: 2012-09-04
I present models of disk evolution around young binary stars. I show that the primary factor in determining circumbinary disk lifetimes is the rate of disk photoevaporation. I also find that photoevaporative clearing leaves a signature on the distribution of circumbinary disk lifetimes, with a sharp increase in disk lifetimes for binary separations a<~0.3-1AU. Observations of young binary stars can therefore be used to test models of disk evolution, and I show that current data set a strong upper limit to the rate of on-going photoevaporation (<1e-9 M_sun/yr). Finally I discuss the implications of these results for planet formation, and suggest that circumbinary planets around close (a<~1AU) binaries should be relatively common.
[34]  oai:arXiv.org:1207.4789  [pdf] - 1124968
A Hubble Space Telescope Survey of H2 Emission in the Circumstellar Environments of Young Stars
Comments: ApJ - accepted. 19 pages, 12 figures
Submitted: 2012-07-19
The formation timescale and final architecture of exoplanetary systems are closely related to the properties of the molecular disks from which they form. Observations of the spatial distribution and lifetime of the molecular gas at planet-forming radii (r < 10 AU) are important for understanding the formation and evolution of exoplanetary systems. Towards this end, we present the largest spectrally resolved survey of H2 emission around low-mass pre-main sequence stars compiled to date. We use a combination of new and archival far-ultraviolet spectra from the COS and STIS instruments on the Hubble Space Telescope to sample 34 T Tauri stars (27 actively accreting CTTSs and 7 non-accreting WTTSs) with ages ranging from roughly 1-10 Myr. We observe fluorescent H2 emission, excited by LyA photons, in 100 of the accreting sources, including all of the transitional disks in our sample (CS Cha, DM Tau, GM Aur, UX Tau A, LkCa15, HD 135344B and TW Hya). The spatial distribution of the emitting gas is inferred from spectrally resolved H2 line profiles. Some of the emitting gas is produced in outflowing material, but the majority of H2 emission appears to originate in a rotating disk. For the disk-dominated targets, the H2 emission originates predominately at r < 3 AU. The emission line-widths and inner molecular radii are found to be roughly consistent with those measured from mid-IR CO spectra.
[35]  oai:arXiv.org:1203.6271  [pdf] - 493692
Circumstellar disks and planets. Science cases for next-generation optical/infrared long-baseline interferometers
Comments: 83 pages; Accepted for publication in "Astronomy and Astrophysics Review"; The final publication is available at http://www.springerlink.com
Submitted: 2012-03-28
We present a review of the interplay between the evolution of circumstellar disks and the formation of planets, both from the perspective of theoretical models and dedicated observations. Based on this, we identify and discuss fundamental questions concerning the formation and evolution of circumstellar disks and planets which can be addressed in the near future with optical and infrared long-baseline interferometers. Furthermore, the importance of complementary observations with long-baseline (sub)millimeter interferometers and high-sensitivity infrared observatories is outlined.
[36]  oai:arXiv.org:1202.5554  [pdf] - 634166
Deserts and pile-ups in the distribution of exoplanets due to photoevaporative disc clearing
Comments: 5 pages, 3 figures. Accepted for publication in MNRAS Letters
Submitted: 2012-02-24
We present models of giant planet migration in evolving protoplanetary discs. We show that disc clearing by EUV photoevaporation can have a strong effect on the distribution of giant planet semi-major axes. During disc clearing planet migration is slowed or accelerated in the region where photoevaporation opens a gap in the disc, resulting in "deserts" where few giant planets are found and corresponding "pile-ups" at smaller and larger radii. However, the precise locations and sizes of these features are strong functions of the efficiency of planetary accretion, and therefore also strongly dependent on planet mass. We suggest that photoevaporative disc clearing may be responsible for the pile-up of ~Jupiter-mass planets at ~1AU seen in exoplanet surveys, and show that observations of the distribution of exoplanet semi-major axes can be used to test models of both planet migration and disc clearing.
[37]  oai:arXiv.org:1110.6312  [pdf] - 1091259
NUV Excess in Slowly Accreting T Tauri Stars: Limits Imposed by Chromospheric Emission
Comments: 30 pages, 9 figures, Accepted to ApJ
Submitted: 2011-10-28
Young stars surrounded by disks with very low mass accretion rates are likely in the final stages of inner disk evolution and therefore particularly interesting to study. We present ultraviolet (UV) observations of the ~5-9 Myr old stars RECX-1 and RECX-11, obtained with the Cosmic Origins Spectrograph (COS) and Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST), as well as optical and near infrared spectroscopic observations. The two stars have similar levels of near UV emission, although spectroscopic evidence indicates that RECX-11 is accreting and RECX-1 is not. The line profiles of Halpha and He I 10830 in RECX-11 show both broad and narrow redshifted absorption components that vary with time, revealing the complexity of the accretion flows. We show that accretion indicators commonly used to measure mass accretion rates, e.g. U band excess luminosity or the Ca II triplet line luminosity are unreliable for low accretors, at least in the middle K spectral range. Using RECX-1 as a template for the intrinsic level of photospheric and chromospheric emission, we determine an upper limit of 3 x 10^-10 solar masses/ year for RECX-11. At this low accretion rate, recent photoevaporation models predict that an inner hole should have developed in the disk. However, the spectral energy distribution of RECX-11 shows fluxes comparable to the median of Taurus in the near infrared, indicating that substantial dust remains. Fluorescent H_2 emission lines formed in the innermost disk are observed in RECX-11, showing that gas is present in the inner disk, along with the dust.
[38]  oai:arXiv.org:1109.4148  [pdf] - 634160
Galactic Centre star formation: the case of the missing gas disc
Comments: 7 pages, 3 figures. Accepted for publication in MNRAS
Submitted: 2011-09-19
We study the dynamical evolution of stars and gas close to the centre of the Milky Way. Any plausible means of forming the young stars observed at the Galactic Centre leaves behind a residual gas disc at ~0.01pc radii. We show that the combined effects of viscous accretion and gravitational interactions with stars do not remove the residual gas efficiently, and that a substantial gas disc, interior to the stellar disc, persists for >10Myr after the stars form. Since no such disc is currently seen at the Galactic Centre we argue that it has been accreted by the super-massive black hole. This scenario offers an attractive connection between nuclear star formation and black hole feeding, and we suggest that the "missing" gas may have been used to power Sgr A*.
[39]  oai:arXiv.org:1109.0276  [pdf] - 1083764
The Relationship between Accretion Disc Age and Stellar Age and its Consequences for Proto-Stellar Discs
Comments: 11 pages, 10 figures, 2 table
Submitted: 2011-09-01
We show that for young stars which are still accreting and for which measurements of stellar age, disc mass and accretion rate are available, nominal disc age (Disc Age = Disc Mass / Accretion Rate) is approximately equal to the stellar age, at least within the considerable observational scatter. We then consider theoretical models of proto-stellar discs through analytic and numerical models. A variety of viscosity prescriptions including empirical power laws, magnetohydrodynamic turbulence and gravitational instability were considered within models describing the disc phenomena of dead zones, photoevaporation and planet formation. These models are generally poor fits to the observational data, showing values of 'Disc Age' which are too high by factors of 3 - 10. We then ask whether a systematic error in the measurement of one of the observational quantities might provide a reasonable explanation for this discrepancy. We show that for the observed systems only disc mass shows a systematic dependence on the value of 'Disc Age / Stellar Age' and we note that a systematic underestimate of the value of disc mass by a factor of around 3 - 5, would account for the discrepancy between theory and observations.
[40]  oai:arXiv.org:1108.3837  [pdf] - 1083479
Disc instability in RS Ophiuchi: a path to Type Ia supernovae?
Comments: 8 pages, 5 figures. Accepted for publication in MNRAS
Submitted: 2011-08-18
We study the stability of disc accretion in the recurrent nova RS Ophiuchi. We construct a one-dimensional time-dependent model of the binary-disc system, which includes viscous heating and radiative cooling and a self-consistent treatment of the binary potential. We find that the extended accretion disc in this system is always unstable to the thermal-viscous instability, and undergoes repeated disc outbursts on ~10-20yr time-scales. This is similar to the recurrence time-scale of observed outbursts in the RS Oph system, but we show that the disc's accretion luminosity during outburst is insufficient to explain the observed outbursts. We explore a range of models, and find that in most cases the accretion rate during outbursts reaches or exceeds the critical accretion rate for stable nuclear burning on the white dwarf surface. Consequently we suggest that a surface nuclear burning triggered by disc instability may be responsible for the observed outbursts. This allows the white dwarf mass to grow over time, and we suggest that disc instability in RS Oph and similar systems may represent a path to Type Ia supernovae.
[41]  oai:arXiv.org:1105.0045  [pdf] - 1076322
The Photoevaporative Wind from the Disk of TW Hya
Comments: Accepted for publication in ApJ
Submitted: 2011-04-30
Photoevaporation driven by the central star is expected to be a ubiquitous and important mechanism to disperse the circumstellar dust and gas from which planets form. Here, we present a detailed study of the circumstellar disk surrounding the nearby star TW Hya and provide observational constraints to its photoevaporative wind. Our new high-resolution (R ~ 30,000) mid-infrared spectroscopy in the [Ne II] 12.81 {\mu}m line confirms that this gas diagnostic traces the unbound wind component within 10AU from the star. From the blueshift and asymmetry in the line profile, we estimate that most (>80%) of the [Ne II] emission arises from disk radii where the midplane is optically thick to the redshifted outflowing gas, meaning beyond the 1 or 4AU dust rim inferred from other observations. We re-analyze high-resolution (R ~ 48, 000) archival optical spectra searching for additional transitions that may trace the photoevaporative flow. Unlike the [Ne II] line, optical forbidden lines from OI, SII, and MgI are centered at the stellar velocity and have symmetric profiles. The only way these lines could trace the photoevaporative flow is if they arise from a disk region physically distinct from that traced by the [Ne II] line, specifically from within the optically thin dust gap. However, the small (~10 km/s) FWHM of these lines suggest that most of the emitting gas traced at optical wavelengths is bound to the system rather than unbound. We discuss the implications of our results for a planet-induced versus a photoevaporation-induced gap.
[42]  oai:arXiv.org:1012.0724  [pdf] - 1042356
Global Gravitational Instabilities in Discs with Infall
Comments: 11 pages, 11 figures, Accepted by MNRAS
Submitted: 2010-12-03
Gravitational instability plays an important role in driving gas accretion in massive protostellar discs. Particularly strong is the global gravitational instability, which arises when the disc mass is of order 0.1 of the mass of the central star and has a characteristic spatial scale much greater than the disc's vertical scale-height. In this paper we use three-dimensional numerical hydrodynamics to study the development of gravitational instabilities in a disc which is embedded in a dense, gaseous envelope. We find that global gravitational instabilities are the dominant mode of angular momentum transport in the disc with infall, in contrast to otherwise identical isolated discs. The accretion torques created by low-order, global modes of the gravitational instability in a disc subject to infall are larger by a factor of several than an isolated disc of the same mass. We show that this global gravitational instability is driven by the strong vertical shear at the interface between the disc and the envelope, and suggest that this process may be an important means of driving accretion on to young stars.
[43]  oai:arXiv.org:0909.4309  [pdf] - 1017673
Radiation-Hydrodynamic Models of X-Ray & EUV Photoevaporating Protoplanetary Discs
Comments: 16 pages, 14 figures, accepted for publication in MNRAS
Submitted: 2009-09-24
(Abridged) We present the first radiation-hydrodynamic model of a protoplanetary disc irradiated with an X-EUV spectrum. In a model where the total ionizing luminosity is divided equally between X-ray and EUV luminosity, we find a photoevaporation rate of 1.4e-8 M_sun/yr, which is two orders of magnitude greater than the case of EUV photoevaporation alone. Thus it is clear that the X-rays are the dominant driving mechanism for photoevaporation. This can be understood inasmuch as X-rays are capable of penetrating much larger columns (~1e22 cm^-2) and can thus effect heating in denser regions and at larger radius than the EUV can. The radial extent of the launching region of the X-ray heated wind is 1-70AU compared with the pure EUV case where the launch region is concentrated around a few AU. When we couple our wind mass-loss rates with models for the disc's viscous evolution, we find that, as in the pure EUV case, there is a photoevaporative switch, such that an inner hole develops at ~ 1 AU at the point that the accretion rate in the disc drops below the wind mass loss rate. At this point, the remaining disc material is quickly removed in the final 15-20% of the disc's lifetime. This is consistent with the 1e5 yr transitional timescale estimated from observations of T-Tauri stars. We also caution that although our mass-loss rates are high compared to some accretion rates observed in young stars, our model has a rather large X-ray luminosity of 2e30 erg/s; further modeling is required in order to investigate the evolutionary implications of the large observed spread of X-ray luminosities in T-Tauri stars.
[44]  oai:arXiv.org:0909.0004  [pdf] - 901554
Giant planet migration, disk evolution, and the origin of transitional disks
Comments: 13 pages, 9 figures. Accepted for publication in ApJ
Submitted: 2009-08-31
We present models of giant planet migration in evolving protoplanetary disks. Our disks evolve subject to viscous transport of angular momentum and photoevaporation, while planets undergo Type II migration. We use a Monte Carlo approach, running large numbers of models with a range in initial conditions. We find that relatively simple models can reproduce both the observed radial distribution of extra-solar giant planets, and the lifetimes and accretion histories of protoplanetary disks. The use of state-of-the-art photoevaporation models results in a degree of coupling between planet formation and disk clearing, which has not been found previously. Some accretion across planetary orbits is necessary if planets are to survive at radii <~1.5AU, and if planets of Jupiter mass or greater are to survive in our models they must be able to form at late times, when the disk surface density in the formation region is low. Our model forms two different types of "transitional" disks, embedded planets and clearing disks, which show markedly different properties. We find that the observable properties of these systems are broadly consistent with current observations, and highlight useful observational diagnostics. We predict that young transition disks are more likely to contain embedded giant planets, while older transition disks are more likely to be undergoing disk clearing.
[45]  oai:arXiv.org:0809.0311  [pdf] - 15865
Massive black hole binary mergers within sub-pc scale gas discs
Comments: 11 pages, 13 figures. MNRAS, accepted. v2: minor changes
Submitted: 2008-09-01, last modified: 2008-11-03
[ABRIDGED] We study supermassive black hole binary mergers driven by angular momentum loss to small-scale gas discs. Such binaries form after major galaxy mergers, but their fate is unclear since hardening through stellar scattering becomes very inefficient at sub-parsec distances. Gas discs may dominate binary dynamics on these scales, and promote mergers. Using numerical simulations, we investigate the evolution of the orbits of binaries embedded within geometrically thin gas discs. Our simulations directly resolve angular momentum transport within the disc, which at the radii of interest is likely dominated by disc self-gravity. We show that the binary decays at a rate which is in good agreement with analytical estimates, while the eccentricity grows. Saturation of eccentricity growth is not observed up to values e > 0.35. Accretion onto the black holes is variable, and is roughly modulated by the binary orbital frequency. Scaling our results, we analytically estimate the maximum rate of binary decay that is possible without fragmentation occuring within the surrounding gas disc, and compare that rate to an estimate of the stellar dynamical hardening rate. For binary masses in the range 10^5 < M/Msun < 10^8 we find that decay due to gas discs may dominate for separations below a ~ 0.01--0.1 pc, where the disc is optically thick. The minimum merger time scale is shorter than the Hubble time for M < 10^7 Msun. This implies that gas discs could commonly attend relatively low mass black hole mergers, and that a significant population of binaries might exist at separations of a few hundredths of a pc, where the combined decay rate is slowest. For more massive binaries, we suggest that scattering of stars formed within a fragmenting gas disc could act as a significant additional sink of binary angular momentum.
[46]  oai:arXiv.org:0809.4137  [pdf] - 16674
Gas Evolution in Protoplanetary Disks
Comments: 9 pages, 6 figures
Submitted: 2008-09-24
This article summarizes a Splinter Session at the Cool Stars XV conference in St. Andrews with 3 review and 4 contributed talks. The speakers have discussed various approaches to understand the structure and evolution of the gas component in protoplanetary disks. These ranged from observational spectroscopy in the UV, infrared and millimeter, through to chemical and hydrodynamical models. The focus was on disks around low-mass stars, ranging from classical T Tauri stars to transitional disks and debris disks. Emphasis was put on water and organic molecules, the relation to planet formation, and the formation of holes and gaps in the inner regions.
[47]  oai:arXiv.org:0809.0316  [pdf] - 15869
[NeII] emission line profiles from photoevaporative disc winds
Comments: 5 pages, 3 figures. Accepted for publication in MNRAS Letters
Submitted: 2008-09-01
I model profiles of the [NeII] forbidden emission line at 12.81um, emitted by photoevaporative winds from discs around young, solar-mass stars. The predicted line luminosities (~ 1E-6 Lsun) are consistent with recent data, and the line profiles vary significantly with disc inclination. Edge-on discs show broad (30-40km/s) double-peaked profiles, due to the rotation of the disc, while in face-on discs the structure of the wind results in a narrower line (~10km/s) and a significant blue-shift (5-10km/s). These results suggest that observations of [NeII] line profiles can provide a direct test of models of protoplanetary disc photoevaporation.
[48]  oai:arXiv.org:0807.1731  [pdf] - 14410
Binary formation and mass function variations in fragmenting discs with short cooling times
Comments: 10 pages, 5 figures. Accepted for publication in MNRAS. Figures 1 and 3 degraded to meet arXiv size limits - version with high resolution figures available at http://www.strw.leidenuniv.nl/~rda/publications.html
Submitted: 2008-07-10
Accretion discs at sub-pc distances around supermassive black holes are likely to cool rapidly enough that self-gravity results in fragmentation. Here, we use high-resolution hydrodynamic simulations of a simplified disc model to study how the outcome of fragmentation depends upon numerical resolution and cooling time, and to investigate the incidence of binary formation within fragmenting discs. We investigate a range of cooling times, from the relatively long cooling time-scales that are marginally unstable to fragmentation down to highly unstable cooling on a time-scale that is shorter than the local dynamical time. The characteristic mass of fragments decreases with reduced cooling time, though the effect is modest and dependent upon details of how rapidly bound clumps radiate. We observe a high incidence of capture binaries, though we are unable to determine their final orbits or probability of survival. The results suggest that faster cooling in the parent disc results in an increased binary fraction, and that a high primordial binary fraction may result from disc fragmentation. We discuss our results in terms of the young massive stars close to the Galactic Centre, and suggest that observations of some stellar binaries close to the Galactic Centre remain consistent with formation in a fragmenting accretion disc.
[49]  oai:arXiv.org:0804.3596  [pdf] - 12003
Stellar dynamical evidence against a cold disc origin for stars in the Galactic Centre
Comments: 5 pages, 6 colour figures. MNRAS Letters in press. (v2: very minor changes)
Submitted: 2008-04-22, last modified: 2008-05-30
Observations of massive stars within the central parsec of the Galaxy show that, while most stars orbit within a well-defined disc, a significant fraction have large eccentricities and / or inclinations with respect to the disc plane. Here, we investigate whether this dynamically hot component could have arisen via scattering from an initially cold disc -- the expected initial condition if the stars formed from the fragmentation of an accretion disc. Using N-body methods, we evolve a variety of flat, cold, stellar systems, and study the effects of initial disc eccentricity, primordial binaries, very massive stars and intermediate mass black holes. We find, consistent with previous results, that a circular disc does not become eccentric enough unless there is a significant population of undetected 100--1000 Msun objects. However, since fragmentation of an eccentric disc can readily yield eccentric stellar orbits, the strongest constraints come from inclinations. We show that_none_ of our initial conditions yield the observed large inclinations, regardless of the initial disc eccentricity or the presence of massive objects. These results imply that the orbits of the young massive stars in the Galactic Centre are largely primordial, and that the stars are unlikely to have formed as a dynamically cold disc.
[50]  oai:arXiv.org:0712.0388  [pdf] - 7686
From discs to planetesimals I: evolution of gas and dust discs
Comments: 28 pages, 9 figures. Refereed review chapter for proceedings of VLTI summer school on "Circumstellar discs and planets at very high angular resolution", to appear in New Astronomy Reviews. See http://www.vlti.org/ for more details
Submitted: 2007-12-03
I review the processes that shape the evolution of protoplanetary discs around young, solar-mass stars. I first discuss observations of protoplanetary discs, and note in particular the constraints these observations place on models of disc evolution. The processes that affect the evolution of gas discs are then discussed, with the focus in particular on viscous accretion and photoevaporation, and recent models which combine the two. I then discuss the dynamics and growth of dust grains in discs, considering models of grain growth, the gas-grain interaction and planetesimal formation, and review recent research in this area. Lastly, I consider the so-called "transitional" discs, which are thought to be observed during disc dispersal. Recent observations and models of these systems are reviewed, and prospects for using statistical surveys to distinguish between the various proposed models are discussed.
[51]  oai:arXiv.org:0711.0759  [pdf] - 314899
Self-gravitating fragmentation of eccentric accretion disks
Comments: 9 pages, 7 figures. Accepted for publication in ApJ
Submitted: 2007-11-05
We consider the effects of eccentricity on the fragmentation of gravitationally unstable accretion disks, using numerical hydrodynamics. We find that eccentricity does not affect the overall stability of the disk against fragmentation, but significantly alters the manner in which such fragments accrete gas. Variable tidal forces around an eccentric orbit slow the accretion process, and suppress the formation of weakly-bound clumps. The "stellar" mass function resulting from the fragmentation of an eccentric disk is found to have a significantly higher characteristic mass than that from a corresponding circular disk. We discuss our results in terms of the disk(s) of massive stars at ~0.1pc from the Galactic Center, and find that the fragmentation of an eccentric accretion disk, due to gravitational instability, is a viable mechanism for the formation of these systems.
[52]  oai:arXiv.org:0706.2325  [pdf] - 2263
The Monitor project: JW 380 -- a 0.26, 0.15 Msol pre main sequence eclipsing binary in the Orion Nebula Cluster
Comments: 11 pages, 9 figures, accepted for publication in MNRAS
Submitted: 2007-06-15
We report the discovery of a low-mass (0.26 +/- 0.02, 0.15 +/- 0.01 Msol) pre-main-sequence eclipsing binary with a 5.3 day orbital period. JW 380 was detected as part of a high-cadence time-resolved photometric survey (the Monitor project) using the 2.5m Isaac Newton Telescope and Wide Field Camera for a survey of a single field in the Orion Nebula Cluster (ONC) region in V and i bands. The star is assigned a 99 per cent membership probability from proper motion measurements, and radial velocity observations indicate a systemic velocity within 1 sigma of that of the ONC. Modelling of the combined light and radial velocity curves of the system gave stellar radii of 1.19 +0.04 -0.18 Rsol and 0.90 +0.17 -0.03 Rsol for the primary and secondary, with a significant third light contribution which is also visible as a third peak in the cross-correlation functions used to derive radial velocities. The masses and radii appear to be consistent with stellar models for 2-3 Myr age from several authors, within the present observational errors. These observations probe an important region of mass-radius parameter space, where there are currently only a handful of known pre-main-sequence eclipsing binary systems with precise measurements available in the literature.
[53]  oai:arXiv.org:astro-ph/0611821  [pdf] - 87207
Dust dynamics during protoplanetary disc clearing
Comments: 13 pages, 5 figures. Accepted for publication in MNRAS
Submitted: 2006-11-27
We consider the dynamics of dust and gas during the clearing of protoplanetary discs. We work within the context of a photoevaporation/viscous model for the evolution of the gas disc, and use a two-fluid model to study the dynamics of dust grains as the gas disc is cleared. Small (<~10um) grains remain well-coupled to the gas, but larger (~1mm) grains are subject to inward migration from large radii (~50AU), suggesting that the time-scale for grain growth in the outer disc is ~10^4-10^5yr. We describe in detail the observable appearance of discs during clearing, and find that pressure gradients in the gas disc result in a strong enhancement of the local dust-to-gas ratio in a ring near to the inner disc edge. Lastly, we consider a simple model of the disc-planet interaction, and suggest that observations of disc masses and accretion rates provide a straightforward means of discriminating between different models of disc clearing.
[54]  oai:arXiv.org:astro-ph/0609812  [pdf] - 316560
Constraints on the Stellar Mass Function from Stellar Dynamics at the Galactic Center
Comments: 9 pages, 5 figures. Accepted for publication in ApJ
Submitted: 2006-09-29
We consider the dynamical evolution of a disk of stars orbiting a central black hole. In particular, we focus on the effect of the stellar mass function on the evolution of the disk, using both analytic arguments and numerical simulations. We apply our model to the ring of massive stars at ~0.1pc from the Galactic Center, assuming that the stars formed in a cold, circular disk, and find that our model requires the presence of a significant population of massive (>100Msun) stars in order to explain the the observed eccentricities of 0.2-0.3. Moreover, in order to limit the damping of the heavier stars' eccentricities, we also require fewer low-mass stars than expected from a Salpeter mass function, giving strong evidence for a significantly ``top-heavy'' mass function in the rings of stars seen near to the Galactic Center. We also note that the maximum possible eccentricities attainable from circular initial conditions at ages of <10Myr are around 0.4-0.5, and suggest that any rings of stars found with higher eccentricities were probably not formed from circular disks.
[55]  oai:arXiv.org:astro-ph/0603253  [pdf] - 80462
Photoevaporation of protoplanetary discs I: hydrodynamic models
Comments: 13 pages, 9 figures. Accepted for publication in MNRAS
Submitted: 2006-03-09
In this paper we consider the effect of the direct ionizing stellar radiation field on the evolution of protoplanetary discs subject to photoevaporative winds. We suggest that models which combine viscous evolution with photoevaporation of the disc (e.g. Clarke, Gendrin & Sotomayor 2001) incorrectly neglect the direct field after the inner disc has drained, at late times in the evolution. We construct models of the photoevaporative wind produced by the direct field, first using simple analytic arguments and later using detailed numerical hydrodynamics. We find that the wind produced by the direct field at late times is much larger than has previously been assumed, and we show that the mass-loss rate scales as $R_{in}^{1/2}$ (where $R_{in}$ is the radius of the instantaneous inner disc edge). We suggest that this result has important consequences for theories of disc evolution, and go on to consider the effects of this result on disc evolution in detail in a companion paper (Alexander, Clarke & Pringle 2006b).
[56]  oai:arXiv.org:astro-ph/0603254  [pdf] - 80463
Photoevaporation of protoplanetary discs II: evolutionary models and observable properties
Comments: 11 pages, 5 figures. Accepted for publication in MNRAS
Submitted: 2006-03-09
We present a new model for protoplanetary disc evolution. This model combines viscous evolution with photoevaporation of the disc, in a manner similar to Clarke, Gendrin & Sotomayor (2001). However in a companion paper (Alexander, Clarke & Pringle 2006a) we have shown that at late times such models must consider the effect of stellar radiation directly incident on the inner disc edge, and here we model the observational implications of this process. We find that the entire disc is dispersed on a time-scale of order $10^5$yr after a disc lifetime of a few Myr, consistent with observations of T Tauri (TT) stars. We use a simple prescription to model the spectral energy distribution of the evolving disc, and demonstrate that the model is consistent with observational data across a wide range of wavelengths. We note also that the model predicts a short ``inner hole'' phase in the evolution of all TT discs, and make predictions for future observations at mid-infrared and millimetre wavelengths.
[57]  oai:arXiv.org:astro-ph/0602059  [pdf] - 79635
The stellar mass-accretion rate relation in T Tauri stars and brown dwarfs
Comments: 4 pages, 2 figures. Accepted for publication in ApJ Letters
Submitted: 2006-02-02
Recent observations show a strong correlation between stellar mass and accretion rate in young stellar and sub-stellar objects, with the scaling $\dot{M}_{acc} \propto M_*^2$ holding over more than four orders of magnitude in accretion rate. We explore the consequences of this correlation in the context of disk evolution models. We note that such a correlation is not expected to arise from variations in disk angular momentum transport efficiency with stellar mass, and suggest that it may reflect a systematic trend in disk initial conditions. In this case we find that brown dwarf disks initially have rather larger radii than those around more massive objects. By considering disk evolution, and invoking a simple parametrization for a shut-off in accretion at the end of the disk lifetime, we show that such models predict that the scatter in the stellar mass-accretion rate relationship should increase with increasing stellar mass, in rough agreement with current observations.
[58]  oai:arXiv.org:astro-ph/0501100  [pdf] - 70212
Constraints on the ionizing flux emitted by T Tauri stars
Comments: 8 pages, 5 figures. Accepted for publication in MNRAS
Submitted: 2005-01-06
We present the results of an analysis of ultraviolet observations of T Tauri Stars (TTS). By analysing emission measures taken from the literature we derive rates of ionizing photons from the chromospheres of 5 classical TTS in the range ~10^41-10^44 photons/s, although these values are subject to large uncertainties. We propose that the HeII/CIV line ratio can be used as a reddening-independent indicator of the hardness of the ultraviolet spectrum emitted by TTS. By studying this line ratio in a much larger sample of objects we find evidence for an ionizing flux which does not decrease, and may even increase, as TTS evolve. This implies that a significant fraction of the ionizing flux from TTS is not powered by the accretion of disc material onto the central object, and we discuss the significance of this result and its implications for models of disc evolution. The presence of a significant ionizing flux in the later stages of circumstellar disc evolution provides an important new constraint on disc photoevaporation models.
[59]  oai:arXiv.org:astro-ph/0406643  [pdf] - 65800
The effects of X-ray photoionization and heating on the structure of circumstellar discs
Comments: 11 pages, 7 figures. Accepted for publication in MNRAS
Submitted: 2004-06-29
We present the results of a theoretical study investigating the effects of photoionization and heating by X-rays on discs around low-mass stars. In particular we address the question of whether or not X-rays can drive a disc wind. First, we construct a 1-dimensional ``quasi-hydrostatic'' model, which solves for the vertical structure introduced by X-ray heating. We consider uniform X-ray illumination of the disc, but the X-ray fluxes required to heat the disc significantly are much greater than those seen by recent observations. When the model is extended to consider heating from a central X-ray source we find that the 1-dimensional model is only valid very close to the star. We extend our analysis to consider a simple 2-dimensional model, treating the disc as a two-layered structure and solving for its density profile self-consistently. For T Tauri stars we are able to set a crude upper limit on the mass-loss rate that can be driven by X-ray photoevaporation, with a value of ~10^-13 g/cm^2/s. Our model is designed to maximise this value, and most likely over-estimates it significantly. However we still find a mass-loss rate which is less than that found in studies of ultraviolet photoevaporation. We conclude that in the presence of a significant UV field, X-ray driven disc winds are unlikely to play a significant role in the evolution of discs around low-mass stars.
[60]  oai:arXiv.org:astro-ph/0311276  [pdf] - 60856
On the origin of ionising photons emitted by T Tauri stars
Comments: 6 pages, 4 figures. Accepted for publication in MNRAS
Submitted: 2003-11-12
We address the issue of the production of Lyman continuum photons by T Tauri stars, in an attempt to provide constraints on theoretical models of disc photoionisation. By treating the accretion shock as a hotspot on the stellar surface we show that Lyman continuum photons are produced at a rate approximately three orders of magnitude lower than that produced by a corresponding black body, and that a strong Lyman continuum is only emitted for high mass accretion rates. When our models are extended to include a column of material accreting on to the hotspot we find that the accretion column is extremely optically thick to Lyman continuum photons. Further, we find that radiative recombination of hydrogen atoms within the column is not an efficient means of producing photons with energies greater than 13.6eV, and find that an accretion column of any conceivable height suppresses the emission of Lyman continuum photons to a level below or comparable to that expected from the stellar photosphere. The photospheric Lyman continuum is itself much too weak to affect disc evolution significantly, and we find that the Lyman continuum emitted by an accretion shock is similarly unable to influence disc evolution significantly. This result has important consequences for models which use photoionisation as a mechanism to drive the dispersal of circumstellar discs, essentially proving that an additional source of Lyman continuum photons must exist if disc photoionisation is to be significant.
[61]  oai:arXiv.org:astro-ph/0302178  [pdf] - 54805
ISOCAM-CVF spectroscopy of the circumstellar environment of young stellar objects
Comments: 15 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics. Version including appendix containing atlas of spectra available at http://www.ast.cam.ac.uk/~rda/cvfpaper.ps.gz
Submitted: 2003-02-10
We present the results of a mid-infrared imaging spectroscopy survey of Young Stellar Objects (YSOs) and their surrounding environment in four low-mass star formation regions: RCrA, Rho Ophiuchi, Serpens and Chamaeleon I. This survey was performed using ISOCAM and its Circular Variable Filters (CVF) and observed 42 YSO candidates. A number of spectral features were measured, most notably the 9.7um silicate feature, the bending modes of both water and CO2 ices at 6.0 and 15.2um respectively and the well-known unidentified feature at 6.8um. The strength of the unidentified feature was observed to correlate very well with that of the water ice bending mode and far less strongly with the CO2 ice bending mode. This suggests, in a manner consistent with previous observations, that the carrier of the unidentified feature is a strongly polar ice. Absorption profiles of the bending mode of CO2 ice are observed to show a significant long wavelength wing, which suggests that a significant fraction of the CO2 ice observed exists in a polar (H2O-rich) phase. The sources observed in RCrA, Rho Oph and Serpens show similar spectral characteristics, whilst the sources observed in Cha I are somewhat anomalous. However this is consistent with previous studies of this region of the Cha I cloud, which contains an unusual cluster of YSOs. We detect an apparent under-abundance of water ice towards the sources in Rho Oph. This may be indicative of differences in chemical composition between the different clouds, or may be due to evaporation. Finally the CO2:H2O ice ratios observed towards the sources in Rho Oph show significantly greater scatter than in the other regions, possibly due to varying local conditions around the YSOs in Rho Oph. (Abridged)