The third transit of snow-line exoplanet Kepler-421b

The Kepler Mission has uncovered a handful of long-period transiting exoplanets that orbit from the cold outer reaches of their systems, despite their low transit probabilities. The atmospheres of these cold gas giant exoplanets are amenable to transit transmission spectroscopy enabling tests of planetary formation and evolution theories. Of particular scientific interest is Kepler-421b, a Neptune-sized exoplanet with a 704-day orbital period residing near the snow-line. Since the Kepler Spacecraft only observed two transits of Kepler-421b, the transit ephemeris is relatively uncertain. We observed Kepler-421 during the anticipated third transit of Kepler-421b in order to constrain the existence and extent of transit timing variations (TTVs). Barring significant TTVs, our visible light, time-series observations from the 4.3-meter Discovery Channel Telescope (DCT) were designed to capture pre-transit baseline and the partial transit of Kepler-421b. We find strong evidence in favor of transit models with no TTVs, suggesting that Kepler-421b is either alone in its system or is only experiencing minor dynamic interactions with an unseen companion. With the combined Kepler and DCT observations, we calculate the timing of future transits and discuss the unique opportunity to characterize the atmosphere of this cold, long-period exoplanet via transmission spectroscopy.http://adsabs.harvard.edu/abs/2016DPS....4812208DPublished versio

Long-term, multiwavelength light curves of ultra-cool dwarfs: II. The evolving light curves of the T2. 5 SIMP 0136 & the uncorrelated light curves of the M9 TVLM 513

We present multiwavelength, multi-telescope, ground-based follow-up photometry of the white dwarf WD 1145+017, that has recently been suggested to be orbited by up to six or more, short-period, low- mass, disintegrating planetesimals. We detect 9 significant dips in flux of between 10% and 30% of the stellar flux from our ground-based photometry. We observe transits deeper than 10% on average every ∼3.6 hr in our photometry. This suggests that WD 1145+017 is indeed being orbited by multiple, short-period objects. Through fits to the multiple asymmetric transits that we observe, we confirm that the transit egress timescale is usually longer than the ingress timescale, and that the transit duration is longer than expected for a solid body at these short periods, all suggesting that these objects have cometary tails streaming behind them. The precise orbital periods of the planetesimals in this system are unclear from the transit-times, but at least one object, and likely more, have orbital periods of ∼4.5 hours. We are otherwise unable to confirm the specific periods that have been reported, bringing into question the long-term stability of these periods. Our high precision photometry also displays low amplitude variations suggesting that dusty material is consistently passing in front of the white dwarf, either from discarded material from these disintegrating planetesimals or from the detected dusty debris disk. For the significant transits we observe, we compare the transit depths in the V- and R-bands of our multiwavelength photometry, and find no significant difference; therefore, for likely compositions the radius of single-size particles in the cometary tails streaming behind the planetesimals in this system must be ∼0.15 μm or larger, or ∼0.06 μm or smaller, with 2σ confidence

Characterizing mid-type M dwarfs in the Kepler field with the Discovery Channel Telescope and WIYN

Planet occurrence rates increase with decreasing stellar mass (later spectral types); therefore, M dwarf systems are our most promising targets in the search for exoplanets. The identification and characterization of stars in the original Kepler field was accomplished using photometry alone, resulting in large uncertainties for late-type stars like M dwarfs. In order to more accurately compute the planet occurrence rate around mid- type M dwarfs, we need to better constrain their stellar radii and masses, properties which strongly correlate with other stellar parameters such as temperature and metallicity. These measurements need to be performed on a statistically significant population of stars including systems with and without planets. Therefore, we have begun to spectroscopically characterize the properties of the 559 probable mid-type M dwarfs in the Kepler field using red optical spectra obtained with the DeVeny Spectrograph on the Discovery Channel Telescope (DCT) and Hydra on the WIYN telescope in order to constrain the planet occurrence rate for such stars. We will be presenting initial results from our DCT and WIYN observations, including new temperature, radius, and mass estimates which we can use in occurrence rate calculations.http://adsabs.harvard.edu/abs/2017AAS...22912608HPublished versio

Characterization of mid-type M dwarfs in the Kepler field

The planet occurrence rate has been found to increase with decreasing stellar mass (later spectral types) in the original Kepler field, and one out of four M dwarfs are expected to host Earth-sized planets within their habitable zones. M dwarf systems are, therefore, our most promising targets in the search for exoplanets. Yet the identification and characterization of M dwarfs in the Kepler field was accomplished using photometry alone and unfortunately this method provides large uncertainties for late-type stars. Notably absent from planet occurrence calculations are single planet mid-type M dwarfs (~M2-M6). In order to make an accurate calculation of the planet occurrence rate around mid-type M dwarfs, we need to constrain stellar radii and masses which depend on other stellar parameters (e.g. temperature and metallicity). We have identified 559 probable mid-type M dwarfs using photometric color selection criteria and have started to gather spectra of these objects in order to better constrain stellar properties and refine planet occurrence rates for this population. Here we outline the methods we are using for stellar classification and characterization and present some results from our initial data.http://adsabs.harvard.edu/abs/2016AAS...22743012HPublished versio

High-resolution broadband spectroscopy using externally dispersed interferometry at the Hale telescope: part 2, photon noise theory

High-resolution broadband spectroscopy at near-infrared (NIR) wavelengths (950 to 2450 nm) has been performed using externally dispersed interferometry (EDI) at the Hale telescope at Mt. Palomar, with the TEDI interferometer mounted within the central hole of the 200-in. primary mirror in series with the comounted TripleSpec NIR echelle spectrograph. These are the first multidelay EDI demonstrations on starlight. We demonstrated very high (10×) resolution boost and dramatic (20× or more) robustness to point spread function wavelength drifts in the native spectrograph. Data analysis, results, and instrument noise are described in a companion paper (part 1). This part 2 describes theoretical photon limited and readout noise limited behaviors, using simulated spectra and instrument model with noise added at the detector. We show that a single interferometer delay can be used to reduce the high frequency noise at the original resolution (1× boost case), and that except for delays much smaller than the native response peak half width, the fringing and nonfringing noises act uncorrelated and add in quadrature. This is due to the frequency shifting of the noise due to the heterodyning effect. We find a sum rule for the noise variance for multiple delays. The multiple delay EDI using a Gaussian distribution of exposure times has noise-to-signal ratio for photon-limited noise similar to a classical spectrograph with reduced slitwidth and reduced flux, proportional to the square root of resolution boost achieved, but without the focal spot limitation and pixel spacing Nyquist limitations. At low boost (∼1×) EDI has ∼1.4× smaller noise than conventional, and at >10× boost, EDI has ∼1.4× larger noise than conventional. Readout noise is minimized by the use of three or four steps instead of 10 of TEDI. Net noise grows as step phases change from symmetrical arrangement with wavenumber across the band. For three (or four) steps, we calculate a multiplicative bandwidth of 1.8:1 (2.3:1), sufficient to handle the visible band (400 to 700 nm, 1.8:1) and most of TripleSpec (2.6:1)

M-Dwarf Fast Rotators and the Detection of Relatively Young Multiple M-Star Systems

We have searched the Kepler light curves of ~3900 M-star targets for evidence of periodicities that indicate, by means of the effects of starspots, rapid stellar rotation. Several analysis techniques, including Fourier transforms, inspection of folded light curves, 'sonograms', and phase tracking of individual modulation cycles, were applied in order to distinguish the periodicities due to rapid rotation from those due to stellar pulsations, eclipsing binaries, or transiting planets. We find 178 Kepler M-star targets with rotation periods, P_rot, of < 2 days, and 110 with P_rot < 1 day. Some 30 of the 178 systems exhibit two or more independent short periods within the same Kepler photometric aperture, while several have three or more short periods. Adaptive optics imaging and modeling of the Kepler pixel response function for a subset of our sample support the conclusion that the targets with multiple periods are highly likely to be relatively young physical binary, triple, and even quadruple M star systems. We explore in detail the one object with four incommensurate periods all less than 1.2 days, and show that two of the periods arise from one of a close pair of stars, while the other two arise from the second star, which itself is probably a visual binary. If most of these M-star systems with multiple periods turn out to be bound M stars, this could prove a valuable way of discovering young hierarchical M-star systems; the same approach may also be applicable to G and K stars. The ~5% occurrence rate of rapid rotation among the ~3900 M star targets is consistent with spin evolution models that include an initial contraction phase followed by magnetic braking, wherein a typical M star can spend several hundred Myr before spinning down to periods longer than 2 days.Comment: 17 pages, 12 figures, 2 tables; accepted for publication in The Astrophysical Journa

DCT astrometry of very low-mass stars

Very low-mass stars (VLMs) are the smallest, least luminous stars in our galaxy, but nonetheless form one of the dominant (baryonic) populations. Precise distances and kinematics of VLMs can provide constraints on the smallest extremes of star formation, as well as important boundary constraints on the star formation process in general. We present preliminary results from a program to measure parallaxes and proper motions for a nearby sample of VLMs using the Large Monolithic Imager at the 4.3m Discovery Channel Telescope. These measurements will complement Gaia observations and allow us to construct high quality luminosity and mass functions, which will help to distinguish between VLM formation scenarios.http://adsabs.harvard.edu/abs/2016csss.confE..36SPublished versio

The mass-radius-luminosity-rotation relationship for M dwarf stars

NASA's future Transiting Exoplanet Survey Satellite (TESS) mission is expected to discover hundreds of terrestrial exoplanets orbiting around M dwarf stars, which will be nearby and amenable to detailed characterization. To accurately measure radii and equilibrium temperatures of these exoplanets, we need to know the host star properties, specifically mass, radius and luminosity, to equal accuracy. However, relationships for M dwarf stellar properties are poorly constrained, which leaves us unprepared to characterize exoplanets to be discovered by the TESS mission. The best way to determine relationships for M dwarf stars is to study mutually eclipsing binaries because the photometric and spectroscopic data empirically determine the physical parameters of the stars. We are conducting an on-going survey to measure infrared eclipses and individual spectra of carefully selected M dwarf eclipsing binary targets. We are using Mimir, a near-infrared wide-field imager, on the 72-inch Perkins Telescope near Flagstaff, Arizona, to determine the J, H, and K band magnitudes of the individual stars, and we are using Keck HIRES to measure the radial velocities of each component. Combining the observations, we determine the masses, radii and the semi-major axes of each component to an accuracy of 1%. We are also using measured parallaxes to determine the individual components' absolute infrared magnitudes and bolometric luminosities. The ultimate goal is to combine the measurements to determine the mass-radius-luminosity-rotation relationship for M dwarf stars. The relationship is critical for choosing the best TESS M dwarf exoplanets for detailed characterization.http://adsabs.harvard.edu/abs/2016AAS...22714221HPublished versio

Proper motions and parallaxes of very low-mass stars using DCT astrometry

Very low-mass stars (VLMs) are the smallest, least luminous stars in our galaxy, but nonetheless form one of the dominant (baryonic) populations. Precise distances and kinematics of VLMs can provide constraints on the smallest extremes of star formation, as well as important boundary constraints on the star formation process in general. However, Gaia will only be ~70% complete at the faint magnitudes of these objects. We present preliminary results from a program to measure parallaxes and proper motions for a nearby sample of 85 VLMs using the Large Monolithic Imager at the 4.3m Discovery Channel Telescope. We present proper motions for the entire sample and preliminary parallaxes for a few sources. These measurements will complement Gaia observations and allow us to construct high quality luminosity and mass functions, which will help to distinguish between VLM formation scenarios.http://adsabs.harvard.edu/abs/2017AAS...22912605SPublished versio

General moments of the inverse real Wishart distribution and orthogonal Weingarten functions

Let $W$ be a random positive definite symmetric matrix distributed according to a real Wishart distribution and let $W^{-1}=(W^{ij})_{i,j}$ be its inverse matrix. We compute general moments $\mathbb{E} [W^{k_1 k_2} W^{k_3 k_4} ... W^{k_{2n-1}k_{2n}}]$ explicitly. To do so, we employ the orthogonal Weingarten function, which was recently introduced in the study for Haar-distributed orthogonal matrices. As applications, we give formulas for moments of traces of a Wishart matrix and its inverse.Comment: 29 pages. The last version differs from the published version, but it includes Appendi