## Two young eclipsing binaries in Orion with temperatures and radii affected by spots and third bodies

Marina Kounkel 1 and Keivan G. Stassun 2

1 Department of Physics and Astronomy, University of North Florida, 1 UNF Dr, Jacksonville, FL, 32224, USA 2 Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA

## ABSTRACT

In this work we present a model of two young eclipsing binaries in the Orion Complex. Both heavily spotted, they present radii and temperatures that are in disagreement with the predictions of standard stellar models. 2M05 -06 consists of two stars with different masses ( ∼ 0.52 and ∼ 0.42 M ⊙ ) but with very similar radii ( ∼ 0.9 R ⊙ ), and with the less massive star having a highly spotted surface that causes it to have a hotter (unspotted) photosphere than the higher-mass star. The other system, 2M05 -00, consists of two stars of very similar masses ( ∼ 0.34 M ), but very different radii ⊙ ( ∼ 0.7 and ∼ 1.0 R ⊙ ), which creates an appearance of the two eclipsing stars being non-coeval. 2M05 -00 appears to have a tertiary companion that could offer an explanation for the unusual properties of the eclipsing stars, as has been seen in some other young triple systems. Comparing the empirically measured properties of these eclipsing binaries to the predictions of stellar models, both standard and magnetic, we find that only the magnetic models correctly predict the observed relationship between mass and effective temperature. However, standard (non-magnetic) models better predict the temperatures of the unspotted photospheres. These observations represent an important step in improving our understanding of pre-main-sequence stellar evolution and the roles of spots and tertiaries on fundamental stellar properties.

## 1. INTRODUCTION

Mass and radius are some of the most fundamental properties of a star. In pre-main sequence (PMS) stars in particular, measuring these parameters is vital in order to understand their evolution. Obtaining accurate, unbiased, and model-independent measurements, however, is difficult to do for most young stars. Typically, these parameters are interpolated from using a combination of photometric or spectroscopic properties relative to the theoretical isochrones. But, these models of isochrones do not perfectly agree with each other and they often offer an imperfect fit to the data, e.g., due to prominent spots on the surface of PMS stars.

aration, orbital period, inclination, and eccentricity, all directly from the data with no model assumptions. Unfortunately, at the current time the census of low mass PMS ESB2s is still quite small, amounting to only ∼ 30 systems, in large part to the great volume of photometric time series data that is needed to observe eclipses, combined with high resolution multi-epoch spectroscopy. A larger, more statistical census of young ESB2s is needed for a more comprehensive analysis.

Thus, a census of young stars with independently derived stellar properties vital to improve the calibration of these models. One of the most direct methods of measuring them is though analyzing the orbits of eclipsing double lined spectroscopic binaries (ESB2s): in them it is possible to solve for masses of the individual stars, their radii, their surface temperatures, the orbital sep-

In this paper we present a characterization of two such system, 2MASS J05351685 -0618158 (hereafter 2M05 -06) and 2MASS J05454167 -0004024 (hereafter 2M05 -00), both associated with the Orion Complex. In Section 2 we describe the light curves and radial velocities used in this work. In Section 3 we describe the fitting procedure to determine the masses, radii, and other stellar and orbital parameters. In Section 4 we present the parameters of both of these systems. In Section 5 we discuss the implications of these parameters on the current evolutionary models, and present comparison with other young EBs. Finally, in Section 6 we conclude our resuts.

Figure 1. TESS light curves for 2M05 -06 (top) and 2M05 -00 (bottom)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-28a979c9-d35f-4222-8f69-349a526d794b.jpg)

APOGEE is a high resolution spectrograph with R ∼ 22 500 covering the wavelength range of 1.5-1.7 , µ mthat is capable of observing up to 300 sources simultaneously in a 2.5 ◦ field of view. It has performed a survey of the Orion Complex, obtaining multi-epoch spectra for ∼ 3500 young stars in this star forming region (Kounkel et al. 2018). In these data, 2M05 -06 and 2M05 -00 have previously been identified as SB2s in (Kounkel et al. 2021), with the full reported complement of the individual RVs for both primary and the secondary in both systems. A total of 7 and 6 epochs are available for both of them respectively; they can be resolved as SB2s in all of them. Through the course of the analysis we had to reassess which RVs were assigned to the primary and which to the secondary to fix mismatched assignments in several epochs in Kounkel et al. (2021) for 2M05 -06.

Through spectral modeling of the APOGEE spectra, 2M05 -06 has been reported to have T eff =3300 ± 25 K and log g =4.17 ± 0.07. For 2M05 -00, T eff =3400 ± 30 K and log g =4.14 ± 0.03 (Sprague et al. 2022); this is in great agreement with their spectral types derived from lower resolution optical spectra of M4 (Fang et al. 2009) and M3 (Hsu et al. 2012) respectively (which convert to 3200 K and 3400 K Pecaut &amp; Mamajek 2013).

Upon examining light curves from TESS Sector 6 full frame images constructed with eleanor (Feinstein et al. 2019), these two systems showed the presence of eclipses, with the respective period of 1.7044 and 1.817138 days(Figure 1). Time series data is also avail- able for them as part of Gaia DR3 (Gaia Collaboration et al. 2022); while they are not labelled as eclipsing in the data release, through being folded onto a correct period, eclipses in all 3 bandpasses (G, BP, RP) are also apparent (Figures 2,3).

TESS light curves for both stars exhibit significant asymmetry in the continuum, indicative of the rotational variability with the timescale comparable to the orbital period. Despite their youth, rotation of both binaries appear to be tidally locked to their respective orbit. In Gaia passband, the effect is less apparent due to the sparcity of the data making it difficult to conclusively observe the shape of the continuum.

## 3. ANALYSIS

## 3.1. Basic Setup

We used PHOEBE (Conroy et al. 2020) to perform the analysis of the light &amp; RV curves for each system. There are a total of 10 orbitals parameters we considered for each system. They include individual radii, R 1 , R 2 , their temperatures T eff , 1 T eff , 2 , orbital period P , time of passage of the superior conjunction t 0 , semi-major axis a , inclination i , mass ratio q , RV of the system v 0 . Eccentricity e and argument of periastron ω were initially considered, but were deemed unnecessary as the systems appear to have circular orbits. We also fit third light contamination l 3 for each of the four light curves: in particular, since TESS has low spatial resolution of 21'/pixel, care needs to be taken by accounting for con-

tamination in denser regions such as star clusters or associations. We also include passband luminosity for all 4 light curves, this is a nuisance parameter corresponding to the scaling of the continuum flux, which is recommended for stability within PHOEBE.

Finally, since these systems do appear to exhibit rotational variability, in both systems we also fit spot parameters on each star: including spot longitude ( ) and l latitude ( b ), and spot radius ( r ). During the fitting the spot temperature contrast was fixed to 0.85, which is typical for spots in stars of these T eff (Berdyugina 2005).

In what follows, we discuss the treatment of T eff and of spots in greater detail, and then summarize the final fitting approach.

## 3.2. T eff from Spectral Energy Distribution Fitting

Eclipsing binaries make it possible to precisely determine the ratio of T eff of the primary and the secondary through considering the differences between the eclipse depths, even when all that is available is just a light curve at a single bandpass (e.g., Prˇa et al. 2008). s However, while it is also possible to determine 'true' T eff of each individual star given sufficient data about the system (such as RVs, as well as light curves at different bandpasses), they are much more difficult to constrain than T eff ratio. In practice, absolute calibration of T eff often requires an independent fit of the spectral energy distribution (SED) to place stronger constraints on the limits of plausible T eff s for the sytem (e.g. Kounkel et al. 2024), as without it, the model may significantly underestimate or overestimate T eff s (while preserving T eff ratio). Unfortunately, however, PHOEBE does enable consideration of single epoch broadband photometry when modeling the system, requiring an independent SED fitting procedure, which may then require iteratively recomputing a model within PHOEBE with a stronger prior in order to converge.

We thus use SEDFit (Kounkel 2023). It downloads all of the available photometry from UV to mid-IR and considers a model photosphere using BT-Settl synthetic grid (Allard et al. 2011). These models are interpolated to T eff and log g of a star, scaled to the appropriate radius, fitting for the distance and extinction. In case of multiple stars, different templates are added together.

## 3.3. Spots

The situation is made more complex when spots are considered. PHOEBE does not report a 'true' effective temperature in spotted stars, rather what it refers to as T eff is the T eff of the unspotted portion of the photosphere. But, since PHOEBE can model a spot independently, with the spot being of a given radius and given temperature contrast, we recalculate T eff in Table 1 with spots in mind across the entire surface of a star as ( ∫ T 4 dA/A ) 0 25 . as per definition. And, for the purposes of SED fitting, to avoid biases, we compute a mesh for each star, each surface in the mesh having an independent estimate of temperature, log g , and surface area; a synthetic model is computed for each surface, and the flux of each of them is added together for a given star.

We note that whether a spot is being fitted or not, the stellar parameters that are determined by PHOEBE are typically not affected; it generally produces an almost identical fit for masses, radii, temperature ratios, third light, and all of the other orbital parameters (Kounkel et al. 2024), which we were able to confirm in the system in this work. Thus, through iterative exploration of the parameter space, considering spots of different sizes alongside with different PHOEBE models and the quality of fit of the resulting SED, we place the prior on the temperature of the photosphere to be 3000-3800 K for 2M05 -06, and 3300-3600 for 2M05 -00.

We further note that the modeling of the spots is somewhat simplistic. We assume that there is only one large spot on the photosphere, as opposed to several small ones, and we assume that both stars are spotted (vs, e.g., one spotted and one unspotted). While the precise configuration of the modeled spots does not play a significant role in the derived spot sizes, contrast ratios, or all of the stellar parameters (Miller et al. 2021), it is nonetheless one of the limitations of the model.

## 3.4. Final Fitting

We utilized emcee module within PHOEBE to solve the systems. In initial set up, we adopt T eff suggested by the APOGEE spectra (Sprague et al. 2022), q and v 0 for the systems derived in Kounkel et al. (2021), and the P as measured by Lomb Scargle periodogram. P required a very narrow initial distribution, as well as a prior, to ensure emcee would accurately converge, and not be pulled, e.g., towards the rotational period of a system. Other parameters were set up with reasonable distributions. The initial distribution guesses are included in Table 1, and PHOEBE is able to improve on these distributions, including going beyond the bounds (e.g. for uniform distribution). An MCMC chain was set up with 250 walkers, and 25,000 steps; full convergence is achieved after ∼ 10,000 steps.

The derived parameters are listed in Table 1. The resulting fit for the light curve and the RV curve are shown in Figures 2 and 3.

## 4. RESULTS

Table 1. Stellar parameters of ESB2s

| Parameter                                                                                             | 2M05 - 06                                              | 2M05 - 06                                                                                                                                                                                           | 2M05 - 00                                               | 2M05 - 00                                                                                                                                                                                                                     |
|-------------------------------------------------------------------------------------------------------|--------------------------------------------------------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|                                                                                                       | Initial                                                |                                                                                                                                                                                                     |                                                         |                                                                                                                                                                                                                               |
| R 1 (R ⊙ ) R 2 (R ⊙ ) M 1 (M ⊙ ) M 2 (M ⊙ ) T phot , 1 T phot , 2 T eff , 1 (K) T eff , 2 (K) log g 1 | 0.6-1.1 0.6-1.1                                        | Fitted 0 . 902 +0 . 029 - 0 . 023 0 . 891 +0 . 022 - 0 . 024 0 . 424 +0 . 023 - 0 . 023 0 . 522 +0 . 029 - 0 . 028 3753 . 0 +47 . 0 - 331 . 0 3630 . 0 +54 . 0 - 316 . 0 3483 . 0 +46 . 0 - 300 . 0 | Initial 0.6-1.1 0.6-1.1 3400 ± 10 3400 ± 10             | Fitted 1 . 011 +0 . 117 - 0 . 093 0 . 725 +0 . 112 - 0 . 077 0 . 348 +0 . 018 - 0 . 021 0 . 339 +0 . 017 - 0 . 019                                                                                                            |
|                                                                                                       |                                                        |                                                                                                                                                                                                     |                                                         | +100 . 0                                                                                                                                                                                                                      |
| (K) (K)                                                                                               | 3400-3500                                              |                                                                                                                                                                                                     |                                                         | 3500 . 0 - 200 . 0 +100 . 0                                                                                                                                                                                                   |
| log g 2 P (d) a (R ⊙ ) i ( ◦ )                                                                        | 3200-3600 1.70439 ± 1e-6 5-7 85-90                     | 3533 . 0 +51 . 0 - 300 . 0 4 . 155 +0 . 026 - 0 . 029 4 . 256 +0 . 025 - 0 . 021 1 . 704404 +4 . 2 e - 07 - 4 . 2 e - 07 5 . 894 +0 . 087 - 0 . 09 88 . 64 +0 . 36 - 0 . 28 +0 . 073                | 1.81719484 ± 1e-6 5-7 85-90                             | 3340 . 0 - 190 . 0 3313 . 0 +153 . 0 - 226 . 0 3124 . 0 +114 . 0 - 194 . 0 3 . 97 +0 . 078 - 0 . 091 4 . 25 +0 . 1 - 0 . 11 1 . 817212 +1 . 8 e - 06 - 1 . 7 e - 06 5 . 532 +0 . 085 - 0 . 103 82 . 3 +1 . 6 - 2 . 2 +0 . 041 |
| q t 0 (MJD) γ (km s - 1 ) r 1 ( ◦ ) l 1 ( ◦ )                                                         | 1.2 ± 0.05 58986.407877 ± 0.0001 20.9 ± 1.5 0-90 0-360 | 1 . 231 - 0 . 067 58986 . 40791 +0 . 00013 - 0 . 00016 21 . 2 +2 . 2 - 2 . 0 89 . 75 +0 . 25 - 2 . 12 255 . 8 +7 . 8 - 8 . 6                                                                        | 0.971 ± 0.05 58986.6671513 ± 0.0001 23.7 ± 1 0-90 0-360 | 0 . 974 - 0 . 035 58986 . 67893 +0 . 00092 - 0 . 00128 27 . 4 +1 . 2 - 1 . 2 75 . 0 +14 . 0 - 17 . 0 0 . 4 +198 . 2 - 7 . 1                                                                                                   |
| b 1 ( ◦ ) r 2 ( ◦ ) ◦                                                                                 | -90-90 0-90                                            | 52 . 1 +6 . 4 - 9 . 3 52 . 1 +6 . 4 - 9 . 3 +8 . 7                                                                                                                                                  | -90-90 0-90                                             | 110 . 0 +11 . 0 - 43 . 0 83 . 7 +6 . 2 - 12 . 2 4 . 9 +45 . 8 - 6 . 0                                                                                                                                                         |
| l 2 ( ) b 2 ( ◦ ) l 3 T                                                                               | 0-360 -90-90 0-0.5                                     | 288 . 5 - 8 . 0 286 . 0 +2 . 4 - 2 . 1 0 . 3212 +0 . 0057 - 0 . 0058 +0 . 00376                                                                                                                     | 0-360 -90-90 0-0.8 0-0.8                                | 160 . 0 +28 . 0 - 44 . 0 0 . 31 +0 . 14 - 0 . 22 0 . 25 +0 . 15                                                                                                                                                               |
| l 3 G l 3 BP                                                                                          | 0-0.05 0-0.05                                          | 0 . 00055 - 0 . 00054 0 . 005 +0 . 0196 - 0 . 005                                                                                                                                                   | 0-0.8                                                   | - 0 . 24 0 . 24 +0 . 19 - 0 . 24                                                                                                                                                                                              |
|                                                                                                       | 0-0.05                                                 |                                                                                                                                                                                                     |                                                         | +0 . 15                                                                                                                                                                                                                       |
| l 3 RP                                                                                                |                                                        | 0 . 016 +0 . 014 - 0 . 014                                                                                                                                                                          |                                                         | 0 . 24 - 0 . 24                                                                                                                                                                                                               |
|                                                                                                       |                                                        |                                                                                                                                                                                                     | 0-0.8                                                   |                                                                                                                                                                                                                               |
| A V (mag)                                                                                             |                                                        | 0.337                                                                                                                                                                                               |                                                         | 0.67                                                                                                                                                                                                                          |
| Dist. (pc)                                                                                            |                                                        |                                                                                                                                                                                                     |                                                         |                                                                                                                                                                                                                               |
|                                                                                                       |                                                        |                                                                                                                                                                                                     |                                                         | 422                                                                                                                                                                                                                           |
|                                                                                                       |                                                        | 334                                                                                                                                                                                                 |                                                         |                                                                                                                                                                                                                               |

## 4.1. 2M05 -06

Overall, the model for 2M05 -06 is able to closely reproduce the observed SED of the system in the optical regime, although in WISE passbands there is a significant infrared excess due to the presence of a protoplanetary disk, although this system does not show evidence of accretion (Hsu et al. 2012; Fang et al. 2013).

as spot coverage of &gt; 50% is not uncommon in pre-main sequence stars (e.g., Gully-Santiago et al. 2017). The spot on the more massive star is smaller. Thus, calculating the true T eff does result in the more massive star in the system being hotter than the less massive star.

The system consists of two stars with comparable luminosity and very similar radii of ∼ 0.9 R ⊙ . However, the masses of them are quite different, with the primary having a mass of ∼ 0 52 M . ⊙ , and the secondary having a mass of ∼ 0 42 . M . ⊙ Interestingly, the temperature of the photosphere in these two stars is backwards in comparison to the expectations, with the more massive star having cooler photosphere. This is explained by the less massive star having a more sizeable spot on its photosphere. PHOEBE estimates that the size of the spot on the secondary is ∼ 90 , i.e, covering 50% of the star ◦ - this is the maximum spot size that is supported by PHOEBE, and the actual spot coverage may be larger,

Such a difference in spot sizes between the stars of different masses is consistent with pre-main sequence stars (in the T eff range of up to 4500 K) having stronger magnetic field strengths with increased mass, and thus having a larger discrepancy between temperatures of spotted and unspotted photospheres (Flores et al. 2022).

## 4.2. 2M05 -00

2M05 -00 consists of two stars with near-equal masses of ∼ 0 34 . M . ⊙ However, they do have very different radii, one with the size of ∼ 1 0 R . ⊙ , the other 0 7 R . ⊙ . The larger one appears to be hotter, and this is true both for the temperature of the photosphere as well as the true T eff , as both stars have comparable spot sizes.

Gaia DR3 parallax of 2M05 -00 is 2.3746 ± 0.0493 mas, translating to a likely distance range of 412.6-430.1 pc.

Figure 2. Model fit of the observational data towards 2M05 -06. Left panel shows the fit of the available lightcurves (in T, G, GP, and RP bandpasses) as well as the residuals (normalized by the typical uncertainties in the fluxes to meaningfully be displayed on the same plot). The blue line shows a flat continuum to demonstrate the significants of spots onto the light curve. Top right panel shows the fit of the radial velocities. Middle right panel shows the model SED (with the SED of the individual stars being shown in yellow, the combined SED in read, the black lines showing the measurements with the uncertainties and the bandpass width, and blue dots showing interpolated model fluxes onto these bandpasses). In this system, the last three fluxes, corresponding to W1, W2, and W3 passbands appear to show IR excess. We also note that the two stars have extremely similar fluxes integrated across the entire surface, thus the two yellow lines are plotted on top of one another. Finally, the bottom right panel shows a mesh grid depicting the surface temperature, size, and the separation of the two stars.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-04d67b61-81ff-4eef-bca2-aa29b7899eda.jpg)

SED fitting using the derived stellar parameters at this distance results in a significant underestimation of the total observed flux towards this system across most of the electromagnetic spectrum.

A better fit can be achieved if the distance is adjusted to 360 pc, however, PHOEBE fit suggests a third light of ∼ 0.31 for TESS passband, and ∼ 0 25 for Gaia pass-. pands. Since Gaia does have very high spatial resolution, it is highly improbable for the third light contamination to originate from a chance alignment of an unrelated star along the same line of sight. Rather, it is likely that 2M05 -00 is a triple system, with a tertiary having a significantly longer orbital period. Indeed, allowing inclusion of a third object while doing an SED fit (at the distance reported by Gaia) favors a presence of a star with T eff ∼ 3130 K and R ∼ 0.92 R ⊙ . This is a similar T eff as a more compact star in the system, and the estimated radius is in between that of either star in the inner binary. Such a star typically has luminosity which is 20-30% that of the combined luminosity of the inner binary, consistent with the estimate for the l 3 that is seen across all passbands.

## 4.3. Degeneracy in the radii

In eclipsing binaries with circular orbits, the radius parameter that is most stringently constrained is the sum of the radii. Although extreme combinations of the radii can be ruled out based on the bottom of the light curve during the eclipse, it is possible to redistribute the radii of the individual stars within a certain range of the sum and still produce a comparable light curve fit. Thus, external confirmation is often useful.

Examining the cross-correlation function of the spectrum consisting of two stars of comparable temperature can be used for such a purpose. The ratio of area under the peak of the CCF is roughly proportional to the ratio of luminosities at the spectral bandpass. We note that this is an approximation that does have its limitations -e.g. in stars of very different temperatures affecting the strength of the cross-correlation depending on the template that is being used (Kounkel et al. 2024).

In a heavily spotted system, such as the ones in this work, this approximation has further complications. The ratio of the luminosities of two stars can vary depending on the phase, depending on the side of a given star that is facing us. Indeed, the ratio of the area under the CCF peaks does vary at different epochs. Unfortunately, while we are able to create a model of spots based on the lightcurves, this model is fairly simplistic, and given that the morphology and the position of spots can vary significantly in time, even if the simplistic model was completely accurate, it may not necessarily hold at the time the spectra was taken, as it was several years prior to the TESS observations. Nonetheless, from the rough estimate of the ratio of the luminosities, we can approximate the ratio of the radii.

The typical ratio of the area under the CCF peaks between the secondary and the primary for is 1.2-1.4 for 2M05 -06, and 0.5-0.7 for 2M05 -00. Using the respective T eff for all of the stars, this results in the estimate of the ratio of the radii of 1.0-1.1 for 2M05 -06, and 0.8-0.9 for 2M05 -00. Combining these ratios with the sum of the radii derived by PHOEBE produces the radii of 0.87+0.91 R ⊙ for 2M05 -06, and 0.91+0.82 R ⊙ for 2M05 -00, which is consistent with the reported uncertainties in Table 1.

## 5. DISCUSSION

## 5.1. Ages

Most commonly used method to determine ages of premain sequence stars is through interpolating their photometry over the isochrones, but ages of the binary stars are commonly underestimated since they have higher total luminosity in comparison to single stars. However, since we have determined all of the fundamental parameters for these stars, we can use them to estimate ages of the individual stars in the system instead. We consider two projections: mass vs radius, as well as T eff (and T phot ) vs radius (Figure 4. We compare our observations to MIST (Choi et al. 2016) isochrones; although we considered other models, they do not have a significant impact on the results.

2M05 -00 appears to be a system of two stars that appear to be mostly coeval, with the age of 6.9-7.1 dex, or 8-12 Myr depending on the precise projection that is used for evaluation. This is broadly consistent with the system being in the foreground of Orion A molecular cloud and being a member of much more evolved Orion D (Kounkel et al. 2018).

However, 2M05 -06 appears to be somewhat peculiar. In this system, the primary favors the age of 6.7 dex ( ∼ 5 Myr), while the secondary favors the age of 7 dex ( ∼ 10 Myr). Such a large apparent difference in age is highly unusual, as typically young EBs do tend to show coevality to within ∼ 20% (Stassun et al. 2014), and indeed this is the expectation from standard binary star formation theory.

Given that 2M05 -06 is a triple star system, it is not impossible to consider that the tertiary companion may have been sufficiently close as to the inner binary so that its orbital motion would have affected the observed ratio of the radial velocities over the 245 day baseline over which RVs have been obtained. This could have

Figure 3. Same as figure 2, but for 2M05 -00. In the middle right panel, SED marked with 1 corresponds to the primary in Table 1, and 2 to the secondary. The green line is the fitted SED of the third star in the system that is needed in order to match the observed flux while keeping the system fixed at the reported distance from Gaia.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-d6ef83e7-ece8-4032-b79b-83560b342bed.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-a7ac9726-d0e7-474e-af3a-102ed69800b8.jpg)

Temperature (K)

3000

3100

3200

3300

3400

3500

3600

3700

3800

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-37b81ad0-7c89-4486-b6ce-35bda8b71f2e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-2da109d1-b0a6-4d59-87ed-b4d8e19aa629.jpg)

Temperature (K)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-8d0529eb-c7e3-4ee3-a171-22e5c0ada713.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-e68b8e70-765a-422e-9b26-e3a9c4e986bc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00886v1-figure-bef75ee6-1d79-49bc-a425-04b9dfd3ced6.jpg)

2M05-06a

2M05-06b

2M05-06a

2M05-06b

2M05-00a

2M05-00b

2M05-00a

2M05-00b

Figure 4. Top: Effective temperature (dashed line) and the temperature of the photosphere (solid line) for 2M05 -00 (red) and 2M05 -06 (blue) versus their radii. The transformation of this parameter space onto mass and age using isochrones is also shown, alongside the masses of the stars in these two systems as a function of radius. Middle: Mass vs radius relative to the isochrones, showing the typical uncertainty in the measurements. Bottom: the difference between the observed masses (from the orbital fitting) and the interpolated model mass from temperatures, divided by the observed mass. Error bars include the uncertainty in the temperature. Left: non-magnetic MIST isochrones (Choi et al. 2016). Right: magnetic SPOT isochrones (Somers et al. 2020), with 34% spot filling fraction.

had an impact on the observed mass ratio, making the determination of the masses to be incorrect.

This, however, would not propagate to T eff s, as T eff ratio is determined purely from the light curves (with the absolute calibration assisted by the SED). Similarly, while the absolute calibration of the individual radii does depend on the orbital speed, the ratio of the radii is driven solely by the eclipse geometry. The apparent factor of 2 difference between the ages of the stars persists in T eff vs radius projection.

2M05 -06 is found in the vicinity of NGC 2068 cluster in Orion B molecular cloud which has a typical age of 1-2 Myr (Kounkel et al. 2018), and although in Orion as a whole there are a number of older stars, there are very few of them near NGC 2068. Thus, of the two derived ages, younger age is significantly favored for this system, and a chance dynamical capture of an isolated older star is highly unlikely. It is unclear what type of activity could be responsible for such a large disparity.

## 5.2. T eff

Surface temperature is strongly correlated with mass. As pre-main-sequence stars evolve down the Hayashi track, there is almost no dependence on age between mass and temperature, resulting in a robust relation.

Spots, however, introduce a significant complication. Masses derived dynamically, in comparison to the masses derived from the temperature, can have 50-100% discrepancy (e.g., David et al. 2019; Flores et al. 2022). We observe this here also. T eff s that have been derived for the systems in these work are generally consistent with spectroscopic estimates (Fang et al. 2009; Hsu et al. 2012; Sprague et al. 2022), pointing to their accuracy. These T eff s significantly underestimate the masses we derive from the full orbital fit by 20-50% when using non-magnetic isochrones such as MIST (Choi et al. 2016), and the difference is especially pronounced in the case of 2M05 -06. Masses derived T phot are a much closer fit to the dynamical masses. On the other hand, comparison made using isochrones that do consider magnetic activity, such as SPOTS tracks (Somers et al. 2020) with 34% spot coverage do produce a much better agreement between the dynamical masses and the observed T eff s. This set of models produced the best agreement between these two parameters; we also examined isochrones with 17% and 51%, with the former underestimating the mass, and the latter overestimating the mass.

## 5.3. Comparison to other systems

While the eclipsing binaries shown in this paper do present some peculiarities, they are not unique. For ex- ample, 2M05 -00 has a remarkable similarity to Parenago 1802, a ∼ 1 Myr system in Orion consisting of two stars of identical masses, but with radii discrepancy of ∼ 7%, T eff discrepancy of ∼ 9%, which results in a luminosity discrepancy of ∼ 60%, which also has an appearance of not being coeval (Stassun et al. 2008). Similarly to 2M05 -00, it appears to be a triple system, and it has been interpreted that the tidal interactions could be responsible for the difference in the apparent luminosity of the binary (G´mez Maqueo Chew et al. 2012). o

On the other hand, 2M05 -06 can be compared to 2MASS J05352184-0546085, a system of eclipsing brown dwarfs in Orion, where the lower mass star in the system appears to be hotter than the primary (Stassun et al. 2007), due to the primary having stronger chromospheric activity than the secondary, suppressing its temperature and inflating its radius (Stassun et al. 2012).

## 6. CONCLUSIONS

In this work we present the analysis of the two young eclipsing binaries 2M05 -06 and 2M05 -00. With ages &lt; 10 Myr, they are valuable benchmarks of testing stellar evolutionary models, as the number of such systems still counts in a few dozen (e.g. Stassun et al. 2014; David et al. 2019). These two systems, however, present interesting properties that differ from the predictions of standard stellar models and that therefore provide insight into the physical mechanisms that can alter the fundamental stellar properties of young stars.

One of them, 2M05 -06, appears to be a mostly coeval system with two stars of different mass ( ∼ 0.52 and ∼ 0.42 M ⊙ ), but very similar radius ( ∼ 0.9 R ⊙ ). With the less massive star being very heavily spotted, it actually has a hotter temperature of the (unspotted) photosphere in comparison to the more massive star.

The other, 2M05 -00 (a likely hierarchical triple system), has the inner binary with very similar masses( ∼ 0.34 M ⊙ ), but very different radii ( ∼ 0.7 and ∼ 1.0 R ⊙ ). This creates the appearance that these stars are not coeval, one with the age of ∼ 4 Myr and the other with the age of ∼ 10 Myr according to standard stellar isochrones, with the population in which they are found favoring the younger age. Similar behavior has been reported previously for at least one other young EB (Par 1802; see Stassun et al. 2008; G´mez Maqueo o Chew et al. 2012), which similarly possesses a tertiary whose dynamical heating has been suggested to possibly explain the anomalous properties of the inner binary (G´ omez Maqueo Chew et al. 2019).

We also find that spots play a significant role in affecting the overall temperature of a star, and that effective temperature appears to significantly underesti-

mate dynamical masses in these EBs when considering non-magnetic isochrones; magnetic tracks offer a much better agreement. Alternatively, the temperature of the unspotted photosphere appears to be in agreement with the non-magnetic isochrones.

As directly measured masses of more pre-mainsequence stars are measured, the improved temperature determination of both spots and photospheres (such as in EBs as in this work, or through spectral decomposition, e.g., Gully-Santiago et al. 2017) will permit a fuller assessment of the role of magnetism in altering the fundamental physical properties of young stars.

MK acknowledges support provided by NASA grant 80NSSC24K0620.

This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www. cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www. cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

Software: PHOEBE (Conroy et al. 2020)

## ACKNOWLEDGMENTS

## REFERENCES

| Allard, F., Homeier, D., & Freytag, B. 2011, in                                                                    |
|--------------------------------------------------------------------------------------------------------------------|
| Astronomical Society of the Pacific Conference Series,                                                             |
| Vol. 448, 16th Cambridge Workshop on Cool Stars,                                                                   |
| Stellar Systems, and the Sun, ed. C. Johns-Krull, M. K.                                                            |
| Browning, & A. A. West, 91,                                                                                        |
| Berdyugina, S. V. 2005, Living Reviews in Solar Physics, 2, 8, doi: 10.12942/lrsp-2005-8                           |
| Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102, doi: 10.3847/0004-637X/823/2/102                     |
| Conroy, K. E., Kochoska, A., Hey, D., et al. 2020, ApJS, 250, 34, doi: 10.3847/1538-4365/abb4e2                    |
| David, T. J., Hillenbrand, L. A., Gillen, E., et al. 2019, ApJ, 872, 161, doi: 10.3847/1538-4357/aafe09            |
| Fang, M., Kim, J. S., van Boekel, R., et al. 2013, ApJS, 207, 5, doi: 10.1088/0067-0049/207/1/5                    |
| Fang, M., van Boekel, R., Wang, W., et al. 2009, A&A, 504, 461, doi: 10.1051/0004-6361/200912468                   |
| Feinstein, A. D., Montet, B. T., Foreman-Mackey, D., et al. 2019, PASP, 131, 094502, doi: 10.1088/1538-3873/ab291c |
| 2022, ApJ, 925, 21, doi: 10.3847/1538-4357/ac37bd Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al.       |
| 2022, arXiv e-prints, arXiv:2208.00211. https://arxiv.org/abs/2208.00211                                           |
| G´mez o Maqueo Chew, Y., Stassun, K. G., Prˇa, s A., et al. 2012, ApJ, 745, 58, doi: 10.1088/0004-637X/745/1/58    |
| G´mez o Maqueo Chew, Y., Hebb, L., Stempels, H. C., et al. 2019, A&A, 623, A23, doi: 10.1051/0004-6361/201833299   |
| Gully-Santiago, M. A., Herczeg, G. J., Czekala, I., et al. 2017, ApJ, 836, 200, doi: 10.3847/1538-4357/836/2/200   |

| Hsu, W.-H., Hartmann, L., Allen, L., et al. 2012, ApJ, 752, 59, doi: 10.1088/0004-637X/752/1/59                              |
|------------------------------------------------------------------------------------------------------------------------------|
| Kounkel, M. 2023, mkounkel/SEDFit: 0.3, 0.3, Zenodo, doi: 10.5281/zenodo.8076501                                             |
| Kounkel, M., Statti, M., Kulkarni, A., Stassun, K. G., & Sun, M. 2024, MNRAS, 527, 3806, doi: 10.1093/mnras/stad3439         |
| Kounkel, M., Covey, K., Su´rez, a G., et al. 2018, AJ, 156, 84, doi: 10.3847/1538-3881/aad1f1                                |
| Kounkel, M., Covey, K. R., Stassun, K. G., et al. 2021, AJ, 162, 184, doi: 10.3847/1538-3881/ac1798                          |
| Miller, A., Kounkel, M., Boggio, C., Covey, K., & Price-Whelan, A. M. 2021, PASP, 133, 044201, doi: 10.1088/1538-3873/abeaf7 |
| Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9, doi: 10.1088/0067-0049/208/1/9                                           |
| Prˇa, s A., Guinan, E. F., Devinney, E. J., et al. 2008, ApJ, 687, 542, doi: 10.1086/591783                                  |
| Somers, G., Cao, L., & Pinsonneault, M. H. 2020, ApJ, 891, 29, doi: 10.3847/1538-4357/ab722e                                 |
| Sprague, D., Culhane, C., Kounkel, M., et al. 2022, AJ, 163, 152, doi: 10.3847/1538-3881/ac4de7                              |
| Stassun, K. G., Feiden, G. A., & Torres, G. 2014, NewAR, 60, 1, doi: 10.1016/j.newar.2014.06.001                             |
| Stassun, K. G., Kratter, K. M., Scholz, A., & Dupuy, T. J. 2012, ApJ, 756, 47, doi: 10.1088/0004-637X/756/1/47               |
| Nature, 453, 1079, doi: 10.1038/nature07069                                                                                  |
| Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2007, ApJ, 664, 1154, doi: 10.1086/519231                                   |
Marina Kounkel Keivan G. Stassun
Two young eclipsing binaries in Orion with temperatures and radii affected by spots and third bodies
In this work we present a model of two young eclipsing binaries in the Orion
Complex. Both heavily spotted, they present radii and temperatures that are in
disagreement with the predictions of standard stellar models. 2M05-06 consists
of two stars with different masses (~0.52 and ~0.42 Msun) but with very similar
radii (~0.9 Rsun), and with the less massive star having a highly spotted
surface that causes it to have a hotter (unspotted) photosphere than the
higher-mass star. The other system, 2M05-00, consists of two stars of very
similar masses (~0.34 Msun), but very different radii (~0.7 and ~1.0 Rsun),
which creates an appearance of the two eclipsing stars being non-coeval.
2M05-00 appears to have a tertiary companion that could offer an explanation
for the unusual properties of the eclipsing stars, as has been seen in some
other young triple systems. Comparing the empirically measured properties of
these eclipsing binaries to the predictions of stellar models, both standard
and magnetic, we find that only the magnetic models correctly predict the
observed relationship between mass and effective temperature. However, standard
(non-magnetic) models better predict the temperatures of the unspotted
photospheres. These observations represent an important step in improving our
understanding of pre-main-sequence stellar evolution and the roles of spots and
tertiaries on fundamental stellar properties.