Source author record

P. J. Armitage

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

17 published item(s)

preprint2022arXiv

Lifetime of the Outer Solar System Nebula From Carbonaceous Chondrites

The evolution and lifetime of protoplanetary disks (PPDs) play a central role in the formation and architecture of planetary systems. Astronomical observations suggest that PPDs evolve in two timescales, accreting onto the star for up to several million years (Myr) followed by gas dissipation within <1 Myr. Because solar nebula magnetic fields are sustained by the gas of the protoplanetary disk, we can use paleomagnetic measurements to infer the lifetime of the solar nebula. Here, we use paleomagnetic measurements of meteorites to constrain this lifetime and investigate whether the solar nebula had a two-timescale evolution. We report on paleomagnetic measurements of bulk subsamples of two CO carbonaceous chondrites: Allan Hills A77307 and Dominion Range 08006. If magnetite in these meteorites can acquire a crystallization remanent magnetization that recorded the ambient field during aqueous alteration, our measurements suggest that the local magnetic field strength at the CO parent body location was <0.9 \muT at some time between 2.7 and 5.1 Myr after the formation of calcium-aluminum-rich inclusions. Coupled with previous paleomagnetic studies, we conclude that the dissipation of the solar nebula in the 3-7 AU region occurred <1.5 Myr after the dissipation of the nebula in the 1-3 AU region, suggesting that protoplanetary disks go through a two-timescale evolution in their lifetime, consistent with dissipation by photoevaporation and/or magnetohydrodynamic winds. We also discuss future directions necessary to obtain robust records of solar nebula fields using bulk chondrites, including obtaining ages from meteorites and experimental work to determine how magnetite acquires magnetization during chondrite parent body alteration.

preprint2014arXiv

Convergence of simulations of self-gravitating accretion discs II: Sensitivity to the implementation of radiative cooling and artificial viscosity

Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.

preprint2014arXiv

Distinguishing an ejected blob from alternative flare models at the Galactic centre with GRAVITY

The black hole at the Galactic centre exhibits regularly flares of radiation, the origin of which is still not understood. In this article, we study the ability of the near-future GRAVITY infrared instrument to constrain the nature of these events. We develop realistic simulations of GRAVITY astrometric data sets for various flare models. We show that the instrument will be able to distinguish an ejected blob from alternative flare models, provided the blob inclination is >= 45deg, the flare brightest magnitude is 14 <= mK <= 15 and the flare duration is >= 1h30.

preprint2011arXiv

Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

preprint2011arXiv

Stability of self-gravitating discs under irradiation

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

preprint2009arXiv

Black hole mergers: the first light

The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.

preprint2009arXiv

Giant planet migration, disk evolution, and the origin of transitional disks

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.

preprint2009arXiv

The role of disc self-gravity in the formation of protostars and protostellar discs

We use time-dependent, one-dimensional disc models to investigate the evolution of protostellar discs that form through the collapse of molecular cloud cores and in which the primary transport mechanism is self-gravity. We assume that these discs settle into a state of thermal equilibrium with Q = 2 and that the strength of the angular momentum transport is set by the cooling rate of the disc. The results suggest that these discs will attain a quasi-steady state that persists for a number of free-fall times and in which most of the mass within 100 au is located inside 10-20 au. This pile-up of mass in the inner disc could result in temperatures that are high enough for the growth of MHD turbulence which could rapidly drain the inner disc and lead to FU Orionis-like outbursts. In all our simulations, the inner regions of the discs (r < 40 au) were stable against fragmentation, while fragmentation was possible in the outer regions (r > 40 au) of discs that formed from cores that had enough initial angular momentum to deposit sufficient mass in these outer regions. The large amounts of mass in these outer regions, however, suggests that fragmentation will lead to the formation of sub-stellar and stellar mass companions, rather than planetary mass objects. Although mass accretion rates were largely consistent with observations, the large disc masses suggest that an additional transport mechanism (such as MRI occuring in the upper layers of the disc) must operate in order to drain the remaining disc material within observed disc lifetimes.

preprint2009arXiv

Time-dependent models of the structure and evolution of self-gravitating protoplanetary discs

Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disk is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (alpha ~ 10^{-3} or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several au at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the \dot{M}-M_* relation for very young systems may allow a test of the model.

preprint2007arXiv

Self-gravitating fragmentation of eccentric accretion disks

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.

preprint2006arXiv

Constraints on the Stellar Mass Function from Stellar Dynamics at the Galactic Center

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.

preprint2003arXiv

Constraints on a planetary origin for the gap in GM Aurigae's protoplanetary disc

The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.

preprint1996arXiv

The ejection of T Tauri stars from molecular clouds and the fate of circumstellar discs

We investigate the evolution of circumstellar discs around T Tauri stars that are ejected from small stellar clusters within molecular clouds. In particular, we study how the interaction that leads to ejection may hasten the transition between Classical and Weak-lined T Tauri status. In our models, ejections of T Tauri stars at velocities of 3-10 km/s truncate the accretion disc at radii between 1 and 10 a.u., reducing the viscous evolution time of the disc so that accretion rapidly ceases. The observational appearance of the resulting systems is then dependent on the presence or absence of a stellar magnetic field. For non-magnetic stars we find that a near-infra red excess should persist due to reprocessing of stellar radiation, but that this is greatly diminished for magnetic T Tauri stars by the presence of a magnetosphere extending to corotation. In either case, there is a period when ejected stars should appear as non-accreting systems with detectable circumstellar material at wavelengths of 5 microns and beyond. We discuss the implications of these results for models in which ejected stars contribute to the halo of pre-main-sequence objects discovered from ROSAT observations of star forming regions and the All-Sky Survey.

preprint1995arXiv

Dynamo driven accretion discs and dwarf nova eruptions

We explore the consequences of a magnetic dynamo origin for the viscosity in accretion discs, for the structure and evolution of discs in dwarf nova systems. We propose that the rapid cooling that sets in at the end of a dwarf nova eruption acts to inhibit the Balbus-Hawley instability, and thereby to quench dynamo action and so reduce disc viscosity. We demonstrate that a modified disc instability model can reproduce the basic properties of dwarf nova eruptions, as well as some properties of quiescent discs. We also discuss some observational consequences of our model.

preprint1995arXiv

Magnetic braking of T Tauri stars

We construct models for the rotation rates of T Tauri stars whose spin is regulated by magnetic linkage between the star and a surrounding accretion disc. Our models utilise a time-dependent disc code to follow the accretion process and include the effects of pre-main-sequence stellar evolution. We find that the initial disc mass controls the evolution of the star-disc system. For sufficiently massive discs, a stellar field of $\sim$ 1 kG is able to regulate the spin rate to the observed values during the classical T Tauri phase. The field then acts to expel the disc and the star spins up at constant angular momentum as a weak-line system. Lower mass discs are ejected at an early epoch and fail to brake the star significantly. We extend the model to close binary systems, and find that the removal of angular momentum from the disc by the secondary significantly prolongs the inner disc lifetime. Such systems should therefore be relatively slow rotators. We also discuss the implications of our model for the spectral energy distributions and variability of T Tauri stars.

preprint1995arXiv

Magnetic cycles and photometric variability of T Tauri stars

We consider the effects of stellar magnetic cycles on disc accretion in T Tauri stars where the inner region of the disc is disrupted by the stellar field. Using time-dependent disc models, we show that order unity variations in field strength can strongly modulate the accretion flow for cycle time-scales of years to decades. Photometric variability is expected to be greatest in the ultraviolet, with smaller but still significant variations probable in the near-infrared. We suggest that variable stellar magnetic fields may contribute to the long-time-scale photometric variability of T Tauri stars. The long-term variability of T Tauri stars may therefore provide a clue to the nature of magnetic dynamo cycles in these stars.

preprint1995arXiv

Magnetically modulated accretion in T Tauri stars

We examine how accretion on to T Tauri stars may be modulated by a time-dependent `magnetic gate' where the inner edge of the accretion disc is disrupted by a varying stellar field. We show that magnetic field variations on time-scales shorter than 10^5 yr can modulate the accretion flow, thus providing a possible mechanism both for the marked photometric variability of T Tauri stars and for the possible conversion of T Tauri stars between classical and weak line status. We thus suggest that archival data relating to the spectrophotometric variability of T Tauri stars may provide an indirect record of magnetic activity cycles in low-mass pre-main-sequence stars.