The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

We report on the first star discovered to host a planet detected by radial velocity (RV) observations obtained within the CARMENES survey for exoplanets around M dwarfs. HD 147379 ($V = 8.9$ mag, $M = 0.58 \pm 0.08$ M$_{\odot}$), a bright M0.0V star at a distance of 10.7 pc, is found to undergo periodic RV variations with a semi-amplitude of $K = 5.1\pm0.4$ m s$^{-1}$ and a period of $P = 86.54\pm0.06$ d. The RV signal is found in our CARMENES data, which were taken between 2016 and 2017, and is supported by HIRES/Keck observations that were obtained since 2000. The RV variations are interpreted as resulting from a planet of minimum mass $m_{\rm p}\sin{i} = 25 \pm 2$ M$_{\oplus}$, 1.5 times the mass of Neptune, with an orbital semi-major axis $a = 0.32$ au and low eccentricity ($e < 0.13$). HD 147379b is orbiting inside the temperate zone around the star, where water could exist in liquid form. The RV time-series and various spectroscopic indicators show additional hints of variations at an approximate period of 21.1d (and its first harmonic), which we attribute to the rotation period of the star.

On March 20, 2015, we obtained 159 spectra of the Sun as a star with the solar telescope and the FTS at the Institut f\"ur Astrophysik G\"ottingen, 76 spectra were taken during partial solar eclipse. We obtained RVs using $I_2$ as wavelength reference and determined the RM curve with a peak-to-peak amplitude of almost 1.4 km s$^{-1}$ at typical RV precision better than 1 m s$^{-1}$. We modeled disk-integrated solar RVs using surface velocities, limb darkening, and information about convective blueshift from 3D magneto-hydrodynamic simulations. We confirm that convective blueshift is crucial to understand solar RVs during eclipse. Our best model reproduced the observations to within a relative precision of 10% with residuals less than 30 m s$^{-1}$. We cross-checked parameterizations of velocity fields using a Dopplergram from the Solar Dynamics Observatory and conclude that disk-integration of the Dopplergram does not provide correct information about convective blueshift necessary for m s$^{-1}$ RV work. As main limitation for modeling RVs during eclipses, we identified limited knowledge about convective blueshift and line shape as functions of solar limb angle. We suspect that our model line profiles are too shallow at limb angles larger than $\mu = 0.6$ resulting in incorrect weighting of the velocities across the solar disk. Alternative explanations cannot be excluded like suppression of convection in magnetic areas and undiscovered systematics during eclipse observations. Accurate observations of solar line profiles across the solar disk are suggested. We publish our RVs taken during solar eclipse as a benchmark curve for codes calculating the RM effect and for models of solar surface velocities and line profiles.

Using a photonic reformatter to eliminate the effects of conventional modal noise could greatly improve the stability of a high resolution spectrograph. However the regimes where this advantage becomes clear are not yet defined. Here we will look at where modal noise becomes a problem in conventional high resolution spectroscopy and what impact photonic spectrographs could have. We will theoretically derive achievable radial velocity measurements to compare photonic instruments and conventional ones. We will discuss the theoretical and experimental investigations that will need to be undertaken to optimize and prove the photonic reformatting concept.

Radial velocity observations of stars are entering the sub-m/s domain revealing fundamental barriers for Doppler precision experiments. Observations of the Sun as a star can easily overcome the m/s photon limit but face other obstacles. We introduce the G\"ottingen Solar Radial Velocity Project with the goal to obtain high precision (cm/s) radial velocity measurements of the Sun as a star with a Fourier Transform Spectrograph. In this first paper, we present the project and first results. The photon limit of our 2 min observations is at the 2 cm/s level but currently limited by strong instrumental systematics. A drift of a few m/s per h is visible in all observing days probably caused by vignetting of the solar disk in our fiber coupled setup, and imperfections of our guiding system adds further offsets in our data. Binning the data into 30 min groups shows m/s stability after correcting for a daily and linear instrumental trend. Our results show the potential of Sun-as-a-star radial velocity measurements that can possibly be achieved after a substantial upgrade of our spectrograph coupling strategy. Sun-as-a-star observations can provide crucial empirical information about the radial velocity signal of convective motion and stellar activity, and on the wavelength dependence of radial velocity signals caused by stellar line profile variations.

We present a new solar flux atlas with the aim to understand wavelength precision and accuracy in solar benchmark data. The atlas covers the wavelength range 405--2300 nm and was observed at the Institut f\"ur Astrophysik, G\"ottingen (IAG) with a Fourier Transform Spectrograph. In contrast to other FTS atlases, the entire visible wavelength range was observed simultaneously using only one spectrograph setting. We compare the wavelength solution of the new atlas to the Kitt Peak solar flux atlases and to the HARPS frequency-comb calibrated solar atlas. Comparison reveals systematics in the two Kitt Peak FTS atlases resulting from their wavelength scale construction, and shows consistency between the IAG and the HARPS atlas. We conclude that the IAG atlas is precise and accurate on the order of $\pm 10$ m s$^{-1}$ in the wavelength range 405--1065 nm while the Kitt Peak atlases show deviations as large as several ten to 100 m s$^{-1}$. We determine absolute convective blueshift across the spectrum from the IAG atlas and report slight differences relative to results from the Kitt Peak atlas that we attribute to the differences between wavelength scales. We conclude that benchmark solar data with accurate wavelength solution are crucial to better understand the effect of convection on stellar RV measurements, which is one of the main limitations of Doppler spectroscopy at m s$^{-1}$ precision.

The performance of highly-multiplexed spectrographs is limited by focal degradation (FRD) in the optical fibres. It has already been shown that this is caused mainly by processes concentrated around the mounting points at the ends of the fibres. We use the thickness of rings produced in the farfield when a fibre is illuminated by a collimated beam, to estimate the size of the region where the FRD is generated. This requires the development of a new model, using features of existing ray-tracing and wave-based models, which fits existing data very well. The results suggest that the amount of FRD is primarily determined by the length of fibre bonded into the supporting ferrule. We point out the implications for the production of future fibre systems.

Optical fibres are essential for many types of highly-multiplexed and precision spectroscopy. The success of the new generation of multifibre instruments under construction to investigate fundamental problems in cosmology, such as the nature of dark energy, requires accurate modellisation of the fibre system to achieve their signal-to-noise goals. Despite their simple construction, fibres exhibit unexpected behaviour including non-conservation of Etendue (Focal Ratio Degradation; FRD) and modal noise. Furthermore, new fibre geometries (non-circular or tapered) have become available to improve the scrambling properties that, together with modal noise, limit the achievable SNR in precision spectroscopy. These issues have often been addressed by extensive tests on candidate fibres and their terminations but these are difficult and time-consuming. Modelling by ray-tracing and wave analysis is possible with commercial software packages but these do not address the more complex features, in particular FRD. We use a phase-tracking ray-tracing method to provide a practical description of FRD derived from our previous experimental work on circular fibres and apply it to non-standard fibres. This allows the relationship between scrambling and FRD to be quantified for the first time. We find that scrambling primarily affects the shape of the near-field pattern but has negligible effect on the barycentre. FRD helps to homogenise the nearfield pattern but does not make it completely uniform. Fibres with polygonal cross-section improve scrambling without amplifying the FRD. Elliptical fibres, in conjunction with tapering, may offer an efficient means of image slicing to improve the product of resolving power and throughput but the result is sensitive to the details of illumination.

Fibre modal noise occurs in high spectral resolution, high signal-to-noise applications. It imposes fundamental limits on the photometric accuracy of state-of-the-art fibre-spectrograph systems. In order to maximize the performance of current and future instruments it is therefore essential to predict fibre modal noise. To attain a predictive model we are using a dual approach, bringing theoretical assumptions in line with the experimental data obtained using a test-bench spectrograph. We show that the task of noise prediction can be reduced to determining the visibility of the modal pattern which can be measured at the detector plane. Subsequently, the visibility-dependence of essential parameters is presented. This work will soon provide a basis for prediction of modal noise limitations in fibre-coupled spectrograph designs.