Probing the circumstellar structures of T Tauri stars and their relationship to those of Herbig stars

We present Hα spectropolarimetry observations of a sample of 10 bright T Tauri stars, supplemented with new Herbig Ae/Be star data. A change in the linear polarization across Hα is detected in most of the T Tauri (9/10) and Herbig Ae (9/11) objects, which we interpret in terms of a compact source of line photons that is scattered off a rotating accretion disc. We find consistency between the position angle (PA) of the polarization and those of imaged disc PAs from infrared and millimetre imaging and interferometry studies, probing much larger scales. For the Herbig Ae stars AB Aur, MWC 480 and CQ Tau, we find the polarization PA to be perpendicular to the imaged disc, which is expected for single scattering. On the other hand, the polarization PA aligns with the outer disc PA for the T Tauri stars DR Tau and SU Aur and FU Ori, conforming to the case of multiple scattering. This difference can be explained if the inner discs of Herbig Ae stars are optically thin, whilst those around our T Tauri stars and FU Ori are optically thick. Furthermore, we develop a novel technique that combines known inclination angles and our recent Monte Carlo models to constrain the inner rim sizes of SU Aur, GW Ori, AB Aur and CQ Tau. Finally, we consider the connection of the inner disc structure with the orientation of the magnetic field in the foreground interstellar medium: for FU Ori and DR Tau, we infer an alignment of the stellar axis and the larger magnetic field direction.     

A search for planets transiting the M-dwarf debris disc host, AU Microscopii

We present high-cadence, high-precision multiband photometry of the young, M1Ve, debris disc star, AU Microscopii. The data were obtained in three continuum filters spanning a wavelength range from 4500 to 6600 Å, plus Hα, over 28 nights in 2005. The light curves show intrinsic stellar variability due to star-spots with an amplitude in the blue band of 0.051 mag and a period of 4.847 d. In addition, three large flares were detected in the data which all occur near the minimum brightness of the star. We remove the intrinsic stellar variability and combine the light curves of all the filters in order to search for transits by possible planetary companions orbiting in the plane of the nearly edge-on debris disc. The combined final light curve has a sampling of 0.35 min and a standard deviation of 6.8 mmag. We performed Monte Carlo simulations by adding fake transits to the observed light curve and find with 95 per cent significance that there are no Jupiter mass planets orbiting in the plane of the debris disc on circular orbits with periods,   P ≤ 5  d. In addition, there are no young Neptune like planets (with radii 2.5 times smaller than the young Jupiter) on circular orbits with periods,   P ≤ 3  d.     

Clues about the star formation history of the M31 disc from WFPC2 photometry

Over the past several years, the Hubble Space Telescope ( HST ) has acquired many broad-band images of various regions in the M31 disc. I have obtained 27 such fields from the HST data archive in order to produce colour–magnitude diagrams (CMDs) of the stellar populations contained within these areas of the disc. I have attempted to reproduce these CMDs using theoretical stellar evolution models in conjunction with statistical tools for determining the star formation history that best fits the observations. The wide range of extinction values within any given field makes the data difficult to reproduce accurately; nevertheless, I have managed to find star formation histories that roughly reproduce the data. These statistically determined star formation histories reveal that, like the disc of the Galaxy, the disc of M31 contains very few old metal-poor stars. The histories also suggest that the star formation rate of the disc as a whole has been low over the past ∼1 Gyr.     

Accretion disc–stellar magnetosphere interaction: field line inflation and the effect on the spin-down torque

We calculate the structure of a force-free magnetosphere which is assumed to corotate with a central star and which interacts with an embedded differentially rotating accretion disc. The magnetic and rotation axes are aligned, and the stellar field is assumed to be a dipole. We concentrate on the case when the amount of field line twisting through the disc–magnetosphere interaction is large , and consider different outer boundary conditions. In general the field line twisting produces field line inflation (e.g. Bardou & Heyvaerts), and in some cases with large twisting many field lines can become open. We calculate the spin-down torque acting between the star and the disc, and we find that it decreases significantly for cases with large field line twisting. This suggests that the oscillating torques observed for some accreting neutron stars could be caused by the magnetosphere varying between states with low and high field line inflation. Calculations of the spin evolution of T Tauri stars may also have to be revised in the light of the significant effect that field line twisting has on the magnetic torque resulting from star–disc interactions.     

Self-gravitating warped discs around supermassive black holes

We consider warped equilibrium configurations for stellar and gaseous discs in the Keplerian force field of a supermassive black hole, assuming that the self-gravity of the disc provides the only acting torques. Modelling the disc as a collection of concentric circular rings and computing the torques in the non-linear regime, we show that stable, strongly warped precessing equilibria are possible. These solutions exist for a wide range of disc-to-black-hole mass ratios   M _d/ M _bh  , can span large warp angles of up to  ±∼120°  , have inner and outer boundaries, and extend over a radial range of a factor of typically two to four. These equilibrium configurations obey a scaling relation such that in good approximation     where     is the (retrograde) precession frequency and Ω is a characteristic orbital frequency in the disc. Stability was determined using linear perturbation theory and, in a few cases, confirmed by numerical integration of the equations of motion. Most of the precessing equilibria are found to be stable, but some are unstable. The main result of this study is that highly warped discs near black holes can persist for long times without any persistent forcing other than by their self-gravity. The possible relevance of this to galactic nuclei is briefly discussed.     

Hα spectropolarimetry of the Herbig Ae star AB Aurigae

We present spectropolarimetric observations, obtained at H α , of the Herbig Ae star AB Aurigae. Changes in linear polarization across the H α line probe structure in the immediate circumstellar environment of the central star, down to scales of the order of one to a few stellar radii. In the case of AB Aurigae the observed polarimetric signature is complex. After applying a correction for foreground continuum polarization, we find that there is a linear-polarized H α emission component intrinsic to the source. Rotation of the angle of polarization through the emission-line profile suggests scattering in a rotating circumstellar disc. The magnetic accretor model commonly applied to T Tauri stars shows promise of explaining these data.     

Light scattering by moving dust grains in the immediate vicinity of young stars

We consider the problem of the distortion of the photospheric spectrum for a young star as its light is scattered in the inner accretion disk in the dust grain evaporation region. In T Tauri stars, this region is at a distance of the order of several stellar radii and is involved in the large-scale motions of matter with velocities of ～100 km s^?1 or higher. The light scattering in such a medium causes the frequency of the scattered radiation to be shifted due to the Doppler effect. We analyze the influence of this effect on the absorption line profiles in the spectra of T Tauri stars using classical results of the theory of radiative transfer. We consider two models of a scattering medium: (i) a homogeneous cylindrical surface and (ii) a cylindrical surface with an azimuth-dependent height (such conditions take place during the accretion of matter onto a star with an oblique magnetic dipole). We show that in the first case, the scattering of the photospheric radiation causes the absorption lines to broaden. If the motion of the circumstellar matter in the dust evaporation region is characterized by two velocity components, then the line profile of the scattered radiation is asymmetric, with the pattern of the asymmetry depending on the direction of the radial velocity. In the second case, the scattered radiation can cause periodic shifts of the absorption line centroid, which can be perceived by an observer as periodic radial-velocity variations in the star. We suggest that precisely this effect is responsible for the low-amplitude radial-velocity variations with periods close to the stellar rotation periods that have recently been found in some of the T Tauri stars.     

Surface inhomogeneities and line variability on DF Tau

We have mapped surface inhomogeneities on the classical T Tauri star DF Tau, using the Li  i doublet at 670.8 nm, the Ca  i lines at 612.2 nm and 643.9 nm and a calcium and iron blend at 646.3 nm. We find compelling evidence that there are hotspots with temperatures of more than 5000 K. Two of the hotspots produce line-profile deformations that can be traced as they move through the cross-correlated profiles. When one of the hotspots crosses the stellar disc, redshifted absorption components appear in the Na D lines. As these redshifted absorption features are usually tracers for mass-infall we interpret this hotspot as an accretion shock close to the stellar surface.   Parts of the surface of DF Tau are covered with a hot chromosphere that is visible in the Ca  ii infrared triplet lines and the narrow component of He  i . We find no correlation between the veiling and the lines that originate from the hot chromosphere, suggesting that the veiling and the chromospheric emission are produced in physically distinct regions.     

The star formation rate in disc galaxies: thresholds and dependence on gas amount

We reassess the applicability of the Toomre criterion in galactic discs and we study the local star formation law in 16 disc galaxies for which abundance gradients are published. The data we use consist of stellar light profiles, atomic and molecular gas (deduced from CO with a metallicity-dependent conversion factor), star formation rates (from Hα emissivities), metallicities, dispersion velocities and rotation curves. We show that the Toomre criterion applies successfully to the case of the Milky Way disc, but it has limited success with the data of our sample; depending on whether or not the stellar component is included in the stability analysis, we find average values for the threshold ratio of the gas surface density to the critical surface density in the range 0.5–0.7. We also test various star formation laws proposed in the literature, i.e. either the simple Schmidt law or modifications of it, that take into account dynamical factors. We find only small differences among them as far as the overall fit to our data is concerned; in particular, we find that all three star formation laws (with parameters derived from the fits to our data) match observations in the Milky Way disc particularly well. In all cases we find that the exponent n of our best-fitting star formation rate has slightly higher values than in other recent works and we suggest several reasons that may cause that discrepancy.     

On the effects of the stellar magnetic field on the structure of T Tauri accretion discs

The structure of accretion discs around magnetic T Tauri stars is calculated numerically using a particle hydrodynamical code, in which magnetic interaction is included in the framework of King's diamagnetic blob accretion model. Setting up the calculation so as to simulate the density structure of a quasi-steady disc in the equatorial plane of a T Tauri star, we find that the central star's magnetic field typically produces a central hole in the disc and spreads out the surface density distribution. We argue that this result suggets a promising mechanism for explaining the unusual flatness (IR excess) of T Tauri accretion disc spectra.     

Alignment and precession of a black hole with a warped accretion disc

We consider the shape of an accretion disc whose outer regions are misaligned with the spin axis of a central black hole and calculate the steady state form of the warped disc in the case where the viscosity and surface densities are power laws in the distance from the central black hole. We discuss the shape of the resulting disc in both the frame of the black hole and that of the outer disc. We note that some parts of the disc and also any companion star maybe shadowed from the central regions by the warp. We compute the torque on the black hole caused by the Lense–Thirring precession, and hence compute the alignment and precession time-scales. We generalize the case with viscosity and hence surface density independent of radius to more realistic density distributions for which the surface density is a decreasing function of radius. We find that the alignment time-scale does not change greatly but the precession time-scale is more sensitive. We also determine the effect on this time-scale if we truncate the disc. For a given truncation radius, the time-scales are less affected for more sharply falling density distributions.     

Surface magnetic fields on two accreting T Tauri stars: CV Cha and CR Cha

We have produced brightness and magnetic field maps of the surfaces of CV Cha and CR Cha: two actively accreting G- and K-type T Tauri stars in the Chamaeleon I star-forming cloud with ages of 3–5 Myr. Our magnetic field maps show evidence for strong, complex multipolar fields similar to those obtained for young rapidly rotating main-sequence stars. Brightness maps indicate the presence of dark polar caps and low-latitude spots – these brightness maps are very similar to those obtained for other pre-main-sequence and rapidly rotating main-sequence stars.
Only two other classical T Tauri stars have been studied using similar techniques so far: V2129 Oph and BP Tau. CV Cha and CR Cha show magnetic field patterns that are significantly more complex than those recovered for BP Tau, a fully convective T Tauri star.
We discuss possible reasons for this difference and suggest that the complexity of the stellar magnetic field is related to the convection zone; with more complex fields being found in T Tauri stars with radiative cores (V2129 Oph, CV Cha and CR Cha). However, it is clearly necessary to conduct magnetic field studies of T Tauri star systems, exploring a wide range of stellar parameters in order to establish how they affect magnetic field generation, and thus how these magnetic fields are likely to affect the evolution of T Tauri star systems as they approach the main sequence.     

Circumstellar structure of RU Lupi down to au scales

We have used the technique of spectro-astrometry to study the milliarcsecond scale structure of the emission lines in the T Tauri star RU Lupi. The wings of the H α emission are found to be displaced from the star towards the south-west (blue wing) and north-east (red wing) with angular scales of 20–30 mas. This structure is consistent with a bipolar outflow from the star. From a study of the variability of the intensity and position spectra, we argue that a combination of magnetically driven bipolar outflow and accreting gas contributes to the H α emission. On the other hand, the [O  i ] and [S  ii ] emission are displaced from the star to the south-west but at much larger distances than the H α , hundreds of milliarcseconds for the high-velocity component (HVC) and down to 30 mas for the low-velocity components (LVCs). The presence of both redshifted and blueshifted outflows in H α but only a blueshifted outflow in the forbidden lines can be explained if the disc obscures the redshifted forbidden line outflow, but a disc gap with outer radius 3–4 au allows the redshifted H α to be seen. This gap could be induced by an unseen companion.     

Black hole spin and radio loudness in a Λ cold dark matter universe

We use a combination of a cosmological N -body simulation of the concordance Λ cold dark matter paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the α model of Shakura & Sunyaev and consider the King et al. warped disc alignment criterion. The orientation of the accretion disc is inferred from the angular momentum of the source of accreted material, which bears a close relationship to the large-scale structure in the simulation. We find that the final BH spin depends almost exclusively on the accretion history and only weakly on the warped disc alignment. The main mechanisms of BH spin-up are found to be gas cooling processes and disc instabilities, a result that is only partially compatible with Monte Carlo models where the main spin-up mechanisms are major mergers and disc instabilities; the latter results are reproduced when implementing randomly oriented accretion discs in our model. Regarding the BH population, we find that more massive BHs, which are hosted by massive ellipticals, have higher spin values than less massive BHs, hosted by spiral galaxies. We analyse whether gas accretion rates and BH spins can be used as tracers of the radio loudness of active galactic nuclei (AGN). We find that the current observational indications of an increasing trend of radio-loud AGN fractions with stellar and BH mass can be easily obtained when placing lower limits on the BH spin, with a minimum influence from limits on the accretion rates; a model with random accretion disc orientations is unable to reproduce this trend. Our results favour a scenario where the BH spin is a key parameter to separate the radio-loud and radio-quiet galaxy populations.     

The rapidly oscillating Ap star HD 99563 and its distorted dipole pulsation mode

We undertook a time-series photometric multisite campaign for the rapidly oscillating Ap (roAp) star HD 99563 and also acquired mean light observations over four seasons. The pulsations of the star, which show flatter light maxima than minima, can be described with a frequency quintuplet centred on 1557.653 μHz and some first harmonics of it. The amplitude of the pulsation is modulated with the rotation period of the star that we determine with 2.91179 ± 0.00007 d from the analysis of the stellar pulsation spectrum and of the mean light data. We break up the distorted oscillation mode into its pure spherical harmonic components and find it is dominated by the ℓ= 1 pulsation, and also has a notable ℓ= 3 contribution, with weak ℓ= 0 and 2 components. The geometrical configuration of the star allows us to see both pulsation poles for about the same amount of time; HD 99563 is only the fourth roAp star for which both pulsation poles are seen and only the third where the distortion of the pulsation modes has been modelled. We point out that HD 99563 is very similar to the well-studied roAp star HR 3831. Finally, we note that the visual companion of HD 99563 is located in the δ Scuti instability strip and may thus show pulsation. We show that if the companion was physical, the roAp star would be a 2.03-M_⊙, object, seen at a rotational inclination of 44°, which then predicts a magnetic obliquity     .     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

The origin of the light distribution in spiral galaxies

We analyse a high-resolution, fully cosmological, hydrodynamical disc galaxy simulation, to study the source of the double-exponential light profiles seen in many stellar discs, and the effects of stellar radial migration upon the spatiotemporal evolution of both the disc age and metallicity distributions. We find a 'break' in the pure exponential stellar surface brightness profile, and trace its origin to a sharp decrease in the star formation per unit surface area, itself produced by a decrease in the gas volume density due to a warping of the gas disc. Star formation in the disc continues well beyond the break. We find that the break is more pronounced in bluer wavebands. By contrast, we find little or no break in the mass density profile. This is, in part, due to the net radial migration of stars towards the external parts of the disc. Beyond the break radius, we find that ∼60 per cent of the resident stars migrated from the inner disc, while ∼25 per cent formed in situ . Our simulated galaxy also has a minimum in the age profile at the break radius but, in disagreement with some previous studies, migration is not the main mechanism producing this shape. In our simulation, the disc metallicity gradient flattens with time, consistent with an 'inside-out' formation scenario. We do not find any difference in the intensity or the position of the break with inclination, suggesting that perhaps the differences found in empirical studies are driven by dust extinction.     

Hydrodynamic simulations of propagating warps and bending waves in accretion discs

We present the results of a study of propagating warp or bending waves in accretion discs. Three-dimensional hydrodynamic simulations were performed using smoothed particle hydrodynamics (SPH), and the results are compared with calculations based on the linear theory of warped discs.
We examine the response of a gaseous disc to an initially imposed warping disturbance under a variety of physical conditions. We consider primarily the physical regime in which the dimensionless viscosity parameter α < H r , where H r is the disc aspect ratio, so that bending waves are expected to propagate. We also performed calculations for disc models in which α > H r , where the warps are expected to evolve diffusively. Small-amplitude (linear) perturbations are studied in both Keplerian and slightly non-Keplerian discs, and we find that the results of the SPH calculations can be reasonably well fitted by those of the linear theory. The main results of these calculations are: (i) the warp in Keplerian discs when α < H r propagates with little dispersion, and damps at a rate expected from estimates of the code viscosity; (ii) warps evolve diffusively when α > H r ; (iii) the slightly non-Keplerian discs lead to a substantially more dispersive behaviour of the warps, which damp at a similar rate to the Keplerian case, when α < H r .
Initially imposed higher amplitude, non-linear warping disturbances were studied in Keplerian discs. The results indicate that non-linear warps can lead to the formation of shocks, and that the evolution of the warp becomes less wave-like and more diffusive in character.
This work is relevant to the study of the warped accretion discs that may occur around Kerr black holes or in misaligned binary systems, and is mainly concerned with discs in which α < H r . The results indicate that SPH can model the hydrodynamics of warped discs, even when using rather modest numbers of particles.     

The effect of radiation stress on gas motions in AGN accretion discs

In this paper we look at one of the effects of irradiation on a warped accretion disc in the context of active galactic nuclei (AGN). A warp will catch a substantial amount of the radiation emitted by the central object. We consider the fluid motions that may arise inside a warped disc when the surface is subject to a radiation stress, and also the net mass flows that result. We find that, to first order, we have a balance of the viscous and Coriolis-type forces. The radial radiation stress causes outward motion of the surface layer, but only the azimuthal Poynting–Robertson drag leads to an increase in the net accretion rate. We investigate the distribution of the velocity perturbations and find them to be significant in determining the local structure of the disc.
An unexpected result is that the picture changes significantly when we take into account the periodic illumination of the warped disc. A type of resonance at the local Keplerian rotation frequency causes a flow that penetrates the whole thickness of the disc; these flows are faster than the flows due to unchanging illumination. They totally dominate the induced flows in terms of sheer mass, but significant impact on disc structure still occurs only near the surface, where velocity perturbations typically go up to some kilometres per second.     

T Tauri star magnetic fields and magnetospheres

Stellar magnetic fields govern key aspects of the evolution of a young star, from controlling accretion to regulating the angular momentum evolution of the system. Spectro‐polarimetric studies of T Tauri stars have revealed a surprising range of magnetic field topologies. Meanwhile multi‐wavelength campaigns have probed T Tauri star systems from stellar photosphere to inner disk, allowing us to study magnetospheric accretion in unprecedented detail. We review recent results and discuss their implications for understanding the evolution of young stars (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)