Source author record

Arpita Roy

Arpita Roy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

68works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

68 published item(s)

preprint2026arXiv

The Lazuli Space Observatory: Architecture & Capabilities

The Lazuli Space Observatory is a 3-meter aperture astronomical facility designed for rapid-response observations and precision astrophysics across visible to near-infrared wavelengths (400-1700 nm bandpass). An off-axis, freeform telescope delivers diffraction-limited image quality (Strehl $>$0.8 at 633 nm) to three instruments across a wide, flat focal plane. The three instruments provide complementary capabilities: a Wide-field Context Camera (WCC) delivers multi-band imaging over a 35' $\times$ 12' footprint with high-cadence photometry; an Integral Field Spectrograph (IFS) provides continuous 400-1700 nm spectroscopy at R $\sim$ 100-500 for stable spectrophotometry; and an ExtraSolar Coronagraph (ESC) enables high-contrast imaging expected to reach raw contrasts of $10^{-8}$ and post-processed contrasts approaching $10^{-9}$. Operating from a 3:1 lunar-resonant orbit, Lazuli will respond to targets of opportunity in under four hours--a programmatic requirement designed to enable routine temporal responsiveness that is unprecedented for a space telescope of this size. Lazuli's technical capabilities are shaped around three broad science areas: (1) time-domain and multi-messenger astronomy, (2) stars and planets, and (3) cosmology. These capabilities enable a potent mix of science spanning gravitational wave counterpart characterization, fast-evolving transients, Type Ia supernova cosmology, high-contrast exoplanet imaging, and spectroscopy of exoplanet atmospheres. While these areas guide the observatory design, Lazuli is conceived as a general-purpose facility capable of supporting a wide range of astrophysical investigations, with open time for the global community. We describe the observatory architecture and capabilities in the preliminary design phase, with science operations anticipated following a rapid development cycle from concept to launch.

preprint2026arXiv

TOI-4495: A Pair of Aligned, Near-Resonant Sub-Neptunes that Likely Experienced Overstable Migration

We report the discovery of a sub-Neptune and a Neptune-like planet ($R_b = 2.48^{+0.14}_{-0.10}\,R_\oplus$, $R_c = 4.03^{+0.23}_{-0.15}\,R_\oplus$) orbiting the F-type star TOI-4495. The planets have orbital periods of 2.567 days and 5.185 days, lying close to a 2:1 mean-motion resonance (MMR). Our photodynamical analysis of the TESS light curves constrains the planetary masses to $M_b = 7.7 \pm 1.4\,M_\oplus$ and $M_c = 23.2 \pm 4.7\,M_\oplus$. The measured masses and radii indicate the presence of volatile-rich gaseous envelopes on both planets. The Rossiter-McLaughlin effect and the Doppler shadow of TOI-4495 c reveal a well-aligned orbit with a projected stellar obliquity of $λ= -2.3^{+8.3}_{-7.8}\,\mathrm{deg}$. Combined with the low mutual inclination constrained by the photodynamical analysis ($ΔI < 8.7\,\mathrm{deg}$), the planetary orbits are likely coplanar and aligned with the host star's spin axis. We show that the planets are near, but not in, the 2:1 MMR, with a circulating resonant angle. We also find substantial free eccentricity for the inner planet, TOI-4495 b ($e_b = 0.078^{+0.020}_{-0.013}$). Given the observed proximity to the 2:1 resonance and the more massive outer planet, TOI-4495 b and c are particularly susceptible to resonant overstability, which can convert resonantly excited eccentricity into free eccentricity. However, additional mechanisms (e.g., planetesimal scattering) may be required to further excite the eccentricity by $\sim 4\%$. To prevent tidal damping from reducing the eccentricity below the observed level over the star's lifetime (1.9 Gyr), the reduced tidal quality factor of TOI-4495 b must be $Q' \gtrsim 10^5$, consistent with the presence of a thick envelope on the planet.

preprint2024arXiv

Giant Outer Transiting Exoplanet Mass (GOT 'EM) Survey. IV. Long-term Doppler Spectroscopy for 11 Stars Thought to Host Cool Giant Exoplanets

Discovering and characterizing exoplanets at the outer edge of the transit method's sensitivity has proven challenging owing to geometric biases and the practical difficulties associated with acquiring long observational baselines. Nonetheless, a sample of giant exoplanets on orbits longer than 100 days has been identified by transit hunting missions. We present long-term Doppler spectroscopy for 11 such systems with observation baselines spanning a few years to a decade. We model these radial velocity observations jointly with transit photometry to provide initial characterizations of these objects and the systems in which they exist. Specifically, we make new precise mass measurements for four long-period giant exoplanets (Kepler-111 c, Kepler-553 c, Kepler-849 b, and PH-2 b), we place new upper limits on mass for four others (Kepler-421 b, KOI-1431.01, Kepler-1513 b, and Kepler-952 b), and we show that several "confirmed" planets are in fact not planetary at all. We present these findings to complement similar efforts focused on closer-in short-period giant planets, and with the hope of inspiring future dedicated studies of cool giant exoplanets.

preprint2022arXiv

A close-in puffy Neptune with hidden friends: The enigma of TOI 620

We present the validation of a transiting low-density exoplanet orbiting the M2.5 dwarf TOI 620 discovered by the NASA TESS mission. We utilize photometric data from both TESS and ground-based follow-up observations to validate the ephemerides of the 5.09-day transiting signal and vet false positive scenarios. High-contrast imaging data are used to resolve the stellar host and exclude stellar companions at separations $\gtrsim 0.2''$. We obtain follow-up spectroscopy and corresponding precise radial velocities (RVs) with multiple PRV spectrographs to confirm the planetary nature of the transiting exoplanet. We calculate a 5$σ$ upper limit of $M_P < 7.1$ M$_\oplus$ and $ρ_P < 0.74$ g cm$^{-3}$, and we identify a non-transiting 17.7-day candidate. We also find evidence for a substellar (1-20 M$_{\rm J}$) companion with a projected separation $\lesssim 20$ au from a combined analysis of Gaia, AO imaging, and RVs. With the discovery of this outer companion, we carry out a detailed exploration of the possibilities that TOI 620 b might instead be a circum-secondary planet or a pair of eclipsing binary stars orbiting the host in a hierarchical triple system. We find, under scrutiny, that we can exclude both of these scenarios from the multi-wavelength transit photometry, thus validating TOI 620 b as a low-density exoplanet transiting the central star in this system. The low density of TOI 620 b makes it one of the most amenable exoplanets for atmospheric characterization, such as with JWST and Ariel, validated or confirmed by the TESS mission to date.

preprint2022arXiv

GJ 3929: High Precision Photometric and Doppler Characterization of an Exo-Venus and its Hot, Mini-Neptune-mass Companion

We detail the follow up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b, (TOI-2013b) and its non-transiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star's Venus-zone (P$_{b}$ = 2.616272 $\pm$ 0.000005 days; S$_{b}$ = 17.3$^{+0.8}_{-0.7}$ S$_{\oplus}$) orbiting a nearby M dwarf. GJ 3929c is most likely a non-transiting sub-Neptune. Using the new, ultra-precise NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak National Observatory, we are able to modify the mass constraints of planet b reported in previous works and consequently improve the significance of the mass measurement to almost 4$σ$ confidence (M$_{b}$ = 1.75 $\pm$ 0.45 M$_{\oplus}$). We further adjust the orbital period of planet c from its alias at 14.30 $\pm$ 0.03 days to the likely true period of 15.04 $\pm$ 0.03 days, and adjust its minimum mass to m$\sin i$ = 5.71 $\pm$ 0.92 M$_{\oplus}$. Using the diffuser-assisted ARCTIC imager on the ARC 3.5 m telescope at Apache Point Observatory, in addition to publicly available TESS and LCOGT photometry, we are able to constrain the radius of planet b to R$_{p}$ = 1.09 $\pm$ 0.04 R$_{\oplus}$. GJ 3929b is a top candidate for transmission spectroscopy in its size regime (TSM = 14 $\pm$ 4), and future atmospheric studies of GJ 3929b stand to shed light on the nature of small planets orbiting M dwarfs.

preprint2022arXiv

Leveraging space-based data from the nearest Solar-type star to better understand stellar activity signatures in radial velocity data

Stellar variability is a key obstacle in reaching the sensitivity required to recover Earth-like exoplanetary signals using the radial velocity (RV) detection method. To explore activity signatures in Sun-like stars, we present SolAster, a publicly-distributed analysis pipeline that allows for comparison of space-based measurements with ground-based disk-integrated RVs. Using high spatial resolution Dopplergrams, magnetograms, and continuum filtergrams from the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO), we estimate 'Sun-as-a-star' disk-integrated RVs due to rotationally modulated flux imbalances and convective blueshift suppression, as well as other observables such as unsigned magnetic flux. Comparing these measurements with ground-based RVs from the NEID instrument, which observes the Sun daily using an automated solar telescope, we find a strong relationship between magnetic activity indicators and RV variation, supporting efforts to examine unsigned magnetic flux as a proxy for stellar activity in slowly rotating stars. Detrending against measured unsigned magnetic flux allows us to improve the NEID RV measurements by ~20\% (~50 cm/s in a quadrature sum), yielding an RMS scatter of ~60 cm/s over five months. We also explore correlations between individual and averaged spectral line shapes in the NEID spectra and SDO-derived magnetic activity indicators, motivating future studies of these observables. Finally, applying SolAster to archival planetary transits of Venus and Mercury, we demonstrate the ability to recover small amplitude (< 50 cm/s) RV variations in the SDO data by directly measuring the Rossiter-McLaughlin (RM) signals.

preprint2022arXiv

NEID Rossiter-McLaughlin Measurement of TOI-1268b: A Young Warm Saturn Aligned with Its Cool Host Star

Close-in gas giants present a surprising range of stellar obliquity, the angle between a planet's orbital axis and its host star's spin axis. It is unclear whether the obliquities reflect the planets' dynamical history (e.g., aligned for in situ formation or disk migration versus misaligned for high-eccentricity tidal migration) or whether other mechanisms (e.g., primordial misalignment or planet-star interactions) are more important in sculpting the obliquity distribution. Here we present the stellar obliquity measurement of TOI-1268 (TIC-142394656, $V_{\rm mag} {\sim} 10.9$), a young K-type dwarf hosting an 8.2-day period, Saturn-sized planet. TOI-1268's lithium abundance and rotation period suggest the system age between the ages of Pleiades cluster (${\sim}120$ Myr) and Praesepe cluster (${\sim}670$ Myr). Using the newly commissioned NEID spectrograph, we constrain the stellar obliquity of TOI-1268 via the Rossiter-McLaughlin (RM) effect from both radial velocity (RV) and Doppler Tomography (DT) signals. The 3$σ$ upper bounds of the projected stellar obliquity $|λ|$ from both models are below 60$^\circ$. The large host star separation ($a/R_\star {\sim} 17$), combined with the system's young age, makes it unlikely that the planet has realigned its host star. The stellar obliquity measurement of TOI-1268 probes the architecture of a young gas giant beyond the reach of tidal realignment ($a/R_\star {\gtrsim} 10$) and reveals an aligned or slightly misaligned system.

preprint2022arXiv

Observing the Sun as a star: Design and early results from the NEID solar feed

Efforts with extreme-precision radial velocity (EPRV) instruments to detect small-amplitude planets are largely limited, on many timescales, by the effects of stellar variability and instrumental systematics. One avenue for investigating these effects is the use of small solar telescopes which direct disk-integrated sunlight to these EPRV instruments, observing the Sun at high cadence over months or years. We have designed and built a solar feed system to carry out "Sun-as-a-star" observations with NEID, a very high precision Doppler spectrometer recently commissioned at the WIYN 3.5m Telescope at Kitt Peak National Observatory. The NEID solar feed has been taking observations nearly every day since December 2020; data is publicly available at the NASA Exoplanet Science Institute (NExScI) NEID Solar Archive: \url{https://neid.ipac.caltech.edu/search_solar.php}. In this paper, we present the design of the NEID solar feed and explanations behind our design intent. We also present early radial velocity (RV) results which demonstrate NEID's RV stability on the Sun over 4 months of commissioning: 0.66~m/s RMS under good sky conditions and improving to 0.41~m/s RMS under best conditions.

preprint2022arXiv

Report on Community Cadence Observing to Maximize the Scientific Output of the Keck Planet Finder

The arrival of the Keck Planet Finder (KPF) in 2022 represents a major advance in the precision radial velocity (PRV) capabilities of the W. M. Keck Observatory. In preparation for KPF science, our committee of PRV experts and WMKO staff studied the current implementation of cadence observing at Keck and other PRV facilities. We find that many of KPF's major science cases are not feasible through Keck's standard allocations of full or half nights to individual PIs. Pooling time among several PIs as is currently done by the California Planet Search (CPS) collaboration with HIRES results in lower quality science results than is possible when KPF is available at higher observational cadence. This strategy also creates barriers to entry, particularly for researchers wishing to lead small proposals. This report makes recommendations for optimizing PRV cadence at Keck subject to the following constraints: preservation of clear boundaries between cadence observations and classically scheduled time; and ensuring fairness and scientific independence of different Keck TACs and different KPF PIs. We recommend establishing a new category of Keck time allocation, "KPF Community Cadence" (KPF-CC). In many ways, KPF-CC will formalize observing strategies provided by CPS, but with higher observational cadence appropriate for KPF science and with universal access to the program for all Keck users. We recommend that KPF-CC time be scheduled classically into blocks as small as a quarter night subject to considerations of bright/dark time, variations in proposal pressure with the seasons, and the needs of non-KPF observing programs. Within KPF-CC time, the Keck Observing Assistants would execute observations generated by a dynamic scheduler. We recommend that Keck staff and a board of PRV experts design and maintain the scheduling software.

preprint2022arXiv

Retrieving the C and O Abundances of HR 7672~AB: a Solar-Type Primary Star with a Benchmark Brown Dwarf

A benchmark brown dwarf (BD) is a BD whose properties (e.g., mass and chemical composition) are precisely and independently measured. Benchmark BDs are valuable in testing theoretical evolutionary tracks, spectral synthesis, and atmospheric retrievals for sub-stellar objects. Here, we report results of atmospheric retrieval on a synthetic spectrum and a benchmark BD -- HR 7672~B -- with \petit. First, we test the retrieval framework on a synthetic PHOENIX BT-Settl spectrum with a solar composition. We show that the retrieved C and O abundances are consistent with solar values, but the retrieved C/O is overestimated by 0.13-0.18, which is $\sim$4 times higher than the formal error bar. Second, we perform retrieval on HR 7672~B using high spectral resolution data (R=35,000) from the Keck Planet Imager and Characterizer (KPIC) and near infrared photometry. We retrieve [C/H], [O/H], and C/O to be $-0.24\pm0.05$, $-0.19\pm0.04$, and $0.52\pm0.02$. These values are consistent with those of HR 7672~A within 1.5-$σ$. As such, HR 7672~B is among only a few benchmark BDs (along with Gl 570~D and HD 3651~B) that have been demonstrated to have consistent elemental abundances with their primary stars. Our work provides a practical procedure of testing and performing atmospheric retrieval, and sheds light on potential systematics of future retrievals using high- and low-resolution data.

preprint2022arXiv

Rotational modulation of spectroscopic Zeeman signatures in low-mass stars

Accurate tracers of the stellar magnetic field and rotation are cornerstones for the study of M dwarfs and for reliable detection and characterization of their exoplanetary companions. Such measurements are particularly challenging for old, slowly rotating, fully convective M dwarfs. To explore the use of new activity and rotation tracers, we examined multi-year near-infrared spectroscopic monitoring of two such stars -- GJ 699 (Barnard's Star) and Teegarden's Star -- carried out with Habitable Zone Planet Finder spectrograph. We detected periodic variations in absorption line widths across the stellar spectrum with higher amplitudes towards longer wavelengths. We also detected similar variations in the strength and width of the 12435.67 Angstrom neutral potassium (K I) line, a known tracer of the photospheric magnetic field. Attributing these variations to rotational modulation, we confirm the known $145\pm15$ d rotation period of GJ 699, and measure the rotation period of Teegarden's Star to be $99.6\pm1.4$ d. Based on simulations of the K I line and the wavelength-dependence of the line width signal, we argue that the observed signals are consistent with varying photospheric magnetic fields and the associated Zeeman effect. These results highlight the value of detailed line profile measurements in the near-infrared for diagnosing stellar magnetic field variability. Such measurements may be pivotal for disentangling activity and exoplanet-related signals in spectroscopic monitoring of old, low-mass stars.

preprint2022arXiv

TESS-Keck Survey IX: Masses of Three Sub-Neptunes Orbiting HD 191939 and the Discovery of a Warm Jovian Plus a Distant Sub-Stellar Companion

Exoplanet systems with multiple transiting planets are natural laboratories for testing planetary astrophysics. One such system is HD 191939 (TOI-1339), a bright (V=9) and Sun-like (G9V) star, which TESS found to host three transiting planets (b, c, and d). The planets have periods of 9, 29, and 38 days each with similar sizes from 3 to 3.4 $R_{\oplus}$. To further characterize the system, we measured the radial velocity (RV) of HD 191939 over 415 days with Keck/HIRES and APF/Levy. We find that $M_b = 10.4 \pm 0.9 M_{\oplus}$ and $M_c = 7.2 \pm 1.4 M_{\oplus}$, which are low compared to most known planets of comparable radii. The RVs yield only an upper-limit on $M_d$ (<5.8 $M_{\oplus}$ at 2$σ$). The RVs further reveal a fourth planet (e) with a minimum mass of $0.34 \pm 0.01 M_{Jup}$ and an orbital period of 101.4 $\pm$ 0.4 days. Despite its non-transiting geometry, secular interactions between planet e and the inner transiting planets indicate that planet e is coplanar with the transiting planets ($Δ$i < 10$^{\circ}$). We identify a second high mass planet (f) with 95% confidence intervals on mass between $2-11 \, M_{Jup}$ and period between 1700-7200 days, based on a joint analysis of RVs and astrometry from $Gaia$ and $Hipparcos$. As a bright star hosting multiple planets with well-measured masses, HD 191939 presents many options for comparative planetary astronomy including characterization with JWST.

preprint2022arXiv

The Aligned Orbit of WASP-148b, the Only Known Hot Jupiter with a Nearby Warm Jupiter Companion, from NEID and HIRES

We present spectroscopic measurements of the Rossiter-McLaughlin effect for WASP-148b, the only known hot Jupiter with a nearby warm-Jupiter companion, from the WIYN/NEID and Keck/HIRES instruments. This is one of the first scientific results reported from the newly commissioned NEID spectrograph, as well as the second obliquity constraint for a hot Jupiter system with a close-in companion, after WASP-47. WASP-148b is consistent with being in alignment with the sky-projected spin axis of the host star, with $λ=-8^{\circ}.2^{{+8^{\circ}.7}}_{-9^{\circ}.7}$. The low obliquity observed in the WASP-148 system is consistent with the orderly-alignment configuration of most compact multi-planet systems around cool stars with obliquity constraints, including our solar system, and may point to an early history for these well-organized systems in which migration and accretion occurred in isolation, with relatively little disturbance. By contrast, previous results have indicated that high-mass and hot stars appear to more commonly host a wide range of misaligned planets: not only single hot Jupiters, but also compact systems with multiple super-Earths. We suggest that, to account for the high rate of spin-orbit misalignments in both compact multi-planet and isolated-hot-Jupiter systems orbiting high-mass and hot stars, spin-orbit misalignments may be caused by distant giant planet perturbers, which are most common around these stellar types.

preprint2022arXiv

The TESS-Keck Survey. VIII. Confirmation of a Transiting Giant Planet on an Eccentric 261 day Orbit with the Automated Planet Finder Telescope

We report the discovery of TOI-2180 b, a 2.8 $M_{\rm J}$ giant planet orbiting a slightly evolved G5 host star. This planet transited only once in Cycle 2 of the primary Transiting Exoplanet Survey Satellite (TESS) mission. Citizen scientists identified the 24 hr single-transit event shortly after the data were released, allowing a Doppler monitoring campaign with the Automated Planet Finder telescope at Lick Observatory to begin promptly. The radial velocity observations refined the orbital period of TOI-2180 b to be 260.8$\pm$0.6 days, revealed an orbital eccentricity of 0.368$\pm$0.007, and discovered long-term acceleration from a more distant massive companion. We conducted ground-based photometry from 14 sites spread around the globe in an attempt to detect another transit. Although we did not make a clear transit detection, the nondetections improved the precision of the orbital period. We predict that TESS will likely detect another transit of TOI-2180 b in Sector 48 of its extended mission. We use giant planet structure models to retrieve the bulk heavy-element content of TOI-2180 b. When considered alongside other giant planets with orbital periods over 100 days, we find tentative evidence that the correlation between planet mass and metal enrichment relative to stellar is dependent on orbital properties. Single-transit discoveries like TOI-2180 b highlight the exciting potential of the TESS mission to find planets with long orbital periods and low irradiation fluxes despite the selection biases associated with the transit method.

preprint2022arXiv

The TESS-Keck Survey. XI. Mass Measurements for Four Transiting sub-Neptunes orbiting K dwarf TOI-1246

Multi-planet systems are valuable arenas for investigating exoplanet architectures and comparing planetary siblings. TOI-1246 is one such system, with a moderately bright K dwarf ($\rm{V=11.6,~K=9.9}$) and four transiting sub-Neptunes identified by TESS with orbital periods of $4.31~\rm{d},~5.90~\rm{d},~18.66~\rm{d}$, and $~37.92~\rm{d}$. We collected 130 radial velocity observations with Keck/HIRES and TNG/HARPS-N to measure planet masses. We refit the 14 sectors of TESS photometry to refine planet radii ($\rm{2.97 \pm 0.06~R_\oplus},\rm{2.47 \pm 0.08~R_\oplus}, \rm{3.46 \pm 0.09~R_\oplus}$, $\rm{3.72 \pm 0.16~R_\oplus}$), and confirm the four planets. We find that TOI-1246 e is substantially more massive than the three inner planets ($\rm{8.1 \pm 1.1 M_\oplus}$, $\rm{8.8 \pm 1.2 M_\oplus}$, $\rm{5.3 \pm 1.7 M_\oplus}$, $\rm{14.8 \pm 2.3 M_\oplus}$). The two outer planets, TOI-1246 d and TOI-1246 e, lie near to the 2:1 resonance ($\rm{P_{e}/P_{d}=2.03}$) and exhibit transit timing variations. TOI-1246 is one of the brightest four-planet systems, making it amenable for continued observations. It is one of only six systems with measured masses and radii for all four transiting planets. The planet densities range from $\rm{0.70 \pm 0.24}$ to $3.21 \pm 0.44 \rm{g/cm^3}$, implying a range of bulk and atmospheric compositions. We also report a fifth planet candidate found in the RV data with a minimum mass of 25.6 $\pm$ 3.6 $\rm{M_\oplus}$. This planet candidate is exterior to TOI-1246 e with a candidate period of 93.8 d, and we discuss the implications if it is confirmed to be planetary in nature.

preprint2022arXiv

The TESS-Keck Survey. XIII. An Eccentric Hot Neptune with a Similar-Mass Outer Companion around TOI-1272

We report the discovery of an eccentric hot Neptune and a non-transiting outer planet around TOI-1272. We identified the eccentricity of the inner planet, with an orbital period of 3.3 d and $R_{\rm p,b} = 4.1 \pm 0.2$ $R_\oplus$, based on a mismatch between the observed transit duration and the expected duration for a circular orbit. Using ground-based radial velocity measurements from the HIRES instrument at the Keck Observatory, we measured the mass of TOI-1272b to be $M_{\rm p,b} = 25 \pm 2$ $M_\oplus$. We also confirmed a high eccentricity of $e_b = 0.34 \pm 0.06$, placing TOI-1272b among the most eccentric well-characterized sub-Jovians. We used these RV measurements to also identify a non-transiting outer companion on an 8.7-d orbit with a similar mass of $M_{\rm p,c}$ sin$i= 27 \pm 3$ $M_\oplus$ and $e_c \lesssim 0.35$. Dynamically stable planet-planet interactions have likely allowed TOI-1272b to avoid tidal eccentricity decay despite the short circularization timescale expected for a close-in eccentric Neptune. TOI-1272b also maintains an envelope mass fraction of $f_{\rm env} \approx 11\%$ despite its high equilibrium temperature, implying that it may currently be undergoing photoevaporation. This planet joins a small population of short-period Neptune-like planets within the "Hot Neptune Desert" with a poorly understood formation pathway.

preprint2022arXiv

The Warm Neptune GJ 3470b has a Polar Orbit

The warm Neptune GJ 3470b transits a nearby ($d=29$pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5m Telescope at Kitt Peak Observatory, we model the classical Rossiter-Mclaughlin effect yielding a sky-projected obliquity of $λ=98_{-12}^{+15\:\circ}$ and a $v \sin i = 0.85_{-0.33}^{+0.27}$km/s. Leveraging information about the rotation period and size of the host star, our analysis yields a true obliquity of $ψ=95_{-8}^{+9\:\circ}$, revealing that GJ 3470b is on a polar orbit. Using radial velocities from HIRES, HARPS and the Habitable-zone Planet Finder, we show that the data are compatible with a long-term RV slope of $\dotγ = -0.0022 \pm 0.0011$m/s/day over a baseline of 12.9 years. If the RV slope is due to acceleration from another companion in the system, we show that such a companion is capable of explaining the polar and mildly eccentric orbit of GJ 3470b using two different secular excitation models. The existence of an outer companion can be further constrained with additional RV observations, Gaia astrometry, and future high-contrast imaging observations. Lastly, we show that tidal heating from GJ 3470b's mild eccentricity has most likely inflated the radius of GJ 3470b by a factor of $\sim$1.5-1.7, which could help account for its evaporating atmosphere.

preprint2022arXiv

TKS V. Twin sub-Neptunes Transiting the Nearby G Star HD 63935

We present the discovery of two nearly identically-sized sub-Neptune transiting planets orbiting HD 63935, a bright ($V=8.6$ mag), sun-like ($T_{eff}=5560K$) star at 49 pc. TESS identified the first planet, HD 63935 b (TOI-509.01), in Sectors 7 and 34. We identified the second signal (HD 63935 c) in Keck HIRES and Lick APF radial velocity data as part of our followup campaign. It was subsequently confirmed with TESS photometry in Sector 34 as TOI-509.02. Our analysis of the photometric and radial velocity data yields a robust detection of both planets with periods of $9.0600 \pm 0.007$ and $21.40 \pm 0.0019$ days, radii of $2.99 \pm 0.14$ and $2.90 \pm 0.13$ $R_\oplus$, and masses of $10.8 \pm 1.8$ and $11.1 \pm 2.4$ $M_\oplus$. We calculate densities for planets b and c consistent with a few percent of the planet mass in hydrogen/helium envelopes. We also describe our survey's efforts to choose the best targets for JWST atmospheric followup. These efforts suggest that HD 63935 b will have the most clearly visible atmosphere of its class. It is the best target for transmission spectroscopy (ranked by Transmission Spectroscopy Metric, a proxy for atmospheric observability) in the so-far uncharacterized parameter space comprising sub-Neptune-sized (2.6 $R_\oplus$ $<$ $R_p$ $<$ 4 $R_\oplus$), moderately-irradiated (100 $F_\oplus$ $<$ $F_p$ $<$ 1000 $F_\oplus$) planets around G-stars. Planet c is also a viable target for transmission spectroscopy, and given the indistinguishable masses and radii of the two planets, the system serves as a natural laboratory for examining the processes that shape the evolution of sub-Neptune planets.

preprint2022arXiv

TOI-1696 and TOI-2136: Constraining the Masses of Two Mini-Neptunes with HPF

We present the validation of two planets orbiting M dwarfs, TOI-1696b and TOI-2136b. Both planets are mini-Neptunes orbiting nearby stars, making them promising prospects for atmospheric characterization with the James Webb Space Telescope. We validated the planetary nature of both candidates using high contrast imaging, ground-based photometry, and near-infrared radial velocities. Adaptive Optics images were taken using the ShARCS camera on the 3 m Shane Telescope. Speckle images were taken using the NN-Explore Exoplanet Stellar Speckle Imager on the WIYN 3.5 m telescope. Radii and orbital ephemerides were refined using a combination of TESS, the diffuser-assisted ARCTIC imager on the 3.5m ARC telescope at Apache Point Observatory, and the 0.6 m telescope at Red Buttes Observatory. We obtained radial velocities using the Habitable-Zone Planet Finder on the 10 m Hobby-Eberly Telescope, which enabled us to place upper limits on the masses of both transiting planets. TOI-1696b (P = 2.5 days; R$_{p}$ = 3.24 R$_{\oplus}$; M$_{p}$ $<$ 56.6 M$_{\oplus}$) falls into a sparsely-populated region of parameter space considering its host star's temperature (T$_{\rm{eff}}$ = 3168 K, M4.5), as planets of its size are quite rare around mid to late M dwarfs. On the other hand, TOI-2136b (P = 7.85 days; R$_{p}$ = 2.09 R$_{\oplus}$; M$_{p}$ $<$ 15.0 M$_{\oplus}$) is an excellent candidate for atmospheric follow-up with JWST.

preprint2022arXiv

TOI-3714 b and TOI-3629 b: Two gas giants transiting M dwarfs confirmed with HPF and NEID

We confirm the planetary nature of two gas giants discovered by TESS to transit M dwarfs. TOI-3714 ($V=15.24,~J=11.74$) is an M2 dwarf hosting a hot Jupiter ($M_p=0.70 \pm 0.03~\mathrm{M_J}$ and $R_p=1.01 \pm 0.03~\mathrm{R_J}$) on an orbital period of $2.154849 \pm 0.000001$ days with a resolved white dwarf companion. TOI-3629 ($V=14.63,~J=11.42$) is an M1 dwarf hosting a hot Jupiter ($M_p=0.26 \pm 0.02~\mathrm{M_J}$ and $R_p=0.74 \pm 0.02~\mathrm{R_J}$) on an orbital period of $3.936551_{-0.000006}^{+0.000005}$ days. We characterize each transiting companion using a combination of ground-based and space-based photometry, speckle imaging, and high-precision velocimetry from the Habitable-zone Planet Finder and the NEID spectrographs. With the discovery of these two systems, there are now nine M dwarfs known to host transiting hot Jupiters. Among this population, TOI-3714 b ($T_{eq}=750\pm20$ K and $\mathrm{TSM}=98\pm7$) and TOI-3629 b ($T_{eq}=690\pm20$ K and $\mathrm{TSM}=80\pm9$) are warm gas giants amenable to additional characterization with transmission spectroscopy to probe atmospheric chemistry and, for TOI-3714, obliquity measurements to probe formation scenarios.

preprint2022arXiv

TOI-3757 b: A low density gas giant orbiting a solar-metallicity M dwarf

We present the discovery of a new Jovian-sized planet, TOI-3757 b, the lowest density planet orbiting an M dwarf (M0V). It orbits a solar-metallicity M dwarf discovered using TESS photometry and confirmed with precise radial velocities (RV) from HPF and NEID. With a planetary radius of $12.0^{+0.4}_{-0.5}$ $R_{\oplus}$ and mass of $85.3^{+8.8}_{-8.7}$ $M_{\oplus}$, not only does this object add to the small sample of gas giants ($\sim 10$) around M dwarfs, but also, its low density ($ρ=$ $0.27^{+0.05}_{-0.04}$ $\textrm{g~cm}^{-3}$) provides an opportunity to test theories of planet formation. We present two hypotheses to explain its low density; first, we posit that the low metallicity of its stellar host ($\sim$ 0.3 dex lower than the median metallicity of M dwarfs hosting gas giants) could have played a role in the delayed formation of a solid core massive enough to initiate runaway accretion. Second, using the eccentricity estimate of $0.14 \pm 0.06$ we determine it is also plausible for tidal heating to at least partially be responsible for inflating the radius of TOI-3757b b. The low density and large scale height of TOI-3757 b makes it an excellent target for transmission spectroscopy studies of atmospheric escape and composition (TSM $\sim$ 190). We use HPF to perform transmission spectroscopy of TOI-3757 b using the helium 10830 Å~ line. Doing this, we place an upper limit of 6.9 \% (with 90\% confidence) on the maximum depth of the absorption from the metastable transition of He at $\sim$ 10830 Å, which can help constraint the atmospheric mass loss rate in this energy limited regime.

preprint2022arXiv

Validation of 13 Hot and Potentially Terrestrial TESS Planets

The James Webb Space Telescope (JWST) will be able to probe the atmospheres and surface properties of hot, terrestrial planets via emission spectroscopy. We identify 18 potentially terrestrial planet candidates detected by the Transiting Exoplanet Survey Satellite (TESS) that would make ideal targets for these observations. These planet candidates cover a broad range of planet radii ($R_{\rm p} \sim 0.6 - 2.0 R_\oplus$) and orbit stars of various magnitudes ($K_s = 5.78 - 10.78$, $V = 8.4 - 15.69$) and effective temperatures ($T_{\rm eff }\sim 3000 - 6000$ K). We use ground-based observations collected through the TESS Follow-up Observing Program (TFOP) and two vetting tools -- DAVE and TRICERATOPS -- to assess the reliabilities of these candidates as planets. We validate 13 planets: TOI-206 b, TOI-500 b, TOI-544 b, TOI-833 b, TOI-1075 b, TOI-1411 b, TOI-1442 b, TOI-1693 b, TOI-1860 b, TOI-2260 b, TOI-2411 b, TOI-2427 b, and TOI-2445 b. Seven of these planets (TOI-206 b, TOI-500 b, TOI-1075 b, TOI-1442 b, TOI-2260 b, TOI-2411 b, and TOI-2445 b) are ultra-short-period planets. TOI-1860 is the youngest ($133 \pm 26$ Myr) solar twin with a known planet to date. TOI-2260 is a young ($321 \pm 96$ Myr) G dwarf that is among the most metal-rich ([Fe/H] = $0.22 \pm 0.06$ dex) stars to host an ultra-short-period planet. With an estimated equilibrium temperature of $\sim 2600$ K, TOI-2260 b is also the fourth hottest known planet with $R_{\rm p} < 2 \, R_\oplus$.

preprint2021arXiv

Gaia 20eae: A newly discovered episodically accreting young star

The Gaia Alert System issued an alert on 2020 August 28, on Gaia 20eae when its light curve showed a $\sim$4.25 magnitude outburst. We present multi-wavelength photometric and spectroscopic follow-up observations of this source since 2020 August and identify it as the newest member of the FUor/EXor family of sources. We find that the present brightening of Gaia 20eae is not due to the dust clearing event but due to an intrinsic change in the spectral energy distribution. The light curve of Gaia 20eae shows a transition stage during which most of its brightness ($\sim$3.4 mag) has occurred at a short timescale of 34 days with a rise-rate of 3 mag/month. Gaia 20eae has now started to decay at a rate of 0.3 mag/month. We have detected a strong P Cygni profile in H$α$ which indicates the presence of winds originating from regions close to the accretion. We find signatures of very strong and turbulent outflow and accretion in Gaia 20eae during this outburst phase. We have also detected a red-shifted absorption component in all the Ca II IR triplet lines consistent with signature of hot in-falling gas in the magnetospheric accretion funnel. This enables us to constrain the viewing angle with respect to the accretion funnel. Our investigation of Gaia 20eae points towards magnetospheric accretion being the phenomenon for the current outburst.

preprint2021arXiv

High resolution near-infrared spectroscopy of a flare around the ultracool dwarf vB 10

We present high-resolution observations of a flaring event in the M8 dwarf vB 10 using the near-infrared Habitable zone Planet Finder (HPF) spectrograph on the Hobby Eberly Telescope (HET). The high stability of HPF enables us to accurately subtract a VB 10 quiescent spectrum from the flare spectrum to isolate the flare contributions, and study the changes in the relative energy of the Ca II infrared triplet (IRT), several Paschen lines, the He 10830 Å~ triplet lines, and select iron and magnesium lines in HPF`s bandpass. Our analysis reveals the presence of a red asymmetry in the He 10830 Å~ triplet; which is similar to signatures of coronal rain in the Sun. Photometry of the flare derived from an acquisition camera before spectroscopic observations, and the ability to extract spectra from up-the-ramp observations with the HPF infrared detector, enables us to perform time-series analysis of part of the flare, and provide coarse constraints on the energy and frequency of such flares. We compare this flare with historical observations of flares around vB 10 and other ultracool M dwarfs, and attempt to place limits on flare-induced atmospheric mass loss for hypothetical planets around vB 10.

preprint2021arXiv

On the origin of nitrogen at low metallicity

Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

preprint2021arXiv

Target Prioritization and Observing Strategies for the NEID Earth Twin Survey

NEID is a high-resolution optical spectrograph on the WIYN 3.5-m telescope at Kitt Peak National Observatory and will soon join the new generation of extreme precision radial velocity instruments in operation around the world. We plan to use the instrument to conduct the NEID Earth Twin Survey (NETS) over the course of the next 5 years, collecting hundreds of observations of some of the nearest and brightest stars in an effort to probe the regime of Earth-mass exoplanets. Even if we take advantage of the extreme instrumental precision conferred by NEID, it will remain difficult to disentangle the weak (~10 cm s$^{-1}$) signals induced by such low-mass, long-period exoplanets from stellar noise for all but the quietest host stars. In this work, we present a set of quantitative selection metrics which we use to identify an initial NETS target list consisting of stars conducive to the detection of exoplanets in the regime of interest. We also outline a set of observing strategies with which we aim to mitigate uncertainty contributions from intrinsic stellar variability and other sources of noise.

preprint2021arXiv

The TESS-Keck Survey IV: A Retrograde, Polar Orbit for the Ultra-Low-Density, Hot Super-Neptune WASP-107b

We measured the Rossiter-McLaughlin effect of WASP-107b during a single transit with Keck/HIRES. We found the sky-projected inclination of WASP-107b's orbit, relative to its host star's rotation axis, to be $|λ| = {118}^{+38}_{-19}$ degrees. This confirms the misaligned/polar orbit that was previously suggested from spot-crossing events and adds WASP-107b to the growing population of hot Neptunes in polar orbits around cool stars. WASP-107b is also the fourth such planet to have a known distant planetary companion. We examined several dynamical pathways by which this companion could have induced such an obliquity in WASP-107b. We find that nodal precession and disk dispersal-driven tilting can both explain the current orbital geometry while Kozai-Lidov cycles are suppressed by general relativity. While each hypothesis requires a mutual inclination between the two planets, nodal precession requires a much larger angle which for WASP-107 is on the threshold of detectability with future Gaia astrometric data. As nodal precession has no stellar type dependence, but disk dispersal-driven tilting does, distinguishing between these two models is best done on the population level. Finding and characterizing more extrasolar systems like WASP-107 will additionally help distinguish whether the distribution of hot-Neptune obliquities is a dichotomy of aligned and polar orbits or if we are uniformly sampling obliquities during nodal precession cycles.

preprint2021arXiv

The TESS-Keck Survey: Science Goals and Target Selection

Space-based transit missions such as Kepler and TESS have demonstrated that planets are ubiquitous. However, the success of these missions heavily depends on ground-based radial velocity (RV) surveys, which combined with transit photometry can yield bulk densities and orbital properties. While most Kepler host stars are too faint for detailed follow-up observations, TESS is detecting planets orbiting nearby bright stars that are more amenable to RV characterization. Here we introduce the TESS-Keck Survey (TKS), an RV program using ~100 nights on Keck/HIRES to study exoplanets identified by TESS. The primary survey aims are investigating the link between stellar properties and the compositions of small planets; studying how the diversity of system architectures depends on dynamical configurations or planet multiplicity; identifying prime candidates for atmospheric studies with JWST; and understanding the role of stellar evolution in shaping planetary systems. We present a fully-automated target selection algorithm, which yielded 103 planets in 86 systems for the final TKS sample. Most TKS hosts are inactive, solar-like, main-sequence stars (4500 K < Teff < 6000 K) at a wide range of metallicities. The selected TKS sample contains 71 small planets (Rp < 4 Re), 11 systems with multiple transiting candidates, 6 sub-day period planets and 3 planets that are in or near the habitable zone of their host star. The target selection described here will facilitate the comparison of measured planet masses, densities, and eccentricities to predictions from planet population models. Our target selection software is publicly available (at https://github.com/ashleychontos/sort-a-survey) and can be adapted for any survey which requires a balance of multiple science interests within a given telescope allocation.

preprint2020arXiv

A Mini-Neptune and a Venus-Zone Planet in the Radius Valley Orbiting the Nearby M2-dwarf TOI-1266: Validation with the Habitable-zone Planet Finder

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.

preprint2020arXiv

A warm Jupiter transiting an M dwarf: A TESS single transit event confirmed with the Habitable-zone Planet Finder

We confirm the planetary nature of a warm Jupiter transiting the early M dwarf TOI-1899, using a combination of available TESS photometry; high-precision, near-infrared spectroscopy with the Habitable-zone Planet Finder; and speckle and adaptive optics imaging. The data reveal a transiting companion on an $\sim29$-day orbit with a mass and radius of $0.66\pm0.07\ \mathrm{M_{J}}$ and $1.15_{-0.05}^{+0.04}\ \mathrm{R_{J}}$, respectively. The star TOI-1899 is the lowest-mass star known to host a transiting warm Jupiter, and we discuss the follow-up opportunities afforded by a warm ($\mathrm{T_{eq}}\sim362$ K) gas giant orbiting an M0 star. Our observations reveal that TOI-1899.01 is a puffy warm Jupiter, and we suggest additional transit observations to both refine the orbit and constrain the true dilution observed in TESS.

preprint2020arXiv

Evidence for He I 10830 Å~ absorption during the transit of a warm Neptune around the M-dwarf GJ 3470 with the Habitable-zone Planet Finder

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 Å~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 Å~in absorption (equivalent width $\sim$ $0.012 \pm 0.002$ Å) 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 Å\, 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.

preprint2020arXiv

Helium and Nitrogen Enrichment in Massive Main Sequence Stars: Mechanisms and Implications for the Origin of WNL Stars

The evolutionary paths taken by massive stars with $M \gtrsim 60 \, \mathrm{M}_\odot$ remain substantially uncertain. They begin their lives as main sequence (MS) O-stars. Depending on their masses, rotation rates, and metallicities, they can then encounter a wide range of evolutionary states with an equally broad set of possible surface compositions and spectral classifications. We present a new grid of calculations for the evolution of such stars that covers a broad range in mass, M/M$_\odot = 60$ to $150$, rotation rate, $v \, / \, v_{\rm crit} = 0$ to $0.6$, metallicity, $[\mathrm{Fe}/\mathrm{H}] = -4$ to $0$, and $α$-element enhancement, $[α/\mathrm{Fe}] = 0$ to $0.4$. We show that rotating stars undergo rotationally-induced dredge-up of nucleosynthetic products, mostly He and N, to their surfaces while still on the MS. Non-rotating metal-rich stars also reveal the products of nucleosynthesis on their surfaces because even modest amounts of mass loss expose their "fossil" convective cores: regions that are no longer convective, but which were part of the convective core at an early stage in the star's evolution. Thus surface enhancement of He and N is expected for rotating stars at all metallicities, and for non-rotating stars if they are relatively metal-rich. We calculate a stellar atmosphere for a representative model from our grid, properly accounting for He- and N-enhancement, and show that the resulting spectrum provides a good match to observed WNL stars, strongly suggesting that the physical mechanisms we have identified are the ultimate cause of the WNL phase.

preprint2020arXiv

Persistent starspot signals on M dwarfs: multi-wavelength Doppler observations with the Habitable-zone Planet Finder and Keck/HIRES

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.

preprint2020arXiv

Solar Contamination in Extreme Precision Radial Velocity Measurements: Deleterious Effects and Prospects for Mitigation

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.

preprint2020arXiv

The multi-planet system TOI-421 -- A warm Neptune and a super puffy mini-Neptune transiting a G9 V star in a visual binary

We report the discovery of a warm Neptune and a hot sub-Neptune transiting TOI-421 (BD-14 1137, TIC 94986319), a bright (V=9.9) G9 dwarf star in a visual binary system observed by the TESS space mission in Sectors 5 and 6. We performed ground-based follow-up observations -- comprised of LCOGT transit photometry, NIRC2 adaptive optics imaging, and FIES, CORALIE, HARPS, HIRES, and PFS high-precision Doppler measurements -- and confirmed the planetary nature of the 16-day transiting candidate announced by the TESS team. We discovered an additional radial velocity signal with a period of 5 days induced by the presence of a second planet in the system, which we also found to transit its host star. We found that the inner mini-Neptune, TOI-421b, has an orbital period of Pb =5.19672 +- 0.00049 days, a mass of Mb = 7.17 +- 0.66 Mearth and a radius of Rb = 2.68+0.19-0.18 Rearth, whereas the outer warm Neptune, TOI-421 c, has a period of Pc =16.06819 +- 0.00035 days, a mass of Mc = 16.42+1.06-1.04 Mearth, a radius of Rc = 5.09+0.16-0.15 Rearth and a density of rho_c =0.685+0.080-0.072 g cm-3 . With its characteristics the inner planet (rho_b=2.05+0.52-0.41 g cm-3) is placed in the intriguing class of the super-puffy mini-Neptunes. TOI-421b and TOI-421c are found to be well suitable for atmospheric characterization. Our atmospheric simulations predict significant Ly-alpha transit absorption, due to strong hydrogen escape in both planets, and the presence of detectable CH_4 in the atmosphere of TOI-421c if equilibrium chemistry is assumed.

preprint2020arXiv

The SDSS-HET Survey of Kepler Eclipsing Binaries. Description of the Survey and First Results

The Kepler mission has provided a treasure trove of eclipsing binaries (EBs), observed at extremely high photometric precision, nearly continuously for several years. We are carrying out a survey of ~100 of these EBs to derive dynamical masses and radii with precisions of 3% or better. We use multiplexed near-infrared H band spectroscopy from the SDSS-III and -IV APOGEE instrument and optical spectroscopy from the HET High-Resolution Spectrograph to derive double-lined spectroscopic orbits and dynamical mass-ratios for the EB sample. This information is combined with Kepler photometry to derive orbital inclination, dynamical masses of the components, radii and temperatures. These measurements are directly applicable for benchmarking stellar models that are integrating the next generation of improvements, such as the magnetic suppression of convection efficiency, updated opacity tables, and fine-tuned equations of state. We selected our EB sample to include systems with low-mass ($M \lt 0.8 \; M_{\odot}$) primary or secondary components, many expected to populate the relatively sparse parameter space below $0.5 \; M_{\odot}$. In this paper, we describe our EB sample and the analysis techniques we are utilizing, and present masses and radii for two systems that inhabit particularly underpopulated regions of mass-radius-period space: KIC 2445134 and KIC 3003991. Our analysis of KIC 2445134 ($q=0.411 \pm 0.001$) yields masses and radii of $M_A = 1.29 \pm 0.03 \; M_{\odot}$, $M_B = 0.53 \pm 0.01 \; M_{\odot}$, $R_A = 1.42 \pm 0.01 \; R_{\odot}$, $R_B = 0.510 \pm 0.004 \; R_{\odot}$, and a temperature ratio of $T_B/T_A = 0.635 \pm 0.001$; our analysis of KIC 3003991 ($q=0.298 \pm 0.006$) yields $M_A = 0.74 \pm \; 0.04 M_{\odot}$, $M_B = 0.222 \pm \; 0.007 M_{\odot}$, $R_A = 0.84 \pm 0.01 \; R_{\odot}$, $R_B = 0.250 \pm 0.004 \; R_{\odot}$, and a temperature ratio of $T_B/T_A= 0.662 \pm 0.001$.

preprint2020arXiv

The TESS-Keck Survey I: A Warm Sub-Saturn-mass Planet and a Caution about Stray Light in TESS Cameras

We report the detection of a Saturn-size exoplanet orbiting HD 332231 (TOI 1456) in light curves from the Transiting Exoplanet Survey Satellite (TESS). HD 332231, an F8 dwarf star with a V-band magnitude of 8.56, was observed by TESS in Sectors 14 and 15. We detect a single-transit event in the Sector 15 presearch data conditioning (PDC) light curve. We obtain spectroscopic follow-up observations of HD 332231 with the Automated Planet Finder, Keck I, and SONG telescopes. The orbital period we infer from the radial velocity (RV) observations leads to the discovery of another transit in Sector 14 that was masked by PDC due to scattered light contamination. A joint analysis of the transit and RV data confirms the planetary nature of HD 332231 b, a Saturn-size ($0.867^{+0.027}_{-0.025} \; R_{\rm J}$), sub-Saturn-mass ($0.244\pm0.021 \; M_{\rm J}$) exoplanet on a 18.71 day circular orbit. The low surface gravity of HD 332231 b and the relatively low stellar flux it receives make it a compelling target for transmission spectroscopy. Also, the stellar obliquity is likely measurable via the Rossiter-McLaughlin effect, an exciting prospect given the 0.14 au orbital separation of HD 332231 b. The spectroscopic observations do not provide substantial evidence for any additional planets in the HD 332231 system, but continued RV monitoring is needed to further characterize this system. We also predict that the frequency and duration of masked data in the PDC light curves for TESS Sectors 14-16 could hide transits of some exoplanets with orbital periods between 10.5 and 17.5 days.

preprint2020arXiv

The TESS-Keck Survey II: An Ultra-Short Period Rocky Planet and its Siblings Transiting the Galactic Thick-Disk Star TOI-561

We report the discovery of TOI-561, a multi-planet system in the galactic thick disk that contains a rocky, ultra-short period planet (USP). This bright ($V=10.2$) star hosts three small transiting planets identified in photometry from the NASA TESS mission: TOI-561 b (TOI-561.02, P=0.44 days, $R_b = 1.45\pm0.11\,R_\oplus$), c (TOI-561.01, P=10.8 days, $R_c=2.90\pm0.13\,R_\oplus$), and d (TOI-561.03, P=16.3 days, $R_d=2.32\pm0.16\,R_\oplus$). The star is chemically ([Fe/H]$=-0.41\pm0.05$, [$α$/H]$=+0.23\pm0.05$) and kinematically consistent with the galactic thick disk population, making TOI-561 one of the oldest ($10\pm3\,$Gyr) and most metal-poor planetary systems discovered yet. We dynamically confirm planets b and c with radial velocities from the W. M. Keck Observatory High Resolution Echelle Spectrometer. Planet b has a mass and density of $3.2\pm0.8\,M_\oplus$ and $5.5^{+2.0}_{-1.6}\,$g$\,$cm$^{-3}$, consistent with a rocky composition. Its lower-than-average density is consistent with an iron-poor composition, although an Earth-like iron-to-silicates ratio is not ruled out. Planet c is $7.0\pm2.3\,M_\oplus$ and $1.6\pm0.6\,$g$\,$cm$^{-3}$, consistent with an interior rocky core overlaid with a low-mass volatile envelope. Several attributes of the photometry for planet d (which we did not detect dynamically) complicate the analysis, but we vet the planet with high-contrast imaging, ground-based photometric follow-up and radial velocities. TOI-561 b is the first rocky world around a galactic thick-disk star confirmed with radial velocities and one of the best rocky planets for thermal emission studies.

preprint2020arXiv

TOI-1235 b: a keystone super-Earth for testing radius valley emergence models around early M dwarfs

Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally-driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/non-rocky transition in period-radius space. Here we present the confirmation of TOI-1235 b ($P=3.44$ days, $r_p=1.738^{+0.087}_{-0.076}$ R$_{\oplus}$), a planet whose size and period are intermediate between the competing model predictions thus making the system an important test case for emergence models of the rocky/non-rocky transition around early M dwarfs ($R_s=0.630\pm 0.015$ R$_{\odot}$, $M_s=0.640\pm 0.016$ M$_{\odot}$). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial-velocities from HARPS-N and HIRES. We measure a planet mass of $6.91^{+0.75}_{-0.85}$ M$_{\oplus}$, which implies an iron core mass fraction of $20^{+15}_{-12}$% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like and we constrain a H/He envelope mass fraction to be $<0.5$% at 90% confidence. Our results are consistent with model predictions from thermally-driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remain efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically-determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin ($P=21.8^{+0.9}_{-0.8}$ days, $m_p\sin{i}=13.0^{+3.8}_{-5.3}$ M$_{\oplus}$) that cannot be firmly ruled out by our data.

preprint2020arXiv

TOI-1728b: The Habitable-zone Planet Finder confirms a warm super Neptune orbiting an M dwarf host

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 Å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 Å~ 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.

preprint2019arXiv

A sub-Neptune sized planet transiting the M2.5-dwarf G 9-40: Validation with the Habitable-zone Planet Finder

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.

preprint2016arXiv

A comprehensive radial velocity error budget for next generation Doppler spectrometers

We describe a detailed radial velocity error budget for the NASA-NSF Extreme Precision Doppler Spectrometer instrument concept NEID (NN-explore Exoplanet Investigations with Doppler spectroscopy). Such an instrument performance budget is a necessity for both identifying the variety of noise sources currently limiting Doppler measurements, and estimating the achievable performance of next generation exoplanet hunting Doppler spectrometers. For these instruments, no single source of instrumental error is expected to set the overall measurement floor. Rather, the overall instrumental measurement precision is set by the contribution of many individual error sources. We use a combination of numerical simulations, educated estimates based on published materials, extrapolations of physical models, results from laboratory measurements of spectroscopic subsystems, and informed upper limits for a variety of error sources to identify likely sources of systematic error and construct our global instrument performance error budget. While natively focused on the performance of the NEID instrument, this modular performance budget is immediately adaptable to a number of current and future instruments. Such an approach is an important step in charting a path towards improving Doppler measurement precisions to the levels necessary for discovering Earth-like planets.

preprint2016arXiv

A Versatile Technique to Enable sub-milli-Kelvin Instrument Stability for Precise Radial Velocity Measurements: Tests with the Habitable-zone Planet Finder

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.

preprint2016arXiv

Detection of a very low mass star in an Eclipsing Binary system

We report the detection of a very low mass star (VLMS) companion to the primary star 1SWASPJ234318.41+295556.5A (J2343+29A), using radial velocity (RV) measurements from the PARAS (PRL Advanced Radial-velocity Abu-sky Search) high resolution echelle spectrograph. The periodicity of the single-lined eclipsing binary (SB1) system, as determined from 20 sets of RV observations from PARAS and 6 supporting sets of observations from SOPHIE data, is found to be 16.953 d as against the 4.24 d period reported from SuperWasp photometry. It is likely that inadequate phase coverage of the transit with SuperWasp photometry led to the incorrect determination of the period for this system. We derive the spectral properties of the primary star from the observed stellar spectra: Teff = 5125 +/- 67 K, [Fe/H] = 0.1 +/- 0.14 and log g = 4.6 +/- 0.14, indicating a K1V primary. Applying the Torres relation to the derived stellar parameters, we estimate a primary mass 0.864 +/- 0.097 M_sun and a radius of 0.854 +/- 0.050 R_sun. We combine RV data with SuperWASP photometry to estimate the mass of the secondary, M_B = 0.098 +/- 0.007 M_sun, and its radius, R_B = 0.127 +/- 0.007 R_sun, with an accuracy of around 7%. Although the observed radius is found to be consistent with the Baraffe's theoretical models, the uncertainties on the mass and radius of the secondary reported here are model dependent and should be used with discretion. Here, we establish this system as a potential benchmark for the study of VLMS objects, worthy of both photometric follow-up and the investment of time on high-resolution spectrographs paired with large-aperture telescopes.

preprint2016arXiv

Evidence for Reflected Light from the Most Eccentric Exoplanet Known

Planets in highly eccentric orbits form a class of objects not seen within our Solar System. The most extreme case known amongst these objects is the planet orbiting HD~20782, with an orbital period of 597~days and an eccentricity of 0.96. Here we present new data and analysis for this system as part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS). We obtained CHIRON spectra to perform an independent estimation of the fundamental stellar parameters. New radial velocities from AAT and PARAS observations during periastron passage greatly improve our knowledge of the eccentric nature of the orbit. The combined analysis of our Keplerian orbital and Hipparcos astrometry show that the inclination of the planetary orbit is $> 1.22\degr$, ruling out stellar masses for the companion. Our long-term robotic photometry show that the star is extremely stable over long timescales. Photometric monitoring of the star during predicted transit and periastron times using MOST rule out a transit of the planet and reveal evidence of phase variations during periastron. These possible photometric phase variations may be caused by reflected light from the planet's atmosphere and the dramatic change in star--planet separation surrounding the periastron passage.

preprint2016arXiv

Molecular outflows in starburst nuclei

Recent observations have detected molecular outflows in a few nearby starburst nuclei. We discuss the physical processes at work in such an environment in order to outline a scenario that can explain the observed parameters of the phenomenon, such as the molecular mass, speed and size of the outflows. We show that outflows triggered by OB associations, with $N_{OB}\ge 10^5$ (corresponding to a star formation rate (SFR)$\ge 1$ M$_{\odot}$ yr$^{-1}$ in the nuclear region), in a stratified disk with mid-plane density $n_0\sim 200\hbox{--}1000$ cm$^{-3}$ and scale height $z_0\ge 200 (n_0/10^2 \, {\rm cm}^{-3})^{-3/5}$ pc, can form molecules in a cool dense and expanding shell. The associated molecular mass is $\ge 10^7$ M$_\odot$ at a distance of a few hundred pc, with a speed of several tens of km s$^{-1}$. We show that a SFR surface density of $10 \le Σ_{SFR} \le 50$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$ favours the production of molecular outflows, consistent with observed values.

preprint2016arXiv

Precision velocimetry planet hunting with PARAS: Current performance and lessons to inform future extreme precision radial velocity instruments

PARAS is a fiber-fed stabilized high-resolution cross-dispersed echelle spectrograph, located on the 1.2 m telescope in Mt. Abu India. Designed for exoplanet detection, PARAS is capable of single-shot spectral coverage of 3800 - 9600 A, and currently achieving radial velocity (RV) precisions approaching ~1 m/s over several months using simultaneous ThAr calibration. As such, it is one of the few dedicated stabilized fiber-fed spectrographs on small (1-2 m) telescopes that are able to fill an important niche in RV follow-up and stellar characterization. The success of ground-based RV surveys is motivating the push into extreme precisions, with goals of ~10 cm/s in the optical and <1 m/s in the near-infrared (NIR). Lessons from existing instruments like PARAS are invaluable in informing hardware design, providing pipeline prototypes, and guiding scientific surveys. Here we present our current precision estimates of PARAS based on observations of bright RV standard stars, and describe the evolution of the data reduction and RV analysis pipeline as instrument characterization progresses and we gather longer baselines of data. Secondly, we discuss how our experience with PARAS is a critical component in the development of future cutting edge instruments like (1) the Habitable Zone Planet Finder (HPF), a NIR spectrograph optimized to look for planets around M dwarfs, scheduled to be commissioned on the Hobby Eberly Telescope in 2017, and (2) the NEID optical spectrograph, designed in response to the NN-EXPLORE call for an extreme precision Doppler spectrometer (EPDS) for the WIYN telescope. In anticipation of instruments like TESS and GAIA, the ground-based RV support system is being reinforced. We emphasize that instruments like PARAS will play an intrinsic role in providing both complementary follow-up and battlefront experience for these next generation of precision velocimeters.

preprint2016arXiv

Proxima Centauri as a Benchmark for Stellar Activity Indicators in the Near Infrared

A new generation of dedicated Doppler spectrographs will attempt to detect low-mass exoplanets around mid-late M stars at near infrared (NIR) wavelengths, where those stars are brightest and have the most Doppler information content. A central requirement for the success of these instruments is to properly measure the component of radial velocity (RV) variability contributed by stellar magnetic activity and to account for it in exoplanet models of RV data. The wavelength coverage for many of these new instruments will not include the Ca II H&K or H-alpha lines, the most frequently used absorption-line tracers of magnetic activity. Thus, it is necessary to define and characterize NIR activity indicators for mid-late M stars in order to provide simultaneous activity metrics for NIR RV data. We have used the high-cadence UVES observations of the M5.5 dwarf Proxima Centauri from Fuhrmeister et al. (2011) to compare the activity sensitivity of 8 NIR atomic lines to that of H-alpha. We find that equivalent width-type measurements of the NIR K I doublet and the Ca II NIR triplet are excellent proxies for the canonical optical tracers. The Ca II triplet will be acquired by most of the new and upcoming NIR Doppler spectrographs, offering a common, reliable indicator of activity.

preprint2015arXiv

'Modal-noise' in single-mode fibers: A cautionary note for high precision radial velocity instruments

Exploring the use of single-mode fibers (SMFs) in high precision Doppler spectrometers has become increasingly attractive since the advent of diffraction-limited adaptive optics systems on large-aperture telescopes. Spectrometers fed with these fibers can be made significantly smaller than typical 'seeing-limited' instruments, greatly reducing cost and overall complexity. Importantly, classical mode interference and speckle issues associated with multi-mode fibers, also known as 'modal noise', are mitigated when using SMFs, which also provide perfect radial and azimuthal image scrambling. However, these fibers do support multiple polarization modes, an issue that is generally ignored for larger-core fibers given the large number of propagation modes. Since diffraction gratings used in most high resolution astronomical instruments have dispersive properties that are sensitive to incident polarization changes, any birefringence variations in the fiber can cause variations in the efficiency profile, degrading illumination stability. Here we present a cautionary note outlining how the polarization properties of SMFs can affect the radial velocity measurement precision of high resolution spectrographs. This work is immediately relevant to the rapidly expanding field of diffraction-limited, extreme precision RV spectrographs that are currently being designed and built by a number of groups.

preprint2015arXiv

An efficient, compact, and versatile fiber double scrambler for high precision radial velocity instruments

We present the design and test results of a compact optical fiber double-scrambler for high-resolution Doppler radial velocity instruments. This device consists of a single optic: a high-index $n$$\sim$2 ball lens that exchanges the near and far fields between two fibers. When used in conjunction with octagonal fibers, this device yields very high scrambling gains and greatly desensitizes the fiber output from any input illumination variations, thereby stabilizing the instrument profile of the spectrograph and improving the Doppler measurement precision. The system is also highly insensitive to input pupil variations, isolating the spectrograph from telescope illumination variations and seeing changes. By selecting the appropriate glass and lens diameter the highest efficiency is achieved when the fibers are practically in contact with the lens surface, greatly simplifying the alignment process when compared to classical double-scrambler systems. This prototype double-scrambler has demonstrated significant performance gains over previous systems, achieving scrambling gains in excess of 10,000 with a throughput of $\sim$87% using uncoated Polymicro octagonal fibers. Adding a circular fiber to the fiber train further increases the scrambling gain to $>$20,000, limited by laboratory measurement error. While this fiber system is designed for the Habitable-zone Planet Finder spectrograph, it is more generally applicable to other instruments in the visible and near-infrared. Given the simplicity and low cost, this fiber scrambler could also easily be multiplexed for large multi-object instruments.

preprint2015arXiv

Narrow escape: how ionizing photons escape from disc galaxies

In this paper we calculate the escape fraction ($f_{\rm esc}$) of ionizing photons from starburst galaxies. Using 2-D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the center of the disk along different angles through the superbubble and the gas disk. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of $\sim 40 ^\circ$, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scale heights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed $\sim [1- \cos (1 \, {\rm radian})] = 0.5$ from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time- and angle-averaged escape fraction on the mid-plane disk gas density (in the range $n_0=0.15-50$ cm $^{-3}$) and the disk scale height (between $z_0=10-600$ pc). We find that the escape fraction is related to the disk parameters (the mid-plane disk density and scale height) roughly so that $f_{\rm esc}^αn_0^2 z_0^3$ (with $α\approx 2.2$) is a constant. For disks with a given WNM temperature, massive disks have lower escape fraction than low mass galaxies. For Milky Way ISM parameters, we find $f_{\rm esc}\sim 5\%$, and it increases to $\approx 10\%$ for a galaxy ten times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the universe.

preprint2015arXiv

Response to Comment on "Stellar activity masquerading as planets in the habitable zone of the M dwarf Gliese 581"

Anglada-Escudé and Tuomi question the statistical rigor of our analysis while ignoring the stellar activity aspects that we present. Although we agree that improvements in multiparametric radial velocity (RV) modeling are necessary for the detection of Earth-mass planets, the key physical points we raised were not challenged. We maintain that activity on Gliese 581 induces RV shifts that were interpreted as exoplanets.

preprint2015arXiv

Stellar activity mimics a habitable-zone planet around Kapteyn's star

Kapteyn's star is an old M subdwarf believed to be a member of the Galactic halo population of stars. A recent study has claimed the existence of two super-Earth planets around the star based on radial velocity (RV) observations. The innermost of these candidate planets--Kapteyn b (P = 48 days)--resides within the circumstellar habitable zone. Given recent progress in understanding the impact of stellar activity in detecting planetary signals, we have analyzed the observed HARPS data for signatures of stellar activity. We find that while Kapteyn's star is photometrically very stable, a suite of spectral activity indices reveals a large-amplitude rotation signal, and we determine the stellar rotation period to be 143 days. The spectral activity tracers are strongly correlated with the purported RV signal of "planet b," and the 48-day period is an integer fraction (1/3) of the stellar rotation period. We conclude that Kapteyn b is not a planet in the Habitable Zone, but an artifact of stellar activity.

preprint2015arXiv

The APOGEE Spectroscopic Survey of Kepler Planet Hosts: Feasibility, Efficiency, and First Results

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

preprint2014arXiv

Determination of mass and orbital parameters of a low-mass star HD 213597B

HD 213597 is an eclipsing binary system which was detected by the STEREO spacecraft and was speculated to host a low-mass stellar companion. We used high-resolution spectroscopy with the 10-m Hobby-Eberly Telescope and the 1.2-m telescope in Mt Abu for radial velocity (RV) measurements of this source. We performed aperture photometry for this star on the STEREO archival data and thereby confirm the transit signature. We also did follow-up ground-based photometry with a 10-inch telescope from Mt Abu. The spectroscopic RV semi-amplitude of the primary (33.39 km s^-1) indicates that the secondary is an M dwarf making the system a short period F+M eclipsing binary. These RVs along with the inclination derived from our combined photometric analysis (i=84.9 degree), enable us to estimate the mass of the secondary as M_B ~ 0.286 M_sun and radius as R_B ~ 0.344 R_sun using an estimated mass M_A ~ 1.3 M_sun and radius R_B ~ 1.97 R_sun of the primary. Our spectral analysis returned the following parameters: T_eff = 6625 +/-121 K, [Fe/H] = -0.095 +\-0.08 and log(g) = 3.72 +/-0.22 for the primary. When log(g) is constrained to a value of 3.96, we derive T_eff = 6753 +/-52 K and [Fe/H] = -0.025+\-0.05.

preprint2014arXiv

Earthshine on a Young Moon: Explaining the Lunar Farside Highlands

The lunar farside highlands problem refers to the curious and unexplained fact that the farside lunar crust is thicker, on average, than the nearside crust. Here we recognize the crucial influence of Earthshine, and propose that it naturally explains this hemispheric dichotomy. Since the accreting Moon rapidly achieved synchronous rotation, a surface and atmospheric thermal gradient was imposed by the proximity of the hot, post-Giant-Impact Earth. This gradient guided condensation of atmospheric and accreting material, preferentially depositing crust-forming refractories on the cooler farside, resulting in a primordial bulk chemical inhomogeneity that seeded the crustal asymmetry. Our model provides a causal solution to the lunar highlands problem: the thermal gradient created by Earthshine produced the chemical gradient responsible for the crust thickness dichotomy that defines the lunar highlands.

preprint2014arXiv

In a hot bubble: why does superbubble feedback work, but isolated supernovae do not?

Using idealized 1-D Eulerian hydrodynamic simulations, we contrast the behavior of isolated supernovae with the superbubbles driven by multiple, collocated supernovae. Continuous energy injection via successive supernovae going off within the hot/dilute bubble maintains a strong termination shock. This strong shock keeps the superbubble over-pressured and drives the outer shock well after it becomes radiative. Isolated supernovae, in contrast, with no further energy injection, become radiative quite early ($\lesssim 0.1$ Myr, 10s of pc), and stall at scales $\lesssim 100$ pc. We show that isolated supernovae lose almost all of their mechanical energy by a Myr, but superbubbles can retain up to $\sim 40\%$ of the input energy in form of mechanical energy over the lifetime of the star cluster (few 10s of Myr). These conclusions hold even in the presence of realistic magnetic fields and thermal conduction. We also compare various recipes for implementing supernova feedback in numerical simulations. For various feedback prescriptions we derive the spatial scale below which the energy needs to be deposited for it to couple to the interstellar medium (ISM). We show that a steady thermal wind within the superbubble appears only for a large number ($\gtrsim 10^4$) of supernovae. For smaller clusters we expect multiple internal shocks instead of a smooth, dense thermalized wind.

preprint2014arXiv

Stellar Activity Masquerading as Planets in the Habitable Zone of the M dwarf Gliese 581

The M dwarf Gliese 581 is believed to host four planets, including one (GJ 581d) near the habitable zone that could possibly support liquid water on its surface if it is a rocky planet. The detection of another habitable-zone planet--GJ 581g--is disputed, as its significance depends on the eccentricity assumed for d. Analyzing stellar activity using the H-alpha line, we measure a stellar rotation period of 130+/-2 days and a correlation for H-alpha modulation with radial velocity. Correcting for activity greatly diminishes the signal of GJ 581d (to 1.5 sigma), while significantly boosting the signals of the other known super-Earth planets. GJ 581d does not exist, but is an artifact of stellar activity which, when incompletely corrected, causes the false detection of planet g.

preprint2014arXiv

The Habitable-zone Planet Finder Calibration System

We present the design concept of the wavelength calibration system for the Habitable-zone Planet Finder instrument (HPF), a precision radial velocity (RV) spectrograph designed to detect terrestrial-mass planets around M-dwarfs. HPF is a stabilized, fiber-fed, R$\sim$50,000 spectrograph operating in the near-infrared (NIR) z/Y/J bands from 0.84 to 1.3 microns. For HPF to achieve 1 m s$^{-1}$ or better measurement precision, a unique calibration system, stable to several times better precision, will be needed to accurately remove instrumental effects at an unprecedented level in the NIR. The primary wavelength calibration source is a laser frequency comb (LFC), currently in development at NIST Boulder, discussed separately in these proceedings. The LFC will be supplemented by a stabilized single-mode fiber Fabry-Perot interferometer reference source and Uranium-Neon lamp. The HPF calibration system will combine several other new technologies developed by the Penn State Optical-Infrared instrumentation group to improve RV measurement precision including a dynamic optical coupling system that significantly reduces modal noise effects. Each component has been thoroughly tested in the laboratory and has demonstrated significant performance gains over previous NIR calibration systems.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2013arXiv

A Cautionary Tale: MARVELS Brown Dwarf Candidate Reveals Itself To Be A Very Long Period, Highly Eccentric Spectroscopic Stellar Binary

We report the discovery of a highly eccentric, double-lined spectroscopic binary star system (TYC 3010-1494-1), comprising two solar-type stars that we had initially identified as a single star with a brown dwarf companion. At the moderate resolving power of the MARVELS spectrograph and the spectrographs used for subsequent radial-velocity (RV) measurements (R ~ <30,000), this particular stellar binary mimics a single-lined binary with an RV signal that would be induced by a brown dwarf companion (Msin(i)~50 M_Jup) to a solar-type primary. At least three properties of this system allow it to masquerade as a single star with a very low-mass companion: its large eccentricity (e~0.8), its relatively long period (P~238 days), and the approximately perpendicular orientation of the semi-major axis with respect to the line of sight (omega~189 degrees). As a result of these properties, for ~95% of the orbit the two sets of stellar spectral lines are completely blended, and the RV measurements based on centroiding on the apparently single-lined spectrum is very well fit by an orbit solution indicative of a brown dwarf companion on a more circular orbit (e~0.3). Only during the ~5% of the orbit near periastron passage does the true, double-lined nature and large RV amplitude of ~15 km/s reveal itself. The discovery of this binary system is an important lesson for RV surveys searching for substellar companions; at a given resolution and observing cadence, a survey will be susceptible to these kinds of astrophysical false positives for a range of orbital parameters. Finally, for surveys like MARVELS that lack the resolution for a useful line bisector analysis, it is imperative to monitor the peak of the cross-correlation function for suspicious changes in width or shape, so that such false positives can be flagged during the candidate vetting process.

preprint2013arXiv

Bar 1-Visibility Drawings of 1-Planar Graphs

A bar 1-visibility drawing of a graph $G$ is a drawing of $G$ where each vertex is drawn as a horizontal line segment called a bar, each edge is drawn as a vertical line segment where the vertical line segment representing an edge must connect the horizontal line segments representing the end vertices and a vertical line segment corresponding to an edge intersects at most one bar which is not an end point of the edge. A graph $G$ is bar 1-visible if $G$ has a bar 1-visibility drawing. A graph $G$ is 1-planar if $G$ has a drawing in a 2-dimensional plane such that an edge crosses at most one other edge. In this paper we give linear-time algorithms to find bar 1-visibility drawings of diagonal grid graphs and maximal outer 1-planar graphs. We also show that recursive quadrangle 1-planar graphs and pseudo double wheel 1-planar graphs are bar 1-visible graphs.

preprint2013arXiv

MARVELS-1: A face-on double-lined binary star masquerading as a resonant planetary system; and consideration of rare false positives in radial velocity planet searches

We have analyzed new and previously published radial velocity observations of MARVELS-1, known to have an ostensibly substellar companion in a ~6- day orbit. We find significant (~100 m/s) residuals to the best-fit model for the companion, and these residuals are naively consistent with an interior giant planet with a P = 1.965d in a nearly perfect 3:1 period commensuribility (|Pb/Pc - 3| < 10^{-4}). We have performed several tests for the reality of such a companion, including a dynamical analysis, a search for photometric variability, and a hunt for contaminating stellar spectra. We find many reasons to be critical of a planetary interpretation, including the fact that most of the three-body dynamical solutions are unstable. We find no evidence for transits, and no evidence of stellar photometric variability. We have discovered two apparent companions to MARVELS-1 with adaptive optics imaging at Keck; both are M dwarfs, one is likely bound, and the other is likely a foreground object. We explore false-alarm scenarios inspired by various curiosities in the data. Ultimately, a line profile and bisector analysis lead us to conclude that the ~100 m/s residuals are an artifact of spectral contamination from a stellar companion contributing ~15-30% of the optical light in the system. We conclude that origin of this contamination is the previously detected radial velocity companion to MARVELS-1, which is not, as previously reported, a brown dwarf, but in fact a G dwarf in a face-on orbit.

preprint2013arXiv

Superbubble breakout and galactic winds from disk galaxies

We study the conditions for disk galaxies to produce superbubbles that can break out of the disk and produce a galactic wind. We argue that the threshold surface density of supernovae rate for seeding a wind depends on the ability of superbubble energetics to compensate for radiative cooling. We first adapt Kompaneets formalism for expanding bubbles in a stratified medium to the case of continuous energy injection and include the effects of radiative cooling in the shell. With the help of hydrodynamic simulations, we then study the evolution of superbubbles evolving in stratified disks with typical disk parameters. We identify two crucial energy injection rates that differ in their effects, the corresponding breakout ranging from being gentle to a vigorous one. (a) Superbubbles that break out of the disk with a Mach number of order 2-3 correspond to an energy injection rate of order 10^{-4} erg cm^{-2} s^{-1}, which is relevant for disk galaxies with synchrotron emitting gas in the extra-planar regions. (b) A larger energy injection threshold, of order 10^{-3} erg cm^{-2} s^{-1}, or equivalently, a star formation surface density of \sim 0.1 solar mass yr^{-1} kpc^{-2}, corresponds to superbubbles with a Mach number \sim 5-10. While the milder superbubbles can be produced by large OB associations, the latter kind requires super-starclusters. These derived conditions compare well with observations of disk galaxies with winds and the existence of multiphase halo gas. Furthermore, we find that contrary to the general belief that superbubbles fragment through Rayleigh-Taylor (RT) instability when they reach a vertical height of order the scale height, the superbubbles are first affected by thermal instability for typical disk parameters and that RT instability takes over when the shells reach a distance of approximately twice the scale height.

preprint2013arXiv

The PRL Stabilized High Resolution Echelle Fiber-fed Spectrograph: Instrument Description & First Radial Velocity Results

We present spectrograph design details and initial radial velocity results from the PRL optical fiber-fed high-resolution cross-dispersed echelle spectrograph (PARAS), which has recently been commissioned at the Mt Abu 1.2 m telescope, in India. Data obtained as part of the post-commissioning tests with PARAS show velocity precision better than 2m/s over a period of several months on bright RV standard stars. For observations of sigma-Dra we report 1.7m/s precision for a period of seven months and 2.1m/s for HD 9407 over a period of 2 months. PARAS is capable of a single-shot spectral coverage of 3800A - 9500A at a resolution of about 67,000. The RV results were obtained between 3800A and 6900A using simultaneous wavelength calibration with a Thorium-Argon (ThAr) hollow cathode lamp. The spectrograph is maintained under stable conditions of temperature with a precision of 0.01 - 0.02C (rms) at 25.55C, and enclosed in a vacuum vessel at pressure of 0.1 +/-0.03 mbar. The blaze peak efficiency of the spectrograph between 5000A and 6500A, including the detector, is 30%; and about 25% with the fiber transmission. The total efficiency, including spectrograph, fiber transmission, focal ratio degradation (FRD), and telescope (with 81% reflectivity) is about 7% in the same wavelength region on a clear night with good seeing conditions.

preprint2012arXiv

The SDSS-HET Survey of Kepler Eclipsing Binaries: Spectroscopic Dynamical Masses of the Kepler-16 Circumbinary Planet Hosts

We have used high-resolution spectroscopy to observe the Kepler-16 eclipsing binary as a double-lined system, and measure precise radial velocities for both stellar components. These velocities yield a dynamical mass-ratio of q=0.2994+-0.0031. When combined with the inclination, i=90.3401+0.0016-0.0019 deg, measured from the Kepler photometric data by Doyle et al. 2011, we derive dynamical masses for the Kepler-16 components of M_A=0.654+-0.017 M_sun and M_B=0.1959+-0.0031 M_sun, a precision of 2.5% and 1.5% respectively. Our results confirm at the ~2% level the mass-ratio derived by Doyle et al. with their photometric-dynamical model, q=0.2937+-0.0006. These are among the most precise spectroscopic dynamical masses ever measured for low-mass stars, and provide an important direct test of the results from the photometric-dynamical modeling technique.

preprint2010arXiv

First Light results from PARAS: The PRL Echelle Spectrograph

We present the first light commissioning results from the Physical Research Laboratory (PRL) optical fiber-fed high resolution cross-dispersed Echelle Spectrograph. It is capable of a single- shot spectral coverage of 3700A to 8600A at R ~ 63,000 and is under very stable conditions of temperature (0.04°C at 23°C). In the very near future pressure control will also be achieved by enclosing the entire spectrograph in a low-pressure vacuum chamber (~0.01mbar). It is attached to a 1.2m telescope using two 50micron core optical fibers (one for the star and another for simultaneous Th-Ar spectral calibration). The 1.2m telescope is located at Mt. Abu, India, and we are guaranteed about 80 to 100 nights a year for observations with the spectrograph. The instrument will be ultimately used for radial-velocity searches of exoplanets around 1000 dwarf stars, brighter than 10th magnitude, for the next 5 years with a precision of 3 to 5m/s using the simultaneous Th-Ar spectral lamp reference technique. The spectrograph has already achieved a stability of 3.7m/s in short-term time scale and in the near future we expect the stability to be at 1m/s once we install the spectrograph inside the vacuum chamber.

preprint2010arXiv

The Pathfinder Testbed: Exploring Techniques for Achieving Precision Radial Velocities in the Near-Infrared

The Penn State Pathfinder is a prototype warm fiber-fed Echelle spectrograph with a Hawaii-1 NIR detector that has already demonstrated 7-10 m/s radial velocity precision on integrated sunlight. The Pathfinder testbed was initially setup for the Gemini PRVS design study to enable a systematic exploration of the challenges of achieving high radial velocity precision in the near-infrared, as well as to test possible solutions to these calibration challenges. The current version of the Pathfinder has an R3 echelle grating, and delivers a resolution of R~50,000 in the Y, J or H bands of the spectrum. We will discuss the on sky-performance of the Pathfinder during an engineering test run at the Hobby Eberly Telescope as well the results of velocity observations of M dwarfs. We will also discuss the unique calibration techniques we have explored, like Uranium-Neon hollow cathode lamps, notch filter, and modal noise mitigation to enable high precision radial velocity observation in the NIR. The Pathfinder is a prototype testbed precursor of a cooled high-resolution NIR spectrograph capable of high radial velocity precision and of finding low mass planets around mid-late M dwarfs.