Recent photometric studies have revealed that surface spots that produce flux variations are present on virtually all L and T dwarfs. Their likely magnetic or dusty nature has been a much-debated problem, the resolution to which has been hindered by paucity of diagnostic multi-wavelength observations. To test for a correlation between magnetic activity and photometric variability, we searched for Ha emission among eight L3-T2 ultra-cool dwarfs with extensive previous photometric monitoring, some of which are known to be variable at 3.6 mu m or 4.5 mu m. We detected Ha only in the non-variable T2 dwarf 2MASS J12545393-0122474. The remaining seven objects do not show Ha emission, even though six of them are known to vary photometrically. Combining our results with those for 86 other L and T dwarfs from the literature show that the detection rate of Ha emission is very high (94%) for spectral types between L0 and L3.5 and much smaller (20%) for spectral types. >= L4, while the detection rate of photometric variability is approximately constant (30%-55%) from L0 to T8 dwarfs. We conclude that chromospheric activity, as evidenced by H alpha emission, and large-amplitude photometric variability are not correlated. Consequently, dust clouds are the dominant driver of the observed variability of ultra-cool dwarfs at spectral types, at least as early as L0.

We report the first detailed chemical abundance analysis of the exoplanet-hosting M-dwarf stars Kepler-138 and Kepler-186 from the analysis of high-resolution (R similar to 22,500) H-band spectra from the SDSS-IV-APOGEE survey. Chemical abundances of 13 elements-C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe-are extracted from the APOGEE spectra of these early M-dwarfs via spectrum syntheses computed with an improved line list that takes into account H2O and FeH lines. This paper demonstrates that APOGEE spectra can be analyzed to determine detailed chemical compositions of M-dwarfs. Both exoplanet-hosting M-dwarfs display modest sub-solar metallicities: [Fe/H](Kepler-138) = -0.09 +/- 0.09 dex and [Fe/H](Kepler-186) = -0.08 +/- 0.10 dex. The measured metallicities resulting from this high-resolution analysis are found to be higher by similar to 0.1-0.2 dex than previous estimates from lower-resolution spectra. The C/O ratios obtained for the two planet-hosting stars are near-solar, with values of 0.55 +/- 0.10 for Kepler-138 and 0.52 +/- 0.12 for Kepler-186. Kepler-186 exhibits a marginally enhanced [Si/Fe] ratio.

Open clusters have been the focus of several exoplanet surveys, but only a few planets have so far been discovered. The Kepler spacecraft revealed an abundance of small planets around small cool stars, therefore, such cluster members are prime targets for exoplanet transit searches. Kepler's new mission, K2, is targeting several open clusters and star-forming regions around the ecliptic to search for transiting planets around their low-mass constituents. Here, we report the discovery of the first transiting planet in the intermediate-age (800 Myr) Beehive cluster (Praesepe). K2-95 is a faint (Kp = 15.5 mag) M3.0 +/- 0.5 dwarf from K2's Campaign 5 with an effective temperature of 3471 +/- 124 K, approximately solar metallicity and a radius of 0.402 +/- 0.050 R-circle dot. We detected a transiting planet with a radius of 3.47(-0.53)(+0.78)R(circle plus) and an orbital period of 10.134 days. We combined photometry, medium/high-resolution spectroscopy, adaptive optics/speckle imaging, and archival survey images to rule out any false-positive detection scenarios, validate the planet, and further characterize the system. The planet's radius is very unusual as M-dwarf field stars rarely have Neptune-sized transiting planets. The comparatively large radius of K2-95b is consistent with the other recently discovered cluster planets K2-25b (Hyades) and K2-33b (Upper Scorpius), indicating systematic differences in their evolutionary states or formation. These discoveries from K2 provide a snapshot of planet formation and evolution in cluster environments and thus make excellent laboratories to test differences between field-star and cluster planet populations.

Condensate clouds fundamentally impact the atmospheric structure and spectra of exoplanets and brown dwarfs, but the connections between surface gravity, cloud structure, dust in the upper atmosphere, and the red colors of some brown dwarfs remain poorly understood. Rotational modulations enable the study of different clouds in the same atmosphere, thereby providing a method to isolate the effects of clouds. Here, we present the discovery of high peak-to-peak amplitude (8%) rotational modulations in a low-gravity, extremely red (J-K-s = 2.55) L6 dwarf WISEP J004701.06+680352.1 (W0047). Using the Hubble Space Telescope (HST) time-resolved grism spectroscopy, we find a best-fit rotational period (13.20 +/- 0.14 hr) with a larger amplitude at 1.1 mu m than at 1.7 mu m. This is the third-largest near-infrared variability amplitude measured in a brown dwarf, demonstrating that large-amplitude variations are not limited to the L/T transition but are present in some extremely red L-type dwarfs. We report a tentative trend between the wavelength dependence of relative amplitude, possibly proxy for small dust grains lofted in the upper atmosphere, and the likelihood of large-amplitude variability. By assuming forsterite as a haze particle, we successfully explain the wavelength-dependent amplitude with submicron-sized haze particle sizes of around 0.4 mu m. W0047 links the earlier spectral and later spectral type brown dwarfs in which rotational modulations have been observed; the large amplitude variations in this object make this a benchmark brown dwarf for the study of cloud properties close to the L/T transition.

Some low-mass stars appear to have larger radii than predicted by standard 1D structure models; prior work has suggested that inefficient convective heat transport, due to rotation and/or magnetism, may ultimately be responsible. In this thesis, we explore this possibility using a combination of 1D stellar models, 2D and 3D simulations, and analytical theory. First, we examine this issue using 1D stellar models constructed using the Modules for Experiments in Stellar Astrophysics (MESA) code. We begin by considering standard models that do not explicitly include rotational/magnetic effects, with convective inhibition modelled by decreasing a depth-independent mixing length theory (MLT) parameter αMLT. We provide formulae linking changes in αMLT to changes in the interior specific entropy, and hence to the stellar radius. Next, we modify the MLT formulation in MESA to mimic explicitly the influence of rotation and magnetism, using formulations suggested by Stevenson (1979) and MacDonald and Mullan (2014) respectively. We find rapid rotation in these models has a negligible impact on stellar structure, primarily because a star’s adiabat, and hence its radius, is predominantly affected by layers near the surface; convection is rapid and largely uninfluenced by rotation there. Magnetic fields, if they influenced convective transport in the manner described by MacDonald and Mullan (2014), could lead to more noticeable radius inflation. Finally, we show that these non-standard effects on stellar structure can be fabricated using a depth-dependent αMLT: a non-magnetic, non-rotating model can be produced that is virtually indistinguishable from one that explicitly parameterises rotation and/or magnetism using the two formulations above. We provide formulae linking the radially-variable αMLT to these putative MLT reformulations. We make further comparisons between MLT and simulations of convection, to establish how heat transport and stellar structure are influenced by rotation and magnetism, by looking at the entropy content of 2D local and 3D global convective calculations. Using 2D “box in a star” simulations, created using the convection code Dedalus, we investigate changes in bulk properties of the specific entropy for increasingly stratified domains. We observe regions stable against convection near the bottom boundary, resulting in the specific entropy in the bulk of the domain exceeding the bottom boundary value: this could be a result of physical effects, such as increased amounts of viscous dissipation for more supercritical, highly stratified cases, but may also be influenced by the artificial boundary conditions imposed by these local simulations. We then turn to 3D global simulations, created using the convection code Rayleigh, and investigate these same properties as a function of rotation rate. We find the average of the shell-averaged specific entropy gradient in the middle third of the domain to scale with rotation rate in a similar fashion to the scaling law derived via MLT arguments in Barker et al. (2014), i.e., |⟨ds/dr⟩| ∝ Ω^4/5.

I investigate the occurrence of extreme mid-infrared (MIR) excesses, a tracer of large amounts of dust orbiting stars, in low-mass stellar systems. Extreme MIR excesses, defined as an excess IR luminosity greater than 1% of the stellar luminosity (L_IR/L∗ > 0.01), have previously only been observed around a small number of solar-mass (M⊙) stars. The origin of this excess has been hypothesized to be massive amounts of orbiting dust, created by collisions between terrestrial planets or large planetesimals. Until recently, there was a dearth of low-mass (M∗ < 0.6M⊙) stars exhibiting extreme MIR excesses, even though low-mass stars are ubiquitous (~70% of all stars), and known to host multiple terrestrial planets (~3 planets per star). I combine the spectroscopic sample of low-mass stars from the Sloan Digital Sky Survey (SDSS) Data Release 7 (70,841 stars) with MIR photometry from the Wide-field Infrared Survey Explorer (WISE), to locate stars exhibiting extreme MIR excesses. I find the occurrence frequency of low-mass field stars (stars with ages > 1 Gyr) exhibiting extreme MIR excesses is much larger than that for higher-mass field stars (0.41 ± 0.03% versus 0.00067 ± 0.00033%, respectively). In addition, I build a larger sample of low-mass stars based on stellar colors and proper motions using SDSS, WISE, and the Two-Micron All-Sky Survey (8,735,004 stars). I also build a galactic model to simulate stellar counts and kinematics to estimate the number of stars missing from my sample. I perform a larger, more complete study of low-mass stars exhibiting extreme MIR excesses, and find a lower occurrence frequency (0.020 ± 0.001%) than found in the spectroscopic sample but that is still orders of magnitude larger than that for higher-mass stars. I find a slight trend for redder stars (lower-mass stars) to exhibit a higher occurrence frequency of extreme MIR excesses, as well as a lower frequency with increased stellar age. Lastly, I use white dwarf and low-mass star binary systems to investigate if the frequency of planetary collisions (traced through extreme MIR excesses) are increased in these environments. I find that these binary systems are more likely to host collisional debris, and therefore exhibit increased excess MIR flux, over single stars. These samples probe important questions into the habitability of worlds discovered around low-mass stars.

This dissertation presents the results of a study to (1) determine the census of the nearby southern M dwarf stellar population via three types of distances and (2) determine the multiplicity rate of nearby M dwarfs using two different search methodologies. The first part of this work reports three types of distance calculations (photographic, photometric, and trigonometric) for 1748 southern M dwarfs. Distances were estimated for 500 red dwarfs using photographic plate BRI magnitudes from SuperCOSMOS, while estimates were made for 667 stars using CCD VRI magnitudes. Both BRI and VRI were combined with 2MASS infrared JHK magnitudes. Distances for an additional 581 red dwarfs were derived from trigonometric parallaxes, 124 of which were measured for the first time during this work. For the second portion of this thesis, an all-sky sample of 1122 M dwarfs, known via trigonometric parallaxes to lie within 25 pc of the Sun, was surveyed for stellar companions at separations of 2" to 600". I-band images using primarily the CTIO/SMARTS 0.9m and the Lowell 42in telescopes were obtained in order to search these systems for companions at separations of 2" to 180". A complementary reconnaissance of wider companions to 600" was also done via blinking SuperCOSMOS BRI images. We find a stellar multiplicity fraction of 27.4 $\pm$ 1.3% for M dwarfs. Using this new gauge of M dwarf multiplicity near the end of the stellar main sequence, we calculate a multiplicity fraction of 30.1% for stellar systems of all types, implying that most systems are single. We find a peak in the separation distribution of the companions at 26 AU, i.e., distances on the scale of our Solar System, with a weak trend of smaller projected separations for lower mass primaries. A hint that M dwarf multiplicity may be a function of age/composition was revealed, with faster moving (and generally older) systems being multiple slightly less often. We calculate that at least 16% of M dwarf mass is made up of the stellar companions of multiple systems. Finally, we show that the mass function for M dwarfs increases to the end of the main sequence.

The response of low energy nuclear recoil in NaI(Tl) is investigated in the following experiment. Such detectors have been used recently to search for evidence of dark matter in the form of weakly interacting massive particles (WIMPs). NaI(Tl) crystal response to nuclear recoil energy deposition is a key element in these searches. I discuss the cosmological and experimental motivations for these experiments, followed by an overview of the physics of direct detection and current relevant WIMP search experiments. With the experiment motivations covered, the details of NaI(Tl) detectors are reviewed. The specifics of our experiment are laid out including the neutron production, neutron beam calibration, shielding optimization, experimental design and setup. Then the crystal response calibration studies and Geant4 simulations are discussed followed by the final quenching factor values and uncertainties. This experiment measured quenching factors for sodium recoils in the energy range of (9 keV-40 keV) of 19%-27% QF. These results are similar to current published measurements. Interesting features of the QF measurements include an increase at low energies and a dip in the values at 30 keV, the iodine K-shell absorption edge. The goal of this experiment was to add valuable measurements of nuclear recoils at low energies that are relevant to low-mass WIMP experiments. Future plans will improve and expand on these measurements in order to better understand the response of NaI(Tl) at low energies.

The aim of this dissertation is to address the topic of distinguishing very low mass stars from brown dwarfs through observational means. To that end, we seek to better characterize both populations and establish mechanisms that facilitate establishing an individual object's membership in either the very low mass star or the brown dwarf populations. The dissertation is composed of three separate observational studies. In the first study we report on our analysis of HST/NICMOS snapshot high resolution images of 255 stars in 201 systems within ~10 parsecs of the Sun. We establish magnitude and separation limits for which companions can be ruled out for each star in the sample, and then perform a comprehensive sensitivity and completeness analysis for the subsample of 138 M dwarfs in 126 systems. We calculate a multiplicity fraction of $0.0-0.0+3.5% for L companions to M dwarfs in the separation range of 5 to 70 AU, and $2.3-0.7+5.0% for L and T companions to M dwarfs in the separation range of 10 to 70 AU. Considering these results and results from several other studies, we argue that the so-called "brown dwarf desert" extends to binary systems with low mass primaries and is largely independent of primary mass, mass ratio, and separation. In the second study we construct a Hertzsprung-Russell diagram for the stellar/substellar boundary based on a sample of 63 objects ranging in spectral type from M6V to L4. We report new VRI photometry for 63 objects and new trigonometric parallaxes for 37 objects. We employ a novel SED fitting algorithm to determine effective temperatures, bolometric luminosities, and radii. We find evidence for the local minimum in the radius-temperature and radius-luminosity trends that may indicate the end of the stellar main sequence and the start of the brown dwarf sequence at $Teff ~2075K, log(L/Lsun) ~ -3.9, and (R/Rsun) ~ 0.086. The third study is a pilot study for future work and part of a long term search for astrometric binaries that have the potential to yield dynamical masses. We report the discovery of five new multiple systems and discuss their potential for determining dynamical masses: LHS 2071AB, GJ 1215 ABC, LTT 7434 AB, LHS 501 AC, and LHS 3738 AB.

abstract: The lowest-mass stars, known as M-dwarfs, form target samples for upcoming exoplanet searches, and together with lower-mass substellar objects known as brown dwarfs, are among prime targets for detailed study with high-contrast adaptive optics (AO) imaging and sub-millimeter interferometry. In this thesis, I describe results from three studies investigating the companion properties and environments of low-mass systems: (1) The 245-star M-dwarfs in Multiples (MinMs) Survey, a volume-limited survey of field M-dwarf companions within 15 pc, (2) the Taurus Boundary of Stellar/Substellar (TBOSS) Survey, an ongoing study of disk properties for low-mass members within the Taurus star-forming region, and (3) spectroscopy of a brown dwarf companion using the Gemini Planet Imager (GPI). Direct imaging of M-dwarfs is a sensitive technique to identify low-mass companions over a wide range of orbital separation, and the high proper motion of nearby M-dwarfs eases confirmation of new multiple stars. Combining AO and wide-field imaging, the MinMs Survey provides new measurements of the companion star fraction (CSF), separation distribution, and mass ratio distribution for the nearest K7-M6 dwarfs. These results demonstrate the closer orbital separations (~6 AU) and lower frequency (~23% CSF) of M-dwarf binaries relative to higher-mass stars. From the TBOSS project, I report 885µm Atacama Large Millimeter/sub-millimeter Array continuum measurements for 24 Taurus members spanning the stellar/substellar boundary (M4-M7.75). Observations of submillimeter emission from dust grains around the lowest-mass hosts show decreasing disk dust mass for decreasing host star mass, consistent with low frequencies of giant planets around M-dwarfs. Compared to the older stellar association of Upper Scorpius, Taurus disks have a factor of four higher mass in submillimeter-sized grains. From the GPI Exoplanet Survey, I describe near-infrared spectroscopy of an unusually red companion orbiting inside the debris disk of an F5V star. As the second brown dwarf discovered within the innermost region of a debris disk, the properties of this system offer important dynamical constraints for companion-disk interaction and a useful benchmark for brown dwarf and giant planet atmospheric study. / Dissertation/Thesis / Doctoral Dissertation Astrophysics 2017