Radial-velocity (RV) signals induce RV variations an order of magnitude larger than the signal created by the orbit of Earth-twins, thus preventing their detection. The goal of this paper is to compare the efficiency of the different methods used to deal with stellar signals to recover extremely low-mass planets despite. However, because observed RV variations at the m/s precision level or below is a combination of signals induced by unresolved orbiting planets, by the star, and by the instrument, performing such a comparison using real data is extremely challenging. To circumvent this problem, we generated simulated RV measurements including realistic stellar and planetary signals. Different teams analyzed blindly those simulated RV measurements, using their own method to recover planetary signals despite stellar RV signals. By comparing the results obtained by the different teams with the planetary and stellar parameters used to generate the simulated RVs, it is therefore possible to compare the efficiency of these different methods. The most efficient methods to recover planetary signals {take into account the different activity indicators,} use red-noise models to account for stellar RV signals and a Bayesian framework to provide model comparison in a robust statistical approach. Using the most efficient methodology, planets can be found down to K/N= K_pl/RV_rms*sqrt{N_obs}=5 with a threshold of K/N=7.5 at the level of 80-90% recovery rate found for a number of methods. These recovery rates drop dramatically for K/N smaller than this threshold. In addition, for the best teams, no false positives with K/N > 7.5 were detected, while a non-negligible fraction of them appear for smaller K/N. A limit of K/N = 7.5 seems therefore a safe threshold to attest the veracity of planetary signals for RV measurements with similar properties to those of the different RV fitting challenge systems.

The Second Workshop on Extreme Precision Radial Velocities defined circa 2015 the state of the art Doppler precision and identified the critical path challenges for reaching 10 cm/s measurement precision. The presentations and discussion of key issues for instrumentation and data analysis and the workshop recommendations for achieving this precision are summarized here. Beginning with the HARPS spectrograph, technological advances for precision radial velocity measurements have focused on building extremely stable instruments. To reach still higher precision, future spectrometers will need to produce even higher fidelity spectra. This should be possible with improved environmental control, greater stability in the illumination of the spectrometer optics, better detectors, more precise wavelength calibration, and broader bandwidth spectra. Key data analysis challenges for the precision radial velocity community include distinguishing center of mass Keplerian motion from photospheric velocities, and the proper treatment of telluric contamination. Success here is coupled to the instrument design, but also requires the implementation of robust statistical and modeling techniques. Center of mass velocities produce Doppler shifts that affect every line identically, while photospheric velocities produce line profile asymmetries with wavelength and temporal dependencies that are different from Keplerian signals. Exoplanets are an important subfield of astronomy and there has been an impressive rate of discovery over the past two decades. Higher precision radial velocity measurements are required to serve as a discovery technique for potentially habitable worlds and to characterize detections from transit missions. The future of exoplanet science has very different trajectories depending on the precision that can ultimately be achieved with Doppler measurements.

A new apodized Keplerian (AK) model is proposed for the analysis of precision radial velocity (RV) data to model both planetary and stellar activity (SA) induced RV signals. A symmetrical Gaussian apodization function with unknown width and center can distinguish planetary signals from SA signals on the basis of the span of the apodization window. The general model for $m$ apodized Keplerian signals includes a linear regression term between RV and the stellar activity diagnostic $\log(R'hk)$, as well as an extra Gaussian noise term with unknown standard deviation. The model parameters are explored using a Bayesian fusion MCMC code. A differential version of the Generalized Lomb-Scargle periodogram that employs a control diagnostic provides an additional way of distinguishing SA signals and helps guide the choice of new periods. Results are reported for a recent international RV blind challenge which included multiple state of the art simulated data sets supported by a variety of stellar activity diagnostics. In the current implementation, the AK method achieved a reduction in SA noise by a factor of approximately 6. Final parameter estimates for the planetary candidates are derived from fits that include AK signals to model the SA components and simple Keplerians to model the planetary candidates. Preliminary results are also reported for AK models augmented by a moving average component that allows for correlations in the residuals.

A re-analysis of Gliese 667C HARPS precision radial velocity data was carried out with a Bayesian multi-planet Kepler periodogram (from 0 to 7 planets) based on a fusion Markov chain Monte Carlo algorithm. The most probable number of signals detected is 6 with a Bayesian false alarm probability of 0.012. The residuals are shown to be consistent with white noise. The 6 signals detected include two previously reported with periods of 7.198 (b) and 28.14 (c) days, plus additional periods of 30.82 (d), 38.82 (e), 53.22, and 91.3 (f) days. The 53 day signal is probably the second harmonic of the stellar rotation period and is likely the result of surface activity. The existence of the additonal Keplerian signals suggest the possibilty of further planets, two of which (d and e) could join Gl 667Cc in the central region of the habitable zone. N-body simulations are required to determine which of these signals are consistent with a stable planetary system. $M \sin i$ values corresponding to signals b, c, d, e, and f are $\sim$ 5.4, 4.8, 3.1, 2.4, and 5.4 M$_{\earth}$, respectively.

A re-analysis of Gliese 581 HARPS and HIRES precision radial velocity data was carried out with a Bayesian multi-planet Kepler periodogram (from 1 to 6 planets) based on a fusion Markov chain Monte Carlo algorithm. In all cases the analysis included an unknown parameterized stellar jitter noise term. For the HARPS data set the most probable number of planetary signals detected is 5 with a Bayesian false alarm probability of 0.01. These include the $3.1498\pm0.0005$, $5.3687\pm0.0002$, $12.927_{-0.004}^{+0.006}$, and $66.9\pm0.2$d periods reported previously plus a $399_{-16}^{+14}$d period. The orbital eccentricities are $0.0_{-0.0}^{+0.2}$, $0.00_{-0.00}^{+0.02}$, $0.10_{-0.10}^{+0.06}$, $0.33_{-0.10}^{+0.09}$, and $0.02_{-0.02}^{+0.30}$, respectively. The semi-major axis and $M sin i$ of the 5 planets are ($0.0285\pm0.0006$ au, $1.9\pm0.3$M$_{\earth}$), ($0.0406\pm0.0009$ au, $15.7\pm0.7$M$_{\earth}$), ($0.073\pm0.002$ au, $5.3\pm0.4$M$_{\earth}$), ($0.218\pm0.005$ au, $6.7\pm0.8$M$_{\earth}$), and ($0.7\pm0.2$ au, $6.6_{-2.7}^{+2.0}$M$_{\earth}$), respectively. The analysis of the HIRES data set yielded a reliable detection of only the strongest 5.37 and 12.9 day periods. The analysis of the combined HIRES/HARPS data again only reliably detected the 5.37 and 12.9d periods. Detection of 4 planetary signals with periods of 3.15, 5.37, 12.9, and 66.9d was only achieved by including an additional unknown but parameterized Gaussian error term added in quadrature to the HIRES quoted errors. The marginal distribution for the sigma of this additional error term has a well defined peak at $1.8\pm0.4$m s$^{-1}$. It is possible that this additional error arises from unidentified systematic effects. We did not find clear evidence for a fifth planetary signal in the combined HIRES/HARPS data set.

A Bayesian analysis of 47 Ursae Majoris (47 UMa) radial velocity data confirms and refines the properties of two previously reported planets with periods of 1079 and 2325 days and finds evidence for an additional long period planet with a period of approximately 10000 days. The three planet model is found to be 10^5 times more probable than the next most probable model which is a two planet model. The nonlinear model fitting is accomplished with a new hybrid Markov chain Monte Carlo (HMCMC) algorithm which incorporates parallel tempering, simulated annealing and genetic crossover operations. Each of these features facilitate the detection of a global minimum in chi-squared. By combining all three, the HMCMC greatly increases the probability of realizing this goal. When applied to the Kepler problem it acts as a powerful multi-planet Kepler periodogram. The measured periods are 1078 \pm 2, 2391{+100}{-87}, and 14002{+4018}{-5095}d, and the corresponding eccentricities are 0.032 \pm 0.014, 0.098{+.047}{-.096}, and 0.16{+.09}{-.16}. The results favor low eccentricity orbits for all three. 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 semi-major axis are (2.53{+.07}{-.06}MJ, 2.10\pm0.02au), (0.54\pm0.07MJ, 3.6\pm0.1au), and (1.6{+0.3}{-0.5}MJ, 11.6{+2.1}{-2.9}au), respectively. We have also characterized a noise induced eccentricity bias and designed a correction filter that can be used as an alternate prior for eccentricity, to enhance the detection of planetary orbits of low or moderate eccentricity.

A Bayesian re-analysis of published radial velocity data sets is providing evidence for additional planetary candidates. The nonlinear model fitting is accomplished with a new hybrid Markov chain Monte Carlo (HMCMC) algorithm which incorporates parallel tempering, simulated annealing and genetic crossover operations. Each of these features facilitate the detection of a global minimum in chi^2. By combining all three, the HMCMC greatly increases the probability of realizing this goal. When applied to the Kepler problem it acts as a powerful multi-planet Kepler periodogram for both parameter estimation and model selection. The HMCMC algorithm is embedded in a unique two stage adaptive control system that automates the tuning of the MCMC proposal distributions through an annealing operation.

A Bayesian multi-planet Kepler periodogram has been developed for the analysis of precision radial velocity data (Gregory 2005b and 2007). The periodogram employs a parallel tempering Markov chain Monte Carlo algorithm. The HD 11964 data (Butler et al. 2006) has been re-analyzed 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+-4 and 1924+44-43 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 semi-major axis are 0.090+0.15-0.14 M_J, 0.253+-0.009 au, 0.21+0.06-0.07 M_J, 1.13+-0.04 au, 0.77+-0.08 M_J, 3.46+-0.13 au, respectively. The small difference (1.3 sigma) between the 360 day period and one year suggests that it might be worth investigating the barycentric correction for the HD 11964 data.

An automatic Bayesian Kepler periodogram has been developed for identifying and characterizing multiple planetary orbits in precision radial velocity data. The periodogram is powered by a parallel tempering MCMC algorithm which is capable of efficiently exploring a multi-planet model parameter space. The periodogram employs an alternative method for converting the time of an observation to true anomaly that enables it to handle much larger data sets without a significant increase in computation time. Improvements in the periodogram and further tests using data from HD 208487 have resulted in the detection of a second planet with a period of $909_{-92}^{82}$d, an eccentricity of $0.37_{-0.20}^{0.26}$, a semi-major axis of $1.87_{-0.14}^{0.13}$ AU and an $M \sin i = 0.45_{-0.11}^{0.13}$ M$_{\rm J}$. The revised parameters of the first planet are period $= 129.8\pm 0.4$d, eccentricity $= 0.20\pm 0.09$, semi-major axis $= 0.51\pm0.02$ AU and $M \sin i = 0.41\pm0.05$ M$_{\rm J}$. Particular attention is paid to several methods for calculating the model marginal likelihood which is used to compare the probabilities of models with different numbers of planets.

Radial velocity (RV) planet searches are increasingly finding planets with small velocity amplitudes, with long orbital periods, or in multiple planet systems. Bayesian inference has the potential to improve the interpretation of existing observations, the planning of future observations and ultimately inferences concerning the overall population of planets. In recent years, the refinement of Markov chain Monte Carlo (MCMC) algorithms has made it practical to accurately characterize orbital parameters and their uncertainties from RV observations of single-planet and weakly interacting multiple-planet systems. Unfortunately, MCMC is not sufficient for Bayesian model selection, i.e., comparing the marginal posterior probability of models, as is necessary to determine how strongly the observational data favor a model with n+1 planets over a model with just n planets. Many of the obvious estimators for the marginal posterior probability suffer from poor convergence properties. We compare several estimators of the marginal likelihood and feature those that display desirable convergence properties based on the analysis of a sample data set for HD 88133b. We find that methods based on importance sampling are most efficient, provided that a good analytic approximation of the posterior probability distribution is available. We present a simple algorithm for using a sample from the posterior to construct a mixture distribution that approximates the posterior and can be used for importance sampling and Bayesian model selection. We conclude with some suggestions for the development and refinement of computationally efficient and robust estimators of marginal posterior probabilities.

Precision radial velocity data for HD 208487 has been re-analyzed using a new Bayesian multi-planet Kepler periodogram. The periodgram employs a parallel tempering Markov chain Monte Carlo algorithm with a novel statistical control system. We confirm the previously reported orbit (Tinney et al. 2005) of 130 days. In addition, we conclude there is strong evidence for a second planet with a period of 998 -62 +57 days, an eccentricity of 0.19 -0.18 +0.05, and an M sin i = 0.46 -0.13 +0.05 of Jupiter's mass.

We report the discovery of three binary millisecond pulsars during the Parkes Multibeam Pulsar Survey of the Galactic Plane. The objects are highly recycled and are in orbits of many tens of days about low-mass white-dwarf companions. The eccentricity of one object, PSR J1853+1303, is more than an order of magnitude lower than predicted by the theory of convective fluctuations during tidal circularization. We demonstrate that, under the assumption that the systems are randomly oriented, current theoretical models of the core-mass--orbital-period relation for the progenitors of these systems likely overestimate the white-dwarf masses, strengthening previous concerns about the match of these models to the data. The new objects allow us to update the limits on violation of relativistic equivalence principles to 95% confidence upper limits of 5.6 x 10^-3 for the Strong Equivalence Principle parameter Delta and 4.0 x 10^-20 for the Lorentz-invariance/momentum-conservation parameter alpha_3.

The Jodrell Bank - VLA Astrometric Survey (JVAS) of flat-spectrum sources yielded a catalog of 2118 compact radio sources in the northern sky. Those sources are being used as phase calibrators for many synthesis arrays. JVAS suffers from a zone of avoidance because Galactic confusion prevented selection of flat-spectrum sources at $\mid b \mid < 2.5\arcdeg$. This confusion problem was overcome by selecting variable GB6 sources in the JVAS zone of avoidance. A subset of these sources were observed at 8.5 GHz with the NRAO VLA in the JVAS style, leading to 29 compact radio sources in the zone of avoidance whose positions have been measured to an rms accuracy of about 10~mas or less. This extension of the JVAS lists to lower Galactic latitudes (1) improves the prospects for VLBA phase-referencing of Galactic targets and (2) strengthens the lists' usefulness for studies of the Galactic interstellar medium.