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

HAL
open science

# Magnetospheric accretion in the intermediate-mass T Tauri star HQ Tauri

K. Pouilly, J. Bouvier, E. Alecian, S. H. P. Alencar, A.-M. Cody, J.-F. Donati,

K. Grankin, G. A. J. Hussain, L. Rebull, C. P. Folsom

► To cite this version:

K. Pouilly, J. Bouvier, E. Alecian, S. H. P. Alencar, A.-M. Cody, et al. Magnetospheric accretion in the intermediate-mass T Tauri star HQ Tauri. Astronomy and Astrophysics - A&A, 2020, 642, pp.A99. 10.1051/0004-6361/202038086. hal-02962089

HAL Id: hal-02962089

https://hal.science/hal-02962089

Submitted on 8 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
# Magnetospheric accretion in the intermediate-mass T Tauri star HQ Tauri*

K. Pouilly<fnref n="1" />, J. Bouvier<fnref n="1" />, E. Alecian<fnref n="1" />, S. H. P. Alencar<fnref n="2" />, A.-M. Cody<fnref n="3" />, J.-F. Donati<fnref n="4" />, K. Grankin<fnref n="5" />, G. A. J. Hussain<fnref n="6" />, L. Rebull<fnref n="7" />, and C. P. Folsom<fnref n="4" />

Received 3 April 2020 / Accepted 3 August 2020

## ABSTRACT

**Context.** Classical T Tauri stars are pre-main sequence stars surrounded by an accretion disk. They host a strong magnetic field, and both magnetospheric accretion and ejection processes develop as the young magnetic star interacts with its disk. Studying this interaction is a major goal toward understanding the properties of young stars and their evolution.

**Aims.** The goal of this study is to investigate the accretion process in the young stellar system HQ Tau, an intermediate-mass T Tauri star (1.9 $M_\odot$).

**Methods.** The time variability of the system is investigated both photometrically, using *Kepler*-K2 and complementary light curves, and from a high-resolution spectropolarimetric time series obtained with ESPaDOnS at CFHT.

**Results.** The quasi-sinusoidal *Kepler*-K2 light curve exhibits a period of 2.424 d, which we ascribe to the rotational period of the star. The radial velocity of the system shows the same periodicity, as expected from the modulation of the photospheric line profiles by surface spots. A similar period is found in the red wing of several emission lines (e.g., HI, CaII, NaI), due to the appearance of inverse P Cygni components, indicative of accretion funnel flows. Signatures of outflows are also seen in the line profiles, some being periodic, others transient. The polarimetric analysis indicates a complex, moderately strong magnetic field which is possibly sufficient to truncate the inner disk close to the corotation radius, $r_{cor} \sim 3.5 R_\star$. Additionally, we report HQ Tau to be a spectroscopic binary candidate whose orbit remains to be determined.

**Conclusions.** The results of this study expand upon those previously reported for low-mass T Tauri stars, as they indicate that the magnetospheric accretion process may still operate in intermediate-mass pre-main sequence stars, such as HQ Tauri.

**Key words.** stars: variables: T Tauri, Herbig Ae/Be – stars: pre-main sequence – accretion, accretion disks – stars: magnetic field – stars: individual: HQ Tau – starspots

## 1. Introduction

Classical T Tauri stars (cTTs) are young stellar objects still surrounded by an accretion disk. They possess a strong magnetic field that truncates the inner disk at a distance of a few stellar radii above the stellar surface and drives accretion through funnel flows, a process referred to as magnetospheric accretion (see reviews in, e.g., Bouvier et al. 2007a; Hartmann et al. 2016). The kinetic energy of the infalling material is dissipated in a shock at the stellar surface, creating a localized hot spot. The star-disk interaction takes place within 0.1 au or less, a scale hardly resolved by current interferometers. An alternative approach to study this compact region is to monitor the variability of the system through contemporaneous photometric, spectroscopic, and spectropolarimetric observing campaigns focused on specific targets. Over the past decade, our group has reported

several such studies (e.g., Bouvier et al. 2007b; Donati et al. 2007, 2019; Alencar et al. 2012, 2018).

This new study focuses on HQ Tau (RA = 04h35m, Dec. = +22°50'), a 1.9 $M_\odot$ T Tauri star located in the Taurus star forming region at a distance of 159 pc (Gaia Collaboration 2018). The system is moderately bright ($V \approx 12.5$; Norton et al. 2007), has a K0-K2 spectral type (Nguyen et al. 2012; Herczeg & Hillenbrand 2014), and is so far considered as single (Pascucci et al. 2015, but see Sect. 4 below). From the SuperWASP campaign, Norton et al. (2007) measured a light curve modulation with a peak-to-peak amplitude of about 2 mag in the V-band and a period of 2.4546 days. Rodriguez et al. (2017) analyzed a 9 yr long KELT light curve of the system and reported several long-duration, non-periodic dimming events with an amplitude of about 1.5 mag and lasting for weeks to months, which led them to classify this object as UXor-like (Grinin 1992). Nguyen et al. (2012) derived a radial velocity of $16.65 \pm 0.11 \text{ km s}^{-1}$, $v \sin i$ of $48 \pm 2 \text{ km s}^{-1}$, and a 10% width of the $H\alpha$ line profile of $442 \pm 93 \text{ km s}^{-1}$, indicative of ongoing accretion onto the star (see also Duchêne et al. 2017). Simon et al. (2016) measured $EW(H\alpha) = 2.22 \AA$, from which they

<footnote n="1">¹ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
e-mail: kim.pouilly@univ-grenoble-alpes.fr</footnote>

<footnote n="2">² Departamento de Fisica – ICEx – UFMG, Av. Antônio Carlos 6627, 30270-901 Belo Horizonte, MG, Brazil</footnote>

<footnote n="3">³ Bay Area Environmental Research Institute, 625 2nd St Ste. 209, Petaluma, CA 94952, USA</footnote>

<footnote n="4">⁴ Univ. de Toulouse, CNRS, IRAP, 14 avenue Belin, 31400 Toulouse, France</footnote>

<footnote n="5">⁵ Crimean Astrophysical Observatory, Nauchny, Crimea 298409, Ukraine</footnote>

<footnote n="6">⁶ European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching bei München, Germany</footnote>

<footnote n="7">⁷ Infrared Science Archive (IRSA), IPAC, California Institute of Technology, 1200 E. California Blvd, MS 100-22, Pasadena, CA 91125, USA</footnote>

<footnote>* Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii.</footnote>
derived a mass accretion rate of $\dot{M}_{acc} = 2 \times 10^{-9} M_{\odot} \text{yr}^{-1}$. More recently, from high-resolution ALMA observations, Long et al. (2019) derived an inclination of $53.8 \pm 3.2^{\circ}$ for the circumstellar disk on a scale of ~25 au. HQ Tau is the cTTs with the faintest disk of their Taurus sample in the mm range, and they report a depletion of dust toward the inner regions, perhaps indicative of an unresolved inner disk cavity. Akeson et al. (2019) similarly derived a low disk mass $M_{d} = 0.4 \times 10^{-3} M_{\odot}$ and mm flux compared to single young stars of similar mass.

We resolved to launch a monitoring campaign on this object as a representative of intermediate-mass T Tauri stars (IMTTs), having a mass at the upper range of cTTs. Moreover, its K2 light curve exhibits a clear and smooth periodicity, which makes this target amenable to the variability approach developed here. Our main goal is to investigate whether the magnetospheric accretion process, which seems to be quite common among low-mass T Tauri stars applies as well to the intermediate-mass range. In Sect. 2, we describe the photometric and spectropolarimetric datasets used in this study. In Sect. 3, we present the results of the analysis, which includes deriving stellar parameters and investigating the photometric, spectral, and polarimetric variability. We discuss the results in Sect. 4, where we provide a framework for their interpretation. We also report long-term radial velocity variations, which indicate that the system is a spectroscopic binary. We conclude on the presence of the magnetospheric accretion process in this young system of the intermediate-mass group in Sect. 5.

## 2. Observations

We briefly describe here the origin of the photometric and spectropolarimetric observations that were obtained during the HQ Tau campaign. We also describe how they were processed.

### 2.1. Photometry

HQ Tau (EPIC 247583818) was observed by *Kepler*-K2 during Campaign 13, which took place over 80 days from March 8, 2017 to May 23, 2017. The observations were performed in a broad band filter (420–900 nm) with measurements taken at a cadence of 30 minutes. The K2 light curve was reduced by A.-M. Cody (Cody & Hillenbrand 2018) and we used the PDC version in this work.

Additional photometric observations of HQ Tau were secured in the Johnson's $VR_IJ$ filters at the Crimean Astrophysical Observatory (CrAO) from December 22, 2016 to February 15, 2018, on the AZT-11 1.25m telescope equipped with the CCD camera ProLine PL23042. CCD images were bias subtracted and flat-field corrected following a standard procedure. We performed differential photometry between HQ Tau and a non-variable, nearby comparison star, 2MASS J04361769+2247125, whose brightness and colors in the Johnson system are 10.16 (V), 1.61 (V − R)$_J$, and 2.98 (V − I)$_J$. A nearby control star of similar brightness, 2MASS J04362344+2252171, was used to verify that the comparison star was not variable. It also provided an estimate of the photometric rms error in each filter, which amounts to 0.010, 0.012, and 0.011 in the $VR_IJ$ bands, respectively.

### 2.2. Spectroscopy and spectropolarimetry

High-resolution optical spectropolarimetry was obtained for HQ Tau using the Echelle SpectroPolarimetric Device for the Observation of Stars (ESPaDOnS; Donati 2003) at the Canada-

<table><caption>**Table 1.** Journal of ESPaDOnS observations.</caption><thead><tr><th>Date<br>(2017)</th><th>HJD<br>(2 450 000+)</th><th>S/N</th><th>S/N<sub>LSD</sub></th><th>φ<sub>rot</sub></th></tr></thead><tbody><tr><td>28 Oct</td><td>8054.93481</td><td>165</td><td>687</td><td>0.349</td></tr><tr><td>29 Oct</td><td>8055.92519</td><td>180</td><td>776</td><td>0.757</td></tr><tr><td>30 Oct</td><td>8056.97501</td><td>156</td><td>656</td><td>1.190</td></tr><tr><td>31 Oct</td><td>8057.91808</td><td>121</td><td>482</td><td>1.579</td></tr><tr><td>01 Nov</td><td>8058.99693</td><td>170</td><td>655</td><td>2.024</td></tr><tr><td>02 Nov</td><td>8060.01605</td><td>180</td><td>784</td><td>2.445</td></tr><tr><td>03 Nov</td><td>8060.98848</td><td>191</td><td>842</td><td>2.846</td></tr><tr><td>04 Nov</td><td>8061.91652</td><td>204</td><td>923</td><td>3.229</td></tr><tr><td>05 Nov</td><td>8063.04627</td><td>169</td><td>724</td><td>3.695</td></tr><tr><td>07 Nov</td><td>8064.96108</td><td>192</td><td>835</td><td>4.485</td></tr><tr><td>08 Nov</td><td>8065.92044</td><td>140</td><td>432</td><td>4.881</td></tr><tr><td>08 Nov</td><td>8065.96526</td><td>153</td><td>609</td><td>4.899</td></tr><tr><td>08 Nov</td><td>8066.08352</td><td>180</td><td>761</td><td>4.948</td></tr><tr><td>09 Nov</td><td>8066.88036</td><td>174</td><td>712</td><td>5.277</td></tr></tbody></table>

Notes. The columns list the date of observation, the heliocentric Julian date, the S/N by spectral resolution element at 731 nm, the effective S/N of the Stokes V LSD profiles, and the rotational phase.

France-Hawaii Telescope (CFHT) between October 28, 2017 and November 9, 2017. We obtained 14 spectra covering the 370–1000 nm range at a resolution of 68 000, reaching an S/N from 120 to 200 at 731 nm. Each observation consists in the combination of 4 individual spectra taken in different polarization states. All the data were reduced using the Libre-ESpRIT package (Donati et al. 1997), which provides an intensity Stokes I spectrum and a polarized Stokes V spectrum for each observation. The reduced spectra were normalized using a polynomial fit over selected continuum points (see Folsom et al. 2016), which produces a flat continuum, and simultaneously accounts for the spectral order overlap. The journal of observations is provided in Table 1.

## 3. Results

We describe in the following subsections the analysis of photometric, spectral, and spectropolarimetric variations observed in the HQ Tau system. We also derive the system's fundamental properties.

### 3.1. Photometric variability

The 80 day-long detrended *Kepler*-K2 light curve of HQ Tau is shown in Fig. 1. It exhibits clear periodic modulation. A Lomb-Scargle periodogram analysis (Scargle 1982; Press & Rybicki 1989) yields a period $P = 2.424$ d, where the period estimate is derived from a Gaussian fit to the periodogram's main peak, with an associated false alarm probability (FAP) of less than $10^{-4}$, computed following the Baluev (2008) approximation. Assuming that the formal error on the period is measured as the $\sigma$ of the Gaussian fit (VanderPlas 2018), it would amount to 0.028 d. However, when the K2 light curve is folded in phase, it is clear that the phase coherency between the first and last minima is lost, meaning that their location differs by at least 0.05 in phase, as soon as the period is varied by more than 0.002 d. We therefore use this empirical estimate as the true error on the period. The light curve folded in phase with this period (see Fig. 1) displays a stable, sinusoidal-like pattern with a nearly constant amplitude, which we ascribe to the modulation of the stellar flux by a cool
<figure><img src="image_3.png" /><figcaption>Fig. 1. Top: HQ Tau's detrended *Kepler* K2 light curve. Low frequency variations have been removed from the original light curve shown in Fig. 2. The spike appearing at J.D. 2.457 890.0 is presumably an instrumental defect. Middle: HQ Tau's light curve folded in phase with P = 2.424 d. The origin of phase is taken as JD 2.457 823.81 to correspond to the epoch of maximum brightness. Bottom: lomb-Scargle periodogram obtained from HQ Tau's *Kepler* K2 light curve.</figcaption></figure>

surface spot (Herbst et al. 1994). We therefore assume the photometric period P = 2.424 ± 0.002 d is the rotational period of the star. This period is consistent with the P1 = 2.42320 d period found by Rebull et al. (2020) whose work also highlighted a secondary period P2 = 5.02495 d. We do see a secondary peak in the periodogram of the detrended light curve at a period of 5.04 ± 0.13 d, with a FAP of less than 10⁻⁴ (see Fig. 1). We argue below that HQ Tau is probably a binary system, and we cannot exclude that the secondary period is to be associated with the rotation period of the companion.

Figure 2 shows the full HQ Tau light curve over 2 epochs, from December 2016 to March 2018, including the K2, before detrending, and CrAO datasets, as well as publicly available photometry from the ASAS-SN and AAVSO surveys. The K2 light curve was rescaled to the 7 CrAO measurements taken contemporaneously. The normalized K2 fluxes were thus converted to magnitudes applying a 12.3 mag zero-point offset. Over the time frame where they overlap, the amplitude of the K2 light curve appears somewhat smaller than that of the V-band CrAO light curve, presumably reflecting the longer effective wavelength of the K2 wide-band filter. Interestingly, the low-frequency part of the K2 light curve, most notably the brightness “bump” seen over its first part, is recovered over the slightly overlapping parts of the CrAO and ASAS-SN light curves, which indicates intrinsic longer term variations superimposed onto the spot modulation.

Photometric data points were also available from the ASAS-SN survey (Jayasinghe et al. 2019). HQ Tau was observed in

both the V-band and g-band filters during this period. Taking advantage of the overlap between the V- and g-band light curves, the g-band measurements were arbitrarily shifted by -0.8 mag in order to match the V-band ones. The resulting ASAS-SN light curve contains a few measurements during the first half of the ESPaDOnS observations and many additional measurements before and after the run was completed. It thus provides an estimate of the amount of variability the source exhibited at the time of the spectroscopic observations, performed about 6 months after the K2 monitoring. Overall, it suggests relatively mild variability at this epoch, amounting to about 0.2 mag in amplitude in the V-band. A significant dimming event, which lasted for about a month just prior to the start of the spectroscopic observations had apparently ceased by the time of the ESPaDOnS measurements. We notice, however, that the mean flux level was changing over the spectroscopic run, with the source being about 0.3 mag brighter toward the end of the run. Both the amplitude of variability and the low-frequency variations are not unlike those seen in the spot-driven K2 light curve obtained several months earlier. This suggests that at the time of the ESPaDOnS observations, the source was in a state of relatively mild periodic variability, and probably not in a drastic dipper state as reported at some other epochs by Rodriguez et al. (2017), in spite of the occurrence of a recent dimming event.

Indeed, we verified that the ASAS-SN g-band and V-band datasets restricted to the 7 days following the ESPaDOnS run (JD 2.458 067.8-2.458 073.7) exhibit a smoothly varying light curve consistent with a period of 2.424 days when folded in phase. Using the same epoch for the origin of phase as for the K2 light curve above (i.e., JD 2.457 823.81), we find that the photometric minimum of the ASAS-SN dataset occurs around phase ~0.56, which we estimated by interpolating the two lowest photometric measurements occurring at phase 0.47 and 0.66, respectively. Hence, the ASAS-SN photometric minimum is not far from phase 0.50 of the K2 photometric minimum, the slight difference being easily accounted for by the 0.002 d uncertainty on the K2 period. The near conservation of phase over the 155 day-long temporal gap stretching between the end of the K2 observations and the beginning of the ESPaDOnS run suggests the modulation results from a long lived, relatively stable spot structure. For the rest of the paper, we thus use the following ephemeris:

$$HJD(d) = 2457823.81 + 2.424E, \quad (1)$$

where E is the rotational phase of the system.

Color variations associated to the brightness changes are shown in Fig. 3. In both the (V-R)j and (V-I)j colors, the source becomes redder when fainter. The color slope of the small scale variations likely results from spot modulation (Herbst et al. 1994; Venuti et al. 2015). A single deep faintening event was recorded on JD 2.458 025.6 in the VRIj filters of the CrAO dataset, and is confirmed by the single-filter ASAS-SN light curve (see Fig. 2). The color plot shows that as the brightness decreased by about 1.2 mag in the V-band, the system became much redder, with a color slope close to that expected for extinction by ISM-like grains. This suggests that the dimming event was caused by circumstellar dust crossing the line of sight.

## 3.2. Stellar parameters

We used the ESPaDOnS high resolution spectra to derive HQ Tau's stellar parameters, namely effective temperature (Teff), and the radial (Vr), rotational (vsini), and microturbulent (vmic)
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** HQ Tau's *V*-band light curve. CrAO measurements are shown as filled grey circles, the K2 light curve as a continuous line, AAVSO measurements as filled triangles, and the ASAS-SN rescaled *g*- and *V*-band datasets as crosses and open circles, respectively. The *upper panel* shows the light curve over 2 observing seasons while the *lower one* displays the photometric variability of the source around the CFHT/ESPaDOnS observations.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Color-magnitude plots from CrAO VR<sub>J</sub>I<sub>J</sub> photometry. Left: (V - R)<sub>J</sub> vs. V. Right: (V - I)<sub>J</sub> vs. V. The source gets redder when fainter. The dashed line is the expected reddening line for ISM extinction. The photometric measurements are color-coded by Julian date to distinguish the 2 observing seasons: lighter points correspond to the first epoch, darker points to the second one (see Fig. 2). We note that the deep faintening event of about 1.2 mag in the *V*-band is accompanied by a strong reddening of the system.</figcaption></figure>

We averaged the 14 ESPaDOnS spectra gathered during the campaign and fit synthetic spectra calculated with the ZEEMAN code (Landstreet 1988; Wade et al. 2001; Folsom et al. 2012) based on MARCS stellar atmosphere grids (Gustafsson et al. 2008), VALD line lists (Ryabchikova et al. 2015), and including the same oscillator strength corrections as those used in Folsom et al. (2016). We explored a range of T<sub>eff</sub>, v sin i, V<sub>r</sub>, and v<sub>mic</sub>, and obtained the best fit to HQ Tau's mean spectrum through a χ<sup>2</sup> minimization procedure using a Levenberg-Marquardt algorithm (see Folsom et al. 2013). The fit was performed on 11 independent spectral windows spanning the range from 422 to 754 nm, each with a width ranging from 4 to 10 nm. Each win-

dow contains well resolved, relatively unblended photospheric lines, and is devoid of emission features. The spectral windows are shown in Fig. 4.

In order to derive the stellar properties from spectral fitting, we first fixed the macroturbulent velocity to 2 km s<sup>-1</sup>, log g to 4.0, and assumed solar metallicity, in other words, values that are typical of low-mass PMS stars in the solar neighborhood (Padgett 1996; James et al. 2006; Santos et al. 2008; Taguchi et al. 2009; D'Orazi et al. 2011). An example of the resulting fit is shown in Fig. 4. Once stellar parameters were obtained for each spectral window, we averaged them and computed the rms dispersion, and then removed the windows that yielded results
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Upper six rows: spectral windows selected from HQ Tau's mean ESPaDOnS spectrum to derive stellar parameters. Bottom row: as an example, HQ Tau's mean spectrum (blue) fitted with a ZEEMAN synthetic spectrum (orange) in the 547–557 nm window.</figcaption></figure>

beyond $1\sigma$ from the mean value. Usually, between 3 or 4 windows were thus rejected. We thus derive $T_{\text{eff}} = 4997 \pm 160$ K, $v \sin i = 53.9 \pm 0.9$ km s$^{-1}$, $V_r = 6.64 \pm 0.71$ km s$^{-1}$, and $v_{\text{mic}} = 1.4 \pm 0.2$ km s$^{-1}$. The $T_{\text{eff}}$ of 4997 K we derive is consistent with a spectral type K0-K1 according to Pecaut & Mamajek (2013) conversion scale.

We next obtained an estimate of HQ Tau's luminosity, based on its 2MASS J-band magnitude ($J = 8.665 \pm 0.024$, Skrutskie et al. 2006), using the Gaia parallax ($\pi = 6.304 \pm 0.213$ mas, Gaia Collaboration 2018). Extinction was computed from the ($V-J$) color excess, using intrinsic colors from Pecaut & Mamajek (2013), which yields $A_V = 2.60 \pm 0.1$ mag and $A_J = 0.73 \pm 0.03$ mag, consistent with the value found by Herczeg & Hillenbrand (2014). We applied bolometric corrections from Pecaut & Mamajek (2013), and obtained $L_\star = 3.91 \pm 0.65 L_\odot$, and $R_\star = 2.64 \pm 0.31 R_\odot$.

We then plotted HQ Tau in an Hertzsprung-Russell (HR) diagram and used a grid of CESTAM evolutionary models (Marques et al. 2013; Villebrun et al. 2019) to derive its mass and internal structure (see Fig. 5). From bilinear interpolation in model grids, we obtain $M_\star = 1.87_{-0.55}^{+0.21} M_\odot$, and an age of ~2 Myr. According to evolutionary models, the star, is partly radiative with $M_{\text{rad}} = 0.55_{-0.55}^{+0.35} M_\star$ and $R_{\text{rad}} = 0.51_{-0.51}^{+0.14} R_\star$, where $M_{\text{rad}}$ and $R_{\text{rad}}$ are the mass and radius of the radiative core, respectively.

We find that HQ Tau is a fast rotator with $v \sin i = 53.9$ km s$^{-1}$. Combining the stellar radius with the star's $v \sin i$

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** HQ Tau's position in the HR diagram. The red cross illustrates HQ Tau's position with corresponding uncertainties. The black curves are the evolutionary tracks from CESTAM PMS models, with the corresponding mass indicated on the left. The blue curves show the position where the radiative core of the star reaches 75%, 50% and 25% of the stellar mass.</figcaption></figure>

and the K2 rotational period, we derive the inclination angle of the rotational axis onto the line of sight, $i = 75_{-17}^{+15}$ deg, indicative of a highly inclined system.

### 3.3. Spectroscopic properties and variability

We investigate the spectral variability of the system by analyzing the 14 ESPaDOnS spectra taken over 12 nights. Table 2 summarizes radial velocity and line equivalent width measurements.

#### 3.3.1. Radial velocity variations

Radial velocity measurements were obtained by cross-correlating each HQ Tau's spectrum with the spectrum of a spectral template. We used the ESPaDOnS observations of the Hyades cluster member Melotte 25-151 as a template, a slowly rotating ($v \sin i = 4.8$ km s$^{-1}$) K2-type star with $V_r = 37.98$ km s$^{-1}$ and $T_{\text{eff}} = 4920$ K (Folsom et al. 2018), reduced and normalized in the same way as HQ Tau. We computed the cross correlation function (CCF) over 5 spectral windows (542–547, 558–563, 585–587, 608–613, and 639–649 nm), and fit a Gaussian profile to derive the radial velocity difference between HQ Tau and the template. We averaged the results over the 5 spectral windows to get a mean value of $V_r$ and its rms uncertainty. The results are listed in Table 2 and the radial velocity curve is shown in Fig. 6.

The radial velocity appears to be modulated and a sinusoidal fit yields a period of $2.48 \pm 0.16$ d using the $1\sigma$ confidence level on $\chi^2$ minimization, consistent with the stellar rotation period within uncertainties. We therefore ascribe this modulation to surface spots. The phased radial velocity curve is shown in Fig. 6. Indeed, the sinusoidal fit indicates that the mean radial velocity ($\langle V_r \rangle = 7.22 \pm 0.27$ km s$^{-1}$) occurs around phase 0.6, which is expected when the spot faces the observer (Vogt & Penrod 1983) and is consistent with the photometric minimum of the contemporaneous ASAS-SN light curve. The amplitude of the $V_r$ modulation amounts to hardly a tenth of the star's vsini, which suggests it is mostly driven by a large area and high latitude cold spot. We also notice a regular downward drift of $V_r$ with
<table><caption>**Table 2.** Photospheric radial velocity and equivalent widths of selected emission lines.</caption><thead><tr><th rowspan="2">HJD<br>(2 450 000+)</th><th rowspan="2">V<sub>r</sub><br>(km s<sup>-1</sup>)</th><th rowspan="2">σV<sub>r</sub></th><th colspan="4">EW (Å)</th></tr><tr><th>Hα</th><th>Hβ</th><th>CaII<br>λ8542</th><th>CaII<br>λ8662</th></tr></thead><tbody><tr><td>8054.93481</td><td>9.77</td><td>0.36</td><td>3.25</td><td>0.96</td><td>1.66</td><td>1.42</td></tr><tr><td>8055.92519</td><td>6.76</td><td>0.53</td><td>2.40</td><td>0.09</td><td>1.16</td><td>0.96</td></tr><tr><td>8056.97501</td><td>8.23</td><td>0.40</td><td>3.18</td><td>0.96</td><td>1.59</td><td>1.36</td></tr><tr><td>8057.91808</td><td>8.48</td><td>0.57</td><td>3.84</td><td>1.30</td><td>1.80</td><td>1.63</td></tr><tr><td>8058.99693</td><td>7.63</td><td>0.41</td><td>1.77</td><td>0.06</td><td>1.09</td><td>0.94</td></tr><tr><td>8060.01605</td><td>8.37</td><td>0.41</td><td>2.88</td><td>0.88</td><td>1.61</td><td>1.46</td></tr><tr><td>8060.98848</td><td>5.58</td><td>0.65</td><td>2.05</td><td>0.12</td><td>1.25</td><td>1.17</td></tr><tr><td>8061.91652</td><td>7.97</td><td>0.33</td><td>2.96</td><td>0.84</td><td>1.52</td><td>1.46</td></tr><tr><td>8063.04627</td><td>7.47</td><td>0.72</td><td>2.64</td><td>0.59</td><td>1.46</td><td>1.31</td></tr><tr><td>8064.96108</td><td>7.35</td><td>0.69</td><td>3.91</td><td>1.03</td><td>1.88</td><td>1.65</td></tr><tr><td>8065.92044</td><td>5.44</td><td>0.67</td><td>1.36</td><td>0.28</td><td>0.92</td><td>0.82</td></tr><tr><td>8065.96526</td><td>4.82</td><td>0.78</td><td>1.37</td><td>0.29</td><td>0.92</td><td>0.91</td></tr><tr><td>8066.08352</td><td>5.20</td><td>0.54</td><td>1.64</td><td>0.06</td><td>1.13</td><td>1.05</td></tr><tr><td>8066.88036</td><td>6.88</td><td>0.34</td><td>3.37</td><td>1.16</td><td>1.60</td><td>1.46</td></tr></tbody></table>

an amplitude of 2–3 km s⁻¹ over a timescale of 10 days. We will come back to this feature in Sect. 4.

### 3.3.2. Emission line profiles

We computed residual emission line profiles by subtracting the rotationally broadened photospheric profiles of the template Melotte 25-151 from HQ Tau’s profiles. Hα, Hβ, and the CaII infrared triplet (IRT) are the only lines in HQ Tau’s spectrum exhibiting significant emission flux. The 3 lines of the CaII IRT are similar in shape but the line at 866.2 nm may be affected by a hydrogen line of the Paschen series and the line at 849.8 nm exhibited a variation far lower than the other two. We therefore focused on the 854.2 nm component. Figure 7 shows the Hα, Hβ, and the selected CaII IRT residual line profiles. The raw line profiles are shown in Fig. A.1.

The Hα line profile exhibits high velocity wings to the blue and red sides, up to about 300 km s⁻¹. The red wing displays large variability, with at times signs of high velocity redshifted absorptions reaching below the continuum, which are Inverse P Cygni (IPC) profiles. The blue wing also displays significant variability. The Hβ line is dominated by deep IPC signatures extending up to about 300 km s⁻¹ seen in nearly half of the observations. The CaII line also exhibits strong variability from the line center to its red wing, with however only marginal IPC components below the continuum level.

We computed periodograms in each velocity bin across the line profiles, the result of that is a 2D periodogram for each line, which is shown in Fig. 7. In all 3 lines, a clear periodicity is seen at a frequency of 0.4 d⁻¹ (P = 2.5 d) with a FAP reaching ~10⁻⁴ and a typical value of ~10⁻², corresponding to the stellar rotational period, extending from the line’s central peak all the way to the red wing. This is a clear indication that the IPC components are modulated by stellar rotation, as expected for funnel flow accretion. We do notice significant power as well at about the same frequency at highly blueshifted velocities (~200 km s⁻¹) in the Hα profile, which might also be present in the CaII line profile. The peak at a frequency of 0.6 d⁻¹ is the 1-day alias of the 0.4 d⁻¹ frequency. It appears clearly here, though at a higher FAP of 0.1, due to the night-to-night sampling of ESPaDOnS spectra.

Figure 8 shows the residual Hα and Hβ profiles ordered by day and by phase. The phase ordering illustrates well the periodicity of IPC components: they appear over the 5 rotational cycles covered by the spectral series at specific phases, from 0.69 to 0.02 in both profiles, with a maximum depth around phase 0.90. The depth of the IPC components appears to slightly vary from one rotational cycle to the next, being stronger at phase 4.90 than at phase 2.85 for instance.

We investigate the relationship between the various components seen in the line profiles by computing correlation matrices (Johns & Basri 1995a; Oliveira et al. 2000; Alencar & Batalha 2002; Kurosawa et al. 2005). Correlation matrices consist of the computation of Pearson’s linear correlation coefficient on line intensity between 2 velocity channels of the same or different line profiles. The coefficient approaches 1 for a strong correlation, 0 for no correlation or −1 for anticorrelated intensity variations. Correlation matrices are shown in Fig. 9.

A strong correlation is seen between velocity channels within the red wing of all profiles, and within the blue wing for Hα and CaII. These correlations appear as bright squares in the correlation matrices, located on each side of the line center and extending to high velocities. As discussed above, the strong correlation within the red side of the line profiles is probably related to the periodic appearance of IPC components, which extends from close to the line center up to +300 km s⁻¹. Strikingly, little correlation is seen between the red and blue wings of the line profiles, suggesting that their variability is driven by unrelated physical processes, presumably accretion for redshifted velocity channels and outflows for blueshifted ones. Yet, an interesting anti-correlation appears between a restricted range of blueshifted velocity channels, from about −130 to −220 km s⁻¹, and the redshifted part of the line profile. This is clearly seen in the Hα correlation matrix and also appears in the CaII matrix. This indicates that, as the IPC components appear on the red side of the line profile, the high velocity part of the blue wing becomes more intense, perhaps the signature of an accretion-driven high velocity outflow (e.g., Johns & Basri 1995b).

### 3.3.3. The NaI D profile

The NaI doublet is sensitive to wind signatures (Mundt 1984), and can thus help to better understand the accretion-ejection connection. The NaI D lines are seen in absorption in HQ Tau’s spectrum and exhibit significant variability (see Fig. 10). As the two lines of the doublet have similar profiles and behavior, we focus here on the 589.0 nm line (NaI D2), leaving aside the 589.6 nm twin line, which lies at the edge of a spectral order.

Two prominent, relatively narrow, and apparently stable absorption components are seen around the line center. The stable component located at +10 km s⁻¹ (V<sub>helio</sub> ≈ 17 km s⁻¹) was previously reported by Pascucci et al. (2015) and ascribed to local interstellar absorption. This narrow and stable absorption component is also present at the same redshifted velocity in all our spectra in the KI 770 nm line profile. We notice however significant variability on both the blue and red sides of this component in the NaI D2 line, over a range of about ±10 km s⁻¹, which indicates an additional source of variable absorption, linked to the stellar system itself. These transient components are also seen in the KI profile. The second narrow absorption component is blueshifted and located at −20 km s⁻¹ (V<sub>helio</sub> ≈ −13 km s⁻¹). This component is also apparent in the NaI D line profile reported in Pascucci et al. (2015), but is not discussed there. CI Tau, which is located only 0.5 degrees away from HQ Tau, does not exhibit this second component, while it
<figure><img src="image_8.png" /></figure>

**Fig. 6.** *From top to bottom:* radial velocity, $H\alpha$ line flux, and average longitudinal magnetic field as a function of Julian date (**left**) and rotational phase (**right**). The sinusoidal fit, forced at the K2 photometric period, to radial velocity variations is shown as a black dash-dotted curve in the **top panels**. The $B_l = 0$ G line is drawn for reference in the **bottom panels**.

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

**Fig. 7.** *Top:* residual $H\alpha$ (**left**), $H\beta$ (**middle**), and CaII 854.2 nm (**right**) emission lines of HQ Tau spectra. *Bottom:* corresponding 2D periodograms with the mean profile superimposed as a white curve. The maximum power of the periodogram is in yellow and the minimum in dark blue. The red and blue wings of $H\alpha$ and the red wing of $H\beta$ and CaII, as well as the center of the 3 line profiles, exhibit a periodic modulation at $f \approx 0.4$ d$^{-1}$, corresponding to the rotational period of the star. The power seen at $f \approx 0.6$ d$^{-1}$ is the 1-day alias induced by the observational window.

does display the redshifted one. This suggests the blueshifted narrow absorption component is not due to interstellar cloud absorption. Yet, its narrow width compared to photospheric lines and its stability suggest it is not related to the inner variable system. This component could conceivably be the signature of the distant, cold part of a stellar wind, as discussed in Mundt (1984).

The variability seen across the NaD line profile is relatively complex, with different components. The 2D periodogram shown in Fig. 10 reveals that the red side of the profile, from about +40 to +150 km s$^{-1}$, varies periodically, at the stellar rotation period. The comparison of successive NaI D2 profiles ordered in phase, shown in Fig. 11, suggests this is due to the appearance of IPC components seen between phase 0.69 and 0.02, as reported above for Balmer and CaII emission lines. In contrast, the variability seen in the blue wing of the line profile is not periodic. Sporadic blueshifted absorption components appear in nearly half of the observations, over velocity channels ranging from about −120 to −30 km s$^{-1}$. Figure 11 shows the occurrence of deep, high velocity blueshifted absorption components from JD 2 458 060 to 2 458 066 covering velocity channels from about −120 to −50 km s$^{-1}$, with evidence for a gradual shift in velocity over this timescale. Lower velocity absorption components

are also seen for a couple of days from JD 2 458 056 to 2 458 057 around −60 to −30 km s$^{-1}$. None of these blueshifted components appear to last over the whole observing period. They appear to be transient phenomena occurring only during specific cycles, are of limited duration, and are not modulated by stellar rotation. These are presumably the signatures of episodic outflows.

Finally, Fig. 9f presents the correlation matrix of the Na D2 line. It resembles the matrices of the other emission lines in the red part of the profile, being similarly modulated by the IPC components. It does not show any correlated variability in the blue wing, as none is expected from the sporadic, independent outflows components analyzed above.

### 3.3.4. Mass accretion rate

We computed the mass accretion rate using the relationship between the accretion luminosity and the luminosity in the $H\alpha$ residual line (Alcalá et al. 2017). Duchêne et al. (2017) notice that the $H\alpha$ emission line of HQ Tau might not be a good accretion indicator, due to its weakness, though its width and double peaked shape favor an accretion origin. In order to confirm the mass accretion rate values found with $H\alpha$, we also computed the
<figure><img src="image_10.png" /><figcaption>Fig. 8. Residual Hα (left panels) and Hβ (right panels) profiles sorted by rotational phase and by day of observation. The successive rotational cycles are shown by different colors, as in Fig. 6.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Fig. 9. Correlation matrices of the line profiles. (a) Hα vs. Hα, (b) Hβ vs. Hβ, (c) Hβ vs. Hα, (d) CaII vs. CaII, (e) Hα vs. CaII, (f) NaD vs. NaD.</figcaption></figure>

mass accretion rate using the Hβ, and CaII IRT residual lines. The line flux was computed from the line equivalent width as $F_{line} = F_0 \cdot EW \times 10^{-0.4m_{\lambda}}$, where $F_{line}$ is the line flux, $F_0$ is the reference flux in the selected filter, EW the line equivalent width, and $m_{\lambda}$ is the extinction-corrected magnitude of HQ Tau in the selected filter, namely $R_C$ for Hα, $R_J$ for the Ca II IRT, and B for Hβ. The line luminosity is then derived from $L_{line} = 4\pi d^2 F_{line}$,

<figure><img src="image_12.png" /><figcaption>Fig. 10. Top: NaI D2 line appearing in the 14 spectra of HQ Tau. The black dash-dotted line is a synthetic NaI D2 line produced using the ZEEMAN code with HQ Tau's parameters. Bottom: 2D periodogram of the NaI D2 line.</figcaption></figure>

where d is the distance to HQ Tau. The accretion luminosity follows from the line luminosity by $log(L_{acc}) = a \log(L_{line}) + b$, where a and b are coefficients taken from Alcalá et al. (2017). The mass accretion rate is obtained from:

$$L_{\text{acc}} = \frac{GM_{\star}\dot{M}_{\text{acc}}}{R_{\star}}\left[1 - \frac{R_{\star}}{R_{l}}\right], \qquad (2)$$

where $R_l$ is the magnetospheric truncation radius, typically $5R_\star$ for cTTs (Bouvier et al. 2007a).

As the equivalent width of Balmer lines is strongly affected by the IPC components, we computed the mean value of $L_{acc}$ and $\dot{M}_{acc}$ by averaging the results excluding the observations showing strong IPC profiles. The mass accretion rates we deduce from the various lines are consistent within $2\sigma$. This yields $L_{acc} = 0.019 \pm 0.005 L_{\odot}$ and $\dot{M}_{acc} = 1.26 \pm 0.35 \times 10^{-9} M_{\odot} \text{yr}^{-1}$ for the system. As EW(Hα) is relatively weak in HQ Tau, the emission line flux might include a non negligible chromospheric contribution. In order to quantify it, we considered a sample
<figure><img src="image_13.png" /><figcaption>Fig. 11. NaD2 line profiles sorted by rotational phase (left) and by day of observation (right). The successive rotational cycles are shown by different colors, as in Fig. 6.</figcaption></figure>

of wTTSs with similar spectral type and v sin i and (or) rota-
tion period as HQ Tau, which were chosen from the list of
James et al. (2006). We thus derived $\langle EW(H\alpha) \rangle = 0.7 \pm 0.1$ Å
for these wTTSs, assumed this value to be the chromospheric
contribution, and subtracted it from HQ Tau’s EW(Hα). We
then recomputed $\dot{M}_{\text{acc}}$ from the different line estimates to obtain
$\dot{M}_{\text{acc}} = 1.16 \pm 0.35 \times 10^{-9} M_{\odot} \text{yr}^{-1}$. The temporal variation of
the Hα line flux, a proxy for $\dot{M}_{\text{acc}}$, is shown in Fig. 6. It is clearly
modulated and a sinusoidal fit yields a period of 2.4 d, consistent
with the stellar rotational period. The peak-to-peak amplitude of
the modulation is about a factor 3, typical of cTTs $\dot{M}_{\text{acc}}$ variabil-
ity on these timescales (Venuti et al. 2014).

Another signature of accretion is the so-called veiling. Veiling is an additional continuum component emitted by the accretion shock, which fills the photospheric lines. In order to compute the amount of veiling in HQ Tau’s spectrum, we used the non-accreting star Melotte 25-151 as a template. Varying the amount of continuum line filling over the spectral window 639–649 nm in the template spectrum, we matched HQ Tau’s spectrum by a $\chi^2$ minimizing method:

$$
I(\lambda) = \frac{I_t(\lambda)(1+r)}{1 + I_t(\lambda)r}, \qquad (3)
$$

where *I*<sub>*t*</sub> is the intensity spectrum of Melotte 25-151, the tem-
plate, and *r* the fractional veiling that is the excess continuum
flux divided by the stellar continuum flux at a given wavelength.
We derived a weak veiling at 640 nm ranging from 0.15 to 0.20
in the 14 HQ Tau spectra, as expected for the modest mass accre-
tion rate and bright photosphere of the system, with uncertainties
of order of 0.20, thus preventing us from detecting any signifi-
cant temporal variations.

3.4. Spectropolarimetry

We analyzed the ESPaDOnS spectropolarimetric data to study
the magnetic field properties of HQ Tau. We used the Least
Square Deconvolution method (LSD) from Donati et al. (1997)

to compute the mean Stokes *I* and *V* photospheric line pro-
files. This method increases the S/N by extracting the Zeeman
signature in many photospheric lines and averaging them. The
parameters used for the LSD computation are the mean wave-
length, intrinsic line depth, and Landé factor, set at 640 nm, 0.2,
and 1.2, respectively, as in Donati et al. (2010). We extracted
a list of spectral lines from VALD atomic database for a star
with HQ Tau’s fundamental parameters, and selected absorption
lines located between 450 and 850 nm. We then removed lines
contaminated with emission, blended with strong broad lines, or
affected by telluric lines. The LSD Stokes *I* and *V* profiles are
shown in Fig. 12. The typical rms is 3.5 × 10<sup>-4</sup> and 1.4 × 10<sup>-4</sup> for
the Stokes *I* and *V* profiles, respectively, computed from more
than 9000 photospheric lines.

The shape of the Stokes I profiles displays a clear modula-
tion over the rotational cycle. A periodogram analysis shows the
same peak at a frequency 0.4 d⁻¹ as for emission lines, rang-
ing from −25 km s⁻¹ to 60 km s⁻¹, with a FAP reaching 10⁻⁴.
This indicates that the modulation is periodic at the stellar rota-
tional period, as expected from a large rotating dark spot at the
stellar surface. The profile is nearly symmetric at phase 3.69,
confirming the spot is facing the observer around phase 0.65, as
defined above from the V<sub>r</sub> curve (see Sect. 3.3.1). The Stokes V
signal is much weaker than the Stokes I profile. Nevertheless, a
clear magnetic signature is seen at several phases, with a com-
plex shape.

The average surface longitudinal magnetic field, $B_l$, was
computed from the Stokes I and V profiles, as:

$$
B_l = -2.14 \times 10^{11} \times \frac{\int vV(v)dv}{\lambda gc \int (1 - I(v))dv}, \quad (4)
$$

where $B_l$ is given in Gauss, $v$ is velocity relative to line center, $\lambda$ in nm is the mean wavelength, and $g$ the Landé factor chosen for the LSD computation (Donati et al. 1997; Wade et al. 2000). We computed the error by propagating the uncertainties in the trapezoidal integration over the range $\pm 60$ km s$^{-1}$.

Figure 6 shows the variation of the longitudinal magnetic field as a function of time and rotational phase. $B_l$ varies from –70 G to +68 G. We notice that the minimum of $B_l$ is reached at $\phi \sim 0.7$, consistent with the rotational phase at which the dark spot that modulates the radial velocity faces the observer, and reaches a maximum around phase 0.9 when the IPC components are best seen in the emission line profiles. Apart from these extrema, there is little evidence for $B_l$ to be modulated, possibly the result of a complex magnetic topology at the stellar surface and (or) the strongest field regions being hidden in dark spots. We ran a similar analysis in the CaII IRT 854.2 nm line profile. However, we did not detect any significant Stokes V signature within the line profile and our data can only provide a $3\sigma$ upper limit of 450 G to the longitudinal magnetic field component in the line.

The full Zeeman Doppler Imaging (ZDI, Donati et al. 2011, 2012) analysis was performed based on the LSD Stokes I and V profiles, using the code described in Folsom et al. (2018). The resulting best fit of the Stokes I and V profiles is shown in Fig. 12.

This first step provides a Doppler image shown in Fig. 13. The fitting process starts from a uniform brightness distribution across the stellar disk and iteratively adds darker and brighter features onto the stellar surface to reproduce the shape of Stokes I profile. For each cell on the stellar surface, the local profile is assumed to be a Voigt profile and the brightness of the pixel is adjusted by minimizing $\chi^2$ and maximizing entropy.
<figure><img src="image_14.png" /><figcaption>**Fig. 12.** HQ Tau's Stokes *I* (left), *V* (middle), and Null (right) profiles (black) ordered by increasing fractional phase, and the fit produced by the ZDI analysis (red). The fractional number next to each profile indicates the rotational phase.</figcaption></figure>

The solutions for brightness and magnetic field distributions are not unique based on $\chi^2$. Therefore, an additional constraint is added using entropy. ZDI uses the algorithm of Skilling & Bryan (1984) for maximizing entropy while minimizing $\chi^2$. We adopted 0.66, 2.8 km s⁻¹, and 1.9 km s⁻¹ for the limb darkening coefficient, the Gaussian and Lorentzian widths of the Voigt profile, respectively. The latter values were chosen based on theoretical values found in Gray (2005) for a K0 spectral type star around 600 nm and adapted to fit the observed profiles. The Doppler image reveals a large cold spot extending over phases 0.4 to 0.8 in longitude, ranging from 90° to about 50° in latitude, and covering about 20% of the stellar visible hemisphere.

ZDI then takes the brightness map as input and fits the Stokes *V* profiles by adjusting the spherical harmonic coefficients that describe the magnetic field (see Donati et al. 2006). It is also possible to fit both intensity and polarisation profiles simultaneously but the resulting maps are essentially the same. Taking advantage of the large $v \sin i$ of the star, which provides spatial resolution on the stellar surface, the spherical harmonic expansion was carried out to the 15th order in *l* (Folsom et al. 2016). Letting the stellar parameters freely vary in the fitting process, we obtain an inclination of $50 \pm 5^\circ$, a period of $2.453_{-0.008}^{+0.006}$ d, and a $v \sin i$ of $51.2 \pm 0.5$ km s⁻¹. All uncertainties were obtained using the 68.27 % confidence level on $\chi^2$. All these values are within $3\sigma$ of those we derived in the previous sections, although the inclination we derive from ZDI is on the lower side of the estimate we obtained in Sect. 3.2 by combining the rotational velocity and period with the stellar radius. This suggests that the stellar radius may be underestimated by about 20%. The solution that minimizes $\chi^2$ did not allowed us to constrain the differential rotation $d\Omega$. We thus fixed $d\Omega = 0.0$ rad d⁻¹ to derive the other parameters, as the derived period was consistent within $3\sigma$ with the photometric one.

The large-scale magnetic reconstruction we recover from the analysis of LSD profiles reveals a mostly poloidal field, with the toroidal component contributing 25% to the total magnetic energy. The main structure of the radial magnetic field extends from phase 0.3 to 0.8, from 90° to 60° in latitude, and reaches a strength of -562 G locally. The dipolar field amounts to only 14% of the poloidal component, is tilted by 31.5° from the rotation axis, and reaches a maximum strength of 63 G at the stellar photosphere.

<figure><img src="image_15.png" /><figcaption>**Fig. 13.** Upper panel: surface brightness map of HQ Tau. The color scale goes from dark brown to light yellow, for dark to bright features. Lower panels, from top to bottom: radial, azimuthal, and meridional surface magnetic field distribution. Ticks on top of the panels indicate the phase of observations.</figcaption></figure>

The quadrupole and octupole components contribute respectively 9.7% and 8.3% to the total magnetic flux.
The complex topology and moderate strength of HQ Tau's magnetic field is reminiscent of those reported by Hussain et al. (2009) and Villebrun et al. (2019) for IMTTs.

We caution that the values reported here from the LSD analysis are lower limits to the strength of the magnetic field that interacts with the disk. Previous similar studies have shown that the ZDI analysis of emission lines arising from the accretion shock located at the foot of funnel flows, most notably HeI 588 nm, yields a much stronger magnetic field strength than that deduced for the large-scale field at the photospheric level from LSD analysis (e.g., Donati et al. 2019). Unfortunately, HQ Tau's spectrum does not show the post-shock HeI line in emission, owing to the modest accretion rate producing a relatively weak shock seen against a bright photosphere.

## 4. Discussion

We selected HQ Tau among many Taurus sources observed by K2 (Rebull et al. 2020) for a spectropolarimetric follow-up campaign in order to investigate the magnetospheric accretion process in a representative member of the class of relatively massive T Tauri stars, the so-called IMTTs (Calvet et al. 2004). With a mass of $M_{\star} = 1.9 M_{\odot}$, a short rotational period of 2.424 d and $v \sin i = 53.9$ km s$^{-1}$, HQ Tau's properties are intermediate between the cool low-mass T Tauri stars and the hotter intermediate-mass Herbig Ae stars. Indeed, PMS models suggest HQ Tau's interior is already partly radiative, with $R_{\text{rad}}/R_{\star} \sim 0.51$, as the star is transiting from the Hayashi to the Henyey track in the HR diagram. The mass accretion rate we derive, $\dot{M}_{\text{acc}} \sim 1.16 \times 10^{-9} M_{\odot} \text{yr}^{-1}$, is however relatively low and more typical of low-mass TTS (e.g., Mendigutía et al. 2011). Only few such IMTTs had been monitored for their magnetic properties so far (e.g., Hussain et al. 2009) and the goal here was to explore the possible extension of the magnetospheric accretion process that is ubiquitous among low-mass T Tauri stars to the higher mass range.

HQ Tau is undoubtedly a member of the Taurus star forming region (Luhman 2018). Yet, the mean radial velocity we measured, $\langle V_r \rangle = 7.22 \pm 0.27$ km s$^{-1}$, is significantly different from the radial velocity distribution of Taurus members, $\langle V_r \rangle = 16.3 \pm 6.43$ km s$^{-1}$ (Bertout & Genova 2006). This led us to suspect that HQ Tau might not be a single source, as usually assumed, and the small but regular drift in $V_r$ seen during the ESPaDOnS run may be additional evidence for that (see Fig. 6). Historically, Simon et al. (1987) reported the system to be a tight binary from lunar occultation, with a separation of $4.9 \pm 0.4$ mas, later revised to $9.0 \pm 2$ mas by Chen et al. (1990). However, the former authors state that “the binary nature of HQ Tau is not obvious” and indeed this result was not confirmed by following studies (Richichi et al. 1994; Simon et al. 1995, 1996; Mason 1996).

We searched for additional radial velocity measurements in the literature. They are summarized in Table 3, together with new measurements, kindly obtained for us recently by L. Hillenbrand and H. Bouy. The results are illustrated in Fig. 14. HQ Tau exhibits clear velocity variations. Measurements obtained from 2006 to 2017 show $V_r$ oscillations between about 7 and $22$ km s$^{-1}$ on a timescale of years. However, the two most recent measurements taken 50 days apart in late 2019, show the same amplitude. This suggests that HQ Tau is a short period type 1 spectral binary (SB1), whose orbit determination awaits additional measurements. We note that the companion must be relatively faint compared to the primary in the optical, as we see no evidence for a double-lined system in the Stokes I LSD profiles.

<table><caption>**Table 3.** $V_r$ values from the literature, and from our and new measurements.</caption><thead><tr><th>HJD<br>(-2450000)</th><th><i>V</i><sub>r</sub><br>km s<sup>-1</sup></th><th>σ<i>V</i><sub>r</sub><br>km s<sup>-1</sup></th><th>Ref.<br>—</th><th>Instrument<br>—</th></tr></thead><tbody><tr><td>3887.77</td><td>16.3</td><td>0.02</td><td>1</td><td>MIKE</td></tr><tr><td>4467.83</td><td>9.2</td><td>1.2</td><td>2</td><td>HIRES</td></tr><tr><td>6262.77</td><td>17.2</td><td>0.2</td><td>3</td><td>HIRES</td></tr><tr><td>8060.00</td><td>7.22</td><td>0.15</td><td>4</td><td>ESPaDOnS</td></tr><tr><td>8816.77</td><td>21.9</td><td>0.5</td><td>2</td><td>HIRES</td></tr><tr><td>8866.44</td><td>7.92</td><td>3.74</td><td>5</td><td>SOPHIE</td></tr></tbody></table>

**References.** (1) Nguyen et al. (2012); (2) L. Hillebrand (priv. comm.); (3) Pascucci et al. (2015); (4) this work; (5) H. Bouy (priv. comm.).

<figure><img src="image_16.png" /><figcaption>**Fig. 14.** Radial velocity of HQ Tau as a function of date of observation. The error bars are indicated, and can be smaller than the symbol size.</figcaption></figure>

Nevertheless, the detection of a secondary period at 5.03 d in the K2 light curve (see Sect. 3.1) suggests that the companion may account for a low-level contribution to the optical flux of the system. We caution that we have implicitly assumed in the above analysis that the companion's contribution to the emission line flux and to the ZDI reconstruction process could be neglected. Although this assumption may not be fully verified, we have currently no way to estimate the flux contribution of the companion to the system. Finally, we suggest that the faint mm flux of HQ Tau's disk and its suspected inner cavity (Long et al. 2019; Akeson et al. 2019) could conceivably be related to a low-mass companion orbiting within the inner disk and, at least partly, clearing it.

On a much shorter timescale of days, the radial velocity of the system is modulated by surface spots with a period consistent with the photometric period of 2.424 d derived from the K2 light curve. From the additional ASAS-SN photometry obtained at the time of the spectroscopic campaign, we find that the photometric minimum occurs around $\phi = 0.54$, which is close to the rotational phase when the spot which modulates radial velocity variations faces the observer (at $\phi = 0.6$ in Fig. 6). We deduce that a large, cool spot is mainly responsible for both the photometric and radial velocity variability on a timescale of a few days. From the flat-bottom shape of Stokes I LSD profiles at phase 0.69, the relatively low amplitude of variability (~5% in both photometric flux and $V_r$), and the sinusoidal shape of the K2 light curve, we further infer that the cool spot is presumably situated at a high latitude on the star. Indeed, the Doppler image reconstruction indicates the presence of a cool spot straddling the stellar pole and covering about 20% of the stellar surface.

In order to check if a cold spot can produce the observed modulation of radial velocity and photometric flux, we used SOAP2 (Dumusque et al. 2014) to generate synthetic variability curves. We used the stellar parameters of HQ Tau reported above and ran a set of simulations for a range of spot latitudes,
spot sizes, and temperature differences between the spotted area and the stellar photosphere. The strong degeneracy between the spot location and its size does not allow us to strictly constrain these parameters. Nevertheless, we found a solution reproducing the amplitude of modulation of both radial velocity and flux with a spot located at a latitude of 60° and whose projected surface covers 30% of the visible hemisphere, qualitatively consistent with the Doppler map. These spot properties are not uncommon for T Tauri stars (e.g., Bouvier & Bertout 1989).

Most of the circumstellar diagnostics also seem to be modulated at the stellar rotation period, consistent with the expectations for magnetospheric accretion. The line profiles exhibit periodically modulated IPC components, whose maximum occurs around phase 0.9, and last for nearly half of the rotational cycle from phase 0.69 to 0.02. This is consistent with magnetospheric accretion funnel flows covering a significant azimuthal extension around the central star, as expected for a tilted large-scale magnetosphere. The shape of the profiles themselves are reminiscent of those computed from magnetospheric accretion models. For instance, the mean $H\beta$ profile of HQ Tau is quite similar to that computed by Muzerolle et al. (2001) for $\dot{M}_{\text{acc}} = 10^{-9} M_{\odot} \text{yr}^{-1}$ at an inclination of 60°, which shows a pronounced IPC component for a funnel flow temperature of 8000 K. The shape of the $H\alpha$ profile is comparable to that computed by Lima et al. (2010) for this mass accretion rate.

The apparent mass accretion rate itself is modulated at the stellar rotational period. As it scales with the equivalent width of the line profiles, this is mostly the result of the periodic appearance of IPC components, which reduce the line flux. Hence, $\dot{M}_{\text{acc}}$ variations are primarily related to the geometric projection of the corotating accretion funnel flow. We do not have evidence for significant intrinsic variations of $\dot{M}_{\text{acc}}$ over a timescale of a couple of weeks, although we note that the depth of the IPC components slightly varies over this timescale. The deepest IPC components occur around $\phi = 0.9$, which is significantly later than the epoch at which the cold spot that modulates the radial velocity faces the observer. If the cold spot marks the magnetic pole where the funnel flow is anchored, it might indicate that the upper part of the funnel flow that produces the redshifted absorptions is trailing. This may occur if the magnetospheric truncation radius is located beyond the disk's corotation radius.

Such a phase delay has been previously reported for T Tauri stars, such as AA Tau (Bouvier et al. 2003) and V2129 Oph (Alencar et al. 2012). Unfortunately, the veiling is too weak in HQ Tau to trace the location of the accretion shock from the modulation of optical excess flux. Indeed, a weak veiling is not unexpected against a relatively bright photosphere ($T_{\text{eff}} = 4997$ K) at such a low $\dot{M}_{\text{acc}}$. Similarly, we do not see a clear modulation of the longitudinal component of the magnetic field, presumably due to its complex topology, which prevents us from assigning a rotational phase to the magnetic pole. It may therefore be that the accretion shock at the base of the funnel flow is located at a slightly different longitude than the cold dark spot around the stellar pole, thus accounting for the IPC components appearing around $\phi = 0.9$ without having to resort to a twisted magnetospheric accretion column.

Together with evidence for magnetospheric funnel flows, some line profiles reveal signatures of outflows. The $H\alpha$ line profile exhibits a high velocity blueshifted absorption component, from about −200 to −150 km s$^{-1}$, which is modulated at the stellar rotational period. Interestingly, this component appears to be anti-correlated with the appearance of IPC features in the line profile. This is reminiscent of the “egg beater” model introduced by Johns & Basri (1995b) for SU Aur, where funnels flow and

mass outflows occur on opposite azimuths at the disk inner edge: while accretion funnel flow is favored where the magnetic pole is inclined toward the disk, the situation is reversed at the opposite azimuth and favors interface or inner disk winds. This interpretation would be consistent with the overall variability reported here for the $H\alpha$ profile of HQ Tau.

Transient mass loss episodes are seen on a timescale of a few days in the NaI D2 line profile. Several deep absorption components are seen over blueshifted velocities ranging from about −40 down to −120 km s$^{-1}$. These absorptions are not rotationally modulated. Instead, they are seen over a few consecutive days and then disappear altogether (see Fig. 11). These episodic mass outflows could conceivably be related to magnetic reconnections at the star-disk interface, following the magnetospheric inflation scenario put forward by Bouvier et al. (2003) for AA Tau: as the magnetospheric field lines inflate under the shear of differential rotation, they eventually open up and reconnect (Goodson et al. 1997). Such an inflationary cycle first induces a phase of reduced mass accretion onto the star during inflation, followed by a transient outflow during reconnection, and the restoration of funnel flow accretion.

A detailed analysis of the line profile variability may support the magnetospheric inflation scenario. The transient blueshifted absorptions are most conspicuous in the NaI D2 line profile over 2 rotational cycles, from J.D. 8060 to 8065 (see Fig. 11), and the central velocity of these components appear to slowly drift toward the line center over this timeframe. The IPC components seen in Balmer lines are the deepest on J.D. 8066 ($\phi = 4.88, 4.90$), meaning right after the transient outflow episode. This would be consistent with the restoration of a magnetospheric accretion funnel after reconnection of the inflated magnetosphere. In contrast, at the start of the transient outflow episode, on J.D. 8061 ($\phi = 2.85$), the IPC profiles are the weakest, consistent with an inflated state for the magnetosphere. As discussed in Bouvier et al. (2003) and Alencar et al. (2018), the timescale for the magnetospheric inflation cycle is expected to be several times the rotational period of the star (Zanni & Ferreira 2013), which is presumably the reason why it is rarely observed in spectroscopic time series extending over only a couple of weeks.

Although we detect only relatively weak Stokes V signatures in the LSD spectra of HQ Tau, they are sufficient to derive the longitudinal component of the magnetic field and reconstruct surface magnetic maps. We do not detect rotational modulation of the longitudinal magnetic field, unlike what is usually seen in T Tauri stars (e.g., Donati et al. 2019, 2020). The rotational modulation of $B_l$ in photospheric lines is often complex, except in the dipolar magnetic field case. It may also be that the weakness of the signatures and the complexity of the field topology at the stellar surface combine to hide the modulation. As HQ Tau has already developed a significant radiative core, we do not expect a surface magnetic field dominated by a dipolar component (Gregory et al. 2012).

The full ZDI analysis suggests a weak dipolar component, inclined by 31.5° from the rotational axis, and amounting to about 63 G. However, cancellation of opposite polarities on small spatial scales at the stellar surface may affect the total magnetic energy of the star. Also, strong field regions may be hidden in dark stellar spots. Therefore, the total magnetic energy we measure with ZDI on photospheric LSD profiles is to be taken as a lower limit (Lavail et al. 2019; Sokal et al. 2020). Indeed, the magnetic field strength derived from Stokes V signatures in the CaII and (or) HeI emission line profiles of T Tauri stars usually exceeds that deduced from LSD profiles by at least a factor of 5 to 6 because they probe small and localised highly magnetised
regions on the surface of the star (Donati et al. 2019, 2020). In HQ Tau, the only emission line available to investigate additional magnetic tracers is CaII IRT. As this line is formed at least partly in the accretion funnel flow, it is a potential probe of the specific field connection between the star and the disk. From the lack of detectable Stokes V signal, we could merely derived a 3σ upper limit of ~450 G for the longitudinal field component in the emission line core.

The two intermediate mass T Tauri stars previously imaged in this way, CR Cha ($P_{\text{rot}}$ = 2.3 d) and CV Cha ($P_{\text{rot}}$ = 4.4 d) by Hussain et al. (2009), have similar spectral types and masses as HQ Tau but have larger radiative cores ($M_{\text{rad}}/M_{\star}$ ≈ 0.65 and 0.92 respectively) and are therefore probably further along in their PMS evolution. In fact, HQ Tau and CR Cha also have very similar rotation periods and mass accretion rates, with $\dot{M}_{\text{acc}}$ = 2 × 10⁻⁹ $M_{\odot}$yr⁻¹ for CR Cha (Nisini et al. 2018). It is interesting to note that there are strong similarities between all three stars in their surface activity, as shown by comparing their brightness maps and their large scale magnetic field distributions. All three stars show a large cool spot near their poles, with CR Cha showing an almost-identical off-centre-polar cap like HQ Tau. The large scale magnetic fields of all three systems are all complex and non-axisymmetric.

Are our results consistent with magnetospheric accretion being at work in HQ Tau? Unfortunately, due to the lack of Zeeman signatures in emission lines, we cannot derive from our data the strength of the dipolar magnetic field in the accretion shock, which would allow us to compute the magnetospheric truncation radius, $r_{\text{mag}}$, according to Bessolaz et al. (2008)'s prescription:

$$ \frac{r_{\text{mag}}}{R_{\star}} = 2m_{s}^{2/7} B_{\star}^{4/7} \dot{M}_{\text{acc}}^{-2/7} M_{\star}^{-1/7} R_{\star}^{5/7}, \quad (5) $$

where $m_s$ ≈ 1, $B_\star$ is the equatorial magnetic field strength, $\dot{M}_{\text{acc}}$ is the mass accretion rate, $M_\star$ the stellar mass, and $R_\star$ the stellar radius, respectively in units of 140 G, 10⁻⁸ $M_{\odot}$yr⁻¹, 0.8 $M_{\odot}$, and $2R_\odot$. Instead, we may assume that the truncation radius is located close to the corotation radius, $r_{\text{cor}}$ = 3.55 ± 0.35 $R_\star$. We then derive from Eq. (5) a magnetic field strength of ~120 G at the equator, which translates to 240 G at the pole for the dipolar component (and to ~3 G at the truncation radius, following Gregory 2011). Although we cannot measure this component directly, we note that the value required at the stellar surface is not inconsistent with the large-scale ZDI magnetic maps derived above nor with the upper limit we set on the CaII line field strength.

An independent support to the magnetospheric accretion scenario comes from the maximum velocity of IPC components measured in the emission line profiles, $v_{\text{IPC}}^{\text{max}}$ ~ 330 ± 19 km s⁻¹. We measured this velocity by fitting the red part of Hβ's IPC components by a straight line. The location where the line joins the continuum yields to $v_{\text{IPC}}^{\text{max}}$. This provides an estimate of the free-fall velocity of the accreted material projected onto the line of sight. Accounting for projection effects that combines the star's inclination (i ~ 60°, as the average between the rotational and ZDI estimates above) and the magnetic obliquity (β = 31.5°), we have $v_{\text{IPC}}^{\text{max}}$ ≈ $v_{\text{ff}} \cos(i - \beta)$, which yields $v_{\text{ff}}$ = 376 km s⁻¹. Assuming free fall accretion from $r_{\text{cor}}$ to $R_\star$, the material would hit the stellar surface at a velocity $v_{\text{ff}}$ = (2GM_\star/R_\star)¹/²(1 - $R_\star/R_{\text{cor}}$)¹/² = 424 km s⁻¹, which is consistent with the estimate obtained from the maximum redshifted velocity of the IPC components, and thus supports a magnetospheric truncation radius being located close to the disk's corotation radius.

Finally, while we observed HQ Tau during a relatively quiescent phase of variability, we note that its behavior on longer

timescales may be more complex. In particular, the deep and long lasting UXOr events reported for the system by Rodriguez et al. (2017) must have another origin than the variability described here. The high inclination we derive for the system, i = 75±15°, could favor transient circumstellar extinction events, such as UXOr and (or) AA Tau type. One of these deep faintening events occurred just before our spectroscopic observations, reaching nearly 1.5 magnitudes and lasting for a month, with significant intra-variability. The system became much redder during the event, consistent with obscuration by circumstellar dust. We may envision two related scenarios to explain this additional component to the system's variability. One is that it might result from a sudden change in the vertical scale-height of the inner disk. As the inner disk edge lies at a distance of only ~3 $R_\star$ from the stellar surface, and the system is seen at high inclination, a puffed-up inner disk could conceivably obscure the star for the duration of the instability (Turner et al. 2010). Another possibility is a change in the location of the magnetospheric truncation radius relative to the dust sublimation radius. From the expression of Monnier & Millan-Gabet (2002), we derive a sublimation radius $r_{\text{sub}}$ = 5.55 ± 1.00 $R_\star$, using $Q_R$ = 1 and $T_s$ = 1500 K. As $r_{\text{sub}}$ is larger than $r_{\text{mag}}$, there is no dust at the magnetospheric truncation radius, which accounts for the lack of dipper-like occultations at the time of our (and K2) observations. However, should the truncation radius increase, following either an increase in the magnetic field strength (by a factor of 3) or a decrease in the mass accretion rate (by a factor of 8), to eventually reach the sublimation radius, a dusty disk warp would result and could induce occultation events in this highly inclined system. We therefore suggest that the episodic UXOr events regularly observed in this system are due to the combination of close-in circumstellar dusty material around the central star and the specific geometry under which the system is seen.

## 5. Conclusions

Following the K2 Taurus campaign, we monitored the young, intermediate mass HQ Tau system with spectropolarimetry in order to investigate the accretion and ejection processes on a timescale of days. The results of this campaign provide clear diagnostics of magnetically mediated accretion occurring in the system. We observed redshifted absorptions in emission line profiles that are periodically modulated at the stellar rotational period, indicative of magnetospheric funnel flows passing through the line of sight. This is expected from a global dipolar topology of the stellar magnetic field at few stellar radii, despite the fact that, at the stellar surface, the Stokes V signatures indicate a more complex topology. While the star hosts a modest dipolar magnetic field component, the reduced mass accretion rate and the star's rapid rotation result in the magnetospheric truncation radius being close to the disk's corotation radius, as observed in most T Tauri systems investigated so far. Concurrent with accretion diagnostics, spectral signatures of outflows are clearly seen in line profiles, some being modulated by stellar rotation, others being transient phenomena, possibly resulting from instabilities at the star-disk magnetospheric boundary. We thus conclude that intermediate-mass pre-main sequence systems may undergo similar accretion and ejection processes as lower mass T Tauri stars, including stable magnetospheric accretion funnel flows, in spite of moderate magnetic field strength. Systems seen at high inclination also experience longer term faintening events, which result from the occultation of the inner system by circumstellar dusty material.

HQ Tau is a fascinating pre-main sequence system, reported here to possibly be a short period spectroscopic binary whose
orbit remains to be determined. It exhibits a mix of variability behaviors, including accretion signatures modulated on a timescale of a few days as seen in lower mass T Tauri stars, as well as longer term UXOR events more typical of Herbig Ae stars. The richness of variability patterns this intermediate-mass young star displays warrants additional multi-wavelength studies of the system and its circumstellar environment on the long term.

**Acknowledgements.** We thank an anonymous referee whose comments improved the content of this paper. We thank Lynne Hillenbrand and Hervé Bouy for taking additional high resolution spectra of HQ Tau that helped us assessing its binary nature. We acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. This paper includes data collected by the *Kepler* mission. Funding for the *Kepler* mission is provided by the NASA Science Mission directorate. This work has made use of the VALD database, operated at Uppsala University, the Institute of Astronomy RAS in Moscow, and the University of Vienna. This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement No 742095; SPIDI: Star-Planets-Inner Disk-Interactions; http://www.spidi-eu.org). This project was funded in part by INSU/CNRS Programme National de Physique Stellaire and Observatoire de Grenoble Labex OSUG2020. S.H.P. Alencar acknowledges financial support from CNPq, CAPES and Fapemig.

## References

Akeson, R. L., Jensen, E. L. N., Carpenter, J., et al. 2019, ApJ, 872, 158

Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20

Alencar, S. H. P., & Batalha, C. 2002, ApJ, 571, 378

Alencar, S. H. P., Bouvier, J., Walter, F. M., et al. 2012, A&A, 541, A116

Alencar, S. H. P., Bouvier, J., Donati, J. F., et al. 2018, A&A, 620, A195

Baluev, R. V. 2008, MNRAS, 385, 1279

Bertout, C., & Genova, F. 2006, A&A, 460, 499

Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155

Bouvier, J., & Bertout, C. 1989, A&A, 211, 99

Bouvier, J., Grankin, K. N., Alencar, S. H. P., et al. 2003, A&A, 409, 169

Bouvier, J., Alencar, S. H. P., Boulelier, T., et al. 2007a, A&A, 463, 1017

Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., & Romanova, M. M. 2007b, Protostars and Planets V, 479

Calvet, N., Muzerolle, J., Briceño, C., et al. 2004, AJ, 128, 1294

Chen, W. P., Simon, M., Longmore, A. J., Howell, R. R., & Benson, J. A. 1990, ApJ, 357, 224

Cody, A. M., & Hillenbrand, L. A. 2018, AJ, 156, 71

Donati, J. F. 2003, ASP Conf. Ser., 307, 41

Donati, J.-F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658

Donati, J.-F., Howarth, I. D., Jardine, M. M., et al. 2006, MNRAS, 370, 629

Donati, J. F., Jardine, M. M., Gregory, S. G., et al. 2007, MNRAS, 380, 1297

Donati, J. F., Skelly, M. B., Bouvier, J., et al. 2010, MNRAS, 409, 1347

Donati, J. F., Gregory, S. G., Alencar, S. H. P., et al. 2011, MNRAS, 417, 472

Donati, J. F., Gregory, S. G., Alencar, S. H. P., et al. 2012, MNRAS, 425, 2948

Donati, J. F., Bouvier, J., Alencar, S. H., et al. 2019, MNRAS, 483, L1

Donati, J. F., Bouvier, J., Alencar, S. H., et al. 2020, MNRAS, 491, 5660

D'Orazi, V., Biazzo, K., & Randich, S. 2011, A&A, 526, A103

Duchêne, G., Becker, A., Yang, Y., et al. 2017, MNRAS, 469, 1783

Dumusque, X., Boisse, I., & Santos, N. C. 2014, ApJ, 796, 132

Folsom, C. P., Bagnulo, S., Wade, G. A., et al. 2012, MNRAS, 422, 2072

Folsom, C. P., Likuski, K., Wade, G. A., et al. 2013, MNRAS, 431, 1513

Folsom, C. P., Petit, P., Bouvier, J., et al. 2016, MNRAS, 457, 580

Folsom, C. P., Bouvier, J., Petit, P., et al. 2018, MNRAS, 474, 4956

Gaia Collaboration 2018, VizieR Online Data Catalog: I/345

Goodson, A. P., Winglee, R. M., & Böhm, K.-H. 1997, ApJ, 489, 199

Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres, 3rd edn. (Cambridge University Press)

Gregory, S. G. 2011, Am. J. Phys., 79, 461

Gregory, S. G., Donati, J. F., Morin, J., et al. 2012, ApJ, 755, 97

Grinin, V. P. 1992, A&AT, 3, 17

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135

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

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

Hussain, G. A. J., Collier Cameron, A., Jardine, M. M., et al. 2009, MNRAS, 398, 189

James, D. J., Melo, C., Santos, N. C., & Bouvier, J. 2006, A&A, 446, 971

Jayasinghe, T., Stanek, K. Z., Kochanek, C. S., et al. 2019, MNRAS, 485, 961

Johns, C. M., & Basri, G. 1995a, AJ, 109, 2800

Johns, C. M., & Basri, G. 1995b, ApJ, 449, 341

Kurosawa, R., Harries, T. J., & Symington, N. H. 2005, MNRAS, 358, 671

Landstreet, J. D. 1988, ApJ, 326, 967

Lavail, A., Kochukhov, O., & Hussain, G. A. J. 2019, A&A, 630, A99

Lima, G. H. R. A., Alencar, S. H. P., Calvet, N., Hartmann, L., & Muzerolle, J. 2010, A&A, 522, A104

Long, F., Herczeg, G. J., Harsono, D., et al. 2019, ApJ, 882, 49

Luhman, K. L. 2018, AJ, 156, 271

Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74

Mason, B. D. 1996, AJ, 112, 2260

Mendigutía, I., Calvet, N., Montesinos, B., et al. 2011, A&A, 535, A99

Monnier, J. D., & Millan-Gabet, R. 2002, ApJ, 579, 694

Mundt, R. 1984, ApJ, 280, 749

Muzerolle, J., Calvet, N., & Hartmann, L. 2001, ApJ, 550, 944

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

Nisini, B., Antoniucci, S., Alcalá, J. M., et al. 2018, A&A, 609, A87

Norton, A. J., Wheatley, P. J., West, R. G., et al. 2007, A&A, 467, 785

Oliveira, J. M., Foing, B. H., van Loon, J. T., & Unruh, Y. C. 2000, A&A, 362, 615

Padgett, D. L. 1996, ApJ, 471, 847

Pascucci, I., Edwards, S., Heyer, M., et al. 2015, ApJ, 814, 14

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

Press, W. H., & Rybicki, G. B. 1989, ApJ, 338, 277

Rebull, L. M., Stauffer, J. R., Cody, A. M., et al. 2020, AJ, 159, 6

Richichi, A., Leinert, C., Jameson, R., & Zinnecker, H. 1994, A&A, 287, 145

Rodriguez, J. E., Ansdell, M., Oelkers, R. J., et al. 2017, ApJ, 848, 97

Ryabchikova, T., Piskunov, N., Kurcz, R. L., et al. 2015, Phys. Scr., 90, 054005

Santos, N. C., Melo, C., James, D. J., et al. 2008, A&A, 480, 889

Scargle, J. D. 1982, ApJ, 263, 835

Simon, M., Howell, R. R., Longmore, A. J., et al. 1987, ApJ, 320, 344

Simon, M., Ghez, A. M., Leinert, C., et al. 1995, ApJ, 443, 625

Simon, M., Holfeltz, S. T., & Taff, L. G. 1996, ApJ, 469, 890

Simon, M. N., Pascucci, I., Edwards, S., et al. 2016, ApJ, 831, 169

Skilling, J., & Bryan, R. K. 1984, MNRAS, 211, 111

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

Sokal, K. R., Johns-Krull, C. M., Mace, G. N., et al. 2020, ApJ, 888, 116

Taguchi, Y., Itoh, Y., & Mukai, T. 2009, PASJ, 61, 251

Turner, N. J., Carballido, A., & Sano, T. 2010, ApJ, 708, 188

VanderPlas, J. T. 2018, ApJS, 236, 16

Venuti, L., Bouvier, J., Flaccomio, E., et al. 2014, A&A, 570, A82

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

Villebrun, F., Alecian, E., Hussain, G., et al. 2019, A&A, 622, A72

Vogt, S. S., & Penrod, G. D. 1983, PASP, 95, 565

Wade, G. A., Donati, J. F., Landstreet, J. D., & Shorlin, S. L. S. 2000, MNRAS, 313, 851

Wade, G. A., Bagnulo, S., Kochukhov, O., et al. 2001, A&A, 374, 265

Zanni, C., & Ferreira, J. 2013, A&A, 550, A99

## Appendix A: Raw emission line profiles

The average $H\alpha$, $H\beta$, and CaII IRT line profiles of HQ Tau are shown in Fig. A.1. These figures show the lines before subtracting the photospheric template profiles.

<figure><img src="image_17.png" /><figcaption>Fig. A.1. From left to right: mean $H\alpha$, $H\beta$, and CaII IRT HQ Tau's line profiles (blue) overplotted on those of the photospheric template Melotte 25-151 (dash dotted black). The template has been corrected for the radial velocity shift and rotationally broadened to the $v\sin i$ of HQ Tau.</figcaption></figure>