Search

We simulate the yield of small (0.5-4.0 R$_\oplus$) transiting exoplanets around single mid-M and ultra-cool dwarfs (UCDs) in the Nancy Grace Roman Space Telescope Galactic Bulge Time Domain Survey. We consider multiple approaches for simulating M3-T9 sources within the survey fields, including scaling local space densities and using Galactic stellar population synthesis models. These approaches independently predict $\sim$100,000 single mid-M dwarfs and UCDs brighter than a Roman F146 magnitude of 21 that are within the survey fields. Assuming planet occurrence statistics previously measured for early-to-mid M dwarfs, we predict that the survey will discover 1347$^{+208}_{-124}$ small transiting planets around these sources, each to a significance of 7.1$\sigma$ or greater. Significant departures from this prediction would test whether the occurrence rates of small planets increase or decrease around mid-M dwarfs and UCDs compared to early-M dwarfs. We predict the detection of 13$^{+4}_{-3}$ habitable, terrestrial planets ($R_p, Comment: 19 pages, 10 figures, accepted to AJ

About 70-80% of stars in our solar and galactic neighborhood are M dwarfs. They span a range of low masses and temperatures relative to solar-type stars, facilitating molecule formation throughout their atmospheres. Standard stellar atmosphere models primarily designed for FGK stars face challenges when characterizing broadband molecular features in spectra of cool stars. Here, we introduce SPHINX--a new 1-D self-consistent radiative-convective thermochemical equilibrium chemistry model grid of atmospheres and spectra for M dwarfs in low-resolution (R~250). We incorporate the latest pre-computed absorption cross-sections with pressure-broadening for key molecules dominant in late-K, early/main-sequence-M stars. We then validate our grid models by acquiring fundamental properties (Teff, log(g), [M/H], radius, and C/O) for 10 benchmark M+G binary stars with known host metallicities and 10 M dwarfs with interferometrically measured angular diameters. Incorporating a Gaussian-process inference tool Starfish, we account for correlated and systematic noise in low-resolution (spectral stitching of SpeX, SNIFS, and STIS) observations and derive robust estimates of fundamental M dwarf atmospheric parameters. Additionally, we assess the influence of photospheric heterogeneity on acquired [M/H] and find that it could explain some deviations from observations. We also probe whether the model-assumed convective mixing-length parameter influences inferred radii, effective temperature, and [M/H] and again find that may explain discrepancies between interferometry observations and model-derived stellar parameters for cooler M dwarfs. Mainly, we show the unique strength in leveraging broadband molecular absorption features occurring in low-resolution M dwarf spectra and demonstrate the ability to improve constraints on fundamental properties of exoplanet hosts and late brown dwarf companions., Comment: Submitted to ApJ

We describe the Perkins INfrared Exosatellite Survey (PINES), a near-infrared photometric search for short-period transiting planets and moons around a sample of 393 spectroscopically confirmed L- and T-type dwarfs. PINES is performed with Boston University's 1.8 m Perkins Telescope Observatory, located on Anderson Mesa, Arizona. We discuss the observational strategy of the survey, which was designed to optimize the number of expected transit detections, and describe custom automated observing procedures for performing PINES observations. We detail the steps of the $\texttt{PINES Analysis Toolkit}$ ($\texttt{PAT}$), software that is used to create light curves from PINES images. We assess the impact of second-order extinction due to changing precipitable water vapor on our observations and find that the magnitude of this effect is minimized in Mauna Kea Observatories $\textit{J}$-band. We demonstrate the validity of $\texttt{PAT}$ through the recovery of a transit of WASP-2 b and known variable brown dwarfs, and use it to identify a new variable L/T transition object: the T2 dwarf WISE J045746.08-020719.2. We report on the measured photometric precision of the survey and use it to estimate our transit detection sensitivity. We find that for our median brightness targets, assuming contributions from white noise only, we are sensitive to the detection of 2.5 $R_\oplus$ planets and larger. PINES will test whether the increase in sub-Neptune-sized planet occurrence with decreasing host mass continues into the L and T dwarf regime., Comment: 25 pages, 15 figures, accepted to AJ

We describe a new transit detection algorithm designed to detect single transit events in discontinuous Perkins INfrared Exosatellite Survey (PINES) observations of L and T dwarfs. We use this algorithm to search for transits in 131 PINES light curves and identify two transit candidates: 2MASS J18212815+1414010 (2MASS J1821+1414) and 2MASS J08350622+1953050 (2MASS J0835+1953). We disfavor 2MASS J1821+1414 as a genuine transit candidate due to the known variability properties of the source. We cannot rule out the planetary nature of 2MASS J0835+1953's candidate event and perform follow-up observations in an attempt to recover a second transit. A repeat event has yet to be observed, but these observations suggest that target variability is an unlikely cause of the candidate transit. We perform a Markov chain Monte Carlo simulation of the light curve and estimate a planet radius ranging from $4.2^{+3.5}_{-1.6}R_\oplus$ to $5.8^{+4.8}_{-2.1}R_\oplus$, depending on the host's age. Finally, we perform an injection and recovery simulation on our light curve sample. We inject planets into our data using measured M dwarf planet occurrence rates and attempt to recover them using our transit search algorithm. Our detection rates suggest that, assuming M dwarf planet occurrence rates, we should have roughly a 1$\%$ chance of detecting a candidate that could cause the transit depth we observe for 2MASS J0835+1953. If 2MASS J0835+1953 b is confirmed, it would suggest an enhancement in the occurrence of short-period planets around L and T dwarfs in comparison to M dwarfs, which would challenge predictions from planet formation models., Comment: 23 pages, 15 figures, accepted to AJ

V471 Tau is a post-common-envelope binary consisting of an eclipsing DA white dwarf and a K-type main-sequence star in the Hyades star cluster. We analyzed publicly available photometry and spectroscopy of V471 Tau to revise the stellar and orbital parameters of the system. We used archival K2 photometry, archival Hubble Space Telescope spectroscopy, and published radial-velocity measurements of the K-type star. Employing Gaussian processes to fit for rotational modulation of the system flux by the main-sequence star, we recovered the transits of the white dwarf in front of the main-sequence star for the first time. The transits are shallower than would be expected from purely geometric occultations owing to gravitational microlensing during transit, which places an additional constraint on the white-dwarf mass. Our revised mass and radius for the main-sequence star is consistent with single-star evolutionary models given the age and metallicity of the Hyades. However, as noted previously in the literature, the white dwarf is too massive and too hot to be the result of single-star evolution given the age of the Hyades, and may be the product of a merger scenario. We independently estimate the conditions of the system at the time of common envelope that would result in the measured orbital parameters today., Comment: 11 pages, 7 figures, Accepted to The Astronomical Journal

Full list of authors: Kemmer, J.; Stock, S.; Kossakowski, D.; Kaminski, A.; Molaverdikhani, K.; Schlecker, M.; Caballero, J. A.; Amado, P. J.; Astudillo-Defru, N.; Bonfils, X.; Ciardi, D.; Collins, K. A.; Espinoza, N.; Fukui, A.; Hirano, T.; Jenkins, J. M.; Latham, D. W.; Matthews, E. C.; Narita, N.; Pallé, E.; Parviainen, H.; Quirrenbach, A.; Reiners, A.; Ribas, I.; Ricker, G.; Schlieder, J. E.; Seager, S.; Vanderspek, R.; Winn, J. N.; Almenara, J. M.; Béjar, V. J. S.; Bluhm, P.; Bouchy, F.; Boyd, P.; Christiansen, J. L.; Cifuentes, C.; Cloutier, R.; Collins, K. I.; Cortés-Contreras, M.; Crossfield, I. J. M.; Crouzet, N.; de Leon, J. P.; Della-Rose, D. D.; Delfosse, X.; Dreizler, S.; Esparza-Borges, E.; Essack, Z.; Forveille, Th.; Figueira, P.; Galadí-Enríquez, D.; Gan, T.; Glidden, A.; Gonzales, E. J.; Guerra, P.; Harakawa, H.; Hatzes, A. P.; Henning, Th.; Herrero, E.; Hodapp, K.; Hori, Y.; Howell, S. B.; Ikoma, M.; Isogai, K.; Jeffers, S. V.; Kürster, M.; Kawauchi, K.; Kimura, T.; Klagyivik, P.; Kotani, T.; Kurokawa, T.; Kusakabe, N.; Kuzuhara, M.; Lafarga, M.; Livingston, J. H.; Luque, R.; Matson, R.; Morales, J. C.; Mori, M.; Muirhead, P. S.; Murgas, F.; Nishikawa, J.; Nishiumi, T.; Omiya, M.; Reffert, S.; Rodríguez López, C.; Santos, N. C.; Schöfer, P.; Schwarz, R. P.; Shiao, B.; Tamura, M.; Terada, Y.; Twicken, J. D.; Ueda, A.; Vievard, S.; Watanabe, N.; Zechmeister, M., We present the confirmation and characterisation of GJ 3473 b (G 50-16, TOI-488.01), a hot Earth-sized planet orbiting an M4 dwarf star, whose transiting signal (P = 1.1980035 ± 0.0000018 d) was first detected by the Transiting Exoplanet Survey Satellite (TESS). Through a joint modelling of follow-up radial velocity observations with CARMENES, IRD, and HARPS together with extensive ground-based photometric follow-up observations with LCOGT, MuSCAT, and MuSCAT2, we determined a precise planetary mass, Mb = 1.86 ± 0.30 M·, and radius, Rb = 1.264 ± 0.050 R·. Additionally, we report the discovery of a second, temperate, non-transiting planet in the system, GJ 3473 c, which has a minimum mass, Mc sin i = 7.41 ± 0.91 M·, and orbital period, Pc = 15.509 ± 0.033 d. The inner planet of the system, GJ 3473 b, is one of the hottest transiting Earth-sized planets known thus far, accompanied by a dynamical mass measurement, which makes it a particularly attractive target for thermal emission spectroscopy. © 2020 ESO., CARMENES is an instrument at the Centro Astronomico Hispano-Aleman de Calar Alto (CAHA, Almeria, Spain). CARMENES was funded by the German Max-Planck-Gesellschaft (MPG), the Spanish Consejo Superior de Investigaciones Cientificas (CSIC), the European Union through FEDER/ERF FICTS-2011-02 funds, and the members of the CARMENES Consortium (Max-Planck-Institut fur Astronomie, Instituto de Astrofisica de Andalucia, Landessternwarte Konigstuhl, Institut de Ciencies de l'Espai, Insitut fur Astrophysik Gottingen, Universidad Complutense de Madrid, Thuringer Landessternwarte Tautenburg, Instituto de Astrofisica de Canarias, Hamburger Sternwarte, Centro de Astrobiologia and Centro Astronomico Hispano-Aleman), with additional contributions by the Spanish Ministry of Economy, the German Science Foundation through the Major Research Instrumentation Programme and Deutsche Forschungsgemeinschaft (DFG) Research Unit FOR2544 "Blue Planets around Red Stars", the Klaus Tschira Stiftung, the states of Baden-Wurttemberg and Niedersachsen, and by the Junta de Andalucia. Part of this work was funded by the Ministerio de Ciencia e Innovacion via projects PID2019-109522GB-C51/2/3/4 and AYA2016-79425C3-1/2/3-P, (MEXT/)JSPS KAKENHI via grants 15H02063, JP17H04574, 18H05442, JP18H01265, JP18H05439, 19J11805, JP19K14783, and 22000005, JST PRESTO via grant JPMJPR1775, FCT/MCTES through national funds (PIDDAC, PTDC/FIS-AST/32113/2017) via grant UID/FIS/04434/2019, FEDER - Fundo Europeu de Desenvolvimento Regional through COMPETE2020 - Programa Operacional Competitividade e Internacionalizacao (POCI-010145-FEDER-032113), the FONDECYT project 3180063, Fundacao para a Ciencia e a Tecnologia through national funds and by FEDER through COMPETE2020 - Programa Operacional Competitividade e Internacionalizacao by these grants: UID/FIS/04434/2019; UIDB/04434/2020; UIDP/04434/2020; PTDC/FIS-AST/32113/2017 & POCI-01-0145-FEDER-032113; PTDC/FISAST/28953/2017 & POCI-01-0145-FEDER-028953. and the International Graduate Program for Excellence in Earth-Space Science. Part of the data analysis was carried out on the Multi-wavelength Data Analysis System operated by the Astronomy Data Center, National Astronomical Observatory of Japan. Funding for the TESS mission is provided by NASA's Science Mission directorate. We acknowledge the use of public TESS Alert data from pipelines at the TESS Science Office and at the TESS Science Processing Operations Center. This research has made use of the Exoplanet Follow-up Observation Program website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. Resources supporting this work were provided by the NASA High-End Computing Program through the NASA Advanced Supercomputing Division at Ames Research Center for the production of the SPOC data products. This paper includes data collected by the TESS mission, which are publicly available from the Mikulski Archive for Space Telescopes. This work makes use of observations from the LCOGT network. The analysis of this work has made use of a wide variety of public available software packages that are not referenced in the manuscript: Exo-Striker (Trifonov 2019), astropy (Astropy Collaboration 2018), scipy (Virtanen et al. 2020), numpy (Oliphant 2006), matplotlib (Hunter 2007), tqdm (da Costa-Luis 2019), pandas (The pandas development team 2020), seaborn (Waskom et al. 2020), lightkurve (Lightkurve Collaboration 2018) and PyFITS (Barrett et al. 2012).

We present $Spitzer$ 4.5$��$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s)., 14 Pages, 5 Figures, 2 Tables, Accepted to AJ

We develop an extension to an open source Python framework called “lightkurve” for high precision time series photometric signal injection and recovery with support for a range of astrophysical transients, including planetary transits, starspot rotation, stellar flares, and supernovae. These tools will provide support for measuring transit detection efficiency in ultracool dwarfs observed by K2. We are interested in our ability to recover various types of transits around dwarfs at the spectral type M/L transition; to investigate this, we inject transits that simulate planets from 0.5 to 3.5 Earth radii and of periods from 0.3 to 26 days and attempt to recover them. We find a threshold of detectability as a relation to orbital period and radius. These detection efficiencies allow us to calculate the probability of seeing planets around these stars as it relates to the planet occurrence rate for a certain type of planet. We can constrain the probability of seeing no planets around a given number of stars and put an upper limit on the planet occurrence rate. We have calculated an upper limit for planet occurrence for each type of planet around M/L dwarfs observed by K2.

We characterized the effects of telluric absorption lines on the radial velocity (RV) precision of stellar spectra taken through an iodine cell. To isolate the effects induced by telluric contamination from other stellar, instrumental, or numerical systematic RV noise, we extracted RVs from simulated iodine calibrated spectra of three RV standard stars regularly observed by Keck/HIRES. We add in water absorption lines according to measured precipitable water vapor (PWV) contents over a one-year period. We concluded that telluric contamination introduces additional RV noise and spurious periodic signals on the level of 10-20 cm/s, consistent with similar previous studies. Our findings show that forward modeling the telluric lines effectively recovers the RV precision and accuracy, with no prior knowledge of the PWV needed. Such a recovery is less effective when the water absorption lines are relatively deep in the stellar template used in the forward modeling. Overall, telluric contamination plays an insignificant role for typical iodine-calibrated RV programs aiming at ~1-2 m/s, but we recommend adding modeling of telluric lines and taking stellar template observations on nights with low humidity for programs aiming to achieve sub-m/s precision., Comment: 17 pages, 9 figures. AJ accepted

Planets occur most frequently around cool dwarfs, but only a handful of specific examples are known to orbit the latest-type M stars. Using TESS photometry, we report the discovery of two planets transiting the low-mass star called LP 791-18 (identified by TESS as TOI 736). This star has spectral type M6V, effective temperature 2960 K, and radius 0.17 R_Sun, making it the third-coolest star known to host planets. The two planets straddle the radius gap seen for smaller exoplanets; they include a 1.1 R_Earth planet on a 0.95 day orbit and a 2.3 R_Earth planet on a 5 day orbit. Because the host star is small the loss of light during these planets' transits is fairly large (0.4% and 1.7%). This has allowed us to detect both planets' transits from ground-based photometry, refining their radii and orbital ephemerides. In the future, radial velocity observations and transmission spectroscopy can both probe these planets' bulk interior and atmospheric compositions, and additional photometric monitoring would be sensitive to even smaller transiting planets., 20 pages, 7 figures, 2 tables, 2 planets, 1 cool star. Submitted to AAS Journals

We present design considerations for a ground-based survey for transiting exoplanets around L and T dwarfs, spectral classes that have yet to be thoroughly probed for planets. We simulate photometry for L and T targets with a variety of red-optical and near-infrared detectors, and compare the scatter in the photometry to anticipated transit depths. Based on these results, we recommend the use of a low-dark-current detector with H-band NIR photometric capabilities. We then investigate the potential for performing a survey for Earth-sized planets for a variety of telescope sizes. We simulate planetary systems around a set of spectroscopically confirmed L and T dwarfs using measured M dwarf planet occurrence rates from $\textit{Kepler}$, and simulate their observation in surveys ranging in duration from 120 to 600 nights, randomly discarding 30% of nights to simulate weather losses. We find that an efficient survey design uses a 2-meter-class telescope with a NIR instrument and 360-480 observing nights, observing multiple L and T targets each night with a dithering strategy. Surveys conducted in such a manner have over an 80% chance of detecting at least one planet, and detect around 2 planets, on average. The number of expected detections depends on the true planet occurrence rate, however, which may in fact be higher for L and T dwarfs than for M dwarfs., 18 pages, 9 figures; accepted for publication in PASP

Low-mass eclipsing binaries show systematically larger radii than model predictions for their mass, metallicity and age. Prominent explanations for the inflation involve enhanced magnetic fields generated by rapid rotation of the star that inhibit convection and/or suppress flux from the star via starspots. However, derived masses and radii for individual eclipsing binary systems often disagree in the literature. In this paper, we continue to investigate low-mass eclipsing binaries (EBs) observed by NASA's {\it Kepler} spacecraft, deriving stellar masses and radii using high-quality space-based light curves and radial velocities from high-resolution infrared spectroscopy. We report masses and radii for three {\it Kepler} EBs, two of which agree with previously published masses and radii (KIC 11922782 and KIC 9821078). For the third EB (KIC 7605600), we report new masses and show the secondary component is likely fully convective ($M_2 = 0.17 \pm 0.01 M_{\sun}$ and $R_2 = 0.199^{+0.001}_{-0.002} R_{\sun}$). Combined with KIC 10935310 from Han et al. (2017), we find that the masses and radii for four low-mass {\it Kepler} EBs are consistent with modern stellar evolutionary models for M dwarf stars and do not require inhibited convection by magnetic fields to account for the stellar radii., Comment: arXiv admin note: text overlap with arXiv:1707.07001

Previous studies of planet occurrence rates largely relied on photometric stellar characterizations. In this paper, we present planet occurrence rates for mid-type M dwarfs using spectroscopy, parallaxes, and photometry to determine stellar characteristics. Our spectroscopic observations have allowed us to constrain spectral type, temperatures, and in some cases metallicities for 337 out of 561 probable mid-type M dwarfs in the primary Kepler field. We use a random forest classifier to assign a spectral type to the remaining 224 stars. Combining our data with Gaia parallaxes, we compute precise ($\sim$3%) stellar radii and masses, which we use to update planet parameters and planet occurrence rates for Kepler mid-type M dwarfs. Within the Kepler field, there are seven M3 V to M5 V stars which host 13 confirmed planets between 0.5 and 2.5 Earth radii and at orbital periods between 0.5 and 10 days. For this population, we compute a planet occurrence rate of $1.19^{+0.70}_{-0.49}$ planets per star. For M3 V, M4 V, and M5 V, we compute planet occurrence rates of $0.86^{+1.32}_{-0.68}$, $1.36^{+2.30}_{-1.02}$, and $3.07^{+5.49}_{-2.49}$ planets per star, respectively., Comment: 35 pages, 10 figures, 4 tables, Accepted for publication in AJ

High-resolution, near-infrared spectra will be the primary tool for finding and characterizing Earth-like planets around low-mass stars. Yet, the properties of exoplanets can not be precisely determined without accurate and precise measurements of the host star. Spectra obtained with the Immersion GRating INfrared Spectrometer (IGRINS) simultaneously provide diagnostics for most stellar parameters, but the first step in any analysis is the determination of the effective temperature. Here we report the calibration of high-resolution H-band spectra to accurately determine effective temperature for stars between 4000-3000 K ($\sim$K8--M5) using absorption line depths of Fe I, OH, and Al I. The field star sample used here contains 254 K and M stars with temperatures derived using BT-Settl synthetic spectra. We use 106 stars with precise temperatures in the literature to calibrate our method with typical errors of about 140 K, and systematic uncertainties less than $\sim$120 K. For the broadest applicability, we present T$_{\rm eff}$--line-depth-ratio relationships, which we test on 12 members of the TW Hydrae Association and at spectral resolving powers between $\sim$10,000--120,000. These ratios offer a simple but accurate measure of effective temperature in cool stars that is distance and reddening independent., Comment: 19 pages, 11 figures and 3 tables. Accepted in ApJ

We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system., 27 pages, 22 figures, AJ accepted

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

We perform a new analysis of the M dwarf-M dwarf eclipsing binary system NSVS 07394765 in order to investigate the reported hyper-inflated radius of one of the component stars. Our analysis is based on archival photometry from the Wide Angle Search for Planets (WASP), new photometry from the 32 cm {Command Module Observatory (CMO) telescope in Arizona and the 70 cm telescope at Thacher Observatory in California}, and new high-resolution infrared spectra obtained with the Immersion Grating Infrared Spectrograph (IGRINS) on the Discovery Channel Telescope. The masses and radii we measure for each component star disagree with previously reported measurements. We show that both stars are early M-type main-sequence stars without evidence for youth or hyper-inflation ($M_1= 0.661\ ^{+0.008}_{-0.036}\ \rm{M_{sun}}$, $M_2= 0.608\ ^{+0.003}_{-0.028}\ \rm{M_{sun}}$, $R_1= 0.599\ ^{+0.032}_{-0.019}\ \rm{R_{sun}}$, $R_2= 0.625\ ^{+0.012}_{-0.027}\ \rm{R_{sun}}$), and we update the orbital period and eclipse ephemerides for the system. We suggest that the likely cause of the initial hyper-inflated result is the use of moderate-resolution spectroscopy for precise radial velocity measurements., Comment: 14 pages, 11 figures. Accepted for publication in The Astronomical Journal

The occurrence of planets orbiting ultracool dwarfs is poorly constrained. We present results from a Guest Observer program on NASA's K2 spacecraft to search for transiting planets orbiting a sample of 827 ultracool dwarfs. Having found no transiting planets in our sample, we determined an upper limit on the occurrence of planets. We simulated planets orbiting our sample for a range of orbital periods and sizes. For the simulated planets that transit their host, we injected the transit light curve into the real K2 light curves, then attempted to recover the injected planets. For a given occurrence rate, we calculated the probability of seeing no planets, and use the results to place an upper limit on planet occurrence as a function of planet radius and orbital period. We find that short period, mini-Neptune- and Jupiter-sized planets are rare around ultracool dwarfs, consistent with results for early- and mid-type M dwarf stars. We constrain the occurrence rate $\eta$ for planets between $0.5$ and $10$ R$_{\oplus}$ with orbital periods between 1 and 26.3 days., Comment: 9 pages, 7 figures, 1 table. Machine-readable table in source. Accepted to the Astronomical Journal

We report the discovery of GJ 1252 b, a planet with a radius of 1.193 $\pm$ 0.074 $R_{\oplus}$ and an orbital period of 0.52 days around an M3-type star (0.381 $\pm$ 0.019 $M_{\odot}$, 0.391 $\pm$ 0.020 $R_{\odot}$) located 20.385 $\pm$ 0.019 pc away. We use TESS data, ground-based photometry and spectroscopy, Gaia astrometry, and high angular resolution imaging to show that the transit signal seen in the TESS data must originate from a transiting planet. We do so by ruling out all false positive scenarios that attempt to explain the transit signal as originating from an eclipsing stellar binary. Precise Doppler monitoring also leads to a tentative mass measurement of 2.09 $\pm$ 0.56 $M_{\oplus}$. The host star proximity, brightness ($V$ = 12.19 mag, $K$ = 7.92 mag), low stellar activity, and the system's short orbital period make this planet an attractive target for detailed characterization, including precise mass measurement, looking for other objects in the system, and planet atmosphere characterization., Accepted to ApJ Letters

In young sun-like stars and field M dwarf stars, chromospheric and coronal magnetic activity indicators such as H$\alpha$, X-ray and radio emission are known to saturate with low Rossby number ($Ro \lesssim 0.1$), defined as the ratio of rotation period to convective turnover time. The mechanism for the saturation is unclear. In this paper, we use photospheric Ti I and Ca I absorption lines in $Y$ band to investigate magnetic field strength in M dwarfs for Rossby numbers between 0.01 and 1.0. The equivalent widths of the lines are magnetically enhanced by photospheric spots, a global field or a combination of the two. The equivalent widths behave qualitatively similar to the chromospheric and coronal indicators: we see increasing equivalent widths (increasing absorption) with decreasing $Ro$ and saturation of the equivalent widths for $Ro \lesssim 0.1$. The majority of M dwarfs in this study are fully convective. The results add to mounting evidence that the magnetic saturation mechanism occurs at or beneath the stellar photosphere., Comment: Accepted to AJ. Machine readable tables in the source files. Reduced spectra will be available in the published version

The physics of molecular vibration causes absorption spectra of atmospheric molecules to be a group of approximately periodic fine lines. This is fortuitous for detecting exoplanet biosignificant molecules, since it approximately matches the periodic sinusoidal transmission of an interferometer. The series addition of a 0.6 cm interferometer with a dispersive spectrograph creates moire patterns. These enhance detection by several orders of magnitude for initially low resolution spectrographs. We simulate the Gemini Planet Imager integral field spectrograph observing a telluric spectrum of native resolutions 40 and 70 for 1.65 and 2 micron bands– too low to resolve the fine lines. The interferometer addition increases the detectability of the molecular signal, relative to photon noise, to a level similar to a R=4400 (at 1.65 micron) or R=3900 (at 2 micron) spectrograph.

We present a catalog of cool dwarf targets ($V-J>2.7$, $T_{\rm eff} \lesssim 4000 K$) and their stellar properties for the upcoming Transiting Exoplanet Survey Satellite (TESS), for the purpose of determining which cool dwarfs should be observed using two-minute observations. TESS has the opportunity to search tens of thousands of nearby, cool, late K and M-type dwarfs for transiting exoplanets, an order of magnitude more than current or previous transiting exoplanet surveys, such as {\it Kepler}, K2 and ground-based programs. This necessitates a new approach to choosing cool dwarf targets. Cool dwarfs were chosen by collating parallax and proper motion catalogs from the literature and subjecting them to a variety of selection criteria. We calculate stellar parameters and TESS magnitudes using the best possible relations from the literature while maintaining uniformity of methods for the sake of reproducibility. We estimate the expected planet yield from TESS observations using statistical results from the Kepler Mission, and use these results to choose the best targets for two-minute observations, optimizing for small planets for which masses can conceivably be measured using follow up Doppler spectroscopy by current and future Doppler spectrometers. The catalog is incorporated into the TESS Input Catalog and TESS Candidate Target List until a more complete and accurate cool dwarf catalog identified by ESA's Gaia Mission can be incorporated., Accepted to The Astronomical Journal. For the full catalog, please contact the corresponding author

We present the spectroscopic orbit of LHS 1610A, a newly discovered single-lined spectroscopic binary with a trigonometric distance placing it at 9.9 pm 0.2 pc. We obtained spectra with the TRES instrument on the 1.5m Tillinghast Reflector at the Fred Lawrence Whipple Observatory located on Mt. Hopkins in AZ. We demonstrate the use of the TiO molecular bands at 7065 -- 7165 Angstroms to measure radial velocities and achieve an average estimated velocity uncertainty of 28 m/s. We measure the orbital period to be 10.6 days and calculate a minimum mass of 44.8 pm 3.2 Jupiter masses for the secondary, indicating that it is likely a brown dwarf. We place an upper limit to 3 sigma of 2500 K on the effective temperature of the companion from infrared spectroscopic observations using IGRINS on the 4.3m Discovery Channel Telescope. In addition, we present a new photometric rotation period of 84.3 days for the primary star using data from the MEarth-South Observatory, with which we show that the system does not eclipse., 10 pages, 5 figures; accepted for publication in the Astronomical Journal

M subdwarfs are low-metallicity M dwarfs that typically inhabit the halo population of the Galaxy. Metallicity controls the opacity of stellar atmospheres; in metal poor stars, hydrostatic equilibrium is reached at a smaller radius, leading to smaller radii for a given effective temperature. We compile a sample of 88 stars that span spectral classes K7 to M6 and include stars with metallicity classes from solar-metallicity dwarf stars to the lowest metallicity ultra-subdwarfs to test how metallicity changes the stellar radius. We fit models to Palomar Double Spectrograph (DBSP) optical spectra to derive effective temperatures ($T_\mathrm{eff}$) and we measure bolometric luminosities ($L_\mathrm{bol}$) by combining broad wavelength-coverage photometry with Gaia parallaxes. Radii are then computed by combining the $T_\mathrm{eff}$ and $L_\mathrm{bol}$ using the Stefan-Boltzman law. We find that for a given temperature, ultra-subdwarfs can be as much as five times smaller than their solar-metallicity counterparts. We present color-radius and color-surface brightness relations that extend down to [Fe/H] of $-$2.0 dex, in order to aid the radius determination of M subdwarfs, which will be especially important for the WFIRST exoplanetary microlensing survey., Comment: Accepted to AJ

M dwarf stars are exciting targets for exoplanet investigations; however, their fundamental stellar properties are difficult to measure. Perhaps the most challenging property to measure is stellar age. Once on the main sequence, M dwarfs change imperceptibly in their temperature and luminosity, necessitating novel statistical techniques for estimating their ages. In this paper, we infer ages for known eccentric-planet-hosting M dwarfs using a combination of kinematics and $\alpha$-element-enrichment, both shown to correlate with age for Sun-like FGK stars. We calibrate our method on FGK stars in a Bayesian context. To measure $\alpha$-enrichment, we use publicly-available spectra from the CARMENES exoplanet survey and a recently developed [Ti/Fe] calibration utilizing individual Ti I and Fe I absorption lines in $Y$ band. Tidal effects are expected to circularize the orbits of short-period planets on short timescales; however, we find a number of mildly eccentric, close-in planets orbiting old ($\sim$8 Gyr) stars. For these systems, we use our ages to constrain the tidal dissipation parameter of the planets, $Q_\mathrm{p}$. For two mini-Neptune planets, GJ 176b and GJ 536b, we find they have $Q_\mathrm{p}$ values more similar to the ice giants than the terrestrial planets in our Solar System. For GJ 436b, we estimate an age of $8.9^{+2.3}_{-2.1}$ Gyr and constrain the $Q_\mathrm{p}$ to be $>10^5$, in good agreement with constraints from its inferred tidal heating. We find that GJ 876d has likely undergone significant orbital evolution over its $8.4^{+2.2}_{-2.0}$ Gyr lifetime, potentially influenced by its three outer companions which orbit in a Laplace resonance., Comment: accepted for publication in ApJ

We explore the radial velocity performance benefits of coupling starlight to a fast-scanning interferometer and a fast-readout spectrometer with zero readout noise. By rapidly scanning an interferometer, we can decouple wavelength calibration errors from precise radial velocity measurements, exploiting the advantages of lock-in amplification. In a Bayesian framework, we investigate the correlation between wavelength calibration errors and resulting radial velocity errors. We construct an end-to-end simulation of this approach to address the feasibility of achieving 10 cm s^(-1) radial velocity precision on a typical Sun-like star using existing, 5 m-class telescopes. We find that such a precision can be reached in a single night, opening up possibilities for ground-based detections of Earth-Sun analog systems.

The Kepler mission has shown that small planets are extremely common. It is likely that nearly every star in the sky hosts at least one rocky planet. We just need to look hard enough - but this requires vast amounts of telescope time. MINERVA (MINiature Exoplanet Radial Velocity Array) is a dedicated exoplanet observatory with the primary goal of discovering rocky, Earth-like planets orbiting in the habitable zone of bright, nearby stars. The MINERVA team is a collaboration among UNSW Australia, Harvard-Smithsonian Center for Astrophysics, Penn State University, University of Montana, and the California Institute of Technology. The four-telescope MINERVA array will be sited at the F.L. Whipple Observatory on Mt Hopkins in Arizona, USA. Full science operations will begin in mid-2015 with all four telescopes and a stabilised spectrograph capable of high-precision Doppler velocity measurements. We will observe 100 of the nearest, brightest, Sun-like stars every night for at least five years. Detailed simulations of the target list and survey strategy lead us to expect 154 new low-mass planets.

The ability to perform detailed chemical analysis of Sun-like F-, G-, and K-type stars is a powerful tool with many applications including studying the chemical evolution of the Galaxy and constraining planet formation theories. Unfortunately, complications in modeling cooler stellar atmospheres hinders similar analysis of M-dwarf stars. Empirically-calibrated methods to measure M dwarf metallicity from moderate-resolution spectra are currently limited to measuring overall metallicity and rely on astrophysical abundance correlations in stellar populations. We present a new, empirical calibration of synthetic M dwarf spectra that can be used to infer effective temperature, Fe abundance, and Ti abundance. We obtained high-resolution (R~25,000), Y-band (~1 micron) spectra of 29 M dwarfs with NIRSPEC on Keck II. Using the PHOENIX stellar atmosphere modeling code (version 15.5), we generated a grid of synthetic spectra covering a range of temperatures, metallicities, and alpha-enhancements. From our observed and synthetic spectra, we measured the equivalent widths of multiple Fe I and Ti I lines and a temperature-sensitive index based on the FeH bandhead. We used abundances measured from widely-separated solar-type companions to empirically calibrate transformations to the observed indices and equivalent widths that force agreement with the models. Our calibration achieves precisions in Teff, [Fe/H], and [Ti/Fe] of 60 K, 0.1 dex, and 0.05 dex, respectively and is calibrated for 3200 K < Teff < 4100 K, -0.7 < [Fe/H] < +0.3, and -0.05 < [Ti/Fe] < +0.3. This work is a step toward detailed chemical analysis of M dwarfs at a similar precision achieved for FGK stars., accepted for publication in ApJ, all synthetic spectra available at http://people.bu.edu/mveyette/phoenix/

Several low-mass eclipsing binary stars show larger than expected radii for their measured mass, metallicity and age. One proposed mechanism for this radius inflation involves inhibited internal convection and starspots caused by strong magnetic fields. One particular eclipsing binary, T-Cyg1-12664, has proven confounding to this scenario. \citet{Cakirli2013a} measured a radius for the secondary component that is twice as large as model predictions for stars with the same mass and age, but a primary mass that is consistent with predictions. \citet[][]{Iglesias2017} independently measured the radii and masses of the component stars and found that the radius of the secondary is not in fact inflated with respect to models, but that the primary is, consistent with the inhibited convection scenario. However, in their mass determinations, \citet[][]{Iglesias2017} lacked independent radial velocity measurements for the secondary component due to the star's faintness at optical wavelengths. The secondary component is especially interesting as its purported mass is near the transition from partially-convective to a fully-convective interior. In this article we independently determined the masses and radii of the component stars of T-Cyg1-12664 using archival {\it Kepler} data and radial velocity measurements of both component stars obtained with IGRINS on the Discovery Channel Telescope and NIRSPEC and HIRES on the Keck Telescopes. We show that neither of the component stars is inflated with respect to models. Our results are broadly consistent with modern stellar evolutionary models for main-sequence M dwarf stars and do not require inhibited convection by magnetic fields to account for the stellar radii.

We report parallaxes and proper motions from the Hawaii Infrared Parallax Program for eight nearby M dwarf stars with transiting exoplanets discovered by Kepler. We combine our directly measured distances with mass-luminosity and radius-luminosity relationships to significantly improve constraints on the host stars' properties. Our astrometry enables the identification of wide stellar companions to the planet hosts. Within our limited sample, all the multi-transiting planet hosts (three of three) appear to be single stars, while nearly all (four of five) of the systems with a single detected planet have wide stellar companions. By applying strict priors on average stellar density from our updated radius and mass in our transit fitting analysis, we measure the eccentricity probability distributions for each transiting planet. Planets in single-star systems tend to have smaller eccentricities than those in binaries, although this difference is not significant in our small sample. In the case of Kepler-42bcd, where the eccentricities are known to be $\simeq$0, we demonstrate that such systems can serve as powerful tests of M dwarf evolutionary models by working in $L_*-\rho_*$ space. The transit-fit density for Kepler-42bcd is inconsistent with model predictions at 2.1$\sigma$ (22%), but matches more empirical estimates at 0.2$\sigma$ (2%), consistent with earlier results showing model radii of M dwarfs are under-inflated. Gaia will provide high precision parallaxes for the entire Kepler M dwarf sample, and TESS will identify more planets transiting nearby, late-type stars, enabling significant improvements in our understanding of the eccentricity distribution of small planets and the parameters of late-type dwarfs., Comment: Accepted to AJ. Updated with corrections from proofs (May 22, 2017)

The Transiting Exoplanet Survey Satellite (TESS) will be conducting a nearly all-sky photometric survey over two years, with a core mission goal to discover small transiting exoplanets orbiting nearby bright stars. It will obtain 30-minute cadence observations of all objects in the TESS fields of view, along with 2-minute cadence observations of 200,000 to 400,000 selected stars. The choice of which stars to observe at the 2-min cadence is driven by the need to detect small transiting planets, which leads to the selection of primarily bright, cool dwarfs. We describe the catalogs assembled and the algorithms used to populate the TESS Input Catalog (TIC). We also describe a ranking system for prioritizing stars according to the smallest transiting planet detectable, and assemble a Candidate Target List (CTL) using that ranking. We discuss additional factors that affect the ability to photometrically detect and dynamically confirm small planets, and we note additional stellar populations of interest that may be added to the final target list. The TIC is available on the STScI MAST server, and an enhanced CTL is available through the Filtergraph data visualization portal system at the URL https://filtergraph.vanderbilt.edu/tess_ctl ., Comment: Published in The Astronomical Journal, 56 pages, 24 figures, 6 appendices. This is the definitive version of the documentation for the TIC and CTL as of TESS launch. This updated version includes a correction to Figure 5 and associated text as an erratum on the first page (erratum will also appear in AJ)

The Minerva-Australis telescope array is a facility dedicated to the follow-up, confirmation, characterisation, and mass measurement of bright transiting planets discovered by the Transiting Exoplanet Survey Satellite (TESS) -- a category in which it is almost unique in the southern hemisphere. It is located at the University of Southern Queensland's Mount Kent Observatory near Toowoomba, Australia. Its flexible design enables multiple 0.7m robotic telescopes to be used both in combination, and independently, for high-resolution spectroscopy and precision photometry of TESS transit planet candidates. Minerva-Australis also enables complementary studies of exoplanet spin-orbit alignments via Doppler observations of the Rossiter-McLaughlin effect, radial velocity searches for non-transiting planets, planet searches using transit timing variations, and ephemeris refinement for TESS planets. In this first paper, we describe the design, photometric instrumentation, software, and science goals of Minerva-Australis, and note key differences from its Northern hemisphere counterpart -- the Minerva array. We use recent transit observations of four planets--WASP-2b, WASP-44b, WASP-45b, and HD 189733b to demonstrate the photometric capabilities of Minerva-Australis., Comment: 33 pages, 7 figures, 5 tables. Published in PASP on November 2019

Aims.L 98-59 (TIC 307210830, TOI-175) is a nearby M3 dwarf around which TESS revealed three small transiting planets (0.80, 1.35, 1.57 Earth radii) in a compact configuration with orbital periods shorter than 7.5 days. Here we aim to measure the masses of the known transiting planets in this system using precise radial velocity (RV) measurements taken with the HARPS spectrograph.Methods.We considered both trained and untrained Gaussian process regression models of stellar activity, which are modeled simultaneously with the planetary signals. Our RV analysis was then supplemented with dynamical simulations to provide strong constraints on the planets’ orbital eccentricities by requiring long-term stability.Results.We measure the planet masses of the two outermost planets to be 2.42 ± 0.35 and 2.31 ± 0.46 Earth masses, which confirms the bulk terrestrial composition of the former and eludes to a significant radius fraction in an extended gaseous envelope for the latter. We are able to place an upper limit on the mass of the smallest, innermost planet of e< 0.1.Conclusions.L 98-59 is likely a compact system of two rocky planets plus a third outer planet with a lower bulk density possibly indicative of the planet having retained a modest atmosphere. The system offers a unique laboratory for studies of planet formation, dynamical stability, and comparative atmospheric planetology as the two outer planets are attractive targets for atmospheric characterization through transmission spectroscopy. Continued RV monitoring will help refine the characterization of the innermost planet and potentially reveal additional planets in the system at wider separations.

WD 1145+017 is a unique white dwarf system that has a heavily polluted atmosphere, an infrared excess from a dust disk, numerous broad absorption lines from circumstellar gas, and changing transit features, likely from fragments of an actively disintegrating asteroid. Here, we present results from a large photometric and spectroscopic campaign with Hubble, Keck , VLT, Spitzer, and many other smaller telescopes from 2015 to 2018. Somewhat surprisingly, but consistent with previous observations in the u' band, the UV transit depths are always shallower than those in the optical. We develop a model that can quantitatively explain the observed "bluing" and the main findings are: I. the transiting objects, circumstellar gas, and white dwarf are all aligned along our line of sight; II. the transiting object is blocking a larger fraction of the circumstellar gas than of the white dwarf itself. Because most circumstellar lines are concentrated in the UV, the UV flux appears to be less blocked compared to the optical during a transit, leading to a shallower UV transit. This scenario is further supported by the strong anti-correlation between optical transit depth and circumstellar line strength. We have yet to detect any wavelength-dependent transits caused by the transiting material around WD 1145+017., Comment: 16 pages, 11 figures, 6 tables, ApJ, in press

We describe demonstrations of remarkable robustness to instrumental noises by using a multiple delay externally dispersed interferometer (EDI) on stellar observations at the Hale telescope. Previous observatory EDI demonstrations used a single delay. The EDI (also called “TEDI”) boosted the 2,700 resolution of the native TripleSpec NIR spectrograph (950-2450 nm) by as much as 10x to 27,000, using 7 overlapping delays up to 3 cm. We observed superb rejection of fixed pattern noises due to bad pixels, since the fringing signal responds only to changes in multiple exposures synchronous to the applied delay dithering. Remarkably, we observed a ~20x reduction of reaction in the output spectrum to PSF shifts of the native spectrograph along the dispersion direction, using our standard processing. This allowed high resolution observations under conditions of severe and irregular PSF drift otherwise not possible without the interferometer. Furthermore, we recently discovered an improved method of weighting and mixing data between pairs of delays that can theoretically further reduce the net reaction to PSF drift to zero. We demonstrate a 350x reduction in reaction to a native PSF shift using a simple simulation. This technique could similarly reduce radial velocity noise for future EDI’s that use two delays overlapped in delay space (or a single delay overlapping the native peak). Finally, we show an extremely high dynamic range EDI measurement of our ThAr lamp compared to a literature ThAr spectrum, observing weak features (~0.001x height of nearest strong line) that occur between the major lines. Because of individuality of each reference lamp, accurate knowledge of its spectrum between the (unfortunately) sparse major lines is important for precision radial velocimetry.

We present an updated optical and mechanical design of NEWS: the Near-infrared Echelle for Wide-band Spectroscopy (formerly called HiJaK: the High-resolution J, H and K spectrometer), a compact, high-resolution, near-infrared spectrometer for 5-meter class telescopes. NEWS provides a spectral resolution of 60,000 and covers the full 0.8-2.5 micron range in 5 modes. We adopt a compact, lightweight, monolithic design and developed NEWS to be mounted to the instrument cube at the Cassegrain focus of the the new 4.3-meter Discovery Channel Telescope., Comment: Proc. SPIE 9908, Ground-based and Airborne Instrumentation for Astronomy VI, 99086M (August 9, 2016)