Stellar spectral classification is a fundamental tool of modern astronomy, providing insight into physical characteristics such as effective temperature, surface gravity, and metallicity. Accurate and fast spectral typing is an integral need for large all-sky spectroscopic surveys like the SDSS and LAMOST. Here, we present the next version of PyHammer, stellar spectral classification software that uses optical spectral templates and spectral line index measurements. PyHammer v2.0 extends the classification power to include carbon (C) stars, DA white dwarf (WD) stars, and also double-lined spectroscopic binaries (SB2). This release also includes a new empirical library of luminosity-normalized spectra that can be used to flux calibrate observed spectra, or to create synthetic SB2 spectra. We have generated physically reasonable SB2 combinations as templates, adding to PyHammer the ability to spectrally type SB2s. We test classification success rates on SB2 spectra, generated from the SDSS, across a wide range of spectral types and signal-to-noise ratios. Within the defined range of pairings described, more than $95\%$ of SB2s are correctly classified.

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

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

Main sequence, fully-convective M dwarfs in eclipsing binaries are observed to be larger than stellar evolutionary models predict by as much as $10-15\%$. A proposed explanation for this discrepancy involves effects from strong magnetic fields, induced by rapid-rotation via the dynamo process. Although, a handful of single, slowly-rotating M dwarfs with radius measurements from interferometry also appear to be larger than models predict, suggesting that rotation or binarity specifically may not be the sole cause of the discrepancy. We test whether single, rapidly rotating, fully convective stars are also larger than expected by measuring their $R \sin i$ distribution. We combine photometric rotation periods from the literature with rotational broadening ($v \sin i$) measurements reported in this work for a sample of 88 rapidly rotating M dwarf stars. Using a Bayesian framework, we find that stellar evolutionary models underestimate the radii by $10-15\% \substack{+3 \\ -2.5}$, but that at higher masses ($0.18<M<0.4 M_{Sun}$) the discrepancy is only about $6\%$ and comparable to results from interferometry and eclipsing binaries. At the lowest masses ($0.08<M<0.18 M_{Sun}$), we find the discrepancy between observations and theory is $13-18\%$, and we argue that the discrepancy is unlikely to be due to effects from age. Furthermore, we find no statistically significant radius discrepancy between our sample and the handful of M dwarfs with interferometric radii. We conclude that neither rotation nor binarity is responsible for the inflated radii of fully convective M dwarfs, and that all fully-convective M dwarfs are larger than models predict.

We present a library of empirical stellar spectra created using spectra from the Sloan Digital Sky Survey's Baryon Oscillation Spectroscopic Survey (BOSS). The templates cover spectral types O5 through L3, are binned by metallicity from -2.0 dex through +1.0 dex and are separated into main sequence (dwarf) stars and giant stars. With recently developed M dwarf metallicity indicators, we are able to extend the metallicity bins down through the spectral subtype M8, making this the first empirical library with this degree of temperature \emph{and} metallicity coverage. The wavelength coverage for the templates is from 3650 Angstroms through 10200 Angstroms at a resolution better than R~2000. Using the templates, we identify trends in color space with metallicity and surface gravity, which will be useful for analyzing large data sets from upcoming missions like LSST. Along with the templates, we are releasing a code for automatically (and/or visually) identifying the spectral type and metallicity of a star.

We present radiation transfer models of rotating young stellar objects (YSOs) with hotspots in their atmospheres, inner disk warps and other 3-D effects in the nearby circumstellar environment. Our models are based on the geometry expected from the magneto-accretion theory, where material moving inward in the disk flows along magnetic field lines to the star and creates stellar hotspots upon impact. Due to rotation of the star and magnetosphere, the disk is variably illuminated. We compare our model light curves to data from the Spitzer YSOVAR project (Morales-Calderon et al. 2014, Cody et al. 2014) to determine if these processes can explain the variability observed at optical and mid-infrared wavelengths in young stars. We focus on those variables exhibiting "dipper" behavior that may be periodic, quasi-periodic, or aperiodic. We find that the stellar hotspot size and temperature affects the optical and near-infrared light curves, while the shape and vertical extent of the inner disk warp affects the mid-IR light curve variations. Clumpy disk distributions with non-uniform fractal density structure produce more stochastic light curves. We conclude that the magneto-accretion theory is consistent with certain aspects of the multi-wavelength photometric variability exhibited by low-mass YSOs. More detailed modeling of individual sources can be used to better determine the stellar hotspot and inner disk geometries of particular sources.