EDP Sciences
Home All issues Volume 582 (October 2015) A&A, 582 (2015) A38 Full HTML
Free Access
Issue
A&A
Volume 582, October 2015
Article Number A38
Number of page(s) 12
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201525771
Published online 02 October 2015

© ESO, 2015

1. Introduction

Observational studies of magnetism in late-type stars have evolved dramatically during the past two decades from the classical chromospheric activity diagnostics (e.g. Mount Wilson H-K project, Baliunas et al. 1995) to spectro-polarimetric snapshot surveys (the BCool project, Marsden et al. 2014) and detailed long-term magnetic monitoring (e.g. Morgenthaler et al. 2012). This has been enabled by the advent of improved instrumentation (e.g. ESPaDOnS at CFHT, Donati 2003), together with advanced data analysis techniques for detection (e.g. least squares deconvolution, Donati et al. 1997) and mapping (e.g. Zeeman Doppler imaging, Donati & Brown 1997) of magnetic fields. These recent studies have opened new possibilities for different areas of astrophysical research, in particular, on dynamo processes and the origin of stellar magnetic fields across the HR diagram (see Donati & Landstreet 2009). In the case of main sequence solar-type stars, complex and relatively weak (≪1 kG) large-scale surface field topologies have been reported. The appearance of ring-like structures of significant (and even dominant) toroidal fields in these stars seems to be connected with the rotation period, hence with the dynamo mechanism behind their generation (Petit et al. 2008).

From the perspective of exoplanet studies, activity-related signatures in late-type stars are known to affect the detection techniques in the form of radial-velocity jitter and photometric flicker (Bastien et al. 2014). By simulating the effects induced by active regions and spots in Sun-like stars, recent tools have been developed to estimate and remove their contribution from the observations (e.g. Dumusque et al. 2014). Other studies have considered the detectability of planets around active cool stars by modelling the stellar activity from recovered starspot and magnetic field maps (Jeffers et al. 2014a; Donati et al. 2014). It is clear that knowing the characteristics of the stellar magnetic field is crucial for addressing the presence of exo-planets in a given system (e.g. the case of ϵ Eridani, Anglada-Escudé & Butler 2012; Jeffers et al. 2014b).

In addition, the stellar magnetic field dominates the environment around late-type stars. This includes transient events, such as flares and coronal mass ejections (Shibata & Magara 2011), and the development of persistent solar-like winds and astrospheres (Wood 2004). These phenomena are known to have a profound impact on the structure of exoplanet atmospheres, a critical factor in the habitability of exo-planetary systems (see Vidotto et al. 2013; Cohen et al. 2011b, 2014). They can erode atmospheres through thermal evaporation or non-thermal processes, such as sputtering and ion pick-up. Significant mass loss has been detected for exoplanets that is driven by the stellar wind (Vidal-Madjar et al. 2003; Linsky et al. 2010; Jensen et al. 2012). Models of this interaction require accurate knowledge of the wind properties (Ekenbäck et al. 2010; Sanz-Forcada et al. 2011) and, therefore, of the host-star surface magnetic field. However, given the observational limitations, robust surface field distributions are known for a very limited number of Sun-like planet-hosting stars (e.g. Fares et al. 2013; Fares 2014). In this context, detailed studies of these systems are very valuable resources not only as direct stellar counterparts of our solar system, but also for the growing interest in finding suitable Earth-like life supporting places in the Universe.

In this article we present our detailed study of one such planet-hosting Sun-like star (HD 1237), in which we investigate the large scale magnetic field and chromospheric activity. This is the first step in characterising the impact of the stellar magnetic field on the circumstellar environment around this system. In particular, the conditions and possible interactions via the magnetically driven stellar wind, with the Jupiter-size exoplanet that comes as close as 0.25 AU in its orbit (Naef et al. 2001). In Sect. 2, we summarise the main properties of the star. Details of the observations and calibration procedures are given in Sect. 3. We present the activity diagnostics and variability in Sect. 4. Section 5 contains a description of the implemented technique for extracting the magnetic field signatures from the spectro-polarimetric data. The required steps for the imaging procedure and the resulting surface field maps are presented in Sect. 6. In Sect. 7, we discuss our findings in the context of previous and ongoing studies of solar-type stars. Our main conclusions are summarised in Sect. 8.

2. HD 1237 stellar properties

HD 1237 (GJ 3021) is a bright (Vmag = 6.58), Sun-like star (G8V) located about 17.5 pc from the Sun in the southern constellation of Hydrus (Koen et al. 2010). This object is a relatively young (~0.88 Gyr), chromospherically active, and confirmed exoplanet host star. Naef et al. (2001) used the enhanced chromospheric activity to explain the large residuals arising from the best Keplerian orbital solution of the planet (Mpsin(i) = 3.37 ± 0.14M♃, Porb = 133.7 ± 0.2 days, e = 0.51 ± 0.02, a = 0.49 AU).

Table 1 contains the basic stellar properties of HD 1237 taken from Ghezzi et al. (2010), Torres et al. (2006), and Saffe et al. (2005). Rotation period (Prot) estimates are sparse, ranging from ~4.0 to 12.6 days, as summarised by (Watson et al. 2010). As presented in Sect. 6.1, we obtain Prot = 7.0 ± 0.7 from our observations (Sect. 3). In addition, we estimated the radial velocity, vR, and the rotational velocity, vsini, using an automatic spectral classification tool (MagIcS, Donati et al. 2012) and the fundamental properties of the star (Table 1). Using several of our observed spectra, we found on average a vsini of 5.3 ± 1.0 km s-1 and a vR of −5.2 ± 0.2 km s-1. Literature values for vsini range between ~4.5 km s-1 and 5.5 km s-1 (Naef et al. 2001; Torres et al. 2006; Schröder et al. 2009). However, as indicated by Naef et al. (2001), the vsini value can be over-estimated for metal-rich stars such as HD 1237. For the subsequent analysis, we adopted a vsini of 5.3 km s-1, which is consistent, within the errors, with the value reported by Naef et al. (2001) and Torres et al. (2006).

Table 1

HD 1237 basic properties.

3. Observational data

We obtained observations using the polarimetric mode (Piskunov et al. 2011) of the HARPS echelle spectrograph (Mayor et al. 2003) at the ESO 3.6 m telescope at La Silla Observatory. The wavelength coverage of the observations range from 378 nm to 691 nm, with a 8 nm gap starting at 526 nm.

Data were reduced using the REDUCE package (Piskunov & Valenti 2002; Makaganiuk et al. 2011), which was modified for the HARPS instrument configuration. This package produces an optimal extraction of the bias-subtracted spectra after flat-fielding corrections and cosmic ray removal have been carried out. The continuum level is determined by masking out the strongest, broadest features (e.g. the Balmer lines) and then fitting a smooth slowly varying function to the envelope of the entire spectrum. Spectra are obtained with resolutions varying from 95 000 to 113 000, depending on the wavelength, with a median value of 106 000. Uncertainties are derived for each pixel assuming photon statistics. The star was observed at two epochs separated by 5 months (July and December) in 2012. A summary of the observations is presented in Table 2.

Table 2

Journal of observations.

The exposure times listed correspond to one circularly polarised spectrum (Stokes V), which results from combining four individual sub-exposures using the ratio method. As explained in Donati et al. (1997), the polarisation signal is obtained by dividing spectra with perpendicular (orthogonal) polarisation states (for HARPSpol and Stokes V: 45°, 135°, 225°, and 315°, using the quarter waveplate). Additionally, a null-polarisation spectrum is constructed to check for possible spurious polarisation contributions in the observations. More details can be found in Bagnulo et al. (2009). Owing to bad weather conditions, two consecutive Stokes V spectra were added together for the night of 2012 July 23.

4. Magnetic activity and variability

To characterise the chromospheric activity level of the star during the observed epochs, we used the CaII H (396.8492 nm) & K (393.3682 nm) lines and the classic Mount Wilson S-index, SMW, defined as (1)Here, H and K represent the fluxes measured in each of the Ca line cores using 0.105 nm wide spectral windows, and R and V are the fluxes measured in the continuum over 2 nm windows centred at 390.1 nm and 400.1 nm, respectively, on both sides of the CaII region.

4.1. Index calibration

To compare the activity level of HD 1237 with other stars, we need to convert the measured HARPS S-index, SH, to the Mount Wilson scale. For this, we require a calibration factor, α, which is an instrument-dependent quantity that linearly relates the values for the classic SMW and the HARPS fluxes H, K, R, V: (2)We estimated α by including a set of stars with previous measurements of chromospheric activity via SMW (Santos et al. 2000) and within the HARPS observations’ database. The spectral type of the reference stars and their reported SMW values are listed in Table 3. The linear relation between SMW and the HARPS fluxes is plotted in Fig. 1. The derived calibration factor, within the 1σ uncertainty, is (3)All spectra were co-aligned using a high-S/N HARPS solar spectrum as reference. We estimate a 5% typical error size in our HARPS flux measurements based on the possible differences in the continuum normalisation, which was performed in the same way for all the stars in the calibration. Therefore, the errors in SH are dominated by the conversion procedure.

Table 3

Stars included in the α calibration.

thumbnail Fig. 1

Linear fit between Santos et al. (2000)SMW values and the HARPS fluxes . A 5% error is estimated in our HARPS flux measurements.
Open with DEXTER

We proceed with the estimation of the activity index SH, with the corresponding indicators RHK (Middelkoop 1982) and (Noyes et al. 1984), which account for colour and photospheric correction, respectively.

4.2. Activity indicators

Figure 2 shows the cores of the CaII H and K lines in the HARPS normalised spectra of HD 1237 for different observations, which are compared with the mean profile (red) derived from the entire dataset. For the averaging procedure we take the slight differences in the wavelength range into account from each observation by an interpolation procedure to match the largest wavelength data points in the observed spectra. A high-S/N HARPS solar spectrum1 is shown as reference.

thumbnail Fig. 2

Core regions of the CaII K (left) and H (right) lines of HD 1237. Four spectra from our sample are plotted with the corresponding dates on the left y-axis and vertically shifted (0.25 units) for visualisation purposes. The black spectrum at the top was taken at a much later epoch (2012 Dec. 07). The red line shows the mean profile for the entire dataset, while the purple line is a HARPS solar spectrum used as reference.
Open with DEXTER

The upper plot shows one observation at a later epoch (2012 Dec 7), where it is possible to observe a variation in the line profile. This is interpreted as a slight change in the chromospheric and photospheric activity of the star in comparison with the mean behaviour of the red line (which is dominated by profiles from the first epoch), especially in the H line (Fig. 2, right panel). The K-line region of the spectrum contained more noise. No change in the activity level of the star is visible in this particular line.

Figure 3 shows the measured HARPS fluxes, (H + K)/(R + V), for each observation starting from 2012 Jul 15 (vertical blue line, HJD = 2 456 123.5). The x-axis units correspond to days after this initial date. The activity of the star showed a marginal variation, in a similar way as in the first epoch. This may be due to the rotation of active regions over the stellar surface. The similarities between the activity levels between both epochs could be an indication of a stable large-scale magnetic field configuration.

The chromospheric activity level of the star can be quantified by using Eq. (2) and the derived value for α in Eq. (3). We estimate an average value of SH ≃ 0.46 ± 0.02 for the available observations of HD 1237. Similar activity levels have been reported for the Sun-like star ξ boo A (Teff = 5600 K, Prot = 6.4 days, Age: ~0.2 Gyr, Mamajek & Hillenbrand 2008; Morgenthaler et al. 2012). As a reference value, the Solar S-index is S⊙ ≃ 0.1783 with a variation of ~0.02 from solar maximum to minimum (Lockwood et al. 2007).

Using the mean derived values of the S-index, we can now apply a transformation to obtain the parameter RHK, which takes the colour of the star in the activity estimation into account (Middelkoop 1982), Here RHK is defined as (4)where CCF is a colour-dependent function. For main sequence stars with 0.3 ≤ (B − V) ≤ 1.6, CCF is given by (5)Noyes et al. (1984) derived an expression in order to include photospheric corrections to the CaII core fluxes, , which is written as (6)with Rphot expressed also as a function of (B − V): (7)and valid in the range of 0.44 < (B − V) ≲ 1.0. Table 4 summarises the activity indicators for the available observations of HD 1237. The errors quoted are the mean measurement errors for these quantities. The B and V magnitudes were taken from the Hipparcos catalogue (B = 7.335 mag, V = 6.578 mag, Koen et al. 2010). High activity levels are related to magnetic fields on the stellar surface. These magnetic signatures are encoded in the polarised spectra of the star and is discussed in the next section.

thumbnail Fig. 3

HARPS fluxes, (H + K)/(R + V), for the available observations of HD 1237. The blue and green vertical lines denote the beginning of each observed epoch. The x-axis contains the number of days since the Heliocentric Julian Date (HJD = 2 456 123.5) of the first observation (2012 Jul. 15). The last two data points of the first epoch are unevenly distributed (after the vertical black line). We estimate a 5% error in our measurements.
Open with DEXTER

Table 4

Average activity indicators for HD 1237.

5. Magnetic field signatures

thumbnail Fig. 4

Calculated Bℓ for the available observations of HD 1237. The blue and green vertical lines denote the beginning of each observed epoch. The x-axis contains the number of days since the Heliocentric Julian Date (HJD = 2 456 123.5) of the first observation (2012 Jul. 15). The last two data points of the first epoch are unevenly distributed (after the black line). Colours indicate the line mask used for the LSD procedure.
Open with DEXTER

The signal-to-noise ratio (S/N) in the observations is not high enough to detect magnetically induced spectro-polarimetric signatures in single lines. However, by applying a multi-line technique, e.g. least squares deconvolution (LSD, Donati et al. 1997), it is possible to increase the S/N by a factor of ~50 − 100, adding-up the signal from thousands of spectral lines over the entire spectral range (for HARPS: 378–691 nm), enhancing our sensitivity to magnetic signatures in the observations (see Kochukhov et al. 2010 for a recent review of the LSD technique). This procedure requires a photospheric model (line list) matching the spectral type of our target star. This is done using an atomic line list database2 (Kupka et al. 2000) and the stellar fundamental properties listed in Table 1. We assumed a micro-turbulence parameter of 1.3 km s-1 (Ghezzi et al. 2010) and solar abundance for the photospheric line list that included ~15 000 lines within the HARPS spectral range.

From this initial photospheric line list, we generated two different masks used with the LSD calculation. In the first mask the strong lines (and lines blended with these lines) that form in the chromosphere or that break the basic assumptions of LSD (e.g. Ca II H&K, Hα) are removed (cleaned mask). This reduced the number of lines included to ~11 000. After the mask cleaning, a numerical routine based on the Levenberg-Marquardt, non-linear least-squares algorithm from the mpfit library (Moré 1978; Markwardt 2009) is applied to fit the line mask to the observed Stokes I spectrum after adjusting the individual depths of the spectral lines (cleaned-tweaked mask). This is performed through the entire HARPS wavelength coverage. This step is more commonly carried out when applying LSD to hot (OB-type) stars that have fewer lines (e.g. Neiner et al. 2012). Finally, LSD was applied to the spectro-polarimetric data using the final masks, generating in this way a single, averaged line profile per observation (LSD Stokes I, V, and diagnostic N profiles). A velocity step Δv = 1.4 km s-1 was used to construct the LSD profiles. This velocity spacing considers two pixels per spectral element of the instrument (in the case of HARPSpol, R = 2.5 km s-1 and 3.4 px per resolution element).

Both procedures lead us to consistent results in the obtained LSD signatures of the star. While no clear change is observed in the null polarisation check, subtle qualitative differences appear in both Stokes profiles, in the sense that the cleaned-tweaked mask seems to get a broader unpolarised profile with a slightly weaker signature in the circular polarised profile, in comparison with the clean line mask. These minimal shape differences in the LSD line profiles can have a much greater effect in hot stars where fewer lines are generally available. On the other hand, the resulting S/N of the LSD profiles for each individual observation was systematically higher using the cleaned-tweaked mask than in the cleaned case, despite the same number of spectral lines in their masks (11 048). On average, a ~7% increase was obtained in the S/N of the LSD profiles with the cleaned-tweaked mask.

5.1. Longitudinal magnetic field

thumbnail Fig. 5

Left: minimisation results for HD 1237 showing the reduced χ2 as a function of Prot. The black segmented lines corresponds to a fourth-order polynomial fit. Right: power spectrum obtained from the Lomb-Scargle periodogram analysis on the temporal variation of Bℓ (Fig. 4).
Open with DEXTER

Using the derived Stokes I and V LSD profiles, it is possible to obtain information about the surface averaged longitudinal magnetic field (Bℓ). Following Donati & Landstreet (2009), an estimate of Bℓ (in G) is given by (8)where the radial velocity shift v (in km s-1) is measured with respect to the average line derived from LSD with central wavelength λ (in μm) and mean Landé factor . As this measurement is an integrated quantity over the visible surface, it cannot provide complete information for stars that host complex large-scale magnetic fields. From multiple measurements of Bℓ taken over a stellar rotation period, it is possible to gain a first insight into inhomogeneities of the disk-integrated magnetic field. It is thus also possible to estimate the stellar rotation period, using the modulation in a time series of Bℓ measurements, provided that they span more than one rotation period.

Figure 4 shows the measurements of Bℓ for the available observations of HD 1237. The integration was centred on the radial velocity of the star (−5.2 km s-1) and covered the entire Stokes V signature (± 12.5 km s-1). The uncertainties were computed from standard error propagation from the spectra. The colours denote the line mask used in the LSD procedure, showing that the slight differences in each of the measurements are consistent within the errors. This implies that either set of LSD profiles can be used to obtain robust longitudinal magnetic field measurements. As the two sets of profiles are so similar, there is only a difference in S/N without any noticeable impact on the maps structure.

Reflecting the additive nature of the chromospheric emission, Bℓ shows a clearer rotational variation in comparison with SH (Fig. 3), in contrast to the mixed polarity effects of the magnetic field. The behaviour of the longitudinal field is consistent throughout the observed epochs with an varying amplitude of roughly ~10 G. This value is somewhat higher than the solar value (Bℓ< 4 G, Daou et al. 2006; Kotov et al. 1998) and than the average value from snapshot observations of other Sun-like stars of the same and different spectral types (3.3, 3.2, and 5.7 G for F, G and K-dwarfs respectively, Marsden et al. 2014). However, caution is advised in these averaged comparisons, given the rotational variability of Bℓ and the nature of a snapshot survey. Similar variations have been observed in the long-term monitoring of the active Sun-like star ξ Boo A (Morgenthaler et al. 2012).

6. Surface magnetic field mapping

6.1. Optimal line profile and stellar parameters

We reconstructed surface magnetic field maps by applying the tomographic inversion technique of Zeeman Doppler imaging (ZDI; Vogt et al. 1987; Semel 1989). As described by Hussain et al. (2000), ZDI has been used to recover magnetic field maps on the surfaces of stars ranging from T Tauri stars to binary systems (e.g. Barnes et al. 2004; Dunstone et al. 2008). The code recovers the magnetic flux distribution across the stellar disk, modulated by the stellar rotation, by using time series of photospheric absorption line profiles (LSD Stokes I) and circularly polarised profiles (LSD Stokes V).

To recover reliable magnetic field maps, it is necessary to properly model the LSD Stokes I and V profiles and their temporal variation. To do this, good constraints should be obtained on the local line profile description and the stellar parameters. Two different synthetic line shapes were tested: a Gaussian profile, which is commonly used in magnetic field studies in Sun-like stars (e.g. Boro Saikia et al. 2015), and a Milne-Eddington profile, fitted to a solar LSD profile derived from a high-S/N HARPS spectrum. This last approach was previously considered in ZDI of accreting T Tauri stars (Donati et al. 2008a) and M dwarfs (Morin et al. 2008). For both cases, we assumed a linear dependence of the continuum limb darkening with the cosine of the limb angle (slope u ≃ 0.65, Sing 2010).

thumbnail Fig. 6

Optimal χ2 selection criteria applied to the ZDI solution set of HD 1237. Each panel contains the results for the July dataset using the Gaussian (left) and Milne-Eddington (right) line profiles. The red symbols show the behaviour of the total entropy content ST as a function of the reduced χ2. Each point corresponds to a converged ZDI solution. Green and blue symbols represent the first and second derivatives as indicated. A fourth-order polynomial fit has been applied (segmented line) to find the value for which the rate of change in the information growth (second derivative) in the ZDI solution set is maximised. The optimal fit level is indicated in each case.
Open with DEXTER

In addition, we estimated the rotational period Prot differential rotation profile (e.g. Ω(l) = Ωeq − dΩsin2(l), see Petit et al. 2002) and inclination angle i of the star. This is done by generating a grid of ZDI models that covers a range of values for the involved quantities and minimising the reduced χ2 from synthetic line profile fitting (see Collier Cameron 1995; Hussain et al. 2009). Figure 5 (left) shows the results of the minimisation analysis over Prot. As mentioned in Sect. 2, the rotation period of HD 1237 is not well known. Our analysis shows the harmonic behaviour of this parameter with a fundamental value of ~7.0 days. We estimate a 10% error for this period determination, given the width of the minima in the left-hand panel of Fig. 4. This estimate is consistent with the high activity level measured in this star (Sect. 4). For completeness, we performed a Lomb-Scargle periodogram, using the first epoch dataset, on the temporal variations of Bℓ (Fig. 4). Figure 5 (right) shows the obtained power spectrum, having a best-fit period of Prot = 6.8 ± 0.2 days with an associated P-value statistic of 5.18 × 10-5 (Zechmeister & Kürster 2009).

No good constraints were obtained by χ2 minimisations for the differential rotation parameters and the inclination angle of the star. Therefore, no differential rotation profile was included in the mapping procedure. For the inclination angle, we considered the expected value from rigid body rotation, i.e. sini = (Prot·vsini) / (2πR∗). Given the uncertainties of the involved quantities (see Table 1), the inclination angle should lie somewhere between ~40–60°. We performed reconstructions for 40, 50, and 60° and find no substantial differences between these magnetic field reconstructions. We present here the maps obtained assuming a 50° inclination angle.

Finally, we compared the optimum stellar and line parameters associated with the LSD profiles produced using both masks (cleaned mask and clean-tweaked mask). No significant difference is found in the optimal parameters. In the subsequent analysis, therefore, we only consider the LSD profiles produced using the cleaned-tweaked mask.

6.2. Optimal fit quality: entropy content in ZDI maps

The last of part of the analysis corresponds to the selection of the optimal fit quality (reduced χ2), of the model with respect to the observations. In principle, the recovered field distribution and associated profiles should try to achieve the lowest possible value of χ2. Still, given the limitations of the observations and the ZDI technique (S/N, spatial resolution, phase-coverage, etc.), the goodness-of-fit level has to be determined carefully to avoid the appearance of numerical artefacts in the final maps. This is particularly important in the case of resolution-limited maps (i.e. slowly rotating stars). However, there are only a few published procedures for establishing a robust stopping criterion, i.e. the point at which noise starts to affect the reconstructed image. Motivated by this, we propose a systematic method for estimating the optimal fit quality for a given set of ZDI magnetic field maps, using HD 1237 as a test-case. It is important to note here that since this procedure is defined a posteriori over the resulting maps themselves (2D images), it does not modify the regularisation functions imposed to ZDI. In this sense, its application to other stellar systems should be straightforward.

We begin with the definition of the entropy S, as an estimate of the information content in an image3. Following Sonka et al. (2007), let P(k) be the probability that the difference between two adjacent pixels is equal to k. The image entropy can be estimated as (9)where log 2 is the base 2 logarithm. This implies that a larger or smaller amount of entropy in the image will depend on the contrast between adjacent pixels. An image that is perfectly constant will have an entropy of zero. For the methodology described below, we are not interested in the absolute values of the entropy, but rather its overall behaviour as a function of the reduced χ2.

For a given converged ZDI solution (i.e. a particular value of reduced χ2), we can calculate the total entropy content (ST) by applying the definition given by Eq. (9) to each one of the recovered maps: (10)where SR, SM, and SA represent the entropy contained in the radial, meridional, and azimuthal field components maps, respectively. Figure 6 shows in red the behaviour of ST as a function of the reduced χ2 for the ZDI solution set of HD 1237.

As expected, by decreasing the reduced χ2, the information content in the resulting ZDI solution increases (field strength and structure). The overall behaviour of the total entropy content and information growth is consistent for both cases. In the case of the Gaussian profile (Fig. 6, left panel), the entropy growth remains fairly constant (close to zero) for high reduced χ2 values, reaching a maximum4 around χ2 ≃ 1.2. However, by that point the concavity of the curve has changed (negative second derivative), suggesting a different regime for the information growth in the ZDI solutions. This is interpreted as a noise signature, reflected as artefacts in the final maps leading to an additional increment in the information growth. For this reason we adopt as optimal fit level the reduced χ2 value for which the rate of change in the information growth (second derivative) in the ZDI solution set is maximised. In this particular case, this occurs around χ2 ≃ 1.4, as indicated by the fourth-order polynomial fit (segmented line) in the left panel of Fig. 6. A similar criterion plot is constructed for the Milne-Eddington line profile (Fig. 6, right panel), where a lower optimal reduced χ2 ≃ 1.1 is achieved in this case5.

6.3. ZDI maps and synthetic Stokes V profiles

We reconstructed the ZDI surface magnetic field maps and the synthetic circularly polarised profiles based on the time series of LSD Stokes V spectra. For this we used the cleaned-tweaked line mask (Sect. 5) and the stellar and line parameters derived in Sect. 6.1. The goodness-of-fit level in each case is selected under the criterion described in the last section.

thumbnail Fig. 7

Results of the ZDI analysis for the first-epoch observations of HD 1237 using the Gaussian line profile. The first three panels show the surface magnetic field components BR, BM, and BA. The colour scale indicates the polarity and the magnitude of the magnetic field component in G, while the phase coverage is indicated by the black ticks in the upper y-axis. The segmented horizontal line indicates the surface visibility limit, imposed by the adopted inclination angle of the star (i = 50°). The last panel shows the comparison between synthetic (red) and observed (black) Stokes V profiles obtained for this particular epoch in each observational phase Φ, where the recovered maps fit the spectro-polarimetric data to an optimal reduced χ2 = 1.4.
Open with DEXTER

thumbnail Fig. 8

Results of the ZDI analysis for the first-epoch observations of HD 1237 using the Milne-Eddington line profile. See caption of Fig. 7 for more details. In this case, the maps fit the spectro-polarimetric data to an optimal reduced χ2 = 1.1.
Open with DEXTER

Figures 7 and 8 show the results of the ZDI procedure for the first-epoch observations (2012 July) of HD 1237. Four vertical panels are presented where the first three correspond to the Mercator-projected magnetic field maps in Gauss (G), decomposed into the radial (BR–top), meridional (BM–middle), and azimuthal (BA–bottom) components. The spatial resolution of the maps is ~18° in longitude at the stellar equator (at the poles, the map has a considerably poorer resolution than the equator). However, the phase coverage also has a significant impact on the level of detail that can be recovered. In this case, the right-hand side would have a slightly better resolution than the left-hand side of the image. The last panel shows the fitted synthetic Stokes V profiles to the spectro-polarimetric observations in each rotational phase (Φ).

The recovered maps show a relatively complex field distribution across the surface. The field is dominated by the azimuthal component, displaying a strong (~+90 G) large-scale ring-like structure around 45° in latitude. In the radial component, two large and moderately strong (~±50 G) regions of opposite polarities are also located at higher latitudes, while weaker (~±25 G) small-scale features of mixed polarities appear close to the equator. These large magnetic features are somewhat preserved with reversed polarities in the meridional magnetic field maps. However, some cross-talk from the radial component may be present in the meridional map (see Donati & Brown 1997).

Although the overall large-scale structure and field strength are consistent between both line profiles, several differences are clearly visible in the characteristics of the field components. First, the maps recovered using the Milne-Eddington line shape include additional small-scale features leading to a more complex field distribution in the surface. These are more prominent in the maps for the radial and meridional components. Second, the recovered magnetic field distribution, in the Milne-Eddington case, seems to be slightly shifted to lower latitudes. These differences can be understood from the fact that the Milne-Eddington line shape is a better representation of the derived LSD Stokes I profile. This is also true for the shape of the circularly polarised profile, leading to a more detailed (additional small-scale structures) and a field distribution shifted to lower latitudes (as a consequence of the sensitivity in the core of the profile). Both elements are translated into the lower optimal reduced χ2 value that can be reliably achieved in this case (see Sect. 6.2), and therefore an improved fit (lower panel in Fig. 8).

The reconstructed map for the second-epoch observations (2012 December) and the corresponding synthetic Stokes V profiles are presented in Appendix A.

7. Summary and discussion

In this paper, we have presented a detailed study that covers two observational epochs of the activity and magnetic field structure of the Sun-like planet-hosting star HD 1237. The chromospheric activity level of HD 1237 estimated from the calibrated SH = 0.46 ± 0.02 and values is similar to other active Sun-like stars (e.g. ξ boo A, Morgenthaler et al. 2012) and considerably higher than the solar case (~3 times higher in terms of the average S index, Lockwood et al. 2007). A much larger difference has been reported for the X-ray activity level, with a log (LX) = 29.02 ± 0.06 for HD 1237 (Kashyap et al. 2008), which is two orders of magnitude higher than the value estimated for the Sun during solar maximum (Peres et al. 2000). The chromospheric activity level remained fairly constant over the near five-month period between the two sets of observations. Our estimate falls between the previous measurements of (Naef et al. 2001) and (Saffe et al. 2005). Given the large uncertainties in the activity indicators, it is difficult to address whether these variations have some correspondence to a magnetic cycle in the star.

To extract magnetic field signatures from spectro-polarimetric data, we applied the LSD multi-line technique to the observations. Two different line masks for the LSD profiles were compared and tested. The standard procedure involves employing a mask “cleaned” of chromospheric and strong NLTE line profiles (e.g. Marsden et al. 2014); we compared the results from this procedure with improving the line list further by “tweaking” line strengths, so that they are tailored to the line depths in the observed spectrum. Both approaches lead to similar results in the obtained LSD profiles, and therefore in the physical quantities inferred (e.g. the longitudinal magnetic field, Bℓ). However, for the same number of spectral lines, the average S/N of the LSD profiles recovered with the aid of the cleaned-tweaked line mask was ~7% higher than in the cleaned case.

The longitudinal magnetic field (Bℓ), estimated from the Stokes I and V LSD profiles, showed a clear rotational modulation with an amplitude of ~10 G. This behaviour was preserved in both observed epochs. Placing these measurements in the context of other stars, similar variations have been observed in the long-term monitoring of the active Sun-like star ξ Boo A (~4–9 G, Age: ~0.2 Gyr, Mamajek & Hillenbrand 2008; Morgenthaler et al. 2012), the K-dwarf exoplanet host ϵ Eri (~10–12 G, Age: ~0.2–0.8 Gyr, Janson et al. 2008; Jeffers et al. 2014b), and more recently for the young solar analogue HN Peg (~14 G, Age: ~0.2 Gyr, Eisenbeiss et al. 2013; Boro Saikia et al. 2015). The chromospheric and X-ray activity levels of these stars are also very similar to HD 1237. The source ξ Boo A has an average S-index of SHK ≃ 0.45 and log (LX) ≃ 28.91 (Gray et al. 1996; Wood & Linsky 2010), and ϵ Eri has strong magnetic activity with a mean SHK ≃ 0.50 and log (LX) ≃ 28.22 (Jeffers et al. 2014b; Poppenhaeger et al. 2011). Similarly, previous reports for HN Peg show SHK ≃ 0.35 and log (LX) ≃ 29.19 (Boro Saikia et al. 2015; Schmitt & Liefke 2004). Despite the various similarities among these systems, the complete relation between the magnetic field and its influence over different layers of the stellar atmosphere (activity) have not yet been fully understood. Similar to these other systems, the chromospheric activity of HD 1237 does not show a clear correlation with Bℓ. This is interpreted as the result of probing different spatial and temporal scales in each of these measured quantities.

We recovered the optimal stellar parameters of HD 1237 using a reduced χ2 minimisation scheme, based on the tomographic inversion technique of ZDI. The analysis yields a rotation period, Prot = 7.0 ± 0.7 days. Literature values of Prot are uncertain, ranging between ~4.1 − 12.6 days (see Watson et al. 2010 and references therein). The 7.0-day value found in this work is consistent with the ZDI analysis and both the chromospheric and coronal activity levels of the star. Similar procedures were applied to estimate the inclination angle and the differential rotation of HD 1237. However, the available observations did not provide enough constraints for a robust determination of these parameters. Therefore, we considered an inclination angle derived using the stellar properties of the star and solid body rotation (i.e. i ~ 50°). No differential rotation profile was included in the reconstruction of the magnetic field maps.

For the surface magnetic field mapping procedure, two different synthetic line shapes (Gaussian/Milne-Eddington) were tested in order to investigate their impact on the ZDI maps for slowly rotating solar-type stars. We showed that both profiles recover robust magnetic field maps. However, the Milne-Eddington line profile yields a better spectro-polarimetric fit (lower optimal reduced χ2), leading to a more detailed structure recovered in the ZDI maps compared to the Gaussian case. In connection to this, a fit-stopping criterion based on the information content (entropy) of the ZDI maps solution set was introduced. This allows identification of the optimal reduced χ2 value, avoiding to some extent the appearance of numerical artefacts in the ZDI maps. The optimal fit level is given by the reduced χ2 value for which the rate of change of the entropy growth in the ZDI solution set is maximised.

The large-scale magnetic field of HD 1237 showed a complex distribution at the stellar surface. The strongest magnetic field features appear at middle latitudes (~45°) in the azimuthal and radial components. The field is dominated by the azimuthal component, displaying latitudinal belts or ring-like structures across the stellar surface. This has been observed in other studies of active Sun-like stars (e.g. Folsom et al. 2014) and in numerical simulations of a rapidly rotating Sun (Brown et al. 2010). As previously suggested, the appearance of significant or even dominant toroidal fields in the surface of these types of stars is connected with the rotation period, with ~12 days as the rotation threshold (Petit et al. 2008).

In comparison with other cool stars, HD 1237 (Prot = 7.0 days, this work, M∗ = 1.0 ± 0.1M⊙, Ghezzi et al. 2010) could be located near ξ Boo A (G8V, Age: ~0.2 Gyr) and GJ 182 (M0V, Age: Gyr, Liu et al. 2004) in the mass-period, large-scale magnetic field diagram (Fig. 3 in Donati & Landstreet 2009) with a mostly toroidal field topology and a considerable contribution from the radial component. With the same spectral type as HD 1237, ξ Boo A displays a slightly weaker field in its surface (~±60 G) with an alternating dominance between the radial and the azimuthal components in the observed long-trend evolution (~4 years, Morgenthaler et al. 2012). During the observed azimuthal-dominated epochs of this star, the field distribution is very similar to the one derived for HD 1237 in this work. A strong uni-directional azimuthal field appears at low latitudes with large mixed polarity regions in the radial field and a minor contribution from the meridional component. However, the magnetic regions of ξ Boo A are much larger and less concentrated than the ones found for HD 1237. This could be related to the small difference in their rotation periods (6.43 days for ξ Boo A and 7.0 days for HD 1237) and/or with the differential rotation that may be occurring in the surface. A similar situation appears in the case of GJ 182 in terms of the recovered magnetic field topology. However, the field strength is much larger in this case (up to 400 G, Donati et al. 2008b), which is mostly connected with the rapid rotation, large latitudinal shear, and the nearly fully convective nature of this star.

Regardless of the low vsini of HD 1237, it is clear that a more complex field distribution is required to fit the observed spectra. Some weaker small-scale regions are recovered closer to the stellar equator. Despite their relative low strength compared to the total surface field, the contribution from these small-scale features to magnetically-related phenomena may be significantly higher. In combination with strong mixed-polarity regions missed by ZDI, these features can influence the quiescent coronal emission (Johnstone et al. 2010), the X-ray modulation (Arzoumanian et al. 2011), and the wind structure around very active stars (Lang et al. 2014) and the Sun (Garraffo et al. 2013).

The large-scale field distribution appears similar at both of the observed epochs, confirming the observed chromospheric activity and longitudinal magnetic field behaviour. Using the final ZDI maps and the derived synthetic Stokes V profiles, we were able to reproduce the variations in Bℓ consistently (for both observed epochs). As an example, the maximum value of Bℓ in Fig. 4 coincides when the large positive polarity region of the radial field is located close to the limb (Day 3). Similarly, the minimum value obtained for Bℓ in Fig. 4 results from the negative polarity region in the radial field map, this time located much closer to the disk centre. This corroborates the robustness of our ZDI magnetic maps for this system.

As mentioned in Sect. 2, HD 1237 has a planetary companion with a mean separation of 0.49 AU and a projected mass of Mpsin(i) = 3.37 ± 0.14M♃. This exoplanet is relatively far out in comparison with other similar systems where ZDI maps of the host star are available (see Fares et al. 2013 and references therein). However, the enhanced activity levels of HD 1237 and the relatively strong and complex surface magnetic field (in connection with its stellar wind) could affect the conditions experienced by the exoplanet through its orbit significantly. Previous parametric studies of this system have predicted a relatively high mass loss rate (~85 ), which could even lead to magnetospheric radio emission from the exoplanet (Stevens 2005). This will be considered in a future study, incorporating the recovered ZDI magnetic field maps into a detailed 3D magneto-hydrodynamics (MHD) code (BATS-R-US, Powell et al. 1999; Tóth et al. 2012), which was originally developed and validated for the solar wind and corona (e.g. Sokolov et al. 2013; van der Holst et al. 2014) and recently applied in the stellar context (e.g. Cohen & Drake 2014; Cohen et al. 2014).

8. Conclusions

We reconstructed magnetic field maps of the young planet-hosting G-type star, HD 1237 using the technique of Zeeman Doppler imaging. As part of this detailed spectro-polarimetric study, we have performed the following.

  • We find that assuming a Milne-Eddington approximation for thelocal line profile produces a better fit to the shape of the observedLSD profiles. This influences the number of magnetic structuresthat can be recovered in the magnetic field maps.

  • We propose a robust method for defining the stopping criterion in ZDI techniques. This allows one to choose the optimum degree of fit, beyond which the model ceases to provide a good fit to the observed dataset owing to the introduction of artefacts into the resulting image. We successfully applied this to two different datasets with vastly different degrees of phase sampling.

  • As part of our optimisation routines, we recovered a rotation period of 7.0 days. This is consistent with the measured chromospheric and coronal activity levels of the star (, this paper; log (LX) = 29.02 ± 0.06, Kashyap et al. 2008).

  • The magnetic field reconstructions for HD 1237 are dominated by a band of strong uni-directional azimuthal field at high latitudes, accompanied by a complex multi-polar radial field distribution. The largest magnetic regions show field strengths of ~90 and 50 G for the azimuthal and the radial field components, respectively.

  • We note that the field topology recovered for HD 1237 is fully commensurate with studies of ξ Boo A and GJ 182, the two other stars that are closest to it in the mass, period diagram (Donati & Landstreet 2009). Larger sample sizes are needed to confirm these trends. If confirmed, it will be possible to predict the global magnetic field topologies and therefore the extended environments of planet-hosting stars at various stages in their evolution based on these fundamental parameters.

We will address the influence of the different magnetic scales on the coronal structure and wind properties of HD 1237 in a follow-up paper, using the ZDI maps presented here as boundary conditions (see Cohen et al. 2010, 2011a). Possible star-planet interactions occurring in the system via transient (e.g. CME events, Cohen et al. 2011b) or quiescent phenomena (e.g. planetary radio emission, Stevens 2005; Vidotto et al. 2012) can also be considered in future work.

1

S/N: 347 at 550 nm, date: 2007 Apr 12 – low activity period.

2

http://vald.astro.uu.se/ – Vienna Atomic Line Database (VALD3).

3

A similar implementation of entropy is commonly used as a regularisation function in ZDI (see Piskunov & Kochukhov 2002).

4

The apparent negative sign in the first derivative is due to the reversed direction of the x-axis (reduced χ2).

5

The same criterion was applied to generate the magnetic field maps in the December dataset (Appendix A). An optimal reduced χ2 ≃ 0.6 was obtained in this case, as expected for a dataset with fewer constraints.

Acknowledgments

Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 089.D-0138.

References

Appendix A: Second-epoch dataset (2012 Dec)

We consider only the Milne-Eddington line profile for the 2012 Dec dataset. Figure A.1 contains the recovered ZDI maps for this epoch with the corresponding synthetic Stokes V profiles. Despite having a lower phase coverage in this case, we were able to recover robust magnetic field maps, fitting the spectro-polarimetric data up to an optimal reduced χ2 = 0.6 (see Sect. 6.2). This low value of reduced χ2 results as a consequence of the fewer constraints available for this dataset.

thumbnail Fig. A.1

Results of the ZDI analysis for the second-epoch observations of HD 1237 using the Milne-Eddington line profile. See caption of Fig. 7 for more details. In this case the maps fit the spectro-polarimetric data to an optimal reduced χ2 = 0.6. In this case, Φ = 0.0 is assigned to the observations acquired at HJD 2 456 265.0.
Open with DEXTER

The field distribution clearly resembles the one obtained for the July dataset with a large contribution from the azimuthal and radial components to the total field. The main large-scale magnetic features are preserved between both observed epochs. This is consistent with the behaviour shown by the activity indicators (Sect. 4.2, Fig. 3) and the longitudinal magnetic field (Sect. 5.1, Fig. 4) in the entire dataset. Some of the smaller-scale structure is not recovered, and the ring of azimuthal field is not clear in this case. These changes are expected from the number of observations and phase coverage in this epoch (e.g. Donati & Brown 1997). The slight shift in longitude is due to the initial phase selection in this case, where Φ = 0.0 is assigned to the observations acquired at HJD 2 456 265.0 (Table 2).

All Tables

Table 1

HD 1237 basic properties.

In the text
Table 2

Journal of observations.

In the text
Table 3

Stars included in the α calibration.

In the text
Table 4

Average activity indicators for HD 1237.

In the text

All Figures

thumbnail Fig. 1

Linear fit between Santos et al. (2000)SMW values and the HARPS fluxes . A 5% error is estimated in our HARPS flux measurements.
Open with DEXTER
In the text
thumbnail Fig. 2

Core regions of the CaII K (left) and H (right) lines of HD 1237. Four spectra from our sample are plotted with the corresponding dates on the left y-axis and vertically shifted (0.25 units) for visualisation purposes. The black spectrum at the top was taken at a much later epoch (2012 Dec. 07). The red line shows the mean profile for the entire dataset, while the purple line is a HARPS solar spectrum used as reference.
Open with DEXTER
In the text
thumbnail Fig. 3

HARPS fluxes, (H + K)/(R + V), for the available observations of HD 1237. The blue and green vertical lines denote the beginning of each observed epoch. The x-axis contains the number of days since the Heliocentric Julian Date (HJD = 2 456 123.5) of the first observation (2012 Jul. 15). The last two data points of the first epoch are unevenly distributed (after the vertical black line). We estimate a 5% error in our measurements.
Open with DEXTER
In the text
thumbnail Fig. 4

Calculated Bℓ for the available observations of HD 1237. The blue and green vertical lines denote the beginning of each observed epoch. The x-axis contains the number of days since the Heliocentric Julian Date (HJD = 2 456 123.5) of the first observation (2012 Jul. 15). The last two data points of the first epoch are unevenly distributed (after the black line). Colours indicate the line mask used for the LSD procedure.
Open with DEXTER
In the text
thumbnail Fig. 5

Left: minimisation results for HD 1237 showing the reduced χ2 as a function of Prot. The black segmented lines corresponds to a fourth-order polynomial fit. Right: power spectrum obtained from the Lomb-Scargle periodogram analysis on the temporal variation of Bℓ (Fig. 4).
Open with DEXTER
In the text
thumbnail Fig. 6

Optimal χ2 selection criteria applied to the ZDI solution set of HD 1237. Each panel contains the results for the July dataset using the Gaussian (left) and Milne-Eddington (right) line profiles. The red symbols show the behaviour of the total entropy content ST as a function of the reduced χ2. Each point corresponds to a converged ZDI solution. Green and blue symbols represent the first and second derivatives as indicated. A fourth-order polynomial fit has been applied (segmented line) to find the value for which the rate of change in the information growth (second derivative) in the ZDI solution set is maximised. The optimal fit level is indicated in each case.
Open with DEXTER
In the text
thumbnail Fig. 7

Results of the ZDI analysis for the first-epoch observations of HD 1237 using the Gaussian line profile. The first three panels show the surface magnetic field components BR, BM, and BA. The colour scale indicates the polarity and the magnitude of the magnetic field component in G, while the phase coverage is indicated by the black ticks in the upper y-axis. The segmented horizontal line indicates the surface visibility limit, imposed by the adopted inclination angle of the star (i = 50°). The last panel shows the comparison between synthetic (red) and observed (black) Stokes V profiles obtained for this particular epoch in each observational phase Φ, where the recovered maps fit the spectro-polarimetric data to an optimal reduced χ2 = 1.4.
Open with DEXTER
In the text
thumbnail Fig. 8

Results of the ZDI analysis for the first-epoch observations of HD 1237 using the Milne-Eddington line profile. See caption of Fig. 7 for more details. In this case, the maps fit the spectro-polarimetric data to an optimal reduced χ2 = 1.1.
Open with DEXTER
In the text
thumbnail Fig. A.1

Results of the ZDI analysis for the second-epoch observations of HD 1237 using the Milne-Eddington line profile. See caption of Fig. 7 for more details. In this case the maps fit the spectro-polarimetric data to an optimal reduced χ2 = 0.6. In this case, Φ = 0.0 is assigned to the observations acquired at HJD 2 456 265.0.
Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.