Roche tomography of cataclysmic variables - III. Star-spots on AE Aqr

We present a Roche tomography reconstruction of the secondary star in the cataclysmic variable AE Aqr. The tomogram reveals several surface inhomogeneities that are due to the presence of large, cool star-spots. In addition to a number of lower latitude spots, the maps also show the presence of a large, high-latitude spot similar to that seen in Doppler images of rapidly rotating isolated stars, and a relative paucity of spots at a latitude of 40 degrees. In total, we estimate that some 18 per cent of the Northern hemisphere of AE Aqr is spotted.

The component masses of the cataclysmic variable V347 Puppis

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variable V347 Pup ( = LB 1800). There is evidence of irradiation on the inner hemisphere of the secondary star, which we correct for using a model to give a secondary-star radial velocity of K R = 198 +/- 5 km s(-1). The rotational velocity of the secondary star in V347 Pup is found to be v sin i = 131 +/- 5 km s(-1) and the system inclination is i = 84degrees.0 +/- 2degrees.3. From these parameters we obtain masses of M-1 = 0.63 +/- 0.04 M for the white dwarf primary and M-2 = 0.52 +/- 0.06 M for the M0.5V secondary star, giving a mass ratio of q = 0.83 +/- 0.05. On the basis of the component masses, and the spectral type and radius of the secondary star in V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows many of the characteristics of the SW Sex stars, exhibiting single-peaked emission lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to be close to the maximum allowed in a pressureless disc.

The masses of the cataclysmic variables AC Cancri and V363 Aurigae

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variables AC Cnc and V363 Aur (=Lanning 10). There is evidence of irradiation on the inner hemisphere of the secondary star in both systems, which we correct for using a model that reproduces the observations remarkably well. We find the radial velocity of the secondary star in AC Cnc to be K-R=176+/-3 km s(-1) and its rotational velocity to be v sin i=135+/-3 km s(-1). From these parameters we obtain masses of M-1=0.76+/-0.03 M-circle dot for the white-dwarf primary and M-2=0.77+/-0.05 M-circle dot for the K2+/-1 V secondary star, giving a mass ratio of q=1.02+/-0.04. We measure the radial and rotational velocities of the G7+/-2V secondary star in V363 Aur to be K-R=168+/-5 km s(-1) and v sin i=143+/-5 km s(-1), respectively. The component masses of V363 Aur are M-1=0.90+/-0.06M(circle dot) and M-2=1.06+/-0.11 M-circle dot giving a mass ratio of q=1.17+/-0.07. The mass ratios for AC Cnc and V363 Aur fall within the theoretical limits for dynamically and thermally stable mass transfer. Both systems are similar to the SW Sex stars, exhibiting single-peaked emission lines with transient absorption features, high-velocity S-wave components and phase-offsets in their radial-velocity curves. The Balmer lines in V363 Aur show a rapid increase in flux around phase 0 followed by a rapid decrease, which we attribute to the eclipse of an optically thick region at the centre of the disc. This model could also account for the behaviour of other SW Sex stars where the Balmer lines show only a shallow eclipse compared to the continuum.

Roche tomography of cataclysmic variables - I. Artefacts and techniques

Roche tomography is a technique used for imaging the Roche-lobe-filling secondary stars in cataclysmic variables (CVs). In order to interpret Roche tomograms correctly, one must determine whether features in the reconstruction are real, or the result of statistical or systematic errors. We explore the effects of systematic errors using reconstructions of simulated data sets, and show that systematic errors result in characteristic distortions of the final reconstructions that can be identified and corrected. In addition, we present a new method of estimating statistical errors on tomographic reconstructions using a Monte Carlo bootstrapping algorithm, and show this method to be much more reliable than Monte Carlo methods which 'jiggle' the data points in accordance with the size of their error bars.

Large molecules in the envelope surrounding IRC+10216

A new chemical model of the circumstellar envelope surrounding the carbon-rich star IRC+10216 is developed that includes carbon-containing molecules with up to 23 carbon atoms. The model consists of 3851 reactions involving 407 gas-phase species. Sizeable abundances of a variety of large molecules - including carbon clusters, unsaturated hydrocarbons and cyanopolyynes - have been calculated. Negative molecular ions of chemical formulae C-n(-) and CnH- (7 less than or equal to n less than or equal to 23) exist in considerable abundance, with peak concentrations at distances from the central star somewhat greater than their neutral counterparts. The negative ions might be detected in radio emission, or even in the optical absorption of background field stars. The calculated radial distributions of the carbon-chain CnH radicals are looked at carefully and compared with interferometric observations.

The evolution of solar systems

The recent detection of extra-solar planets and the increasing ability of modern instruments to image discs around young stars has brought a renewed interest in the formation of solar systems. In this article, I shall briefly review what we know about extra-solar planets and the physical structure of protostellar discs. One of the most powerful means of studying these disc systems is to observe the rotational line emission from molecules which can give information on physics and dynamics. At present, the observations are relatively crude but future instruments should be able to resolve molecular structures in the disc around nearby stars. As a prelude to these observations, I discuss some conceptually simple, although numerically challenging, models of the physical and chemical processes involved in determining the molecular distributions.

WASP-10b: a 3M(J), gas-giant planet transiting a late-type K star

We report the discovery of WASP-10b, a new transiting extrasolar planet (ESP) discovered by the Wide Angle Search for Planets ( WASP) Consortium and confirmed using Nordic Optical Telescope FIbre-fed Echelle Spectrograph and SOPHIE radial velocity data. A 3.09-d period, 29 mmag transit depth and 2.36 h duration are derived for WASP-10b using WASP and high-precision photometric observations. Simultaneous fitting to the photometric and radial velocity data using a Markov Chain Monte Carlo procedure leads to a planet radius of 1.28R(J), a mass of 2.96M(J) and eccentricity of approximate to 0.06. WASP-10b is one of the more massive transiting ESPs, and we compare its characteristics to the current sample of transiting ESP, where there is currently little information for masses greater than approximate to 2M(J) and non-zero eccentricities. WASP-10's host star, GSC 2752-00114 (USNO-B1.0 1214-0586164) is among the fainter stars in the WASP sample, with V = 12.7 and a spectral type of K5. This result shows promise for future late-type dwarf star surveys.

Polarimetric signatures of infrared jets in X-ray binaries

We present near-infrared linear spectropolarimetry of a sample of persistent X-ray binaries, Sco X-1, Cyg X-2 and GRS 1915+105. For Sco X-1 and Cyg X-2, the polarization levels at 2.4 µm are 1.3+/-0.10% and 5.4+/-0.7%, respectively, which is greater than the polarization level at 1.65 µm. This cannot be explained by interstellar polarization or electron scattering in the anisotropic environment of the accretion flow. We propose that the most likely explanation is that this is the polarimetric signature of synchrotron emission arising from close to the base of the jet. For Sco X-1 the position angle of the radio jet on the sky is approximately perpendicular to the near-infrared position angle (electric vector), suggesting that the magnetic field is aligned with the jet. These observations may be a first step towards probing the ordering, alignment, and variability of the outflow magnetic field, in a region closer to the central accreting object than is observed in the radio band.

A SuperWASP search for additional transiting planets in 24 known systems

We present results from a search for additional transiting planets in 24 systems already known to contain a transiting planet. We model the transits due to the known planet in each system and subtract these models from light curves obtained with the SuperWASP (Wide Angle Search for Planets) survey instruments. These residual light curves are then searched for evidence of additional periodic transit events. Although we do not find any evidence for additional planets in any of the planetary systems studied, we are able to characterize our ability to find such planets by means of Monte Carlo simulations. Artificially generated transit signals corresponding to planets with a range of sizes and orbital periods were injected into the SuperWASP photometry and the resulting light curves searched for planets. As a result, the detection efficiency as a function of both the radius and orbital period of any second planet is calculated. We determine that there is a good (>50 per cent) chance of detecting additional, Saturn-sized planets in P ~ 10 d orbits around planet-hosting stars that have several seasons of SuperWASP photometry. Additionally, we confirm previous evidence of the rotational stellar variability of WASP-10, and refine the period of rotation. We find that the period of the rotation is 11.91 +/- 0.05 d, and the false alarm probability for this period is extremely low (~10-13).

An orbital period of 0.94days for the hot-Jupiter planet WASP-18b

The `hot Jupiters' that abound in lists of known extrasolar planets are thought to have formed far from their host stars, but migrate inwards through interactions with the proto-planetary disk from which they were born, or by an alternative mechanism such as planet-planet scattering. The hot Jupiters closest to their parent stars, at orbital distances of only ~0.02 astronomical units, have strong tidal interactions, and systems such as OGLE-TR-56 have been suggested as tests of tidal dissipation theory. Here we report the discovery of planet WASP-18b with an orbital period of 0.94days and a mass of ten Jupiter masses (10MJup), resulting in a tidal interaction an order of magnitude stronger than that of planet OGLE-TR-56b. Under the assumption that the tidal-dissipation parameter Q of the host star is of the order of 106, as measured for Solar System bodies and binary stars and as often applied to extrasolar planets, WASP-18b will be spiralling inwards on a timescale less than a thousandth that of the lifetime of its host star. Therefore either WASP-18 is in a rare, exceptionally short-lived state, or the tidal dissipation in this system (and possibly other hot-Jupiter systems) must be much weaker than in the Solar System.

The main-sequence rotation-colour relation in the Coma Berenices open cluster

We present the results of a photometric survey of rotation rates in the Coma Berenices (Melotte 111) open cluster, using data obtained as part of the SuperWASP exoplanetary transit-search programme. The goal of the Coma survey was to measure precise rotation periods for main-sequence F, G and K dwarfs in this intermediate-age (~600 Myr) cluster, and to determine the extent to which magnetic braking has caused the stellar spin periods to converge. We find a tight, almost linear relationship between rotation period and J - K colour with an rms scatter of only 2 per cent. The relation is similar to that seen among F, G and K stars in the Hyades. Such strong convergence can only be explained if angular momentum is not at present being transferred from a reservoir in the deep stellar interiors to the surface layers. We conclude that the coupling time-scale for angular momentum transport from a rapidly spinning radiative core to the outer convective zone must be substantially shorter than the cluster age, and that from the age of Coma onwards stars rotate effectively as solid bodies. The existence of a tight relationship between stellar mass and rotation period at a given age supports the use of stellar rotation period as an age indicator in F, G and K stars of Hyades age and older. We demonstrate that individual stellar ages can be determined within the Coma population with an internal precision of the order of 9 per cent (rms), using a standard magnetic braking law in which rotation period increases with the square root of stellar age. We find that a slight modification to the magnetic-braking power law, P ~ t0.56, yields rotational and asteroseismological ages in good agreement for the Sun and other stars of solar age for which p-mode studies and photometric rotation periods have been published.

Atmospheric parameters and rotational velocities for a sample of Galactic B-type supergiants

High-resolution optical spectra of 57 Galactic B-type supergiant stars have been analysed to determine their rotational and macroturbulent velocities. In addition, their atmospheric parameters (effective temperature, surface gravity and microturbulent velocity) and surface nitrogen abundances have been estimated using a non-local thermodynamic equilibrium grid of model atmospheres. Comparisons of the projected rotational velocities have been made with the predictions of stellar evolutionary models and in general good agreement was found. However, for a small number of targets, their observed rotational velocities were significantly larger than predicted, although their nitrogen abundances were consistent with the rest of the sample. We conclude that binarity may have played a role in generating their large rotational velocities. No correlation was found between nitrogen abundances and the current projected rotational velocities. However, a correlation was found with the inferred projected rotational velocities of the main-sequence precursors of our supergiant sample. This correlation is again in agreement with the predictions of single star evolutionary models that incorporate rotational mixing. The origin of the macroturbulence and microturbulent velocity fields is discussed and our results support previous theoretical studies that link the former to subphotospheric convection and the latter to non-radial gravity mode oscillations. In addition, we have attempted to identify differential rotation in our most rapidly rotating targets.

Li-6 in the atmosphere of GJ 117 revisited

Detection of Li-6 has been shown for energetic solar events, one chromospherically active binary, and several dwarf halo stars. We had previously found a Li-6/Li-7 = 0.03 +/- 0.01 for active K dwarf GJ 117 using VLT UVES observations. Here we present high signal-to-noise (> 1000) high spectral resolution observations taken with the McDonald Observatory's 2.7 m and echelle spectrometer of GJ 117. We have used the solar spectrum and template stars to eliminate possible blends, such as Ti I, in the Li-6 spectral region. Our new analysis, using an updated PHOENIX model atmosphere, finds Li-6/Li-7 = 0.05 +/- 0.02. In addition, bisector analysis showed no significant red asymmetries that would affect the lithium line profile. No changes above the statistical uncertainties are found between the VLT and McDonald data. The amount of Li-6 derived for GJ 117 is consistent with creation in spallation reactions on the stellar surface, but we caution that uncertainties in the continuum level may cause additional uncertainty in the Li-6 fraction.

Line-profile tomography of exoplanet transits - II. A gas-giant planet transiting a rapidly rotating A5 star

Most of our knowledge of extrasolar planets rests on precise radial-velocity measurements, either for direct detection or for confirmation of the planetary origin of photometric transit signals. This has limited our exploration of the parameter space of exoplanet hosts to solar- and later-type, sharp-lined stars. Here we extend the realm of stars with known planetary companions to include hot, fast-rotating stars. Planet-like transits have previously been reported in the light curve obtained by the SuperWASP survey of the A5 star HD15082 (WASP-33 V = 8.3, v sini = 86 km s-1). Here we report further photometry and time-series spectroscopy through three separate transits, which we use to confirm the existence of a gas-giant planet with an orbital period of 1.22d in orbit around HD15082. From the photometry and the properties of the planet signal travelling through the spectral line profiles during the transit, we directly derive the size of the planet, the inclination and obliquity of its orbital plane and its retrograde orbital motion relative to the spin of the star. This kind of analysis opens the way to studying the formation of planets around a whole new class of young, early-type stars, hence under different physical conditions and generally in an earlier stage of formation than in sharp-lined late-type stars. The reflex orbital motion of the star caused by the transiting planet is small, yielding an upper mass limit of 4.1MJupiter on the planet. We also find evidence of a third body of substellar mass in the system, which may explain the unusual orbit of the transiting planet. In HD 15082, the stellar line profiles also show evidence of non-radial pulsations, clearly distinct from the planetary transit signal. This raises the intriguing possibility that tides raised by the close-in planet may excite or amplify the pulsations in such stars.

Extensive optical and near-infrared observations of the nearby, narrow-lined type Ic SN 2007gr: days 5 to 415

We present photometric and spectroscopic observations at optical and near-infrared wavelengths of the nearby type Ic supernova 2007gr. These represent the most extensive data-set to date of any supernova of this sub-type, with frequent coverage from shortly after discovery to more than one year post-explosion. We deduce a rise time to B-band maximum of 11.5 +/- 2.7 d. We find a peak B-band magnitude of M-B = -16.8, and light curves which are remarkably similar to the so-called "hypernova" SN 2002ap. In contrast, the spectra of SNe 2007gr and 2002ap show marked differences, not least in their respective expansion velocities. We attribute these differences primarily to the density profiles of their progenitor stars at the time of explosion i.e. a more compact star for SN 2007gr compared to SN 2002ap. From the quasi-bolometric light curve of SN 2007gr, we estimate that 0.076 +/- 0.010 M-circle dot of Ni-56 was produced in the explosion. Our near-infrared (IR) spectra clearly show the onset and disappearance of the first overtone of carbon monoxide (CO) between similar to 70 to 175 d relative to B-band maximum. The detection of the CO molecule implies that ionised He was not microscopically mixed within the carbon/oxygen layers. From the optical spectra, near-IR light curves, and colour evolution, we find no evidence for dust condensation in the ejecta out to about +400 d. Given the combination of unprecedented temporal coverage, and high signal-to-noise data, we suggest that SN 2007gr could be used as a template object for supernovae of this sub-class.