This report models the refractive index above a telescope site by a homogeneous atmosphere with exponential decay of the refractive index (susceptibility) as a function of altitude. The air is represented as a spherical hull with a symmetry center at the Earth center. We compute (i) the differential zenith angle -- the difference between the actual zenith angle at arrival of rays at the telescope and a hypothetical zenith angle without the atmosphere -- and (ii) the optical path length distribution of rays at arrival in the entrance pupil as a function of the distance to the pupil center. The key technique of the semi-numerical mathematics in this work is to expand some integrals -- that depend on the refractive index profile along the curved path of each ray from some virtual plane in the direction of the star up to the entrance pupil -- in power series of small parameters.

This works presents a program which computes velocities of an Earth-bound observatory in the reference frame of the barycenter of the solar system. It feeds from ephemerides files of the Jet Propulsion Laboratory to extract the velocity of the geocenter, optionally with corrections from Earth rotation data of the International Earth Rotation Service, takes a datum (time and geodetic location) of the observer as parameters, and processes these data with the program library of the working group `Standards of Fundamental Astronomy' of the International Astronomical Union. The prospective application of the computed velocities is to subtract their projection onto a pointing direction from observed velocities in a step of data reduction of astronomic radial velocities.

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.

The displacement of star images by atmospheric refraction observed by an Earth-bound telescope is dominated by a familiar term proportional to the product of the tangent of the zenith angle by the refractivity at the ground. The manuscript focuses on the torsion of the ray path through the atmosphere in a model of atmospheric layers above the ellipsoidal Earth surface, induced by the two slightly different principal curvatures along N--S and E--W pointing directions. This breaking of the azimuthal symmetry effects apparent places at the sub-milliarcsecond scale at optical and infrared wavelengths.

The polarization effect of the 31 reflections within the mirror train of an auxiliary telescope of the Very Large Telescope Interferometer is calculated as a function of pointing direction and rotator angle. With a rough estimate of the mean complex index of refraction of the reflecting surfaces, their Jones matrices are concatenated while tracing a ray from M1 up to the feeding mirror in front of an instrument. The net effect is summarized in terms of the axis ratio of the polarization ellipse of star light that was circularly polarized above the primary mirror.

The phase structure function accumulated by two parallel rays after transmission through a layer of turbulent air is best known by a proportionality to the 5/3rd power of the lateral distance in the aperture, derived from an isotropic Kolmogorov spectrum of the refractive index. For a von-Karman spectrum of the refractive index, a dependence involving a modified Bessel function of the ratio of the distance over the outer scale is also known. A further standard proposition is a proportionality to the path length through the atmospheric layer. The manuscript modifies this factor through a refined calculation of an integral representation of the structure function. The correction establishes a sub-linearity as the lateral distance grows in proportion to the layer thickness; it is more important for large than for small outer scales.

Estimators of outer scales of atmospheric turbulence usually fit the phase screen snapshots derived from local wave front sensors to a Zernike basis, and then compare the spectrum of expansion coefficients in this basis with a narrowing associated with decreasing outer scales. This manuscript discusses aspects that arise if the Zernike basis is exchanged for a Karhunen-Loeve basis of statistically independent modes. Data acquisition turns out to be more demanding because sensing the tip-tilt mode is required. The data reduction methodology is replaced by fitting of variance ratios of an entire (long exposure) data set. Statistical testing of hypotheses on outer scale models can be applied to a set of modes - supposing other noise originating from detector readout and the optical train can be disentangled.

Phase screens above a telescope pupil represent the variation of the phase of the electromagnetic field induced by atmospheric turbulence. Instances drawn from such statistics are represented by a vector of random phase amplitudes which are coefficients of a linear superposition of two-dimensional basis functions across the pupil. This work shortly reviews Fried's analysis of this modal decomposition for the case of Kolmogorov statistics of the phase covariance as a function of separation in the pupil plane. We focus on the numerical synthesis of phase screens. The statistically independent modes are transformed into the eigen-modes of a gradient matrix as time-dependence is introduced such that on short time scales the instances of the phase screens are rigidly shifted into a direction imposed by some wind velocity - known as the Taylor frozen screen approximation. This simple technique factorizes spatial and temporal variables and aims at binding the time dependence of the phase screens to the few expansion coefficients of the basis functions that obey a stochastic time-dependent differential equation.

Baseline calibration of a stellar interferometer is a prerequisite to data reduction of astrometric operations. This technique of astrometry is triangulation of star positions. Since angles are deduced from the baseline and delay side of these triangles, length and pointing direction (in the celestial sphere) of the baseline vector at the time of observation are key input data. We assume that calibration follows from reverse astrometry; a set of calibrator stars with well-known positions is observed and inaccuracies in these positions are leveled by observing many of them for a common best fit. The errors in baseline length and orientation angles drop proportional to the inverse square roots of the number of independent data taken, proportional to the errors in the individual snapshots of the delay, and proportional to the errors in the apparent positions of the calibrators. Scheduling becomes important if the baseline components are reconstructed from the sinusoidal delay of a single calibrator as a function of time.

Adaptive optics will almost completely remove the effects of atmospheric turbulence at 10 microns on the Extremely Large Telescope (ELT) generation of telescopes. In this paper, we observationally confirm that the next most important limitation to image quality is atmospheric dispersion, rather than telescope diffraction. By using the 6.5 meter MMT with its unique mid-IR adaptive optics system, we measure atmospheric dispersion in the N-band with the newly commissioned spectroscopic mode on MIRAC4-BLINC. Our results indicate that atmospheric dispersion is generally linear in the N-band, although there is some residual curvature. We compare our measurements to theory, and make predictions for ELT Strehls and image FHWM with and without an atmospheric dispersion corrector (ADC). We find that for many mid-IR applications, an ADC will be necessary on ELTs.

The basic geometry of a stellar interferometer with two telescopes consists of a baseline vector and a direction to a star. Two derived vectors are the delay vector, and the projected baseline vector in the plane of the wavefronts of the stellar light. The manuscript deals with the trigonometry of projecting the baseline further outwards onto the celestial sphere. The position angle of the projected baseline is defined, measured in a plane tangential to the celestial sphere, tangent point at the position of the star. This angle represents two orthogonal directions on the sky, differential star positions which are aligned with or orthogonal to the gradient of the delay recorded in the u-v plane. The North Celestial Pole is chosen as the reference direction of the projected baseline angle, adapted to the common definition of the "parallactic" angle.

We discuss the effect of atmospheric dispersion on the performance of a mid-infrared adaptive optics assisted instrument on an extremely large telescope (ELT). Dispersion and atmospheric chromaticity is generally considered to be negligible in this wavelength regime. It is shown here, however, that with the much-reduced diffraction limit size on an ELT and the need for diffraction-limited performance, refractivity phenomena should be carefully considered in the design and operation of such an instrument. We include an overview of the theory of refractivity, and the influence of infrared resonances caused by the presence of water vapour and other constituents in the atmosphere. `Traditional' atmospheric dispersion is likely to cause a loss of Strehl only at the shortest wavelengths (L-band). A more likely source of error is the difference in wavelengths at which the wavefront is sensed and corrected, leading to pointing offsets between wavefront sensor and science instrument that evolve with time over a long exposure. Infrared radiation is also subject to additional turbulence caused by the presence of water vapour in the atmosphere not seen by visible wavefront sensors, whose effect is poorly understood. We make use of information obtained at radio wavelengths to make a first-order estimate of its effect on the performance of a mid-IR ground-based instrument. The calculations in this paper are performed using parameters from two different sites, one `standard good site' and one `high and dry site' to illustrate the importance of the choice of site for an ELT.

The statistically independent Karhunen-Loeve modes of refractive indices with isotropic Kolmogorov spectrum of the covariance are calculated in a sphere of given radius, rendered as series of 3D Zernike functions. Many of the symmetry arguments of the associated 2D problem for the circular input pupil remain valid. The technique of efficient diagonalization of the eigenvalue problem in wavenumber space is founded on the Fourier representation of the 3D Zernike basis.

Coupling efficiencies of an electromagnetic field with a Kolmogorov phase statistics into a step-index fiber in its monomode regime of wavelengths are computed from the overlap integral between the phase screens and the far-field of the monomode at infrared wavelengths. The phase screens are composed from Karhunen-Loeve basis functions, optionally cutting off some of the eigenmodes of largest eigenvalue as if Adaptive Optics had corrected for some of the perturbations. The examples are given for telescope diameters of 1 and 1.8 m, and Fried parameters of 10 and 20 cm. The wavelength of the stellar light is in the J, H, or K band of atmospheric transmission, where the fiber core diameter is tailored to move the cutoff wavelength of the monomode regime to the edges of these bands.

Long-baseline interferometry detects fringes created by superposition of two beams of light collected by two telescopes pointing into a common direction. The external path difference generated by pointing away from the zenith is commonly compensated by adding a variable optical path length (delay) through air for one beam such that the optical path difference between the beams remains close to zero near the detector. The ABCD formula assigns a (wrapped) phase to the amplitudes A to D of an interference pattern shifted by multiples of 90 degrees in phase. We study the interplay between a wide band pass of the optics and the dispersion of the air in the compensating delay, which leads to small deviations between the ABCD phase and the reduced, monochromatic group-delay representation of the wave packets. In essence, this adds dispersion to the effects that have been discussed for evacuated interferometers (telescopes in space) before [J. Opt. Soc. Am. A 22 (2005) 2774].

The basic effect of the earth's atmospheric refraction on telescope operation is the reduction of the true zenith angle to the apparent zenith angle, associated with prismatic aberrations due to the dispersion in air. If one attempts coherent superposition of star images in ground-based interferometry, one is in addition interested in the optical path length associated with the refracted rays. In a model of a flat earth, the optical path difference between these is not concerned as the translational symmetry of the setup means no net effect remains. Here, I evaluate these interferometric integrals in the more realistic arrangement of two telescopes located on the surface of a common earth sphere and point to a star through an atmosphere which also possesses spherical symmetry. Some focus is put on working out series expansions in terms of the small ratio of the baseline over the earth radius, which allows to bypass some numerics which otherwise is challenged by strong cancellation effects in building the optical path difference.