Search

Gas has been detected in many exoplanetary systems ($>$10 Myr), thought to be released in the destruction of volatile-rich planetesimals orbiting in exo-Kuiper belts. In this letter, we aim to explore whether gas is also expected in the Kuiper belt (KB) in our Solar System. To quantify the gas release in our Solar System, we use models for gas release that have been applied to extrasolar planetary systems, as well as a physical model that accounts for gas released due to the progressive internal warming of large planetesimals. We find that only bodies larger than about 4 km can still contain CO ice after 4.6 Gyr of evolution. This finding may provide a clue as to why Jupiter-family comets, thought to originate in the Kuiper belt, are deficient in CO compared to Oort-clouds comets. We predict that gas is still produced in the KB right now at a rate of $2 \times 10^{-8}$ M$_\oplus$/Myr for CO and orders of magnitude more when the Sun was younger. Once released, the gas is quickly pushed out by the Solar wind. Therefore, we predict a gas wind in our Solar System starting at the KB location and extending far beyond with regards to the heliosphere with a current total CO mass of $\sim 2 \times 10^{-12}$ M$_\oplus$. We also predict the existence of a slightly more massive atomic gas wind made of carbon and oxygen (neutral and ionized) with a mass of $\sim 10^{-11}$ M$_\oplus$. We predict that gas is currently present in our Solar System beyond the Kuiper belt and that although it cannot be detected with current instrumentation, it could be observed in the future with an in situ mission using an instrument similar to Alice on New Horizons with larger detectors. Our model of gas release due to slow heating may also work for exoplanetary systems and provide the first real physical mechanism for the gas observations., accepted for publication as a Letter to the editor in A&A; abstract shortened; 15 pages, 5 figures, 4 tables

Be stars are rapidly rotating B stars with Balmer emission lines that indicate the presence of a Keplerian, rotationally supported, circumstellar gas disc. Current disc models, referred to as "decretion discs", make use of the zero torque inner boundary condition typically applied to accretion discs, with the 'decretion' modelled by adding mass to the disc at a radius of about two per cent larger than the inner disc boundary. We point out that, in this model, the rates at which mass and energy need to be added to the disc are implausibly large. What is required is that the disc has not only a source of mass but also a continuing source of angular momentum. We argue that the disc evolution may be more physically modelled by application of the non-zero torque inner boundary condition of Nixon & Pringle (2020), which determines the torque applied at the boundary as a fraction of the advected angular momentum flux there and approaches the accretion and decretion disc cases in the appropriate limits. We provide supporting arguments for the suggestion that the origin of the disc material is small-scale magnetic flaring events on the stellar surface, which, when combined with rapid rotation, can provide sufficient mass to form, and sufficient angular momentum to maintain, a Keplerian Be star disc. We discuss the origin of such small-scale magnetic fields in radiative stars with differential rotation. We conclude that small-scale magnetic fields on the stellar surface, may be able to provide the necessary mass flux and the necessary time-dependent torque on the disc inner regions to drive the observed disc evolution., 6 pages, no figures, accepted for publication in The Astrophysical Journal Letters

Magnetic fields are observed in star forming regions. However simulations of the late stages of star formation that do not include magnetic fields provide a good fit to the properties of young stars including the initial mass function (IMF) and the multiplicity. We argue here that the simulations that do include magnetic fields are unable to capture the correct physics, in particular the high value of the magnetic Prandtl number, and the low value of the magnetic diffusivity. The artificially high (numerical and uncontrolled) magnetic diffusivity leads to a large magnetic flux pervading the star forming region. We argue further that in reality the dynamics of high magnetic Prandtl number turbulence may lead to local regions of magnetic energy dissipation through reconnection, meaning that the regions of molecular clouds which are forming stars might be essentially free of magnetic fields. Thus the simulations that ignore magnetic fields on the scales on which the properties of stellar masses, stellar multiplicities and planet-forming discs are determined, may be closer to reality than those which include magnetic fields, but can only do so in an unrealistic parameter regime., Comment: 28 pages, 0 figures, accepted for publication in New Astronomy

The rapidly evolving dust and gas extinction observed towards WD 1145+017 has opened a real-time window onto the mechanisms for destruction-accretion of planetary bodies onto white dwarf stars, and has served to underline the importance of considering the dynamics of dust particles around such objects. Here it is argued that the interaction between (charged) dust grains and the stellar magnetic field is an important ingredient in understanding the physical distribution of infrared emitting particles in the vicinity of such white dwarfs. These ideas are used to suggest a possible model for WD 1145+017 in which the unusual transit shapes are caused by opaque clouds of dust trapped in the stellar magnetosphere. The model can account for the observed transit periodicities if the stellar rotation is near 4.5 h, as the clouds of trapped dust are then located near or within the co-rotation radius. The model requires the surface magnetic field to be at least around some tens of kG. In contrast to the eccentric orbits expected for large planetesimals undergoing tidal disintegration, the orbits of magnetospherically-trapped dust clouds are essentially circular, consistent with the observations., Comment: 5 pages, accepted to MNRAS Letters

Motivated by recent observations of short-timescale variations in the infrared emission of circumstellar disks, we propose that coronal mass ejections can remove dust grains on timescales as short as a few days. Continuous monitoring of stellar activity, coupled to infrared observations, can place meaningful constraints on the proposed mechanism., Comment: accepted to Astrophysical Journal Letters

Standard, planar accretion discs operate through a dissipative mechanism, usually thought to be turbulent, and often modelled as a viscosity. This acts to take energy from the radial shear, enabling the flow of mass and angular momentum in the radial direction. In a previous paper we discussed observational evidence for the magnitude of this viscosity, and pointed out discrepancies between these values and those obtained in numerical simulations. In this paper we discuss the observational evidence for the magnitude of the dissipative effects which act in non--planar discs, both to transfer and to eliminate the non--planarity. Estimates based on the model by Ogilvie (1999), which assumes a small--scale, isotropic viscosity, give alignment timescales for fully ionized discs which are apparently too short by a factor of a few compared with observations, although we emphasize that more detailed computations as well as tighter observational constraints are required to verify this conclusion. For discs with low temperature and conductivity, we find that the timescales for disc alignment based on isotropic viscosity are too short by around two orders of magnitude. This large discrepancy suggests that our understanding of viscosity in quiescent discs is currently inadequate., Comment: 7 pages, no figures. MNRAS, in press

We propose two new means of identifying the main class of progenitors of Type Ia supernovae--single or double degenerate: (i) If the range of supernova properties is significantly determined by the range of viewing angles of non-spherically symmetric explosions, then the nature of the correlation between polarization and another property (for example, the velocity gradient) can be used to determine the geometry of the asymmetry and hence the nature of the progenitor, and (ii) in the double- but not in the single-degenerate case, the range in the observed properties (e.g., velocity gradients) is likely to increase with the amount of carbon seen in the ejecta., Comment: 10 pages, 2 figures; ApJL in press

Gas at intermediate temperature between the hot X-ray emitting coronal gas in galaxies at the centers of galaxy clusters, and the much cooler optical line emitting filaments, yields information on transport processes and plausible scenarios for the relationship between X-ray cool cores and other galactic phenomena such as mergers or the onset of an active galactic nucleus. Hitherto, detection of intermediate temperature gas has proven elusive. Here, we present FUV imaging of the "low excitation" emission filaments of M87 and show strong evidence for the presence of CIV 1549 A emission which arises in gas at temperature ~10^5K co-located with Halpha+[NII] emission from cooler ~10^4K gas. We infer that the hot and cool phases are in thermal communication, and show that quantitatively the emission strength is consistent with thermal conduction, which in turn may account for many of the observed characteristics of cool core galaxy clusters., Comment: Accepted for publication in ApJ Letters

Sokoloski et al (2008) have recently reported evidence that the recurrent nova RS Ophiuchi produced a pair of highly collimated radio jets within days of its 2006 outburst. This suggests that an accretion disc must be present during the outburst. However in the standard picture of recurrent novae as thermonuclear events, any such disc must be expelled from the white dwarf vicinity, as the nuclear energy yield greatly exceeds its binding energy. We suggest instead that the outbursts of RS Oph are thermal--viscous instabilities in a disc irradiated by the central accreting white dwarf. The distinctive feature of RS Oph is the very large size of its accretion disc. Given this, it fits naturally into a consistent picture of systems with unstable accretion discs. This picture explains the presence and speed of the jets, the brightness and duration of the outburst, and its rise time and linear decay, as well as the faintness of the quiescence. By contrast, the hitherto standard picture of recurrent thermonuclear explosions has a number of severe difficulties. These include the presence of jets, the faintness of quiescence, and the fact the the accretion disc must be unstable unless it is far smaller than any reasonable estimate., Comment: MNRAS, in press

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are thought to be caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion, $\sigma_k \approx 15\rm km s^{-1}$, is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically $i \ge 25^\circ$. Alternatively a higher velocity kick distribution, $\sigma_k = 265 \rm km s^{-1}$, is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularisation of the binary orbits has taken place. We discuss briefly how this might be achieved., Comment: Accepted for publication in MNRAS

The microquasar GRO J1655-40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady state disc warped by Lense-Thirring precession and find the timescale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically fifty times higher in the latter case but we find that, in both cases, the lifetime of the mass transfer state is at most a few times the alignment timescale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter-aligned after its supernova kick but that it is most likely to be close to alignment rather than counteralignment now., Comment: Accepted for publication in MNRAS

We suggest that most nearby active galactic nuclei are fed by a series of small--scale, randomly--oriented accretion events. Outside a certain radius these events promote rapid star formation, while within it they fuel the supermassive black hole. We show that the events have a characteristic time evolution. This picture agrees with several observational facts. The expected luminosity function is broadly in agreement with that observed for moderate--mass black holes. The spin of the black hole is low, and aligns with the inner disc in each individual feeding event. This implies radio jets aligned with the axis of the obscuring torus, and uncorrelated with the large--scale structure of the host galaxy. The ring of young stars observed about the Galactic Centre are close to where our picture predicts that star formation should occur., Comment: MNRAS, in press

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)., 13 pages, 9 figures. Accepted for publication in MNRAS

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., Comment: 8 pages, 5 figures. Accepted for publication in MNRAS

In this paper we consider the evolution of small planetesimals in marginally stable, self-gravitating protoplanetary discs. The drag force between the disc gas and the embedded planetesimals generally causes the planetesimals to drift inwards through the disc at a rate that depends on the particle size. In a marginally stable, self-gravitating disc, however, the planetesimals are significantly influenced by the non-axisymmetric spiral structures resulting from the growth of the gravitational instability. For small particles, that are strongly coupled to the disc gas, and for large particles, that have essentially decoupled from the disc gas, the effect is not particularly significant. Intermediate sized particles, which would generally have the largest radial drift rates, do, however, become significantly concentrated at the peaks of the spiral arms. These high density regions may persist for, of order, an orbital period and may attain densities comparable to that of the disc gas. Although at the end of the simulation only ~ 25 % of the planetesimal particles lie in regions of enhanced density, during the course of the simulation at least 75 % of the planetesimal particles have at some stage been in a such a region. We find that the concentration of particles in the spiral arms results in an increased collision rate, an effect that could significantly accelerate planetesimal growth. The density enhancements may also be sufficient for the growth of planetesimals through direct gravitational collapse. The interaction between small planetesimals and self-gravitating spiral structures may therefore play an important role in the formation of large planetesimals that will ultimately coagulate to form terrestrial planets or the cores of gas/ice giant planets., Comment: 11 pages, 9 figures, Accepted for publication in MNRAS

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic ($\gamma = 5/3$). Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance ${\cal M}_{\rm r} \ga 0.01$ undergo shocks within a distance of order the resonance radius., Comment: 28 pages, 21 figures - 8 as GIF, 13 embedded postscript, Accepted for publication in MNRAS. Full postscript version available from http://www.astro.ex.ac.uk/people/mbate

We investigate the migration of massive extrasolar planets due to gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 AU at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log (a), where a is the orbital radius. Between 1 AU and 3 AU, the expected number of planets per logarithmic interval in orbital radius roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, about 15 percent of stars targeted by existing radial velocity searches are predicted to possess planets with masses 0.3 M_Jupiter < M_p sin (i) < 10 M_Jupiter, and radii 0.1 AU < a < 5 AU. A third of these planets (around 5 percent of the target stars) lie at the radii most amenable to detection via microlensing. A further 5-10 percent of stars could have planets at radii of 5 AU < a < 8 AU that have migrated outwards. We discuss the probability of forming a system (akin to the Solar System) in which significant radial migration of the most massive planet does not occur. About 10-15 percent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low mass planets than high mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function., MNRAS, in press

We demonstrate that, by combining optical, radio and X-ray observations of a Seyfert, it is possible to provide a direct measurement of the angle $\beta$ between the direction of the radio jet and the normal to the plane of the spiral host galaxy. To do so, we make the assumptions that the inner radio jet is perpendicular to the X-ray observed inner accretion disk, and that the observed jet (or the stronger component, if the jet is two-sided) is physically closer to Earth than the plane of the galaxy. We draw attention to the possibility of measurement producing a result which is not self-consistent, in which case for that galaxy, one of the assumptions must fail., Comment: 11 pages, 1 figure, accepted for publication in The Astrophysical Journal Letters

We present a new method for the statistical investigation into the distributions of the angle beta between the radio axis and the normal to the galactic disk for a sample of Seyfert galaxies. We discuss how further observations of the sample galaxies can strengthen the conclusions. Our data are consistent with the hypothesis that AGN jets are oriented randomly in space, independent of the position of the plane of the galaxy. By making the simple assumption that the Standard Model of AGN holds, with a universal opening angle of the thick torus of phi_c, we demonstrate a statistical method to obtain an estimate of phi_c. Our data are not consistent with the simple-minded idea that Seyfert 1s and Seyfert 2s are differentiated solely by whether or not our line of sight lies within some fixed angle of the jet axis. Our result is significant on the 2 sigma level and can thus be considered only suggestive, not conclusive. A complete sample of Seyfert galaxies selected on an isotropic property is required to obtain a conclusive result., Comment: 13 pages, Tex, 5 Postscript figures. Accepted ApJ