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.

Hot Jupiters seem to get rarer with decreasing stellar mass. The goal of the Pan-Planets transit survey was the detection of such planets and a statistical characterization of their frequency. Here, we announce the discovery and validation of two planets found in that survey, Wendelstein-1b and Wendelstein-2b, which are two short-period hot Jupiters that orbit late K host stars. We validated them both by the traditional method of radial velocity measurements with the HIgh Resolution Echelle Spectrometer (HIRES) and the Habitable-zone Planet Finder (HPF) instruments and then by their Transit Color Signature (TraCS). We observed the targets in the wavelength range of $4000 - 24000$ Angstr\"om and performed a simultaneous multiband transit fit and additionally determined their thermal emission via secondary eclipse observations. Wendelstein-1b is a hot Jupiter with a radius of $1.0314_{-0.0061}^{+0.0061}$ $R_J$ and mass of $0.592_{-0.129}^{+0.165}$ $M_J$, orbiting a K7V dwarf star at a period of $2.66$ d, and has an estimated surface temperature of about $1727_{-90}^{+78}$ K. Wendelstein-2b is a hot Jupiter with a radius of $1.1592_{-0.0210}^{+0.0204}$ $R_J$ and a mass of $0.731_{-0.311}^{+0.541}$ $M_J$, orbiting a K6V dwarf star at a period of $1.75$ d, and has an estimated surface temperature of about $1852_{-140}^{+120}$ K. With this, we demonstrate that multiband photometry is an effective way of validating transiting exoplanets, in particular for fainter targets since radial velocity (RV) follow-up becomes more and more costly for those targets.

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.

We report the detection of a hot Jupiter ($M_{p}=1.75_{-0.17}^{+0.14}\ M_{J}$, $R_{p}=1.38\pm0.04\ R_{J}$) orbiting a middle-aged star ($\log g=4.152^{+0.030}_{-0.043}$) in the Transiting Exoplanet Survey Satellite (TESS) southern continuous viewing zone ($\beta=-79.59^{\circ}$). We confirm the planetary nature of the candidate TOI-150.01 using radial velocity observations from the APOGEE-2 South spectrograph and the Carnegie Planet Finder Spectrograph, ground-based photometric observations from the robotic Three-hundred MilliMeter Telescope at Las Campanas Observatory, and Gaia distance estimates. Large-scale spectroscopic surveys, such as APOGEE/APOGEE-2, now have sufficient radial velocity precision to directly confirm the signature of giant exoplanets, making such data sets valuable tools in the TESS era. Continual monitoring of TOI-150 by TESS can reveal additional planets and subsequent observations can provide insights into planetary system architectures involving a hot Jupiter around a star about halfway through its main-sequence life.

We present the characterization of two engineered diffusers mounted on the 2.5 meter Nordic Optical Telescope, located at Roque de Los Muchachos, Spain. To assess the reliability and the efficiency of the diffusers, we carried out several test observations of two photometric standard stars, along with observations of one primary transit observation of TrES-3b in the red (R-band), one of CoRoT-1b in the blue (B-band), and three secondary eclipses of WASP-12b in V-band. The achieved photometric precision is in all cases within the sub-millimagnitude level for exposures between 25 and 180 seconds. Along a detailed analysis of the functionality of the diffusers, we add a new transit depth measurement in the blue (B-band) to the already observed transmission spectrum of CoRoT-1b, disfavouring a Rayleigh slope. We also report variability of the eclipse depth of WASP-12b in the V-band. For the WASP-12b secondary eclipses, we observe a secondary-depth deviation of about 5-sigma, and a difference of 6-sigma and 2.5-sigma when compared to the values reported by other authors in similar wavelength range determined from Hubble Space Telescope data. We further speculate about the potential physical processes or causes responsible for this observed variability

Teledyne's H2RG detector images suffer from cross-hatch like patterns which arises from sub-pixel quantum efficiency (QE) variation. In this paper we present our measurements of this sub-pixel QE variation in the Habitable-Zone Planet Finder's H2RG detector. We present a simple model to estimate the impact of sub-pixel QE variations on the radial velocity, and how a first order correction can be implemented to correct for the artifact in the spectrum. We also present how the HPF's future upgraded laser frequency comb will enable us to implement this correction.

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.

The Transiting Exoplanet Survey Satellite (TESS, launched early 2018) is expected to find a multitude of new transiting planet candidates around the nearest and brightest stars. Timely high-precision follow-up observations from the ground are essential in confirming and further characterizing the planet candidates that TESS will find. However, achieving extreme photometric precisions from the ground is challenging, as ground-based telescopes are subject to numerous deleterious atmospheric effects. Beam-shaping diffusers are emerging as a low-cost technology to achieve hitherto unachievable differential photometric precisions from the ground. These 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. In this paper, we expand on our previous work (Stefansson et al. 2017; Stefansson et al. 2018 [submitted]), providing a further detailed discussion of key guidelines when sizing a diffuser for use on a telescope. Furthermore, we present our open source Python package iDiffuse which can calculate the expected PSF size of a diffuser in a telescope system, along with its expected on-sky diffuser-assisted photometric precision for a host star of a given magnitude. We use iDiffuse to show that most ($\sim$80\%) of the planet hosts that TESS will find will be scintillation limited in transit observations from the ground. Although iDiffuse has primarily been developed to plan challenging transit observations using the diffuser on the ARCTIC imager on the ARC 3.5m Telescope at Apache Point observatory, iDiffuse is modular and can be easily extended to calculate the expected diffuser-assisted photometric precisions on other telescopes.

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 present precision transit observations of the Neptune-sized planets K2-28b and K2-100b, using the Engineered Diffuser on the ARCTIC imager on the ARC 3.5m Telescope at Apache Point Observatory. K2-28b is a $R_{p} = 2.56 R_\oplus$ mini-Neptune transiting a bright (J=11.7) metal-rich M4 dwarf, offering compelling prospects for future atmospheric characterization. K2-100b is a $R_{p} = 3.45 R_\oplus$ Neptune in the Praesepe Cluster and is one of few planets known in a cluster transiting a host star bright enough ($V=10.5$) for precision radial velocity observations. Using the precision photometric capabilities of the diffuser/ARCTIC system, allows us to achieve a precision of $105^{+87}_{-37}$ppm, and $38^{+21}_{-11}$ppm in 30 minute bins for K2-28b, and K2-100b, respectively. Our joint-fits to the K2 and ground-based light-curves give an order of magnitude improvement in the orbital ephemeris for both planets, yielding a timing precision of 2min in the JWST era. Although we show that the currently available broad-band measurements of K2-28b's radius are currently too imprecise to place useful constraints on K2-28b's atmosphere, we demonstrate that JWST/NIRISS will be able to discern between a cloudy/clear atmosphere in a modest number of transit observations. Our light-curve of K2-100b marks the first transit follow-up observation of this challenging-to-observe transit, where we obtain a transit depth of $819 \pm 50 \mathrm{ppm}$ in the SDSS $i^\prime$ band. We conclude that diffuser-assisted photometry can play an important role in the TESS era to perform timely and precise follow-up of the expected bounty of TESS planet candidates.

We present a single candidate transit event around Proxima Centauri, found during a blind transit search using a robotic 30\,cm telescope at Las Campanas Observatory. The event lasted 1 hour, with an estimated depth of 5\,mmag, and is inconsistent with the transit window predicted for the recently discovered planet b. We modeled the lightcurve under the assumption that the event was caused by a transiting exoplanet, and our model predicts the planet has a radius $R \sim 1 R_{\oplus}$. We encourage continued monitoring of Proxima to elucidate the origin of this event.

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.

Insufficient instrument thermo-mechanical stability is one of the many roadblocks for achieving 10cm/s Doppler radial velocity (RV) precision, the precision needed to detect Earth-twins orbiting Solar-type stars. Highly temperature and pressure stabilized spectrographs allow us to better calibrate out instrumental drifts, thereby helping in distinguishing instrumental noise from astrophysical stellar signals. We present the design and performance of the Environmental Control System (ECS) for the Habitable-zone Planet Finder (HPF), a high-resolution (R=50,000) fiber-fed near infrared (NIR) spectrograph for the 10m Hobby Eberly Telescope at McDonald Observatory. HPF will operate at 180K, driven by the choice of an H2RG NIR detector array with a 1.7micron cutoff. This ECS has demonstrated 0.6mK RMS stability over 15 days at both 180K and 300K, and maintained high quality vacuum (<$10^{-7}$Torr) over months, during long-term stability tests conducted without a planned passive thermal enclosure surrounding the vacuum chamber. This control scheme is versatile and can be applied as a blueprint to stabilize future NIR and optical high precision Doppler instruments over a wide temperature range from ~77K to elevated room temperatures. A similar ECS is being implemented to stabilize NEID, the NASA/NSF NN-EXPLORE spectrograph for the 3.5m WIYN telescope at Kitt Peak, operating at 300K. A full SolidWorks 3D-CAD model and a comprehensive parts list of the HPF ECS are included with this manuscript to facilitate the adaptation of this versatile environmental control scheme in the broader astronomical community.