We confirm the planetary nature of TOI-1728b using a combination of ground-based photometry, near-infrared Doppler velocimetry and spectroscopy with the Habitable-zone Planet Finder.TOI-1728 is an old, inactive M0 star with \teff{} $= 3980^{+31}_{-32}$ K, which hosts a transiting super Neptune at an orbital period of $\sim$ 3.49 days. Joint fitting of the radial velocities and TESS and ground-based transits yields a planetary radius of $5.05_{-0.17}^{+0.16}$ R$_{\oplus}$, mass $26.78_{-5.13}^{+5.43}$ M$_{\oplus}$ and eccentricity $0.057_{-0.039}^{+0.054}$. We estimate the stellar properties, and perform a search for He 10830 \AA absorption during the transit of this planet and claim a null detection with an upper limit of 1.1$\%$ with 90\% confidence. A deeper level of He 10830 \AA ~ absorption has been detected in the planet atmosphere of GJ 3470b, a comparable gaseous planet. TOI-1728b is the largest super Neptune -- the intermediate subclass of planets between Neptune and the more massive gas-giant planets -- discovered around an M dwarf. With its relatively large mass and radius, TOI-1728 represents a valuable datapoint in the M-dwarf exoplanet mass-radius diagram, bridging the gap between the lighter Neptune-sized planets and the heavier Jovian planets known to orbit M-dwarfs. With a low bulk density of $1.14_{-0.24}^{+0.26}$ g/cm$^3$, and orbiting a bright host star (J $\sim 9.6$, V $\sim 12.4$), TOI-1728b is also a promising candidate for transmission spectroscopy both from the ground and from space, which can be used to constrain planet formation and evolutionary models.

We report on the validation of two planets orbiting the nearby (36pc) M2 dwarf TOI-1266 observed by the TESS mission. The inner planet is sub-Neptune-sized ($R=2.46 \pm 0.08 R_\oplus$) with an orbital period of 10.9 days. The outer planet has a radius of $1.67_{-0.11}^{+0.09} R_\oplus$ and resides in the exoplanet Radius Valley---the transition region between rocky and gaseous planets. With an orbital period of 18.8 days, the outer planet receives an insolation flux of 2.4 times that of Earth, similar to the insolation of Venus. Using precision near-infrared radial velocities with the Habitable-zone Planet Finder Spectrograph, we place upper mass limits of $15.9 M_\oplus$ and $6.4 M_\oplus$ at 95% confidence for the inner and outer planet, respectively. A more precise mass constraint of planet c, achievable with current RV instruments given the host star brightness (V=12.9, J=9.7), will yield further insights into the dominant processes sculpting the exoplanet Radius Valley.

Young, rapidly-rotating M dwarfs exhibit prominent starspots, which create quasiperiodic signals in their photometric and Doppler spectroscopic measurements. The periodic Doppler signals can mimic radial velocity (RV) changes expected from orbiting exoplanets. Exoplanets can be distinguished from activity-induced false positives by the chromaticity and long-term incoherence of starspot signals, but these qualities are poorly constrained for fully-convective M stars. Coherent photometric starspot signals on M dwarfs may persist for hundreds of rotations, and the wavelength dependence of starspot RV signals may not be consistent between stars due to differences in their magnetic fields and active regions. We obtained precise multi-wavelength RVs of four rapidly-rotating M dwarfs (AD Leo, G 227-22, GJ 1245B, GJ 3959) using the near-infrared (NIR) Habitable-zone Planet Finder, and the optical Keck/HIRES spectrometer. Our RVs are complemented by photometry from Kepler, TESS, and the Las Cumbres Observatory (LCO) network of telescopes. We found that all four stars exhibit large spot-induced Doppler signals at their rotation periods, and investigated the longevity and optical-to-NIR chromaticity for these signals. The phase curves remain coherent much longer than is typical for Sunlike stars. Their chromaticity varies, and one star (GJ 3959) exhibits optical and NIR RV modulation consistent in both phase and amplitude. In general, though, we find that the NIR amplitudes are lower than their optical counterparts. We conclude that starspot modulation for rapidly-rotating M stars frequently remains coherent for hundreds of stellar rotations, and gives rise to Doppler signals that, due to this coherence, may be mistaken for exoplanets.

Understanding the dynamics and kinematics of out-flowing atmospheres of hot and warm exoplanets is crucial to understanding the origins and evolutionary history of the exoplanets near the evaporation desert. Recently, ground based measurements of the meta-stable Helium atom's resonant absorption at 10830 \AA~has become a powerful probe of the base environment which is driving the outflow of exoplanet atmospheres. We report evidence for the He I 10830 \AA~in absorption (equivalent width $\sim$ $0.012 \pm 0.002$ \AA) in the exosphere of a warm Neptune orbiting the M-dwarf GJ 3470, during three transits using the Habitable Zone Planet Finder (HPF) near infrared spectrograph. This marks the first reported evidence for He I 10830 \AA\, atmospheric absorption for a planet orbiting an M-dwarf. Our detected absorption is broad and its blueshifted wing extends to -36 km/sec, the largest reported in the literature to date. We modelled the state of Helium atoms in the exosphere of GJ3470b based on assumptions on the UV and X-ray flux of GJ 3470, and found our measurement of flux-weighted column density of meta-stable state Helium $(N_{He^2_3S} = 2.4 \times 10^{10} \mathrm{cm^{-2}})$, derived from our transit observations, to be consistent with model, within its uncertainties. The methodology developed here will be useful to study and constrain the atmospheric outflow models of other exoplanets like GJ 3470b which are near the edge of the evaporation desert.

Solar contamination, due to moonlight and atmospheric scattering of sunlight, can cause systematic errors in stellar radial velocity (RV) measurements that significantly detract from the ~10cm/s sensitivity required for the detection and characterization of terrestrial exoplanets in or near Habitable Zones of Sun-like stars. The addition of low-level spectral contamination at variable effective velocity offsets introduces systematic noise when measuring velocities using classical mask-based or template-based cross-correlation techniques. Here we present simulations estimating the range of RV measurement error induced by uncorrected scattered sunlight contamination. We explore potential correction techniques, using both simultaneous spectrometer sky fibers and broadband imaging via coherent fiber imaging bundles, that could reliably reduce this source of error to below the photon-noise limit of typical stellar observations. We discuss the limitations of these simulations, the underlying assumptions, and mitigation mechanisms. We also present and discuss the components designed and built into the NEID precision RV instrument for the WIYN 3.5m telescope, to serve as an ongoing resource for the community to explore and evaluate correction techniques. We emphasize that while "bright time" has been traditionally adequate for RV science, the goal of 10cm/s precision on the most interesting exoplanetary systems may necessitate access to darker skies for these next-generation instruments.

We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distance of $d=27.9\mathrm{pc}$, G 9-40b is the second closest transiting planet discovered by K2 to date. The planet's large transit depth ($\sim$3500ppm), combined with the proximity and brightness of the host star at NIR wavelengths (J=10, K=9.2) makes G 9-40b one of the most favorable sub-Neptune-sized planet orbiting an M-dwarf for transmission spectroscopy with JWST, ARIEL, and the upcoming Extremely Large Telescopes. The star is relatively inactive with a rotation period of $\sim$29 days determined from the K2 photometry. To estimate spectroscopic stellar parameters, we describe our implementation of an empirical spectral matching algorithm using the high-resolution NIR HPF spectra. Using this algorithm, we obtain an effective temperature of $T_{\mathrm{eff}}=3404\pm73$K, and metallicity of $\mathrm{[Fe/H]}=-0.08\pm0.13$. Our RVs, when coupled with the orbital parameters derived from the transit photometry, exclude planet masses above $11.7 M_\oplus$ with 99.7% confidence assuming a circular orbit. From its radius, we predict a mass of $M=5.0^{+3.8}_{-1.9} M_\oplus$ and an RV semi-amplitude of $K=4.1^{+3.1}_{-1.6}\mathrm{m\:s^{-1}}$, making its mass measurable with current RV facilities. We urge further RV follow-up observations to precisely measure its mass, to enable precise transmission spectroscopic measurements in the future.

M dwarfs, though the most abundant star in the galaxy, form only a small subset of stellar hosts with exoplanets with measured radii and masses. In this paper we analyze the Mass-Radius (M-R) relationship of planets around M dwarfs using M-R measurements for 24 exoplanets. In particular, we apply both parametric and nonparametric models and compare the two different fitting methods. We also use these methods to compare the results of the M dwarf M-R relationship with that from the Kepler sample. Using the nonparametric method, we find that the predicted masses for the smallest and largest planets around M dwarfs are smaller than that given by the relation fit to the Kepler data, but that the distribution of masses for 3 Earth Radii planets does not substantially differ between the two datasets. With future additions to the M dwarf M-R relation from TESS and instruments like the Habitable Zone Planet Finder, we will be able to characterize these differences in more detail, which will help illuminate the process of planetary formation and evolution around these stars. We release a publicly available Python code called MRExo which uses the nonparametric algorithm introduced by Ning et al. (2018) to fit the M-R relationship. Such a nonparametric fit does not assume an underlying power law fit to the measurements and hence can be used on samples spanning large mass and radii ranges. Therefore by not assuming a functional form, the fit is less biased. MRExo also offers a tool to predict mass from radius posteriors, and vice versa. This functionality can help inform observational strategies for radial velocity campaigns, such as TESS follow-up studies, as well as predict radii with microlensing planet masses.

Two key areas of emphasis in contemporary experimental exoplanet science are the detailed characterization of transiting terrestrial planets, and the search for Earth analog planets to be targeted by future imaging missions. Both of these pursuits are dependent on an order-of-magnitude improvement in the measurement of stellar radial velocities (RV), setting a requirement on single-measurement instrumental uncertainty of order 10 cm/s. Achieving such extraordinary precision on a high-resolution spectrometer requires thermo-mechanically stabilizing the instrument to unprecedented levels. Here, we describe the Environment Control System (ECS) of the NEID Spectrometer, which will be commissioned on the 3.5 m WIYN Telescope at Kitt Peak National Observatory in 2019, and has a performance specification of on-sky RV precision < 50 cm/s. Because NEID's optical table and mounts are made from aluminum, which has a high coefficient of thermal expansion, sub-milliKelvin temperature control is especially critical. NEID inherits its ECS from that of the Habitable-zone Planet Finder (HPF), but with modifications for improved performance and operation near room temperature. Our full-system stability test shows the NEID system exceeds the already impressive performance of HPF, maintaining vacuum pressures below $10^{-6}$ Torr and an RMS temperature stability better than 0.4 mK over 30 days. Our ECS design is fully open-source; the design of our temperature-controlled vacuum chamber has already been made public, and here we release the electrical schematics for our custom Temperature Monitoring and Control (TMC) system.

The discovery and characterization of exoplanets around nearby stars is driven by profound scientific questions about the uniqueness of Earth and our Solar System, and the conditions under which life could exist elsewhere in our Galaxy. Doppler spectroscopy, or the radial velocity (RV) technique, has been used extensively to identify hundreds of exoplanets, but with notable challenges in detecting terrestrial mass planets orbiting within habitable zones. We describe infrared RV spectroscopy at the 10 m Hobby-Eberly telescope that leverages a 30 GHz electro-optic laser frequency comb with nanophotonic supercontinuum to calibrate the Habitable Zone Planet Finder spectrograph. Demonstrated instrument precision <10 cm/s and stellar RVs approaching 1 m/s open the path to discovery and confirmation of habitable zone planets around M-dwarfs, the most ubiquitous type of stars in our Galaxy.

Many articulations of the Fermi Paradox have as a premise, implicitly or explicitly, that humanity has searched for signs of extraterrestrial radio transmissions and concluded that there are few or no obvious ones to be found. Tarter et al. (2010) and others have argued strongly to the contrary: bright and obvious radio beacons might be quite common in the sky, but we would not know it yet because our search completeness to date is so low, akin to having searched a drinking glass's worth of seawater for evidence of fish in all of Earth's oceans. Here, we develop the metaphor of the multidimensional "Cosmic Haystack" through which SETI hunts for alien "needles" into a quantitative, eight-dimensional model and perform an analytic integral to compute the fraction of this haystack that several large radio SETI programs have collectively examined. Although this model haystack has many qualitative differences from the Tarter et al. (2010) haystack, we conclude that the fraction of it searched to date is also very small: similar to the ratio of the volume of a large hot tub or small swimming pool to that of the Earth's oceans. With this article we provide a Python script to calculate haystack volumes for future searches and for similar haystacks with different boundaries. We hope this formalism will aid in the development of a common parameter space for the computation of upper limits and completeness fractions of search programs for radio and other technosignatures.

The Habitable-zone Planet Finder (HPF) is a highly stabilized fiber fed precision radial velocity (RV) spectrograph working in the Near Infrared (NIR): 810 - 1280 nm . In this paper we present an overview of the preparation of the optical fibers for HPF. The entire fiber train from the telescope focus down to the cryostat is detailed. We also discuss the fiber polishing, splicing and its integration into the instrument using a fused silica puck. HPF was designed to be able to operate in two modes, High Resolution (HR- the only mode mode currently commissioned) and High Efficiency (HE). We discuss these fiber heads and the procedure we adopted to attach the slit on to the HR fibers.

We announce barycorrpy (BCPy) , a Python implementation to calculate precise barycentric corrections well below the 1 cm/s level, following the algorithm of Wright and Eastman (2014). This level of precision is required in the search for 1 Earth mass planets in the Habitable Zones of Sun-like stars by the Radial Velocity (RV) method, where the maximum semi-amplitude is about 9 cm/s. We have developed BCPy to be used in the pipeline for the next generation Doppler Spectrometers - Habitable-zone Planet Finder (HPF) and NEID. In this work, we also develop an automated leap second management routine to improve upon the one available in Astropy. It checks for and downloads a new leap second file before converting from the UT time scale to TDB.

We demonstrate a path to hitherto unachievable differential photometric precisions from the ground, both in the optical and near-infrared (NIR), using custom-fabricated beam-shaping diffusers produced using specialized nanofabrication techniques. Such diffusers mold the focal plane image of a star into a broad and stable top-hat shape, minimizing photometric errors due to non-uniform pixel response, atmospheric seeing effects, imperfect guiding, and telescope-induced variable aberrations seen in defocusing. This PSF reshaping significantly increases the achievable dynamic range of our observations, increasing our observing efficiency and thus better averages over scintillation. Diffusers work in both collimated and converging beams. We present diffuser-assisted optical observations demonstrating $62^{+26}_{-16}$ppm precision in 30 minute bins on a nearby bright star 16-Cygni A (V=5.95) using the ARC 3.5m telescope---within a factor of $\sim$2 of Kepler's photometric precision on the same star. We also show a transit of WASP-85-Ab (V=11.2) and TRES-3b (V=12.4), where the residuals bin down to $180^{+66}_{-41}$ppm in 30 minute bins for WASP-85-Ab---a factor of $\sim$4 of the precision achieved by the K2 mission on this target---and to 101ppm for TRES-3b. In the NIR, where diffusers may provide even more significant improvements over the current state of the art, our preliminary tests have demonstrated $137^{+64}_{-36}$ppm precision for a $K_S =10.8$ star on the 200" Hale Telescope. These photometric precisions match or surpass the expected photometric precisions of TESS for the same magnitude range. This technology is inexpensive, scalable, easily adaptable, and can have an important and immediate impact on the observations of transits and secondary eclipses of exoplanets.