The large-scale magnetic field of the eccentric pre-main-sequence binary system V1878 Ori

We report time-resolved, high-resolution optical spectropolarimetric observations of the young double-lined spectroscopic binary V1878 Ori. Our observations were collected with the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope through the BinaMIcS large programme. V1878 Ori A and B are partially convective intermediate mass weak-line T Tauri stars on an eccentric and asynchronous orbit. We also acquired X-ray observations at periastron and outside periastron. Using the least-squares deconvolution technique (LSD) to combine information from many spectral lines, we clearly detected circular polarization signals in both components throughout the orbit. We refined the orbital solution for the system and obtained disentangled spectra for the primary and secondary components. The disentangled spectra were then employed to determine atmospheric parameters of the two components using spectrum synthesis. Applying our Zeeman Doppler imaging code to composite Stokes $IV$ LSD profiles, we reconstructed brightness maps and the global magnetic field topologies of the two components. We find that V1878 Ori A and B have strikingly different global magnetic field topologies and mean field strengths. The global magnetic field of the primary is predominantly poloidal and non-axisymmetric (with a mean field strength of 180 G). while the secondary has a mostly toroidal and axisymmetric global field (mean strength of 310 G). These findings confirm that stars with very similar parameters can exhibit radically different global magnetic field characteristics. The analysis of the X-ray data shows no sign of enhanced activity at periastron, suggesting the lack of strong magnetospheric interaction at this epoch.


INTRODUCTION
Magnetic fields play very important roles throughout the entire stellar evolution. They directly influence many physical processes, both within stars and in their immediate surroundings. During the early stages of stellar formation and on the pre-main-sequence (PMS), magnetic fields are particularly important, influencing for instance the collapse of molecular clouds (Hennebelle & Inutsuka 2019), the accretion mechanism (Hartmann et al. 2016), the formation and ⋆ E-mail: alexis.lavail@physics.uu.se collimation of jets and outflows (Pudritz & Ray 2019), and stellar spin down (Bouvier 2013). The development of dedicated instrumentation, such as high-resolution spectropolarimeters, and tomographic surface mapping techniques has allowed to measure and map surface magnetic fields in a plethora of cool and hot stars (Donati & Landstreet 2009;Reiners 2012). Yet, we do not have a complete theory that can consistently reproduce the stellar magnetic fields characteristics across the Hertzsprung-Russell diagram (Brun & Browning 2017).
A promising avenue to disentangle the influence of magnetic fields from other effects is the study of binary stars. For such systems, one can assume that the stars were formed in the same conditions, with the same material, and at the same time, which allows to separate the effects of magnetism from the initial conditions. This avenue is being explored by "The Binarity and Magnetic Interactions in various classes of stars" (BinaMIcS) project (Alecian et al. 2015). BinaMIcs aims to study the interplay between binarity and magnetism. This project was allocated two large programmes that ran between 2013 and 2017 with the twin spectropolarimeters ESPaDOnS and Narval, respectively at the Canada-France-Hawaii-Telescope and the Télescope Bernard Lyot at the Pic du Midi observatory. Until now, magnetic fields of premain-sequence double-lined binary systems have been seldom studied through tomographic techniques such as Zeeman Doppler imaging (ZDI). The magnetic field maps of two systems, HD 155555 (M 1,2 ≈ 1M ⊙ ) and V4046 Sgr (M ≈ 0.9M ⊙ ), were published by Dunstone et al. (2008) and Donati et al. (2011), respectively.
V1878 Ori (= Parenago 523 = RX J0530.7-0434 = 1RXS J053043.1-043453) is a weak-line T Tauri star binary system with an eccentric orbit and two nearly equal-mass components. The spectroscopic binary nature of V1878 Ori was first reported by Alcalá et al. (2000). Shortly afterwards, the object was studied by Covino et al. (2001b) using time series of high-resolution optical spectra. They measured variation of the radial velocities (RVs) of each component with the cross-correlation technique. With the help of these measurements, they determined orbital parameters and reported an orbital period of P orb = 40.5738±0.0047 days and a relatively high eccentricity with e = 0.3181±0.002. These authors inferred nearly equal mass and luminosity for the two components, with M A /M B = 1.002 and L A /L B = 1.0, and determined that both components have a K2-K3 spectral type and a projected rational velocity v e sin i = 13 ± 2 km s −1 . Photometric data from Covino et al. (2001a) showed variability, possibly tracing a rotation period of 13.5 days, which amounts to about 1/3 of the system's orbital period. Through a follow-up photometric monitoring, Marilli et al. (2007) reported a photometric rotation period of 12.9 days, continuing to hint towards a non-synchronous rotation.
V1878 Ori was also observed in the radio at millimeter wavelengths (Kóspál et al. 2011) and in X-rays (Getman et al. 2016) to hunt for signs of flares or enhanced emissions at periastron where the magnetospheres could interact and magnetic reconection occur. However, no millimeter wavelength signal was detected from V1878 Ori by Kóspál et al. (2011), and the modelling of the X-ray observations shows only a non-significant increase of the X-ray luminosity at periastron (Getman et al. 2016).
According to the orbital solution by Covino et al. (2001b), confirmed in this work (see Table 2), both components of V1878 Ori have a mass of at least 1.5-1.6 M ⊙ . This means that it is the first binary system of intermediatemass T Tauri stars (IMTTS) that can be studied with ZDI. IMTTS, with masses 1M ⊙ < M < 4M ⊙ , are at a particularly interesting stage of stellar evolution. These stars are believed to be the precursors of PMS Herbig Ae/Be stars and A-and B-type stars on the main-sequence (MS). While all low-mass T Tauri stars (M < 1M ⊙ ) harbour a convective envelope and a priori generate a magnetic field, only 5-10% of Herbig Ae/Be stars and A/B stars show a detectable magnetic field (Alecian et al. 2013; Sikora et al. 2019), which has quite dif-ferent characteristics than for low-mass stars. Studying magnetic fields of IMTTS is therefore a means to understand the origin of magnetism in intermediate-mass stars Villebrun et al. 2019).
Our paper is structured as follows. We present our spectropolarimetric and X-ray observations in Sect. 2. In Sect. 3 we describe our application of the least-squares deconvolution (LSD) technique to the observed data and detection of the Zeeman signatures in Stokes V LSD profiles. We discuss in Sect. 4 our spectral disentangling procedure, yielding a refined orbital solution as well as a disentangled spectra for each component. The disentangled spectra are analyzed to determine stellar parameters as described in Sect. 5. Finally, we present in Sect. 6 simultaneous reconstruction of the global magnetic field maps for the two components using ZDI, the analysis of the X-ray data in Sect. 7, and discuss our results in Sect. 8.

OBSERVATIONS
Our study of V1878 Ori is based on the high-resolution optical spectropolarimetric observations obtained in the context of the BinaMIcS large programme (Alecian et al. 2015) between 2014 and 2016. The data were acquired with the Echelle SpectroPolarimetric Device for the Observation of Stars (ESPaDOnS) instrument mounted at the 3.6-metre Canada-France-Hawaii Telescope (CFHT). ESPaDOnS (Donati 2003;Donati et al. 2006) is a highresolution (R ≡ λ/∆λ ≈ 65000) cross-dispersed spectropolarimeter with a spectral coverage from 370 to 1000 nm that can record circularly and linearly polarized spectra as well as intensity spectra. Each polarized observation consists of a set of four sub-exposures during which the configuration of the polarimeter is changed to swap the position of the two orthogonal polarization beams through the instrument optics and on the detector. This beam-switching technique (Semel et al. 1993) efficiently removes instrumental polarization and effects arising from inhomogeneities between the two optical paths.
We obtained a total of 22 Stokes IV observations taken with exposure times of either 4 × 1715 s (for the observing date 2014-02-21) or 4 × 1355 s (all other observations). The main dataset is comprised of 21 observations obtained over 16 consecutive nights in January 2016. Our data were reduced on the fly at the telescope with the Upena pipeline 1 that runs the LIBRE-ESPRIT routines (Donati et al. 1997). The median S/N of the reduced data per spectral pixel at λ = 566 nm is 128. Information on individual observations can be found in the observing log given in Table 1.
We complemented our study of the magnetic properties of V1878 Ori by analyzing also a long XMM-Newton observation. We observed V1878 Ori with XMM-Newton for a total of ∼ 220 ks. The aim of the observation was to monitor its X-ray emission during and outside periastron. Therefore the observation was divided into two different segments: the first exposure of ∼ 110 ks (ObsID 0763720401) was performed on 2016 March 1 with the system out of periastron; the second exposure of ∼ 110 ks (ObsID 0763720301) was performed on 2016 March 21 near the periastron passage. Neither the spatial nor spectral resolution of XMM-Newton allows to separate the two components. X-ray data were therefore analyzed for the whole binary system.
We analyzed X-ray data gathered with the three EPIC instruments (PN, M1, and M2). EPIC data were first processed using the SAS V15.0 standard tasks. X-ray source events were extracted from a circle with a radius of 37.5 ′′ centered on the target position. A similar nearby region, free from other X-ray sources, was used to extract background events. We considered only events with energy ranging between 0.3 and 7.9 keV. During the two observations some time intervals were affected by significantly high background levels. We therefore discarded these high-background intervals, slightly reducing the actual exposure times.

LEAST-SQUARES DECONVOLUTION
Cool stars generally exhibit very weak Zeeman polarization signals in individual spectral lines. Combining the signal from a large set of spectral lines is a widely used solution in order to obtain a high S/N mean line profile. Here we applied the least squares deconvolution technique (LSD, Donati et al. 1997;Kochukhov et al. 2010) to compute high S/N mean Stokes IV profiles for each observation. LSD effectively performs a weighted co-addition of a set of spectral lines, defined by a line mask, into a mean profile. To construct the line mask, we retrieved a line list from the VALD3 database (Ryabchikova et al. 2015) using a MARCS stellar atmosphere (Gustafsson et al. 2008) with T eff = 4750 K and log g = 4.0. After having removed spectral regions where tellurics or particularly strong and broad spectral lines were present, we obtained a line mask containing 5304 atomic lines deeper than 20% of the continuum. The LSD Stokes IV profiles were then computed with the LSD code described by Kochukhov et al. (2010). The mean parameters used to scale the LSD profiles were the mean wavelength λ 0 = 5643Å and the mean Landé factor z 0 = 1.221. The LSD procedure provided a median S/N gain of 66 with a median S/N of the LSD Stokes V profile of 8263.
The resulting Stokes IV profiles are shown in Fig. 1. The circular polarization signatures are detected unambiguously in both components throughout the part of the orbit sampled by our observations (between orbital phases 0.232 and 0.678). The Stokes I profiles show that the primary component has somewhat deeper spectral lines. As our data mostly samples orbital phases where the two components are not clearly separated, we do not study variation of the mean longitudinal magnetic field for each component.

SPECTRUM DISENTANGLING AND ORBITAL SOLUTION
V1878 Ori is a double-lined spectroscopic binary for which time-dependent superposition of the component spectra greatly complicates interpretation of the observed spectroscopic time series. In order to isolate the contribution from individual stars, we used different spectral disentangling methods. This allowed us to (a) measure accurate radial velocities for each of the two components at each epoch, (b) derive an orbital solution for the binary system, (c) obtain individual disentangled stellar spectra for each of the two components.
In this work, we used the disentangling algorithm described in Folsom et al. (2010) and subsequently used by e.g. Kochukhov et al. (2018); Rosén et al. (2018); Kochukhov & Shulyak (2019). First, we modelled our series of the Stokes I LSD profiles to measure the RV variation of the two components. The algorithm retrieves a mean Stokes I profile for each component as well as a set of RV mesurements. It assumes that, at each epoch, the observed Stokes I LSD profile is a combination of two profiles which stay constant throughout the time series and are shifted by a Covino et al. (2001) This work  Table 1. Log of observations. Columns 1-2 list the UT date and heliocentric Julian date (HJD) of each observation. Our RV measurements from the disentangling of the Stokes I LSD profiles are reported in columns 3-4. The estimated uncertainty of RV measurements is 0.9 km s −1 . Column 5 contains the orbital phase, columns 6-7 respectively list the S/N of the reduced observed spectra per spectral pixel at λ = 566 nm and the S/N of the Stokes V LSD profiles relative to the unpolarized continuum level.   Covino et al. (2001b) and in this work. Covino et al. (2001b) This work We then fitted the RV measurements to derive an orbital solution. We used both our own measurements, which sample around half of the orbit, and the data from Covino et al. (2001b), which, though less precise, has the benefit of covering the entire orbit. We present our orbital solution in Table 2 alongside the previously determined solution by Covino et al. (2001b). The orbital parameters derived from the RV fit include the orbital period P orb , the time of periastron passage T 0 , the eccentricity e, the radial velocity of the mass centre of the system γ, the radial velocity semi-amplitudes K A,B , and the longitude of periastron ω. We also determined the minimum stellar masses M A,B sin 3 i and the projected semi-major axes a A,B sin i. The RV measurements and the best-fitting model are shown in Fig. 2. From the scatter of our RV points around the model, we determined that our measurements have a typical precision of 0.9 km s −1 .
Finally, we applied the disentangling algorithm to the spectra themselves, working iteratively on 10-50 nm-long regions of the spectrum, to eventually cover a large part of the spectrum between 450 and 900 nm. In the disentangling procedure, we kept the previously found RV solution fixed, and recovered the individual spectra for each component, assuming that they exhibit no intrinsic variation. The disentangled spectra S A,B of each components were then rescaled in order to correct for the dilution of the continuum using the formula where r L is the flux ratio of the two components, here r L = 1.0. The spectra were finally individually renormalized by fitting the continuum level with a low-order polynomial. Fig. 3 shows our set of observed spectra as well as the disentangled spectra for the A and B components in a narrow wavelength interval around 601 nm. Table 3. Parameters of V1878 Ori A and B inferred with SME.

STELLAR PARAMETERS
We used the disentangled stellar spectra, corrected for the continuum dilution, to determine stellar parameters for the two components by fitting synthetic spectra to observations. To calculate theoretical stellar spectra and optimise free parameters, we used the Spectroscopy Made Easy (SME) spectrum synthesis code (Piskunov & Valenti 2017). Adopting the solar abundances from Asplund et al. (2009), we downloaded absorption line lists from the VALD database (Ryabchikova et al. 2015) using the Extract Stellar mode with T eff = 4750 K and log g = 4.0. We used 1D model atmospheres from the MARCS grid (Gustafsson et al. 2008). We fitted a set of wavelength intervals sensitive to different stellar parameters, mostly in the wavelength regions 516-519 nm and 600-620 nm, as used in Valenti & Fischer (2005). The same fitting windows around lines of interest were used for both components. We adopted a fixed macroturbulence v mac = 0 km s −1 due to its degeneracy with v e sin i. The instrumental broadening was modelled with a gaussian profile corresponding to a spectral resolution of R = 65000. The free parameters in the fit were the effective temperature T eff , the surface gravity log g, the projected rotational velocity v e sin i, the microturbulence v mic and the overall metallicity [M/H]. The inferred stellar parameters are listed in Table 3. The error bars reported in this table were derived with the analysis of cumulative distribution functions, implemented in recent versions of SME (Piskunov & Valenti 2017). Comparison of the best-fitting synthetic spectra and observations of the A and B components are shown in Fig. 4 for one of the wavelength regions used. Based on the CESAM PMS evolutionary tracks (Morel & Lebreton 2008;Marques et al. 2013) computed in Villebrun et al. (2019), our best-fit stellar parameters correspond to masses of 1.7 ± 0.3M ⊙ and 1.6 ± 0.3M ⊙ for the primary and secondary, respectively and a young age of approximately 0.85 Myr. This suggests that both stars are beginning to develop radiative cores (∼ 10% by mass and ∼ 25% by radius). It is important to note however, that there are large uncertainties in this region of the PMS HR diagram, and that another set of stellar evolutionary tracks such as the Yale-Potsdam Stellar Isochrones (YaPSI; Spada et al. 2017) would place our stars just across the fully-convective limit. Our spectra of V1878 Ori A and B confirm previous observations that these stars have an abnormally high Li abundance (Covino et al. 2001b), suggesting indeed a very young age. This however, is in contrast with the absence of high veiling or other signs of accretion onto the stars which would be expected at ages lower than 1 Myr. It is possible that binarity has impacted the accretion process very early in that case.

ZEEMAN DOPPLER IMAGING
We carried out reconstruction of the large-scale magnetic field geometry of the two components simultaneously using the binary Zeeman Doppler imaging code InversLSDB code described by Rosén et al. (2018), which is based on the single star ZDI code InversLSD code (Kochukhov et al. 2014). InversLSDB has previously been used in the studies of the double-lined spectroscopic binary systems σ 2 CrB (Rosén et al. 2018) and YY Gem (Kochukhov & Shulyak 2019). Only the main dataset acquired in 2016 was used for ZDI inversions. The single observation from 2014 was discarded due to lack of precise rotation periods necessary for phasing together data taken two years apart and a high likelihood that magnetic field topologies have evolved during that time period.
In our analysis, we computed local Stokes IV profiles using the Unno-Rachkovsky solution of the polarised radiative transfer equation, assuming a Milne-Eddington atmosphere.
We also assumed that the LSD profiles behave as a single spectral line, with mean parameters, such as central wavelength and effective Landé factor, adopted from the LSD line mask applied to the observations. The local line profile strength was adjusted to reproduce the Stokes I profiles. Additionally, the local profiles of the primary were scaled by a factor of 1.3 to account for the difference in line depths between the primary and secondary.
Moreover, we assumed that the stars have a spherical geometry and orbital parameters were taken from Table 2. As the V1878 Ori system has an eccentric and asynchronous orbit, it was necessary to adopt individual values of v e sin i, rotation periods P A,B , and inclination angles i A,B for each component in our analysis. We adopted the v e sin i values from Table 3 that were determined using spectrum synthesis. The individual rotation periods of the primary and the secondary were found by running ZDI inversions for a grid of varying rotation periods for the two stars using a step of 0.25 days. For each inversion, we computed the deviation between the observed and reconstructed Stokes V profiles, as illustrated in Fig. 5. The best combination of periods, P A = 12.81 d and P B = 13.24 d, was determined by fitting a 2-D spline function to the deviation and finding the location of the minimum. Similarly, the inclination of the two stars was constrained -more coarsely -by running inversions for a grid of inclinations for the two stars with a step of 20 degrees. We adopted the best-fitting inclination values within that grid, namely i A = i B = 60 • . This inclination value, together with the rotation periods and v e sin i that we adopted, imply stellar radii of R A = R B = 4.5 R ⊙ . This is larger than R A = R B = 3.4-3.6 R ⊙ corresponding to the parameters obtained by Covino et al. (2001b) and Marilli et al. (2007). These inclination values imply stellar masses of M A ∼ 2.4M ⊙ and M B ∼ 2.3M ⊙ using the M A,B sin 3 i values from Table 2. These are in disagreement with the stellar masses estimated from PMS evolutionary tracks in Sect. 5 (M A = 1.7 ± 0.3M ⊙ and M B = 1.6 ± 0.3M ⊙ ). While we believe the projected masses M A,B sin 3 i infered from the RV solution to be accurate, the discrepancy can stem from the inclination values which are not well-constrained by the ZDI inversion or from the masses derived with the PMS models which can potentially suffer from systematic uncertainties.
The brightness mapping was carried out modelling the Stokes I LSD profiles with our inversion code inversLSDB. This reconstruction employed Tikhonov regularization, following the approach described in Rosén et al. (2018), in order to obtain a solution with minimum brightness contrast.
The magnetic field of each star was parametrized as a combination of poloidal and toroidal fields, each decomposed using spherical harmonics as described in Kochukhov et al. (2014). The spherical harmonic coefficients α ℓ,m and β ℓ,m represent respectively the contributions of the radial and horizontal poloidal components, and γ ℓ,m the horizontal toroidal component. The parameters ℓ and m are, respectively, the angular degree and the azimuthal order of the spherical harmonic function. Here, we set the maximum degree of the spherical harmonic coefficients to ℓ max = 10. This choice of ℓ max does not restrict the complexity of the recovered field given the moderate v e sin iof the studied stars, as we let this role be played by the regularization in our analysis. Indeed, a ZDI reconstruction using Stokes V data alone is an ill-posed inverse problem, and one needs to use regularization in order to obtain a stable and unique solution. We adopted the following regularization function (Kochukhov et al. 2014) where Λ is the regularization parameter. This function penalises the presence of high order modes, hence favouring simpler magnetic field topologies. The optimal value of the regularisation parameter Λ was determined following Kochukhov (2017), by running inversions with a grid of successively decreasing regularization parameter, selecting the value that offers a good fit to the data without fitting the noise. The recovered brightness and magnetic field maps are shown in Fig. 6 while the comparison between observed and model LSD Stokes IV profiles is presented in Fig. 7. The summary of the magnetic field characteristics for each component is gathered in Table 4. The observed Stokes I and V LSD profiles are reproduced very well at all orbital phases. The presence of spots explains the apparent reversal of the line depth ratio between the primary and secondary components, observed at the rotational phases φ A between 0.529 and 0.683 (corresponding to the orbital phases 0.505-0.531 in Figs. 1 and 3). As indicated by the dashed line in Fig. 7, this feature is unexplained if one assumes a homogeneous brightness distribution, resulting in a noticeably inferior line profile fit.
The large-scale magnetic fields of the two stars have rather different topologies. On the one hand, the magnetic field of the primary is mostly poloidal (80% of the magnetic energy is concentrated in poloidal components) and non-axisymmetric (91% in the modes with |m| > ℓ/2). On the other hand, the field of the secondary is mostly toroidal (87% of the field energy is in the toroidal modes) and axisymmetric (87%). Additionaly, there is also a large difference in the average magnetic field strength of the two components. The primary has the weaker global magnetic field of the pair, with a mean field strength of 180 G and a maximum strength of 410 G. The secondary has a roughly twice stronger field, with a mean field strength of 320 G and a maximum strength of 810 G. Fig. 6 suggests that this difference is due to a strong axisymmetric toroidal field, dominating the azimuthal field map of the secondary but absent in the primary. As the two components have similar stellar parameters and rotation periods, it could be expected -if a dynamo mechanism generating their magnetic field was strictly a function of the stellar parameters -that their magnetic field would be similar. But here, as already observed before in fully-convective stars ) and rapidly rotating Sun-like stars (Rosén et al. 2018), two components of a binary system with similar stellar parameters exhibit significantly different field configurations. However, a trivial interpretation of this observation is that one or both of the stars have variable magnetic activity and the observations analysed here correspond to one particular phase of activity cycle, not necessarily representative of the time-averaged activity levels of the two components.

X-RAY DATA ANALYSIS
To study the X-ray emission of V1878 Ori observed by XMM-Newton during non-periastron (companion separation ∼ 21 R ⋆ ) and periastron (companion separation ∼ 11 R ⋆ ) we first inspected the X-ray light curves, obtained by adding all the events registered by the three EPIC instruments (Fig. 8). During both segments, the X-ray emission of V1878 Ori shows significant variability, likely due to constant flaring activity. For the strongest flare, occurring during the final 30 ks of the non-periastron segment, the X-ray count rate increases up to 30%. In addition to this short timescale variability, the average X-ray emission level dur-  ing the two observing intervals, separated by ∼ 20 d, was different: during the non-periastron phase the average Xray count rate was 0.601 ± 0.003 cts s −1 , significantly higher than 0.477 ± 0.003 cts s −1 , observed during the periastron.
To derive the characteristics of the coronal plasma we inspected the PN, M1, and M2 X-ray spectra collected during non-periastron and periastron, respectively. To increase the signal-to-noise ratio we rebinned each spectra to obtain at least 30 counts in each channel. To better constrain the Xray emitting plasma properties we simultaneously fitted the spectra of the three EPIC instruments. The spectral analysis was performed using XSPEC V12.10. We adopted a model of an optically-thin plasma emission (using the APED atomic database, Smith et al. 2001) with two isothermal components, and taking into account also the interstellar absorption. We left as free parameters the abundances of O, Ne, S, and Fe, while the abundances of the other elements were linked to the iron one.
The observed and predicted spectra are shown in Fig. 9. The values of the best-fit parameters are reported in Table 5. In both of the observing segments the temperature of the hot component is unconstrained, and only a lower limit could be obtained. The plasma properties during the two segments are very similar, with the non-periastron X-ray emission showing only a slightly larger amount of plasma emission measure (E M) in both the thermal components. From the models that we derived for V1878 Ori, we obtain X-ray luminosities of 1.6×10 31 erg s −1 and 1.2×10 31 erg s −1 in the 0.5 − 8.0 keV band, for the non-periastron and periastron phases, respectively.
A hot and X-ray bright corona, such as that of V1878 Ori, is typical for young late-type stars (e.g. Preibisch et al. 2005). In particular, V1878 Ori displays average plasma temperatures, luminosity, and abundances comparable to those of other young intermediate-mass stars with convective envelopes (e.g. Principe et al. 2017;Jensen et al. 2009).
The general picture of coronal physics of late type stars, which implies that hot coronal plasma is heated and con-

Brightness Brightness
Radial field Radial field

Meridional field Meridional field
Azimuthal field Azimuthal field fined by stellar magnetic fields, is robust. However some fundamental aspects, e.g. the extent of the magnetic structure, are still debated. Monitoring the coronal emission of eccentric binaries offers the opportunity to probe the magnetospheric extents, by investigating whether and how the two magnetospheres interact. Constraining and understanding the possible magnetospheric interactions in binary stars is also important for the understanding of star-planet magnetic interactions (e.g. Cohen et al. 2009;Fares et al. 2010;Poppenhaeger 2019).
Evidence of enhanced X-ray emission at periastron in eccentric binaries, suggesting interacting magnetospheres, was obtained for DQ Tau (Getman et al. 2011) and, tentatively, for a few other stars (Getman et al. 2016). In that study, X-ray snapshot observations indicated a 50% increase in the periastron flux relative to outside periastron but as only one non-periastron epoch was observed robust statistics could not be ascertained.
The two long X-ray observations of V1878 Ori obtained at and outside periastron indicate that, for this target and during this epoch, periastron passage was not associated with enhanced X-ray emission. Therefore the magnetospheric interaction of the two companions at periastron was not strong enough to affect the coronal activity. This lack of magnetospheric interaction in V1878 Ori, if compared to the previous results obtained for DQ Tau, may be explained by the smaller periastron separation of DQ Tau with respect to V1878 Ori (8 R ⋆ vs. 11 R ⋆ ), different field geometries, and, possibly, by an even stronger magnetic field in the cooler young low-mass component stars of DQ Tau, although these have not yet been measured for DQ Tau. The upper limit inferred for the coronal extent of V1878 Ori agrees with the value of ∼ 2 R ⋆ inferred for the coronal struc- ture dimension, obtained studying the supersaturation effects of coronal activity in young intermediate-mass stars (Argiroffi et al. 2016).

SUMMARY AND DISCUSSION
This work presented the spectropolarimetric analysis of the double-lined spectroscopic binary system V1878 Ori, consisting of two intermediate-mass weak-line T Tauri stars. We used time-resolved high-resolution spectropolarimetric observations to derive high signal-to-noise ratio LSD profiles, measure the radial velocities of the two components and refine the system's orbital solution. We obtained disentangled intensity spectra for the two stars and inferred stellar parameters with spectrum synthesis modelling. Finally, we applied the ZDI technique to the time-series of composite Stokes I and V LSD profiles to simultaneously reconstruct the brightness as well as vector magnetic field maps of V1878 Ori A and B.
We found that the secondary has a stronger global field than the primary. These fields also turned out to have different topological properties. The global magnetic geometry of the primary is predominantly poloidal and nonaxisymmetric while the field of the secondary is mostly toroidal and axisymmetric. Individually, both stars agree well with the trend that strong toroidal fields appear predominantly axisymmetric, as found by See et al. (2015) based on ensemble analysis of single-star ZDI studies.
This is the first ZDI study of a binary system consisting of two intermediate-mass T Tauri stars. The overall properties of the reconstructed magnetic fields are in broad agreement with previous studies of single T Tauri stars: the mean large-scale magnetic field (i.e. recovered from spectropolarimetric measurements subject to flux cancelation) is stronger than for hotter fully-radiative IMTTS on which the surface magnetic field is too weak and/or too complex to be detected (Villebrun et al. 2019) but weaker than the simple strong axisymmetric fields recovered in cool fully convective T Tauri stars Hill et al. (2019). Also in agreement with these studies, for T eff > 4300 K the large-scale surface magnetic topology is more complex than the simple superposition of a dipole plus an octupole. A previous ZDI study of two higher-mass intermediate-mass T Tauri stars in the Chamaeleon cluster revealed even more complex large-scale magnetic field distributions than those recovered for V1878 Ori. As those stars were similar masses but more evolved than the V1878 Ori primary and secondary stars, they have much larger radiative cores (> 60%R ⊙ ; this indicates that the field complexity in V1878 Ori will likely also increase as its radiative core developes further (Hussain et al. 2009). This study hints that the magnetic fields of intermediate-mass stars can change a lot around the convective boundary, but we have relatively few magnetic studies of stars in this region (see Villebrun et al. 2019) thus far. More magnetic studies of candidates covering this region of the PMS are needed to better understand the evolution of the stellar dynamo in young stars.
The fact that the two V1878 Ori stars, which have nearly identical stellar parameters and the same evolutionary history, have radically different magnetic field characteristics is striking. For fully convective stars, different magnetic fields at similar stellar parameters (as observed in e.g. GJ65 A and B,  can be tentatively explained by a bi-stability of the dynamo process (Gastine et al. 2013). For partially convective binary stars showing the same phenomenon, such as σ 2 CrB (Rosén et al. 2018) and V1878 Ori, such an explanation does not hold. Rather, our results may provide a clue that the spread in the ratio of toroidal to poloidal magnetic energy of partially convective stars visible in Figure 2 of See et al. (2015) is intrinsic to dynamo action in these stars and not related to a hidden parameter (such as age, composition, or initial conditions). Here, the only parameter that appeared significantly different between the two components of V1878 Ori in the spectrum synthesis is the microturbulence.
We now know that young rapidly rotating partially convective stars can exhibit a wide range of toroidal magnetic energy fractions Folsom et al. (2016Folsom et al. ( , 2018. Future observa- Table 5. Best fit parameters of the X-ray spectra of V1878 Ori collected during non-periastron and periastron phases. The plasma emission measure E M is defined as n e n H V , where n e and n H are the electron and hydrogen densities, and V is the plasma volume. Abundances refer to the solar photospheric values of Anders & Grevesse (1989). Uncertainties correspond to 1σ.  Figure 8. Background-subtracted X-ray light curves of V1878 Ori observed during non-periastron and periastron observing segments. Count rates were obtained by adding the photons registered by the three EPIC instruments and adopting a time bin of 2 ks. Error bars correspond to 1σ. Figure 9. X-ray spectra of V1878 Ori observed during non-periastron and periastron observing segments. Dots with error bars indicate the observed spectra (black, red, and blue mark PN, M1, and M2, respectively), continuous histograms represent predicted spectra. Error bars correspond to 1σ.
tions monitoring systems such as V1878 Ori for a period of years to decades will allow us to devise to which extent the differences observed in the magnetic topologies of V1878 Ori A and B can be attributed to intrinsic dynamo-related variability, whether cyclic or chaotic. Indeed, for main sequence partially convective stars, the fraction of toroidal magnetic energy as well as the degree of axisymmetry of the poloidal component can display dramatic variations along the magnetic cycle (Dunstone et al. 2008;Boro Saikia et al. 2018).
It is also important to note that the strength of the field recovered using ZDI is a small fraction of the total magnetic field strength at the surface of the star (Kochukhov & Shulyak 2019;Lavail et al. 2019;Kochukhov et al. 2020). This fraction depends on the magnetic field complexity, as (potentially strong) local tangled fields can cancel out in circular polarization. Analysis of Zee-man broadening in intensity spectra would be needed to obtain a reliable estimate of the total magnetic field strength but was not carried out with this dataset as Zeeman broadening is relatively small in the optical. Self-consistent Zeeman broadening measurements and ZDI inversions, as well as spot coverage measurements using molecular bands are becoming more feasible using high-resolution near-infrared spectropolarimeters such as SPIRou at CFHT or the upcoming CRIRES+ instrument at VLT. These will certainly help to obtain a more complete picture of the magnetism of these objects and understanding of the dynamo processes taking place in the interiors of pre-main sequence stars such as these. Canada-France-Hawaii Telescope were performed with care and respect from the summit of Maunakea which is a significant cultural and historic site.