A Bayesian periodogram finds evidence for three planets in HD 11964

A Bayesian multiplanet Kepler periodogram has been developed for the analysis of precision radial velocity data. The periodogram employs a parallel tempering Markov chain Monte Carlo algorithm. The HD 11964 data have been re-analysed using 1, 2, 3 and 4 planet models. Assuming that all the models are equally probable a priori, the three planet model is found to be ≥600 times more probable than the next most probable model which is a two planet model. The most probable model exhibits three periods of  38.02^+0.06_−0.05, 360⁺⁴₋₄ and 1924⁺⁴⁴₋₄₃ d  , and eccentricities of  0.22^+0.11_−0.22, 0.63^+0.34_−0.17 and 0.05^+0.03_−0.05  , respectively. Assuming the three signals (each one consistent with a Keplerian orbit) are caused by planets, the corresponding limits on planetary mass ( M sin  i ) and semimajor axis are     respectively. The small difference (1.3σ) between the 360-d period and one year suggests that it might be worth investigating the barycentric correction for the HD 11964 data.     

WASP-14b: 7.3 MJ transiting planet in an eccentric orbit

We report the discovery of a 7.3 M _J exoplanet WASP-14b, one of the most massive transiting exoplanets observed to date. The planet orbits the 10th-magnitude F5V star USNO-B1 11118−0262485 with a period of 2.243 752 d and orbital eccentricity   e = 0.09  . A simultaneous fit of the transit light curve and radial velocity measurements yields a planetary mass of 7.3 ± 0.5 M _J and a radius of 1.28 ± 0.08 R _J. This leads to a mean density of about 4.6 g cm⁻³ making it the densest transiting exoplanets yet found at an orbital period less than 3 d. We estimate this system to be at a distance of  160 ± 20  pc. Spectral analysis of the host star reveals a temperature of  6475 ± 100 K, log  g = 4.07  cm s⁻² and   v sin  i = 4.9 ± 1.0  km s⁻¹, and also a high lithium abundance,  log  N (Li) = 2.84 ± 0.05  . The stellar density, effective temperature and rotation rate suggest an age for the system of about 0.5–1.0 Gyr.     

WASP-3b: a strongly irradiated transiting gas-giant planet

We report the discovery of WASP-3b, the third transiting exoplanet to be discovered by the WASP and SOPHIE collaboration. WASP-3b transits its host star USNO-B1.0 1256−0285133 every  1.846 834 ± 0.000 002  d. Our high-precision radial velocity measurements present a variation with amplitude characteristic of a planetary-mass companion and in phase with the light curve. Adaptive optics imaging shows no evidence for nearby stellar companions, and line-bisector analysis excludes faint, unresolved binarity and stellar activity as the cause of the radial velocity variations. We make a preliminary spectroscopic analysis of the host star and find it to have   T _eff= 6400 ± 100 K  and  log   g = 4.25 ± 0.05  which suggests it is most likely an unevolved main-sequence star of spectral type F7-8V. Our simultaneous modelling of the transit photometry and reflex motion of the host leads us to derive a mass of  1.76^+0.08_−0.14 M _J  and radius  1.31^+0.07_−0.14 R _J  for WASP-3b. The proximity and relative temperature of the host star suggests that WASP-3b is one of the hottest exoplanets known, and thus has the potential to place stringent constraints on exoplanet atmospheric models.     

Debris disc stirring by secular perturbations from giant planets

Detectable debris discs are thought to require dynamical excitation ('stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semimajor axes to precess at different rates, and after some time   t _cross  initially non-intersecting orbits begin to cross. We show that   t _cross∝ a ^9/2_disc/( m _pl e _pl a ³_pl)  , where   m _pl, e _pl  and   a _pl  are the mass, eccentricity and semimajor axis of the planet, and   a _disc  is the semimajor axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced among planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100 s of au. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star ε Eridani, we find that the known radial velocity planet can excite the planetesimal belt at 60 au sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.     

Expected planets in globular clusters

We argue that all transient searches for planets in globular clusters have a very low detection probability. Planets of low-metallicity stars typically do not reside at small orbital separations. The dependence of planetary system properties on metallicity is clearly seen when the quantity   I _e ≡ M _p[ a (1 − e )]²  is considered;   M _p, a   and e are the planet mass, semimajor axis and eccentricity, respectively. In high-metallicity systems, there is a concentration of systems at high and low values of I _e , with a low-populated gap near   I _e ∼ 0.3 M _J au²  , where M _J is Jupiter's mass. In low-metallicity systems, the concentration is only at the higher range of I _e , with a tail to low values of I _e . Therefore, it is still possible that planets exist around main-sequence stars in globular clusters, although at small numbers because of the low metallicity, and at orbital periods of ≳10 d. We discuss the implications of our conclusions on the role that companions can play in the evolution of their parent stars in globular clusters, for example, influencing the distribution of horizontal branch stars on the Hertzsprung–Russell diagram of some globular clusters, and in forming low-mass white dwarfs.     

Why are there so few hot Jupiters?

We use numerical simulations to model the migration of massive planets at small radii and compare the results with the known properties of 'hot Jupiters' (extrasolar planets with semimajor axes   a < 0.1  au). For planet masses   M _pl sin  i > 0.5 M _J  , the evidence for any 'pile-up' at small radii is weak (statistically insignificant), and although the mass function of hot Jupiters is deficient in high-mass planets as compared to a reference sample located further out, the small sample size precludes definitive conclusions. We suggest that these properties are consistent with disc migration followed by entry into a magnetospheric cavity close to the star. Entry into the cavity results in a slowing of migration, accompanied by a growth in orbital eccentricity. For planet masses in excess of 1 Jupiter mass we find eccentricity growth time-scales of a few ×10⁵ yr, suggesting that these planets may often be rapidly destroyed. Eccentricity growth appears to be faster for more massive planets which may explain changes in the planetary mass function at small radii and may also predict a pile-up of lower mass planets, the sample of which is still incomplete.     

Spectroscopic solution of the star QX Cas

High-resolution spectroscopic observations around the Hα line of the binary star QX Cas covering the whole orbital period are presented. Our radial velocity solution, the first ever determined, requires an eccentric orbit with the following orbital parameters: eccentricity,   e = 0.22 ± 0.01  ; longitude of periastron,  ω= 45°± 5°  ; semi-amplitudes of the radial velocity curves of the primary and secondary stars,   K ₁ sin  i = 125.8 ± 0.9 km s⁻¹  and   K ₂ sin  i = 144.8 ± 1.1 km s⁻¹  ; gamma velocity,   V ₀= 65.1 ± 0.5 km s⁻¹  ; and mass ratio,   q = 0.869 ± 0.013  . The corresponding lower limits of the masses of the components and their separation are         , and   a sin  i = 31.34 ± 0.48 R_⊙  .     

A dynamical analysis of the 14 Herculis planetary system

Precision radial velocity measurements of the Sun-like dwarf 14 Herculis published by Naef et al., Butler et al. and Wittenmyer, Endl & Cochran reveal a Jovian planet in a 1760-d orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant planet. We derive dynamical limits to the orbital parameters of the putative outer Jovian companion in an orbit within ∼13 au. In this case, the mutual interactions between the Jovian planets are important for the long-term stability of the system. The best self-consistent and stable Newtonian fit to an edge-on configuration of Jovian planets has the outer planet in 9-au orbit with a moderate eccentricity of ∼0.2 and confined to a zone spanned by the low-order mean motion resonances 5:1 and 6:1. This solution lies in a shallow minimum of (χ²_ν)^1/2 and persists over a wide range of the system inclination. Other stable configurations within 1σ confidence interval of the best fit are possible for the semimajor axis of the outer planet in the range of (6,13) au and the eccentricity in the range of (0, 0.3). The orbital inclination cannot yet be determined but when it decreases, both planetary masses approach ∼10 m _J and for i ∼ 30° the hierarchy of the masses is reversed.     

An iterative filter to reconstruct planetary transit signals in the presence of stellar variability

We report the identification, from a photometric, astrometric and spectroscopic study, of a massive white dwarf member of the nearby, approximately solar metallicity, Coma Berenices open star cluster (Melotte 111). We find the optical to near-infrared energy distribution of WD 1216+260 to be entirely consistent with that of an isolated DA and determine the effective temperature and surface gravity of this object to be   T _eff= 15 739⁺¹⁹⁷₋₁₉₆ K  and  log  g = 8.46^+0.03_−0.02  . We set tight limits on the mass of a putative cool companion,   M ≳ 0.036 M_⊙  (spatially unresolved) and   M ≳ 0.034 M_⊙  (spatially resolved and   a ≲ 2500 au  ). Based on the predictions of CO core, thick H layer evolutionary models we determine the mass and cooling time of WD 1216+260 to be   M _WD= 0.90 ± 0.04 M_⊙  and  τ_cool= 363⁺⁴⁶₋₄₁ Myr  , respectively. For an adopted cluster age of  τ= 500 ± 100 Myr  we infer the mass of its progenitor star to be   M _init= 4.77^+5.37_−0.97 M_⊙  . We briefly discuss this result in the context of the form of the stellar initial mass–final mass relation.     

The mass ratio and the orbital parameters of the sdOB binary AA Doradus

The time sequence of 105 spectra covering one full orbital period of AA Dor has been analysed. Direct determination of   V  sin  i   for the sdOB component from 97 spectra outside of the eclipse for the lines Mg  ii 4481 Å and Si  iv 4089 Å clearly indicated a substantially smaller value than estimated before. Detailed modelling of line-profile variations for eight spectra during the eclipse for the Mg  ii 4481 Å line, combined with the out-of-eclipse fits, gave   V  sin  i = 31.8 ± 1.8 km s⁻¹  . The previous determinations of   V  sin  i   , based on the He  ii 4686 Å line, appear to be invalid because of the large natural broadening of the line. With the assumption of the solid-body, synchronous rotation of the sdOB primary, the measured values of the semi-amplitude K ₁ and   V  sin  i   lead to the mass ratio   q = 0.213 ± 0.013  which in turn gives K ₂ and thus the masses and radii of both components. The sdOB component appears to be less massive than assumed before,   M ₁= 0.25 ± 0.05 M_⊙  , but the secondary has its mass–radius parameters close to theoretically predicted for a brown dwarf,   M ₂= 0.054 ± 0.010 M_⊙  and   R ₂= 0.089 ± 0.005 R_⊙  . Our results do not agree with the recent determination of Vŭcković et al. based on a K ₂ estimate from line-profile asymmetries.     

Orbital parameters, masses and distance to β Centauri determined with the Sydney University Stellar Interferometer and high-resolution spectroscopy

The bright southern binary star β Centauri (HR 5267) has been observed with the Sydney University Stellar Interferometer (SUSI) and spectroscopically with the European Southern Observatory Coude Auxiliary Telescope and Swiss Euler telescope at La Silla. The interferometric observations have confirmed the binary nature of the primary component and have enabled the determination of the orbital parameters of the system. At the observing wavelength of 442 nm the two components of the primary system have a magnitude difference of  0.15 ± 0.02  . The combination of interferometric and spectroscopic data gives the following results: orbital period  357.00 ± 0.07 d  , semimajor axis  25.30 ± 0.19 mas  , inclination  674 ± 03  , eccentricity  0.821 ± 0.003  , distance  102.3 ± 1.7 pc  , primary and secondary masses   M ₁= M ₂= 9.1 ± 0.3 M_⊙  and absolute visual magnitudes of the primary and secondary   M _{1 V} =−3.85 ± 0.05  and   M _{2 V} =−3.70 ± 0.05  , respectively. The high degree of accuracy of the results offers a fruitful starting point for future asteroseismic modelling of the pulsating binary components.     

Spectroscopic and photometric observations of the eclipsing star UV Leo

High-resolution spectroscopic observations around the Hα line and BVRI photometry of the eclipsing short-period RS CVn star UV Leo are presented. The simultaneous light-curve solution and radial velocity-curve solution led to the following values of the global parameters of the binary: temperatures   T ₁= 6000 ± 100 K  and   T ₂= 5970 ± 20 K  ; masses   M ₁= 0.976 ± 0.067 M_⊙  and   M ₂= 0.931 ± 0.052 M_⊙  ; separation   a = 3.716 ± 0.048 R_⊙  ; orbital inclination     ; radii   R ₁= 1.115 ± 0.052 R_⊙  and   R ₂= 1.078 ± 0.051 R_⊙  ; equatorial velocities   V ₁= 98.8 ± 2.3 km s⁻¹  and   V ₂= 89.6 ± 2.7 km s⁻¹  . These results lead to the conclusion that the two components of UV Leo are slightly oversized for their masses and lie within the main-sequence band on the mass–radius diagram, close to the isochrone 9 × 10¹⁰ yr.     

Microlensing in dark matter haloes

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.     

Photospheric phosphorus in the FUSE spectra of GD71 and two similar DA white dwarfs

We report the detection, from the Far Ultraviolet Spectroscopic Explorer (FUSE) data, of phosphorus in the atmospheres of GD71 and two similar DA white dwarfs. This is the first detection of a trace metal in the photosphere of the spectrophotometric standard star GD71. Collectively, these objects represent the coolest DA white dwarfs in which photospheric phosphorus has been observed. We use a grid of homogeneous non-local thermodynamic equilibrium synthetic spectra to measure abundances of  [P/H]=−8.57^+0.09_−0.13, −8.70^+0.23_−0.37  and  −8.36^+0.14_−0.19  in GD71, RE J1918+595 and RE J0605−482 respectively. At the observed level we find that phosphorus has no significant impact on the overall energy distribution of GD71. We explore possible mechanisms responsible for the presence of this element in these stars, concluding that the most likely is an interplay between radiative levitation and gravitational settling, possibly modified by weak mass loss.     

A survey for planetary transits in the field of NGC 7789

We present results from 30 nights of observations of the open cluster NGC 7789 with the Wide Field Camera on the Isaac Newton Telescope, La Palma. From ∼900 epochs, we obtained light curves and Sloan   r '− i '  colours for ∼33 000 stars, with ∼2400 stars having better than 1 per cent precision. We expected to detect ∼2 transiting hot Jupiter planets if 1 per cent of stars host such a companion and a typical hot Jupiter radius is  ∼1.2  R _J  . We find 24 transit candidates, 14 of which we can assign a period. We rule out the transiting planet model for 21 of these candidates using various robust arguments. For two candidates, we are unable to decide on their nature, although it seems most likely that they are eclipsing binaries as well. We have one candidate exhibiting a single eclipse, for which we derive a radius of  1.81^+0.09_−0.00  R _J  . Three candidates remain that require follow-up observations in order to determine their nature.     

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⁻¹  . The rotational velocity of the secondary star in V347 Pup is found to be   v sin  i = 131 ± 5 km s⁻¹  and the system inclination is   i = 840 ± 23  . From these parameters we obtain masses of   M ₁= 0.63 ± 0.04 M_⊙  for the white dwarf primary and   M ₂= 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.     

Orbital solution and fundamental parameters of σ Scorpii

The first orbital solution for the spectroscopic pair in the multiple star system σ Scorpii, determined from measurements with the Sydney University Stellar Interferometer, is presented. The primary component is of β Cephei variable type and has been one of the most intensively studied examples of its class. The orbital solution, when combined with radial velocity results found in the literature, yields a distance of  174⁺²³₋₁₈ pc  , which is consistent with, but more accurate than the Hipparcos value. For the primary component we determine  18.4 ± 5.4 M_⊙, −4.12 ± 0.34 mag  and  12.7 ± 1.8 R_⊙  for the mass, absolute visual magnitude and radius, respectively. A B1 dwarf spectral type and luminosity class for the secondary is proposed from the mass determination of  11.9 ± 3.1 M_⊙  and the estimated system age of 10 Myr.     

Maximum-likelihood method for estimating the mass and period distributions of extrasolar planets

We investigate the distribution of mass M and orbital period P of extrasolar planets, taking account of selection effects caused by the limited velocity precision and duration of existing surveys. We fit the data on 72 planets to a power-law distribution of the form  d n = CM ^−α P ^−β(d M / M )(d P / P )  , and find  α= 0.11 ± 0.10  ,  β=−0.27 ± 0.06  for   M ≲ 10 M _J  , where   M _J  is the mass of Jupiter. The correlation coefficient between these two exponents is −0.31, indicating that uncertainties in the two distributions are coupled. We estimate that 4 per cent of solar-type stars have companions in the range  1 M _J < M < 10 M _J  ,  2_d < P < 10 yr  .     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Optical spectroscopy and photometry of SAX J1808.4−3658 in outburst

We present phase resolved optical spectroscopy and photometry of V4580 Sagittarii, the optical counterpart to the accretion powered millisecond pulsar SAX J1808.4−3658, obtained during the 2008 September/October outburst. Doppler tomography of the N  iii λ4640.64 Bowen blend emission line reveals a focused spot of emission at a location consistent with the secondary star. The velocity of this emission occurs at  324 ± 15 km s⁻¹  ; applying a ' K -correction', we find the velocity of the secondary star projected on to the line of sight to be  370 ± 40 km s⁻¹  . Based on existing pulse timing measurements, this constrains the mass ratio of the system to be  0.044^+0.005_−0.004  , and the mass function for the pulsar to be  0.44^+0.16_−0.13 M_⊙  . Combining this mass function with various inclination estimates from other authors, we find no evidence to suggest that the neutron star in SAX J1808.4−3658 is more massive than the canonical value of  1.4 M_⊙  . Our optical light curves exhibit a possible superhump modulation, expected for a system with such a low mass ratio. The equivalent width of the Ca  ii H and K interstellar absorption lines suggest that the distance to the source is ∼2.5 kpc. This is consistent with previous distance estimates based on type-I X-ray bursts which assume cosmic abundances of hydrogen, but lower than more recent estimates which assume helium-rich bursts.