21 April 2015 Miniature Exoplanet Radial Velocity Array I: design, commissioning, and early photometric results
Author Affiliations +
J. of Astronomical Telescopes, Instruments, and Systems, 1(2), 027002 (2015). doi:10.1117/1.JATIS.1.2.027002
The Miniature Exoplanet Radial Velocity Array (MINERVA) is a U.S.-based observational facility dedicated to the discovery and characterization of exoplanets around a nearby sample of bright stars. MINERVA employs a robotic array of four 0.7-m telescopes outfitted for both high-resolution spectroscopy and photometry, and is designed for completely autonomous operation. The primary science program is a dedicated radial velocity survey and the secondary science objective is to obtain high-precision transit light curves. The modular design of the facility and the flexibility of our hardware allows for both science programs to be pursued simultaneously, while the robotic control software provides a robust and efficient means to carry out nightly observations. We describe the design of MINERVA, including major hardware components, software, and science goals. The telescopes and photometry cameras are characterized at our test facility on the Caltech campus in Pasadena, California, and their on-sky performance is validated. The design and simulated performance of the spectrograph is briefly discussed as we await its completion. New observations from our test facility demonstrate sub-mmag photometric precision of one of our radial velocity survey targets, and we present new transit observations and fits of WASP-52b—a known hot-Jupiter with an inflated radius and misaligned orbit. The process of relocating the MINERVA hardware to its final destination at the Fred Lawrence Whipple Observatory in southern Arizona has begun, and science operations are expected to commence in 2015.
Swift, Bottom, Johnson, Wright, McCrady, Wittenmyer, Plavchan, Riddle, Muirhead, Herzig, Myles, Blake, Eastman, Beatty, Barnes, Gibson, Lin, Zhao, Gardner, Falco, Criswell, Nava, Robinson, Sliski, Hedrick, Ivarsen, Hjelstrom, de Vera, and Szentgyorgyi: Miniature Exoplanet Radial Velocity Array I: design, commissioning, and early photometric results



The field of exoplanetary science has experienced rapid growth since the discoveries of the first planetary-mass companions more than two decades ago. These initial discoveries spawned myriad observational efforts that have expanded our view of planetary systems from a single example—our own Solar System—to a diverse statistical ensemble containing hundreds of confirmed systems and thousands of candidates.78.9.10 Modern techniques for discovering and characterizing exoplanets include transits,1112.13.14.15 microlensing,1617.18 and direct detection.19,20 Our understanding of planet formation and evolution and the possibility for other intelligent life in the cosmos has been transformed by this swift progress that continues to accelerate in exciting directions.

At the time of the launch of the Kepler Mission in 2009,21 radial velocity (RV) surveys had discovered more than 400 planets orbiting nearby stars by detecting the minute, periodic Doppler shifts in stellar spectra induced by orbiting planetary companions. The diverse collection of RV-detected planets revealed many important correlations between planet occurrence and stellar properties,2223.24 as well as detailed information about the physical and orbital characteristics of planets outside our Solar System (for a few interesting examples, see Refs. 2526.27). We now know that planets throughout the Galaxy have a wider range of masses, radii, and internal structures than the planets of our own Solar System. Of particular interest is a new class of planet—the so-called super-Earths—with masses and radii intermediate to the Solar System terrestrial planets and the ice-giants Neptune and Uranus (see Refs. 22 and 2829.30.31 for some early examples).

From the standpoint of searching for low-mass exoplanets that resemble those of the inner Solar System, one of the most exciting statistical results from RV-detected planets is a planetary mass function that rises steeply toward Earth-mass planets. Howard et al.32 analyzed the planet discoveries and detection efficiency of the NASA/UC η Survey conducted at Keck Observatory using the HIRES spectrograph and found the number of planets per interval in MminMpsini scales as dN/dlogMminMmin0.48. Extrapolation of this relationship to terrestrial minimum masses (0.5<Mpsini/M<2.0) and periods P<50 days led to the remarkable prediction that 23% of Sun-like stars harbor an Earth-mass planet. This number agrees well with results from the HARPS surveys of chromospherically quiet FGK dwarfs33,34 and M dwarfs.35 Small planets are far more common than large planets throughout the Galaxy, which bodes well for the search for Earth-like planets in the Solar Neighborhood.

In the past four years, >4000 planet candidates36 detected by the Kepler Space Mission3738.39.40 have strengthened the scientific results from RV surveys and greatly expanded our knowledge of exoplanet properties down to sizes comparable to, and even below the Earth.41,42 Consistent with the results of RV surveys, the Kepler discoveries strikingly illustrate that the number of planets increases rapidly with decreasing planet radius.4344.45 There is still uncertainty about how much this trend softens or even turns over toward the smallest detected planets.46,47 However, it is clear that there are more planets smaller than 4 R than larger ones in the Galaxy. Estimates of the occurrence of Earth-size planets around Sun-like stars range from 10% to 15%,45,48 and the occurrence of Earth-size planets in Earth-like orbits is estimated to be between 2% and 6%.47,48

Less than a decade ago, the only known terrestrial planets orbiting main sequence stars resided in our Solar System, and there was little expectation that, even if they existed around other stars, they would be presently discovered. Now Kepler has shown us that they are very common and may even constitute the dominant population of exoplanets, particularly if one considers recent planet occurrence estimates around low-mass, M-type dwarfs, which are the most numerous stars in the Galaxy.46,49,50 While the characteristic distance of a Kepler target is 1kpc, the statistical results from the Kepler Mission should extend to the Solar Neighborhood, thereby informing us about the closest stars to Earth (but see Ref. 51). It is now clear that the night sky is teeming with unseen terrestrial-mass planets and super-Earths. The proximity of these low-mass planets in the Solar Neighborhood will facilitate follow-up studies that would be difficult or impossible with Kepler stars, and their physical properties will inform the search for life outside of the Solar System. Indeed, a small fraction of these nearby planetary systems have already been discovered.5253.54.55

Exoplanet transit searches require a near-perfect alignment of the orbital plane along the line of sight, necessitating fairly large sample sizes to ensure a detection. For example, the transit probability of a super-Earth orbiting a Sun-like star in an Earth-like orbit is 0.5%. This means a sample size of at least several hundred would be needed to ensure a single detection of one such planet. The limited number of local stars thus makes this method unfavorable for discovering large numbers of nearby planets.

RV surveys are more promising, as the detection probability is less sensitive to orbital inclination. However, the velocity precision required for detecting small planets is just at or beyond the limits of most current instruments. Equally important to the success of such an RV survey is the tremendous cadence and phase coverage needed to densely sample a planet’s full orbit and a range of stellar noise sources.54,56 Attaining this cadence for more than a small handful of stars is not realistic within the framework of shared telescope time allocation. Expanding our planetary census to dozens of stars in the Solar Neighborhood requires a dedicated observatory capable of highly precise RV measurements.

Some fraction of nearby planets discovered via the RV method will transit their host star, with an increase in transit probability for planets in shorter orbital periods. The RV data can, therefore, be used to guide searches of the transit windows of low-mass planets with precise photometry. This RV-first, transit-second method has proven to be a powerful observational technique producing the detection of the first transiting planet in 2000 (HD 209458 b)14,15 and providing us with the brightest transiting systems known.57,58 These bright (V<8) systems are important because they are the most amenable to follow-up science using space and ground-based facilities.

In addition to the short-period (P10 days) transiting planets around nearby stars, there should also exist a large population of super-Earths that lie within their respective habitable zones.59,60 There have been a few examples from this population recently discovered that support this claim.55,61,62 Locating and characterizing these planets from the ground will populate the target lists of future space-based missions designed to produce direct images and spectra of Earth-like planets. In this way, the high-precision RV surveys of today will be an important stepping stone toward discovering a true analog to our own Earth.

We aim to address the need for new Doppler-based planet detection facilities by building a dedicated ground-based observatory for the detection of small planets in our Solar Neighborhood called the Miniature Exoplanet Radial Velocity Array (MINERVA). The philosophy and design of MINERVA are presented in Sec. 2, which gives an overview of the project concept followed by a recent status update. The basic hardware that will be used to carry out MINERVA science is described in Sec. 3, which includes the telescopes, their enclosures, and the photometry cameras. A brief presentation of the MINERVA spectrograph, KiwiSpec-MINERVA, is presented in this section; the design details and on-sky performance will be presented in a forthcoming publication. The software that will run the MINERVA array is being adapted from the robotic brain of the Robo-AO system,63 and we summarize its basic architecture and functionality in Sec. 5. In Sec. 6, we present results validating the expected performance of our telescopes and cameras with photometric observations of one of our RV survey target stars, and we also present new observations and model fits to the transiting hot-Jupiter WASP-52 b in Sec. 7. Last, in Sec. 8, we briefly summarize this publication and offer a look toward future MINERVA opportunities.


Dedicated Facility for Exoplanet Detection and Characterization

MINERVA is designed around the primary science goal of detecting super-Earths within the habitable zones of nearby stars, as well as terrestrial-mass planets in close-in orbits. For this latter class, the planet candidate orbiting α Cen B is the prototype.54 Our strategy to achieve this goal is to build a dedicated observatory and perform a precise RV survey with a nominal three-year time baseline. We will target solar-type stars in the optical part of the spectrum based on optimizations performed by our team64 and the expertise of the project principle investigators. The relatively small number of targets, their brightness, their distribution on the sky, and our signal-to-noise ratio requirements to achieve our target precision of 0.8ms1 along with estimated cost and time line for construction are the primary considerations that determine the nature of the telescope(s) required for this program. Multiplexing is not feasible since the stars will, in general, be widely spaced on the sky, and to achieve our target RV precision, it would be inadvisable to risk crosstalk between multiple stellar spectra in a single exposure.

Large-aperture telescopes (D3m) have a relatively low duty cycle for these bright stars (3V8) because exposure times are short compared to the overheads incurred from CCD readout, slew time, and source acquisition. For example, the exposure time needed to achieve optimal signal-to-noise per pixel for a star with V=5 using the HIRES spectrometer on the Keck I 10-m telescope is 5 to 10 s, while the readout time is 42 s and a typical source acquisition time is 1 to 2 min. In addition, mitigating the effects of p-mode oscillations is best accomplished by matching exposure times to multiple oscillation periods, which in these stars is on the order of minutes.65 Indeed, early RV planet searches worked effectively with apertures from 0.6 to 3 m (two prominent examples are described in Refs. 6 and 66).

Another benefit of using small telescopes is cost. We have conducted an informal survey of telescope costs by obtaining list prices for high-end, commercially available small telescopes and cost estimates for larger, professional, custom telescopes up to 2.4 m through private consultation (the survey was conducted circa 2010 and may not reflect the current state of the market). The data shown in Fig. 1 reveal that the cost of small telescopes scales roughly with the aperture. There is a significant discontinuity in this scaling, however, between the largest commercial telescopes (1m) and slightly larger professional telescopes, of about a factor of four, presumably reflecting the economies of scale and less stringent engineering requirements for the amateur market. Interestingly, extrapolation of the custom scaling to very large telescopes approximately describes the cost of the Keck 10-m telescopes ($100million each) and the Thirty-Meter Telescope ($1billion). We also note that there appears to be a separate, lower track for truly mass-produced small telescopes (not shown in the figure).

Fig. 1

Cost curve for the purchase of telescopes of a specified aperture. The cost of amateur telescopes seems to scale with aperture diameter (red points and curve). There is a discontinuity near an aperture of 1 m between the largest commercial telescopes and slightly larger professional telescopes that seem to follow a cost curve that scales as collecting area (blue points and line). The cost and effective apertures corresponding to successive numbers of CDK-700 0.7 m telescopes are shown in green.


Since cost scales with aperture, there is no reduction in hardware cost in buying a single telescope versus several small telescopes. We can take advantage of this fact to put MINERVA’s cost onto the high-end amateur track by using multiple commercial 0.7-m telescopes whose light can be combined to create a larger effective aperture. This allows us to construct a large light bucket by feeding our spectrograph with multiple, smaller light buckets. The exact model of our 0.7-m telescopes was chosen based on the specific features offered, such as two instrument ports per telescope and fast slew speed as well as the proximity of the manufacturer to our test facility (see Sec. 3.1).

In addition to the factor of four decrease in hardware cost over a single custom telescope, this design choice offers several advantages over using a single telescope. They can be purchased off-the-shelf, complete with control software, allowing a quicker path to on-sky commissioning and diminished development risks, which offset the increased complexity of organizing a telescope array. The smaller étendue of the optical systems translates to a smaller spectrograph that is easier to stabilize both mechanically and thermally. This reduces the cost of both the spectrograph and the facility needed to adequately stabilize the instrument. Last, the modularity of the MINERVA design offers several benefits, such as redundancy, the ability to change the scope of the project, and flexibility with our observing strategy.

The survey target list is drawn from the NASA/UC η sample made up of 166 nearby, chromospherically inactive stars currently monitored by Keck/HIRES for orbiting exoplanets.67 The projected yield from this target list dictates the minimum effective aperture required for MINERVA. We choose from this list the maximum number of targets that can be observed from the final location of the array in southern Arizona to the precision necessary to detect planets with Mpsini=3M in their respective habitable zones with three years of observations.68 Figure 2 shows the η target list with the required integration time per night. As we do not have commissioning data for our spectrograph in hand, we perform simulations assuming a hard RV precision limit of 0.8ms1 and a photon noise model. We assume a total system throughput of 10% (see Secs. 3.4, 6.1, and 6.4).

Fig. 2

The required integration time per night for the MINERVA array to detect three M planets within the habitable zones of each star in the η sample according to a photon limited noise model as a function of V magnitude. Data points are colored according to their effective temperatures, determined using stellar masses from Ref. 32.


Our simulations account for stellar jitter, target observability, and weather losses based on historical weather records for southern Arizona. The RV stellar jitter is modeled to match the spot-induced photometric variations observed in quiet G and K dwarfs by Kepler.69 We find that a spot model with four to seven spot pairs per star with sizes ranging from 1.4% to 1.8% of the stellar radius reproduce the observed photometric variations of between 3×104 and 7×104 well. We populate the surfaces of our simulated stars with spots having lifetimes and latitudes following the Solar butterfly diagram and assign rotation periods based on the period distribution seen by Kepler.70 The spot-induced pseudo-Doppler shifts are added to the simulated dynamical shifts caused by the planet. We then add white noise to the simulated RVs, recording the reduction in detectability of the simulated planet. The maximum amount of white noise added such that the planet was detected in 99% of the realizations was taken as the required per-night precision per star. Our final target list is chosen in consideration of the length of nights, declination of target stars, calibration observations, 10% overhead for secondary science and other programs, and a 25 s telescope slew and source acquisition time.

The projected exoplanet yield from the MINERVA Doppler survey is estimated using the statistical results from the Kepler Mission.71 We extrapolated the reported occurrence rates out to periods of 400 days for planets above 2 R and assume the same frequencies per log bin as those for 85-day periods. The latter is likely a conservative assumption, as the frequency of planets appears to rise beyond 50-day periods. For each target in the MINERVA sample, we randomly drew planets in the radius-period grid based on the extrapolated frequency surface and then converted the exoplanet radius to a mass using a density relation, MR2.29 (Ref. 71). Signals >3ms1 would have already been detected by the η program and are not included in the yield. Multiplanet extractions were not performed in this simulation and may delay confirmation of the largest RV signals in multiplanet systems.

The final yield results are obtained from the results of repeating the simulation 1000 times. Using an effective aperture of 1.4 m (four CDK-700 telescopes), we are able to observe the 82 brightest stars from the η list, and we find a mean yield of 15±4 planets with amplitudes between 0.8 and 3ms1. Improvements to these estimates are currently being pursued with the use of a more realistic stellar activity model and optimized observing strategies based on our recovery methods. Modeling of the stellar activity cycles will be necessary for the lowest-mass planets, and we recognize that the efficacy of these algorithms, which have yet to be quantified, will affect our estimated yield.

However, the ability to measure Doppler shifts of our targets at this precision every night (weather permitting) is unique to MINERVA and will be a key factor in recovering the RV signals of low-mass planets in their respective habitable zones. The observing cadence achievable with MINERVA allows us to account for the stellar variability of our sample of stars in a way that current facilities cannot. Based on our simulations and our expectations from a more refined treatment of stellar variability, we find a 1.28 m aperture to be the minimum required for ensured success of the project. We, therefore, design MINERVA around the use of four 0.7-m telescopes.

The light from each telescope will be fed into 50-μm octagonal fibers using a custom focal plane unit (see Ref. 72 and Sec. 3.3). These four fibers will then form a pseudo-slit at the entrance of the MINERVA spectrograph. KiwiSpec-MINERVA is a highly stabilized, bench mounted spectrograph covering wavelengths from 500 to 630 nm with a resolving power of R80,000 (see Sec. 3.4). This is an optimal spectral range for RV surveys of solar-type stars64 that can also be wavelength calibrated with an iodine cell.73

Built in to the design of MINERVA is the capability for flexible scheduling and simultaneous science and education programs. Of the 15±4 simulated detections described above, 1.0±0.8 are expected to transit their host star. Currently, there are 16 RV detected planets with declinations δ>20deg, periods <30 days, Msini<50 M, and V<10. With an additional 10 from MINERVA, the total transit yield is expected to exceed unity. This motivates the secondary science objective of MINERVA to search for transits of super-Earths among its RV discovered planets and to further characterize known transiting planets with multiband light curves. This requires a broadband photometry precision of <1mmag in the optical. We demonstrate a comparable level of photometric precision from our commissioning site on the Caltech campus in Pasadena, California (see Sec. 6). On the educational front, students in lab courses can use one of the telescopes to conduct their course assignment and gain valuable observing experience—a community need recently expressed by Ref. 74—while the other telescopes simultaneously conduct the primary science program.

During the early stages of the project, MINERVA will be used to follow up newly found Jupiter- and Neptune-sized planets from surveys like HATNet,75 WASP,76 and NGTS77 in addition to providing long-term monitoring for some TERMS long period planets.78 MINERVA photometry can also be used to follow up space-based discoveries. The Kepler K2 mission79 will produce thousands of transit discoveries, but will only monitor each target field for 75 days. At later times, the number of potential targets will grow considerably as the Transiting Exoplanet Survey Satellite (TESS) Mission,80 set to launch in 2017, will yield many more detections over the whole sky, but with continuous monitoring of only 27 days for each nonoverlapping field.


Project Status and Approximate Timeline

As of the writing of this manuscript, MINERVA telescopes 1 and 2 are located on the Caltech campus in the first Aqawan (see Sec. 3.2), where we have been developing the MINERVA Robotic Software (MRS) and validating the performance of the various hardware components. The hardware on-site has been fully tested, and now the goal of this facility is to achieve coordinated, automated control of both telescopes and Aqawan 1. This is expected to be complete by early 2015, and once this goal has been reached, the entire facility will be moved to Mt. Hopkins.

The infrastructure at Mt. Hopkins necessary for the relocation of the Aqawans and telescopes has been completed. Aqawan 2 has been constructed in California and was transported to Mt. Hopkins on December 9, 2014. MINERVA telescopes 3 and 4 have had their performance validated and were relocated to Aqawan 2 at Mt. Hopkins on December 15, 2014.

The custom room designed for the spectrograph is currently under construction. The outermost layer in a three-stage environmental control scheme will be a 100 K clean room temperature stabilized to ±1°C, peak-to-peak. A second, interior room will then be erected inside which the spectrograph will be mounted with its critical optical elements inside a vacuum chamber. The manufacturing of the spectrograph has been completed by Callaghan Innovation in New Zealand. Work is now underway on the input optics as well as the iodine cell mount. It is undergoing lab tests and awaits the completion of the spectrograph room at Mt. Hopkins. Delivery of the spectrograph is expected by the end of the first quarter of 2015 when on-sky commissioning will begin.

Allowing for minor unforeseen delays, fully automated, robotic control of the array and spectrograph is expected by mid-year 2015. The primary survey is projected to begin by the end of the year.


MINERVA Hardware


CDK-700 by PlaneWave

The PlaneWave CDK-700 is a 0.7-m, altitude/azimuth mounted telescope system.81,82 It has a compact design, standing just <8ft tall when pointing at the zenith, with a 5-ft radius of maximum extent when pointing horizontally. The telescopes use a corrected Dall-Kirkham (CDK) optical setup consisting of an elliptical primary mirror, a spherical secondary mirror, and a pair of correcting lenses to remove off-axis coma, astigmatism, and field curvature. This results in a flatter, more coma and astigmatism-free field than the Ritchey-Chretien design, with the added benefit that the spherical secondary mirror makes alignment forgiving compared to the hyperbolic secondary of the Ritchey-Chretien design. The CDK-700 has dual Nasmyth ports with output beams at f/6.5 accessed with a rotating tertiary mirror. The CDK-700 specifications are summarized in Table 1.

Table 1

CDK-700 specifications.

Optical system
Optical designCorrected Dall-Kirkham (CDK)
Aperture700 mm (27.56 in.)
Focal length4540 mm
Focal ratio6.5
Central obscuration47% primary diameter
Back focus305 mm from mounting surface
Focus positionNasmyth (dual)
Dimensions93.73 arcsec⁡ H×43.25 arcsec  W× 39 arcmin⁡ D
Weight1200 lbs
Optical performance1.8  μm
Image scale22  μm
Optimal field of view70 mm (0.86 deg)
Fully baffled field60 mm
Mechanical structure
ForkMonolithic U-shaped fork arm
Azimuth bearing20 in. diameter thrust bearing
Altitude bearing2×8.5  in.
Optical tubeDual truss structure
Motion control
MotorsDirect drive, three phase axial flux torque motor
EncodersStainless steel encoder tape with 81 mas resolution
Motor torque∼35  ft-lbs
Slew rate15  deg s−1
System performance
Pointing accuracy10 arcsec RMS
Pointing precision2 arcsec RMS
Tracking accuracy1 arcsec RMS over 3 min
Field derotator3  μm


These values were validated during the commissioning. Other values are the manufacturer specifications.

The telescope pointing is controlled by two direct-drive electromagnetic motors with an encoder resolution of 81 mas resulting in a pointing accuracy of 10 arcsec RMS, a pointing precision of 2 arcsec, and a tracking accuracy of 1 arcsec over a 3-min period. The slew rate is 15degs1, which keeps slew times between any two points in the sky to <10s. The focus mechanism and image derotator are combined into a single, motor-controlled unit that can be remotely adjusted. Cooling fans and temperature sensors are used to keep the primary mirror in thermal equilibrium, and the control software is built to automatically correct for perturbations, such as wind gusts.

Section 6 presents the performance validation for our telescopes, including control, pointing, guiding, and throughput. The commissioning performed over the past two years at Caltech has proven the feasibility of using high-end, off-the-shelf hardware for professional astronomical research. Thereby, MINERVA has retired many of the risks involved in developing an observational facility from the ground up. Figure 3 shows an image of the MINERVA telescopes at the Caltech test site.

Fig. 3

Telescope 2, a PlaneWave CDK-700, is shown inside the MINERVA Aqawan telescope enclosure at the Caltech commissioning site.



Aqawan Telescope Enclosures

An Aqawan (Chumash native American word for “to be dry”) is a telescope enclosure developed by Las Cumbres Observatory Global Telescope for their 0.4-m telescopes specifically designed for remote, robotic operations around the world.83 The design offers full access to the sky, limiting the effects of dome seeing, and eliminates the need to coordinate dome slit positioning while maintaining a relatively small footprint. We have purchased two custom-built Aqawans with longer sides that can each accommodate two CDK-700s without any possibility of collision and can close safely with the telescopes pointed in any orientation. Stronger motors with a higher gear ratio were also installed to handle the heavier roof panels. Figure 4 shows the design and realization of our first Aqawan, which has been delivered to the Caltech campus for commissioning.

Fig. 4

(a) Design drawing of our custom Aqawan telescope enclosure with two PlaneWave CDK-700 telescopes inside. (b) MINERVA commissioning site on the Caltech campus showing the open Aqawan and telescopes 1 and 2 inside.


The Aqawan receives 208V/30 A three phase power that is converted to 24 V dc within its control panel. This power runs through an internal uninterruptible power supply (UPS) and powers a programmable automated controller (PAC) that controls all the functionality of the enclosure. In addition to basic opening and closing, the Aqawan has many auxiliary features, including a web camera that provides 360 deg coverage, a temperature and humidity sensor, fans to promote temperature equalization, fluorescent lighting, and a smoke alarm. Communication to the Aqawan PAC is established via TCP/IP, with commands consisting of ASCII strings. The Aqawan firmware is designed such that if a “heartbeat” command is not issued each minute, the roof automatically closes. This feature along with built-in backup power and the ability for the Aqawan to close with the telescopes in any configuration offers safety against power and connectivity failure modes.


Fiber and Fiber Coupling

Astronomical light collected from MINERVA’s four 0.7-m telescopes will feed a stabilized spectrograph via fiber optic cables. Using a fiber feed offers many advantages and some challenges. Single-mode fibers offer superior control of the instrumental profile. However, they are diffraction-limited by nature and require a high-performance adaptive optics system for efficient coupling of starlight. Multimodal fibers couple to starlight much more easily, but they have a near- and far-field output that is variable. The near-field variations are due to the interference of modes at the output of the fiber, and when imaged on the detector, these variations can limit the signal-to-noise ratio of the observations. This effect can be mitigated by use of an octagonal fiber, which mixes the modes as they traverse the fiber.84,85 Physical agitation of the fiber enhances this effect.86

The far-field uniformity and stability are also important as the far field is incident on the echelle grating. Different grooves will be illuminated if the spatial intensity distribution changes, which can introduce spurious wavelength shifts as the grooves are not identical. One way to improve this performance is to introduce a double scrambler, which inverts the field and angle distribution, at a cost of reduced throughput.87 This is an option that MINERVA will explore further if necessary.

We will use 50-μm octagonal fibers with a 94μm cladding diameter and a numerical aperture of 0.22 to feed light collected from each of the four 0.7-m telescopes to our stabilized spectrograph. Starlight will be coupled to our fibers at the native f/6.5 of the CDK-700 bypassing the need for small optics for the purpose of better throughput. KiwiSpec-MINERVA is designed to accommodate an f/5.0 beam or larger, permitting up to 20% focal ratio degradation before impacting overall throughput. Therefore, our fibers project onto the sky with a 2.27 arcsec diameter.

Our fiber coupling system consists of a fiber acquisition unit (FAU) and control software that can interface with the unit and provide closed-loop guiding with the telescope control. Besides KiwiSpec-MINERVA, this is one of the only truly custom hardware components of MINERVA. The details of this system are presented in Ref. 72, and here we present only a brief review.

The FAU has three accessible optical paths that are based upon the design presented in Ref. 88. In the primary path, the telescope beam is fed directly into the fiber. Separately, there is an optical path that relays a portion of the light to a guiding camera via a pellicle. Finally, there is a set of relay lenses and a corner retroreflector. If the fiber is illuminated from the exit end near the spectrograph, the input end of the fiber tip will be imaged on the guide camera. The corner retroreflector guarantees that misalignments of the optics do not affect the image position; this allows the determination of the pixel position on the guider that corresponds to the fiber tip and, hence, a setpoint for guiding. An annotated schematic of the FAU design is shown in Fig. 5, and a picture of one of the FAUs mounted on a CDK-700 can be seen in Fig. 6.

Fig. 5

Ray-trace illustration of the fiber acquisition unit (FAU), which measures 305 mm (12in.) across. Light enters from the right, delivered by the CDK-700 with an f/6.5 beam. A pellicle reflects 2% of the beam to an SBIG ST-i guide camera with a V-band filter, matching the wavelengths used for radial velocity observations. The remaining 98% of the beam comes to a focus on the fiber tip. We include an achromatic lens and corner cube on the opposite side of the pellicle such that the fiber tip can be imaged onto the guide camera via backillumination. This allows for quickly determining the location of the fiber tip on the guide camera for precise guiding.


Fig. 6

The MINERVA FAU mounted on one of the MINERVA telescopes with all major components labeled. The focuser/derotator will be disengaged for standard spectroscopic observations.




KiwiSpec-MINERVA is an adaptation of an existing spectrograph design89,90 that leverages the successes of existing facilities9192.93 and techniques.73,94 This critical component of our facility will be described in detail in a forthcoming publication. Here, we provide a brief description of KiwiSpec-MINERVA for the sake of completeness of this publication.

It is a bench mounted, fiber-fed spectrograph of asymmetric, white pupil design. Primary dispersion is achieved using an R4 echelle, and a volume phase holographic grism is used for cross-dispersion. Four distinct traces will be imaged on a 2k×2k detector covering a spectral range from 500 to 630 nm over 26 echelle orders. The resolving power of the spectrograph is R80,000. Two additional calibration fibers bracket the four science fibers and provide stable wavelength calibration with the use of a stabilized etalon wavelength source (as simulated in Fig. 7), or Thorium-Argon light. A subsection of a simulated echellogram created using ray-tracing techniques95 is presented in Fig. 7. Each science trace corresponds to the input from an independent telescope and fiber system.

Fig. 7

A subsection of a simulated MINERVA echellogram showing four science spectra bracketed by a simulated stabilized etalon wavelength source. Each science fiber trace represents the input from one of the four MINERVA telescopes that will be observing the same astronomical target.


The estimated throughput of the prototype spectrograph from which the KiwiSpec-MINERVA has been designed is 30%. This has been validated with on-sky measurements taken at Mt. John Observatory in New Zealand.90 For these tests, the mirrors for the prototype spectrograph were bare aluminum, the lenses had only single layer antireflection coatings, and several off-the-shelf optical components were used, all of which contribute to a suboptimal throughput. Improvements have now been implemented, including high-efficiency coatings and the installation of a new custom camera such that the throughput is expected to increase to 45%. The throughput of our fiber system is expected to be 70%, and the total throughput of our system up to the entrance slit of the spectrograph is expected to be 50% (see Secs. 6.1 and 6.4). Therefore, the total throughput of our system including losses due to sky extinction at the Fred Lawrence Whipple Observatory (FLWO) is expected to be 10% or greater.

The line spread function will be sampled with 3 pixels and a minimum of 4 null pixels will lie between each trace. The spectral extraction and data reduction pipelines are currently being developed. The strategy is to model the two-dimensional spectrum directly, i.e., all orders will be modeled simultaneously such that cross-contamination between orders and scattered light will be accounted for in the model. This approach is likely only possible in the high signal-to-noise regime in which MINERVA will be working. The fact that all telescopes will be observing the same target simultaneously will also help to mitigate the effects of cross-contamination between orders.

The spectrograph will be placed inside a purpose-built, two-stage room and the critical components will reside inside a vacuum chamber stabilized to ±0.01°C. An iodine cell will be mounted off the optical bench and can be inserted or removed from the optical path, allowing the option of simultaneous wavelength calibration of the science traces in the echellogram.


Cameras and Filters

Each MINERVA telescope will also be equipped with a wide-field CCD camera on one of the two Nasmyth ports available on the CDK-700. The array will incorporate three different camera models for its four telescopes. Two of the MINERVA telescopes will be equipped with identical Andor iKON-L cameras.96 These cameras have a back-illuminated sensor with wide band (BV) coating and 2048×2048 square 13.5μm pixels for a total chip size of 27.6 mm corresponding to a 20.9 arcmin field. Figure 8 shows an image of one of these cameras mounted on a MINERVA telescope. A third telescope will be equipped with an additional Andor iKON-L that is identical to the two described above except that it will contain a deep depletion sensor with fringe suppression (BR-DD). This camera is sensitive to light at wavelengths out to 1μm and can be used for precision photometry in the near-infrared (e.g., i, z, and Y). The iKON-L cameras come with a four-stage thermoelectric cooling system that can achieve operating temperatures 80°C, nearly eliminating dark current. The third camera model will be an Apogee Aspen CG23097 with 2048×2048 15×15μm2pixels for a chip size of 30.7 mm corresponding to a 23.2 arcmin field.

Fig. 8

The optics chain for MINERVA photometry with major components along the light path to the camera labeled.


Each MINERVA telescope has a filter wheel on the photometry port, and we employ two different Apogee filter wheel models, each with 50 mm square slots. Three of the MINERVA telescopes have the AFW50-7s filter wheel, while one telescope has a custom double filter wheel comprising two AFW50-10s wheels mounted back to back. One open slot in each 10-slot wheel allows access to the other 18 slots. Currently, the filters available to the MINERVA system are the Johnson U, B, V, R, and I; second-generation Sloan g, r, i, and z; narrow band H α, [Sii], [Oiii], and amateur filters, L, R, G, and B. The Andor cameras are mated to the filter wheels through a custom adaptor plate designed by Paul Gardner of Caltech Optical Observatories and implemented by Andor.


MINERVA Site: Fred Lawrence Whipple Observatory

The performance validation and preliminary commissioning of all MINERVA components aside from the spectrograph have been done at the commissioning site on the Caltech campus in Pasadena, California (see Fig. 4), or at the PlaneWave warehouse in Rancho Dominguez, California. Once the primary commissioning tasks are complete and fully automated observations have been achieved, the entire facility will be moved to its final location at FLWO on Mt. Hopkins outside of Amado, Arizona (see Sec. 2.1).

The FLWO site was chosen after a site selection study conducted over the summer of 2013. We accumulated historical data and visited with staff and personnel from Mt. Wilson Observatory in Los Angeles, California, McDonald Observatory in Jeff Davis County, Texas, and San Pedro Mártir Observatory in Baja California, Mexico, in addition to FLWO. Mt. Hopkins was determined to be the optimal choice for MINERVA based on good overall weather and seeing conditions, the existing infrastructure available for use, full-time support staff, and financial considerations in setting up and maintaining the site. The weather data for Mt. Hopkins consisted of a compilation of observing logs from the HAT-Net project75 located on the Mt. Hopkins Ridge 130m to the North, and data from the MEarth project98 located 430m to the North and at a slightly higher elevation than MINERVA (see Fig. 9). We also incorporated information presented in Ref. 99 as the relative conditions at the Ridge correlate well with the conditions at the summit where MMT Observatory is located. From these data, we anticipate 271 nights per year that we will be able to observe for 6.5 h or more with a median seeing of 1.2 arcsec.

Fig. 9

Locations of current and future astronomical facilities on the Mt. Hopkins Ridge. The MINERVA site and building are labeled at the bottom of the image.


The MINERVA site is located at (ϕ,λ)=31°4049.4N, 110°52′44.6′′W at an elevation of 7816 ft. Figure 9 shows the MINERVA telescope and building location in relation to the rest of the astronomical facilities on the Mt. Hopkins Ridge. The MINERVA telescopes will all be placed on the flat and smooth area approximately where the decommissioned FLWO 10-m gamma-ray telescope was located. The horizon limits have been measured to be between 15 and 20 deg over the full 360 deg azimuth range except for an 80-ft lightning rod northwest of the array, which covers a small total solid angle, but must be considered when tracking high declination sources across the sky, and the top corner of the MINERVA building, which reaches an elevation of 30deg north of the MINERVA telescope locations. The foundations for two Aqawans oriented with their long axes East-West are 2-ft-thick reinforced concrete that have been properly grounded. The telescope piers are 36 in. in diameter and extend 5ft below the ground surface. One of the piers is anchored into the underground remains of the 10-m pier.

MINERVA will use the building that was used for the FLWO 10-m gamma-ray telescope that is now decommissioned and whose function was replaced by VERITAS,100 located at the base of Mt. Hopkins. There are several rooms in the building, three of which are instrumental to the operation of MINERVA: (1) the spectrograph room is the southeastern-most room that is being converted into a class 100,000 clean room within which the KiwiSpec-MINERVA will be mounted; (2) the UPS room is located immediately northwest of the spectrograph room, where the power for the entire MINERVA facility will be routed through a facility grade UPS; and (3) the control room adjacent to the spectrograph room, which will house control computers and network equipment.


MINERVA Robotic Software

MINERVA will be a completely autonomous facility. MRS is being adapted from the Robo-AO software63 that has been successfully operating with a laser guide star adaptive optics system on the 60-in. telescope at Palomar Observatory for the past three years. Robo-AO has already completed the largest adaptive optics surveys to date with high observing efficiency and robust operation.101 The Robo-AO software was developed in a modular way such that it can be easily replicated and used for the robotic operation of laser adaptive optics imaging on other 1 to 2 m class telescopes. However, this design also allows a straightforward adaptation of the Robo-AO brain to the MINERVA system. The MRS architecture is shown in Fig. 10. The operation of MINERVA is carried out by six separate computers, main, telcom1, telcom2, telcom3, telcom4, and spec. All computers use Ubuntu 12.04.2 as the base operating system, and all source code is written in C++. Communications between the subsystems of MRS use the custom TCP/IP protocol developed for Robo-AO. This protocol is used to pass commands and exchange telemetry between each of the subsystems. This system will detect when one of the control daemons for a subsystem dies and will restart the subsystem automatically.

Fig. 10

MINERVA Robotic Software architecture. Blue boxes are subsystem control daemons, gray boxes are control or oversight daemons that control more than one daemon, and red boxes represent data file storage. The red dashed lines with arrows signify the paths for telemetry through the system, black lines denote the command paths, and the blue are the data paths (the three asterisks indicate three other data pathways from T1, T2, and T3 that we leave off for the sake of neatness). The green lines signify the pathways for the status information that is passed to the system monitoring software.


The main computer, main, coordinates operation of the facility while monitoring the performance and state of all system components as well as environmental conditions. The three control or oversight daemons running on this computer are the Robo System Control, the Watchdog, and the System Monitor. The Robo System Control provides the high-level control of MINERVA, coordinating the execution of all operations.

The System Monitor manages the information flow of the status of the subsystems to the entire robotic system. This part of the software regularly examines the status of each of the other software elements for their state of operation. It detects when one of the software subsystems has an error, crashes, or has other problems that might hinder the proper operation of the system. Issues are flagged and the operation of the automated system observations is stopped until the subsystem daemon can clear the issue. If the subsystem cannot correct the error, the automation system can take steps, up to and including restarting subsystems, in an attempt to continue operations. If it is unable to restart the system, a message will be sent for human assistance, and an attempt can be made to continue operation without the failed subsystem, e.g., if one telescope fails, the other three telescopes can continue with the primary science program. If operations cannot continue without the failed subsystem (e.g., the spectrograph), the software will close the Aqawans and shut everything down, leaving the system in a safe state.

The sole function of the Watchdog system will be to make sure the observing system remains in operation. If any system crashes or stops working properly, it will attempt to restart it, as well as stop telescope operations. If the process cannot be restarted and is not essential, the system will continue without the failed subsystem. If the failed subsystem is essential, the shutdown sequence will be initiated followed by a request for human help.

Each of the four telescopes and their suite of instruments are run with a single computer on which the station daemons are run, T1 through T4. Each station daemon reports telemetry information and handles the operation of the subsystems mounted on each telescope, including the telescope control system (TCS), guide cameras, imaging science camera, and filter wheel.

The spec computer is dedicated to the control of the spectrograph and runs the spectrograph daemon, which coordinates the operation of the spectrograph and monitoring of the spectrograph environment. This daemon also gathers data from the System Monitor and outputs the science data with extensive header information regarding the state of the facility during the time over which data were collected.

The software control of each hardware subsystem consists of a set of individual modules. Each interface module handles configuration file interactions, initialization, and error control. These modules are stacked together into larger modules, which are then managed by other facets of the robotic control system. The subsystems are run as daemons in the operating system; each separately manages the hardware under its control and runs a status monitor to sample subsystem performance and register errors that occur. Each of these subsystems is composed of many separate functions that initialize the hardware, monitor its function, and manage the operation of the hardware to achieve successful scientific output. In essence, each of the subsystem daemons are individual robotic programs that manage their hardware and operate according to external commands. The subsystem daemons communicate their state through the TCP/IP protocol to the System Monitor, and the System Monitor then parses and relays this information to the Watchdog and Robo System Control.

MRS will also employ a queue system modeled after the Robo-AO system built to read and organize program and target files in XML format.102 While the elimination and weighting criteria for the MINERVA program will be slightly different than for Robo-AO, the approach will be essentially the same for the primary MINERVA spectroscopic survey. MINERVA will operate primarily in a coordinated fashion. Targets will be observed by all telescopes simultaneously, vastly simplifying operations. On the occasion when one or more of the telescopes is needed for another program, e.g., transit photometry, a natural break in operations will be identified when the needed telescope(s) can be seamlessly left out of the next observing sequence. The telescope(s) will then be released from robotic control and can be operated remotely.

Development of a robotic system from scratch is an involved process, and can take years, especially for a system controlling multiple enclosures and telescopes. Using the Robo-AO software as a base, the MINERVA software development time has been cut by at least 50%. Most subsystems are under computer control and debugged; work continues on the TCS interface, and detailed spectrograph software development is awaiting delivery of the spectrograph hardware. Software to control the individual stations and the overall robotic software control system (including the queue scheduler) is currently undergoing final development, with robotic operations, under queue control, expected at the test facility in early 2015.



The commissioning of all four PlaneWave CDK-700 telescopes was performed from either the Caltech commissioning site (34°08′10.0′′N, 118°07′34.5′′W; elevation 800ft) or from the PlaneWave warehouse in Rancho Dominguez, California (33°52′14.1′′N, 118°14′49.4′′W; elevation 100ft). The basic functionality of the telescopes including software control via the PlaneWave Interface software (PWI) has been validated through frequent use beginning April 12, 2013, when the first telescope was delivered to the Caltech campus, and extending through the writing of this publication. In the following sections, we present the procedures and results for tests of telescope throughput and vignetting as well as fiber throughput and guiding.


Telescope Throughput

The equation for the photoelectron detection rate is


where S is the rate of detection of photoelectrons in astronomical data unit (ADU)/s for a source with apparent magnitude m relative to a zeropoint flux scale of Fλ(0) at an airmass X. The quantity Aeff is the effective aperture, 2998.33cm2, for the CDK-700 corresponding to 22% central obscuration by the secondary, f(λ) is the filter transmission function, η(λ) is the camera quantum efficiency, and g is the camera gain in units of e/ADU. The α factor is directly proportional to the atmospheric extinction with a proportionality constant of 0.4ln(10), and we use X=sec(z) to estimate airmass, which is sufficiently accurate for the range of zenith angles under consideration. We explicitly show the wavelength dependence of all quantities, and h and c are Planck’s constant and the speed of light, respectively. The efficiency term, τ(λ), characterizes all sources of attenuation not accounted for by the other terms in Eq. (1), such as the reflectivity of the mirrors and the transmission of the lenses.

Camera gain was directly measured from flat field images, and the theoretical transmission and reflectance of our optical elements were supplied by the manufacturers and are shown in Fig. 11. Multiplying these curves together (including the Nasmyth lens curve four times to account for the four surfaces on two lenses) and then integrating over the the V filter band gives a theoretical upper limit to the telescope throughput of τV(X=0)=84%.

Fig. 11

Transmission and reflectance of our optical elements as a function of wavelength. The reflectance of the tertiary mirror was measured at a 45 deg incident angle, and the Nasmyth lens transmission is for a single lens surface. The quantum efficiency of the Andor camera (BV chip) is shown for reference, as is the Johnson V bandpass, which roughly corresponds to the spectral range of KiwiSpec-MINERVA.


To calibrate the throughput of our telescopes, we observed bright (V<12) standard stars.103,104 The standard stars were observed in sequences of between 3 and 10 images per pointing with integration times adjusted to give high counts within the linear regime of the CCD. The signal-to-noise achieved at each observation was typically a few hundred. We then cycled through a list of standard stars for a given night between three and five times producing measurements over a range of airmasses and at different azimuthal angles.

The total counts measured within the integration time of each observation were converted into a V-band averaged throughput for a calculated airmass by rearranging Eq. (1):



We were able to collect data over airmasses between 1 and 2 with limited horizons at both the Caltech and Rancho Dominguez sites. Observations were spread as widely in the azimuthal angle as possible on a given night, and we repeated observations several times to average over a large variance in the individual measurements that we attribute to a highly variable aerosol content in the Los Angeles basin atmosphere. Figure 12 shows the cumulative results of our throughput measurements using different telescopes and different cameras over the course of five separate nights. The overall best fit for the mean telescope throughput is 70% with a high formal error of 30% due to the covariance between the atmospheric extinction—which is poorly constrained by our data—and the throughput scale. However, the best fit value for the V band atmospheric extinction, 0.27mags/airmass, is close to what we would expect and lends credence to the derived throughput of 70%. Fits from individual nights and sources using different telescopes and cameras also agree with this value within errors.

Fig. 12

Summary and fit of our throughput observations. The total throughput of the CDK-700 optics is estimated to be 70% in agreement with expectations.




To characterize the level of vignetting along our optical path, we perform astronomical observations and compare with the optical model of the telescope. The optical model of the telescope is summarized in Fig. 13, where the vignetting percentage with and without mirror baffles and RMS spot size is plotted as a function of the distance from the optical axis in the focal plane. For the f/6.5 optics of CDK-700, this corresponds to 45.3arcsec/mm.

Fig. 13

Optical model of the CDK-700 provided by PlaneWave Inc. showing the expected level of vignetting in percent as a function of off-axis distance (blue dots and line) and the RMS spot size in microns (gray dots and line). The on-sky measurement of Section 6.2 is shown in gold.


The vignetting experiments used a bright standard star at high elevation, SA 111 773 (V=8.97), observed alternately at the center of the CCD image and near to the four corners of the chip, 15.5 arcmin off the field center. We use standard aperture photometry to derive the flux of the star in each of the positions and then fit a polynomial function to the measured flux values of the standard star at the center of the chip to account for varying atmospheric conditions over the course of the observations. These variations were at the 1% level. The flux values normalized to the polynomial fit reveal the relative flux decrement observed with the star at the corners of the CCD. The results from our observations using the fourth MINERVA telescope performed on September 9, 2014, are shown in Fig. 14. The average vignetting measured at the off-axis positions is 4.9%±0.3%. This value is 0.8% above the expected level, which is a statistically significant discrepancy given our measurement accuracy. However, this level of vignetting at the edge of our photometric field can be calibrated straightforwardly and is not expected to adversely affect our science goals.

Fig. 14

Results from our on-sky vignetting test. A bright standard star was observed on September 9, 2014 alternately between the field center and the four corners of the CCD chip with telescope baffles in place. After correcting for zeropoint drift due to atmospheric changes, we derive a 4.9%±0.3% relative vignetting at an off-axis distance of 15.5 arcmin.



Pointing, Guiding, and Source Acquisition

The pointing and guiding performance of CDK-700 is dependent on a pointing model that converts astronomical coordinates into altitude and azimuth positions. The MINERVA commissioning site has limited sky access due to buildings and foliage preventing a full sky pointing model. Despite this limitation, we have found that the pointing and guiding of our telescopes typically achieve or exceed the specifications of Table 1. This level of performance surpasses our requirements to place our RV target stars within the 2 arcmin field of view of our active guiding cameras.

The median seeing at the Mt. Hopkins Ridge site is about a factor of two better than our fiber diameters (1.2 arcsec versus 2.3 arcsec), ensuring that minimal flux will be lost at the wings of the seeing disk given adequate guiding. The CDK-700 telescope open-loop tracking is accurate to a couple of arc seconds over the typical integration times of our primary program. Therefore, we have implemented a modified positional proportional-integral-derivative-type controller to correct for drifts and other inaccuracies. The controller input is the star position on the camera, and the output is the altitude and azimuth offsets to the telescope mount. It typically sends corrections once every few seconds. The optical design of the telescope allows us to perform a one-time calibration of the camera field rotation if the telescope derotator is turned off, which is ideal for high-cadence observing.

Once a target star arrives within the field of view of our guide cameras, the controller actively guides the telescope such that the star is placed on the calibrated pixel location of the fiber tip. The fiber tip location can be calibrated as frequently as needed to minimize losses due to offsets between the guiding center and the true center of the fiber. We have tested the stability of the fiber position on the guide camera and have found no measurable drifts or systematic offsets on day timescales. We have not measured the temporal stability over timescales from days to weeks. However, we do not anticipate that frequent calibrations will be needed as the optics are rigid and the fiber is strain relieved. Also, for these bright stars, reflections off the fiber cladding during observations can be detected and used as a secondary check for optical alignment.

Currently, the controller operates near the optimal level, showing an RMS pointing precision of 0.2arcsec, dominated by uncertainties in the measurement from seeing variations. Simulations of the required pointing accuracy indicate that coupling penalties <5% are incurred for a pointing accuracy of 0.2 arcsec RMS at any seeing from 0.5 arcsec to 2.5 arcsec. It is thus unlikely that the control system is contributing to any major loss of throughput in the system. The most convincing evidence that the controller is not adding significant noise to the pointing system is provided by the amplitude spectral density of the pointing errors; the error is comparable or lower at all sensed frequencies when the telescope is guiding (see Fig. 2 of Ref. 72). Typical guide camera exposure times are 0.1 s for stars from 4 to 6 mag, allowing sufficient sensitivity to successfully guide on the dimmest targets in our target list.

Results from a guiding test are shown in Fig. 15. For this test, first presented in Ref. 72, the location of a bright star was tracked on the FAU guide camera up until 49s, after which active guiding was initiated. The time to acquire the source on the chosen pixel (here x,y=160,150) was 20s. Although this controller has not been optimized to minimize the acquisition time, we use this result as the basis for total source acquisition time including telescope slew time.

Fig. 15

Results from the FAU guiding test,72 here showing a zoom-in of the source acquistion. The active guiding system was initiated around 49 s, and the position of the star stabilized on the target pixel in 20s.



Fiber Coupling and Throughput

The theoretical throughput curve for our 50-μm octagonal fibers has been calculated as a function of astronomical seeing based on the vendor supplied transmission specifications for the pellicle, the fiber transmission, and the expected transmission calculated for the input and output reflectance (see Fig. 3 of Ref. 72). On-sky tests of the fiber coupling device performed from the MINERVA commissioning site on the Caltech campus in 2arcsec seeing conditions show very good agreement with these expected values in consideration of fiber losses, reflection losses, and coupling efficiencies.

On-sky throughput observations performed on the Caltech commissioning site have confirmed the expected performance of the FAU design achieving 50% measured efficiency (45% throughput). Further tests will be needed at FLWO to validate the performance at more optimal conditions. The final version of this instrument will be deployed later this year and will incorporate minor modifications, such as a customized pellicle for slightly higher transmission (98 versus 92%). These results suggest that the throughput of our fiber system will be roughly 70% at FLWO where the median seeing is 1.2 arcsec. With these results and those of Sec. 6.1, we expect to lose 50% of the astronomical light from our telescope and fiber systems excluding losses from the atmosphere.


First Science Results

While the fair weather of the Los Angeles basin allows for routine commissioning operations, the Caltech site is a challenging place from which to obtain science grade astronomical data. Despite difficulties involving a highly variable atmosphere, significant obstruction, and copious stray light, we have been able to surpass our lower limit for photometric precision required by our secondary science objectives from this location. We also present new observations of WASP-52b, a transiting hot-Jupiter, as an end-to-end test of the MINERVA photometry pipeline.


High-Precision Photometry of 16 Cygnus

The secondary science goal of MINERVA is to search the transit windows of known and newly discovered super-Earths detected by the RV technique, including potential detections from the MINERVA target list. The transit of a 3 R planet around a Sun-like star (0.8M/M1.2) produces a decrement of light on the order of 1 mmag. This is an approximate upper bound of what would be considered a super-Earth and, therefore, represents a lower limit to the precision that must be achieved with MINERVA for our secondary science program to be viable. The timescale for this precision is also important. A super-Earth in its respective habitable zone of a Sun-like star will transit with a duration of 13h and have an ingress/egress time of 20min. Planets closer to their host star are easier to detect, have higher transit probabilities, and shorter transit durations. Therefore, this level of precision should be attained on timescales <20min.

On the evening of UT September 18, 2014, we observed one of the η targets that will be part of our RV survey, 16 Cyg AB (V=5.95/6.20, also HD 186,408/186,427). The first MINERVA telescope was equipped with an Andor iKON-L camera, an SBIG ST-i guide camera, and a seven-slot filter wheel. The telescope was controlled through the PWI interface, while the camera and active guiding were controlled through Maxim DL.105 A series of flats were taken during twilight. There were some clouds in the East, but the area of sky where we were observing remained clear throughout our observations.

Once on our target field, we aggressively defocused the camera while keeping our guide camera focused on a suitable guide star 20arcmin off-axis. The defocusing allows for the ample flux from our target to be spread over many more pixels mitigating the photometric error contributions from pixel-to-pixel variations and Poisson noise from the target itself. A second and important benefit from defocusing is to allow for long enough integration times to preclude shutter effects. To this second end, we also chose to observe in z band where the quantum efficiency of our camera is a factor of 2 below peak. A series of bias and dark frames were taken with the same setup following our observations.

Using 7 s integrations and full frame readouts, we obtained 208 images of the 16 Cyg field spanning almost exactly an hour from UT 0352 to UT 0453. We actively guided throughout the course of these observations, but we did not achieve optimal guiding results. A drift of a few pixels was seen over the course of these observations, and there was one episode of a fairly large (7pixels or 4 arcsec) guiding excursion that took place over 1.5min. The frames were bias and dark subtracted, and divided by our calibrated median twilight flats to correct for pixel-to-pixel variations.

To extract the photometry from the calibrated science images, we used the multiaperture mode of AstroImageJ,106 which uses simple aperture photometry and sky-background subtraction. For all of the stars for which we measured a lightcurve, we used a constant aperture size of 30 pixels (18.4 arcsec) and a sky annulus with an inner radius of 90 pixels and an outer radius of 100 pixels. The rather large sky annulus was necessary so that the background annuli centered on each of the stars in 16 Cyg did not include the other member of the binary. On each science image, we recentered the apertures on the stellar centroids using the center-of-light method.107

We used a set of five nearby comparison stars to remove systematics in the light curves of 16 Cyg A and B. The set of comparison stars included the corresponding other member of the 16 Cyg binary, which effectively provided all of the comparison information, as the next brightest comparison star was 3.5 magnitudes fainter. To remove any lingering systematic trends in the data, we then performed a linear detrending against airmass.

The detrended photometric time series of both 16 Cyg A and B are shown in Fig. 16, where the individual points have been binned into 1-min intervals. The Allan Variance of 16 Cyg A is plotted in Fig. 17. The RMS of the unbinned photometry is 2.7 mmag, while we achieved sub-mmag precision on 3- to 5-min timescales. The stability of the atmosphere on Mt. Hopkins is considerably better than in Pasadena. Therefore, these first results from our commissioning site support the prospect of routinely achieving sub-mmag photometric precision from FLWO.

Fig. 16

Detrended photometric time series of the 16 Cyg A and B observations performed on UT September 18, 2014 from the MINERVA commissioning site in Pasadena, California. The individual 7-s exposures have been binned into 1-min intervals.


Fig. 17

Allan Variance plot of the photometric time series of 16 Cygnus A performed on UT September 18, 2014 from the MINERVA commissioning site in Pasadena, California. On the 17-s duty cycle of our observations, a 2.7-mmag precision was achieved that bins down to <1mmag on 3- to 5-min timescales.



WASP-52b: New Transit Observations and Modeling

WASP-52b is an inflated hot-Jupiter with M=0.5MJ and R=1.3RJ in a slightly misaligned, 1.75 d orbit.108 The transits of WASP-52b were first observed by the SuperWASP survey76 in 2008 and 2009, and the most recent observations in the literature are precision light curves obtained in September 2011.108 The host star is reported to have a mass of 0.87 M and a rotation period of 11.8 days, suggesting a gyrochronological age of 0.4 Gyr. The age, small orbital distance, and obliquity of this planet may have important implications for the mechanism by which it formed and for the formation of hot-Jupiter systems in general.109

We observed WASP-52 on the evening of UT September 18, 2014, from the MINERVA commissioning site on the Caltech campus in Pasadena, California. It has a V magnitude of 12, >200 times fainter than 16 Cygnus. To maximize the signal-to-noise ratio of our observations of WASP-52, we again aggressively defocused. The observations were performed in an r band, and we used an integration time of 120 s to build up high signal-to-noise on the defocused star image. Active offset guiding was used throughout the observations, but the guide star was faint and the target star drifted by 15pixels or 9 arcsec. This drift was little more than half the size of our defocused star images. The target was tracked for 4 h 19 min starting 13min. before the start of ingress.

We used the same AstroImageJ106 reduction pipeline as in our 16 Cyg observations. We bias and dark subtracted our raw science images, before using a median twilight flat to remove image inhomogeneities. We then conducted simply aperture photometry with sky background subtraction on our calibrated images and extracted light curves for WASP-52 and 12 other nearby comparison stars. For all of the stars for which we extracted photometry, we used a fixed 20-pixel aperture radius (12.3 arcsec) and a sky annulus with an inner radius of 30 pixels and an outer radius of 50 pixels. We recentered the apertures on the individual stellar centroids in each calibrated image using the center-of-light method.107

Figure 18 shows the calibrated and detrended photometry of WASP-52. We achieved 3 mmag precision on the 131-s duty cycle of these observations, which binned down to 1mmag on 30-min timescales. The Allan Variance for this photometric time series can be seen in Fig. 19.

Fig. 18

Normalized and detrended light curve data for WASP-52 with best fit transit model overlaid. The transit center, TC, is observed to be 6918.79085 BJDTDB.


Fig. 19

Allan Variance plot of the photometric time series of WASP-52 performed on UT September 18, 2014 from the MINERVA commissioning site in Pasadena, California. On the 131-s duty cycle of our observations, a 2.7-mmag precision was achieved that bins down to 1 mmag on 30-min timescales.


The midpoint Julian Date in Coordinated Universal Time (JDUTC) of each integration was recorded for each observation, which we convert to BJDTDB.110 The light curve of WASP-52 was then fit using EXOFAST.111 We employ priors on the planet period, P; the stellar age, Age; metallicity, [Fe/H]; and effective temperature, Teff; all taken from Ref. 108. The transit parameters are tied to the Yonsei-Yale stellar models112 through the stellar density as determined by the scaled semimajor axis, a/R. This, in turn, informs a prior for a quadratic stellar limb darkening model parameterized by u1 and u2.113 The remaining free parameters of the fit are the baseline flux, F0; transit time, TC; cosi; RP/R; logM; an additional noise term added in quadrature to the transit errors, σr2; and an error scaling for the uncertainties in the photometric time series, TranScale. While the differential photometry technique accounts for the majority of the airmass effects on our photometry measurements, additional drift is noticed that we suspect is due to scattered light effects. We, therefore, include an additional nuisance parameter in the fit that accounts for a linear drift with airmass.

The differential evolution Markov Chain Monte Carlo (MCMC) sampler required 2080 burn-in steps and 13,074 steps after the burn-in for adequate mixing.111 Figure 18 shows the best fit transit model overlaid on the detrended light curve data. The median parameter values and 1σ errors are reported in Table 2.

Table 2

Median values and 68% confidence interval for WASP-52b.

Stellar parameters:
M⋆Mass (0.842−0.045+0.041
R⋆Radius (0.732−0.042+0.040
L⋆Luminosity (0.305−0.051+0.056
ρ⋆Density (cgs)3.03−0.33+0.39
log(g⋆)Surface gravity (cgs)4.634−0.026+0.028
AgeAge (Gyr)0.45−0.24+0.27
TeffEffective temperature (K)5020±100
Planetary parameters:
PPeriod (days)1.7497798±0.0000012
aSemimajor axis (AU)0.02682−0.00049+0.00043
RPRadius (1.166−0.088+0.084
TeqEquilibrium temperature (K)1264−46+45
〈F〉Incident flux (0.579−0.080+0.086
Radial velocity parameters:
TATime of ascending node2456918.35341±0.00036
TDTime of descending node2456919.22830±0.00036
Primary transit parameters:
TCTime of transit (2456918.79085±0.00036
RP/R⋆Radius of planet in stellar radii0.1637−0.0039+0.0035
a/R⋆Semimajor axis in stellar radii7.88−0.29+0.33
u1Linear limb-darkening coeff0.560±0.056
u2Quadratic limb-darkening coeff0.174−0.054+0.053
iInclination (degrees)86.52−0.63+0.80
bImpact parameter0.479−0.095+0.066
δTransit depth0.0268−0.0013+0.0012
TFWHMFWHM duration (days)0.0620−0.0013+0.0014
τIngress/egress duration (days)0.0134±0.0014
T14Total duration (days)0.0754−0.0012+0.0013
PTA priori0.1061−0.0039+0.0038
PT,GA priori0.1477−0.0063+0.0060
DepthFlux decrement at mid transit0.03155±0.00065
d/R⋆Separation at mid transit7.88−0.29+0.33
σr2Variance of transit red noise−0.0000034−0.0000100+0.0000083
TranScaleScaling of transit errors1.14−0.66+0.49
σrTransit red noise0.0000−0.00+0.0022
F0Baseline flux1.00003−0.00030+0.00029
Secondary eclipse parameters:
TSTime of eclipse (2456919.66574±0.00036

Our parameter distributions are in broad agreement with those of Ref. 108. However, 1 to 1.5σ discrepancies are seen between the scaled semimajor axis of the orbit, a/R, and associated parameters. This discrepancy may arise from the difference in parametrization of the stellar limb darkening (Hebrard et al. use a four-parameter limb darkening model, we use two) or from the MCMC chains stepping in different parameters, which introduces different priors.

The mid-transit time we derive for WASP-52b is TC=23±31s advanced from the predicted mid-transit time, 643 transits from the literature ephemeris. This is consistent with the published ephemeris and extends the time baseline of WASP-52b transit photometry to approximately six years. The possibility of a third companion was mentioned in Ref. 108 based on a potential acceleration of 30ms1year1 seen in the RV data. If real, this could be due to a massive giant planet at a few to several astronomical unit (AU) or a stellar companion with mass <0.8 M further out.114 At those large orbital distances, the light travel time effect would dominate a transit timing variation (TTV) signal.115 Using the lack of TTVs in our data as a constraint, we can rule out the existence of brown dwarfs, M10 to 15MJup, out to 6AU.

We also examine the photometry data from Ref. 108 and find no significant variation of the transit depths in the r band for three epochs with full transit coverage. Independent fits to two epochs of data taken with the EulerCam on the 1.2-m Euler-Swiss telescope in La Silla, Chile, and one epoch of data from the 0.94-m James Gregory Telescope in St. Andrews, Scotland, agree with our transit depth to within 1σ. This suggests that WASP-52b may be a suitable target for follow-up observations of occultation events for the purpose of atmospheric studies.116


Concluding Remarks

The statistics from both RV and transit surveys of exoplanets have informed us that planets are very common throughout the Galaxy and that the most common type of planet may be of a variety unrepresented in the Solar System—super-Earths—with masses between that of the Earth and Neptune. The prevalence of these planets inferred from survey data imply that the stars in our local Solar Neighborhood should harbor many of these, with some fraction orbiting in their respective habitable zones. The RV technique is the most promising method to date for detecting these planets. However, the detection of such small planets via the reflex motion of their host stars presents a significant challenge that requires cutting-edge instrumentation and a significant amount of observing time as the planet signals lie at or below the level of systematic noise generated from the stellar surface.

MINERVA is an innovative facility designed to address these demands in a modular, cost-effective manner. By employing four small-aperture telescopes from a commercial vendor, we obtain a 1.4 m effective aperture for a fraction of the cost of a single telescope. The small étendue translates to a smaller spectrograph that is easier to stabilize and requires less infrastructure. Once the facility is complete, the array of telescopes will undertake an automated survey of a fixed target list over the course of several years that is expected to result in exoplanet detections of high scientific value.

We have presented the design and major components of MINERVA herein, and the on-sky performance of our equipment up to and including the feeding of starlight into our fibers has been validated. Meanwhile, the procedures and data reduction pipeline for our secondary science goal of detecting transit events around nearby stars have been demonstrated with observations from our test facility. There are, however, several more milestones the project will need to reach before science operations commence.


Future Prospects: MINERVA-South and MINERVA-Red

In addition to the primary survey presented in this article, the modularity of the MINERVA design offers several opportunities for expansion that are already being pursued. The abundant yield of transiting exoplanets around bright stars in both the Northern and Southern hemispheres expected from the Kepler K2 Ecliptic Mission79 and the upcoming TESS80 will require intensive RV follow-up. Planets with periods up to 100d will be detected within the continuous viewing zone of TESS near the ecliptic poles that will require long-term RV monitoring. As has been made clear from the Kepler prime mission, the most interesting planets will require a substantial amount of dedicated telescope time.117118.119 Hence, there are significant opportunities for a MINERVA-like facility in the Southern hemisphere. MINERVA-South will take advantage of the same economies as described in Sec. 2, with the added advantage that our team’s investment in software/hardware development and operational expertise will easily translate to the Southern facility. We expect MINERVA-South to be of essentially the same design as MINERVA, sited in Australia or Chile, and operational by 2018 to capitalize on the coming flood of planet candidates from K2 and TESS.

In the near term, we are expanding the reach of the MINERVA project to include nearby M stars with a second instrument specifically designed for gathering precise RVs of the closest low-mass stars to the Sun. Recent results from Kepler and ground based surveys indicate that compact systems of planets in orbit around low-mass stars may be extremely common.35,49,50 Statistically speaking, we expect some of the closest stars to the Sun to host systems of planets. However, these small, cool stars are often too faint to observe at the optical wavelengths where most precision RV instruments operate, including KiwiSpec-MINERVA. MINERVA-Red, a parallel effort to the main MINERVA survey, will specifically target a small sample of nearby mid- to late-M stars.

The MINERVA-Red instrument is a fiber-fed echelle spectrograph housed in a vacuum chamber and optimized to cover the 800 to 900 nm spectral region at a resolution of R55,000. The instrument is designed around a single-mode fiber input, which allows the instrument to be very compact and stable, and also eliminates modal noise as a source of RV error. While the single-mode fiber reduces the possible coupling efficiency of starlight into the fiber, special attention has been paid to maximizing the optical throughput of the rest of the instrument. To that end, the spectrograph will operate with a pair of PlaneWave CDK-700 telescopes that have gold mirror coatings, significantly enhancing the telescope efficiency in the wavelength range of interest. These will be in addition to the four CDK-700s that will be used for the primary survey. The instrument optics are all optimized for this relatively narrow spectral range, and a deep depletion detector will ensure high quantum efficiency and low fringing. The MINERVA-Red instrument is currently under construction and is slated to be deployed with the first of its two telescopes at the Mt. Hopkins site by mid-year 2015.


This work was partially supported by funding from the Center for Exoplanets and Habitable Worlds. The Center for Exoplanets and Habitable Worlds is supported by the Pennsylvania State University, the Eberly College of Science, and the Pennsylvania Space Grant Consortium. MINERVA hardware has been partially funded by the Australian Research Council’s Linkage, Infrastructure, Equipment and Facilities funding scheme (project LE140100050). We are grateful to the Mt. Cuba Astronomical Foundation and the David and Lucile Packard Foundation for their generous funding of MINERVA hardware and personnel. C.H.B. is supported by a NASA Nancy Grace Roman Technology Fellowship.



B. Campbell, G. A. H. Walker and S. Yang, “A search for substellar companions to solar-type stars,” Astrophys. J. 331, 902–921 (1988).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/166608Google Scholar


D. W. Latham et al., “The unseen companion of HD114762: a probable brown dwarf,” Nature 339, 38–40 (1989).NATUAS0028-0836http://dx.doi.org/10.1038/339038a0Google Scholar


A. Wolszczan and D. A. Frail, “A planetary system around the millisecond pulsar PSR1257 + 12,” Nature 355, 145–147 (1992).NATUAS0028-0836http://dx.doi.org/10.1038/355145a0Google Scholar


A. Wolszczan, “Confirmation of Earth-mass planets orbiting the millisecond pulsar PSR B1257 + 12,” Science 264, 538–542 (1994).SCIEAS0036-8075http://dx.doi.org/10.1126/science.264.5158.538Google Scholar


M. Mayor and D. Queloz, “A Jupiter-mass companion to a solar-type star,” Nature 378, 355–359 (1995).NATUAS0028-0836http://dx.doi.org/10.1038/378355a0Google Scholar


G. W. Marcy and R. P. Butler, “A planetary companion to 70 Virginis,” Astrophys. J. Lett. 464, L147 (1996).AJLEEY0004-637Xhttp://dx.doi.org/10.1086/310096Google Scholar


J. T. Wright et al., “The Exoplanet Orbit Database,” Publ. Astron. Soc. Pac. 123, 412–422 (2011).PASPAU0004-6280http://dx.doi.org/10.1086/659427Google Scholar


E. Han et al., “Exoplanet Orbit Database. II. Updates to Exoplanets.org,” Publ. Astron. Soc. Pac. 126, 827–837 (2014).PASPAU0004-6280http://dx.doi.org/10.1086/678447Google Scholar


R. L. Akeson et al., “The NASA Exoplanet Archive: data and tools for exoplanet research,” PASP 125, 989–999 (2013).PASPAU0004-6280http://dx.doi.org/10.1086/672273Google Scholar



F. Rosenblatt, “A two-color photometric method for detection of extra solar planetary systems,” Icarus 14, 71–93 (1971).ICRSA50019-1035http://dx.doi.org/10.1016/0019-1035(71)90103-5Google Scholar


W. J. Borucki and A. L. Summers, “The photometric method of detecting other planetary systems,” Icarus 58, 121–134 (1984).ICRSA50019-1035http://dx.doi.org/10.1016/0019-1035(84)90102-7Google Scholar


W. J. Borucki, J. D. Scargle and H. S. Hudson, “Detectability of extrasolar planetary transits,” ApJ 291, 852–854 (1985).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/163121Google Scholar


G. W. Henry et al., “A transiting ‘51 Peg-like’ planet,” Astrophys. J. Lett. 529, L41–L44 (2000).AJLEEY0004-637Xhttp://dx.doi.org/10.1086/312458Google Scholar


D. Charbonneau et al., “Detection of planetary transits across a Sun-like star,” Astrophys. J. Lett. 529, L45–L48 (2000).AJLEEY0004-637Xhttp://dx.doi.org/10.1086/312457Google Scholar


S. Mao and B. Paczynski, “Gravitational microlensing by double stars and planetary systems,” ApJ Lett. 374, L37–L40 (1991).AJLEEY0004-637Xhttp://dx.doi.org/10.1086/186066Google Scholar


A. Gould and A. Loeb, “Discovering planetary systems through gravitational microlenses,” ApJ 396, 104–114 (1992).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/171700Google Scholar


I. A. Bond et al., “OGLE 2003-BLG-235/MOA 2003-BLG-53: a planetary microlensing event,” Astrophys. J. Lett. 606, L155–L158 (2004).AJLEEY0004-637Xhttp://dx.doi.org/10.1086/420928Google Scholar


C. Marois et al., “Direct imaging of multiple planets orbiting the star HR 8799,” Science 322, 1348–1352 (2008).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1166585Google Scholar


P. Kalas et al., “Optical images of an exosolar planet 25 light-years from Earth,” Science 322, 1345–1348 (2008).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1166609Google Scholar


W. J. Borucki et al., “Kepler planet-detection mission: introduction and first results,” Science 327, 977–980 (2010).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1185402Google Scholar


N. C. Santos et al., “The HARPS survey for southern extra-solar planets. II. A 14 Earth-masses exoplanet around μ Arae,” A&A 426, L19–L23 (2004).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361:200400076Google Scholar


D. A. Fischer and J. Valenti, “The planet-metallicity correlation,” Astrophys. J. 622, 1102–1117 (2005).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/apj.2005.622.issue-2Google Scholar


J. A. Johnson et al., “Giant planet occurrence in the stellar mass-metallicity plane,” Publ. Astron. Soc. Pac. 122, 905–915 (2010).PASPAU0004-6280http://dx.doi.org/10.1086/653442Google Scholar


E. B. Ford, V. Lystad and F. A. Rasio, “Planet-planet scattering in the upsilon Andromedae system,” Nature 434, 873–876 (2005).NATUAS0028-0836http://dx.doi.org/10.1038/nature03427Google Scholar


A. C. M. Correia et al., “The HARPS search for southern extra-solar planets. XVI. HD 45364, a pair of planets in a 32 mean motion resonance,” Astron. Astrophys. 496, 521–526 (2009).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361:200810774Google Scholar


E. J. Rivera et al., “The Lick-Carnegie exoplanet survey: a Uranus-Mass fourth planet for GJ 876 in an extrasolar Laplace configuration,” Astrophys. J. 719, 890–899 (2010).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/719/1/890Google Scholar


E. J. Rivera et al., “A 7.5 M planet orbiting the nearby star, GJ 876,” Astrophys. J. 634, 625–640 (2005).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/apj.2005.634.issue-1Google Scholar


D. Valencia, R. J. O’Connell and D. Sasselov, “Internal structure of massive terrestrial planets,” Icarus 181, 545–554 (2006).ICRSA50019-1035http://dx.doi.org/10.1016/j.icarus.2005.11.021Google Scholar


D. Valencia, D. D. Sasselov and R. J. O’Connell, “Radius and structure models of the first super-Earth planet,” ApJ 656, 545–551 (2007).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/509279Google Scholar


S. Udry et al., “The HARPS search for southern extra-solar planets. XI. Super-Earths (5 and 8 M) in a 3-planet system,” A&A 469, L43–L47 (2007).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361:20077612Google Scholar


A. W. Howard et al., “The occurrence and mass distribution of close-in super-Earths, Neptunes, and Jupiters,” Science 330, 653–655 (2010).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1194854Google Scholar


C. Lovis et al., “Towards the characterization of the hot Neptune/super-Earth population around nearby bright stars,” in IAU Symp., Vol. 253, pp. 502–505, Cambridge University Press, New York (2009).Google Scholar


M. Mayor et al., “The HARPS search for southern extra-solar planets XXXIV. Occurrence, mass distribution and orbital properties of super-Earths and Neptune-mass planets,” arXiv:1109.2497, http://adsabs.harvard.edu/abs/2011arXiv1109.2497M (2011).Google Scholar


X. Bonfils et al., “The HARPS search for southern extra-solar planets. XXXI. The M-dwarf sample,” Astron. Astrophys. 549, A109 (2013).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361/201014704Google Scholar



W. J. Borucki et al., “Characteristics of Kepler planetary candidates based on the first data set,” Astrophys. J. 728, 117 (2011).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/728/2/117Google Scholar


W. J. Borucki et al., “Characteristics of planetary candidates observed by Kepler. II. Analysis of the first four months of data,” Astrophys. J. 736, 19 (2011).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/736/1/19Google Scholar


N. M. Batalha et al., “Planetary candidates observed by Kepler. III. Analysis of the first 16 months of data,” Astrophys. J. Suppl. 204, 24 (2013).APJSA20067-0049http://dx.doi.org/10.1088/0067-0049/204/2/24Google Scholar


C. J. Burke et al., “Planetary candidates observed by Kepler IV: planet sample from Q1–Q8 (22 months),” Astrophys. J. Suppl. 210, 19 (2014).APJSA20067-0049http://dx.doi.org/10.1088/0067-0049/210/2/19Google Scholar


P. S. Muirhead et al., “Characterizing the cool KOIs. III. KOI 961: a small star with large proper motion and three small planets,” Astrophys. J. 747, 144 (2012).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/747/2/144Google Scholar


T. Barclay et al., “A sub-Mercury-sized exoplanet,” Nature 494, 452–454 (2013).NATUAS0028-0836http://dx.doi.org/10.1038/nature11914Google Scholar


A. W. Howard et al., “The California Planet Survey. I. Four new giant exoplanets,” Astrophys. J. 721, 1467–1481 (2010).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/721/2/1467Google Scholar


A. Cassan et al., “One or more bound planets per Milky Way star from microlensing observations,” Nature 481, 167–169 (2012).NATUAS0028-0836http://dx.doi.org/10.1038/nature10684Google Scholar


F. Fressin et al., “The false positive rate of Kepler and the occurrence of planets,” Astrophys. J. 766, 81 (2013).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/766/2/81Google Scholar


T. D. Morton and J. Swift, “The radius distribution of planets around cool stars,” Astrophys. J. 791, 10 (2014).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/791/1/10Google Scholar


D. Foreman-Mackey, D. W. Hogg and T. D. Morton, “Exoplanet population inference and the abundance of Earth analogs from noisy, incomplete catalogs,” arXiv:1406.3020 (2014).Google Scholar


E. A. Petigura, A. W. Howard and G. W. Marcy, “Prevalence of Earth-size planets orbiting Sun-like stars,” Proc. Natl. Acad. Sci. 110, 19273–19278 (2013).PNASA60027-8424http://dx.doi.org/10.1073/pnas.1319909110Google Scholar


J. J. Swift et al., “Characterizing the cool KOIs. IV. Kepler-32 as a prototype for the formation of compact planetary systems throughout the galaxy,” Astrophys. J. 764, 105 (2013).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/764/1/105Google Scholar


C. D. Dressing and D. Charbonneau, “The occurrence rate of small planets around small stars,” Astrophys. J. 767, 95 (2013).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/767/1/95Google Scholar


J. T. Wright et al., “The frequency of hot Jupiters orbiting nearby solar-type stars,” Astrophys. J. 753, 160 (2012).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/753/2/160Google Scholar


X. Bonfils et al., “The HARPS search for southern extra-solar planets. X. A m sin i=11M planet around the nearby spotted M dwarf GJ 674,” Astron. Astrophys. 474, 293–299 (2007).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361:20077068Google Scholar


J. N. Winn et al., “A super-Earth transiting a naked-eye star,” Astrophys. J. Lett. 737, L18 (2011).AJLEEY0004-637Xhttp://dx.doi.org/10.1088/2041-8205/737/1/L18Google Scholar


X. Dumusque et al., “An Earth-mass planet orbiting α Centauri B,” Nature 491, 207–211 (2012).NATUAS0028-0836http://dx.doi.org/10.1038/nature11572Google Scholar


R. A. Wittenmyer et al., “GJ 832c: a super-Earth in the habitable zone,” Astrophys. J. 791, 114 (2014).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/791/2/114Google Scholar


M. Endl et al., “The Mt John University Observatory search for Earth-mass planets in the habitable zone of Alpha Centauri,” Int. J. Astrobio. 14, 305–312 (2014).http://dx.doi.org/10.1017/S1473550414000081Google Scholar


B. Sato et al., “The N2K Consortium. II. A transiting hot Saturn around HD 149026 with a large dense core,” Astrophys. J. 633, 465–473 (2005).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/apj.2005.633.issue-1Google Scholar


J. N. Winn et al., “A super-Earth transiting a naked-eye star,” Astrophys. J. Lett. 737, L18 (2011).AJLEEY0004-637Xhttp://dx.doi.org/10.1088/2041-8205/737/1/L18Google Scholar


J. F. Kasting, D. P. Whitmire and R. T. Reynolds, “Habitable zones around main sequence stars,” Icarus 101, 108–128 (1993).ICRSA50019-1035http://dx.doi.org/10.1006/icar.1993.1010Google Scholar


R. K. Kopparapu et al., “Habitable zones around main-sequence stars: new estimates,” Astrophys. J. 765, 131 (2013).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/765/2/131Google Scholar


G. Anglada-Escudé et al., “A planetary system around the nearby M Dwarf GJ 667C with at least one super-Earth in its habitable zone,” Astrophys. J. Lett. 751, L16 (2012).AJLEEY0004-637Xhttp://dx.doi.org/10.1088/2041-8205/751/1/L16Google Scholar


P. Robertson and S. Mahadevan, “Disentangling planets and stellar activity for Gliese 667C,” Astrophys. J. Lett. 793, L24 (2014).http://dx.doi.org/10.1088/2041-8205/793/2/L24Google Scholar


R. L. Riddle et al., “The Robo-AO software: fully autonomous operation of a laser guide star adaptive optics and science system,” Proc. SPIE 8447, 84472O (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.925475Google Scholar


M. Bottom et al., “Optimizing Doppler surveys for planet yield,” Publ. Astron. Soc. Pac. 125, 240–251 (2013).PASPAU0004-6280http://dx.doi.org/10.1086/670174Google Scholar


X. Dumusque et al., “Planetary detection limits taking into account stellar noise. I. Observational strategies to reduce stellar oscillation and granulation effects,” Astron. Astrophys. 525, A140 (2011).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361/201014097Google Scholar


A. Baranne et al., “ELODIE: a spectrograph for accurate radial velocity measurements,” Astron. Astrophys. Suppl. 119, 373–390 (1996).AAESB90365-0138http://dx.doi.org/10.1051/aas:1996251Google Scholar


A. W. Howard et al., “The NASA-UC Eta-Earth Program. I. A super-Earth orbiting HD 7924,” Astrophys. J. 696, 75–83 (2009).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/696/1/75Google Scholar


N. McCrady and C. Nava, “Minerva exoplanet detection sensitivity from simulated observations,” 2014, http://extras.aas.org/meeting/aas223/348/348.03 (21 March 2015).Google Scholar


G. Basri et al., “Photometric variability in Kepler target stars. II. An overview of amplitude, periodicity, and rotation in first quarter data,” Astron. J. 141, 20 (2011).ANJOAA0004-6256http://dx.doi.org/10.1088/0004-6256/141/1/20Google Scholar


A. McQuillan, S. Aigrain and T. Mazeh, “Measuring the rotation period distribution of field M dwarfs with Kepler,” Mon. Not. R. Astron. Soc. 432, 1203–1216 (2013).MNRAA40035-8711http://dx.doi.org/10.1093/mnras/stt536Google Scholar


A. W. Howard et al., “Planet occurrence within 0.25 AU of solar-type stars from Kepler,” Astrophys. J. Suppl. 201, 15 (2012).APJSA20067-0049http://dx.doi.org/10.1088/0067-0049/201/2/15Google Scholar


M. Bottom et al., “Design, motivation, and on-sky tests of an efficient fiber coupling unit for 1-meter class telescopes” Proc. SPIE 9147, 91472E (2014).http://dx.doi.org/10.1117/12.2055605Google Scholar


G. W. Marcy and R. P. Butler, “Precision radial velocities with an iodine absorption cell,” Publ. Astron. Soc. Pac. 104, 270–277 (1992).PASPAU0004-6280http://dx.doi.org/10.1086/132989Google Scholar


S. Tuttle et al., “Builders instead of consumers: training astronomers in instrumentation & observation,” http://adsabs.harvard.edu/abs/2014arXiv1410.2856T (2014).Google Scholar


G. Á. Bakos et al., “System description and first light curves of the Hungarian automated telescope, an autonomous observatory for variability search,” Publ. Astron. Soc. Pac. 114, 974–987 (2002).PASPAU0004-6280http://dx.doi.org/10.1086/pasp.2002.114.issue-799Google Scholar


D. L. Pollacco et al., “The WASP project and the SuperWASP cameras,” Publ. Astron. Soc. Pac. 118, 1407–1418 (2006).PASPAU0004-6280http://dx.doi.org/10.1086/pasp.2006.118.issue-848Google Scholar


P. J. Wheatley et al., “The Next Generation Transit Survey (NGTS),” 2013, http://www.epj-conferences.org/articles/epjconf/abs/2013/08/epjconf_hpcs2012_13002/epjconf_hpcs2012_13002.html (21 March 2015).Google Scholar


D. Dragomir et al., “TERMS photometry of known transiting exoplanets,” Astron. J. 142, 115 (2011).ANJOAA0004-6256http://dx.doi.org/10.1088/0004-6256/142/4/115Google Scholar


S. B. Howell et al., “The K2 mission: characterization and early results,” Publ. Astron. Soc. Pac. 126, 398–408 (2014).PASPAU0004-6280http://dx.doi.org/10.1086/676406Google Scholar


G. R. Ricker et al., “Transiting exoplanet survey satellite (TESS),” Proc. SPIE 9143, 914320 (2014).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.2063489Google Scholar



R. L. Hedrick et al., “New paradigms for producing high-performing meter class ground-based telescopes,” Proc. SPIE 7739, 77391Q (2010).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.857496Google Scholar


T. M. Brown et al., “Las Cumbres Observatory global telescope network,” Publ. Astron. Soc. Pac. 125, 1031–1055 (2013).PASPAU0004-6280http://dx.doi.org/10.1086/673168Google Scholar


B. Chazelas et al., “New scramblers for precision radial velocity: square and octagonal fibers,” Proc. SPIE 7739, 773947 (2010).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.856874Google Scholar


F. Bouchy et al., “SOPHIE+: First results of an octagonal-section fiber for high-precision radial velocity measurements,” Astron. Astrophys. 549, A49 (2013).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361/201219979Google Scholar


J. Baudrand and G. A. H. Walker, “Modal noise in high-resolution, fiber-fed spectra: a study and simple cure,” Publ. Astron. Soc. Pac. 113, 851–858 (2001).PASPAU0004-6280http://dx.doi.org/10.1086/pasp.2001.113.issue-785Google Scholar


S. I. Barnes and P. J. MacQueen, “A high-efficiency fibre double-scrambler prototype,” Proc. SPIE 7735, 773567 (2010).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.857945Google Scholar


P. P. Plavchan et al., “Precision near-infrared radial velocity instrumentation II: noncircular core fiber scrambler,” Proc. SPIE 8864, 88640G (2013).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.2023696Google Scholar


S. I. Barnes et al., “KiwiSpec—an advanced spectrograph for high resolution spectroscopy: optical design and variations,” Proc. SPIE 8446, 844688 (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.926527Google Scholar


S. Gibson et al., “KiwiSpec—an advanced spectrograph for high resolution spectroscopy: prototype design and performance,” Proc. SPIE 8446, 844648 (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.925501Google Scholar


M. Mayor et al., “Setting new standards with HARPS,” Messenger 114, 20–24 (2003).Google Scholar


R. Cosentino et al., “Harps-N: the new planet hunter at TNG,” Proc. SPIE 8446, 84461V (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.925738Google Scholar


J. D. Crane et al., “The Carnegie Planet Finder Spectrograph: integration and commissioning,” Proc. SPIE 7735, 773553 (2010).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.857792Google Scholar


R. P. Butler et al., “Attaining Doppler precision of 3 M s-1,” Publ. Astron. Soc. Pac. 108, 500 (1996).PASPAU0004-6280http://dx.doi.org/10.1086/133755Google Scholar


S. I. Barnes, “EchMod: a MATLAB toolbox for modeling high resolution échelle spectrographs,” Proc. SPIE 8550, 85501T (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.981563Google Scholar




J. Irwin et al., “The MEarth project: searching for transiting habitable super-Earths around nearby M dwarfs,” in IAU Symp., Vol. 253, pp. 37–43, Cambridge University Press, New York (2009).Google Scholar


J. D. Gibson, G. G. Williams and T. Trebisky, “Seeing trends from deployable Shack-Hartmann wavefront sensors, MMT Observatory, Arizona, USA,” Proc. SPIE 8444, 844432 (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.924897Google Scholar


Very Energetic Radiation Imaging Telescope Array System (VERITAS), http://veritas.sao.arizona.edu/.Google Scholar


C. Baranec et al., “High-efficiency autonomous laser adaptive optics,” Astrophys. J. Lett. 790, L8 (2014).AJLEEY0004-637Xhttp://dx.doi.org/10.1088/2041-8205/790/1/L8Google Scholar


R. L. Riddle et al., “The Robo-AO automated intelligent queue system,” Proc. SPIE 9152, 91521E (2014).http://dx.doi.org/10.1117/12.2056534Google Scholar


A. U. Landolt, “UBVRI photometric standard stars in the magnitude range 11.5–16.0 around the celestial equator,” Astron. J. 104, 340–371 (1992).ANJOAA0004-6256http://dx.doi.org/10.1086/116242Google Scholar


A. U. Landolt, “UBVRI photometric standard stars around the sky at +50deg declination,” Astron. J. 146, 131 (2013).ANJOAA0004-6256http://dx.doi.org/10.1088/0004-6256/146/5/131Google Scholar



K. Collins and J. Kielkopf, “AstroImageJ: ImageJ for astronomy,” Astrophys. Sour. Code Libr. 1, 09001 (2013).Google Scholar


S. B. Howell, Handbook of CCD Astronomy, Cambridge University Press, New York (2006).Google Scholar


G. Hébrard et al., “WASP-52b, WASP-58b, WASP-59b, and WASP-60b: four new transiting close-in giant planets,” Astron. Astrophys. 549, A134 (2013).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361/201220363Google Scholar


F. Valsecchi and F. A. Rasio, “Planets on the edge,” Astrophys. J. Lett. 787, L9 (2014).AJLEEY0004-637Xhttp://dx.doi.org/10.1088/2041-8205/787/1/L9Google Scholar


J. Eastman, R. Siverd and B. S. Gaudi, “Achieving better than 1 minute accuracy in the heliocentric and barycentric Julian dates,” Publ. Astron. Soc. Pac. 122, 935–946 (2010).PASPAU0004-6280http://dx.doi.org/10.1086/653442Google Scholar


J. Eastman, B. S. Gaudi and E. Agol, “EXOFAST: a fast exoplanetary fitting suite in IDL,” Publ. Astron. Soc. Pac. 125, 83–112 (2013).PASPAU0004-6280http://dx.doi.org/10.1086/669497Google Scholar


P. Demarque et al., “Y2 isochrones with an improved core overshoot treatment,” Astrophys. J. Suppl. 155, 667–674 (2004).APJSA20067-0049http://dx.doi.org/10.1086/apjs.2004.155.issue-2Google Scholar


A. Claret and S. Bloemen, “Gravity and limb-darkening coefficients for the Kepler, CoRoT, Spitzer, uvby, UBVRIJHK, and Sloan photometric systems,” Astron. Astrophys. 529, A75 (2011).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361/201116451Google Scholar


G. Torres, “Substellar companion masses from minimal radial velocity or astrometric information: a Monte Carlo approach,” Publ. Astron. Soc. Pac. 111, 169–176 (1999).PASPAU0004-6280http://dx.doi.org/10.1086/pasp.1999.111.issue-756Google Scholar


T. Borkovits et al., “Transit timing variations in eccentric hierarchical triple exoplanetary systems. I. Perturbations on the time scale of the orbital period of the perturber,” Astron. Astrophys. 528, A53 (2011).AAEJAF0004-6361http://dx.doi.org/10.1051/0004-6361/201015867Google Scholar


Z. K. Berta et al., “The flat transmission spectrum of the super-Earth GJ1214b from wide field camera 3 on the Hubble Space Telescope,” Astrophys. J. 747, 35 (2012).ASJOAB0004-637Xhttp://dx.doi.org/10.1088/0004-637X/747/1/35Google Scholar


A. W. Howard et al., “A rocky composition for an Earth-sized exoplanet,” Nature 503, 381–384 (2013).NATUAS0028-0836http://dx.doi.org/10.1038/nature12767Google Scholar


F. Pepe et al., “An Earth-sized planet with an Earth-like density,” Nature 503, 377–380 (2013).NATUAS0028-0836http://dx.doi.org/10.1038/nature12768Google Scholar


G. W. Marcy et al., “Masses, radii, and orbits of small Kepler planets: the transition from gaseous to rocky planets,” Astrophys. J. Suppl. 210, 20 (2014).APJSA20067-0049http://dx.doi.org/10.1088/0067-0049/210/2/20Google Scholar


Jonathan Swift obtained a PhD from UC Berkeley in 2006 and served as the Submillimeter Array Postdoctoral Fellow at the Institute for Astronomy in Honolulu, HI, before taking on project management of MINERVA in 2012. He now has a joint position in the mathematics and science departments at the Thacher School in Ojai, CA, where he will also be directing the Thacher Observatory.

Biographies of the other authors are not available.

Jonathan J. Swift, Michael Bottom, John A. Johnson, Jason T. Wright, Nate McCrady, Rob Wittenmeyer, Peter P. Plavchan, Reed L. Riddle, Philip S. Muirhead, Erich Herzig, Justin Myles, Cullen H. Blake, Jason D. Eastman, Thomas G. Beatty, Stuart I. Barnes, Steven R. Gibson, Brian Lin, Ming Zhao, Paul Gardner, Emilio E. Falco, Stephen Criswell, Chantanelle Nava, Connor Robinson, David H. Sliski, Richard Hedrick, Kevin Ivarsen, Annie Hjelstrom, Jon de Vera, Andrew Szentgyorgyi, "Miniature Exoplanet Radial Velocity Array I: design, commissioning, and early photometric results," Journal of Astronomical Telescopes, Instruments, and Systems 1(2), 027002 (21 April 2015). http://dx.doi.org/10.1117/1.JATIS.1.2.027002







Optical fibers


Back to Top