The nearest exoplanets to the Sun are our best possibilities for detailed characterization. We report the discovery of a compact multi-planet system of super-Earths orbiting the nearby red dwarf GJ 887, using radial velocity measurements. The planets have orbital periods of 9.3 and 21.8~days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8 day planet is approx 350 K; which is interior, but close to the inner edge, of the liquid-water habitable zone. We also detect a further unconfirmed signal with a period of 50 days which could correspond to a third super-Earth in a more temperate orbit. GJ 887 is an unusually magnetically quiet red dwarf with a photometric variability below 500 parts-per-million, making its planets amenable to phase-resolved photometric characterization.

AU Mic~b is a Neptune size planet on a 8.47-day orbit around the nearest pre-main sequence ($\sim$20 Myr) star to the Sun, the bright (V=8.81) M dwarf AU Mic. The planet was preliminary detected in Doppler radial velocity time series and recently confirmed to be transiting with data from the TESS mission. AU Mic~b is likely to be cooling and contracting and might be accompanied by a second, more massive planet, in an outer orbit. Here, we present the observations of the transit of AU Mic~b using ESPRESSO on the VLT. We obtained a high-resolution time series of spectra to measure the Rossiter-McLaughlin effect and constrain the spin-orbit alignment of the star and planet, and simultaneously attempt to retrieve the planet's atmospheric transmission spectrum. These observations allow us to study for the first time the early phases of the dynamical evolution of young systems. We apply different methodologies to derive the spin-orbit angle of AU Mic~b, and all of them retrieve values consistent with the planet being aligned with the rotation plane of the star. We determine a conservative spin-orbit angle $\lambda$ value of $-2.96^{+10.44}_{-10.30}$, indicative that the formation and migration of the planets of the AU Mic system occurred within the disk. Unfortunately, and despite the large SNR of our measurements, the degree of stellar activity prevented us from detecting any features from the planetary atmosphere. In fact, our results suggest that transmission spectroscopy for recently formed planets around active young stars is going to remain very challenging, if at all possible, for the near future.

The ability to make independent detections of the signatures of exoplanets with complementary telescopes and instruments brings a new potential for robust identification of exoplanets and precision characterization. We introduce PEXO, a package for Precise EXOplanetology to facilitate the efficient modeling of timing, astrometry, and radial velocity data, which will benefit not only exoplanet science but also various astrophysical studies in general. PEXO is general enough to account for binary motion and stellar reflex motions induced by planetary companions and is precise enough to treat various relativistic effects both in the solar system and in the target system. We also model the post-Newtonian barycentric motion for future tests of general relativity in extrasolar systems. We benchmark PEXO with the pulsar timing package TEMPO2 and find that PEXO produces numerically similar results with timing precision of about 1 ns, space-based astrometry to a precision of 1 $\mu$as, and radial velocity of 1 $\mu$m/s and improves on TEMPO2 for decade-long timing data of nearby targets, due to its consideration of third-order terms of Roemer delay. PEXO is able to avoid the bias introduced by decoupling the target system and the solar system and to account for the atmospheric effects which set a practical limit for ground-based radial velocities close to 1 cm/s. Considering the various caveats in barycentric correction and ancillary data required to realize cm/s modeling, we recommend the preservation of original observational data. The PEXO modeling package is available at GitHub (https://github.com/phillippro/pexo).

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

We report individual dynamical masses for the brown dwarfs Epsilon Indi B and C, which have spectral types of T1.5 and T6, respectively, measured from astrometric orbit mapping. Our measurements are based on a joint analysis of astrometric data from the Carnegie Astrometric Planet Search and the Cerro Tololo Inter-American Observatory Parallax Investigation as well as archival high resolution imaging, and use a Markov Chain Monte Carlo method. We find dynamical masses of 75.0+-0.82 Mjup for the T1.5 B component and 70.1+-0.68 Mjup for the T6 C component. These masses are surprisingly high for substellar objects and challenge our understanding of substellar structure and evolution. We discuss several evolutionary scenarios proposed in the literature and find that while none of them can provide conclusive explanations for the high substellar masses, evolutionary models incorporating lower atmospheric opacities come closer to approximating our results. We discuss the details of our astrometric model, its algorithm implementation, and how we determine parameter values via Markov Chain Monte Carlo Bayesian inference.

Exoplanet Proxima b will be an important laboratory for the search for extraterrestrial life for the decades ahead. Here we discuss the prospects of detecting carbon dioxide at 15 um using a spectral filtering technique with the Medium Resolution Spectrograph (MRS) mode of the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST). At superior conjunction, the planet is expected to show a contrast of up to 100 ppm with respect to the star. At a spectral resolving power of R=1790-2640, about 100 spectral CO2 features are visible within the 13.2-15.8 um (3B) band, which can be combined to boost the planet atmospheric signal by a factor 3-4, depending on the atmospheric temperature structure and CO2 abundance. If atmospheric conditions are favorable (assuming an Earth-like atmosphere), with this new application to the cross-correlation technique carbon dioxide can be detected within a few days of JWST observations. However, this can only be achieved if both the instrumental spectral response and the stellar spectrum can be determined to a relative precision of <= 1e-4 between adjacent spectral channels. Absolute flux calibration is not required, and the method is insensitive to the strong broadband variability of the host star. Precise calibration of the spectral features of the host star may only be attainable by obtaining deep observations of the system during inferior conjunction that serve as a reference. The high-pass filter spectroscopic technique with the MIRI MRS can be tested on warm Jupiters, Neptunes, and super-Earths with significantly higher planet/star contrast ratios than the Proxima system.

High precision measurements of stellar spectroscopic line profiles and their changes over time contain very valuable information about the physics of the stellar photosphere (stellar activity) and can be used to characterize extrasolar planets via the Rossiter-McLaughlin effect or from reflected light from the planet. In this paper we present a new method for measuring small changes in the mean line profile of a spectrum by performing what we call differential Least Squares Deconvolution (dLSD). The method consists in finding the convolution function (or kernel) required to transform a high signal-to-noise ratio template of the star into each observed spectrum. Compared to similar techniques, the method presented here does not require any assumptions on the template spectrum (eg. no line-list or cross-correlation mask required). We show that our implementation of dLSD is able to perform -at least- as good as other techniques by applying it to star-planet obliquity measurements of exoplanet HD183799 during its transit. Among other things, the method should enable model independent detection of light reflected by an exoplanet.

Transit photometry of the M8V dwarf star TRAPPIST-1 (2MASS J23062928-0502285) has revealed the presence of at least seven planets with masses and radii similar to that of Earth orbiting at distances that might allow liquid water to be present on their surfaces. We have been following TRAPPIST-1 since 2011 with the CAPSCam astrometric camera on the 2.5-m du Pont telescope at the Las Campanas Observatory in Chile. In 2016 we noted that TRAPPIST-1 lies slightly farther away than previously thought, at 12.49 pc, rather than 12.1 pc. Here we examine fifteen epochs of CAPSCam observations of TRAPPIST-1, spanning the five years from 2011 to 2016, and obtain a revised trigonometric distance of $12.56 \pm 0.12$ pc. The astrometric data analysis pipeline shows no evidence for a long-period astrometric wobble of TRAPPIST-1. After proper motion and parallax are removed, residuals at the level of $\pm 1.3$ millarcsec (mas) remain. The amplitude of these residuals constrains the masses of any long-period gas giant planets in the TRAPPIST-1 system: no planet more massive than $\sim 4.6 M_{Jup}$ orbits with a 1 yr period, and no planet more massive than $\sim 1.6 M_{Jup}$ orbits with a 5 yr period. Further refinement of the CAPSCam data analysis pipeline, combined with continued CAPSCam observations, should either detect any long-period planets, or put an even tighter constraint on these mass upper limits.

The removal of noise typically correlated in time and wavelength is one of the main challenges for using the radial velocity method to detect Earth analogues. We analyze radial velocity data of tau Ceti and find robust evidence for wavelength dependent noise. We find this noise can be modeled by a combination of moving average models and "differential radial velocities". We apply this noise model to various radial velocity data sets for tau Ceti, and find four periodic signals at 20.0, 49.3, 160 and 642 d which we interpret as planets. We identify two new signals with orbital periods of 20.0 and 49.3 d while the other two previously suspected signals around 160 and 600 d are quantified to a higher precision. The 20.0 d candidate is independently detected in KECK data. All planets detected in this work have minimum masses less than 4$M_\oplus$ with the two long period ones located around the inner and outer edges of the habitable zone, respectively. We find that the instrumental noise gives rise to a precision limit of the HARPS around 0.2 m/s. We also find correlation between the HARPS data and the central moments of the spectral line profile at around 0.5 m/s level, although these central moments may contain both noise and signals. The signals detected in this work have semi-amplitudes as low as 0.3 m/s, demonstrating the ability of the radial velocity technique to detect relatively weak signals.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We present a study of four different formation scenarios that may be able to explain the origin of the recently announced planet (`Proxima b') orbiting Proxima Centauri. The aim is to examine how the formation scenarios differ in their predictions for the multiplicity of the Proxima system, the water/volatile content of Proxima b and its orbital eccentricity, so that these can be tested by future observations. A scenario of in situ formation via giant impacts from a locally enhanced disc of planetary embryos/planetesimals predicts that Proxima b will be in a multiplanet system with a measurably finite eccentricity. Assuming that the local solid enhancement needed to form a Proxima b analogue arises because of the inwards drift of solids in the form of small planetesimals/boulders, this scenario also results in Proxima b analogues being only moderately endowed with water/volatiles. A scenario in which multiple embryos form, migrate and mutually collide within a gas disc results in Proxima b being a member of a multiple system, possibly displaying mean motion resonances, but where the constituent members are Ocean planets due to accretion occurring mainly outside of the snowline. A scenario in which a single accreting embryo forms outside the snowline, and migrates inwards while accreting planetesimals/pebbles results in Proxima b being an isolated Ocean planet on a circular orbit. A scenario in which Proxima b formed via pebble accretion interior to the snowline produces a dry planet on a circular orbit. Future observations that characterise the physical and orbital properties of Proxima b, and any additional planets in the system, will provide valuable insights into the formation history of this neighbouring planetary system.

We report the discovery of a second long-period giant planet orbiting HD 30177, a star previously known to host a massive Jupiter analog (HD 30177b: a=3.8$\pm$0.1 au, m sin $i=9.7\pm$0.5 Mjup). HD 30177c can be regarded as a massive Saturn analog in this system, with a=9.9$\pm$1.0 au and m sin $i=7.6\pm$3.1 Mjup. The formal best fit solution slightly favours a closer-in planet at $a\sim$7 au, but detailed n-body dynamical simulations show that configuration to be unstable. A shallow local minimum of longer-period, lower-eccentricity solutions was found to be dynamically stable, and hence we adopt the longer period in this work. The proposed $\sim$32 year orbit remains incomplete; further monitoring of this and other stars is necessary to reveal the population of distant gas giant planets with orbital separations $a\sim$10 au, analogous to that of Saturn.

Accounting for stellar activity is a crucial component of the search for ever-smaller planets orbiting stars of all spectral types. We use Doppler imaging methods to demonstrate that starspot induced radial velocity variability can be effectively reduced for moderately rotating, fully convective stars. Using starspot distributions extrapolated from sunspot observations, we adopt typical M dwarf starspot distributions with low contrast spots to synthesise line profile distortions. The distortions are recovered using maximum entropy regularised fitting and the corresponding stellar radial velocities are measured. The procedure is demonstrated for a late-M star harbouring an orbiting planet in the habitable zone. The technique is effective for stars with vsini = 1-10 km/s, reducing the stellar noise contribution by factors of nearly an order of magnitude. With a carefully chosen observing strategy, the technique can be used to determine the stellar rotation period and is robust to uncertainties such as unknown stellar inclination. While demonstrated for late-type M stars, the procedure is applicable to all spectral types.

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

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.

HD 11112 is an old, Sun-like star that has a long-term radial velocity (RV) trend indicative of a massive companion on a wide orbit. Here we present direct images of the source responsible for the trend using the Magellan Adaptive Optics system. We detect the object (HD 11112B) at a separation of 2\fasec 2 (100 AU) at multiple wavelengths spanning 0.6-4 \microns ~and show that it is most likely a gravitationally-bound cool white dwarf. Modeling its spectral energy distribution (SED) suggests that its mass is 0.9-1.1 \msun, which corresponds to very high-eccentricity, near edge-on orbits from Markov chain Monte Carlo analysis of the RV and imaging data together. The total age of the white dwarf is $>2\sigma$ discrepant with that of the primary star under most assumptions. The problem can be resolved if the white dwarf progenitor was initially a double white dwarf binary that then merged into the observed high-mass white dwarf. HD 11112B is a unique and intriguing benchmark object that can be used to calibrate atmospheric and evolutionary models of cool white dwarfs and should thus continue to be monitored by RV and direct imaging over the coming years.

We present the results of two 2.3 micron near-infrared radial velocity surveys to detect exoplanets around 36 nearby and young M dwarfs. We use the CSHELL spectrograph (R ~46,000) at the NASA InfraRed Telescope Facility, combined with an isotopic methane absorption gas cell for common optical path relative wavelength calibration. We have developed a sophisticated RV forward modeling code that accounts for fringing and other instrumental artifacts present in the spectra. With a spectral grasp of only 5 nm, we are able to reach long-term radial velocity dispersions of ~20-30 m/s on our survey targets.

We present the results of a precise near-infrared (NIR) radial velocity (RV) survey of 32 low-mass stars with spectral types K2-M4 using CSHELL at the NASA IRTF in the $K$-band with an isotopologue methane gas cell to achieve wavelength calibration and a novel iterative RV extraction method. We surveyed 14 members of young ($\approx$ 25-150 Myr) moving groups, the young field star $\varepsilon$ Eridani as well as 18 nearby ($<$ 25 pc) low-mass stars and achieved typical single-measurement precisions of 8-15 m s$^{-1}$ with a long-term stability of 15-50 m s$^{-1}$. We obtain the best NIR RV constraints to date on 27 targets in our sample, 19 of which were never followed by high-precision RV surveys. Our results indicate that very active stars can display long-term RV variations as low as $\sim$ 25-50 m s$^{-1}$ at $\approx$ 2.3125 $\mu$m, thus constraining the effect of jitter at these wavelengths. We provide the first multi-wavelength confirmation of GJ 876 bc and independently retrieve orbital parameters consistent with previous studies. We recovered RV variability for HD 160934 AB and GJ 725 AB that are consistent with their known binary orbits, and nine other targets are candidate RV variables with a statistical significance of 3-5$\sigma$. Our method combined with the new iSHELL spectrograph will yield long-term RV precisions of $\lesssim$ 5 m s$^{-1}$ in the NIR, which will allow the detection of Super-Earths near the habitable zone of mid-M dwarfs.

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.

We present high-contrast Magellan adaptive optics (MagAO) images of HD 7449, a Sun-like star with one planet and a long-term radial velocity (RV) trend. We unambiguously detect the source of the long-term trend from 0.6-2.15 \microns ~at a separation of \about 0\fasec 54. We use the object's colors and spectral energy distribution to show that it is most likely an M4-M5 dwarf (mass \about 0.1-0.2 \msun) at the same distance as the primary and is therefore likely bound. We also present new RVs measured with the Magellan/MIKE and PFS spectrometers and compile these with archival data from CORALIE and HARPS. We use a new Markov chain Monte Carlo procedure to constrain both the mass ($> 0.17$ \msun ~at 99$\%$ confidence) and semimajor axis (\about 18 AU) of the M dwarf companion (HD 7449B). We also refine the parameters of the known massive planet (HD 7449Ab), finding that its minimum mass is $1.09^{+0.52}_{-0.19}$ \mj, its semimajor axis is $2.33^{+0.01}_{-0.02}$ AU, and its eccentricity is $0.8^{+0.08}_{-0.06}$. We use N-body simulations to constrain the eccentricity of HD 7449B to $\lesssim$ 0.5. The M dwarf may be inducing Kozai oscillations on the planet, explaining its high eccentricity. If this is the case and its orbit was initially circular, the mass of the planet would need to be $\lesssim$ 1.5 \mj. This demonstrates that strong constraints on known planets can be made using direct observations of otherwise undetectable long-period companions.

CRIRES, the ESO high resolution infrared spectrometer, is a unique instrument which allows astronomers to access a parameter space which up to now was largely uncharted. In its current setup, it consists of a single-order spectrograph providing long-slit, single-order spectroscopy with resolving power up to R=100,000 over a quite narrow spectral range. This has resulted in sub-optimal efficiency and use of telescope time for all the scientific programs requiring broad spectral coverage of compact objects (e.g. chemical abundances of stars and intergalactic medium, search and characterization of extra-solar planets). To overcome these limitations, a consortium was set-up for upgrading CRIRES to a cross-dispersed spectrometer, called CRIRES+. This paper presents the updated optical design of the crossdispersion module for CRIRES+. This new module can be mounted in place of the current pre-disperser unit. The new system yields a factor of >10 increase in simultaneous spectral coverage and maintains a quite long slit (10"), ideal for observations of extended sources and for precise sky-background subtraction.

We report the detection of GJ 832c, a super-Earth orbiting near the inner edge of the habitable zone of GJ 832, an M dwarf previously known to host a Jupiter analog in a nearly-circular 9.4-year orbit. The combination of precise radial-velocity measurements from three telescopes reveals the presence of a planet with a period of 35.68+/-0.03 days and minimum mass (m sin i) of 5.4+/-1.0 Earth masses. GJ 832c moves on a low-eccentricity orbit (e=0.18+/-0.13) towards the inner edge of the habitable zone. However, given the large mass of the planet, it seems likely that it would possess a massive atmosphere, which may well render the planet inhospitable. Indeed, it is perhaps more likely that GJ 832c is a "super-Venus," featuring significant greenhouse forcing. With an outer giant planet and an interior, potentially rocky planet, the GJ 832 planetary system can be thought of as a miniature version of our own Solar system.

Variations related to stellar activity and correlated noise can prevent the detections of low-amplitude signals in radial velocity data if not accounted for. This can be seen as the greatest obstacle in detecting Earth-like planets orbiting nearby stars with Doppler spectroscopy regardless of developments in instrumentation and rapidly accumulating amounts of data. We use a statistical model that is not sensitive to aperiodic and/or quasiperiodic variability of stellar origin. We demonstrate the performance of our model by re-analysing the radial velocities of the moderately active star CoRoT-7 ($\log R_{\rm HK} = -4.61$) with a transiting planet whose Doppler signal has proven rather difficult to detect. We find that the signal of the transiting planet can be robustly detected together with signals of two other planet candidates. Our results suggest that rotation periods of moderately active stars can be filtered out of the radial velocity noise, which enables the detections of low-mass planets orbiting such stars.

We present radial velocity measurements of a sample of M5V-M9V stars from our Red-Optical Planet Survey, ROPS, operating at 0.65-1.025 micron. Radial velocities for 15 stars, with r.m.s. precision down to 2.5 m/s over a week long time scale are achieved using Thorium-Argon reference spectra. We are sensitive to planets with m_psin(i) >= 1.5 MEarth (3 MEarth at 2-sigma) in the classical habitable zone and our observations currently rule out planets with m_psin(i) > 0.5 MJup at 0.03 AU for all our targets. A total of 9 of the 15 targets exhibit r.m.s. < 16 m/s, which enables us to rule out the presence of planets with m_psin(i) > 10 MEarth in 0.03 AU orbits. Since the mean rotation velocity is of order 8 km/s for an M6V star and 15 km/s by M9V, we avoid observing only slow rotators that would introduce a bias towards low axial inclination i << 90 deg systems, which are unfavourable for planet detection. Our targets with the highest vsini values exhibit radial velocities significantly above the photon-noise limited precision, even after accounting for vsini. We therefore monitored stellar activity via chromospheric emission from the Halpha and Ca II infrared triplet lines. A clear trend of log_10(Halpha/Lbol) with radial velocity r.m.s. is seen, implying that significant starspot activity is responsible for the observed radial velocity precision floor. The implication that most late M dwarfs are significantly spotted, and hence exhibit time varying line distortions, indicates that observations to detect orbiting planets need strategies to reliably mitigate against the effects of activity induced radial velocity variations.

We have built and commissioned a prototype agitated non-circular core fiber scrambler for precision spectroscopic radial velocity measurements in the near-infrared H band. We have collected the first on-sky performance and modal noise tests of these novel fibers in the near-infrared at H and K bands using the CSHELL spectrograph at the NASA InfraRed Telescope Facility (IRTF). We discuss the design behind our novel reverse injection of a red laser for co-alignment of star-light with the fiber tip via a corner cube and visible camera. We summarize the practical details involved in the construction of the fiber scrambler, and the mechanical agitation of the fiber at the telescope. We present radial velocity measurements of a bright standard star taken with and without the fiber scrambler to quantify the relative improvement in the obtainable blaze function stability, the line spread function stability, and the resulting radial velocity precision. We assess the feasibility of applying this illumination stabilization technique to the next generation of near-infrared spectrographs such as iSHELL on IRTF and an upgraded NIRSPEC at Keck. Our results may also be applied in the visible for smaller core diameter fibers where fiber modal noise is a significant factor, such as behind an adaptive optics system or on a small < 1 meter class telescope such as is being pursued by the MINERVA and LCOGT collaborations.

We have built and commissioned gas absorption cells for precision spectroscopic radial velocity measurements in the near-infrared in the H and K bands. We describe the construction and installation of three such cells filled with 13CH4, 12CH3D, and 14NH3 for the CSHELL spectrograph at the NASA Infrared Telescope Facility (IRTF). We have obtained their high-resolution laboratory Fourier Transform spectra, which can have other practical uses. We summarize the practical details involved in the construction of the three cells, and the thermal and mechanical control. In all cases, the construction of the cells is very affordable. We are carrying out a pilot survey with the 13CH4 methane gas cell on the CSHELL spectrograph at the IRTF to detect exoplanets around low mass and young stars. We discuss the current status of our survey, with the aim of photon-noise limited radial velocity precision. For adequately bright targets, we are able to probe a noise floor of 7 m/s with the gas cell with CSHELL at cassegrain focus. Our results demonstrate the feasibility of using a gas cell on the next generation of near-infrared spectrographs such as iSHELL on IRTF, iGRINS, and an upgraded NIRSPEC at Keck.

Exoplanet Doppler surveys are currently the most efficient means to detect low-mass companions to nearby stars. Among these stars, the light M dwarfs provide the highest sensitivity to detect low-mass exoplanet candidates. Evidence is accumulating that a substantial fraction of these low-mass planets are found in high-multiplicity planetary systems. GJ 163 is a nearby inactive M dwarf with abundant public observations obtained using the HARPS spectrograph. We obtain and analyse radial velocities from the HARPS public spectra of GJ 163 and investigate the presence of a planetary companions orbiting it. The number of planet candidates detected might depend on some prior assumptions. Since the impact of prior choice has not been investigated throughly previously, we study the effects of different prior densities on the detectability of planet candidates around GJ 163. We use Bayesian tools, i.e. posterior samplings and model comparisons, when analysing the GJ 163 velocities. We consider models accounting for the possible correlations of subsequent measurements. We also search for activity-related counterparts of the signals we observe and test the dynamical stability of the planetary systems corresponding to our solutions using direct numerical integrations of the orbits. We find that there are at least three planet candidates orbiting GJ 163. The existence of a fourth planet is supported by the data but the evidence in favor of the corresponding model is not yet conclusive. The second innermost planet candidate in the system with an orbital period of 25.6 days and a minimum mass of 8.7 Me is inside the liquid-water habitable zone of the star.

Stellar activity signatures such as spots and plage can significantly limit the search for extrasolar planets. Current models of activity-induced radial velocity (RV) signals focused on the impact of temperature contrast in spots predicting the signal to diminish toward longer wavelengths. On the other hand, the relative importance of the Zeeman effect on RV measurements should grow with wavelength, because the Zeeman displacement itself grows with \lambda, and because a magnetic and cool spot contributes more to the total flux at longer wavelengths. We model the impact of active regions on stellar RV measurements including both temperature contrast in spots and Zeeman line broadening. We calculate stellar line profiles using polarized radiative transfer models including atomic and molecular Zeeman splitting from 0.5 to 2.3\mum. Our results show that the amplitude of the RV signal caused by the Zeeman effect alone can be comparable to that caused by temperature contrast. Furthermore, the RV signal caused by cool and magnetic spots increases with wavelength. We also calculate the RV signal due to variations in average magnetic field strength from one observation to the next, but find it unlikely that this can significantly influence the search for extrasolar planets. As an example, we derive the RV amplitude of the active M dwarf AD Leo as a function of wavelength using data from the HARPS spectrograph. The RV signal does not diminish at longer wavelengths but shows evidence for the opposite behavior. We conclude that the RV signal of active stars does not vanish at longer wavelength but sensitively depends on the combination of spot temperature and magnetic field; in active low-mass stars, it is even likely to grow with wavelength.

The K2.5 dwarf HD 40307 has been reported to host three super-Earths. The system lacks massive planets and is therefore a potential candidate for having additional low-mass planetary companions. We re-derive Doppler measurements from public HARPS spectra of HD 40307 to confirm the significance of the reported signals using independent data analysis methods. We also investigate these measurements for additional low-amplitude signals. We used Bayesian analysis of our radial velocities to estimate the probability densities of different model parameters. We also estimated the relative probabilities of models with differing numbers of Keplerian signals and verified their significance using periodogram analyses. We investigated the relation of the detected signals with the chromospheric emission of the star. As previously reported for other objects, we found that radial velocity signals correlated with the S-index are strongly wavelength dependent. We identify two additional clear signals with periods of 34 and 51 days, both corresponding to planet candidates with minimum masses a few times that of the Earth. An additional sixth candidate is initially found at a period of 320 days. However, this signal correlates strongly with the chromospheric emission from the star and is also strongly wavelength dependent. When analysing the red half of the spectra only, the five putative planetary signals are recovered together with a very significant periodicity at about 200 days. This signal has a similar amplitude as the other new signals reported in the current work and corresponds to a planet candidate with M sin i = 7 Me (HD 40307 g). ...

We have conducted a kinematic study of 165 young M dwarfs with ages of <300 Myr. Our sample is composed of stars and brown dwarfs with spectral types ranging from K7 to L0, detected by ROSAT and with photometric distances of <25 pc assuming the stars are single and on the main-sequence. In order to find stars kinematically linked to known young moving groups (YMGs), we measured radial velocities for the complete sample with Keck and CFHT optical spectroscopy and trigonometric parallaxes for 75 of the M dwarfs with the CAPSCam instrument on the du Pont 2.5-m Telescope. Due to their youthful overluminosity and unresolved binarity, the original photometric distances for our sample underestimated the distances by 70% on average, excluding two extremely young (<3 Myr) objects found to have distances beyond a few hundred parsecs. We searched for kinematic matches to 14 reported YMGs and identified 9 new members of the AB Dor YMG and 2 of the Ursa Majoris group. Additional possible candidates include 6 Castor, 4 Ursa Majoris, 2 AB Dor members, and 1 member each of the Her-Lyr and beta Pic groups. Our sample also contains 27 young low-mass stars and 4 brown dwarfs with ages <150 Myr which are not associated with any known YMG. We identified an additional 15 stars which are kinematic matches to one of the YMGs, but the ages from spectroscopic diagnostics and/or the positions on the sky do not match. These warn against grouping stars together based only on kinematics and that a confluence of evidence is required to claim that a group of stars originated from the same star-forming event.

We present new radial velocity measurements of the ultra-cool dwarf VB 10, which was recently announced to host a giant planet detected with astrometry. The new observations were obtained using optical spectrographs(MIKE/Magellan and ESPaDOnS/CHFT) and cover a 63% of the reported period of 270 days. We apply Least-squares periodograms to identify the most significant signals and evaluate their corresponding False Alarm Probabilities. We show that this method is the proper generalization to astrometric data because (1) it mitigates the coupling of the orbital parameters with the parallax and proper motion, and (2) it permits a direct generalization to include non-linear Keplerian parameters in a combined fit to astrometry and radial velocity data. In fact, our analysis of the astrometry alone uncovers the reported 270 d period and an even stronger signal at 50 days. We estimate the uncertainties in the parameters using a Markov Chain Monte Carlo approach. The nominal precision of the new Doppler measurements is about 150 s$^{-1}$ while their standard deviation is 250 ms$^{-1}$. However, the best fit solutions still have RMS of 200 ms$^{-1}$ indicating that the excess in variability is due to uncontrolled systematic errors rather than the candidate companions detected in the astrometry. Although the new data alone cannot rule-out the presence of a candidate, when combined with published radial velocity measurements, the False Alarm Probabilities of the best solutions grow to unacceptable levels strongly suggesting that the observed astrometric wobble is not due to an unseen companion.

The Doppler technique measures the reflex radial motion of a star induced by the presence of companions and is the most successful method to detect exoplanets. If several planets are present, their signals will appear combined in the radial motion of the star, leading to potential misinterpretations of the data. Specifically, two planets in 2:1 resonant orbits can mimic the signal of a single planet in an eccentric orbit. We quantify the implications of this statistical degeneracy for a representative sample of the reported single exoplanets with available datasets, finding that 1) around 35 percent of the published eccentric one-planet solutions are statistically indistinguishable from planetary systems in 2:1 orbital resonance, 2) another 40 percent cannot be statistically distinguished from a circular orbital solution and 3) planets with masses comparable to Earth could be hidden in known orbital solutions of eccentric super-Earths and Neptune mass planets.

We are undertaking an astrometric search for gas giant planets and brown dwarfs orbiting nearby low mass dwarf stars with the 2.5-m du Pont telescope at the Las Campanas Observatory in Chile. We have built two specialized astrometric cameras, the Carnegie Astrometric Planet Search Cameras (CAPSCam-S and CAPSCam-N), using two Teledyne Hawaii-2RG HyViSI arrays, with the cameras' design having been optimized for high accuracy astrometry of M dwarf stars. We describe two independent CAPSCam data reduction approaches and present a detailed analysis of the observations to date of one of our target stars, NLTT 48256. Observations of NLTT 48256 taken since July 2007 with CAPSCam-S imply that astrometric accuracies of around 0.3 milliarcsec per hour are achievable, sufficient to detect a Jupiter-mass companion orbiting 1 AU from a late M dwarf 10 pc away with a signal-to-noise ratio of about 4. We plan to follow about 100 nearby (primarily within about 10 pc) low mass stars, principally late M, L, and T dwarfs, for 10 years or more, in order to detect very low mass companions with orbital periods long enough to permit the existence of habitable, Earth-like planets on shorter-period orbits. These stars are generally too faint and red to be included in ground-based Doppler planet surveys, which are often optimized for FGK dwarfs. The smaller masses of late M dwarfs also yield correspondingly larger astrometric signals for a given mass planet. Our search will help to determine whether gas giant planets form primarily by core accretion or by disk instability around late M dwarf stars.

According to current plans, the SIM/NASA mission will be launched just after the end of operations for the Gaia/ESA mission. This is a new situation which enables long term astrometric projects that could not be achieved by either mission alone. Using the well-known perspective acceleration effect on astrometric measurements, the true heliocentric radial velocity of a nearby star can be measured with great precision if the time baseline of the astrometric measurements is long enough. Since white dwarfs are compact objects, the gravitational redshift can be quite large (40-80 km/s), and is the predominant source of any shift in wavelength. The mismatch of the true radial velocity with the spectroscopic shift thus leads to a direct measure of the Mass--Radius relation for such objects. Using available catalog information about the known nearby white dwarfs, we estimate how many masses/gravitational redshift measurements can be obtained with an accuracy better than 2%. Nearby white dwarfs are relatively faint objects (10 < V < 15), which can be easily observed by both missions. We also briefly discuss how the presence of a long period planet can mask the astrometric signal of perspective acceleration.

Context. High accuracy astrometric instruments like Gaia aiming at an accuracy of 1 microarcsecond cannot be considered as point-like observers in the framework of relativistic modelling of observable quantities. Aims. Special-relativistic effects on the imaging by a non-point-like arbitrarily moving optical instrument are discussed. Methods. A special-relativistic reflection law for a mirror of arbitrary shape and motion is derived in the limit of geometrical optics. The aberration patterns are computed with ray tracing using a full pecial-relativistic model for two simple rotating optical instrument. Results. It was found that the effect of special-relativistic reflection law on the photocenters of aberration patterns of an optical system rotating with a moderate angular velocity of $60 \arcsec/{\rm s}$ may be at the level of 1 microarcsecond if the system involves mirrors significantly inclined relative to the optical axis. Conclusions. Special-relativistic optical modelling of the future astrometric instrument is generally speaking indispensable if a level of a few microarcseconds is envisaged.