The MINiature Exoplanet Radial Velocity Array (MINERVA) is a dedicated observatory of four 0.7m robotic telescopes fiber-fed to a KiwiSpec spectrograph. The MINERVA mission is to discover super-Earths in the habitable zones of nearby stars. This can be accomplished with MINERVA's unique combination of high precision and high cadence over long time periods. In this work, we detail changes to the MINERVA facility that have occurred since our previous paper. We then describe MINERVA's robotic control software, the process by which we perform 1D spectral extraction, and our forward modeling Doppler pipeline. In the process of improving our forward modeling procedure, we found that our spectrograph's intrinsic instrumental profile is stable for at least nine months. Because of that, we characterized our instrumental profile with a time-independent, cubic spline function based on the profile in the cross dispersion direction, with which we achieved a radial velocity precision similar to using a conventional "sum-of-Gaussians" instrumental profile: 1.8 m s$^{-1}$ over 1.5 months on the RV standard star HD 122064. Therefore, we conclude that the instrumental profile need not be perfectly accurate as long as it is stable. In addition, we observed 51 Peg and our results are consistent with the literature, confirming our spectrograph and Doppler pipeline are producing accurate and precise radial velocities.

We present the result of calculations to optimize the search for molecular oxygen (O2) in Earth analogs transiting around nearby, low-mass stars using ground-based, high-resolution, Doppler shift techniques. We investigate a series of parameters, namely spectral resolution, wavelength coverage of the observations, and sky coordinates and systemic velocity of the exoplanetary systems, to find the values that optimize detectability of O2. We find that increasing the spectral resolution of observations to R = 300,000 - 400,000 from the typical R ~ 100,000, more than doubles the average depth of O2 lines in planets with atmospheres similar to Earth's. Resolutions higher than about 500,000 do not produce significant gains in the depths of the O2 lines. We confirm that observations in the O2 A-band are the most efficient except for M9V host stars, for which observations in the O2 NIR-band are more efficient. Combining observations in the O2 A, B, and NIR -bands can reduce the number of transits needed to produce a detection of O2 by about 1/3 in the case of white noise limited observations. However, that advantage disappears in the presence of typical levels of red noise. Therefore, combining observations in more than one band produces no significant gains versus observing only in the A-band, unless red-noise can be significantly reduced. Blending between the exoplanet's O2 lines and telluric O2 lines is a known problem. We find that problem can be alleviated by increasing the resolution of the observations, and by giving preference to targets near the ecliptic.

We present a novel implementation for extremely high resolution spectroscopy using custom-designed Fabry Perot Interferometer (FPI) arrays. For a given telescope aperture at the seeing limited case, these arrays can achieve resolutions well in excess of ${\rm R\sim10^5}$ using optical elements orders of magnitude smaller in size than standard echelle spectrographs of similar resolution. We apply this method specifically to the search for molecular oxygen in exoplanetary atmospheres using the ${\rm O_2}$ A-band at 0.76 ${\rm \mu m}$, and show how a FPI array composed of $\sim10$ etalons with parameters optimized for this science case can record ${\rm R=3-5\,\cdot10^5}$ spectra covering the full ${\rm O_2}$ A-band. Using simulated observations of the atmosphere of a transiting nearby Earth-like planet, we show how observations with a FPI array coupled to a long-slit spectrograph can reduce the number of transit observations needed to produce a ${\rm 3\sigma}$ detection of ${\rm O_2}$ by $\sim30\%$ compared to observations with a ${\rm R=10^5}$ echelle spectrograph. This, in turn, leads to a decrease in an observing program duration of several years. The number of transits needed for a ${\rm 3\sigma}$ detection can be further reduced by increasing the efficiency of FPI arrays using dualons (an etalon with a buried reflective layer), and by coupling the FPI array to a dedicated spectrograph optimized for the ${\rm O_2}$ A-band.