<figure><img src="image_1.png" /></figure>

# Placing the Spotted T Tauri Star LkCa 4 on an HR Diagram

Michael A. Gully-Santiago¹, Gregory J. Herczeg¹, Ian Czekala², Garrett Somers³, Konstantin Grankin⁴, Kevin R. Covey⁵,
J. F. Donati⁶,⁷, Silvia H. P. Alencar⁸, Gaitee A. J. Hussain⁹,¹⁰, Benjamin J. Shappee¹¹,<fnref n="17" />, Gregory N. Mace¹², Jae-Joon Lee¹³,
T. W.-S. Holoien¹⁴,¹⁵, Jessy Jose¹, and Chun-Fan Liu¹,¹⁶

¹ Kavli Institute for Astronomy and Astrophysics, Peking University, Yi He Yuan Lu 5, Haidian Qu, 100871 Beijing, People's Republic of China

² KIPAC, Physics and Astronomy, Stanford University, Stanford, CA, USA

³ Department of Physics and Astronomy, Vanderbilt University, 6301 Stevenson Center, Nashville, TN, 37235, USA

⁴ Crimean Astrophysical Observatory, Nauchny, Crimea 298409, Russia

⁵ Department of Physics & Astronomy, Western Washington University, Bellingham WA 98225-9164, USA

⁶ CNRS, IRAP/UMR 5277, Toulouse, 14 avenue E. Belin, F-31400 France

⁷ Université de Toulouse, UPS-OMP, IRAP, 14 avenue E. Belin, Toulouse, F-31400 France

⁸ Departamento de Física—ICEx—UFMG, Av. Antônio Carlos, 6627, 30270-901 Belo Horizonte, MG, Brazil

⁹ European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei München, Germany

¹⁰ Institut de Recherche en Astrophysique et Planétologie, Université de Toulouse, UPS-OMP, F-31400 Toulouse, France

¹¹ Carnegie Observatories, 813 Santa Barbara Street, Pasadena, CA 91101, USA

¹² Department of Astronomy, The University of Texas at Austin, 2515 Speedway Street, Austin, TX, USA

¹³ Korea Astronomy and Space Science Institute, 776 Daedeok-daero, Yuseong-gu, Daejeon 305-348, Korea

¹⁴ Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA

¹⁵ Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University, 191 W. Woodruff Avenue, Columbus, OH 43210, USA

¹⁶ Institute of Astronomy and Astrophysics, Academia Sinica (ASIAA), P.O. Box 23-141, Taipei 10641, Taiwan

Received 2016 September 27; revised 2017 January 13; accepted 2017 January 17; published 2017 February 21

## Abstract

Ages and masses of young stars are often estimated by comparing their luminosities and effective temperatures to pre-main-sequence stellar evolution tracks, but magnetic fields and starspots complicate both the observations and evolution. To understand their influence, we study the heavily spotted weak-lined *T-Tauri* star LkCa 4 by searching for spectral signatures of radiation originating from the starspot or starspot groups. We introduce a new methodology for constraining both the starspot filling factor and the spot temperature by fitting two-temperature stellar atmosphere models constructed from PHOENIX synthetic spectra to a high-resolution near-IR IGRINS spectrum. Clearly discernable spectral features arise from both a hot photospheric component $T_{\text{hot}} \sim 4100$ K and a cool component $T_{\text{cool}} \sim 2700-3000$ K, which covers $\sim 80\%$ of the visible surface. This mix of hot and cool emission is supported by analyses of the spectral energy distribution, rotational modulation of colors and of TiO band strengths, and features in low-resolution optical/near-IR spectroscopy. Although the revised effective temperature and luminosity make LkCa 4 appear to be much younger and of much lower mass than previous estimates from unspotted stellar evolution models, appropriate estimates will require the production and adoption of spotted evolutionary models. Biases from starspots likely afflict most fully convective young stars and contribute to uncertainties in ages and age spreads of open clusters. In some spectral regions, starspots act as a featureless “veiling” continuum owing to high rotational broadening and heavy line blanketing in cool star spectra. Some evidence is also found for an anticorrelation between the velocities of the warm and cool components.

**Key words:** stars: fundamental parameters – stars: individual (LkCa 4) – stars: low-mass – stars: statistics

## 1. Introduction

The causes of age uncertainties and of large luminosity spreads in young clusters are controversial (see reviews by Preibisch 2012; Soderblom et al. 2014). Absolute ages are usually inferred from comparisons between observed stellar properties and pre-main-sequence evolutionary models, which are sensitive to stellar birthlines and the input physics that controls the contraction and internal heating. Within clusters, any real differences in ages between stars may be masked by observational uncertainties in temperature and luminosity measurements (e.g., Hartmann 2001; Slesnick et al. 2008).

Magnetic activity is a likely source of significant uncertainty in both the models and observations of young low-mass stars. Convection at these young ages generates strong magnetic fields, as measured by Zeeman broadening and polarimetry (e.g., Johns-Krull 2007; Donati & Landstreet 2009) and as seen

in starspots (e.g., Stauffer et al. 2003; Grankin et al. 2008). Evolutionary models are just now starting to implement prescriptions for how spots and magnetic fields affect the interior structure (MacDonald & Mullan 2013; Jackson & Jeffries 2014; Somers & Pinsonneault 2015a; Feiden 2016). Stellar evolution models that include these effects can make a coeval 10 Myr population exhibit apparent age spreads of 3–10 Myr. Mass estimates may also be biased to lower masses, as measured for some eclipsing binaries (e.g., Stassun et al. 2014; Rizzuto et al. 2016). Observationally, starspots may also be responsible for biases in stellar effective temperatures derived by different methods. For example, the effective temperatures for 3493 young stars measured using the APOGEE spectrograph (1.5–1.70 µm at R = 22,500 Wilson et al. 2010) are offset from and usually cooler than prior measurements made at optical wavelengths (Cottaar et al. 2014). Similarly, near-IR spectral types of young stars have sometimes been evaluated as later than optical spectral types (Bouvier & Appenzeller 1992; Vacca & Sandell 2011).

<footnote n="17">¹⁷ Hubble, Carnegie-Princeton Fellow.</footnote>
Rotational modulation from starspots is understood as the
cause of cyclic variations in weak-lined T Tauri stars (WTTSs)
and some classical T Tauri stars (e.g., Vrba et al. 1986;
Herbst et al. 1994). Spots have recently been observed with
exceptional photometric precision from monitoring surveys
targeting exoplanet transits (e.g., Harrison et al. 2011; Daven-
port et al. 2015). The single band light curves from such
surveys cannot separate the relative areas and temperatures of
the emitting regions, and therefore cannot yield estimates for
the filling factor of spots without making assumptions about the
relative temperature contrast. Contemporaneous or near-
contemporaneous panchromatic photometric monitoring
(Herbst et al. 1994; Petrov et al. 1994; Bouvier et al. 1995;
Grankin et al. 2007; Cody et al. 2014) can break the
degeneracy between starspot filling factor and starspot temp-
erature, but these measurements are still limited to detecting
large longitudinally asymmetric starspots. Transiting planets or
eclipsing binaries passing in front of starspots (Desert
et al. 2011) can break geometric degeneracies, providing
estimates for the characteristic sizes and lifetimes of starspots
or starspot groups, but the number of such transiting systems
and demands of high photometric precision limit wider
application. Zeeman Doppler Imaging (Donati et al. 2014,
ZDI) produces the latitudinal and longitudinal distribution of
starspots in reconstructed brightness maps based on models for
spectro-polarimetric line shapes but would be unable to identify
uniform distributions of small spots. All of these techniques
would tend to underestimate the fractional coverage of
starspots distributed isotropically on the surface, stars viewed
pole-on, or stars with longitudinally symmetric bands of
starspots.

Spectroscopic strategies have the power to measure both the
starspot areal coverage and temperature in any geometric
configuration. For WTTSs, SED modeling indicated the
presence of cool spots (Wolk & Walter 1996), while spot
temperatures have been estimated from rotational modulation
of optical colors (e.g., Grankin 1998; Venuti et al. 2015;
Koen 2016). The TiO bands in optical spectra have been used
to measure spots on WTTSs (Petrov et al. 1994), stars with
ages of 125 Myr in the Pleiades (Fang et al. 2016), and evolved
stars (e.g., Neff et al. 1995; O'Neal et al. 2001, 2004). Recent
spectroscopic detections of spots on young stars have focused
on TW Hya, an active accretor, and DQ Tau, a close binary in
which both components are accreting (Debes et al. 2013; Bary
& Petersen 2014).

The star LkCa 4, a WTTS member of the Taurus Molecular
Cloud (Herbig et al. 1986; Downes & Keyes 1988; Strom
et al. 1989a, 1989b) is an ideal exemplar for a spotted pre-MS
star because it does not have any mid-IR or millimeter excess
(e.g., Andrews & Williams 2005; Furlan et al. 2006; Buckle
et al. 2015) and is not actively accreting (e.g., Edwards
et al. 2006; Cauley et al. 2012). LkCa 4 has no evidence for a
close companion from either direct imaging searches (Karr
et al. 2010; Kraus et al. 2011; Daemgen et al. 2015)<fnref n="18" /> or from
radial velocity (RV) searches (Nguyen et al. 2012) and high-
precision ZDI (Donati et al. 2014). This single star demon-
strates a large amplitude of sinusoidal photometric variability
indicative of large spots (Grankin et al. 2008; Xiao et al. 2012).
The variability amplitude cannot arise from an eclipsing binary

since the large RV variations would have been seen in
spectroscopic monitoring. All of these observations indicate
that the spectrum of LkCa 4 should be devoid of complicating
factors like near-IR excess veiling, accretion excess, or close
binarity.

LkCa 4 has recently been examined with ZDI (Donati et al. 2014), revealing a complex distribution of cool and warm spots in brightness map reconstructions. Dark polar spots extend to about 30° from the pole, with about five appendages reaching down to about 60° from the pole. There is also evidence for a hot spot in the reconstructed ZDI map. These spot maps demonstrate that multiple temperature components contribute to the broad-band emission from the star, which may explain the large differences between the 4100 K effective temperature measurement obtained from high-resolution optical spectra (Donati et al. 2014) compared to 3600 K obtained from an analysis dominated by TiO bands (Herczeg & Hillenbrand 2014).

In this work, we measure the starspot properties of LkCa 4 with four complementary techniques.<fnref n="19" /> Section 2 describes the full spectroscopic and photometric data set and quantifies the optical variability of LkCa 4 over the last 31 years based on all available photometric monitoring; a collection of spectral observations at various phases of variability is introduced. Section 3 and Appendix A describe a spectral inference technique aimed at generalizing line depth ratio analysis, which is applied to a high-resolution *H* and *K* near-IR spectrum of LkCa 4 in Section 3.5. Section 4 shows results from SED fitting, polychromatic photometric monitoring, and optical TiO line-depth ratio fitting. All lines of evidence are ultimately combined to build a consistent picture of the spectral and temporal evolution of LkCa 4. Section 5 describes how the two-temperature fit affects the placement of LkCa 4 on the HR diagram and what can be understood about pre-main-sequence stellar evolution from this exemplar. Finally, Section 6 describes how our results should be interpreted, the implications for analyses of spots in high-resolution spectra, and methodological limitations.

**2. Observations**

2.1. IGRINS Spectroscopy

We acquired observations of LkCa 4 with the Immersion
Grating Infrared Spectrograph, IGRINS (Park et al. 2014) on
the Harlan J. Smith Telescope at McDonald Observatory on
2015 November 18 09<sup>h</sup> UTC. IGRINS is a high-resolution
near-infrared echelle spectrograph providing simultaneous
*R* ≃ 45,000 spectra over 1.48–2.48 μm. The spectrograph has
two arms with 28 orders in *H* band and 25 orders in *K* band.
The data were reduced with the Pipeline Package (Lee 2015).
Telluric correction was performed by dividing the spectrum of
LkCa 4 by a spectrum of HR 1237, an A0V star observed
immediately before LkCa 4, both at airmass 1.1. The broad
hydrogen lines in the A0V star produced excess residuals in the
spectrum of LkCa 4. No effort was made to mitigate these flux
excess residuals, which affect several spectral orders. These
orders were ultimately excluded from further analyses.

<footnote n="18">¹⁸ The status for wide companions is less clear; see Stauffer et al. (1991), Itoh et al. (2008), Kraus & Hillenbrand (2009), Kraus et al. (2011), and Herczeg & Hillenbrand (2014).</footnote>

<footnote n="19">¹⁹ Some data, Python code, Starfish configuration files, and other documents related to this paper are available on the project’s GitHub repository at https://github.com/BrownDwarf/welter.</footnote>
## 2.2. ESPaDOnS Spectroscopy

We used ESPaDOnS on CFHT to obtain 12 high-resolution optical spectra of LkCa 4 from 2014 January 8–21 as part of the MaTYSSE Large Program. These spectra cover 3900–10000 Å at R ~ 68,000 and were obtained in the spectropolarimetry mode. The ZDI obtained from these observations were analyzed by Donati et al. (2014). In this paper, we concentrate on the intensity spectrum. To calculate TiO indices, the relative flux calibration was obtained from an approximate instrument sensitivity calculated from observations of 72 Tau. The spectral shapes match the flux-calibrated low-resolution spectra described in Section 2.3.

## 2.3. Low-resolution Optical and Near-IR Spectra

A flux-calibrated low-resolution optical spectrum of LkCa 4 was obtained using Palomar/DBSP (DBSP, Oke & Gunn 1982) on 2008 December 30 at 06:30 UT. These spectra cover 3200–8700 Å at a spectral resolution of R ~ 500. The data reduction and a spectral analysis are described by Herczeg & Hillenbrand (2014) within the context of a larger survey.

Near-infrared spectra of LkCa 4 were obtained at 01:30 UT on 2008 December 30, a few hours earlier than the DBSP spectra, using the TripleSpec spectrograph on the Astrophysical Research Consortium (ARC) 3.5 meter telescope at Apache Point Observatory in sunspot, New Mexico. LkCa 4 was observed for a total integration time of 90 s using an ABBA dither sequence along the slit to remove sky emission by differencing sequential exposures. With the 1"/1 slit used for these observations, TripleSpec provides nearly contiguous coverage from 0.95 to 2.46 µm at a resolution of R ~ 3200 (Wilson et al. 2004). Spectra were reduced with a variant of the SpeXTool pipeline, originally developed by Cushing et al. (2004) and modified for use with TripleSpec data. Spectra were differenced, flattened, extracted, and wavelength calibrated prior to telluric correction and flux calibration; the latter two corrections were performed using the XTELLCOR IDL package (Vacca et al. 2003) and a spectrum of HD 25175, a nearby A0V star observed directly after LkCa 4 at a similar airmass.

## 2.4. Photometric Monitoring

Photometry in B, V, and/or R bands in 1216 distinct epochs was assembled from published and unpublished data sources spanning 31 years in 22 observing seasons. Published data targeting LkCa 4 is composed of 29 BVR visits from 1985 to 1986 (Vrba et al. 1993), 26 BVR epochs from 1990 to 1991 (Bouvier et al. 1993), 284 visits in V band (278 with B and 268 with R) from 1992 August to 2004 October (Grankin et al. 2008), and 10 BVR visits from 2013 (Donati et al. 2014). ASAS3 (Pojmański 2004) acquired 63 V-band measurements from 2002 to 2004. Integral-OMC obtained 138 V-band measurements from 2006 to 2008 (Alfonso-Garzón et al. 2012). The AAVSO archive (Kafka 2016) includes 385 V, 23 B, and 10 R measurements from 2013 to 2016. Unpublished data from the ASAS-SN survey (Shappee et al. 2014) were obtained from 2012 January to 2016 March in 186 visits. Finally, recent photometry from the ongoing Crimean Astrophysical Observatory (CrAO) ROTOR project (Grankin et al. 2008) is presented here for the first time. The CrAO data includes 43 BVR visits from 2015 August to 2016 March. Figure 1 displays all the V-band photometry over the

interval 1985–2016. HJD, BJD, and JD are simply referred to as JD, which is accurate to the precision of available data.

Figure 2 shows all available V-band data grouped by the 22 observing seasons and phase-folded by the period P = 3.375 days obtained from multiterm Lomb–Scargle periodograms (Ivezić et al. 2014). The general appearance of the phase-folded light curves does not change with perturbations to the period on the scale of 0.003 days. The estimated BVR magnitudes at the time of the spectral observations are determined from regularized multiterm fits (VanderPlas & Ivezić 2015, i.e., Fourier series truncated to the first ~4 components) shown as the solid blue lines in Figure 2. Table 1 lists the estimated BVR photometry and the observing epoch for data from IGRINS, ESPaDOnS, 2MASS (Skrutskie et al. 2006), DBSP, and TripleSpec.

# 3. Fits to High-resolution Spectra

Starspots possess a spectrum distinct from the ambient photosphere. Approaches like ZDI and line-bisector analysis (e.g., Prato et al. 2008; Donati et al. 2014) target the modulation of stellar photospheric lines as starspots enter and exit the observable stellar disk. In contrast, the technique described in this section measures the spectrum that emerges from both the hot and the cool surfaces of the photosphere.

Starspot line depth ratio analysis has traditionally been limited to isolated portions of the spectrum that possess absorption lines attributable only to cool photospheric components and spectral lines arising from the warmer ambient photosphere (e.g., Neff et al. 1995; O’Neal et al. 2001). The apparent veiling of the lines and their respective temperature dependences can be combined to solve for the cool photosphere and ambient photosphere temperatures and relative areal coverages. This line depth ratio method suffers from the need to identify portions of the spectrum that possess easily identifiable strong lines, and from the need to assemble large atlases of observed spotless spectral templates to which the spotted star spectrum can be compared. In this section, we introduce a generalization of the line depth ratio analysis, which employs pixel-by-pixel modeling of all echelle spectral orders. By using all the spectral data, this strategy has the power to constrain starspot properties at relatively low filling factors and to identify weak lines that originate from the cool photosphere.

We took two approaches to characterizing the photospheric temperature from a spotted star spectrum. First, using the spectral fitting approach introduced in Section 3.1, we separately fit the optical (ESPaDOnS) and near-IR (IGRINS) spectra with distinct, single photospheric component models; these fits are described in Sections 3.2 and 3.3, respectively, and then compared in Section 3.4. Second, we use a two-temperature mixture model to perform a fit to the near-IR spectra in Section 3.5. The near-IR spectral range is preferred over the optical, since the cool photosphere emits most of its flux in the near-IR, enhancing the likelihood of direct detection of emission from the cool component.

## 3.1. Methodology

Czekala et al. (2015, hereafter C15) developed a modular framework<fnref n="20" /> to infer stellar properties from high-resolution spectra. The C15 technique forward models observed spectra with synthetic spectra from pre-computed model grids. The

<footnote n="20">²⁰ The open source codebase and its full revision history is available at https://github.com/iancze/Starfish. The experimental fork discussed in this paper is at https://github.com/gully/Starfish.</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Overview of LkCa 4 V-band photometric monitoring from 1986 to 2016. The vertical lines denote the observing epochs of 2MASS, IGRINS, ESPaDOnS, DBSP, and TripleSpec. The near-contemporaneous DBSP and TripleSpec epochs lay on top of each other on this scale, as do the 12 ESPaDOnS epochs. The abscissa range is equal to the current lifespan of the first author of this paper.</figcaption></figure>

intra-grid-point spectra are “emulated” in a process similar but superior to interpolation since it seamlessly quantifies the uncertainty attributable to the coarsely sampled stellar intrinsic parameters (see the Appendix of C15). The forward model includes calibration parameters, line spread functions, and a Gaussian process noise model to account for correlations in the residual spectrum. We employed the PHOENIX grid of pre-computed synthetic stellar spectra, which are provided at high spectral resolution over a wide wavelength range at 100 K intervals in stellar photospheric temperature<fnref n="21" /> in our range of interest (Husser et al. 2013). The modular framework was altered to accommodate starspot measurements in two ways. First, the single photospheric component was updated to include two photospheric components. Second, the MCMC sampling strategy was altered to accommodate the additional free parameters added by the two-component model.

The stellar photosphere is characterized as two photospheric components with a cool temperature $T_{\text{cool}}$ and hot temperature $T_{\text{hot}}$, with scalar solid angular coverages $\Omega_{\text{cool}}$ and $\Omega_{\text{hot}}$, respectively.<fnref n="22" /> The choice of nomenclature depends on the physical origin of the emission region, which is not directly observable without ancillary information. For many applications, it could be reasonably assumed that the cool photospheric component corresponds to a sunspot or sunspot group, and the hot component corresponds to the ambient photosphere. However, a cool ambient photosphere possessing localized hot spots—like plages on the Sun—could yield the same or similar two-temperature composite spectrum. In principle, the physical origin could be assessed with detailed comparison of average surface magnetic-field strengths through Zeeman line broadening or other magnetic-field sensitive techniques. For now, we index the two photospheric components in the most general way possible, simply cool and hot.

The two photospheric components share the same average intrinsic and extrinsic stellar parameters $v \sin i$, $\log g$, [Fe/H], $v_z$. The composite mixture model for observed flux density is

$$S_{\text{mix}} = \Omega_{\text{hot}}B(T_{\text{hot}}) + \Omega_{\text{cool}}B(T_{\text{cool}}) \quad (1)$$

where $B(T)$ is the spectral radiance from model spectra. The filling factor of the starspots is

$$f_{\text{cool}} = \frac{\Omega_{\text{cool}}}{\Omega_{\text{hot}} + \Omega_{\text{cool}}} \quad (2)$$

The term $f_{\text{cool}}$ represents an instantaneous, observational fill factor seen on one projected hemisphere, not $f_{\text{spot}}$, the “ratio of the spotted surface to the total surface areas” (Somers & Pinsonneault 2015a). In the limit of homogeneously distributed starspots and assuming the cool photosphere arises from starspots, $f_{\text{cool}}$ tends to $f_{\text{spot}}$, but in general $f_{\text{spot}}$ is not a direct observable since some circumpolar regions of inclined stars will always face away from Earth; $f_{\text{cool}}$ can vary cyclically through rotational modulation, whereas $f_{\text{spot}}$ changes secularly through starspot evolution that is not yet understood.

The two-component mixture model includes two more free parameters than the standard C15 model, namely $T_{\text{cool}}$ and $f_{\text{cool}}$. The addition of these two parameters makes the MCMC sampling much more correlated than it was for single component models. This challenge motivated a second technical change to the spectral inference framework, switching from Gibbs sampling (see C15) to ensemble sampling using emcee (Foreman-Mackey et al. 2013). The practical effect of this switch is that all 14 stellar and nuisance parameters are fit simultaneously in a single spectral order, making stellar parameter estimates $\theta_o$ unique for each spectral order $o$, whereas the C15 strategy had the power to provide a single set of stellar parameters $\theta$ that was based on all $N_{\text{ord}}$ spectral orders.

Interpreting the two-temperature mixture model as arising from the same photosphere requires single stars, or at least double stars with extremely large (>100) luminosity ratios or wide separations. Unresolved single-lined or double-lined binaries could mimic a signal that would be misattributed to starspots. As pointed out earlier, LkCa 4 has no detected

<footnote n="21">²¹ Single component PHOENIX photospheres are indexed by the label “effective temperature,” or $T_{\text{eff}}$, which presupposes that the emergent spectrum is comprised of only a single photospheric component. In this work, we avoid labeling the photospheric components by $T_{\text{eff}}$, and instead, refer to them as simply “photospheric temperatures.”</footnote>

<footnote n="22">²² For flux-calibrated spectra, the total solid angle $\Omega$ can be constrained, but for typical echelle spectrographs, only relative values of $\Omega_{\text{cool}}$ and $\Omega_{\text{hot}}$ can be inferred.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Phase-folded light curves constructed assuming *P* = 3.375 days for all 22 observing seasons. The blue solid lines show a “multiterm” regularized periodic fit, that is, keeping the first *M*<sub>max</sub> = 4 Fourier components (VanderPlas & Ivezić 2015). The vertical lines show the epochs at which spectra or ancillary photometry were obtained, with the same line styles and colors as Figure 1. LkCa 4 shows secular changess in its light-curve morphology. The IGRINS spectrum was acquired near the median flux level of the season.</figcaption></figure>
<table><caption>**Table 1**
Estimated BVR Magnitudes</caption><thead><tr><th>JD-2450000</th><th>Instrument</th><th><i>B</i></th><th><i>V</i></th><th><i>R</i></th></tr></thead><tbody><tr><td>781.7106</td><td>2MASS</td><td>14.02</td><td>12.61</td><td>11.25</td></tr><tr><td>4830.5637</td><td>TripleSpec</td><td>...</td><td>13.02</td><td>...</td></tr><tr><td>4830.7719</td><td>DoubleSpec</td><td>...</td><td>12.94</td><td>...</td></tr><tr><td>6665.7204</td><td>ESPaDOnS</td><td>14.36</td><td>12.87</td><td>11.71</td></tr><tr><td>6666.8505</td><td>ESPaDOnS</td><td>14.15</td><td>12.69</td><td>11.65</td></tr><tr><td>6667.7727</td><td>ESPaDOnS</td><td>14.31</td><td>12.87</td><td>11.63</td></tr><tr><td>6668.8699</td><td>ESPaDOnS</td><td>14.44</td><td>12.94</td><td>11.75</td></tr><tr><td>6672.8995</td><td>ESPaDOnS</td><td>14.25</td><td>12.80</td><td>11.56</td></tr><tr><td>6673.8408</td><td>ESPaDOnS</td><td>14.12</td><td>12.67</td><td>11.62</td></tr><tr><td>6674.7746</td><td>ESPaDOnS</td><td>14.46</td><td>12.99</td><td>11.85</td></tr><tr><td>6675.7396</td><td>ESPaDOnS</td><td>14.40</td><td>12.90</td><td>11.73</td></tr><tr><td>6676.7954</td><td>ESPaDOnS</td><td>14.18</td><td>12.72</td><td>11.65</td></tr><tr><td>6677.8699</td><td>ESPaDOnS</td><td>14.29</td><td>12.86</td><td>11.61</td></tr><tr><td>6678.7419</td><td>ESPaDOnS</td><td>14.53</td><td>13.02</td><td>11.80</td></tr><tr><td>6678.8950</td><td>ESPaDOnS</td><td>14.47</td><td>12.97</td><td>11.75</td></tr><tr><td>7344.8610</td><td>IGRINS</td><td>14.24</td><td>12.84</td><td>11.38</td></tr></tbody></table>

companion from direct imaging, and its RV variations are consistent with starspot modulation, so any binary companion would have to have an exceptionally large luminosity (or mass) ratio, rendering it undetectable with our spectral inference methodology with the finite signal-to-noise available in the IGRINS data. The absence of a detectable disk around LkCa 4 means that the spectral modeling requires no further parameters like veiling or accretion.

Further details of the spectral inference methodology are described in Appendix A.

## 3.2. Single Temperature Fitting to the ESPaDOnS Spectrum

We performed spectral fitting on an ESPaDOnS spectrum acquired on 2014 January 11 (JD 2456668.9). The spectrum was split into subsets of $N = 26$ chunks, corresponding approximately to spectral order boundaries. Fitting was performed separately on each spectral chunk. The spectral emulator was trained on stellar parameters in the range of $\log g \in [3.5, 4.0]$, $[\text{Fe/H}] \in [-0.5, 0.5]$, and photospheric temperatures of $T \in [3500, 4200]$ K.

Standard MCMC convergence criteria were assessed. Typical spectral chunks show photospheric temperature point estimates in the range $4000 \pm 130$ K. Figure 3 displays the point estimates with 5th and 95th percentile error bars placed at the central wavelength of each spectral chunk. The spectral chunk from 8473–8707 Å is an outlier, showing an estimated photospheric temperature of ~3500 K.

## 3.3. Single Temperature Fitting to the IGRINS Spectrum

We selected a subset of 32 of the 54 available IGRINS spectral orders<fnref n="23" /> with low telluric absorption artifacts. We performed unique spectral fitting on each order, with free stellar parameters $\theta = (T_{\text{eff}}, \log g, [\text{Fe/H}])$. The same analysis procedure described for the ESPaDOnS spectra is used here, though in the K band the temperature search range is expanded to $T \in [3000, 4200]$ K.

We find a larger dispersion in the point estimates for the stellar parameters derived from the IGRINS data than those

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Top panel: the best-fit temperature measured from our full spectrum fitting to each of 58 spectral orders in the optical and infrared portions of the spectrum, assuming a single component photosphere. The best-fit temperature measured in the K-band is about 800 K lower than that derived in optical. The triangles in the H-band show upper limits with best fits near the edge of our grid. Bottom panel: the flux ratio between the cool and hot components increases to the near-IR, which explains the observed trend in the best-fit temperature.</figcaption></figure>

derived in the optical spectrum. The most conspicuous trend is in the derived photospheric temperature as a function of wavelength shown in Figure 3. The photospheric temperature peaks at values of ~4200 K in the short-wavelength end of the H band and saturates at <3500 K at the long-wavelength end of the H band. The K band shows even lower derived photospheric temperatures of ~3300, or 700 K cooler than estimated from optical. No single photospheric temperature can describe all the spectral lines present in the high-resolution optical and IR spectra. Sources with such discrepancies have been seen previously, for example, in Figures 4 and 5 in Bouvier & Appenzeller (1992). These discrepancies are circumstantial evidence for the detection of spectral features attributable to a second photospheric component.

## 3.4. Heightened Sensitivity to Starspot Spectral Lines in the Infrared

The heightened sensitivity to starspot spectral lines as a function of wavelength can be understood in the following way. Starspots are cooler than their surrounding photosphere and will, therefore, have a longer wavelength of peak emission. The average ratio of flux density between starspot and bulk photosphere will increase with wavelength until asymptotically approaching a fixed value on the Rayleigh–Jeans tail (Wolk & Walter 1996). The bottom panel of Figure 3 shows the flux density ratio for patches of 2800 K and 4100 K photospheres with equal areas (50% filling factor). The blackbody ratios predict smooth flux ratio increases with wavelength, while the ratio of PHOENIX models shows wavelength regions with heightened sensitivity to starspot fluxes, for example, in the

<footnote n="23">²³ The prohibitive computational cost limited us from running all the spectral orders.</footnote>
J band. The visible portion of the spectrum will be dominated by the warm patches. The infrared spectrum will have a lower overall spot contrast and will have more net starspot flux than optical emission.

Some care should be taken when directly comparing results between the ESPaDOnS and IGRINS spectra since they were not taken at the same time. The particular ESPaDOnS spectrum used in Section 3.2 was obtained when LkCa 4 had an estimated V-band brightness of 12.94 mag, compared to 12.84 mag at the time of the IGRINS spectrum. The fainter magnitude during the ESPaDOnS spectrum acquisition implies a greater coverage fraction of cool material than during the IGRINS spectrum acquisition. The photospheric temperature derived in the optical bands and short-wavelength end of H-band yield similar values of ~4100 K, suggesting that the bulk<fnref n="24" /> spectral features are broadly consistent with emission from a single temperature component. A single ESPaDOnS order surrounding the TiO bands shows an exceptionally low estimated $T_{\text{eff}}$. The long-wavelength portion of the H band and all of the K band are more sensitive to starspot spectral signatures than the shorter wavelength portions.

### 3.5. Two-temperature Fitting to IGRINS Spectra

Full spectrum fitting was performed for 48 of the 54 available IGRINS spectral orders, omitting only the orders with the most pathological telluric spectral artifacts. We applied the mixture model as described in Section 3.1 and Appendix A. Standard MCMC convergence tests were evaluated. The stellar parameter ranges were $\log g \in [3.5, 4.0]$ and $[\text{Fe/H}] \in [-0.5, 0.5]$; the PHOENIX model spectrum temperature range was $T_{\text{hot}}$, $T_{\text{cool}} \in [2700, 4500]$. H-band fits had normal distribution priors in place: solar metallicity to ±0.05 dex, $\log g = 3.8 \pm 0.1$, and $v \sin i = 29 \pm 5 \text{ km s}^{-1}$.

The fit quality is first assessed by examining the consistency of the point estimates of $v_z$ and $v \sin i$ across the spectral orders. The distribution of $v_z$ and $v \sin i$ exposed extremely poor fits in two orders ($o = 91$ and 94), with all other orders demonstrating $v_z = 12.4 \pm 2.6 \text{ km s}^{-1}$ and $v \sin i = 28.8 \pm 2.0 \text{ km s}^{-1}$.

Overplotting forward-modeled spectra with the observed IGRINS spectra yielded insights on why some spectral orders perform better than others in assessing stellar properties. Orders with extremely poor telluric correction residuals, large spectral line outliers, and uncorrected H line residuals from A0V standard division, are all excluded from our final stellar parameter compilation. Several orders in the K band were rejected due to poor spectral fits influenced by deep metal lines. Metal lines can be biased by Zeeman broadening, which could explain poorer fit quality with wavelength (Johns-Krull 2007; Deen 2013). Some orders were also discarded because their fits were mediocre or uninformative. The remaining orders are relatively devoid of spectral line outliers and include the most information rich portion of the spectrum.

The IGRINS spectrum demonstrates some features that are present only in the hot photosphere model, and some features that are present only in the cool photosphere model. Figure 4 shows a selection of six such spectral features for a range of

plausible fill factors. In some cases (e.g., the top two panels), featureless cool photosphere spectra veil isolated spectral lines predicted in the hot photosphere models. In other cases (e.g., the middle row and lower left panel), the hot photosphere model veils sequences of shallow spectral features predicted in the cool photosphere model. The cool spectral models predict shallow features because line blanketing from multiple indistinct molecular bands overlaps from rotational broadening; any feature of interest will be biased to non-zero veiling.

The combination of line blanketing and high $v \sin i$ has probably hampered efforts to identify isolated spectral features suitable for line-depth-ratio analysis in the spectra of rotationally broadened young stars. One such feature, the OH 1.563 µm line analyzed by O'Neal et al. (2001), shows a clear pattern of three lines in our data, with the central line exceeding the depth of the adjacent two lines. The PHOENIX models predict a non-negligible hot photosphere contribution to the middle line for our range of hot and cool temperatures. In comparison, the forward-modeling technique described in this work thrives in the presence of long sequences of indistinct yet predictably correlated spectral features. The level of veiling of hot photospheric lines is set by the cool photosphere filling factor and temperature, while the level of veiling of the cool lines is set by the filling factor and temperature of the hot photosphere lines. Similar forward-modeling strategies have successfully identified patterns of weak metal lines in the line-blanketed spectra of low metallicity stars (Kirby 2011; Kirby et al. 2015). Figures 14 and 15 in the Appendix show 42 of the H- and K-band spectra on a log scale with a single random composite model spectrum overplotted.

Figure 5 shows the distribution of $T_{\text{hot}}$ (blue shading), $T_{\text{cool}}$ (red shading), and $f_{\text{cool}}$ (yellow shading) for all of the spectral orders, with the rejected orders grayed out, and the reliable order subset shown in bold. The point estimates for each spectral order are listed in Table 3 in the Appendix.

Best-fit values for the stellar parameters are listed in Table 2. Remarkably, the filling factor of cool photosphere exceeds 50%, with a best-fit value of $f_{\text{cool}} = 80\% \pm 5\%$.

## 4. The Two-temperature Fit to the SED and Stellar Rotation

In the previous section, we established that the high-resolution optical and near-IR spectra of LkCa 4 may be explained by a two-temperature photosphere. In this section, we test the two-component fit by comparing the combined spectrum to observed spectral features and rotational modulation. We adopt the best-fit two-temperature model to the IGRINS spectrum, with components of 2750 K covering 80% of the visible stellar surface and 4100 K covering 20% of the stellar surface, obtained when LkCa 4 had an estimated brightness of $V = 12.84$ mag. The cool and hot temperatures and the filling factor are blindly adopted without any adjustments to attempt to improve fits to the SED or broadband spectra.

Figures 6 and 7 compare synthetic spectra from the two-component photosphere to the SED and flux-calibrated spectra of LkCa 4. The only free parameters in this comparison, the luminosity and extinction, are both scaled to match the spectrum (see Sections 5.1 and 5.2). The optical-IR SED is obtained by estimating photometry from Grankin et al. (2008) during the 2MASS epoch, with $V \sim 12.61$ mag (see Table 1). The SED also includes mid-IR photometry from Spitzer/IRAC

<footnote n="24">²⁴ The bulk appearance and disappearance of spectral lines is mostly controlled by the temperature of the emitting region of the local photosphere. In other words, most of the variance in a spectrum is attributable to temperature, assuming near-solar metallicity. Starspots imbue dearths of flux in the line profiles of optical spectra, but these are secondary to the mere presence of temperature-sensitive lines.</footnote>
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 4.** Examples of spectral features in the observed IGRINS spectrum (light gray line). The composite spectrum model (purple thin line) is consistent with the observed spectrum for a range of fill factors, with examples of the median fill factor (middle panel of triptych) and ±2σ fill factors demarcated on the spectral postage stamps. The upper right triptych shows a Zeeman-sensitive Mg I line that is modeled with no attention to the magnetic field, which may bias estimates of T<sub>hot</sub> and/or f<sub>cool</sub> for individual spectral orders.</figcaption></figure>

(Hartmann et al. 2005), without adjusting for epoch. In the spectroscopic comparison, some minor discrepancies between the optical and near-IR spectra may be introduced because the spot coverage changed in the ~5 hrs between observations (ΔV = 0.08 mag). The synthetic spectrum is obtained from the Phoenix models, as in Section 3, and is extended beyond the longest wavelength (5 µm) to calculate the bolometric luminosity.

The synthetic models match the full spectrum reasonably well. When scaled separately to the red-optical and near-IR spectra, the synthetic models match even better (bottom panels of Figure 7), while single temperature photospheres fail to reproduce large spectral features. The 3900 K component fits the TiO bands at <7400 Å but fails to reproduce molecular bands at longer wavelengths. The 3500 K component suffers from the opposite problem, yielding molecular bands at short
<figure><img src="image_11.png" /><figcaption>**Figure 5.** Marginal probability distributions mirrored through the vertical axis (Waskom et al. 2014, "Violin plot") for 48 IGRINS orders for $T_{\text{hot}}$ (blue shade in the top panel), $T_{\text{cool}}$ (red shade, top), and fill factor $f_{\text{cool}}$ (yellow shade, bottom). The stellar parameters are derived independently in each spectral order. Spectral orders show differing levels of constraint on the cool photosphere and hot photosphere properties, including some orders ($\sigma = 104, 102, 100, 88, 83$) that yield especially tight cool spot filling factors. The starspot temperature is consistent with values even lower than 2700 K, the lower limit of the temperature range used. Many of the K-band orders that show lower estimates for the hot photosphere are unreliable due to spectral line outliers and uncorrected telluric residuals (light shaded distributions).</figcaption></figure>

<table><caption>**Table 2**

Adopted Values of LkCa 4 from IGRINS Spectra</caption><thead><tr><th><i>T</i><sub>hot</sub><br>(K)</th><th><i>T</i><sub>cool</sub><br>(K)</th><th><i>f</i><sub>σ</sub><br>(%)</th></tr></thead><tbody><tr><td>4100 ± 100</td><td>2750<sup>+250</sup><sub>-50</sub></td><td>80 ± 5</td></tr></tbody></table>

**Note.** These are the values at the epoch of the IGRINS spectrum acquisition.

wavelengths that are too deep. The two-temperature photosphere with the parameters from the best fit calculated in Section 3 reproduces the TiO bands, the hump in the H band, and the jump in flux at the long-wavelength end of the J band. This latter feature is caused by opacity in H$_2$O bands and is only seen when the fraction of emission from the cool component is high (Figure 8).

These spectral features should be rotationally modulated as the filling factors of the visible hot and cool components change. In Section 4.1, we first use our fit to the IGRINS spectrum to convert the V-band brightness to the filling factor of the two components. We subsequently calculate the color changes expected based on

<figure><img src="image_12.png" /><figcaption>**Figure 6.** SED of LkCa 4 (red) compared with the two-temperature spectrum obtained from the best fit to the IGRINS spectra (black), scaled to $A_V = 0.3$ and $\log L/L_\odot = 0.20$ at the epoch of the 2MASS observation (see Section 5).</figcaption></figure>

this spot coverage and compare these results to historical data (Section 4.2), and confirm our basic approach by demonstrating that TiO band depths vary with rotation (Section 4.3).
<figure><img src="image_13.png" /><figcaption>**Figure 7.** Top: the low-resolution optical/near-IR spectrum of LkCa 4 obtained from Palomar/DBSP and APO/Triplespec on 2008 December 30 (black), compared to a synthetic spectrum of a two-temperature photosphere (4100 K with 20% fill factor and 2750 K with 80% fill factor, reddened by $A_V$ = 0.3 mag; purple lines), a 3900 K spectrum reddened by 1.3 mag (blue lines), and a 3500 K spectrum with no reddening. All synthetic spectra are fit to the *J*-band spectrum of LkCa 4. The observations in gray and the single temperature photospheres are scaled by factors of three and nine for visual purposes. Bottom: same as the top panel, but the synthetic spectra are scaled separately to the optical spectrum at 0.75 µm and to the near-IR spectrum at 1.5 µm. Warm photospheres accurately reproduce molecular bands at 0.7 µm but fail to fit the spectral features at longer wavelengths. Cooler photospheres predict molecular bands at <0.7 µm that are much deeper than observed. The two-temperature photosphere accurately fits spectral features in both optical and near-IR wavelengths.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 8.** Constraint on the size and temperature of the cool spot from the H₂O absorption band depth between *J* and *H* bands. The contours show the flux ratio at 1.30 and 1.34 µm from the Phoenix models. Comparisons to the APO/TripleSpec spectrum, shown in the inset and normalized at 1.4 µm, yield best-fit solutions (solid purple line) and the range of acceptable fits (shaded cyan region, with contours of flux ratio in dashed black lines) for an epoch when the cool spot fraction was ~84%. The range of acceptable fits is based on an estimate uncertainty of ~2% in the observed flux ratio, which is dominated by uncertainty in telluric correction of strong H₂O absorption.</figcaption></figure>

### 4.1. The V-band Light Curve and Spot Coverage

The V-band brightness reflects the instantaneous filling factor of the cool and hot components. In the two-temperature photosphere, the V-band emission is dominated by the hotter component and is, therefore, a good proxy for the visible surface area of the hot component. Figure 9 shows the winter 2015–2016 ASAS-SN light curve, converted into a cool spot filling factor. During this period, the brightest and faintest phase corresponds to cool component fill factors of 74% and 86%, respectively. Since 1992 the V-band photometry has been as bright as 12.3 mag, which corresponds to a cool component fill factor of 65%, and as faint as 13.2 mag, or a cool component fill factor of 87%. This drastic change in visible area of the hot spot (35%–13%) is needed to explain the full ~1 mag range in brightness, assuming no temperature change in either component.

### 4.2. Rotational Modulation of Colors

The relative contributions from the hot and cool components change as LkCa 4 rotates, leading to color changes. We model this rotational modulation using the main-sequence colors (compiled by Pecaut & Mamajek 2013), which should have less contribution from spots than pre-main-sequence colors.
<figure><img src="image_15.png" /><figcaption>**Figure 9.** V-band magnitude in 2014–2015, converted into fill factor for the cool component. The optical brightness depends mostly on the hot component. If we fix a 75% filling factor, as measured in the IGRINS spectrum, to $V = 12.83$ at the time of the observation, then the V-band amplitude corresponds to filling factors of 67%–83%. The factor of ~2 change in visible surface area of the hot component, from 33% to 17%, is required to produce the $\Delta V = 0.6$.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 10.** Observed optical colors of LkCa 4, mostly from Grankin et al. (2008), compared with predictions from the rotation of a two-temperature photosphere model. The model is constructed instead by converting V-band brightness to a cool spot filling factor and subsequently calculating colors from main-sequence colors and bolometric corrections of Pecaut & Mamajek (2013; purple lines).</figcaption></figure>

The Cousins *R* photometry from Pecaut & Mamajek (2013) is converted to Johnson *R* for comparison with the Grankin et al. (2008) photometry, following the color transformation prescribed by Landolt (1983) and applied by Grankin et al. Photometry from synthetic spectra produced by Phoenix and BT-Settl models cannot be used here because the predicted *B*–*V* colors are much bluer than observed. Offsets between model and observed colors are at least in part attributed to extinction and are discussed in Section 5.1.

Limb darkening effects also become important if starspots are distributed in a finite number of large patches, not many small patches. The variations in the mean projection angle cause confounding color trends in the photometric modulations.

The optical emission is dominated by the hotter component, while both components contribute equally to the infrared emission (see the inset in the top panel of Figure 7). Figure 10

demonstrates that the *B*–*V* color from Grankin et al. (2008) does not depend on *V*, which is consistent with expectations for a single hot component with no contributions from the cooler component. The standard deviation of 0.03 mag in *B*–*V* color (which includes a ~0.01 mag uncertainty in both *V* and *B*) indicates that the spot temperature is stable to <50 K over decades. The stability of this temperature implies that the hotter component may be the ambient temperature of the photosphere.

The correlation between *V* and *V*–*R* indicates that the star becomes redder when the cool spot has a higher filling factor. Our simple model predicts that the correlation should be much weaker than is measured. Most likely, our two-temperature model for the photosphere is overly simplistic, since the *V*–*R* color could be reproduced with some contribution from intermediate temperatures.

Rotational modulation is expected to lead to much smaller brightness changes at near-IR wavelengths, with amplitudes of 0.15–0.2 mag, than the 0.6–0.8 mag amplitudes seen at optical
<figure><img src="image_17.png" /><figcaption>**Figure 11.** Variability in TiO bands measured with ESPaDOnS spectra. The top panel compares spectra between when LkCa 4 was bright (blue) and faint (red). The V-band emission is estimated from fits to the ASAS-SN light curve obtained during the same period. The bottom left and middle panels show correlations between V-band magnitude and individual TiO band indices, while the bottom right panel shows a similar correlation with the average of the TiO-6200, CaH-6800, TiO-7140, TiO-7600, and TiO-8500 indices.</figcaption></figure>

wavelengths. The smaller amplitude of near-IR brightness is also seen in a few spotted stars in the optical/IR monitoring of NGC 2264 (Cody et al. 2014), though those stars have a much smaller V-band amplitude than LkCa 4.

### 4.3. Rotational Modulation of TiO Band Depths

The depths of TiO and other molecular bands are common diagnostics of spectral type and therefore temperature in optical spectra (e.g., Kirkpatrick et al. 1991). In a two-temperature photosphere cool enough for molecules to form, the depth of observed TiO bands will depend on the fractional coverage of the components. The TiO band depths vary with spot coverage using 12 epochs of CFHT/ESPaDOnS spectra obtained over 14 days (~4 rotation periods) by Donati et al. (2014).

The spectra obtained near the estimated maximum and minimum optical brightness are shown in the top panel of Figure 11. Small but significant changes are seen in the TiO bands. On the other hand, the blue spectrum does not change. The TiO band depths in these spectra are measured from spectral indices as described in Herczeg & Hillenbrand (2014). The TiO indices are also combined into a single TiO index by calculating the number of standard deviations each point is from the median value of the index and then averaging the standard deviation.

The TiO bands trace the cool component. When the hotter component has a larger filling fraction, the TiO bands are veiled by the hotter component and should, therefore, be shallower. As expected, the TiO indices all correlate strongly

with the optical brightness. The TiO-7140 band depths are accurately reproduced by varying the two-temperature fit with rotation. Other bands have small offsets (see the middle panel of Figure 11), likely caused by either contributions from intermediate temperatures or by small errors in the synthetic spectra. In the spectral type scheme derived by Herczeg & Hillenbrand (2014), the change in the TiO-7140 index corresponds to a range from M1.5 (3640 K) to M2.2 (3530 K).

## 5. Placing LkCa 4 on the HR Diagram

### 5.1. Extinction Estimates to LkCa 4

An accurate extinction estimate from the full SED is challenging because of a strong dependence on the distribution of temperatures in the photosphere, which in our work is simplified to a two-temperature fit. However, isolating specific colors constrains the possible range of extinctions. Based on the following analysis, an $A_V = 0.3$ mag is adopted throughout this paper, though some uncertainty is introduced by confusion in the choice of template colors.

The B-V color is stable and not affected by cooler components, and may therefore be compared to the expected colors for the hotter (4100 K) temperature. The average $B-V = 1.41$ mag color from Grankin et al. (2008) has a color excess of $E(B-V) = 0.09$ mag relative to the empirical color expected for a main-sequence star and $E(B-V) = 0.20$ mag for colors of pre-main-sequence stars, based on the temperature-color tables of Pecaut & Mamajek (2013). The main-sequence colors are probably most
appropriate in this analysis because they are freer of spots. The B-V color excess leads to an extinction of $A_V = 0.28 \pm 0.15$ mag relative to main-sequence colors, or $0.62 \pm 0.15$ mag relative to pre-main-sequence colors. The listed uncertainties are calculated from extinction estimates relative to colors of 4000 and 4200 K photospheres. The $V-R_J$ color excess leads to $A_V = 0.48$ mag, which should be less affected by gravity (based on $V-I_C$ colors in Pecaut et al.) but is likely to be overestimated because the cool component contributes to the $R_J$-band emission, as seen in Figure 10. In the near-IR, the estimated $V = 12.61$ mag at the time of the 2MASS observations leads to a spot-weighted template (main-sequence) color $J-K_s = 0.87$ mag, with a range of ~0.05 mag for acceptable solutions with other temperatures and covering fractions. The observed $J-K = 0.93$ mag leads to $A_V = 0.35$ mag, with a range of 0.29 mag to account for the range of other acceptable photospheric solutions. The template near-IR color is more robust than optical colors to differences between the main-sequence and pre-main-sequence temperature scales because the pre-main-sequence $J-K$ color is redder at 4100 K and bluer at 2750 K.

The near-IR extinction to LkCa 4, and perhaps to other TTSs, will be overestimated if emission is generated from a multi-temperature spectrum but template colors are based on only the hotter component. This mismatch likely explains why near-IR extinctions are often higher than optical extinctions for young stars (Herczeg & Hillenbrand 2014), IR colors that are often redder than expected from models (Tottle & Mohanty 2015), and general challenges in near-IR templates of T Tauri stars (e.g., Espaillat et al. 2010; Fischer et al. 2011). For LkCa 4, the excess $E(J-K)$ color leads to an extinction of 0.84 mag for a 4100 K photosphere and Pecaut & Mamajek (2013) pre-main-sequence colors. This extinction is consistent with the color and SED analysis of Furlan et al. (2006), where a K7 spectral type was adopted, and is higher than extinctions of $A_V = 0.69$ mag measured from optical photometry (Kenyon & Hartmann 1995) and 0.35 mag measured by comparing optical spectra to young spectral templates (Herczeg & Hillenbrand 2014).

Our results tentatively support the lower extinction obtained at optical wavelengths, though the template colors used to calculate extinction are uncertain and will need revision. Differences between main-sequence and pre-main-sequence temperature scales and colors in the Pecaut & Mamajek (2013) tables are introduced by a combination of spots, in addition to any color differences caused by gravity differences. Uncertainty in colors may be solved with self-consistent comparisons to young stars with similar spectral types if spot coverage is also similar. Discrepancies in colors, spectral type/effective temperatures, and extinctions for LkCa 4, and perhaps more broadly for TTSs, are dominated by ignoring spots and not by spectral type-effective temperature scales. Fundamental problems in molecular opacities and photospheric models are likely important secondary effects.

## 5.2. Spot-corrected HR Diagram Placement

The spectrum of LkCa 4 is well reproduced by a cool component of 2750 K covering 80% of the surface and a 4100 K component covering 20% of the surface. Scaling this composite spectrum of LkCa 4 to the observed 2MASS photometry (~74% spot covering fraction) yields $R = 2.3 R_\odot$

and $\log L/L_\odot = -0.20$, adopting $A_V = 0.3$ mag and a distance of $131 \pm 2$ pc<fnref n="25" /> from the range of parallaxes of nearby Taurus members measured by Torres et al. (2012). At the median 80% spot covering fraction from 2015, the luminosity would be $\log L/L_\odot = -0.26$.

The effective temperature $T_{\text{eff}}$ of these two components may be calculated as

$$T_{\text{eff}} = [T_{\text{hot}}^4 (1 - f_{\text{cool}}) + T_{\text{cool}}^4 f_{\text{cool}}]^{0.25} \quad (3)$$

For the median fill factor of 80%, the observed effective temperature would be 3180 K, with a range of 3070–3280 for 86% and 74% spot fill factors. If the cool spot is 3000 K, at the limit of the acceptable temperature range, then these values would increase by 150 K. The apparent effective temperature in any single epoch depends on the relative fractions of the hot and cool spot. When accounting for the visible contributions of both hot and cool components, this wide range of possible observed effective temperatures in any single epoch underscores the difficulty in placing young stars on HR diagrams.

Figure 12 and Table 4 show that this effective temperature is much cooler than previous measurements from optical spectra, which are clustered around 4000 K. The luminosity does not change significantly because the broadband SED yields similar luminosities. When the photosphere is assumed to be a single temperature, the effective temperature depends on the wavelength where the measurement is made. For LkCa 4, the K-band temperature is closest to the effective temperature of the entire photosphere. These discrepancies may also explain the large temperature, systematic offsets identified when comparing temperatures from H-band APOGEE spectra to those obtained from other (usually optical) methods (Cottaar et al. 2014). Two-component photospheres would also explain deriving different spectral types in the optical and infrared for the same source. On the other hand, many GAIA-ESO spectra cover only a narrow wavelength around $H\alpha$ (e.g., Jeffries et al. 2014; Frasca et al. 2015), which would yield an effective temperature of ~4000 K. Spectral types and effective temperatures measured from TiO depths (e.g., Herczeg & Hillenbrand 2014) are cooler than those estimated from blue spectra but still overestimate the effective temperature. Temperature uncertainties for young K and M dwarfs are larger than any formal uncertainties in individual measurements.

In principle, the improved characterization of the photospheric emission and radius of LkCa 4 should lead to more accurate estimates of mass and age from pre-main-sequence tracks. However, most standard pre-main-sequence stellar evolution models do not take into account the effects of starspots, and those that do (Somers & Pinsonneault 2015a) limit cool spot areal coverage fractions to less than 50%. Comparison of LkCa 4 to unspotted models is inherently flawed, but for the purposes of gauging the effect size of spots on derived parameters, our effective temperature and luminosity yield a mass of 0.15–0.3 $M_\odot$ and an age of ~0.5 Myr in the Baraffe et al. (2015) evolutionary models.<fnref n="26" /> This age is much younger and this mass is much lower than previous estimates

<footnote n="25">²⁵ This distance is consistent with the Gaia distance of 130 ± 4 pc measured to the nearby Taurus members V410 Tau, V773 Tau, and BP Tau (Gaia Collaboration et al. 2016).</footnote>

<footnote n="26">²⁶ The Baraffe et al. (2015) evolutionary tracks were calculated using BT-Settl photospheric spectra (Allard 2014), which may introduce minor uncertainties when compared to results from the PHOENIX models used here.</footnote>
<figure><img src="image_18.png" /></figure>

**Figure 12.** Locations for LkCa 4 on an HR diagram, compared with models of pre-main-sequence evolution calculated by Baraffe et al. (2015) with isochrones (black lines) and evolution models of a single mass (dashed blue lines) as marked. The measured effective temperature and luminosity from this paper, based on the two-component fit and a median V-band magnitude, is shown as the black asterisk. The yellow shaded region corresponds to the range of apparent effective temperatures that would be measured as the hot component rotates into and out of our view for a cool temperature of 2750 K, while the green shaded region shows the same range for a cool temperature of 3000 K. The blue circles correspond to the measurement at blue-optical wavelengths (the large circle and error bar is from Donati et al. 2014), the purple square corresponds to the measurement from low-resolution optical spectra, biased to TiO band strengths, by Herczeg & Hillenbrand (2014, biased to TiO bands), while the red diamond corresponds to what we would measure from the K-band spectrum and 2MASS J-band magnitude.

made assuming the warm component solely accounted for the effective temperature.

Even though LkCa 4 is only a single data point, the cool spot coverage fraction may be evidence that strong magnetic fields are inhibiting convection in LkCa 4. As a result, the surface is cooler and releases less radiation, slowing the contraction rate. LkCa 4 then appears more luminous and cooler than expected for a star of its genuine age from evolutionary tracks that do not consider magnetic fields. This shift to lower temperature may also help to resolve some of the age discrepancies between intermediate-mass and low-mass young stars (e.g., Herczeg & Hillenbrand 2015). The corrected placement of spotted stars on the HR diagram reflects the accurate luminosity, radius, and effective temperature, but evolutionary models are not yet equipped to interpret such heavily spotted stars (Somers & Pinsonneault 2015a). In an HR diagram, spotted evolutionary model tracks at a given mass and age shift to higher luminosity and lower effective temperature. Such a spot-corrected HR diagram position would alter the age and mass estimates to >0.2–0.3 $M_\odot$ and 0.5 Myr; in other words, back in the direction of previous naïve estimates. However, the full magnitude of the correction is an outstanding question.

The IGRINS and ESPaDOnS spectra indicate $v \sin i \sim 28$ km s$^{-1}$ and a rotation period of 3.375 days. Combining these numbers gives $R \sin i \sim 1.87 R_\odot$. Our HR diagram analysis gives $R \sim 2.3 R_\odot$, which suggests $\sin i \sim 0.813$, or an inclination of about 35° from edge on. These values show broad consistency between rotational properties, spectral fitting values, and our interpretation of a tilted star revealing a circumpolar region with large areal coverage of starspots.

As seen in Figures 1 and 2, the areal coverage of spots appears to wax and wane secularly through the observing seasons of the last 31 years. LkCa 4 is currently in a relatively faint phase, 0.2 mag fainter in V than it was through much of the late 1980s and early 1990s. It is conceivable that LkCa 4 happens to be going through a short-lived episode of relatively high coverage of cool photosphere analogous to—albeit distinct from—solar maximum. The intensity and duration of such episodes would further confound the interpretation of any instantaneous position on the HR diagram.

## 6. Discussion

In this paper, we have demonstrated that the IGRINS spectra of LkCa 4 are well fit with two temperatures, a hot photosphere of 4100 K and large cool patches of 2700–3000 K that cover 80% of the stellar surface. This solution reproduces the high-resolution HK-band spectra obtained with IGRINS, the broad-band SED, and rotational modulation of spectral features. Such a high covering fraction for a spot seems extraordinary, but may be common among young stars. In this section, we describe that the need for a large cool spot is robust, despite some simplifying limitations of our spectral inference technique. We also discuss how such a large spot could be interpreted, and suggest that even though LkCa 4 appears extreme, other young sources have similar characteristics.

### 6.1. An Absolute Lower Limit of the Starspot Filling Factor

This simple thought experiment demonstrates that the high starspot filling factor must be approximately correct. A lower
limit to the starspot filling factor can be set by the size of a
black spot that can induce such large ΔV. Consider the
following extremely conservative assumptions.

1. The starspot emits zero flux. This statement is equivalent to the starspot being at absolute zero.

2. The stellar disk shows no starspots at maximum brightness.

3. The starspot transits the center of the stellar disk.

4. The star shows no feature brighter than the quiet photosphere at the visible surface

For LkCa 4 in 2015, the ΔV ~ 0.5 mag amplitude corresponds to an absolute minimum fill factor of f > 0.37. At least 37% of the stellar disk is covered in starspots at minimum light. In the 2004 season, LkCa 4 exhibited ΔV ~ 0.8, corresponding to a minimum fill factor of f = 0.52. If any of assumptions 1–3 are relaxed, the minimum fill factor increases. If assumption 4 is relaxed, the minimum fill factor decreases. The large amplitude photometric variations in LkCa 4 are therefore evidence for either a hot spot in excess of the ambient photosphere or a large coverage fraction of cool starspots.

## 6.2. Interpretation of the Large Filling Factor of Cool Starspots

We have shown that LkCa 4 exposes between 74% and 86% cool photosphere toward our viewing direction in winter 2015–2016, the most recent season for which there is data. In previous seasons, like fall 2004 (see Figures 2 and 9), the filling factor of cool spots could have exceeded 90%. This exceptionally large fraction of surface area covered by spots defies the conventional wisdom that starspots are the minority constituent of the stellar surface. The evidence for this derived fill factor has been shown in the SED, high-resolution spectral inference, direct spectral detection, and TiO variability modeling. Perhaps the strongest evidence for the large fill factor of the cool photosphere is the large amplitude of V-band modulation, which cannot be explained with small spots, even if a spot temperature of absolute zero is assigned to the cool regions.

One conceivable interpretation of the 74%-86% coverage fraction of cool photosphere is a stellar axis inclined toward our observing direction, with a large circumpolar spot covering much of the observable hemisphere, as described by Donati et al. (2014). The ambient hot photosphere would then be close to the equator, which would be limb-darkened due to projection effects. Such limb darkening would cause a nonlinear transformation between our observed fill factor $f_{\text{cool}}$, and the surface area coverage fraction of spots. It is conceivable, and probably likely, that the value we report for $f_{\text{cool}}$ is an overestimate of the surface area actually covered by spots. In this circumpolar spot scenario, local fingers of cool spots could reach down to near the equator, entering into and out of the field of view and giving rise to the large amplitude photometric variation seen on LkCa 4. This scenario is consistent with the inclination derived from relating $v \sin i$, rotation period, and estimates for the stellar radius, as previously noted. Such estimates are limited by estimates of the radius since uncertainties in effective temperature and luminosity cause a wide range of radius estimates.

The starspot geometry is also consistent with morphologies observed in interferometric maps of $\zeta$ Andromedae: large areal

coverage of permanent polar spots, and smaller, transient networks of varying darkness, which thread across the whole surface (Roettenbacher et al. 2016).

## 6.3. Is LkCa 4 an Extreme Source?

The large amplitude of cyclical optical variability of LkCa 4 exceeds all other counterparts from the long-term monitoring of WTTSs by Grankin et al. (2008). From a statistical perspective, it might not be a surprise that the source with the largest photometric variability also possesses among the largest coverage fraction of cool spots. When compared with spot covering fractions of 15%-60% for evolved giants and subgiants (Chugainov 1976; Berdyugina 2005), the 80% spot coverage of LkCa 4 seems high. However, photometric and spectroscopic evidence demonstrates that many other young low-mass stars also have large spot covering fractions.

Several WTTS in the Taurus–Auriga star-forming region have large amplitudes in periodic light variations, achieving 0.4–0.8 mag in the V band: LkCa 4, V410 Tau (Herbst 1989), LkCa 7, V827 Tau, V830 Tau, and V819 Tau (Grankin et al. 2008). In a sample that is not biased to variable stars, WTTSs observed in NGC 2264 show photometric variations similar to—but of lower average amplitude than—LkCa 4 (Cody et al. 2014; Venuti et al. 2015). Modeling of polychromatic photometric monitoring of the locus of WTTSs indicates a lower limit to the true coverage fraction $f_{\text{eq}} \sim 10\%–30\%$ and a difference $T_{\text{hot}} – T_{\text{cool}} \sim 500$ K.

Large amplitudes of the light variation indicate the existence of very extended spotted regions on the stellar surface. To estimate the total area and mean temperature of the spots in the visible stellar hemisphere, Grankin used a simple non-parametric model for analysis of V410 Tau, V819 Tau, V827 Tau, and V830 Tau (Grankin 1998, 1999). For the stars with the highest light-curve amplitudes, this model shows that the spotted regions cover from 17% to 73% of the visible stellar hemisphere and that the mean temperature of the cool regions is 500–1400 K lower than the ambient photosphere. An analysis of the calculated spot parameters of V410 Tau revealed a correlation between the amplitude of periodicity and the spot distribution nonuniformity, which is defined as the difference of spot coverage in the visible hemispheres of the star at maximum and minimum light. As the amplitude increases from 0.39 to 0.63 mag, the degree of nonuniformity in the spot distribution increases from 21% to 37%. An alternative analysis of TiO bands indicate a spot coverage of ~41% for V410 Tau (Petrov et al. 1994). Besides the high amplitude of variability, these objects show the phenomenon of long-term stability of a brightness minimum in the interval from 5 to 19 years.

The Pleiades provides an older (~125 Myr, Stauffer et al. 1998) sample of young stars, many of which have starspots that cover large fractions of the stellar surface. For example, anomalous colors of K dwarfs in the Pleiades can be explained by ~50% fill factors of spots (Stauffer et al. 2003), while LAMOST spectra of the TiO band indicates that ~30% spot coverages may be common for 3500–5000 K stars (Fang et al. 2016). Recent photometric monitoring of the Pleiades indicates that the anomalous colors scale with rotation rate, strengthening the case for magnetic-field induced starspots as the cause of the anomaly (Covey et al. 2016). However, the Covey et al. (2016) Pleiades rotation rates do not scale with the ΔV photometric amplitude at the precision of the PTF
data, suggesting that while the overall starspot coverage
increases with rotation rate, the longitudinally asymmetric
component does not. LkCa 4 represents an example of such an
architecture, with more than 50% of the visible stellar surface
covered in spots at all observational phases. The long-
itudinally symmetric component evades photometric modula-
tion detection but contributes to otherwise anomalous color
offsets.

Spots are also likely important in estimating the properties of
generic low-mass young stars. As an example, the well-studied
CTTS TW Hya has a ~4100 K photosphere measured from
high-resolution optical spectra (e.g., Yang et al. 2005). How-
ever, the near-IR spectrum is consistent with ~M2, when
compared with main-sequence dwarf stars (Vacca & San-
dell 2011). Debes et al. (2013) reconciled this discrepancy by
invoking spots, though the spectral contribution from accretion
and disk continuum emission was uncertain. McClure et al.
(2013) found that the TW Hya spectrum was well reproduced
with near-IR spectra of stars with similar (~K7) optical spectral
type, as long as the comparison spectra were young. Herczeg &
Hillenbrand (2014) found an intermediate spectral type for TW
Hya and other TTSs by focusing on TiO bands at redder
wavelengths than previous spectra, and also by using young
stars as photospheric templates. Similar discrepancies have
been found between optical and near-IR spectral types (e.g.,
Bouvier & Appenzeller 1992; Bary & Petersen 2014; Cottaar
et al. 2014) and in comparisons of spectral types measurement
based on TiO bands and those based on high-resolution blue
spectra (see the discussion in Herczeg & Hillenbrand 2014).

The mismatches between near-IR spectra of dwarf stars and young stars (along with some of the color discrepancies, see Section 5.1) is also likely explained by the prevalence and importance of cool spots in low-mass stars. For classical T Tauri stars, the veiling at ~1 µm is usually larger than can be explained by either most accretion models<fnref n="27" /> or warm disk emission (e.g., Basri & Batalha 1990; Cieza et al. 2005; Fischer et al. 2011; McClure et al. 2013). Similarly, Cottaar et al. (2014) found that most cool stars in the Pleiades are veiled in the H band, with no possibility for any disk or accretion emission. Differences in spot temperatures and coverages likely contribute to and may explain entirely this excess veiling.

6.4. *The Limitation of a Two-temperature Model*

We assume only two-temperature components are present
in the photosphere. In reality, starspots are probably described
by a range of spot temperatures analogous to umbra,
penumbra, and plages on the solar surface. A model
possessing three unique temperatures and two unique fill
factors would better approximate the minor constituents of the
emergent spectrum than our two-temperature model can, but
additional components will always appear to fit data better,
while leading to possible overfitting. The fit quality of the
existing two-temperature fits seems strong through much of
H-band and many K-band orders, so the contribution of minor
photospheric constituents is relatively low. Interestingly, the
K band has several spectral orders that converge toward
hybrid temperatures in the vicinity of ~3500 K. Most of these
orders are flagged as erroneous due to strong metal line
spectral outliers or uncorrected calibration artifacts. A few

<figure><img src="image_19.png" /></figure>

**Figure 13.** Gray points are all Grankin et al. (2008) photometric data, while black points are the subset of observations from the 1999 season. The purple line represents a two-component photometric model described by Grankin et al. (2008), comprising a photospheric temperature of 4040 K (for K7), and a temperature of spotted area of 3030 K (for an M5). In this model, the spot coverage of the observed hemisphere varies from ~36% at maximum brightness, and ~82% at minimum brightness.

orders just show mediocre fits, where no permutation of hot
and cool two-temperature models could describe all of the
variance in the spectrum. These K-band orders could be
examined for a penumbra-like temperature component. We
caution that systematic errors in the models would ultimately
limit the interpretability of models with more than two-
temperature components. Space-based high-precision plane-
tary transit spectrophotometry could offer an avenue to
measure the morphologies and temperature distributions of
starspots for the subset of transiting exoplanet host stars or
eclipsing binaries.

Analyses of optical emission also suggest that a two-
temperature model may be an oversimplification. Figure 13
shows that the V − R colors from Grankin et al. (2008) are
better reproduced if the cooler component is ~3050 K, thereby
contributing more flux at optical wavelengths. The ZDI
imaging of Donati et al. (2014) also indicates a complex
structure of spots that are hotter than what is interpreted here as
the ambient photosphere, while also reproducing the optical
photometry. Despite the need for more complicated emission
structures to simultaneously reproduce all observations, the
broadband SED, spectral shape, and rotational modulation of
spectral indices are reasonably well fit with the two
components presented here. The high- and low-resolution
near-IR spectra, including the H₂O band absorption between
the J and H bands (see Figure 8), require a cool (2700–3000 K)
component, while the global photometry requires a large
covering fraction for this cool component.

6.5. *Limitations of Assumptions in the Spectral Inference*
*Methodology*

We assumed that the PHOENIX synthetic spectra are a good
representation of the observed IGRINS spectrum of LkCa 4, a
young star that possesses a strong magnetic field (Donati et al.
2014). Cool sunspots are associated with locally heightened
magnetic fields on the Sun, so synthetic spectral models
employing negligible magnetic fields will fail to fit magnetic

<footnote n="27">²⁷ Ingleby et al. (2013) invoked multiple accretion flows with a range of densities to explain this red excess.</footnote>
sensitive features. Magnetic fields will have a few effects on the emergent spectrum-line broadening attributable to the Zeeman effect, and pressure broadening attributable to seeing deeper into the photosphere in the line wings than in the line center. The former effect will preferentially impact lines with high magnetic sensitivity (e.g., Johns-Krull et al. 1999; Deen 2013). The latter effect could masquerade as an apparent shift in surface gravity, as seen in mixture models of spotted stars with empirical templates composed of dwarfs and giants (O'Neal et al. 1996). Zeeman amplification can also be present (Basri et al. 1992).

Our tacit assumptions have been that these line-broadening effects are (1) relatively small compared to the gross appearance and disappearance of spectral lines at disparate temperatures, (2) many orders and lines are not affected at all, and (3) the Gaussian process noise model provides some resilience against unknown residual spectrum correlations induced by such under-fitting. We see tentative evidence for physics outside of the pre-computed model grids in the poor fits to the observed spectra in the *K* band. Strong metal lines demonstrate large departures in their best fit [Fe/H], as the naïve fitting strategy futilely attempts to minimize the residuals from Zeeman-broadened lines. No pre-computed model grids including a range of surface magnetic fields are available at the wavelength range and spectral sampling needed to approach these questions within this spectral inference framework.

Uncertainties in the PHOENIX synthetic spectra continuum opacity flow down to uncertainties in the absolute value of the fill factor. Erroneous bulk opacities could be responsible for some of the scatter of fill factor distributions seen across different spectral orders, as in Figures 5 and 4. The large number of spectral orders over which the fill factors are averaged offers resilience against these anomalies.

## 6.6. Evidence for Intra-spectrum RV Jitter

Starspots cause significant jitter in RV surveys for young planets (e.g., Donati et al. 2014; Robertson et al. 2014). The two-temperature model assumed that the hot and cool components shared the same $v \sin i$ and $v_z$. This assumption is accurate for either homogeneously distributed spots, azimuthally symmetric bands of spots, or polar spots. LkCa 4 has clear optical line-profile variations (Nguyen et al. 2012; Donati et al. 2014). We should be able to detect an anticorrelation between the line center positions of spectral lines arising from the ambient photosphere and lines arising from cool starspots, since there is a zero-sum competition for solid-angle on the recessional and advancing stellar limbs. We found tentative, albeit weak, evidence for such anticorrelation in $v_z$ across spectral orders dominated by cool photosphere and those dominated by hot photosphere. For example, the cool photosphere dominated spectral orders 100, 102, and 104 (1.7–1.8 µm) show the largest median $v_z \sim 14.5–17.0 \text{ km s}^{-1}$, while the other orders that possess a mix of primarily hot- and some cool photosphere features demonstrate $v_z \sim 9–13 \text{ km s}^{-1}$. An independent analysis of the RV of the IGRINS spectrum yields a barycentric RV of $14.6 \pm 0.2 \text{ km s}^{-1}$. The direct emission of starspots could constitute a heretofore unaccounted source of correlations in conventional RV signal processing.

RV planet searches in the infrared—where the modulation of RV signals is lessened (e.g., Prato et al. 2008; Crockett et al. 2012)—could potentially reduce RV jitter by characterizing and constraining the effect of direct emission from

starspots, in addition to the flux deficits imbed in optical line profiles. Direct spectral emission from cool starspots is likely much weaker in sources with less starspot coverage than LkCa 4. As demands on RV precision increase, the characterization and modeling of direct emission from starspots should, in principle, decrease RV jitter attributable to starspots in the near-IR spectra of any spotted star.

## 7. Conclusions

This paper establishes that the optical and near-IR spectra of the WTTS LkCa 4 are produced by multiple temperature components. When the temperature distribution is simplified to two temperatures, the photospheric model that best fits the spectrum includes a hot (4100 K) photospheric region that covers 20% of the stellar surface and a cool (2700–3000 K) photospheric that covers 80% of the stellar surface. The previously established V-band brightness variability of 0.8 mag is caused by the rotational modulation of the covering fraction of the hot photospheric component. This spectral modeling leads to an improved effective temperature measurement of LkCa 4. For LkCa 4 and other young low-mass stars, the standard approach of inferring masses and ages of young stars is inappropriate and will lead to large uncertainties, both because effective temperatures are systematically too hot and because most models of pre-main-sequence evolution neglect spots.

New results from Fang et al. (2016) show that large starspot coverage fractions of 50% are common even among stars in the Pleiades, which at ~125 Myr are much older than the 1–10 Myr old stars in the Taurus Molecular Cloud. LkCa 4 breaks the apparent ceiling at the 50% coverage fraction of cool starspots, though it is not clear how to classify stars for which the minority of the surface is purportedly the ambient photosphere. A high fraction of the stellar surface may be covered in spots without inducing any photometric modulation, if the spots are located in circumpolar regions that always face toward our line of sight, in longitudinally symmetric bands, or in many isotropic small spots.

Coupled with the results of Fang et al. (2016) and Covey et al. (2016), the results shown here suggest that the contribution of starspots to pre-main-sequence stellar evolution have been systematically underestimated. Some cool stars in the Pleiades have a large excess in Li abundance, some have inflated radii, some are fast rotators, and some have spots (e.g., Somers & Pinsonneault 2015b). While the source overlap for these properties is unclear, these anomalies may all be explained if the pre-main-sequence evolution of some stars is strongly affected by spots (Somers & Pinsonneault 2015a). LkCa 4 and other young stars with similar properties may help to connect the properties of anomalous Pleiades stars with earlier stages of stellar evolution.

The spectral inference technique described in this paper lays the foundation for future studies that seek accurate constraints on multiple photospheric components. The near-IR has been shown to be the most sensitive region to search for cool emission from starspots. Starspots and disk or accretion veiling can be modeled in this framework, and combining the two will enable studies of the large number of classical T-Tauri stars. Applying this technique to large ensembles of T-Tauri stars across a wide range of mass and age will offer an observational picture of starspot evolution that can revise the current understanding of pre-main-sequence stellar evolution.
We thank the anonymous referee, Chris Kochanek, and Daniel T Jaffe for valuable comments on the draft. M.G.-S. and G.J.H. also thank Steve Saar, Lisa Prato, and Lynne Hillenbrand for valuable discussions about spots.

Portions of the data, Python code, Jupyter notebooks, Starfish configuration files, MCMC samples, and revision history of this manuscript are available on the project's GitHub repository at https://github.com/BrownDwarf/welter.

M.G.-S. and G.J.H. are supported by general grant 11473005 awarded by the National Science Foundation of China. The ESPaDOnS observations are supported by the contribution to the MaTYSSE Large Project on CFHT obtained through the Telescope Access Program (TAP), which has been funded by the “the Strategic Priority Research Program—The Emergence of Cosmological Structures” of the Chinese Academy of Sciences (Grant No.11 XDB09000000) and the Special Fund for Astronomy from the Ministry of Finance. TW-SH is supported by the DOE Computational Science Graduate Fellowship, grant number DE-FG02-97ER25308.

We thank LCOGT and its staff for their continued support of ASAS-SN. ASAS-SN is supported by NSF grant AST-1515927. Development of ASAS-SN has been supported by NSF grant AST-0908816, the Center for Cosmology and AstroParticle Physics at the Ohio State University, the Mt. Cuba Astronomical Foundation, and by George Skestos.

B.J.S. is supported by NASA through Hubble Fellowship grant HST-HF-51348.001 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS 5-26555.

This work used the Immersion Grating Infrared Spectrograph (IGRINS) that was developed under a collaboration between the University of Texas at Austin and the Korea Astronomy and Space Science Institute (KASI) with the financial support of the US National Science Foundation under grant ASTR1229522, of the University of Texas at Austin, and of the Korean GMT Project of KASI.

This research has made use of NASA's Astrophysics Data System. The research is based on data from the OMC Archive at CAB (INTA-CSIC), pre-processed by ISDC. We acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

*Facilities:* Smith (IGRINS), AAVSO, CFHT (ESPaDOnS), INTEGRAL (OMC), ASAS, CrAO:1.25m, ARC (TripleSpec), Hale (DBSP), Gaia.

*Software:* pandas (McKinney 2010), emcee (Foreman-Mackey et al. 2013), matplotlib (Hunter 2007), numpy (Walt et al. 2011), scipy (Jones et al. 2001), ipython (Pérez & Granger 2007), gatspy (Vanderplas 2015), starfish (Czekala et al. 2015), seaborn (Waskom et al. 2014).

# Appendix A
## Further Methodological Details

When the synthetic stellar models used in Starfish are in absolute flux units, our mixture model is complete as stated in

Equation (1). However, C15 employ standardized $f_\lambda$ when constructing the forward model M:

$$ \bar{f}_{\lambda} = \frac{f_{\lambda}}{\int_{0}^{\infty} f_{\lambda} d\lambda} = \frac{f_{\lambda}}{f} \qquad (4) $$

The reason for standardizing the flux is a matter of practicality: the number of PCA eigenspectra components in the spectral emulator scales steeply with the pixel-to-pixel variance of $f_\lambda$ ($\{\theta_*\}^{\text{grid}}$), increasing the computational cost of spectral emulation. The choice to standardize fluxes makes no difference for modeling a single photospheric component. However, for two-component photosphere models, the relative flux of the two model spectra needs to be accounted for to get an accurate estimate of the areal coverage fraction of the cool spots, $f_{\text{cool}} \equiv \Omega_{\text{cool}}/(\Omega_{\text{hot}} + \Omega_{\text{cool}})$. Thus we scale the mixture model in the following way:

$$ f_{\lambda, \text{mix}} = f_{\text{hot}} \bar{f}_{\lambda, \text{hot}} \times \Omega_{\text{hot}} + f_{\text{cool}} \bar{f}_{\lambda, \text{cool}} \times \Omega_{\text{cool}}, \qquad (5) $$

$$ q = q(\theta_*) \equiv \frac{f_{\text{cool}}(T_{\text{cool}}...)}{f_{\text{hot}}(T_{\text{hot}}...)} \qquad (6) $$

$$ f'_{\lambda, \text{mix}} = \bar{f}_{\lambda, \text{hot}} \times \Omega_{\text{hot}} + q \bar{f}_{\lambda, \text{cool}} \times \Omega_{\text{cool}} \qquad (7) $$

where the prime symbol in the final line indicates a re-standardized mixture model flux, where the relative fluxes of the model components are now correctly scaled.

We are then tasked with computing an estimator, $\hat{f}(\theta_*)$, for estimating the scale factor $q$ in-between model grid points. One robust approach would be to follow C15 by training a Gaussian process on the $f(\{\theta_*\}^{\text{grid}})$. Instead, we simply linearly interpolated between the model grid points. Interpolation can cause pile-up near model grid points, as noted in Cottaar et al. (2014), which motivated the spectral emulation procedure in C15. We assume the interpolation of $\hat{f}$ is smooth enough that we will not see such pileups and, even if we did see pileups, they would mostly be discernable as kinks in the distribution of samples in the starspot coverage fraction $f_{\text{cool}}$. One drawback of our estimator method compared to the Gaussian Process regression method is that we do not propagate the uncertainty associated with the absolute flux ratio interpolation into our estimate of $f_{\text{cool}}$. We assume this uncertainty is relatively small and can be ignored.

The mixture model is a linear operation with all the same stellar extrinsic parameters, so we can re-use all of the same post-processed eigenspectra $\tilde{\Sigma}$, mean spectra $\tilde{\xi}_\mu$, and variance spectra $\tilde{\xi}_\sigma$, with the tildes representing all post processing except the $\Omega$ scaling. We calculate two sets of eigenspectra weights $w_{\in(\text{hot,cool})}$, and their associated mean and covariances following the appendix of C15, and yielding

$$ M'_{\text{mix}} = \Omega_{\text{hot}}(\tilde{\xi}_{\mu} + X\mu_{w,\text{hot}}) + q\Omega_{\text{cool}}(\tilde{\xi}_{\mu} + X\mu_{w,\text{cool}}), \qquad (8) $$

$$ C'_{\text{mix}} = \Omega_{\text{hot}}^2 X \Sigma_{w,\text{hot}} X^T + q \Omega_{\text{cool}}^2 X \Sigma_{w,\text{cool}} X^T + C. \qquad (9) $$

# Appendix B
## Further Examples of Full Spectrum Fitting

Figures 14 and 15 show spectra from 42 IGRINS orders.
<figure><img src="image_20.png" /><figcaption>**Figure 14.** IGRINS Orders 122–100, with panels arranged with the shortest wavelength in the upper left corner with central wavelength decreasing toward the bottom of the leftmost column, then decreasing through the subsequent columns. The y-axis is on a logarithmic scale. The red line is the cool photosphere while the blue line is the hot photosphere. The purple line is the composite mixture model.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>**Figure 15.** IGRINS Orders 93–73, with the same layout, ordering, and colors as Figure 14. The y-axis is on a logarithmic scale.</figcaption></figure>

## Appendix C

### Full Table of IGRINS Best Fits

Table 3 lists all the IGRINS spectral order results.

## Appendix D

### Previous Work

Table 4 lists measurements of LkCa 4 from previous studies, compared to this work.
<table><caption>**Table 3**

Results from One- and Two-component IGRINS Spectral Fitting</caption><thead><tr><th rowspan="2">Order</th><th rowspan="2">Instrument</th><th rowspan="2">λ<sub>1</sub> - λ<sub>2</sub><br>(Å)</th><th colspan="1">One Comp.</th><th colspan="3">Two Components</th></tr><tr><th>T<sub>eff</sub><br>(K)</th><th>T<sub>hot</sub><br>(K)</th><th>T<sub>cool</sub><br>(K)</th><th>f<sub>cool</sub></th></tr></thead><tbody><tr><td>0</td><td>ESPaDoNs</td><td>5161-5200</td><td>3734<sup>+183</sup><sub>-63</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>1</td><td>ESPaDoNs</td><td>5199-5263</td><td>3991<sup>+70</sup><sub>-81</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>2</td><td>ESPaDoNs</td><td>5262-5328</td><td>3990<sup>+196</sup><sub>-140</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>3</td><td>ESPaDoNs</td><td>5327-5395</td><td>3923<sup>+72</sup><sub>-77</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>4</td><td>ESPaDoNs</td><td>5394-5465</td><td>3812<sup>+54</sup><sub>-85</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>5</td><td>ESPaDoNs</td><td>5464-5537</td><td>4165<sup>+59</sup><sub>-61</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>6</td><td>ESPaDoNs</td><td>5536-5612</td><td>4129<sup>+37</sup><sub>-38</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>7</td><td>ESPaDoNs</td><td>5611-5689</td><td>4057<sup>+64</sup><sub>-52</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>8</td><td>ESPaDoNs</td><td>5688-5770</td><td>4038<sup>+50</sup><sub>-44</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>9</td><td>ESPaDoNs</td><td>5769-5854</td><td>4212<sup>+71</sup><sub>-64</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>10</td><td>ESPaDoNs</td><td>5853-5941</td><td>4119<sup>+52</sup><sub>-46</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>11</td><td>ESPaDoNs</td><td>5940-6032</td><td>4091<sup>+83</sup><sub>-63</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>12</td><td>ESPaDoNs</td><td>6031-6127</td><td>4108<sup>+46</sup><sub>-42</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>13</td><td>ESPaDoNs</td><td>6126-6225</td><td>3941<sup>+28</sup><sub>-24</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>14</td><td>ESPaDoNs</td><td>6225-6329</td><td>4012<sup>+33</sup><sub>-42</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>15</td><td>ESPaDoNs</td><td>6328-6436</td><td>4238<sup>+49</sup><sub>-40</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>16</td><td>ESPaDoNs</td><td>6435-6549</td><td>3880<sup>+61</sup><sub>-76</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>17</td><td>ESPaDoNs</td><td>6548-6667</td><td>3874<sup>+72</sup><sub>-70</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>18</td><td>ESPaDoNs</td><td>6666-6791</td><td>4054<sup>+42</sup><sub>-32</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>20</td><td>ESPaDoNs</td><td>6920-7057</td><td>3907<sup>+32</sup><sub>-20</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>21</td><td>ESPaDoNs</td><td>7056-7201</td><td>3910<sup>+25</sup><sub>-20</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>22</td><td>ESPaDoNs</td><td>7200-7352</td><td>3878<sup>+28</sup><sub>-15</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>23</td><td>ESPaDoNs</td><td>7351-7512</td><td>4085<sup>+51</sup><sub>-66</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>26</td><td>ESPaDoNs</td><td>7858-8051</td><td>3852<sup>+24</sup><sub>-24</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>29</td><td>ESPaDoNs</td><td>8473-8707</td><td>3504<sup>+49</sup><sub>-46</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>30</td><td>ESPaDoNs</td><td>8706-8954</td><td>3821<sup>+33</sup><sub>-29</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>123</td><td>IGRINS</td><td>14519-14721</td><td>...</td><td>4182<sup>+107</sup><sub>-75</sub></td><td>2773<sup>+254</sup><sub>-68</sub></td><td>0.82<sup>+0.03</sup><sub>-0.72</sub></td></tr><tr><td>122</td><td>IGRINS</td><td>14634-14837</td><td>...</td><td>4185<sup>+187</sup><sub>-137</sub></td><td>2819<sup>+209</sup><sub>-105</sub></td><td>0.78<sup>+0.04</sup><sub>-0.22</sub></td></tr><tr><td>121</td><td>IGRINS</td><td>14750-14956</td><td>...</td><td>4084<sup>+153</sup><sub>-144</sub></td><td>2762<sup>+83</sup><sub>-58</sub></td><td>0.79<sup>+0.04</sup><sub>-0.05</sub></td></tr><tr><td>120</td><td>IGRINS</td><td>14868-15076</td><td>...</td><td>4027<sup>+76</sup><sub>-112</sub></td><td>2739<sup>+94</sup><sub>-37</sub></td><td>0.69<sup>+0.05</sup><sub>-0.08</sub></td></tr><tr><td>119</td><td>IGRINS</td><td>14988-15198</td><td>...</td><td>4093<sup>+123</sup><sub>-132</sub></td><td>2719<sup>+67</sup><sub>-18</sub></td><td>0.76<sup>+0.04</sup><sub>-0.05</sub></td></tr><tr><td>118</td><td>IGRINS</td><td>15111-15322</td><td>...</td><td>4124<sup>+142</sup><sub>-131</sub></td><td>2730<sup>+83</sup><sub>-30</sub></td><td>0.74<sup>+0.08</sup><sub>-0.08</sub></td></tr><tr><td>117</td><td>IGRINS</td><td>15235-15448</td><td>4034<sup>+147</sup><sub>-122</sub></td><td>4105<sup>+63</sup><sub>-73</sub></td><td>2731<sup>+70</sup><sub>-29</sub></td><td>0.73<sup>+0.05</sup><sub>-0.06</sub></td></tr><tr><td>116</td><td>IGRINS</td><td>15362-15577</td><td>4239<sup>+55</sup><sub>-73</sub></td><td>4121<sup>+147</sup><sub>-146</sub></td><td>2795<sup>+169</sup><sub>-85</sub></td><td>0.80<sup>+0.04</sup><sub>-0.09</sub></td></tr><tr><td>115</td><td>IGRINS</td><td>15490-15707</td><td>3718<sup>+94</sup><sub>-130</sub></td><td>3914<sup>+176</sup><sub>-151</sub></td><td>2769<sup>+298</sup><sub>-64</sub></td><td>0.55<sup>+0.14</sup><sub>-0.44</sub></td></tr><tr><td>114</td><td>IGRINS</td><td>15622-15841</td><td>...</td><td>4021<sup>+159</sup><sub>-122</sub></td><td>2751<sup>+147</sup><sub>-43</sub></td><td>0.64<sup>+0.10</sup><sub>-0.13</sub></td></tr><tr><td>113</td><td>IGRINS</td><td>15755-15976</td><td>...</td><td>3874<sup>+265</sup><sub>-170</sub></td><td>2779<sup>+193</sup><sub>-75</sub></td><td>0.77<sup>+0.08</sup><sub>-0.11</sub></td></tr><tr><td>112</td><td>IGRINS</td><td>15891-16114</td><td>...</td><td>3986<sup>+268</sup><sub>-206</sub></td><td>2943<sup>+425</sup><sub>-198</sub></td><td>0.53<sup>+0.16</sup><sub>-0.17</sub></td></tr><tr><td>111</td><td>IGRINS</td><td>16030-16255</td><td>4101<sup>+100</sup><sub>-116</sub></td><td>...</td><td>...</td><td>...</td></tr><tr><td>110</td><td>IGRINS</td><td>16171-16398</td><td>4233<sup>+59</sup><sub>-104</sub></td><td>4180<sup>+264</sup><sub>-193</sub></td><td>2777<sup>+184</sup><sub>-68</sub></td><td>0.83<sup>+0.06</sup><sub>-0.17</sub></td></tr><tr><td>109</td><td>IGRINS</td><td>16314-16543</td><td>3933<sup>+126</sup><sub>-97</sub></td><td>0<sup>+0</sup><sub>-0</sub></td><td>0<sup>+0</sup><sub>-0</sub></td><td>0.50<sup>+0.00</sup><sub>-0.00</sub></td></tr><tr><td>108</td><td>IGRINS</td><td>16461-16692</td><td>4150<sup>+100</sup><sub>-176</sub></td><td>4248<sup>+183</sup><sub>-193</sub></td><td>2802<sup>+339</sup><sub>-84</sub></td><td>0.71<sup>+0.12</sup><sub>-0.35</sub></td></tr><tr><td>107</td><td>IGRINS</td><td>16610-16843</td><td>3750<sup>+186</sup><sub>-112</sub></td><td>3813<sup>+521</sup><sub>-177</sub></td><td>3237<sup>+437</sup><sub>-486</sub></td><td>0.38<sup>+0.37</sup><sub>-0.36</sub></td></tr><tr><td>106</td><td>IGRINS</td><td>16762-16997</td><td>4225<sup>+68</sup><sub>-92</sub></td><td>4324<sup>+161</sup><sub>-168</sub></td><td>2882<sup>+1227</sup><sub>-175</sub></td><td>0.53<sup>+0.19</sup><sub>-0.53</sub></td></tr><tr><td>105</td><td>IGRINS</td><td>16917-17155</td><td>&lt;3711</td><td>3930<sup>+209</sup><sub>-152</sub></td><td>2763<sup>+212</sup><sub>-59</sub></td><td>0.66<sup>+0.06</sup><sub>-0.10</sub></td></tr><tr><td>104</td><td>IGRINS</td><td>17075-17315</td><td>&lt;3776</td><td>4300<sup>+177</sup><sub>-245</sub></td><td>2870<sup>+115</sup><sub>-130</sub></td><td>0.83<sup>+0.05</sup><sub>-0.05</sub></td></tr><tr><td>103</td><td>IGRINS</td><td>17236-17478</td><td>&lt;3700</td><td>3848<sup>+312</sup><sub>-334</sub></td><td>2972<sup>+199</sup><sub>-234</sub></td><td>0.65<sup>+0.18</sup><sub>-0.30</sub></td></tr><tr><td>102</td><td>IGRINS</td><td>17400-17645</td><td>&lt;3629</td><td>4336<sup>+136</sup><sub>-159</sub></td><td>2751<sup>+113</sup><sub>-48</sub></td><td>0.82<sup>+0.04</sup><sub>-0.05</sub></td></tr><tr><td>101</td><td>IGRINS</td><td>17568-17815</td><td>&lt;3543</td><td>4199<sup>+264</sup><sub>-574</sub></td><td>2916<sup>+158</sup><sub>-158</sub></td><td>0.79<sup>+0.05</sup><sub>-0.23</sub></td></tr><tr><td>100</td><td>IGRINS</td><td>17739-17988</td><td>&lt;3598</td><td>4312<sup>+161</sup><sub>-193</sub></td><td>2765<sup>+135</sup><sub>-61</sub></td><td>0.79<sup>+0.04</sup><sub>-0.06</sub></td></tr><tr><td>99</td><td>IGRINS</td><td>17914-18165</td><td>&lt;3552</td><td>...</td><td>...</td><td>...</td></tr><tr><td>94</td><td>IGRINS</td><td>18855-19117</td><td>...</td><td>4293<sup>+185</sup><sub>-1132</sub></td><td>2899<sup>+304</sup><sub>-150</sub></td><td>0.81<sup>+0.07</sup><sub>-0.39</sub></td></tr><tr><td>93</td><td>IGRINS</td><td>19053-19318</td><td>...</td><td>4044<sup>+353</sup><sub>-378</sub></td><td>2740<sup>+119</sup><sub>-36</sub></td><td>0.81<sup>+0.08</sup><sub>-0.19</sub></td></tr><tr><td>92</td><td>IGRINS</td><td>19256-19524</td><td>...</td><td>3770<sup>+155</sup><sub>-183</sub></td><td>2859<sup>+510</sup><sub>-149</sub></td><td>0.56<sup>+0.18</sup><sub>-0.21</sub></td></tr><tr><td>91</td><td>IGRINS</td><td>19462-19734</td><td>3464<sup>+70</sup><sub>-67</sub></td><td>3962<sup>+389</sup><sub>-426</sub></td><td>3126<sup>+857</sup><sub>-325</sub></td><td>0.68<sup>+0.18</sup><sub>-0.56</sub></td></tr><tr><td>90</td><td>IGRINS</td><td>19674-19948</td><td>3594<sup>+75</sup><sub>-70</sub></td><td>3907<sup>+251</sup><sub>-211</sub></td><td>3463<sup>+90</sup><sub>-161</sub></td><td>0.75<sup>+0.10</sup><sub>-0.28</sub></td></tr></tbody></table>
<table><caption>**Table 3**
(Continued)</caption><thead><tr><th rowspan="2">Order</th><th rowspan="2">Instrument</th><th rowspan="2">λ<sub>1</sub> - λ<sub>2</sub><br>(Å)</th><th colspan="1">One Comp.</th><th colspan="3">Two Components</th></tr><tr><th>T<sub>eff</sub><br>(K)</th><th>T<sub>hot</sub><br>(K)</th><th>T<sub>cool</sub><br>(K)</th><th>f<sub>cool</sub></th></tr></thead><tbody><tr><td>89</td><td>IGRINS</td><td>19890-20168</td><td>3557<sup>+82</sup><sub>-72</sub></td><td>3706<sup>+254</sup><sub>-151</sub></td><td>3412<sup>+183</sup><sub>-617</sub></td><td>0.45<sup>+0.20</sup><sub>-0.23</sub></td></tr><tr><td>88</td><td>IGRINS</td><td>20112-20392</td><td>3500<sup>+61</sup><sub>-60</sub></td><td>3917<sup>+208</sup><sub>-129</sub></td><td>2812<sup>+250</sup><sub>-99</sub></td><td>0.75<sup>+0.09</sup><sub>-0.06</sub></td></tr><tr><td>87</td><td>IGRINS</td><td>20338-20622</td><td>3314<sup>+77</sup><sub>-90</sub></td><td>3709<sup>+305</sup><sub>-207</sub></td><td>2848<sup>+167</sup><sub>-119</sub></td><td>0.75<sup>+0.09</sup><sub>-0.24</sub></td></tr><tr><td>86</td><td>IGRINS</td><td>20570-20857</td><td>3544<sup>+51</sup><sub>-48</sub></td><td>3915<sup>+475</sup><sub>-275</sub></td><td>3135<sup>+281</sup><sub>-345</sub></td><td>0.68<sup>+0.15</sup><sub>-0.18</sub></td></tr><tr><td>85</td><td>IGRINS</td><td>20808-21098</td><td>3315<sup>+99</sup><sub>-107</sub></td><td>4220<sup>+252</sup><sub>-322</sub></td><td>2828<sup>+205</sup><sub>-112</sub></td><td>0.85<sup>+0.05</sup><sub>-0.10</sub></td></tr><tr><td>84</td><td>IGRINS</td><td>21051-21345</td><td>3475<sup>+62</sup><sub>-56</sub></td><td>3681<sup>+444</sup><sub>-191</sub></td><td>3414<sup>+125</sup><sub>-130</sub></td><td>0.68<sup>+0.17</sup><sub>-0.12</sub></td></tr><tr><td>83</td><td>IGRINS</td><td>21300-21597</td><td>3234<sup>+81</sup><sub>-108</sub></td><td>4334<sup>+144</sup><sub>-235</sub></td><td>2896<sup>+166</sup><sub>-137</sub></td><td>0.90<sup>+0.03</sup><sub>-0.03</sub></td></tr><tr><td>82</td><td>IGRINS</td><td>21555-21856</td><td>3589<sup>+2</sup><sub>-8</sub></td><td>3535<sup>+582</sup><sub>-258</sub></td><td>2964<sup>+297</sup><sub>-241</sub></td><td>0.62<sup>+0.21</sup><sub>-0.38</sub></td></tr><tr><td>81</td><td>IGRINS</td><td>21817-22121</td><td>3466<sup>+117</sup><sub>-87</sub></td><td>3536<sup>+146</sup><sub>-124</sub></td><td>3373<sup>+152</sup><sub>-150</sub></td><td>0.34<sup>+0.16</sup><sub>-0.10</sub></td></tr><tr><td>80</td><td>IGRINS</td><td>22085-22393</td><td>3157<sup>+132</sup><sub>-124</sub></td><td>3543<sup>+404</sup><sub>-273</sub></td><td>2797<sup>+336</sup><sub>-88</sub></td><td>0.69<sup>+0.16</sup><sub>-0.21</sub></td></tr><tr><td>79</td><td>IGRINS</td><td>22360-22671</td><td>3437<sup>+79</sup><sub>-60</sub></td><td>3511<sup>+138</sup><sub>-92</sub></td><td>3383<sup>+122</sup><sub>-120</sub></td><td>0.30<sup>+0.13</sup><sub>-0.10</sub></td></tr><tr><td>78</td><td>IGRINS</td><td>22643-22957</td><td>3197<sup>+43</sup><sub>-43</sub></td><td>3808<sup>+321</sup><sub>-231</sub></td><td>2727<sup>+159</sup><sub>-26</sub></td><td>0.78<sup>+0.07</sup><sub>-0.06</sub></td></tr><tr><td>77</td><td>IGRINS</td><td>22932-23251</td><td>3273<sup>+52</sup><sub>-49</sub></td><td>4221<sup>+254</sup><sub>-392</sub></td><td>2716<sup>+273</sup><sub>-15</sub></td><td>0.85<sup>+0.03</sup><sub>-0.12</sub></td></tr><tr><td>76</td><td>IGRINS</td><td>23230-23552</td><td>3321<sup>+54</sup><sub>-53</sub></td><td>3428<sup>+317</sup><sub>-116</sub></td><td>3123<sup>+193</sup><sub>-362</sub></td><td>0.46<sup>+0.22</sup><sub>-0.21</sub></td></tr><tr><td>75</td><td>IGRINS</td><td>23535-23861</td><td>...</td><td>3649<sup>+528</sup><sub>-218</sub></td><td>2956<sup>+394</sup><sub>-222</sub></td><td>0.62<sup>+0.17</sup><sub>-0.17</sub></td></tr><tr><td>74</td><td>IGRINS</td><td>23849-24179</td><td>...</td><td>3523<sup>+133</sup><sub>-94</sub></td><td>3053<sup>+388</sup><sub>-308</sub></td><td>0.34<sup>+0.20</sup><sub>-0.19</sub></td></tr><tr><td>73</td><td>IGRINS</td><td>24172-24505</td><td>...</td><td>4119<sup>+345</sup><sub>-433</sub></td><td>3029<sup>+244</sup><sub>-294</sub></td><td>0.80<sup>+0.08</sup><sub>-0.13</sub></td></tr><tr><td>72</td><td>IGRINS</td><td>24503-24840</td><td>...</td><td>3897<sup>+398</sup><sub>-184</sub></td><td>2881<sup>+425</sup><sub>-157</sub></td><td>0.74<sup>+0.11</sup><sub>-0.26</sub></td></tr></tbody></table>

Note. The uncertainties are the 5th and 95th percentiles of the marginalized samples, corresponding roughly to “2σ” error bars.

<table><caption>**Table 4**
Previous Studies of LkCa 4</caption><thead><tr><th>Reference</th><th>Band(s)</th><th>R<br>(λ/δλ)</th><th>SpT</th><th>T<sub>eff</sub><br>(K)</th><th>log g</th><th>A<sub>V</sub></th><th>ν sin i<br>(km s<sup>-1</sup>)</th><th>ν<sub>z</sub><br>(km s<sup>-1</sup>)</th><th>log L/L<sub>☉</sub></th></tr></thead><tbody><tr><td>Herbig et al. (1986)</td><td>V</td><td>100,000</td><td>K7 V</td><td>...</td><td>...</td><td>0.68</td><td>26.1 ± 2.4</td><td>13 ± 4</td><td>0.02</td></tr><tr><td>Hartmann et al. (1987)</td><td>U</td><td>100,000</td><td>...</td><td>...</td><td>...</td><td>...</td><td>26.1 ± 2.4</td><td>16.9 ± 2.6</td><td>...</td></tr><tr><td>Downes & Keyes (1988)</td><td>V</td><td>13 Å</td><td>Me</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>Strom et al. (1989a)</td><td>IRAS</td><td>...</td><td>K7:V</td><td>...</td><td>...</td><td>0.95</td><td>...</td><td>...</td><td>0.04</td></tr><tr><td>Strom et al. (1989b)</td><td>V</td><td>3500</td><td>K7:V</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>Schauffer et al. (1991)</td><td>V R I</td><td>500</td><td>M1</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr><tr><td>Strom & Strom (1994)</td><td>V</td><td>...</td><td>K7</td><td>4000</td><td>...</td><td>1.25</td><td>...</td><td>...</td><td>0.06</td></tr><tr><td>Kenyon & Hartmann (1995)</td><td>U-IRAS</td><td>...</td><td>K7</td><td>4060</td><td>...</td><td>0.69</td><td>...</td><td>...</td><td>-0.08</td></tr><tr><td>Hartigan et al. (1995)</td><td>V</td><td>25,000</td><td>K7</td><td>4000</td><td>...</td><td>0.68</td><td>...</td><td>...</td><td>-0.07</td></tr><tr><td>White & Ghez (2001)</td><td>V</td><td>...</td><td>K7</td><td>...</td><td>...</td><td>1.21</td><td>...</td><td>...</td><td>-0.02</td></tr><tr><td>Nguyen et al. (2009)</td><td>BV R I</td><td>60,000</td><td>...</td><td>...</td><td>...</td><td>...</td><td>30 ± 2</td><td>...</td><td>...</td></tr><tr><td>Nguyen et al. (2012)</td><td>BV R I</td><td>60,000</td><td>...</td><td>...</td><td>...</td><td>...</td><td>30 ± 2</td><td>16.0 ± 4.0</td><td>...</td></tr><tr><td>Grankin (2013)</td><td>BV R I</td><td>...</td><td>K7</td><td>4040</td><td>...</td><td>0.54</td><td>26.1</td><td>...</td><td>-0.13</td></tr><tr><td>Donati et al. (2014)</td><td>U BV R I</td><td>68,000</td><td>...</td><td>4100 ± 50</td><td>3.8 ± 0.1</td><td>0.68 ± 0.15</td><td>28.0 ± 0.5</td><td>16.8 ± 0.1</td><td>-0.04 ± 0.11</td></tr><tr><td>Herczeg & Hillenbrand (2014)</td><td>U BV R I</td><td>700</td><td>M1.3</td><td>3670</td><td>...</td><td>0.35</td><td>...</td><td>...</td><td>-0.29</td></tr><tr><td>This Work: Best</td><td>H K</td><td>45,000</td><td>...</td><td>3180</td><td>~3.8</td><td>0.3</td><td>29 ± 3</td><td>13 ± 3</td><td>-0.26</td></tr><tr><td>This Work: Alternate</td><td>H K</td><td>45,000</td><td>...</td><td>3320</td><td>~3.8</td><td>...</td><td>29 ± 3</td><td>13 ± 3</td><td>~0.26</td></tr></tbody></table>

References

Alfonso-Garzón, J., Domingo, A., Mas-Hesse, J. M., & Giménez, A. 2012, A&A, 548, A79

Allard, F. 2014, in IAU Symp. 299, Exploring the Formation and Evolution of Planetary Systems, ed. M. Booth, B. C. Matthews, & J. R. Graham (Cambridge: Cambridge Univ. Press), 271

Andrews, S. M., & Williams, J. P. 2005, ApJ, 631, 1134

Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42

Bary, J. S., & Petersen, M. S. 2014, ApJ, 792, 64

Basri, G., & Batalha, C. 1990, ApJ, 363, 654

Basri, G., Marcy, G. W., & Valenti, J. A. 1992, ApJ, 390, 622

Berdyugina, S. V. 2005, LRSP, 2, 8

Bouvier, J., & Appenzeller, I. 1992, A&AS, 92, 481

Bouvier, J., Cabrit, S., Fernandez, M., Martin, E. L., & Matthews, J. M. 1993, A&AS, 101, 485

Bouvier, J., Covino, E., Kovo, O., et al. 1995, A&A, 299, 89

Buckle, J. V., Drabek-Maunder, E., Greaves, J., et al. 2015, MNRAS, 449, 2472

Cauley, P. W., Johns-Krull, C. M., Hamilton, C. M., & Lockhart, K. 2012, ApJ, 756, 68

Chugainov, P. F. 1976, IzKry, 55, 94

Cieza, L. A., Kessler-Silacci, J. E., Jaffe, D. T., Harvey, P. M., & Evans, N. J., II 2005, ApJ, 635, 422

Cody, A. M., Stauffer, J., Baglin, A., et al. 2014, AJ, 147, 82

Cottaar, M., Covey, K. R., Meyer, M. R., et al. 2014, ApJ, 794, 125
Covey, K. R., Agüeros, M. A., Law, N. M., et al. 2016, ApJ, 822, 81

Crockett, C. J., Mahmud, N. I., Prato, L., et al. 2012, ApJ, 761, 164

Cushing, M. C., Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 362

Czekala, I., Andrews, S. M., Mandel, K. S., Hogg, D. W., & Green, G. M. 2015, ApJ, 812, 128

Daemgen, S., Bonavita, M., Jayawardhana, R., Lafrenière, D., & Janson, M. 2015, ApJ, 799, 155

Davenport, J. R. A., Hebb, L., & Hawley, S. L. 2015, ApJ, 806, 212

Debes, J. H., Jang-Condell, H., Weinberger, A. J., Roberge, A., & Schneider, G. 2013, ApJ, 771, 45

Deen, C. P. 2013, AJ, 146, 51

Desert, J.-M., Charbonneau, D., Demory, B.-O., et al. 2011, ApJS, 197, 14

Donati, J.-F., Hébrard, E., Hussain, G., et al. 2014, MNRAS, 444, 3220

Donati, J.-F., & Landstreet, J. D. 2009, ARA&A, 47, 333

Downes, R. A., & Keyes, C. D. 1988, AJ, 96, 777

Edwards, S., Fischer, W., Hillenbrand, L., & Kwan, J. 2006, ApJ, 646, 319

Espaillat, C., D'Alessio, P., Hernández, J., et al. 2010, ApJ, 717, 441

Fang, X.-S., Zhao, G., Zhao, J.-K., Chen, Y.-Q., & Bharat Kumar, Y. 2016, arXiv:1608.05452

Feiden, G. A. 2016, arXiv:1604.08036

Fischer, W., Edwards, S., Hillenbrand, L., & Kwan, J. 2011, ApJ, 730, 73

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, M. 2013, PASP, 125, 306

Frasca, A., Biazzo, K., Lanzafame, A. C., et al. 2015, A&A, 575, A4

Furlan, E., Hartmann, L., Calvet, N., et al. 2006, ApJS, 165, 568

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016, arXiv:1609.04172

Grankin, K. N. 1998, AstL, 24, 497

Grankin, K. N. 1999, AstL, 25, 526

Grankin, K. N. 2013, AstL, 39, 251

Grankin, K. N., Bouvier, J., Herbst, W., & Melnikov, S. Y. 2008, A&A, 479, 827

Grankin, K. N., Melnikov, S. Y., Bouvier, J., Herbst, W., & Shevchenko, V. S. 2007, A&A, 461, 183

Harrison, T. E., Coughlin, J. L., Ule, N. M., & Lopez-Morales, M. 2011, AJ, 143, 4

Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736

Hartmann, L. 2001, AJ, 121, 1030

Hartmann, L., Megeath, S. T., Allen, L., et al. 2005, ApJ, 629, 881

Hartmann, L. W., Soderblom, D. R., & Stauffer, J. R. 1987, AJ, 93, 907

Herbig, G. H., Vrba, F. J., & Rydgren, A. E. 1986, AJ, 91, 575

Herbst, W. 1989, AJ, 98, 2268

Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 1994, AJ, 108, 1906

Herczeg, G. J., & Hillenbrand, L. A. 2014, ApJ, 786, 97

Herczeg, G. J., & Hillenbrand, L. A. 2015, ApJ, 808, 23

Hunter, J. D. 2007, CSE, 9, 90

Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6

Ingleby, L., Calvet, N., Herczeg, G., et al. 2013, ApJ, 767, 112

Itoh, Y., Tamura, M., Hayashi, M., et al. 2008, PASJ, 60, 209

Ivezić, Ž., Connelly, A. J., VanderPlas, J. T., & Gray, A. 2014, Statistics, Data Mining, and Machine Learning in Astronomy (Princeton, NJ: Princeton Univ. Press)

Jackson, R. J., & Jeffries, R. D. 2014, MNRAS, 441, 2111

Jeffries, R. D., Jackson, R. J., Cottaar, M., et al. 2014, A&A, 563, A94

Johns-Krull, C. M. 2007, ApJ, 664, 975

Johns-Krull, C. M., Valenti, J. A., & Koresko, C. 1999, ApJ, 516, 900

Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open source scientific tools for Python, (Online; accessed 2016-08-02)

Kafka, S. 2016, Observations from the AAVSO International Database, http://www.aavso.org, (Online; accessed 2016 July)

Karr, J. L., Ohashi, N., Kudo, T., & Tamura, M. 2010, AJ, 139, 1015

Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117

Kirby, E. N. 2011, PASP, 123, 531

Kirby, E. N., Guo, M., Zhang, A. J., et al. 2015, ApJ, 801, 125

Kirkpatrick, J. D., Henry, T. J., & McCarthy, D. W., Jr. 1991, ApJS, 77, 417

Koen, C. 2016, MNRAS, 463, 4383

Kraus, A. L., & Hillenbrand, L. A. 2009, ApJ, 703, 1511

Kraus, A. L., Ireland, M. J., Martinache, F., & Hillenbrand, L. A. 2011, ApJ, 731, 8

Landolt, A. U. 1983, AJ, 88, 439

Lee, J.-J. 2015, plp: Version 2.0, Zenodo, doi:10.5281/zenodo.18579

MacDonald, J., & Mullan, D. J. 2013, ApJ, 765, 126

McClure, M. K., Calvet, N., Espaillat, C., et al. 2013, ApJ, 769, 73

McKinney, W. 2010, in Proc. 9th Python in Science Conf., ed. S. van der Walt & J. Millman, 51

Neff, J. E., O'Neal, D., & Saar, S. H. 1995, ApJ, 452, 879

Nguyen, D. C., Brandeker, A., van Kerkwijk, M. H., & Jayawardhana, R. 2012, ApJ, 745, 119

Nguyen, D. C., Jayawardhana, R., van Kerkwijk, M. H., et al. 2009, ApJ, 695, 1648

O'Neal, D., Neff, J. E., Saar, S. H., & Cuntz, M. 2004, AJ, 128, 1802

O'Neal, D., Neff, J. E., Saar, S. H., & Mines, J. K. 2001, AJ, 122, 1954

O'Neal, D., Saar, S. H., & Neff, J. E. 1996, ApJ, 463, 766

Oke, J. B., & Gunn, J. E. 1982, PASP, 94, 586

Park, C., Jaffe, D. T., Yuk, I.-S., et al. 2014, Proc. SPIE, 9147, 1

Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9

Pérez, F., & Granger, B. E. 2007, CSE, 9, 21

Petrov, P. P., Shcherbakov, V. A., Berdyugina, S. V., et al. 1994, A&AS, 107, 9

Pojmański, G. 2004, AN, 325, 553

Prato, L., Huerta, M., Johns-Krull, C. M., et al. 2008, ApJL, 687, L103

Preibisch, T. 2012, RAA, 12, 1

Rizzuto, A. C., Ireland, M. J., Dupuy, T. J., & Kraus, A. L. 2016, ApJ, 817, 164

Robertson, P., Mahadevan, S., Endl, M., & Roy, A. 2014, Sci, 345, 440

Roettenbacher, R. M., Monnier, J. D., Korhonen, H., et al. 2016, Natur, 533, 217

Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Slesnick, C. L., Hillenbrand, L. A., & Carpenter, J. M. 2008, ApJ, 688, 377

Soderblom, D. R., Hillenbrand, L. A., Jeffries, R. D., Mamajek, E. E., & Naylor, T. 2014, in Protostars and Planets VI, ed. H. Beuther et al. (Tucson, AZ: Univ. of Arizona Press), 219

Somers, G., & Pinsonneault, M. H. 2015a, ApJ, 807, 174

Somers, G., & Pinsonneault, M. H. 2015b, MNRAS, 449, 4131

Stassun, K. G., Feiden, G. A., & Torres, G. 2014, NewAR, 60, 1

Stauffer, J., Herter, T., Hamilton, D., et al. 1991, ApJL, 367, L23

Stauffer, J. R., Jones, B. F., Backman, D., et al. 2003, AJ, 126, 833

Stauffer, J. R., Schultz, G., & Kirkpatrick, J. D. 1998, ApJL, 499, L199

Strom, K. M., & Strom, S. E. 1994, ApJ, 424, 237

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989a, AJ, 97, 1451

Strom, K. M., Wilkin, F. P., Strom, S. E., & Seaman, R. L. 1989b, AJ, 98, 1444

Torres, R. M., Loinard, L., Mioduszewski, A. J., et al. 2012, ApJ, 747, 18

Tottle, J., & Mohanty, S. 2015, ApJ, 805, 57

Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389

Vacca, W. D., & Sandell, G. 2011, ApJ, 732, 8

Vanderplas, J. 2015, gatspy: General tools for Astronomical Time Series in Python, Zenodo, doi:10.5281/zenodo.14833

VanderPlas, J. T., & Ivezić, Z. 2015, ApJ, 812, 18

Venuti, L., Bouvier, J., Irwin, J., et al. 2015, A&A, 581, A66

Vrba, F. J., Chugainov, P. F., Weaver, W. B., & Stauffer, J. S. 1993, AJ, 106, 1608

Vrba, F. J., Rydgren, A. E., Chugainov, P. F., Shakovskaia, N. I., & Zak, D. S. 1986, ApJ, 306, 199

Walt, S. v. d., Colbert, S. C., & Varoquaux, G. 2011, CSE, 13, 22

Waskom, M., Botvinnik, O., Hobson, P., et al. 2014, seaborn: v0.5.0, Zenodo, doi:10.5281/zenodo.12710

White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265

Wilson, J. C., Hearty, F., Skrutskie, M. F., et al. 2010, Proc. SPIE, 7735, 77351C

Wilson, J. C., Henderson, C. P., Herter, T. L., et al. 2004, Proc. SPIE, 5492, 1295

Wolk, S. J., & Walter, F. M. 1996, AJ, 111, 2066

Xiao, H. Y., Covey, K. R., Rebull, L., et al. 2012, ApJS, 202, 7

Yang, H., Johns-Krull, C. M., & Valenti, J. A. 2005, ApJ, 635, 466