Three Neptunes

Three Temperate Neptunes Orbiting Nearby Stars11affiliation: Based on observations obtained at the W. M. Keck Observatory, which is operated jointly by the University of California and the California Institute of Technology. Keck time was granted for this project by the University of Hawai‘i, the University of California, and NASA.

Benjamin J. Fulton1 2 3 , Andrew W. Howard1 , Lauren M. Weiss4 5 , Evan Sinukoff1 6 , Erik A. Petigura7 8 , Howard Isaacson4 , Lea Hirsch4 , Geoffrey W. Marcy4 , Gregory W. Henry9 , Samuel K. Grunblatt1 , Daniel Huber10 11 12 , Kaspar von Braun13 , Tabetha S. Boyajian14 , Stephen R. Kane15 , Justin Wittrock15 , Elliott P. Horch16 , David R. Ciardi17 , Steve B. Howell18 , Jason T. Wright19 20 , Eric B. Ford19 20
2affiliation: Institute for Astronomy, University of Hawai‘i, 2680 Woodlawn Drive, Honolulu, HI 96822, USA
17affiliation: National Science Foundation Graduate Research Fellow
*affiliation: bfulton@hawaii.edu
3affiliation: Department of Astronomy, University of California, Berkeley, CA 94720, USA
18affiliation: Ken and Gloria Levy Graduate Student Research Fellow
19affiliation: Natural Sciences and Engineering Research Council of Canada Graduate Student Fellow
5affiliation: California Institute of Technology, Pasadena, California, U.S.A.
20affiliation: Hubble Fellow
4affiliation: Center of Excellence in Information Systems, Tennessee State University, 3500 John A. Merritt Blvd., Box 9501, Nashville, TN 37209
7affiliation: SETI Institute, 189 Bernardo Avenue, Mountain View, CA 94043, USA
9affiliation: Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006
10affiliation: Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark
6affiliation: Lowell Observatory, Flagstaff, AZ 86001, USA
8affiliation: Yale University, New Haven, CT 06520, USA
11affiliation: Department of Physics and Astronomy, San Francisco State University, 1600 Holloway Avenue, San Francisco, CA 94132, USA
12affiliation: Department of Physics, Southern Connecticut State University, New Haven, CT 06515, USA
13affiliation: NASA Exoplanet Science Institute, Caltech, MS 100-22, 770 South Wilson Avenue, Pasadena, CA 91125, USA
14affiliation: NASA Ames Research Center, Moffett Field, CA 94035, US
15affiliation: Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA 16802
16affiliation: Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802
Abstract

We present the discovery of three modestly-irradiated, roughly Neptune-mass planets orbiting three nearby Solar-type stars. HD 42618 b has a minimum mass of , a semi-major axis of 0.55 AU, an equilibrium temperature of 337 K, and is the first planet discovered to orbit the solar analogue host star, HD 42618. We also discover new planets orbiting the known exoplanet host stars HD 164922 and HD 143761 ( CrB). The new planet orbiting HD 164922 has a minimum mass of  and orbits interior to the previously known Jovian mass planet orbiting at 2.1 AU. HD 164922 c has a semi-major axis of 0.34 AU and an equilibrium temperature of 418 K. HD 143761 c orbits with a semi-major axis of 0.44 AU, has a minimum mass of , and is the warmest of the three new planets with an equilibrium temperature of 445 K. It orbits exterior to the previously known warm Jupiter in the system. A transit search using space-based CoRoT data and ground-based photometry from the Automated Photometric Telescopes (APTs) at Fairborn Observatory failed to detect any transits, but the precise, high-cadence APT photometry helped to disentangle planetary-reflex motion from stellar activity. These planets were discovered as part of an ongoing radial velocity survey of bright, nearby, chromospherically-inactive stars using the Automated Planet Finder (APF) telescope at Lick Observatory. The high-cadence APF data combined with nearly two decades of radial velocity data from Keck Observatory and gives unprecedented sensitivity to both short-period low-mass, and long-period intermediate-mass planets.

Subject headings:
planetary systems — stars: individual (HD 42618) stars: individual (HD 164922) stars: individual (HD 143761)

1. Introduction

The mass function of extrasolar planets potentially offers rich clues to the processes that shape their growth and evolution. This mass function is known to rise with decreasing planet mass based on the discovery and characterization of planets orbiting nearby stars (Howard et al. 2010; Mayor et al. 2011). Although the Kepler mission discovered thousands of transiting exoplanets allowing for detailed characterization of the planet radius distribution (Howard et al. 2012; Fressin et al. 2013; Petigura et al. 2013a), we have a much sparser sample of planets orbiting nearby stars. We only know of 17 confirmed planets with measured minimum masses  and orbital periods days111Based on a query of exoplanets.org on 28 June 2016 (Han et al. 2014). These low-mass, temperate planets reside in an region of parameter space that must be explored in order to understand the formation of the extremely abundant population of close-in super-Earths.

The discovery that the occurrence rate of Jovian planets increases at orbital distances of 1–3 AU (Cumming et al. 2008) has been suggested to be a sign that the ice line is important to the formation of Jovian planets. The increased abundance of solids in the protoplanetary disk beyond the ice-line is expected to speed up the coagulation of planetesimals and  cores that can undergo runaway gas accretion before the gas in the disk is dissipated (Ida & Lin 2008).

The formation of close-in intermediate mass planets known as super-Earths or mini-Neptunes presents some challenges for planet formation theory. It was initially thought that such planets should either remain small terrestrial planets, or if they grow to large cores then they should quickly accrete substantial nebular gas and grow to be gas giants (Ida & Lin 2004; Mordasini et al. 2009). Alternatively, planets that form in situ likely require either a protoplanetary disk that is much more massive than the minimum mass solar nebula (Hansen & Murray 2012; Chiang & Laughlin 2013), an extremely metal rich disk, or fine tuning of the formation timescales (Lee et al. 2014). New models featuring gas-drag driven “pebble accretion” which offers a mechanism to transport solids in the disc (Chatterjee & Tan 2014), or the delayed formation of super-Earth cores until the gas disk begins to dissipate (Lee & Chiang 2015) could facilitate in situ formation of super-Earths and/or mini Neptunes.

Alternatively, migration from further out in the disk where there is plenty of material to form massive planet cores or super-Earths could explain the presence of these planets near their host stars. We would expect to see multi-planet systems in resonant chains if a slow, smooth migration were the dominant mechanism and this is not observed (Veras & Ford 2012; Fabrycky et al. 2014). Another indication that migration is at play would be an increased occurrence rate of super-Earths at large orbital separations. The occurrence rates of super-Earths as a function of orbital period appears fairly flat to periods as long as 200 days (Petigura et al. 2013a), however the sensitivity to long-period planets from transit surveys is very low. In addition, the planet radius distribution can be influenced by a variety of different factors including stellar irradiation, and thermal evolution (Lopez 2014). Only a very small percentage by mass of volatiles can significantly inflate a planet’s radius and hide the fundamental properties of the planet that encode information about the formation mechanism (Lopez & Fortney 2014). RV surveys are now starting to discover a statistically useful sample of super-Earth to Neptune mass planets at larger orbital separations and lower stellar irradiance that will help to map out the details of the mass function for long-period super-Earth to Neptune mass planets.

The Eta-Earth RV survey of nearby stars (Howard et al. 2010) was conducted using the HIRES spectrograph at Keck observatory (Vogt et al. 1994). They searched for planets in a volume-limited sample of 166 nearby G and K dwarfs. With the catalogue of planets detected in this survey (Howard et al. 2009, 2011a, 2011b, 2014), and the completeness limits calculated for each star they were able to measure the occurrence rate of small, short-period planets and show that planets with masses of 3-30  are much more common than planets larger than 30 . Although these low-mass planets are common, they are still very difficult to detect, requiring 100 measurements per star with precision. The number of nearby stars for which large, high-precision radial velocity datasets exist is increasing thanks to the proliferation of dedicated and robotic radial velocity facilities such as the APF and MINERVA (Swift et al. 2015) and the ongoing long-term surveys from Keck, HARPS-N (Motalebi et al. 2015), and HARPS (Pepe et al. 2004).

We are currently using the APF telescope to conduct a RV survey of 51 of the brightest, and least chromospherically active, stars from the Eta-Earth survey. We capitalize on the robotic nature of the telescope to monitor the stars at high cadence for the entire four year duration of the survey. This survey builds on the Eta-Earth Survey, but with improved Doppler precision due to the high observing cadence and larger number of measurements. We will measure the occurrence rate and mass function of small planets in our local neighborhood using the new planets discovered by the APF-50 survey and the set of planets already known to orbit stars in our sample. With a larger sample of planets with measured masses in the 3-30  range we will measure the mass function for small planets with higher mass resolution. Combining the mass function from this survey with the size distribution from Kepler, we will probe the density and core mass properties of super-Earths to inform formation theories of the galaxy’s most abundant planets.

Nearby G and K dwarf stars are observationally advantageous. High SNR spectroscopy can often be obtained with relatively short exposure times, facilitating many precise time series RV measurements with sensitivity to planets with small Doppler amplitudes. The stars can be characterized precisely using spectroscopy and asteroseismology (when available) that improve estimates of the star and planet properties. The advantages of nearby, bright targets are critical for characterization of the planets atmospheres transmission, emission, and direct spectroscopy.

Here we present the discovery of three roughly Neptune mass planets orbiting bright stars within 25 pc. This paper is structured as follows. Our observational setup and RV measurements are described in Section 2. In Section 3 we discuss our derived stellar properties for each of the three stars and compare with previous literature studies. We describe our methods used to discover these planets in Section 4. We describe our modeling procedure used to obtain the final adopted parameters and their associated uncertainties and our various tests to ensure that the signals are planetary in nature in Section 4.2. We analyze photometry of each of the three systems in Section 5, and conclude with a summary and discussion in Section 6.

2. Radial Velocity Measurements

We collected 571 RV measurements of HD 42618 using Keck/HIRES (Vogt et al. 1994) and 35 measurements using the Levy spectrograph on the Automated Planet Finder (APF, Vogt et al. 2014; Radovan et al. 2014) over the past 19 years starting in 1996. For each RV measurement the starlight is passed through a cell of gaseous iodine that imprints a dense forest of molecular absorption lines onto the stellar spectrum and serves as both a wavelength and point spread function (PSF) reference. We also collected a single set of iodine-free observations of this star that was deconvolved with the instrumental PSF and used as a model of the intrinsic stellar spectrum. Each observation was forward modeled as the intrinsic stellar spectrum doppler shifted by an arbitrary amount, then multiplied by the transmission of iodine, and convolved with the instrumental PSF modeled as a sum of 13 Gaussians with fixed widths and positions but heights free to vary (Butler et al. 1996). The Levy slit-fed spectrograph also relies on an iodine cell for precise RVs and our observational setup is described in detail in Fulton et al. (2015).

Our setup was identical for the three stars. We collected 328 Keck/HIRES measurements and 73 APF measurements for HD 164922 over the past 19 years. All of the APF measurements and 244 of the Keck measurements are new since the publication of Butler et al. (2006). For HD 143761 we obtained 519 RV measurements using Keck/HIRES and 157 measurements using the Levy spectrograph on the APF. The Keck observations of HD 143761 started in 2006 and the total observational baseline is 8 years. We do not include the Lick 3.0 m data from Noyes et al. (1997) and Butler et al. (2006) in our analysis. We find that including the Lick data does not significantly improve the uncertainties on any of the orbital parameters and it adds an additional source of systematic uncertainty that is less well characterized and understood.

Our Doppler pipeline has been tuned in small ways over the years to improve RV precision. Here we describe a new pipeline improvement that decorrelates the measured RVs with nuisance parameters from the model and spectrometer state parameters. This decorrelation offers modest improvements in Doppler precision (1 m s) and is only applied to time series RV of stars for which the number of spectra greatly exceeds the number of potential decorrelation parameters. The nuisance parameters in Doppler analysis include descriptions of the point spread function (PSF) over the spectral format and the wavelength solution. The PSF is parameterized as a sum Gaussians with fixed widths and centers, but variable amplitudes. We also have a wealth of information about the weather and environment inside and outside the spectrograph extracted from the FITS headers of the raw APF spectra. Some environmental information is available in the FITS headers for the Keck data, but we have not yet implemented a system to extract these values for the tens of thousands of Keck spectra taken over the last 20 years. For Keck data we only include the nuisance parameters that are part of the forward modeling process. We clean the RVs of systematic trends by removing any correlations that these parameters show with the final RVs. We search for significant correlations of the RVs with all of the PSF parameters by calculating the Spearman rank correlation coefficient (Spearman 1904) and flag any parameters that show correlation coefficients greater than 0.1. The flagged parameters are included in a multivariate ordinary least squares linear regression using the STATSMODELS222https://pypi.python.org/pypi/statsmodels package in Python. This model for RV as a function of all parameters included in the fit is then subtracted from the raw RVs. This process is done blindly in the planet discovery/identification phase but once planet candidates are identified in a given dataset we first model the system to find the best fit N-planet Keplerian model then perform the detrending procedure on the residuals to this model. The detection of all newly discovered planets in this work does not depend on this detrending and they are easily identified at high significance in either the detrended or non-detrended datasets.

We reject measurements with low signal-to-noise ratios (SNR per pixel) and/or uncertainties more then 9 larger then the median uncertainty, which results in the omission of 1% of the data for each of the three stars. Since these stars are exceptionally bright we almost always collect three consecutive measurements in order to average out RV shifts caused by p-mode oscillations (Dumusque et al. 2011). The three measurements are then binned together before the stellar jitter is added in quadrature during the modeling process (see Section 4.2). This effectively reduces the weight of the three measurements to that of a single measurement, but averages out some of the astrophysical noise in the process and prevents time-correlated instrumental systematic noise from biasing the results. We also extract the Ca ii H & K activity index () using the technique of Isaacson & Fischer (2010) for every RV measurement on both Keck and APF, however there may be an arbitrary zero point offset in the  values between the Keck and APF values. The uncertainties for the  measurements are systematically limited to 0.002 for Keck and 0.004 for APF. This was estimated by measuring the standard deviation of all measurements of the extremely chromospherically quiet star, HD 10700. All RV measurements and the associated  values can be found in Table 1. We include only the detrended velocities in Table 1 but the full set of environmental and PSF parameters for each observation along with the non-detrended velocities can be downloaded from https://github.com/bjfultn/three_neptunes.

HD RV11Zero point offsets between instruments have not been removed and must be fit as free parameters when analyzing this dataset Unc. Inst.22k = pre-upgrade Keck/HIRES, j = post-upgrade Keck/HIRES, a = APF 33Uncertainties on  are 0.002 for all Keck measurements and 0.004 for all APF measurements
(– 2440000) (m s) (m s)
42618 2450366.126333 02.85 1.12 k
42618 2453694.093412 01.17 1.13 j 0.161
164922 2454777.744397 09.19 1.04 j 0.152
164922 2457267.662041 04.43 1.90 a 0.145
143761 2455455.762546 41.98 1.27 j 0.149
143761 2457292.685768 08.08 2.92 a 0.135
footnotetext: (This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.)
Table 1Radial Velocities

3. Stellar Properties

3.1. Hd 42618

HD 42618, also known as HIP 29432 and Gl 3387, is a well studied solar analogue located at a distance of 23.5 pc (van Leeuwen 2007). The star was not previously known to host any exoplanets. We analyzed 5 high SNR Keck-HIRES spectra (described below) using SpecMatch (Petigura 2015) to obtain the mean spectroscopic parameters listed in Table 2. SpecMatch uses trilinear interpolation to synthesize high resolution model spectra from the Coelho (2014) grid of models for any set of arbitrary stellar parameters (, , , and ) that are contained within the limits of the model grid. The interpolated models are then compared to the observed spectrum. We maximize the likelihood () to determine the optimal stellar parameters, where is summed over the extracted spectral pixels and normalized by the flux uncertainties.

Our spectral analysis is consistent with the results of Valenti & Fischer (2005) who extracted spectroscopic parameters using Spectroscopy Made Easy (SME) and found K, , and . HD 42618 is chromospherically quiet with (Isaacson & Fischer 2010). It was deemed to be a good solar-analog based on a very similar chemical abundance pattern to the Sun (Morel et al. 2013a). Those authors also derive K, and , consistent with our SpecMatch results. Ramírez et al. (2014) measured the fundamental parameters of HD 42618 differentially relative to the Sun which allowed them to obtain highly precise values for K, , and that show good agreement with our results. It has also been noted that HD 42618 shows a low lithium abundance of (Ramírez et al. 2012) similar to that of the Sun ( Asplund et al. 2009). Our adopted stellar mass and radius for HD 42618 are based on the relations of Torres et al. (2010) using our spectroscopic constraints on , , and . HD 42618 was also target of the CoRoT mission (Baglin et al. 2009), with a preliminary detection of solar-like oscillations presented by Barban et al. (2013). We performed an independent asteroseismic analysis of the CoRoT photometry (see Section 5.1.1), which yielded a mass and radius in agreement with our adopted values.

Parameter HD 42618 HD 164922 HD 143761
Spectral type G4V11MCMC jump parameters that were modified from the physical parameters in order to speed convergence and avoid biasing parameters that must physically be finite and positive. G9V66footnotemark: G0V88footnotemark:
(mag) 0.65722footnotemark: 0.80077footnotemark: 0.600
(mag) 6.83922footnotemark: 6.9977footnotemark: 5.4188footnotemark:
(mag) 5.701 0.02333footnotemark: 5.553 0.02633footnotemark: 4.0933footnotemark:
(mag) 5.385 0.02433footnotemark: 5.203 0.01733footnotemark: 3.9933footnotemark:
(mag) 5.301 0.02033footnotemark: 5.113 0.02033footnotemark: 3.89 0.0588footnotemark:
Distance (pc) 23.500.3044footnotemark: 22.130.2744footnotemark: 17.2360.02444footnotemark:
(K) 5727 60 5293 32 5627 541111footnotemark:
log  (cgs) 4.44 0.07 4.387 0.014 4.121 0.018
 (dex) 0.09 0.04 0.16 0.05 0.31 0.05
 sin  (km s)
() 0.98 0.17 0.703 0.017 1.706 0.0421111footnotemark:
 () 1212footnotemark: 0.874 0.012 0.889 0.030
 () 0.999 0.0871212footnotemark: 0.999 0.017 1111footnotemark:
5.0155footnotemark: 5.0655footnotemark: 5.0555footnotemark:
0.15755footnotemark: 0.15455footnotemark: 0.15055footnotemark:
footnotetext: References - (1) Medhi et al. (2007); (2) Koen et al. (2010); (3) Cutri et al. (2003); (4) van Leeuwen (2007); (5) Isaacson & Fischer (2010); (6) Gray et al. (2003); (7) Butler et al. (2006); (8) van Belle & von Braun (2009); (9) Noyes et al. (1997); (10) van Belle & von Braun (2009); (11) von Braun et al. (2014); (12) Torres et al. (2010)
Table 2Adopted Stellar Properties

3.2. Hd 164922

HD 164922, also known as HIP 88348 and Gl 9613, is a bright, chromospherically inactive (, Isaacson & Fischer 2010) G9 V dwarf located 22.1 pc away (van Leeuwen 2007). It was previously known to host a single Saturn-mass planet orbiting with a semi-major axis of 2.1 AU (Butler et al. 2006). This target was one of several selected for more intensive long-term RV monitoring by Keck/HIRES based on both the stellar properties, and the mass and orbit of the previously detected planet making the system particularly well-suited for detecting additional low-mass planets. It was also on the Eta-Earth target list as part of a deep Doppler survey for low-mass planets (Howard et al. 2010).

We measured the stellar radius for HD 164922 using the CHARA Array. Interferometric observations of HD 164922 were taken on 2012 May 13 and 14 using the Pavo beam combiner (ten Brummelaar et al. 2005; Ireland et al. 2008). Observations of the science target were interleaved with the calibrator stars HD 164900, HD 161019, and HD 165373 (Bonneau et al. 2006, 2011). The data were reduced and calibrated using the standard data reduction pipeline (for details see White et al. 2013). We use the -band limb-darkening coefficient from Claret & Bloemen (2011), , to determine a limb-darkened angular diameter  mas (Figure 1).

Figure 1.— Observed squared visibility versus spatial frequency for HD 164922 (blue diamonds). The red line shows the best fit limb-darkened model. See Section 3.2 for a discussion of the implications of this figure.

Combining the angular diameter with the parallax yields a stellar radius of  R. We determine a stellar bolometric flux of  erg/s/cm by fitting a spectral template from Pickles (1998) to flux calibrated photometry after applying revised filter profiles from Mann & von Braun (2015). This translates to a luminosity  L. Lastly, we use our measured angular diameter with the star’s bolometric flux to derive an empirical effective temperature  K.

We then use the interferometrically determined parameters to inform a SpecMatch analysis of a stack of 5 high SNR APF spectra using an iterative technique. An initial uninformed SpecMatch analysis of the APF spectra (all priors uniform) gives K, and . We use this spectroscopically measured  combined with the  and  determined from the CHARA data along with existing ,, and photometry (Cutri et al. 2003) as Gaussian priors in a fit to Dartmouth isochrones (Dotter et al. 2008) using the isochrones package (Morton 2015)333https://github.com/timothydmorton/isochrones. This gives the first estimate of the full set of stellar parameters. We then re-run SpecMatch with Gaussian priors applied to  and  from the isochrones output. The full likelihood with the Gaussian priors is

(1)

In this case the changes to stellar parameters become negligible after two iterations and the process is halted. The resulting stellar parameters are listed in Table 2.

Since the star is bright and already a known planet host it has been the subject of many spectroscopic studies. Santos et al. (2013) find K, , and , all consistent with our analysis to within 1 . Valenti & Fischer (2005) find a significantly higher value for and K but a consistent metallicity value of . Ghezzi et al. (2010) measure K, , and also using high-resolution spectroscopy which, except for , is also consistent with our analysis to within 1 . As seen by the range of  values obtained by these various studies, it can be difficult to pin down the stellar gravity from high resolution spectra alone. Since our value is constrained via the direct measurement of the stellar radii from interferometry we adopt the new slightly lower value for .

3.3. Hd 143761

HD 143761 is the closest and brightest star of the three studied in this work. The star is also known as Corona Borealis, HIP 78459, and Gl 9537. It is a slightly evolved naked eye () G0 V star (van Belle & von Braun 2009) located at a distance of 17.236 pc (van Leeuwen 2007). It was previously known to host a warm Jupiter-mass planet with an orbital period of 39 days (Noyes et al. 1997). This star was also part of the Eta-Earth survey and was independently selected for intensive long-term RV monitoring based on both the stellar properties, and the mass and orbit of the previously detected planet. Like HD 42618 and HD 164922 this star is chromosherically quiet with (Isaacson & Fischer 2010). As with HD 164922 we performed an iterative interferometric+spectroscopic analysis using the CHARA results from von Braun et al. (2014) and a stack of 5 high SNR APF spectra. The likelihood was essentially the same as Equation 1 but with the  and  values for HD 143761 substituted in the last two terms.

Valenti & Fischer (2005) measure K, and Fuhrmann et al. (1998) measure K both using high resolution spectroscopy. Both of these values are significantly hotter then our adopted value of K from von Braun et al. (2014) . We chose to adopt the value from von Braun et al. (2014) in order to maintain self-consistency with the interferometrically measured stellar radius and luminosity. Our metallicity value is also consistent within 1 to that of Fuhrmann et al. (1998) and but is significantly lower than that of Valenti & Fischer (2005). Valenti & Fischer (2005) again measure a significantly higher value for , but our value is consistent with from Fuhrmann et al. (1998).

4. Keplerian Analysis

4.1. Discovery

We discovered each of the three new planets and re-discovered the previously known planets using a technique essentially identical to that of Fulton et al. (2015). In brief, we calculate a two-dimensional Lomb-Scargle periodogram (2DKLS, O’Toole et al. 2009) to look for significant periodic signals that are well fit by a Keplerian orbital model. Our implementation of the 2DKLS periodogram incorporates arbitrary zero-point offsets between each instrument. The periodogram power () represents the improvement to the statistic relative to that of a baseline fit. When searching for the first planet in a system the baseline fit is simply a flat line or linear trend. If any significant signals are found after the first iteration the baseline model then becomes the single planet Keplerian model and we calculate the improvement to when a second planet is added, without subtracting the first. We repeat this process until no more significant peaks are found in the 2DKLS periodogram. We start the search assuming no known signals in order to ensure that the previously published planets can be automatically detected using our pipeline. An initial jitter term of 2.0 m sis added in quadrature with the RVs before starting the 2DKLS search in order to ensure fair weighting between the Keck and APF data sets. We expect both data sets to be limited by instrumental/astrophysical systematics rather than photon noise.

The discovery pipeline is completely automated in order to facilitate injection recovery tests that will allow us to characterize the pipeline completeness for future occurrence analysis(Howard & Fulton 2016). We calculate an empirical periodogram false alarm probability (eFAP) by fitting a power law to the distribution of periodogram values between the 50th and 97th quartiles. This fit provides an estimate of the significance of periodogram peaks of a given value. When multiplied by the number of independent test periods, the fit gives the approximate probability that we would find a peak of a given value within any particular periodogram. Any periodogram peak with an eFAP below 0.1% is automatically considered a viable candidate and the search is continued until no more periodogram peaks fall above the 0.1% eFAP threshold. Further details of the automated planet detection pipeline can be found in Howard & Fulton (2016). We note that the eFAP metric is used simply to automatically identify candidates. The significance of the corresponding periodogram peaks are checked using the bootstrapping technique described in Section 4.3. Each of the previously known planets were re-discovered with eFAPs much less than those of the new planets announced in this work.

We discover two significant signals in the RV time series of HD 42618. One with a long period of 4850 days, and a second at a period of 149.6 days. Upon inspection of the Ca ii H & K activity index time series we notice that this index shows a periodicity with a period very similar to that of the long period RV signal. The period of 4850 days is also very similar to the period of the sun’s magnetic activity cycle. We conclude that this is likely the signature of the stellar magnetic activity cycle and not the signature of an orbiting planet. We include this long-period signal as an additional eccentric Keplerian in all further modeling. HD 42618 b is easily detected in the Keck data alone and the combined Keck+APF dataset but we don’t yet have enough APF measurements to detect it in the APF data alone. Figure 2 shows the most likely model from the posterior of the two-Keplerian model, the 2DKLS periodogram used to discover HD 42618 b, and the RVs phased to the orbital period of planet b.

Wright et al. (2007) mentioned a candidate planetary signal with a period of 75.8 d and m sorbiting HD 164922 but did not have sufficient data to claim a significant detection. With seven years of additional Keck data, and 2 years of APF data we can firmly establish this signal as being coherent and persistent as expected for the Doppler motion caused by an orbiting planet. The short period planet is easily detected in either the APF or Keck data alone. The long-period planet can only be detected in the long baseline Keck data but we do observe a linear RV trend that emerged during the most recent APF observing season which is a result of the massive outer planet. Figures 3 and 4 show the most likely model from the posteriors for the two planet Keplerian model and the 2DKLS periodograms used to discover/re-discover each of the two planets.

We discover a super-Neptune mass planet orbiting HD 143761 exterior to the known Jupiter mass planet that has an orbital period of 39 days. The new planet has an orbital period of 102 days and a semi-amplitude of 3.7 m s. Each planet is discovered with very high significance in both the Keck and Keck+APF datasets individually. The most likely model from the posteriors for the two planet model and the 2DKLS detection periodograms are shown in Figure 5.

Figure 2.— One-planet Keplerian orbital model plus one additional long-period Keplerian to model the stellar magnetic activity cycle for HD 42618. The most likely model is plotted but the orbital parameters annotated on the figure and listed in Tables 3 and 6 are the median values of the posterior distributions. The process used to find the orbital solution is described in 4.2. a) Full binned RV time series. Open black squares indicate pre-upgrade Keck/HIRES data (see §2), open black circles are post-upgrade Keck/HIRES data, and filled green diamonds are APF measurements. The thin blue line is the most probable 1-planet plus stellar activity model. We add in quadrature the RV jitter term listed in Table 3 with the measurement uncertainties for all RVs. b) Residuals to the most probable 1-planet plus stellar activity model. c) Binned RVs phase-folded to the ephemeris of planet b. The long-period stellar activity signal has been subtracted. The small point colors and symbols are the same as in panel a. For visual clarity, we also bin the velocities in 0.08 units of orbital phase (red circles). The phase-folded model for planet b is shown as the blue line. d) 2DKLS periodogram showing the improvement to for a model including the long period activity signal and a single planet compared to a model that only includes the activity signal.
Figure 3.— Two-planet Keplerian orbital model posterior distributions for HD 164922. The most likely model from the posterior distribution is plotted while the orbital parameters annotated on the figure and listed in Tables 4 and 6 are the median values of the posterior distributions. The process used to find the orbital solution is described in 4.2. a) Full binned RV time series. Open black squares indicate pre-upgrade Keck/HIRES data (see §2), open black circles are post-upgrade Keck/HIRES data, and filled green diamonds are APF measurements. The thin blue line is the most probable 2-planet model. We add in quadrature the RV jitter term listed in Table 4 with the measurement uncertainties for all RVs. b) Residuals to the most probable 2-planet model. c) Binned RVs phase-folded to the ephemeris of planet b. The Keplerian orbital model for planet c has been subtracted. The small point colors and symbols are the same as in panel a. For visual clarity, we also bin the velocities in 0.08 units of orbital phase (red circles). The phase-folded model for planet b is shown as the blue line. d) 2DKLS periodogram comparing a model including only the long period planet to the two planet model. Panels e) and f) are the same as panels c) and d) but for planet HD 164922 c. The shaded region of panel a is re-plotted in Figure 4.
Figure 4.— Recent RVs for HD 164922 highlighted in the grey box of Figure 3 panel a) featuring the high cadence APF observations collected during the most recent observational season.
Figure 5.— Same as Figure 3 but for planets HD 143761 b and c.

4.2. Characterization

We estimated orbital parameters and their associated uncertainties using the ExoPy Differential Evolution Markov Chain Monte Carlo (DE-MCMC, Ter Braak 2006) modeling code and a technique identical to that of Fulton et al. (2013), Knutson et al. (2014), and Fulton et al. (2015).

Our multi-planet RV model is a sum of Keplerian single-planet models over all planets in the system. For each single-planet Keplerian model (indexed by ) we compute posteriors for the orbital period (), time of inferior conjunction (), eccentricity (), argument of periastron of the star’s orbit (), velocity semi-amplitude (). We also compute posteriors for the offsets between pre-upgrade, post-upgrade, and the APF datasets (), and an RV jitter term with the specific prior described in Fulton et al. (2015) and Johnson et al. (2011). In order to speed convergence, we choose to re-parameterize some of the physical parameters as , , and . The full likelihood for this Keplerian model is

(2)

summed over RV measurements (v) with associated uncertainties . is the Keplerian model for observation . We assign uniform priors to , , , , , and for each instrument. We follow the prescription of Eastman et al. (2013) checking for convergence at regular intervals during the MCMC runs by calculating the Gelmin-Rubin statistic and the number of independent draws (, Ford 2006). We consider the chains well mixed and halt the MCMC run when the Gelmin-Rubin statistic is within 1% of unity and for all free parameters. All of the adopted median values and 68% confidence intervals of the posterior distributions are listed in Tables 35.

Parameter Value Units
Modified DE-MCMC Step Parameters11MCMC jump parameters that were modified from the physical parameters in order to speed convergence and avoid biasing parameters that must physically be finite and positive.
log() 2.17497 log(days)
log() 0.276 log(m s)
Model Parameters
149.61 days
2456670.2
0.19
101 degrees
1.89 m s
0.61 m s
0.71 m s
0.0 m sday
0.0 m sday
2.34 m s
Table 3Orbital Parameters for HD 42618
Parameter Value Units
Modified DE-MCMC Step Parameters11MCMC jump parameters that were modified from the physical parameters in order to speed convergence and avoid biasing parameters that must physically be finite and positive.
log() 3.08 log(days)
log() 0.854 m s
log() 1.87947 log(days)
log() 0.346 log(m s)
Model Parameters
1201.1 days
2456778
0.126
129 degrees
7.15 m s
75.765 days
2456277.6
0.22
81 degrees
2.22 m s
0.23 m s
1.02 m s
0.2 m s
0.0 m sday
0.0 m sday
2.63 m s
Table 4Orbital Parameters for HD 164922
Parameter Value Units
Modified DE-MCMC Step Parameters11MCMC jump parameters that were modified from the physical parameters in order to speed convergence and avoid biasing parameters that must physically be finite and positive.
log() 1.600382 log(days)
log() 1.8279 m s
log() 2.01091 log(days)
log() 0.573 log(m s)
Model Parameters
39.8458 days
2455759.091
0.0373
270.6 degrees
67.28 m s
102.54 days
2455822
0.052
190 degrees
3.74 m s
-0.6 m s
-0.7 m s
0.0 m sday
0.0 m sday
2.57 m s
Table 5Orbital Parameters for HD 143761
Parameter Value Units
HD 42618
-0.009
0.14
0.554 AU
119.1
14.4
**Stellar irradiance received at the planet relative to the Earth. 3.16
****Assuming a bond albedo of 0.32 (Demory 2014). 337 K
HD 164922
-0.025
0.032
2.115 AU
454.9
107.6
**Stellar irradiance received at the planet relative to the Earth. 0.1578
****Assuming a bond albedo of 0.32 (Demory 2014). 159.4 K
0.003
0.079
0.3351 AU
72.1
12.9
**Stellar irradiance received at the planet relative to the Earth. 6.29
****Assuming a bond albedo of 0.32 (Demory 2014). 400.5 K
HD 143761
7e-05
-0.0072
0.2196 AU
34.66
332.1
**Stellar irradiance received at the planet relative to the Earth. 34.7
****Assuming a bond albedo of 0.32 (Demory 2014). 614.0 K
0.0001
-0.0001
0.4123 AU
65.1
25
**Stellar irradiance received at the planet relative to the Earth. 9.85
****Assuming a bond albedo of 0.32 (Demory 2014). 448.1 K
Table 6Derived Planet Properties

4.3. Bootstrap False Alarm Assessment

We conduct a bootstrap false alarm assessment to verify and double check that the periodogram peaks with low eFAPs are indeed statistically significant periodic signals and not caused by random fluctuations of noise. For all three stars we scramble the RV time series 1000 times and recalculate the 2DKLS periodogram searching for N+1 planets where N is the number of previously published planets in the system. We record the highest periodogram value from each trial and plot the distribution of periodogram peak heights relative to the periodogram peak values corresponding to the newly discovered planets. These distributions are plotted in Figure 6. The periodogram peak heights corresponding to each of the new planets are well separated from the distribution of peaks in the scrambled RV trials. This indicates that the probability that random noise could conspire to create the periodogram peaks used to detect the planets are . However, a visual inspection of the distribution of periodogram peak heights in Figure 6 suggests that the FAPs are likely much lower.

Figure 6.— Graphical representation of the bootstrap false alarm tests described in Section 4.3. Top: Distribution of maximum periodogram peak heights for 1000 2DKLS periodograms of scrambled RV time series for HD 42618. The long period activity signal was subtracted before scrambling the data set. The vertical dashed blue line marks the height of the original periodogram peak for planet b which is clearly separated from the distribution of peaks caused by random fluctuations. Middle: Same as the top panel for planet HD 164922 c. Bottom: Same as the top panel for planet HD 143761 c. In each case, the scrambled RVs generate peaks with significantly lower power than the power observed from the new planetary signals.

4.4. Brown Dwarf Companion to HD 143761?

HD 143761 b was one of the first exoplanets discovered (Noyes et al. 1997). Gatewood et al. (2001) later reported that they detected the signature of HD 143761 b in astrometric data from the Multichannel Astrometric Photometer and Hipparcos. Zucker & Mazeh (2001) were quick to point out that the statistical significance of this astrometric detection is only 2 . Over a decade after the discovery of HD 143761 b, Reffert & Quirrenbach (2011) claimed to detect the same astrometric signal of the warm Jupiter in a re-reduction of Hipparcos data. The astrometric orbit suggests that the system is nearly face on with an inclination between 0.4 and 0.7 degrees. This would imply that, after correcting for the viewing angle, HD 143761 b is not a planet but instead a low mass M star with M.

4.4.1 Interferometry of HD 143761

Long-baseline interferometry is sensitive to some stellar binaries. HD 143761 was observed interferometrically with the CHARA Array (von Braun et al. 2014). If we assume a face-on orbit of HD 143761 b at a distance of 17.2 pc, then the angular distance between it and the principal component is 14 milliarcseconds independent of phase angle. This is detectable as a separated fringe packet at -band with CHARA (Farrington et al. 2010), provided the brightness contrast is not larger than (Farrington et al. 2010; Raghavan et al. 2012).

Using the mass constraints from the astrometric orbit, a comparison with the Dartmouth isochrones for metallicities spanning the 1  uncertainty of HD 143761 from Table 2 corresponds to an apparent magnitude at a distance of 17.2 pc of around 10. This implies a , which is significantly below what could be detected as a separate component in CHARA data (cf. equation A5 in Boyajian et al. 2008). Since component c’s apparent magnitude is much fainter than component b’s, its detection is impossible in the CHARA data. As expected, an inspection of the CHARA data used in von Braun et al. (2014) did not yield any indication of additional fringe packets.

We also obtained imaging observations of CrB using the Differential Speckle Survey Instrument (DSSI) on Gemini-North during the nights of July 19, 24, and 25 in 2014. The DSSI camera is a dual-channel speckle imaging system, expounded upon in more detail by Horch et al. (2009, 2011). Observations were acquired using red and near-infrared filters centered on 692 nm and 880 nm respectively. Our instrument setup is the same as that described in Horch et al. (2012) and our analysis methodology is outlined by Kane et al. (2014). Briefly, we estimate the limiting magnitude (difference between local image maxima and minima) as a function of target separation resulting in a 5 detection curve. More details on the derivation of the DSSI detection limits can be found in Howell et al. (2011). All of our DSSI CrB observations show no evidence of a stellar companion to the host star. Figure 7 shows the detection curve from the 880 nm image acquired for CrB on the night of July 25 2014. The dashed curve is the cubic spline interpolation of the 5 detection limit from 01 to 12. The results exclude companions with and at separations of 01 and 14 respectively. Given the distance of CrB of 17.236 pc, these angular separations correspond to a physical exclusion range of 1.7–24.1 AU. We can thus rule out stellar companions in close proximity to the host star, supporting the evidence that the system is not a face-on triple star system, but instead a multi-planet system viewed at moderate to high inclination.

Figure 7.— Limiting magnitude as a function of separation from the CrB. Also shown are a cubic spline interpolation of the 5 detection limit (dashed line) and limiting magnitudes for 01 and 02.

4.4.2 Stability of the HD 143761 System

Two planet systems on circular orbits are likely to be unstable if Gladman (1993) for all mutual inclinations, where and R is the Hill radius,

(3)

Sky-projected inclinations smaller than 4 degrees combined with the  constraints listed in Table 6 imply large companion masses and push below the stability threshold. If we assume that the system is dynamically stable then the presence of HD 143761 c rules out the low inclination orbit found by Gatewood et al. (2001), and Reffert & Quirrenbach (2011).

HD 164922 b and c are widely separated and intuitively we would expect them to be in a stable configuration. For completeness, we calculate for the minimum masses which suggests that this system configuration is likely to be dynamically stable for a long time.

4.5. Additional Planet Candidates

There is an additional significant periodic signal in the RV data for HD 42618 at a period of 388 days and a velocity semi-amplitude of 2 m/s. This would be a = 22 planet orbiting just outside 1 AU. The periodogram peak for this candidate falls above the eFAP threshold automatically calculated by our discovery pipeline. However, we do not believe that we have enough evidence to claim a concrete detection of a bona fide planet due to the proximity of this period to 1 year and incomplete phase coverage of the orbit. Telluric contamination of the template or problems with the barycentric correction could inject a false signal with a period near 1 year (Wright & Eastman 2014; Fischer et al. 2016, S. Wang personal communication). The amplitude and periodicity of this signal depends on the outcomes of modeling HD 42618 b, and MCMC runs for models including this candidate fail to converge. We will continue to monitor this star intensively with both Keck and the APF to confirm or refute this planet candidate in the upcoming years.

We find a candidate periodicity in the HD 164922 system with a period of 41.7 d and an amplitude of 1.9 m s. The eFAP of this 2DKLS periodogram peak is 0.00098 and falls just above our 0.1% eFAP threshold. However, we do not consider this to be a viable planet candidate due to its marginal detection and proximity to the expected rotation period for this star (44 days, Isaacson & Fischer 2010). Further monitoring and a detailed analysis that includes the effects of rotational modulation of starspots is needed to determine the nature of this signal.

There is no evidence for significant periodic signals from other candidates in the periodograms for HD 143761. However, visually there appears to be some long-period structure in the residuals to our most probable model (see Figure 5). This marginal variability, if real, likely has a period of years and an amplitude of only a few m sand it appears to be at a shorter period than the stellar magnetic activity cycle as seen in the  values for this star. Long-term monitoring of this target is required to determine if this signal is real and the signature of a planetary companion.

4.6. Chromospheric Activity

These stars were all selected to be part of the APF-50 survey of nearby stars due, in part, to their extremely low mean chromospheric activity of . However, in the case of HD 42618 we do detect significant long-period variability in the  values that is strongly correlated with the RVs (see Figure 8) that is likely the signature of the stellar magnetic activity cycle. We do not find any significant periodic signals in the  values after removal of this long-period trend that might be the signature of rotation. However, we clearly identify the rotation period of the star to be 16.9 days in CoRoT photometry (see Section 5.1). We account for the activity cycle in the RV data of HD 42618 by including an additional long-period Keplerian signal in the model.

HD 164922 shows only a linear trend in the  values but we do not detect the effect of this change in chromospheric activity in the residuals to the two planet fit. There is also a very weak peak in the Lomb-Scargle periodogram (L-S, Lomb 1976; Scargle 1982) of the  values at 37.8 days (see Figure 9). This may be the signature of stellar rotation since this is near the expected period for a star of this type and age (Isaacson & Fischer 2010). However, this rotation period is well separated from the orbital periods of the two planets and does not influence our two planet fits.

We do not detect any long term variability in the  values of HD 143761 but we see a clear peak in a periodogram of the  values at 18.5 days that is likely caused by the rotational modulation of star spots (see Figure 9). Since the rotation period is well separated from the periods of either of the planets orbiting HD 143761 this does not affect our Keplerian modeling and is likely absorbed into the stellar jitter term.

Figure 8.— Velocity-activity correlation for HD 42618. A discussion of the chromospheric activity of each of the three stars can be found in Section 4.6. Top: Binned RV time series of the post-upgrade Keck data with planet b. Middle: Binned  time series of the post-upgrade Keck data only. Note the similarities between the variability in the top and middle panels. Bottom: Spearman rank correlation test of the velocities with  values (Spearman 1904). We do not subtract this correlation from the RVs of HD 42618 but instead model the magnetic activity cycle as an additional long-period Keplerian (see Section 4.2).
Figure 9.— Lomb-Scargle periodograms of  chromospheric activity. In each panel the period of the planet announced in this work is marked by the blue dashed line and the power corresponding to an analytical false alarm probability of 1% is marked by the red dotted line (Schwarzenberg-Czerny 1998).  values measured from spectra with SNR40 per pixel or exposure times 25% longer then the median exposure time (due to clouds and/or seeing) can be badly contaminated by the solar spectrum and cause our  extraction pipeline to produce large outliers. These measurements were excluded before calculating the periodograms. No significant periodicity is detected in any of the stars at the orbital periods of the new planets. Top: Periodogram of  values for HD 42618. Middle: Periodogram of  values for HD 143761. Bottom: Periodogram of  values for HD 164922.

5. Photometry

5.1. CoRoT Photometry of HD 42618

HD 42618 was the target of high cadence, high precision, continuous photometric monitoring for 0.5 years with the purpose of detecting solar like oscillations (Baglin et al. 2012). We perform a simple polynomial detrending of the space-based photometry. After removing large ramp-shaped features at the start of two long observing campaigns we then fit an 8th order polynomial to all continuous segments of the data. These segments are 2-20 days in length.

We use the detrended photometry to look for periodic photometric variability that might be caused by rotationally modulated star spots. We detect significant variability with a period of 16.9 days but with a broad distribution of periodogram power around the highest peak (see Figure 10). This is a clear signature of stellar rotation with slightly changing phase and/or differential rotation which creates a broad distribution of increased power in Fourier space near the true rotation period. This star is very similar to the Sun in mass, age, and chemical abundance so the fact that the rotation period is also similar to that of the Sun (26 days) is not surprising. However, we note that the precise location of the highest peak in the Lomb-Scargle periodogram of HD 42618 depends on the polynomial order used to detrend the CoRoT photometry. We also tried high-pass filtering the CoRoT photometry using running median filters with window widths of 20-50 days and only analyzing continuous segments of data longer than 20 days. We found that the period of highest power is somewhat variable but always falls between 12-18 days. Since the periodogram period is dependent on the detrending algorithm we can’t determine the rotation period of HD 42618 precisely, but we estimate that it falls within the range of 12-18 days. We do not detect any significant periodic signal in the RV data near the photometric period.

Figure 10.— CoRoT photometry of HD 42618 discussed in Section 5.1. Top: Detrended light curve. Middle: Lomb-scargle periodogram of the light curve. Bottom: Photometry phase-folded to the period corresponding to the highest peak in the Lomb-Scargle periodogram (16.9 days). We also bin the photometry with bin widths of 0.05 units of phase (red circles).
Figure 11.— CoRoT photometry of HD 42618 phase-folded to the orbital period of planet b. The transit search for HD 42618 b is discussed in Section 5.1. Top: Photometry over the full orbital phase of planet b. The red circles are binned photometric measurements with bin widths of 0.04 units of orbital phase. The grey shaded region shows the 1 uncertainty on the time of inferior conjunction derived from the RV modeling. Bottom: Same as top panel with the x-axis zoomed-in near the time of inferior conjunction. In this panel we only plot the measurements binned with bin widths of 0.002 units of orbital phase. Again, the shaded region represents the 1 uncertainty on the time of inferior conjunction. The black transit model shows the predicted transit depth for a solid iron planet using the mass-radius relation of Weiss & Marcy (2014).

We searched through the detrended CoRoT light curve using the TERRA planet detection algorithm (Petigura et al. 2013a, b). We did not find any periodic box-shaped dimmings with signal-to-noise ratios (SNR) greater than 7. Searches of Kepler photometry commonly require SNR (Jenkins et al. 2010) or SNR (Petigura et al. 2013b), though the SNR threshold depends on the noise structure of the photometry. We conclude that there are no transiting planets having periods between 0.5 and 60 days with transits that are detectable above Poisson, stellar, and instrumental noise. Given the photometric noise properties of HD 42618, we can rule out planets with transits deeper than 150 ppm (1.3 ) at 5 day orbital periods and transits deeper than 300 ppm (1.9 ) for 50 day orbital periods. The a priori transit probability for HD 42618 b is only 0.8% so it is not surprising that we do not detect transits.

5.1.1 Asteroseismic Mass Determination

Convection in the outer layers of a star excites stochastic oscillations, which can be observed on the stellar surface. In the case of main sequence stars, these oscillations manifest themselves as periodic variations on the order of cm/s in radial velocity data or ppm in photometric data. Photometric space telescopes such as CoRoT proved to be quite effective for measuring and characterizing these oscillations, which can be used to derive global stellar properties (such as radius, mass and age), as well as to constrain the stellar interior (e.g. Michel et al. 2008; Chaplin & Miglio 2013).

We measured the mass of HD42618 from CoRoT photometry obtained during two long observing runs spanning 79 and 94 days, respectively. Through a Fourier analysis of the CoRoT lightcurve, we produced the power spectral density function shown in Figure 12. We then stacked the power spectrum in equally sized pieces to create an echelle diagram, revealing the distinct =0, 1, and 2 latitudinal modes of oscillation. We then collapsed this echelle diagram, effectively creating a binned power spectrum, and fit the power excess with a Gaussian to measure a maximum oscillation power frequency of 3.16 0.10 mHz. We then collapsed the echelle diagram along the perpendicular axis to preserve the frequency spacing, computed its autocorrelation, and fit the autocorrelation with a Gaussian to measure a large oscillation frequency spacing of 141.6 0.8 Hz. Using scaling relations (Christensen-Dalsgaard & Frandsen 1983; Kjeldsen & Bedding 1995; Kallinger et al. 2010; Huber et al. 2011), solar parameters taken from Huber et al. (2011), and an effective temperature equivalent to the Sun’s within errors (Morel et al. 2013b) we measure an asteroseismic radius of 0.95 0.05 R, and an asteroseismic mass of 0.93 0.13 M. This is in agreement with a previous CoRoT asteroseismic analysis (Barban et al. 2013) and our estimate of the stellar mass and radius of HD42618 using our spectroscopic constraints and the Torres et al. (2010) relations. The precision on the asteroseismic mass is lower compared to our spectroscopic+isochrone mass but it is much less model-dependent. If we were to fit the spectroscopic parameters to isochrones derived using different input physics we may find that the error on the spectroscopic mass is much larger. We adopt the higher precision, spectroscopic mass for all calculations of planet minimum masses and orbital separations.

Figure 12.— Smoothed one-dimensional power spectrum of HD 42618 from the CoRoT data. The comb of peaks in the power spectrum near a frequency of 3 mHz is the signature of solar-like asteroseismic oscillations. Our asteroseismic analysis of HD 42618 is described in Section 5.1.1.

5.2. APT Photometry

Long-term photometric observations of HD 42618, HD 143761, and HD 164922 were collected with Tennessee State University’s T11 0.80 m, T4 0.75 m, and T12 0.80 m APTs at Fairborn Observatory. These three stars are among a collection of more than 300 being observed by the APTs to study magnetic cycles in solar-type stars (e.g., Lockwood et al. (2013) and references therein) and have APT observational histories between 15 and 23 years. At the beginning of the APF survey, the vast majority of the target stars were already being observed by the APTs. The remaining few have been added so that all 51 stars in the APF survey are also being observed nightly by the APTs.

The APTs are equipped with two-channel precision photometers that use a dichroic filter and two EMI 9124QB bi-alkali photomultiplier tubes to measure the Strömgren and pass bands simultaneously. The APTs are programmed to make differential brightness measurements of a program star with respect to three comparison stars. For the APF project, we use the two best comparison stars (C1 and C2) and compute the differential magnitudes , , and , correct them for atmospheric extinction, and transform them to the Strömgren system. To maximize the precision of the nightly observations, we combine the differential and observations into a single “passband” and also compute the differential magnitudes of the program star against the mean brightness of the two comparison stars. The resulting precision of the individual differential magnitudes ranges between mag and mag on good nights. Further details of our automatic telescopes, precision photometers, and observing and data reduction procedures can be found in Henry (1999), Eaton et al. (2003), and Henry et al. (2013).

Figure 13.— Long-term photometric observations of the planetary candidate host stars HD 42618 (top), HD 143761 (middle), and HD 164922 (bottom) acquired with TSU’s T11 0.80 m, T4 0.75 m, and T12 0.80 m APTs at Fairborn Observatory in southern Arizona. All three stars are plotted with identical x and y scales. The horizonal line in each panel marks the mean of each data set. The APT photometry and analysis are described in Section 5.2.

5.2.1 APT Photometry of HD 42618

We collected 2241 relative flux measurements of HD 42618 over the past 15 years. We search for photometric variability on short timescales by first subtracting the mean magnitude from each observing season to remove seasonal offsets. This removes all astrophysical and systematic instrumental variability on timescales longer then one year. A L-S period search returns a very weak periodicity with a period of 16.5 days and an amplitude of 0.3 mmag. This may be the same signature of stellar rotation as detected in the CoRoT data but it is too close to the precision limit of the ground-based dataset to be certain.

Photometric variability on the timescale of the orbital period may indicate that the RV fluctuations are the result of rotational modulation of star spots (Queloz et al. 2001). We find no evidence of photometric variability at the orbital period of the planet to the limit of our photometric precision. A least-squares sine fit on the orbital period of HD 42618 gives a semi-amplitude of just 0.000037 mag, showing the complete absence of any surface activity that could affect the radial velocities. Figure 13 shows the full photometric dataset and Figure 14 shows the photometry phase folded to the orbital period of planet b. The lack of variability at the orbital period is consistent with the results of the CoRoT analysis and strengthens our claim that the RV fluctuations are caused by a Neptune-mass planet orbiting HD 42618. We also find no evidence of the transit of HD 42618 b in the APT data.

Figure 14.— Top: Fifteen years of photometric observations of HD 42618 from the top panel of Fig. 13 plotted against the 149.6-day planetary orbital period and time of conjunction derived from the radial velocity observations. A least-squares sine fit on the radial velocity period gives a semi-amplitude of just mag, firmly establishing the lack of stellar activity on the radial velocity period and thus confirming the presence of stellar reflex motion caused by an orbiting planet. Bottom: Closeup of the observations near the time of planetary conjunction at phase 0.0. The solid line shows a toy model transit of a sphere of constant 1.0 g cm density and radius determined by the relation of Weiss & Marcy (2014). The vertical lines mark the uncertainty in the predicted transit times. Our current photometric observations provide no evidence for transits.
Figure 15.— Top: Yearly means of the Mt. Wilson-calibrated S values acquired along with the radial velocity measurements. Bottom three panels: Yearly means of HD 42618’s , , and differential magnitudes. The horizontal dotted lines designate the grand means of the observations while the numbers in the lower-left and lower-right give the total range and standard deviation of each data set, respectively. It is evident that we have resolved low-level brightness variability in HD 42618 compared to the two comparison stars, and . Low-amplitude cycles of roughly 0.001 mag over 5 years are seen in both the and light curves. There appears to be little or no correlation of S values with photometric brightness.

Figure 15 shows the mean  values from Keck and APF and the long-term photometric variability of HD 42618 by plotting the seasonal means of both the S values and APT photometry. For old solar type stars we expect a positive correlation of chromospheric activity as measured by the  values with the mean brightness since the number of bright faculae regions on the star increases during more active periods. However, in this case we see no correlation of mean brightness with . Young stars typically show a negative correlation of brightness with  because their photometric variations are spot-dominated instead of faculae dominated. While somewhat unusual, this behavior is not unprecedented among similar stars (Hall et al. 2009).

5.2.2 APT Photometry of HD 164922

We collected a total of 1095 photometric measurements for HD 164922 over the past 11 observing seasons from 2005 to 2015. As with HD 42618 we remove seasonal offsets from the photometry to search for short period variability and search for transits of HD164922 b and c. We find no significant periodic variability with a period between 1 and 100 days and do not detect the rotation period of the star. We also find no evidence of transits for either planet b or c or periodic photometric variability at the orbital period of either planet (Figure 16).

Figure 16.— Same as Figure 14 but for HD 164922 b and c.

We study the long-term photometric variability of HD 164922 by comparing the mean brightness of the star to the  activity index (Figure 17). In contrast to the results for HD 42618, in this case we see a clear positive correlation of the brightness of HD 164922 with the  index. It is interesting that we do not see a RV vs.  correlation for HD 164922, but we do find that the RV is strongly correlated with  for HD 42618 where the photometry is not. In other cases we have seen a correlation in both the photometry and RV data (e.g. Fulton et al. 2015).

Figure 17.— Same as Figure 15 but for HD 164922. In this case we see a positive correlation of the brightness of HD 164922 with the  index.

5.2.3 APT Photometry of HD 143761

We collected 1586 photometric measurements of HD 143761 over the past 18 observing seasons from 1997 to 2015 (Figure 13). Our reduction and analysis techniques are the same as for HD 42618 and HD 164922 discussed in the previous two sections. We find no evidence of the photometric signature of rotationally modulated star spots or photometric variability at the orbital periods of HD 143761 b or c. There is no evidence of transits of either planet b or c (Figure 18), however shallow transits of a rocky planet c can not be ruled out by this dataset.

Figure 18.— Same as Figure 14 but for HD 143761 b and c.

The mean photometric brightness binned by observing season is well correlated with the  values measured using Keck and APF as expected for an old solar type star. As with HD 164922, we do not see a correlation of  with RV but there is a positive correlation of  with mean brightness (Figure 19).

Figure 19.— Same as Figure 15 but for HD 143761. In this case we again see a positive correlation of the brightness of HD 143761 with the  index.

6. Discussion & Summary

We present the discovery of three approximately Neptune mass planets orbiting three bright, nearby stars. The planet orbiting HD 42618 has a minimum mass of  =  and is the first discovered to orbit this star. There has been some discussion in the literature that stellar abundance patterns similar to the Sun might be evidence of the formation of terrestrial planets similar to those that exist in our solar system (e.g. Meléndez et al. 2009; González Hernández et al. 2010). While we can not rule out the existence of terrestrial planets in the HD 42618 system, the presence of a temperate Neptune mass planet orbiting at 0.554 AU with an orbital period of 149 days shows that, at the present time, this system is not a close analogue to our own solar system. We cannot determine the initial planetary architecture of this system because migration may have played an important role to sculpt the current configuration.

We also detect the signature of the stellar magnetic activity cycle with a period of 12 years. This activity cycle manifests as a 3.1 m samplitude signal in the RV time series. We identify the rotation period of the star to be 17 days using public data from the CoRoT space telescope. Transits of HD 42618 b are expected to be extremely unlikely and we do not find any evidence for transits of this planet in the CoRoT data. This is a temperate planet receiving only 3.1 times the radiation that the Earth receives from the Sun. The planet’s equilibrium temperature, assuming a bond albedo of 0.32 (Demory 2014), is 337 K. We perform an asteroseismic study of HD 42618 to detect solar like oscillations and measure a precise stellar radius and mass.

HD 164922 c is the second planet in a system previously known to host one Jupiter mass planet orbiting at 2.1 AU. The new planet announced in this work is a sub-Neptune mass planet with  =   orbiting at a distance of AU and an orbital period of 75 days. This planet is also temperate with an equilibrium temperature of 401 K and receiving 6.3 times the flux received by the Earth from the Sun.

HD 143761 c is the second planet in a system previously known to host a warm Jupiter mass planet orbiting with a period of 39 days. The new planet is a super-Neptune with  =   orbiting with a period of 102 days. This planet is the warmest of the three with a stellar irradiance 9.6 times that of the Earth-Sun system and an equilibrium temperature of 445 K. We find that the previous low inclination orbit for HD 143761 detected in Hipparcos astrometry can not be stable with the presence of HD 143761 c.

These three planets are some of the nearest long period Neptune mass planets yet discovered. They demonstrate the capabilities of the combined KeckAPF-50 survey and are the beginning of a complete census of small planets in the local neighborhood.

Facilities: Automated Planet Finder (Levy), Keck:I (HIRES), CoRoT

We thank the many observers who contributed to the measurements reported here. We thank Kyle Lanclos, Matt Radovan, Will Deich and the rest of the UCO Lick staff for their invaluable help shepherding, planning, and executing observations, in addition to writing the low-level software that made the automated APF observations possible. We thank Gail Schaefer for her help with the calculations related to the CHARA interferometric observations. We thank Debra Fischer, Jason Wright, and John Johnson for their many nights of observing that contributed to the Keck data presented in this work. We gratefully acknowledge the efforts and dedication of the Keck Observatory staff, especially Scott Dahm, Greg Doppman, Hien Tran, and Grant Hill for support of HIRES and Greg Wirth for support of remote observing. We are grateful to the time assignment committees of the University of Hawai‘i, the University of California, and NASA for their generous allocations of observing time. Without their long-term commitment to RV monitoring, these planets would likely remain unknown. We acknowledge R. Paul Butler and S. S. Vogt for many years of contributing to the data presented here. A. W. H. acknowledges NSF grant AST-1517655 and NASA grant NNX12AJ23G. L. M. W. gratefully acknowledges support from Ken and Gloria Levy. D.H. acknowledges support by the Australian Research Council’s Discovery Projects funding scheme (project number DE140101364) and support by the National Aeronautics and Space Administration under Grant NNX14AB92G issued through the Kepler Participating Scientist Program. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. 2014184874. Any opinion, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation. 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. J.T.W. and E.B.F. acknowledge support from multiple NASA Keck PI Data Awards, administered by the NASA Exoplanet Science Institute, to follow multiple exoplanets systems including HD 164922 and HD 143761 from semester 2010A to through 2012B (semester 2010B excluded). J.T.W. acknowledges support from NSF grant AST-1211441 G. W. H. acknowledges support from NASA, NSF, Tennessee State University, and the State of Tennessee through its Centers of Excellence program. This research made use of the Exoplanet Orbit Database and the Exoplanet Data Explorer at exoplanets.org. Based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnolog a e Innovaci n Productiva (Argentina), and Minist rio da Ci ncia, Tecnologia e Inova o (Brazil). This work made use of the SIMBAD database (operated at CDS, Strasbourg, France), and NASA’s Astrophysics Data System Bibliographic Services. Finally, the authors wish to extend special thanks to those of Hawai‘ian ancestry on whose sacred mountain of Maunakea we are privileged to be guests. Without their generous hospitality, the Keck observations presented herein would not have been possible.

References

  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
  • Baglin et al. (2009) Baglin, A., Auvergne, M., Barge, P., Deleuil, M., Michel, E., & CoRoT Exoplanet Science Team. 2009, in IAU Symposium, Vol. 253, IAU Symposium, ed. F. Pont, D. Sasselov, & M. J. Holman, 71–81
  • Baglin et al. (2012) Baglin, A., Michel, E., & CoRoT Team. 2012, in Astronomical Society of the Pacific Conference Series, Vol. 462, Progress in Solar/Stellar Physics with Helio- and Asteroseismology, ed. H. Shibahashi, M. Takata, & A. E. Lynas-Gray, 492
  • Barban et al. (2013) Barban, C., et al. 2013, Journal of Physics Conference Series, 440, 012030
  • Bonneau et al. (2011) Bonneau, D., Delfosse, X., Mourard, D., Lafrasse, S., Mella, G., Cetre, S., Clausse, J.-M., & Zins, G. 2011, A&A, 535, A53
  • Bonneau et al. (2006) Bonneau, D., et al. 2006, A&A, 456, 789
  • Boyajian et al. (2008) Boyajian, T. S., et al. 2008, ApJ, 683, 424
  • Butler et al. (1996) Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C., Dosanjh, P., & Vogt, S. S. 1996, PASP, 108, 500
  • Butler et al. (2006) Butler, R. P., et al. 2006, ApJ, 646, 505
  • Chaplin & Miglio (2013) Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353
  • Chatterjee & Tan (2014) Chatterjee, S., & Tan, J. C. 2014, ApJ, 780, 53
  • Chiang & Laughlin (2013) Chiang, E., & Laughlin, G. 2013, MNRAS, 431, 3444
  • Christensen-Dalsgaard & Frandsen (1983) Christensen-Dalsgaard, J., & Frandsen, S. 1983, Sol. Phys., 82, 469
  • Claret & Bloemen (2011) Claret, A., & Bloemen, S. 2011, A&A, 529, A75
  • Coelho (2014) Coelho, P. R. T. 2014, MNRAS, 440, 1027
  • Cumming et al. (2008) Cumming, A., Butler, R. P., Marcy, G. W., Vogt, S. S., Wright, J. T., & Fischer, D. A. 2008, PASP, 120, 531
  • Cutri et al. (2003) Cutri, R. M., et al. 2003, VizieR Online Data Catalog, 2246, 0
  • Demory (2014) Demory, B.-O. 2014, The Astrophysical Journal, 789, L20
  • Dotter et al. (2008) Dotter, A., Chaboyer, B., Jevremović, D., Kostov, V., Baron, E., & Ferguson, J. W. 2008, ApJS, 178, 89
  • Dumusque et al. (2011) Dumusque, X., Udry, S., Lovis, C., Santos, N. C., & Monteiro, M. J. P. F. G. 2011, A&A, 525, A140
  • Eastman et al. (2013) Eastman, J., Gaudi, B. S., & Agol, E. 2013, PASP, 125, 83
  • Eaton et al. (2003) Eaton, J. A., Henry, G. W., & Fekel, F. C. 2003, in Astrophysics and Space Science Library, Vol. 288, Astrophysics and Space Science Library, ed. T. D. Oswalt, 189
  • Fabrycky et al. (2014) Fabrycky, D. C., et al. 2014, ApJ, 790, 146
  • Farrington et al. (2010) Farrington, C. D., et al. 2010, AJ, 139, 2308
  • Fischer et al. (2016) Fischer, D. A., et al. 2016, PASP, 128, 066001
  • Ford (2006) Ford, E. B. 2006, ApJ, 642, 505
  • Fressin et al. (2013) Fressin, F., et al. 2013, ApJ, 766, 81
  • Fuhrmann et al. (1998) Fuhrmann, K., Pfeiffer, M. J., & Bernkopf, J. 1998, A&A, 336, 942
  • Fulton et al. (2013) Fulton, B. J., et al. 2013, ApJ, 772, 80
  • Fulton et al. (2015) —. 2015, ApJ, 805, 175
  • Gatewood et al. (2001) Gatewood, G., Han, I., & Black, D. C. 2001, ApJ, 548, L61
  • Ghezzi et al. (2010) Ghezzi, L., Cunha, K., Smith, V. V., de Araújo, F. X., Schuler, S. C., & de la Reza, R. 2010, ApJ, 720, 1290
  • Gladman (1993) Gladman, B. 1993, Icarus, 106, 247
  • González Hernández et al. (2010) González Hernández, J. I., Israelian, G., Santos, N. C., Sousa, S., Delgado-Mena, E., Neves, V., & Udry, S. 2010, ApJ, 720, 1592
  • Gray et al. (2003) Gray, R. O., Corbally, C. J., Garrison, R. F., McFadden, M. T., & Robinson, P. E. 2003, AJ, 126, 2048
  • Hall et al. (2009) Hall, J. C., Henry, G. W., Lockwood, G. W., Skiff, B. A., & Saar, S. H. 2009, AJ, 138, 312
  • Han et al. (2014) Han, E., Wang, S. X., Wright, J. T., Feng, Y. K., Zhao, M., Fakhouri, O., Brown, J. I., & Hancock, C. 2014, PASP, 126, 827
  • Hansen & Murray (2012) Hansen, B. M. S., & Murray, N. 2012, ApJ, 751, 158
  • Henry (1999) Henry, G. W. 1999, PASP, 111, 845
  • Henry et al. (2013) Henry, G. W., et al. 2013, ApJ, 768, 155
  • Horch et al. (2011) Horch, E. P., Gomez, S. C., Sherry, W. H., Howell, S. B., Ciardi, D. R., Anderson, L. M., & van Altena, W. F. 2011, AJ, 141, 45
  • Horch et al. (2012) Horch, E. P., Howell, S. B., Everett, M. E., & Ciardi, D. R. 2012, AJ, 144, 165
  • Horch et al. (2009) Horch, E. P., Veillette, D. R., Baena Gallé, R., Shah, S. C., O’Rielly, G. V., & van Altena, W. F. 2009, AJ, 137, 5057
  • Howard & Fulton (2016) Howard, A. W., & Fulton, B. J. 2016, ArXiv e-prints
  • Howard et al. (2009) Howard, A. W., et al. 2009, ApJ, 696, 75
  • Howard et al. (2010) Howard, A. W., et al. 2010, Science, 330, 653
  • Howard et al. (2011a) Howard, A. W., et al. 2011a, ApJ, 726, 73
  • Howard et al. (2011b) —. 2011b, ApJ, 730, 10
  • Howard et al. (2012) —. 2012, ApJS, 201, 15
  • Howard et al. (2014) —. 2014, ApJ, 794, 51
  • Howell et al. (2011) Howell, S. B., Everett, M. E., Sherry, W., Horch, E., & Ciardi, D. R. 2011, AJ, 142, 19
  • Huber et al. (2011) Huber, D., et al. 2011, ApJ, 743, 143
  • Ida & Lin (2004) Ida, S., & Lin, D. N. C. 2004, ApJ, 604, 388
  • Ida & Lin (2008) —. 2008, ApJ, 673, 487
  • Ireland et al. (2008) Ireland, M. J., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7013, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 24
  • Isaacson & Fischer (2010) Isaacson, H., & Fischer, D. 2010, The Astrophysical Journal, 725, 875
  • Jenkins et al. (2010) Jenkins, J. M., et al. 2010, in Proc. SPIE, Vol. 7740, Software and Cyberinfrastructure for Astronomy, 77400D
  • Johnson et al. (2011) Johnson, J. A., et al. 2011, ApJS, 197, 26
  • Kallinger et al. (2010) Kallinger, T., et al. 2010, A&A, 509, A77
  • Kane et al. (2014) Kane, S. R., et al. 2014, ApJ, 785, 93
  • Kjeldsen & Bedding (1995) Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87
  • Knutson et al. (2014) Knutson, H. A., et al. 2014, ApJ, 785, 126
  • Koen et al. (2010) Koen, C., Kilkenny, D., van Wyk, F., & Marang, F. 2010, MNRAS, 403, 1949
  • Lee & Chiang (2015) Lee, E. J., & Chiang, E. 2015, ArXiv e-prints
  • Lee et al. (2014) Lee, E. J., Chiang, E., & Ormel, C. W. 2014, ApJ, 797, 95
  • Lockwood et al. (2013) Lockwood, G. W., Henry, G. W., Hall, J. C., & Radick, R. R. 2013, in Astronomical Society of the Pacific Conference Series, Vol. 472, New Quests in Stellar Astrophysics III: A Panchromatic View of Solar-Like Stars, With and Without Planets, ed. M. Chavez, E. Bertone, O. Vega, & V. De la Luz, 203
  • Lomb (1976) Lomb, N. R. 1976, Ap&SS, 39, 447
  • Lopez (2014) Lopez, E. D. 2014, PhD thesis, University of California, Santa Cruz
  • Lopez & Fortney (2014) Lopez, E. D., & Fortney, J. J. 2014, ApJ, 792, 1
  • Mann & von Braun (2015) Mann, A. W., & von Braun, K. 2015, PASP, 127, 102
  • Mayor et al. (2011) Mayor, M., et al. 2011, arXiv:1109.2497
  • Medhi et al. (2007) Medhi, B. J., Messina, S., Parihar, P. S., Pagano, I., Muneer, S., & Duorah, K. 2007, A&A, 469, 713
  • Meléndez et al. (2009) Meléndez, J., Asplund, M., Gustafsson, B., & Yong, D. 2009, ApJ, 704, L66
  • Michel et al. (2008) Michel, E., et al. 2008, Science, 322, 558
  • Mordasini et al. (2009) Mordasini, C., Alibert, Y., & Benz, W. 2009, A&A, 501, 1139
  • Morel et al. (2013a) Morel, T., Rainer, M., Poretti, E., Barban, C., & Boumier, P. 2013a, A&A, 552, A42
  • Morel et al. (2013b) —. 2013b, A&A, 552, A42
  • Morton (2015) Morton, T. D. 2015, isochrones: Stellar model grid package, Astrophysics Source Code Library
  • Motalebi et al. (2015) Motalebi, F., et al. 2015, A&A, 584, A72
  • Noyes et al. (1997) Noyes, R. W., Jha, S., Korzennik, S. G., Krockenberger, M., Nisenson, P., Brown, T. M., Kennelly, E. J., & Horner, S. D. 1997, ApJ, 483, L111
  • O’Toole et al. (2009) O’Toole, S. J., Jones, H. R. A., Tinney, C. G., Butler, R. P., Marcy, G. W., Carter, B., Bailey, J., & Wittenmyer, R. A. 2009, ApJ, 701, 1732
  • Pepe et al. (2004) Pepe, F., et al. 2004, A&A, 423, 385
  • Petigura (2015) Petigura, E. A. 2015, arXiv:1510.03902
  • Petigura et al. (2013a) Petigura, E. A., Howard, A. W., & Marcy, G. W. 2013a, Proceedings of the National Academy of Science, 110, 19273
  • Petigura et al. (2013b) Petigura, E. A., Marcy, G. W., & Howard, A. W. 2013b, ApJ, 770, 69
  • Pickles (1998) Pickles, A. J. 1998, PASP, 110, 863
  • Queloz et al. (2001) Queloz, D., et al. 2001, A&A, 379, 279
  • Radovan et al. (2014) Radovan, M. V., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9145, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 2
  • Raghavan et al. (2012) Raghavan, D., Farrington, C. D., ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., Sturmann, L., Sturmann, J., & Turner, N. H. 2012, ApJ, 745, 24
  • Ramírez et al. (2012) Ramírez, I., Fish, J. R., Lambert, D. L., & Allende Prieto, C. 2012, ApJ, 756, 46
  • Ramírez et al. (2014) Ramírez, I., et al. 2014, A&A, 572, A48
  • Reffert & Quirrenbach (2011) Reffert, S., & Quirrenbach, A. 2011, A&A, 527, A140
  • Santos et al. (2013) Santos, N. C., et al. 2013, A&A, 556, A150
  • Scargle (1982) Scargle, J. D. 1982, ApJ, 263, 835
  • Schwarzenberg-Czerny (1998) Schwarzenberg-Czerny, A. 1998, MNRAS, 301, 831
  • Spearman (1904) Spearman, C. 1904, The American Journal of Psychology, 15, pp. 72
  • Swift et al. (2015) Swift, J. J., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 027002
  • ten Brummelaar et al. (2005) ten Brummelaar, T. A., et al. 2005, ApJ, 628, 453
  • Ter Braak (2006) Ter Braak, C. 2006, Statistics and Computing
  • Torres et al. (2010) Torres, G., Andersen, J., & Giménez, A. 2010, A&A Rev., 18, 67
  • Valenti & Fischer (2005) Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141
  • van Belle & von Braun (2009) van Belle, G. T., & von Braun, K. 2009, ApJ, 694, 1085
  • van Leeuwen (2007) van Leeuwen, F. 2007, A&A, 474, 653
  • Veras & Ford (2012) Veras, D., & Ford, E. B. 2012, MNRAS, 420, L23
  • Vogt et al. (1994) Vogt, S. S., et al. 1994, in Proc. SPIE Instrumentation in Astronomy VIII, David L. Crawford; Eric R. Craine; Eds., 2198, 362
  • Vogt et al. (2014) Vogt, S. S., et al. 2014, PASP, 126, 359
  • von Braun et al. (2014) von Braun, K., et al. 2014, MNRAS, 438, 2413
  • Weiss & Marcy (2014) Weiss, L. M., & Marcy, G. W. 2014, ApJ, 783, L6
  • White et al. (2013) White, T. R., et al. 2013, MNRAS, 433, 1262
  • Wright & Eastman (2014) Wright, J. T., & Eastman, J. D. 2014, PASP, 126, 838
  • Wright et al. (2007) Wright, J. T., et al. 2007, ApJ, 657, 533
  • Zucker & Mazeh (2001) Zucker, S., & Mazeh, T. 2001, arXiv, astro-ph/0104098
\enddocument
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
116123
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description