Doppler spectroscopy is a powerful tool for discovering and characterizing exoplanets. For decades, the standard approach to extracting radial velocities (RVs) has been to cross-correlate observed spectra with a weighted template mask. While still widely used, this approach is known to suffer numerous drawbacks, and so in recent years increasing attention has been paid to developing new and improved ways of extracting RVs. In this proof-of-concept paper we present a simple yet powerful approach to RV extraction. We use Gaussian processes to model and align all pairs of spectra with each other; without constructing a template, we combine pairwise RV shifts to produce accurate differential stellar RVs. Doing this on a highly-localized basis enables a data-driven approach to identifying and mitigating spectral contamination, even without the input of any prior astrophysical knowledge. We show that a crude implementation of this method applied to an inactive standard star yields RVs with comparable precision to and significantly lower rms variation than RVs from industry-standard pipelines. Though amenable to numerous improvements, even in its basic form presented here our method could facilitate the study of smaller planets around a wider variety of stars than has previously been possible.

Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) with M_* = 0.892 +/- (0.051, 0.035) M_sun and R_* = 0.8717 +/- (0.0064, 0.0061) R_sun. Kepler-538b is a 2.215 +/- (0.040, 0.034) R_earth sub-Neptune with a period of P = 81.73778 +/- 0.00013 d. It is the only known planet in the system. We collected radial velocity (RV) observations with HIRES on Keck I and HARPS-N on the TNG. We characterized stellar activity by a Gaussian process with a quasi-periodic kernel applied to our RV and cross correlation function full width at half maximum (FWHM) observations. By simultaneously modeling Kepler photometry, RV, and FWHM observations, we found a semi-amplitude of K = 1.68 +/- (0.39, 0.38) m s^-1 and a planet mass of M_p = 10.6 +/- (2.5, 2.4) M_earth. Kepler-538b is the smallest planet beyond P = 50 d with an RV mass measurement. The planet likely consists of a significant fraction of ices (dominated by water ice), in addition to rocks/metals, and a small amount of gas. Sophisticated modeling techniques such as those used in this paper, combined with future spectrographs with ultra high-precision and stability will be vital for yielding more mass measurements in this poorly understood exoplanet regime. This in turn will improve our understanding of the relationship between planet composition and insolation flux and how the rocky to gaseous transition depends on planetary equilibrium temperature.

We present results from a data challenge posed to the radial velocity (RV) community: namely, to quantify the Bayesian "evidence" for n={0,1,2,3} planets in a set of synthetically generated RV datasets containing a range of planet signals. Participating teams were provided the same likelihood function and set of priors to use in their analysis. They applied a variety of methods to estimate Z, the marginal likelihood for each n-planet model, including cross-validation, the Laplace approximation, importance sampling, and nested sampling. We found the dispersion in Z across different methods grew with increasing n-planet models: ~3 for 0-planets, ~10 for 1-planet, ~100-1000 for 2-planets, and >10,000 for 3-planets. Most internal estimates of uncertainty in Z for individual methods significantly underestimated the observed dispersion across all methods. Methods that adopted a Monte Carlo approach by comparing estimates from multiple runs yielded plausible uncertainties. Finally, two classes of numerical algorithms (those based on importance and nested samplers) arrived at similar conclusions regarding the ratio of Zs for n and (n+1)-planet models. One analytic method (the Laplace approximation) demonstrated comparable performance. We express both optimism and caution: we demonstrate that it is practical to perform rigorous Bayesian model comparison for <=3-planet models, yet robust planet discoveries require researchers to better understand the uncertainty in Z and its connections to model selection.

We present a detailed analysis of HARPS-N radial velocity observations of K2-100, a young and active star in the Praesepe cluster, which hosts a transiting planet with a period of 1.7 days. We model the activity-induced radial velocity variations of the host star with a multi-dimensional Gaussian Process framework and detect a planetary signal of $10.6 \pm 3.0 {\rm m\,s^{-1}}$, which matches the transit ephemeris, and translates to a planet mass of $21.8 \pm 6.2 M_\oplus$. We perform a suite of validation tests to confirm that our detected signal is genuine. This is the first mass measurement for a transiting planet in a young open cluster. The relatively low density of the planet, $2.04^{+0.66}_{-0.61} {\rm g\,cm^{-3}}$, implies that K2-100b retains a significant volatile envelope. We estimate that the planet is losing its atmosphere at a rate of $10^{11}-10^{12}\,{\rm g\,s^{-1}}$ due to the high level of radiation it receives from its host star.

This paper reports on the validation and mass measurement of K2-263b, a sub-Neptune orbiting a quiet G9V star. Using K2 data from campaigns C5 and C16, we find this planet to have a period of $50.818947\pm 0.000094$ days and a radius of $2.41\pm0.12$ R$_{\oplus}$. We followed this system with HARPS-N to obtain 67 precise radial velocities. A combined fit of the transit and radial velocity data reveals that K2-263b has a mass of $14.8\pm3.1$ M$_{\oplus}$. Its bulk density ($5.7_{-1.4}^{+1.6}$ g cm$^{-3}$) implies that this planet has a significant envelope of water or other volatiles around a rocky core. EPIC211682544b likely formed in a similar way as the cores of the four giant planets in our own Solar System, but for some reason, did not accrete much gas. The planetary mass was confirmed by an independent Gaussian process-based fit to both the radial velocities and the spectroscopic activity indicators. K2-263b belongs to only a handful of confirmed K2 exoplanets with periods longer than 40 days. It is among the longest periods for a small planet with a precisely determined mass using radial velocities.

Ultra-short period (USP) planets are a class of low mass planets with periods shorter than one day. Their origin is still unknown, with photo-evaporation of mini-Neptunes and in-situ formation being the most credited hypotheses. Formation scenarios differ radically in the predicted composition of USP planets, it is therefore extremely important to increase the still limited sample of USP planets with precise and accurate mass and density measurements. We report here the characterization of an USP planet with a period of 0.28 days around K2-141 (EPIC 246393474), and the validation of an outer planet with a period of 7.7 days in a grazing transit configuration. We derived the radii of the planets from the K2 light curve and used high-precision radial velocities gathered with the HARPS-N spectrograph for mass measurements. For K2-141b we thus inferred a radius of $1.51\pm0.05~R_\oplus$ and a mass of $5.08\pm0.41~M_\oplus$, consistent with a rocky composition and lack of a thick atmosphere. K2-141c is likely a Neptune-like planet, although due to the grazing transits and the non-detection in the RV dataset, we were not able to put a strong constraint on its density. We also report the detection of secondary eclipses and phase curve variations for K2-141b. The phase variation can be modeled either by a planet with a geometric albedo of $0.30 \pm 0.06$ in the Kepler bandpass, or by thermal emission from the surface of the planet at $\sim$3000K. Only follow-up observations at longer wavelengths will allow us to distinguish between these two scenarios.

Initial RV characterisation of the enigmatic planet Kepler-10c suggested a mass of $\sim17$ M$_\oplus$, which was remarkably high for a planet with radius $2.32$ R$_\oplus$; further observations and subsequent analysis hinted at a (possibly much) lower mass, but masses derived using RVs from two different spectrographs (HARPS-N and HIRES) were incompatible at a $3\sigma$-level. We demonstrate here how such mass discrepancies may readily arise from sub-optimal sampling and/or neglecting to model even a single coherent signal (stellar, planetary, or otherwise) that may be present in RVs. We then present a plausible resolution of the mass discrepancy, and ultimately characterise Kepler-10c as having mass $7.37_{-1.19}^{+1.32}$ M$_\oplus$, and mean density $3.14^{+0.63}_{-0.55}$ g cm$^{-3}$.

Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasi-periodic --- spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalising over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine-fitting periodogram and an autocorrelation function method. We also demonstrate that it works well on real data, by inferring rotation periods for 275 Kepler stars with previously measured periods. We provide a table of rotation periods for these 1132 Kepler objects of interest and their posterior probability density function samples. Because this method delivers posterior probability density functions, it will enable hierarchical studies involving stellar rotation, particularly those involving population modelling, such as inferring stellar ages, obliquities in exoplanet systems, or characterising star-planet interactions. The code used to implement this method is available online.

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.

Context. The presence of a small-mass planet (M$_p<$0.1\,M$_{Jup}$) seems, to date, not to depend on metallicity, however, theoretical simulations have shown that stars with subsolar metallicities may be favoured for harbouring smaller planets. A large, dedicated survey of metal-poor stars with the HARPS spectrograph has thus been carried out to search for Neptunes and super-Earths. Aims. In this paper, we present the analysis of \object{HD175607}, an old G6 star with metallicity [Fe/H] = -0.62. We gathered 119 radial velocity measurements in 110 nights over a time span of more than nine years. Methods. The radial velocities were analysed using Lomb-Scargle periodograms, a genetic algorithm, a Markov chain Monte Carlo analysis, and a Gaussian processes analysis. The spectra were also used to derive stellar properties. Several activity indicators were analysed to study the effect of stellar activity on the radial velocities. Results. We find evidence for the presence of a small Neptune-mass planet (M$_{p}\sin i = 8.98\pm1.10$\,M$_{\oplus}$) orbiting this star with an orbital period $P = 29.01\pm0.02$\, days in a slightly eccentric orbit ($e=0.11\pm0.08$). The period of this Neptune is close to the estimated rotational period of the star. However, from a detailed analysis of the radial velocities together with the stellar activity, we conclude that the best explanation of the signal is indeed the presence of a planetary companion rather than stellar related. An additional longer period signal ($P\sim 1400$\,d) is present in the data, for which more measurements are needed to constrain its nature and its properties. Conclusions. HD\,175607 is the most metal-poor FGK dwarf with a detected low-mass planet amongst the currently known planet hosts. This discovery may thus have important consequences for planet formation and evolution theories.

We re-analyse the publicly available radial velocity (RV) measurements for Alpha Cen B, a star hosting an Earth-mass planet candidate, Alpha Cen Bb, with 3.24 day orbital period. We demonstrate that the 3.24 d signal observed in the Alpha Cen B data almost certainly arises from the window function (time sampling) of the original data. We show that when stellar activity signals are removed from the RV variations, other significant peaks in the power spectrum of the window function are coincidentally suppressed, leaving behind a spurious yet apparently-significant 'ghost' of a signal that was present in the window function's power spectrum to begin with. Even when fitting synthetic data with time sampling identical to the original data, but devoid of any genuine periodicities close to that of the planet candidate, the original model used to infer the presence of Alpha Cen Bb leads to identical conclusions: viz., the 3$\sigma$ detection of a half-a-metre-per-second signal with 3.236 day period. Our analysis underscores the difficulty of detecting weak planetary signals in RV data, and the importance of understanding in detail how every component of an RV data set, including its time sampling, influences final statistical inference.

To date, the radial velocity (RV) method has been one of the most productive techniques for detecting and confirming extrasolar planetary candidates. Unfortunately, stellar activity can induce RV variations which can drown out or even mimic planetary signals - and it is notoriously difficult to model and thus mitigate the effects of these activity-induced nuisance signals. This is expected to be a major obstacle to using next-generation spectrographs to detect lower mass planets, planets with longer periods, and planets around more active stars. Enter Gaussian processes (GPs) which, we note, have a number of attractive features that make them very well suited to disentangling stellar activity signals from planetary signals. We present here a GP framework we developed to model RV time series jointly with ancillary activity indicators (e.g. bisector velocity spans, line widths, chromospheric activity indices), allowing the activity component of RV time series to be constrained and disentangled from e.g. planetary components. We discuss the mathematical details of our GP framework, and present results illustrating its encouraging performance on both synthetic and real RV datasets, including the publicly-available Alpha Centauri B dataset.

The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4 GHz continuum emission traced by the Very Large Array (VLA) and the Gemini/T-ReCS 8.7 micron emission, as well as the Pa-alpha line emission, show remarkable morphological similarities within the star-forming ring, suggesting that the underlying emission mechanisms are tightly related. We used an HST/NICMOS Pa-alpha map of similar resolution to our radio maps to disentangle the thermal free-free and non-thermal synchrotron radio emission, from which we obtained the intrinsic synchrotron power-law for each individual region within the central kpc of NGC1614. The radio ring surrounds a relatively faint, steep-spectrum source at the very center of the galaxy, suggesting that the central source is not powered by an AGN, but rather by a compact (r < 90 pc) starburst. Chandra X-ray data also show that the central kpc region is dominated by starburst activity, without requiring the existence of an AGN. We also used publicly available infrared data to model-fit the spectral energy distribution of both the starburst ring and a putative AGN in NGC1614. In summary, we conclude that there is no need to invoke an AGN to explain the observed bolometric properties of the galaxy.

A new algorithm developed to perform autonomous fitting of gravitational microlensing lightcurves is presented. The new algorithm is conceptually simple, versatile and robust, and parallelises trivially; it combines features of extant evolutionary algorithms with some novel ones, and fares well on the problem of fitting binary-lens microlensing lightcurves, as well as on a number of other difficult optimisation problems. Success rates in excess of 90% are achieved when fitting synthetic though noisy binary-lens lightcurves, allowing no more than 20 minutes per fit on a desktop computer; this success rate is shown to compare very favourably with that of both a conventional (iterated simplex) algorithm, and a more state-of-the-art, artificial neural network-based approach. As such, this work provides proof of concept for the use of an evolutionary algorithm as the basis for real-time, autonomous modelling of microlensing events. Further work is required to investigate how the algorithm will fare when faced with more complex and realistic microlensing modelling problems; it is, however, argued here that the use of parallel computing platforms, such as inexpensive graphics processing units, should allow fitting times to be constrained to under an hour, even when dealing with complicated microlensing models. In any event, it is hoped that this work might stimulate some interest in evolutionary algorithms, and that the algorithm described here might prove useful for solving microlensing and/or more general model-fitting problems.

We study star formation (SF) in very active environments, in luminous IR galaxies, which are often interacting. A variety of phenomena are detected, such as central starbursts, circumnuclear SF, obscured SNe tracing the history of recent SF, massive super star clusters, and sites of strong off-nuclear SF. All of these can be ultimately used to define the sequence of triggering and propagation of star-formation and interplay with nuclear activity in the lives of gas rich galaxy interactions and mergers. In this paper we present analysis of high-spatial resolution integral field spectroscopy of central regions of two interacting LIRGs. We detect a nuclear 3.3 um PAH ring around the core of NGC 1614 with thermal-IR IFU observations. The ring's characteristics and relation to the strong star-forming ring detected in recombination lines are presented, as well as a scenario of an outward expanding starburst likely initiated with a (minor) companion detected within a tidal feature. We then present NIR IFU observations of IRAS 19115-2124, aka the Bird, which is an intriguing triple encounter. The third component is a minor one, but, nevertheless, is the source of 3/4 of the SFR of the whole system. Gas inflows and outflows are detected at the locations of the nuclei. Finally, we briefly report on our on-going NIR adaptive optics imaging survey of several dozen LIRGs. We have detected highly obscured core-collapse SNe in the central kpc, and discuss the statistics of "missing SNe" due to dust extinction. We are also determining the characteristics of hundreds of super star clusters in and around the core regions of LIRGs, as a function of host-galaxy properties.

Genetic algorithms (GAs) emulate the process of biological evolution, in a computational setting, in order to generate good solutions to difficult search and optimisation problems. GA-based optimisers tend to be extremely robust and versatile compared to most traditional techniques used to solve optimisation problems. This review paper provides a very brief introduction to GAs and outlines their utility in astronomy and astrophysics.

It is important to understand the interplay between nuclear star-formation and nuclear activity when studying the evolution of gas rich galaxy mergers. We present here new spatially resolved L-band integral field unit observations of the inner kpc of the luminous IR galaxy NGC 1614. A broad ring of 3.3 um PAH emission is found at a distance of approximately 200 pc from the core. This ring overlaps with a previously established star-forming ring detected with Pa_alpha and radio continuum observations, but peaks outside it, especially if determined using the PAH equivalent width. Using the characteristics of the PAH emission and the ionised gas emission we argue that NGC 1614 features an outward propagating ring of star formation, where the equivalent width of the PAH emission localises the regions where the current star formation is just expanding into the molecular gas outward of the nucleus. The core itself shows a highly luminous, slightly resolved (at ~80 pc) L-band continuum source. We find no evidence of AGN activity and rule out the presence of an obscured AGN using L-band diagnostics. Furthermore, we detect the likely companion galaxy from archival HST/ACS imaging. The star formation and dynamical characteristics of the system are consistent with a relatively major merger just after its second passage. An outstanding question is how a gas-rich advanced merger such as this one, with strong LIRG level nuclear starburst and major merger-like tidal features, has not yet developed an active nucleus.