Showing 1–38 of 38 results for author: Hackman, T

Surface distribution of small-scale magnetic field on the active cool star LQ Hydrae

Authors: O. Kochukhov, T. Hackman, J. J. Lehtinen

Abstract: It is well known that small-scale magnetism dominates the surface magnetic field topologies of active late-type stars. However, little information is available on the spatial distribution of this key magnetic field component. Here, we take advantage of the recently developed magnetic field diagnostic procedure relying on the magnetic intensification of iron atomic lines in the optical. We extend t… ▽ More It is well known that small-scale magnetism dominates the surface magnetic field topologies of active late-type stars. However, little information is available on the spatial distribution of this key magnetic field component. Here, we take advantage of the recently developed magnetic field diagnostic procedure relying on the magnetic intensification of iron atomic lines in the optical. We extend this methodology from measuring a single average field strength value to simultaneous Doppler imaging reconstruction of the two-dimensional maps of temperature and magnetic field strength. We applied this novel surface mapping approach to two spectroscopic data sets of the young active Sun-like star LQ Hya. For both epochs, we found a fairly uniform field strength distribution, apart from a latitudinal trend of the field strength increasing from 1.5-2.0 kG at low latitudes to 3.0-3.5 kG, close to the rotational poles. This distribution of the small-scale field does not display a clear correlation with the locations of temperature spots or the global magnetic field structure reconstructed for the same epochs. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 7 pages, 6 figures; accepted for publication in A&A

From convective stellar dynamo simulations to Zeeman-Doppler images

Authors: Thomas Hackman, Oleg Kochukhov, Mariangela Viviani, Jörn Warnecke, Maarit J. Korpi-Lagg, Jyri J. Lehtinen

Abstract: Zeeman-Doppler imaging (ZDI) is used to reconstruct the surface magnetic field of late-type stars from high resolution spectropolarimetric observations. The results are usually described in terms of characteristics of the field topology, i.e. poloidality vs. toroidality and axi-symmetry vs. non-axisymmetry in addition to the field strength. We want to test how well these characteristics are preser… ▽ More Zeeman-Doppler imaging (ZDI) is used to reconstruct the surface magnetic field of late-type stars from high resolution spectropolarimetric observations. The results are usually described in terms of characteristics of the field topology, i.e. poloidality vs. toroidality and axi-symmetry vs. non-axisymmetry in addition to the field strength. We want to test how well these characteristics are preserved when applying the ZDI method on simulated data, i.e. how accurately the field topology is preserved and to what extent stellar parameters influence the reconstruction. We use published magnetic field data from direct numerical MHD simulations. These have variable rotation rates, and hence represent different levels of activity, of an otherwise Sun-like setup. Our ZDI reconstruction is based on spherical harmonics expansion. By comparing the original values to those of the reconstructed images, we study the ability to reconstruct the surface magnetic field in terms of various characteristics of the field. The main large-scale features are reasonably well recovered, but the strength of the recovered magnetic field is just a fraction of the original input. The quality of the reconstruction shows clear correlations with the data quality. Furthermore, there are some spurious dependencies between stellar parameters and the characteristics of the field. Our study uncovers some limits of ZDI. Firstly, the recovered field strength will generally be lower than the "real" value as smaller structures with opposite polarities will be blurred in the inversion. Secondly, the axisymmetry is overestimated. The poloidality vs. toroidality is better recovered. The reconstruction works better for a stronger field and faster rotation velocity. Still, the ZDI method works surprisingly well even for a weaker field and slow rotation, provided the data has a high signal-to-noise and good rotation phase coverage. △ Less

Submitted 22 January, 2024; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 21 figures, A&A, in press, abstract shortened

V889 Her: abrupt changes in the magnetic field or differential rotation?

Authors: Teemu Willamo, Thomas Hackman, Jyri J. Lehtinen, Maarit Korpi-Lagg, Oleg Kochukhov

Abstract: We have applied Zeeman-Doppler imaging (ZDI) to an extensive spectropolarimetric HARPSpol data set of the magnetically active young solar analogue V889 Her, covering 35 spectra obtained during six nights in May 2011. The data set allows us to study Stokes V profiles of the star at almost identical rotational phases, separated by one or more stellar rotations. We use these data to study if the line… ▽ More We have applied Zeeman-Doppler imaging (ZDI) to an extensive spectropolarimetric HARPSpol data set of the magnetically active young solar analogue V889 Her, covering 35 spectra obtained during six nights in May 2011. The data set allows us to study Stokes V profiles of the star at almost identical rotational phases, separated by one or more stellar rotations. We use these data to study if the line profiles evolve from one rotation to the next, and find that some evolution does indeed occur. We consider two possible explanations for this: abrupt changes in the large-scale magnetic field or differential rotation. We find it quite difficult to distinguish between the two alternatives using ZDI alone. A strong differential rotation could, however, explain the changes in the line profiles, so we conclude that it must be present, and the abrupt magnetic field evolution is left uncertain. Commonly, rapidly rotating stars are assumed to have only weak differential rotation. If the strong differential rotation of V889 Her is indeed present, as has been found in other studies as well, it could indicate that the theoretical and numerical results of differential rotation still need to be revised. The rapid changes that may have occurred in the magnetic field indicate that one should be quite cautious when interpreting ZDI maps constructed from data over long time intervals. △ Less

Submitted 14 July, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Zeeman-Doppler imaging of five young solar-type stars

Authors: T. Willamo, J. J. Lehtinen, T. Hackman, M. J. Käpylä, O. Kochukhov, S. V. Jeffers, H. Korhonen, S. C. Marsden

Abstract: The magnetic activity of the Sun changes with the solar cycle. Similar cycles are found in other stars as well, but their details are not known to a similar degree. Characterising stellar magnetic cycles is important for the understanding of the stellar and solar dynamos that are driving the magnetic activity. We present spectropolarimetric observations of five young, solar-type stars, and compare… ▽ More The magnetic activity of the Sun changes with the solar cycle. Similar cycles are found in other stars as well, but their details are not known to a similar degree. Characterising stellar magnetic cycles is important for the understanding of the stellar and solar dynamos that are driving the magnetic activity. We present spectropolarimetric observations of five young, solar-type stars, and compare them to previous observations, with the aim to identify and characterise stellar equivalents of the solar cycle. We use Zeeman-Doppler imaging (ZDI) to map the surface magnetic field and brightness of our targets. The magnetic field is decomposed into spherical harmonic expansions, from which we report the strengths of the axisymmetric vs. non-axisymmetric, and poloidal vs. toroidal components, and compare them to the Rossby numbers of the stars. We present five new ZDI-maps of young, solar-type stars from Dec 2017. Of special interest is the case of V1358 Ori, that has gone through a polarity reversal between our observations and earlier ones. A less evident polarity reversal might also have occurred in HD 35296. There is a preference for more axisymmetric field, and possibly more toroidal field, for the more active stars with lower Rossby number, but a larger sample should be studied to draw any strong conclusions from this. For most of the individual stars, the amounts of toroidal and poloidal field have stayed on similar levels as in earlier observations. We find evidence for a magnetic polarity reversal having occurred in V1358 Ori. \c{hi}1 Ori could be an interesting target for future observations, with a possible short magnetic cycle of a few years. The correlation between the brightness maps and the magnetic field is mostly poor, which could indicate the presence of small-scale magnetic features of different polarities, that cancel each other out, and are not resolved in our maps. △ Less

Submitted 10 May, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Journal ref: A&A 659, A71 (2022)

Hidden magnetic fields of young suns

Authors: O. Kochukhov, T. Hackman, J. J. Lehtinen, A. Wehrhahn

Abstract: Global magnetic fields of active solar-like stars are nowadays routinely detected with spectropolarimetric measurements and are mapped with Zeeman-Doppler imaging (ZDI). However, due to the cancellation of opposite field polarities, polarimetry captures only a tiny fraction of the magnetic flux and cannot assess the overall stellar surface magnetic field if it is dominated by a small-scale compone… ▽ More Global magnetic fields of active solar-like stars are nowadays routinely detected with spectropolarimetric measurements and are mapped with Zeeman-Doppler imaging (ZDI). However, due to the cancellation of opposite field polarities, polarimetry captures only a tiny fraction of the magnetic flux and cannot assess the overall stellar surface magnetic field if it is dominated by a small-scale component. Analysis of Zeeman broadening in high-resolution intensity spectra can reveal these hidden complex magnetic fields. Historically, there were very few attempts to obtain such measurements for G dwarf stars due to the difficulty of disentangling Zeeman effect from other broadening mechanisms affecting spectral lines. Here we developed a new magnetic field diagnostic method based on relative Zeeman intensification of optical atomic lines with different magnetic sensitivity. Using this technique we obtained 78 field strength measurements for 15 Sun-like stars, including some of the best-studied young solar twins. We find that the average magnetic field strength $Bf$ drops from 1.3-2.0 kG in stars younger than about 120 Myr to 0.2-0.8 kG in older stars. The mean field strength shows a clear correlation with the Rossby number and with the coronal and chromospheric emission indicators. Our results suggest that magnetic regions have roughly the same local field strength $B\approx3.2$ kG in all stars, with the filling factor $f$ of these regions systematically increasing with stellar activity. Comparing our results with the spectropolarimetric analyses of global magnetic fields in the same stars, we find that ZDI recovers about 1% of the total magnetic field energy in the most active stars. This figure drops to just 0.01% for the least active targets. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Comments: 25 pages, 26 figures; accepted for publication in A&A

Journal ref: A&A 635, A142 (2020)

Shapes of stellar activity cycles

Authors: T. Willamo, T. Hackman, J. J. Lehtinen, M. J. Käpylä, N. Olspert, M. Viviani, J. Warnecke

Abstract: Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other ste… ▽ More Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other stellar parameters. We also study the differences in the shape of the solar cycle, depending on what activity indicator is used. The observations are also compared to simulated activity cycles. We use the chromospheric Ca II H&K data from the Mount Wilson Observatory HK Project. From this data set we identify 47 individual cycles from 18 stars. We use the statistical skewness of a cycle as a measure of its asymmetry, and compare this to other stellar parameters. A similar analysis has been done to magnetic cycles extracted from direct numerical magnetohydrodynamic simulations of solar-type convection zones. The shape of the solar cycle (fast rise and slower decline) is common in other stars as well, although the Sun has particularly asymmetric cycles. Cycle-to-cycle variations are large, but the average shape of a cycle is still fairly well represented by a sinusoid. We find only slight correlations between the cycle asymmetry and other stellar parameters. There are large differences in the shape of the solar cycle, depending on what activity indicator is used. In the simulated cycles, there is a difference in the symmetry of global simulations covering the full longitudinal range, hence capable of exciting non-axisymmetric dynamo modes, versus wedge simulations covering a partial extent in longitude, where only axisymmetric modes are possible. The former produce preferentially positive skewness, while the latter a negative one. △ Less

Submitted 29 April, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: Accepted to A&A; 12 pages, 12 figures

Journal ref: A&A 638, A69 (2020)

Topological changes in the magnetic field of LQ Hya during an activity minimum

Authors: Jyri J. Lehtinen, Maarit J. Käpylä, Thomas Hackman, Oleg Kochukhov, Teemu Willamo, Stephen C. Marsden, Sandra V. Jeffers, Gregory W. Henry, Lauri Jetsu

Abstract: Previous studies have related surface temperature maps, obtained with the Doppler imaging (DI) technique, of LQ Hya with long-term photometry. We compare surface magnetic field maps, obtained with the Zeeman Doppler imaging (ZDI) technique, with contemporaneous photometry, with the aim of quantifying the star's magnetic cycle characteristics. We inverted Stokes IV spectropolarimetry into magnetic… ▽ More Previous studies have related surface temperature maps, obtained with the Doppler imaging (DI) technique, of LQ Hya with long-term photometry. We compare surface magnetic field maps, obtained with the Zeeman Doppler imaging (ZDI) technique, with contemporaneous photometry, with the aim of quantifying the star's magnetic cycle characteristics. We inverted Stokes IV spectropolarimetry into magnetic field and surface brightness maps using a tomographic inversion code that models high signal-to-noise ratio mean line profiles produced by the least squares deconvolution (LSD) technique. The magnetic field and surface brightness maps reveal similar patterns to previous DI and ZDI studies: non-axisymmetric polar magnetic field structure, void of fields at mid-latitudes, and a complex structure in the equatorial regions. There is a weak but clear tendency of the polar structures to be linked with strong radial field and the equatorial ones with the azimuthal. We find a polarity reversal in the radial field between 2016 and 2017 coincident with an activity minimum seen in the long-term photometry. The inverted field strengths cannot easily be related with the observed spottedness, but we find that they are partially connected with the retrieved field complexity. Comparing to global magnetoconvection models for rapidly rotating young Suns, this field topology and dominance of the poloidal field component could be explained by a turbulent dynamo, where differential rotation does not play a major role (so called alpha^2 Omega or alpha^2 dynamos), and axi- and non-axisymmetric modes are excited simultaneously. The complex equatorial magnetic field structure could arise from the twisted (helical) wreaths often seen in these simulations, while the polar feature would be connected to the mostly poloidal non-axisymmetric component having a smooth spatial structure. △ Less

Submitted 15 December, 2021; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: 13 pages, 8 figures, 7 tables

Journal ref: A&A 660, A141 (2022)

Differences in the solar cycle variability of simple and complex active regions during 1996-2018

Authors: Shabnam Nikbakhsh, Eija Tanskanen, Maarit Käpylä, Thomas Hackman

Abstract: Aims. Our aim is to examine the solar cycle variability of magnetically simple and complex active region. Methods. We studied simple ($α$ and $β$) and complex ($βγ$ and $βγδ$) active regions based on the Mount Wilson magnetic classification by applying our newly developed daily approach. We analyzed the daily number of the simple active regions (SARs) and compared that to the abundance of the comp… ▽ More Aims. Our aim is to examine the solar cycle variability of magnetically simple and complex active region. Methods. We studied simple ($α$ and $β$) and complex ($βγ$ and $βγδ$) active regions based on the Mount Wilson magnetic classification by applying our newly developed daily approach. We analyzed the daily number of the simple active regions (SARs) and compared that to the abundance of the complex active regions (CARs) over the entire solar cycle 23 and cycle 24 until December 2018. Results. We show that CARs evolve differently over the solar cycle from SARs. The time evolution of SARs and CARs on different hemispheres also shows differences, even though on average their latitudinal distributions are shown to be similar. The time evolution of SARs closely follows that of the sunspot number, and their maximum abundance was observed to occur during the early maximum phase, while that of the CARs was seen roughly two years later. We furthermore found that the peak of CARs was reached before the latitudinal width of the activity band starts to decease. Conclusions. Our results suggest that the active region formation process is a competition between the large-scale dynamo (LSD) and the small-scale dynamo (SSD) near the surface, the former varying cyclically and the latter being independent of the solar cycle. During solar maximum, LSD is dominant, giving a preference to SARs, while during the declining phase the relative role of SSD increases. Therefore, a preference for CARs is seen due to the influence of the SSD on the emerging flux. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 9 pages, 10 figures, A&A journal

Journal ref: A&A 629, A45 (2019)

Spot evolution on LQ Hya during 2006--2017: temperature maps based on SOFIN and FIES data

Authors: Elizabeth M. Cole-Kodikara, Maarit J. Käpylä, Jyri J. Lehtinen, Thomas Hackman, Ilya V. Ilyin, Nikolai Piskunov, Oleg Kochukhov

Abstract: Context. LQ Hya is one of the most frequently studied young solar analogue stars. Recently, it has been observed to show intriguing behaviour by analysing long-term photometry: during 2003--2009, a coherent spot structure migrating in the rotational frame has been reported by various authors, but since that time the star has entered a chaotic state where coherent structures seem to have disappeare… ▽ More Context. LQ Hya is one of the most frequently studied young solar analogue stars. Recently, it has been observed to show intriguing behaviour by analysing long-term photometry: during 2003--2009, a coherent spot structure migrating in the rotational frame has been reported by various authors, but since that time the star has entered a chaotic state where coherent structures seem to have disappeared and rapid phase jumps of the photometric minima occur irregularly over time. Aims. LQ Hya is one of the stars included in the SOFIN/FIES long-term monitoring campaign extending over 25 years. Here we publish new temperature maps for the star during 2006--2017, covering the chaotic state of the star. Methods. We use a Doppler imaging technique to derive surface temperature maps from high-resolution spectra. Results. From the mean temperatures of the Doppler maps we see a weak but systematic increase in the surface temperature of the star. This is consistent with the simultaneously increasing photometric magnitude. During nearly all observing seasons we see a high-latitude spot structure which is clearly nonaxisymmetric. The phase behaviour of this structure is very chaotic but agrees reasonably well with the photometry. Equatorial spots are also frequently seen, but many of them we interpret to be artefacts due to the poor to moderate phase coverage. Conclusions. Even during the chaotic phase of the star, the spot topology has remained very similar to the higher activity epochs with more coherent and long-lived spot structures: we see high-latitude and equatorial spot activity, the mid latitude range still being most often void of spots. We interpret the erratic jumps and drifts in phase of the photometric minima to be caused by changes in the high-latitude spot structure rather than the equatorial spots. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 9 pages, 5 figures

Journal ref: A&A 629, A120 (2019)

Starspot activity of HD 199178. Doppler images from 1994--2017

Authors: Thomas Hackman, Iya Ilyin, Jyri J. Lehtinen, Oleg Kochukhov, Maarit J. Käpylä, Nikolai Piskunov, Teemu Willamo

Abstract: Context. Studying the spots of late-type stars is crucial for distinguishing between the various proposed dynamo mechanisms believed to be the main cause of starspot activity. For this research it is important to collect observation time series that are long enough to unravel both long- and short-term spot evolution. Doppler imaging is a very efficient method for studying spots of stars that canno… ▽ More Context. Studying the spots of late-type stars is crucial for distinguishing between the various proposed dynamo mechanisms believed to be the main cause of starspot activity. For this research it is important to collect observation time series that are long enough to unravel both long- and short-term spot evolution. Doppler imaging is a very efficient method for studying spots of stars that cannot be angularly resolved. Aims. High-resolution spectral observations during 1994-2017 are analysed in order to reveal long- and short-term changes in the spot activity of the FK Comae-type subgiant HD 199178. Methods. Most of the observations were collected with the Nordic Optical Telescope. The Doppler imaging temperature maps were calculated using an inversion technique based on Tikhonov regularisation and utilising multiple spectral lines. Results. We present a unique series of 41 temperature maps spanning more than 23 years. All reliable images show a large cool spot region centred near the visible rotation pole. Some lower latitude cool features are also recovered, although the reliability of these is questionable. There is an expected anti-correlation between the mean surface temperature and the spot coverage. Using the Doppler images, we construct the equivalent of a solar butterfly diagram for HD 199178. Conclusions. HD 199178 clearly has a long-term large and cool spot structure at the rotational pole. This spot structure dominated the spot activity during the years 1994-2017. The size and position of the structure has evolved with time, with a gradual increase during the last years. The lack of lower latitude features prevents the determination of a possible differential rotation. △ Less

Submitted 7 June, 2019; v1 submitted 5 December, 2018; originally announced December 2018.

Comments: 11 pages, 7 figures, Astronomy and Astrophysics, parallel version of published article

Journal ref: A&A 625, A79 (2019)

Long-term spot monitoring of the young solar analogue V889 Her

Authors: T. Willamo, T. Hackman, J. J. Lehtinen, M. J. Käpylä, I. Ilyin, G. W. Henry, L. Jetsu, O. Kochukhov, N. Piskunov

Abstract: Starspots are important manifestations of stellar magnetic activity. By studying their behaviour in young solar analogues, we can unravel the properties of their magnetic cycles. This gives crucial information of the underlying dynamo process. Comparisons with the solar cycle enable us to infer knowledge about how the solar dynamo has evolved during the Sun's lifetime. Here we study the correlatio… ▽ More Starspots are important manifestations of stellar magnetic activity. By studying their behaviour in young solar analogues, we can unravel the properties of their magnetic cycles. This gives crucial information of the underlying dynamo process. Comparisons with the solar cycle enable us to infer knowledge about how the solar dynamo has evolved during the Sun's lifetime. Here we study the correlation between photometric brightness variations, spottedness and mean temperature in V889 Her, a young solar analogue. Our data covers 18 years of spectroscopic and 25 years of photometric observations. We use Doppler imaging to derive temperature maps from high-resolution spectra. We use the Continuous Period Search method to retrieve mean V-magnitudes from photometric data. Our Doppler imaging maps show a persistent polar spot structure varying in strength. This structure is centered slightly off the rotational pole. The mean temperature derived from the maps shows an overall decreasing trend, as does the photometric mean brightness, until it reaches its minimum around 2017. The filling factor of cool spots, however, shows only a weak tendency to anti-correlate with the decreasing mean brightness. We interpret V889 Her to have entered into a grand maximum in its activity. The clear relation between the mean temperature of the Doppler imaging surface maps and the mean magnitude supports the reliability of the Doppler images. The lack of correlation between the mean magnitude and the spottedness may indicate that bright features in the Doppler images are real. △ Less

Submitted 27 February, 2019; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: 9 pages, 7 figures, 4 tables

Journal ref: A&A 622, A170 (2019)

Contemporaneous Imaging Comparisons of the Spotted Giant $σ$ Geminorum Using Interferometric, Spectroscopic, and Photometric Data

Authors: Rachael M. Roettenbacher, John D. Monnier, Heidi Korhonen, Robert O. Harmon, Fabien Baron, Thomas Hackman, Gregory W. Henry, Gail H. Schaefer, Klaus G. Strassmeier, Michael Weber, Theo A. ten Brummelaar

Abstract: Nearby, active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary $σ$ Geminorum ($σ$ Gem) was imaged for two epochs of interferometric, h… ▽ More Nearby, active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary $σ$ Geminorum ($σ$ Gem) was imaged for two epochs of interferometric, high-resolution spectroscopic, and photometric observations. The light curves from the reconstructions show good agreement with the observed light curves, supported by the longitudinally-consistent spot features on the different maps. However, there is strong disagreement in the spot latitudes across the methods. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: Accepted to ApJ, 19 pages, 16 figures, 5 tables

Magnetic fields of young solar twins

Authors: L. Rosén, O. Kochukhov, T. Hackman, J. Lehtinen

Abstract: The goal of this work is to study the magnetic fields of six young solar-analogue stars both individually, and collectively, to search for possible magnetic field trends with age. If such trends are found, they can be used to understand magnetism in the context of stellar evolution of solar-like stars and, the past of the Sun and the solar system. This is also important for the atmospheric evoluti… ▽ More The goal of this work is to study the magnetic fields of six young solar-analogue stars both individually, and collectively, to search for possible magnetic field trends with age. If such trends are found, they can be used to understand magnetism in the context of stellar evolution of solar-like stars and, the past of the Sun and the solar system. This is also important for the atmospheric evolution of the inner planets, Earth in particular. We used Stokes IV data from two different spectropolarimeters, NARVAL and HARPSpol. The least-squares deconvolution multi-line technique was used to increase the signal-to-noise ratio of the data. We then applied a modern Zeeman-Doppler imaging code in order to reconstruct the magnetic topology of all stars and the brightness distribution of one of our studied stars. Our results show a significant decrease in the magnetic field strength and energy as the stellar age increases from 100Myr to 250Myr while there is no significant age dependence of the mean magnetic field strength for stars with ages 250-650Myr. The spread in the mean field strength between different stars is comparable to the scatter between different observations of individual stars. The meridional field component has the weakest strength compared to the radial and azimuthal field components in 15 out of the 16 magnetic maps. It turns out that 89-97% of the magnetic field energy is contained in l=1-3. There is also no clear trend with age and distribution of field energy into poloidal/toroidal and axisymmetric/non-axisymmetric components within the sample. The two oldest stars in this study do show a twice as strong octupole component compared to the quadrupole component. This is only seen in one out of 13 maps of the younger stars. One star, chi1 Ori displays two field polarity switches during almost 5 years of observations suggesting a magnetic cycle length of either 2, 6 or 8 years. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 25 pages, 15 figures, 4 tables; Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 593, A35 (2016)

Time series analysis of long-term photometry of BM Canum Venaticorum

Authors: L. Siltala, L. Jetsu, T. Hackman, G. W. Henry, L. Immonen, P. Kajatkari, J. Lankinen, J. Lehtinen, S. Monira, S. Nikbakhsh, A. Viitanen, J. Viuho, T. Willamo

Abstract: Studying RS CVn binaries is challenging, because in addition to spot activity, other effects such as mass transfer between the components and gravitational distortion of their spherical forms may distort their light curves. Such effects can, however, be removed from the data by subtracting a mean light curve phased with the orbital period. We study a quarter of a century of standard Johnson differ… ▽ More Studying RS CVn binaries is challenging, because in addition to spot activity, other effects such as mass transfer between the components and gravitational distortion of their spherical forms may distort their light curves. Such effects can, however, be removed from the data by subtracting a mean light curve phased with the orbital period. We study a quarter of a century of standard Johnson differential V photometry of the RS CVn binary BM CVn. Our main aims are to determine the activity cycles, the rate of surface differential rotation and the rotation period of the active longitudes of BM CVn. The Continuous Period Search (CPS) is applied to the photometry. The changes of the mean and amplitude of the light curves are used to search for activity cycles. The rotation period changes give an estimate of the rate of surface differential rotation. The Kuiper method is applied to the epochs of the primary and secondary minima to search for active longitudes. The photometry reveals the presence of a stable mean light curve (MLC) connected to the orbital period P_orb = 20d.6252 of this binary. We remove this MLC from the original V magnitudes which gives us the corrected V' magnitudes. These two samples of V and V' data are analysed separately with CPS. The fraction of unreliable CPS models decreases when the MLC is removed. The same significant activity cycle of approximately 12.5 years is detected in both V and V' samples. The estimate for the surface differential rotation coefficient, k >= 0.10, is the same for both samples, but the number of unrealistic period estimates decreases after removing the MLC. The same active longitude period of P_al = 20d.511 +- 0d.005 is detected in the V and V' magnitudes. This long-term regularity in the epochs of primary and secondary minima of the light curves is not caused by the MLC. On the contrary, the MLC hampers the detection of active longitudes. △ Less

Submitted 13 September, 2016; v1 submitted 14 January, 2016; originally announced January 2016.

Comments: 11 pages, 13 figures

Activity trends in young solar-type stars

Authors: Jyri Lehtinen, Lauri Jetsu, Thomas Hackman, Perttu Kajatkari, Gregory W. Henry

Abstract: We study a sample of 21 young and active solar-type stars with spectral types ranging from late F to mid K and characterize the behaviour of their activity. We apply the continuous period search (CPS) time series analysis method on 16 to 27 years of photometry to estimate the surface differential rotation and determine the existence and behaviour of active longitudes and activity cycles on the sta… ▽ More We study a sample of 21 young and active solar-type stars with spectral types ranging from late F to mid K and characterize the behaviour of their activity. We apply the continuous period search (CPS) time series analysis method on 16 to 27 years of photometry to estimate the surface differential rotation and determine the existence and behaviour of active longitudes and activity cycles on the stars. We supplement the time series results by calculating new $\log{R'_{\rm HK}}$ emission indices for the stars from high resolution spectroscopy. The photometric rotation period variations reveal a positive correlation between the relative differential rotation coefficient and the rotation period, but do not reveal any dependence of the differential rotation on the effective temperature of the stars. Secondary period searches reveal activity cycles in 18 of the stars and temporary or persistent active longitudes in 11 of them. The activity cycles fall into specific activity branches. We find a new split into sub-branches, indicating multiple simultaneously present cycle modes. Active longitudes appear to be present only on the more active stars. There is a sharp break at approximately $\log{R'_{\rm HK}}=-4.46$ separating the less active stars with long-term axisymmetric spot distributions from the more active ones with non-axisymmetric configurations. In seven out of eleven of our stars with clearly detected long-term non-axisymmetric spot activity the estimated active longitude periods are significantly shorter than the mean photometric rotation periods. This systematic trend can be interpreted either as a sign of the active longitudes being sustained from a deeper level in the stellar interior than the individual spots or as azimuthal dynamo waves exhibiting prograde propagation. △ Less

Submitted 17 March, 2016; v1 submitted 22 September, 2015; originally announced September 2015.

Comments: 24 pages, 14 figures

Journal ref: A&A 588, A38 (2016)

Zeeman-Doppler imaging of active young solar type stars

Authors: Thomas Hackman, Jyri Lehtinen, Lisa Rosén, Oleg Kochukhov, Maarit J. Käpylä

Abstract: By studying young magnetically active late-type stars, i.e. analogues to the young Sun, one can draw conclusions on the evolution of the solar dynamo. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric spots in three young late-type stars. High-resolution spectropolarimetry of the targets were obtained with the HARPSpol i… ▽ More By studying young magnetically active late-type stars, i.e. analogues to the young Sun, one can draw conclusions on the evolution of the solar dynamo. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric spots in three young late-type stars. High-resolution spectropolarimetry of the targets were obtained with the HARPSpol instrument mounted at the ESO 3.6 m telescope. The signal-to-noise ratio of the Stokes IV measurements were boosted by combining the signal from a large number of spectroscopic absorption lines through the least squares deconvolution technique. Surface brightness and magnetic field maps were calculated using the Zeeman-Doppler imaging technique. All the three targets show clear signs of both magnetic fields and cool spots. Only one of the targets, namely V1358 Ori, shows evidence of the dominance of non-axisymmetric modes. In two of the targets, the poloidal field is significantly stronger than the toroidal one, indicative of an $α^2$-type of a dynamo, in which convective turbulence effects dominate over the weak differential rotation. In two of the cases there is a slight anti-correlation between the cool spots and the strength of the radial magnetic field. However, even in these cases the correlation is much weaker than in the case of sunspots. The weak correlation between the measured radial magnetic field and cool spots may indicate a more complex magnetic field structure in the spots or spot groups involving mixed magnetic polarities. Comparison with a previously published magnetic field map shows that on one of the stars, HD 29615, the underlying magnetic field has changed its polarity between 2009 and 2013. △ Less

Submitted 1 March, 2016; v1 submitted 8 September, 2015; originally announced September 2015.

Comments: 7 pages, 6 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 587, A28 (2016)

Study of FK Comae Berenices: VII. Correlating photospheric and chromospheric activity

Authors: K. Vida, H. Korhonen, I. V. Ilyin, K. Oláh, M. I. Andersen, T. Hackman

Abstract: We study the connection between the chromospheric and photospheric behaviour of the active late-type star FK Comae. We use spot temperature modelling, light curve inversion based on narrow- and wide-band photometric measurements, Halpha observations from 1997-2010, and Doppler maps from 2004-2010 to compare the behaviour of chromospheric and photospheric features. Investigating low-resolution Halp… ▽ More We study the connection between the chromospheric and photospheric behaviour of the active late-type star FK Comae. We use spot temperature modelling, light curve inversion based on narrow- and wide-band photometric measurements, Halpha observations from 1997-2010, and Doppler maps from 2004-2010 to compare the behaviour of chromospheric and photospheric features. Investigating low-resolution Halpha spectra we find that the changes in the chromosphere seem to happen mainly on a time scale longer than a few hours, but shorter variations were also observed. According to the Halpha measurements prominences are often found in the chromosphere that reach to more than a stellar radius and are stable for weeks, and which seem to be often, but not every time connected with dark photospheric spots. The rotational modulation of the Halpha emission seems to typically be anticorrelated with the light curve, but we did not find convincing evidence of a clear connection in the long-term trends of the Halpha emission and the brightness of the star. In addition, FK Com seems to be in an unusually quiet state in 2009-2010 with very little chromospheric activity and low spot contrast, that might indicate the long-term decrease of activity. △ Less

Submitted 26 May, 2015; originally announced May 2015.

Comments: 14 pages, 14 figures, 1 table; accepted for publication in A&A

Journal ref: A&A 580, A64 (2015)

Doppler Imaging of LQ Hya for 1998-2002

Authors: E. M. Cole, T. Hackman, M. J. Käpylä, I. Ilyin, O. Kochukhov, N. Piskunov

Abstract: We study the spot distribution on the surface of LQ~Hya during the observing seasons October 1998 -- November 2002. We look for persistent active longitudes, trends in the level of spot activity and compare to photometric data, specifically to the derived time epochs of the lightcurve minima. We apply the Doppler Imaging technique on photospheric spectral lines using an inversion code to retrieve… ▽ More We study the spot distribution on the surface of LQ~Hya during the observing seasons October 1998 -- November 2002. We look for persistent active longitudes, trends in the level of spot activity and compare to photometric data, specifically to the derived time epochs of the lightcurve minima. We apply the Doppler Imaging technique on photospheric spectral lines using an inversion code to retrieve images of the surface temperature. We present new temperature maps using multiple spectral lines for a total of 7 seasons. We calculate spot coverage fractions from each map, and as a result we find a general trend that is interpreted to be an indication of a spot cycle. There is a minimum during the observing season of March 1999. After this the activity increases until November 2000, followed by a general decrease in activity again. We find no evidence for active longitudes persisting over multiple observing seasons. The spot activity appears to be concentrated to two latitude regions. The high latitude spots are particularly strong when the spot coverage is at a maximum. Using the currently accepted rotation period, we find spot structures to show a trend in the phase-time plot, indicative of a need for a longer period. We conclude that the long-term activity of LQ~Hya is more chaotic than that of some magnetically active binary stars analyzed with similar methods, but still with clear indications of an activity cycle. △ Less

Submitted 14 April, 2015; originally announced April 2015.

Comments: 9 pages, 10 figures

Journal ref: A&A 581, A69 (2015)

Stellar activity as noise in exoplanet detection I. Methods and application to solar-like stars and activity cycles

Authors: H. Korhonen, J. M. Andersen, N. Piskunov, T. Hackman, D. Juncher, S. P. Jarvinen, U. G. Joergensen

Abstract: The detection of exoplanets using any method is prone to confusion due to the intrinsic variability of the host star. We investigate the effect of cool starspots on the detectability of the exoplanets around solar-like stars using the radial velocity method. For investigating this activity-caused "jitter" we calculate synthetic spectra using radiative transfer, known stellar atomic and molecular l… ▽ More The detection of exoplanets using any method is prone to confusion due to the intrinsic variability of the host star. We investigate the effect of cool starspots on the detectability of the exoplanets around solar-like stars using the radial velocity method. For investigating this activity-caused "jitter" we calculate synthetic spectra using radiative transfer, known stellar atomic and molecular lines, different surface spot configurations, and an added planetary signal. Here, the methods are described in detail, tested and compared to previously published studies. The methods are also applied to investigate the activity jitter in old and young solar-like stars, and over a solar-like activity cycles. We find that the mean full jitter amplitude obtained from the spot surfaces mimicking the solar activity varies during the cycle approximately between 1 m/s and 9 m/s. With a realistic observing frequency a Neptune mass planet on a one year orbit can be reliably recovered. On the other hand, the recovery of an Earth mass planet on a similar orbit is not feasible with high significance. The methods developed in this study have a great potential for doing statistical studies of planet detectability, and also for investigating the effect of stellar activity on recovered planetary parameters. △ Less

Submitted 5 January, 2015; originally announced January 2015.

Comments: Accepted to MNRAS

Doppler images and the underlying dynamo. The case of AF Leporis

Authors: S. P. Järvinen, R. Arlt, T. Hackman, S. C. Marsden, M. Küker, I. V. Ilyin, S. V. Berdyugina, K. G. Strassmeier, I. A. Waite

Abstract: The (Zeeman-)Doppler imaging studies of solar-type stars very often reveal large high-latitude spots. This also includes F stars that possess relatively shallow convection zones, indicating that the dynamo operating in these stars differs from the solar dynamo. We aim to determine whether mean-field dynamo models of late-F type dwarf stars can reproduce the surface features recovered in Doppler ma… ▽ More The (Zeeman-)Doppler imaging studies of solar-type stars very often reveal large high-latitude spots. This also includes F stars that possess relatively shallow convection zones, indicating that the dynamo operating in these stars differs from the solar dynamo. We aim to determine whether mean-field dynamo models of late-F type dwarf stars can reproduce the surface features recovered in Doppler maps. In particular, we wish to test whether the models can reproduce the high-latitude spots observed on some F dwarfs. The photometric inversions and the surface temperature maps of AF Lep were obtained using the Occamian-approach inversion technique. Low signal-to-noise spectroscopic data were improved by applying the least-squares deconvolution method. The locations of strong magnetic flux in the stellar tachocline as well as the surface fields obtained from mean-field dynamo solutions were compared with the observed surface temperature maps. The photometric record of AF Lep reveals both long- and short-term variability. However, the current data set is too short for cycle-length estimates. From the photometry, we have determined the rotation period of the star to be 0.9660+-0.0023 days. The surface temperature maps show a dominant, but evolving, high-latitude (around +65 degrees) spot. Detailed study of the photometry reveals that sometimes the spot coverage varies only marginally over a long time, and at other times it varies rapidly. Of a suite of dynamo models, the model with a radiative interior rotating as fast as the convection zone at the equator delivered the highest compatibility with the obtained Doppler images. △ Less

Submitted 9 December, 2014; originally announced December 2014.

Comments: accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 574, A25 (2015)

Multiperiodicity, modulations and flip-flops in variable star light curves III. Carrier fit analysis of LQ Hya photometry for 1982-2014

Authors: Nigul Olspert, Maarit J. Käpylä, Jaan Pelt, Elizabeth M. Cole, Thomas Hackman, Jyri Lehtinen, Gregory W. Henry

Abstract: We study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the perio… ▽ More We study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the period that is found to be the optimal one. The harmonic least-squares analysis of all the available data reveals a short period close to 1.6 days as a limiting value for a set of significant frequencies. We interpret this as the rotation period of the spots near the equatorial region. In addition, the distribution of the significant periods is found to be bimodal, hinting of a longer-term modulating period, which we set out to study with a two-harmonic CF model. Weak modulation signal is, indeed retrieved, with a period of roughly 6.9 years. The phase dispersion analysis gives a clear symmetric minimum for coherence times lower than and around 100 days. We interpret this as the mean rotation period of the spots (1.60514 days), and this value is chosen to be used as the carrier period for the CF analysis. With the CF method we seek for any systematic trends in the spot distribution in the global time frame, and locally look for abrupt phase changes earlier reported in rapidly rotating objects. During 2005-2008 the global CF reveals a coherent structure rotating with a period of 1.6037 days, while during most other times the spot distribution appears rather random in phase. The evolution of the spot distribution of the object is found to be very chaotic, with no clear signs of an azimuthal dynamo wave that would persist over longer time scales, although the short-lived coherent structures observed occasionally do not rotate with the same speed as the mean spot distribution. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: 14 pages, 10 figures, submitted to Astronomy and Astrophysics

Journal ref: A&A 577, A120 (2015)

Spot activity of the RS CVn star σ Geminorum

Authors: P. Kajatkari, T. Hackman, L. Jetsu, J. Lehtinen, G. W. Henry

Abstract: We model the photometry of RS CVn star $σ$ Geminorum to obtain new information on the changes of the surface starspot distribution, i.e., activity cycles, differential rotation and active longitudes. We use the previously published Continuous Periods Search-method (CPS) to analyse V-band differential photometry obtained between the years 1987 and 2010 with the T3 0.4 m Automated Telescope at the F… ▽ More We model the photometry of RS CVn star $σ$ Geminorum to obtain new information on the changes of the surface starspot distribution, i.e., activity cycles, differential rotation and active longitudes. We use the previously published Continuous Periods Search-method (CPS) to analyse V-band differential photometry obtained between the years 1987 and 2010 with the T3 0.4 m Automated Telescope at the Fairborn Observatory. The CPS-method divides data into short subsets and then models the light curves with Fourier-models of variable orders and provides estimates of the mean magnitude, amplitude, period and light curve minima. These light curve parameters are then analysed for signs of activity cycles, differential rotation and active longitudes. We confirm the presence of two previously found stable active longitudes, synchronised with the orbital period $P_{\rm{orb}}=19.60$d and find eight events where the active longitudes are disrupted. The epochs of the primary light curve minima rotate with a shorter period $P_{\rm{min,1}}=19.47$d than the orbital motion. If the variations in the photometric rotation period were to be caused by differential rotation, this would give a differential rotation coefficient of $α\ge 0.103$. The presence of two slightly different periods of active regions may indicate a superposition of two dynamo modes, one stationary in the orbital frame and the other one propagating in the azimuthal direction. Our estimate of the differential rotation is much higher than previous results. However, simulations show that this can be caused by insufficient sampling in our data. △ Less

Submitted 5 November, 2013; originally announced November 2013.

Comments: 10 pages, 6 figures. Submitted to A&A

Periodicity in some light curves of the solar analogue V352 CMa

Authors: P. Kajatkari, L. Jetsu, E. Cole, T. Hackman, G. W. Henry, S-L. Joutsiniemi, J. Lehtinen, V. Mäkelä, S. Porceddu, K. Ryynänen, V. Şolea

Abstract: Our aim was to study the light curve periodicity of the solar analogue V352 CMa and in particular show that the presence or absence of periodicity in low amplitude light curves can be modelled with the Continuous Period Search (CPS) method. We applied the CPS method to 14 years of V-band photometry of V352 CMa and obtained estimates for the mean, amplitude, period and minima of the light curves in… ▽ More Our aim was to study the light curve periodicity of the solar analogue V352 CMa and in particular show that the presence or absence of periodicity in low amplitude light curves can be modelled with the Continuous Period Search (CPS) method. We applied the CPS method to 14 years of V-band photometry of V352 CMa and obtained estimates for the mean, amplitude, period and minima of the light curves in the selected datasets. We also applied the Power Spectrum Method (PSM) to these datasets and compared the performance of this frequently applied method to that of CPS. We detected signs of a $11.7 \pm 0.5$ year cycle in in the mean brightness. The long-term average photometric rotation period was $7.24 \pm 0.22$ days. The lower limit for the differential rotation coefficient would be $|k| > 0.12$, if the law of solar surface differential rotation were valid for V352 CMa and the period changes traced this phenomenon. Signs of stable active longitudes rotating with a period of $7.157 \pm 0.002$ days were detected from the epochs of the light minima with the Kuiper method. CPS performed better than the traditional PSM, because the latter always used a sinusoidal model for the data even when this was clearly not the correct model. △ Less

Submitted 16 October, 2013; originally announced October 2013.

Comments: 5 pages, 3 figures, Submitted to A&A

Multiperiodicity, modulations and flip-flops in variable star light curves II. Analysis of II Peg photometry during 1979-2010

Authors: Marjaana Lindborg, Maarit J. Mantere, Nigul Olspert, Jaan Pelt, Thomas Hackman, Gregory W. Henry, Lauri Jetsu, Klaus G. Strassmeier

Abstract: According to earlier Doppler images of the magnetically active primary giant component of the RS CVn binary II Peg, the surface of the star was dominated by one single active longitude that was clearly drifting in the rotational frame of the binary system during 1994-2002; later imaging for 2004-2010, however, showed decreased and chaotic spot activity, with no signs of the drift pattern. Here we… ▽ More According to earlier Doppler images of the magnetically active primary giant component of the RS CVn binary II Peg, the surface of the star was dominated by one single active longitude that was clearly drifting in the rotational frame of the binary system during 1994-2002; later imaging for 2004-2010, however, showed decreased and chaotic spot activity, with no signs of the drift pattern. Here we set out to investigate from a more extensive photometric dataset whether such a drift is a persistent phenomenon, in which case it could be due to either an azimuthal dynamo wave or an indication of the binary system orbital synchronization still being incomplete. We analyse the datasets using the Carrier Fit method (hereafter CF), especially suitable for analyzing time series in which a fast clocking frequency (such as the rotation of the star) is modulated with a slower process (such as the stellar activity cycle). We combine all collected photometric data into one single data set, and analyze it with the CF method. As a result, we confirm the earlier results of the spot activity having been dominated by one primary spotted region almost through the entire data set, and the existence of a persistent, nearly linear drift. Disruptions of the linear trend and complicated phase behavior are also seen, but the period analysis reveals a rather stable periodicity with P(spot)=6.71054d plus/minus 0.00005d. After the linear trend is removed from the data, we identify several abrupt phase jumps, three of which are analyzed in more detail with the CF method. These phase jumps closely resemble what is called flip-flop event, but the new spot configurations do not, in most cases, persist for longer than a few months. △ Less

Submitted 15 April, 2013; originally announced April 2013.

Comments: 9 pages, 7 figures, submitted to Astronomy and Astrophysics

Chemical surface inhomogeneities in late B-type stars with Hg and Mn peculiarity I Spot evolution in HD 11753 on short and long time scales

Authors: H. Korhonen, J. F. Gonzalez, M. Briquet, M. Flores Soriano, S. Hubrig, I. Savanov, T. Hackman, I. V. Ilyin, E. Eulaers, W. Pessemier

Abstract: Aims: Time series of high-resolution spectra of the late B-type star HD 11753 exhibiting HgMn chemical peculiarity are used to study the surface distribution of different chemical elements and their temporal evolution. Methods: High-resolution and high signal-to-noise ratio spectra were obtained using the CORALIE spectrograph at La Silla in 2000, 2009, and 2010. Surface maps of YII, SrII, TiII, an… ▽ More Aims: Time series of high-resolution spectra of the late B-type star HD 11753 exhibiting HgMn chemical peculiarity are used to study the surface distribution of different chemical elements and their temporal evolution. Methods: High-resolution and high signal-to-noise ratio spectra were obtained using the CORALIE spectrograph at La Silla in 2000, 2009, and 2010. Surface maps of YII, SrII, TiII, and CrII were calculated using the Doppler imaging technique. The results were also compared to equivalent width measurements. The evolution of chemical spots both on short and long time scales were investigated. Results: We determine the binary orbit of HD 11753 and fine-tune the rotation period of the primary. The earlier discovered fast evolution of the chemical spots is confirmed by an analysis using both the chemical spot maps and equivalent width measurements. In addition, a long-term decrease in the overall YII and SrII abundances is discovered. A detailed analysis of the chemical spot configurations reveals some possible evidence that a very weak differential rotation is operating in HD 11753. △ Less

Submitted 20 February, 2013; originally announced February 2013.

Comments: 21 pages, 23 figures, accepted for publication in A&A

Magnetic field topology of the RS CVn star II Pegasi

Authors: O. Kochukhov, M. J. Mantere, T. Hackman, I. Ilyin

Abstract: The dynamo processes in cool active stars generate complex magnetic fields responsible for prominent surface stellar activity and variability at different time scales. For a small number of cool stars magnetic field topologies were reconstructed from the time series of spectropolarimetric observations using the Zeeman Doppler imaging (ZDI) method. In this study we follow a long-term evolution of t… ▽ More The dynamo processes in cool active stars generate complex magnetic fields responsible for prominent surface stellar activity and variability at different time scales. For a small number of cool stars magnetic field topologies were reconstructed from the time series of spectropolarimetric observations using the Zeeman Doppler imaging (ZDI) method. In this study we follow a long-term evolution of the magnetic field topology of the RS CVn binary star II Peg. We collected high-resolution circular polarisation observations of II Peg using the SOFIN spectropolarimeter at the Nordic Optical Telescope. These data cover 12 epochs spread over 7 years. A multi-line diagnostic technique in combination with a new ZDI code is applied to interpret these observations. Magnetic inversions using these data reveals evolving magnetic fields with typical local strengths of 0.5-1.0 kG and complex topologies. Despite using a self-consistent magnetic and temperature mapping technique, we do not find a clear correlation between magnetic and temperature features in the ZDI maps. Neither do we confirm the presence of persistent azimuthal field rings found in other RS CVn stars. Reconstruction of the magnetic field topology of II Peg reveals significant evolution of both the surface magnetic field structure and the extended magnetospheric field geometry. From 2004 to 2010 the total field energy drastically declined and the field became less axisymmetric. This also coincided with the transition from predominantly poloidal to mainly toroidal field topology. A qualitative comparison of the ZDI maps of II Peg with the prediction of dynamo theory suggests that the magnetic field in this star is produced mainly by the turbulent alpha^2 dynamo rather than the solar alphaOmega dynamo. Our results do not show a clear active longitude system, nor is there an evidence of the presence of an azimuthal dynamo wave. △ Less

Submitted 8 January, 2013; originally announced January 2013.

Comments: 20 pages, 10 figures; accepted for publication in Astronomy & Astrophysics

Flip-flops of FK Comae Berenices

Authors: Thomas Hackman, Jaan Pelt, Maarit J. Mantere, Lauri Jetsu, Heidi Korhonen, Thomas Granzer, Perttu Kajatkari, Jyri Lehtinen, Klaus G. Strassmeier

Abstract: FK Comae is a rapidly rotating magnetically active star, the light curve of which is modulated by cool spots on its surface. It was the first star where the "flip-flop" phenomenon was discovered. Since then, flip-flops in the spot activity have been reported in many other stars. Therefore, it is of interest to perform a more thorough study of the evolution of the spot activity in FK Com. In this… ▽ More FK Comae is a rapidly rotating magnetically active star, the light curve of which is modulated by cool spots on its surface. It was the first star where the "flip-flop" phenomenon was discovered. Since then, flip-flops in the spot activity have been reported in many other stars. Therefore, it is of interest to perform a more thorough study of the evolution of the spot activity in FK Com. In this study, we analyse 15 years of photometric observations with two different time series analysis methods, with a special emphasis on detecting flip-flop type events from the data. We apply the continuous period search and carrier fit methods on long-term standard Johnson-Cousins V-observations from the years 1995--2010. The observations were carried out with two automated photometric telescopes, Phoenix-10 and Amadeus T7 located in Arizona. We identify complex phase behaviour in 6 of the 15 analysed data segments. We identify five flip-flop events and two cases of phase jumps, where the phase shift is Δφ< 0.4. In addition we see two mergers of spot regions and two cases where the apparent phase shifts are caused by spot regions drifting with respect to each other. Furthermore we detect variations in the rotation period corresponding to a differential rotation coefficient of |k|>0.031. The flip-flop cannot be interpreted as a single phenomenon, where the main activity jumps from one active longitude to another. In some of our cases the phase shifts can be explained by differential rotation: Two spot regions move with different angular velocity and even pass each other. Comparison between the methods show that the carrier fit utility is better in retrieving slow evolution especially from a low amplitude light curve, while the continuous period search is more sensitive in case of rapid changes. △ Less

Submitted 22 March, 2013; v1 submitted 5 November, 2012; originally announced November 2012.

Comments: 14 pages, 11 figures, 4 tables, accepted for publication in Astron. Astrophys

Spot activity of LQ Hya from photometry between 1988 and 2011

Authors: J. Lehtinen, L. Jetsu, T. Hackman, P. Kajatkari, G. W. Henry

Abstract: We investigate the spot activity of the young chromospherically active main sequence star LQ Hya. Our aims are to identify possible active longitudes, estimate the differential rotation and study long and short term changes in the activity. Our analysis is based on 24 years of Johnson V-band photometry. We use the previously published Continuous Period Search (CPS) method to model the evolution of… ▽ More We investigate the spot activity of the young chromospherically active main sequence star LQ Hya. Our aims are to identify possible active longitudes, estimate the differential rotation and study long and short term changes in the activity. Our analysis is based on 24 years of Johnson V-band photometry. We use the previously published Continuous Period Search (CPS) method to model the evolution of the light curve of LQ Hya. The CPS fits a Fourier series model to short overlapping subsets of data. This enables us to monitor the spot configuration of the star with a higher time resolution. We find seasonal variability in the mean level and amplitude of the light curve of LQ Hya. The variability of the light curve amplitude seems not to be cyclic, but the long-term variations in the mean magnitude could be explained by an approximately 13 year cycle. Because of the limited length of the observed time series, it is not yet possible to determine whether this structure really is periodic and represents an activity cycle. We estimate the differential rotation of the star to be small, and the star is potentially very close to a rigid rotator. We search for active longitudes and find that on time scales up to six months there are typically one or two relatively stable active areas on the star with limited phase migration. On time scales longer than one year, no stable active longitudes have been present except for the period between 2003 and 2009 and possibly also some time before 1995. We find any signs of flip-flops with a regular period. The mean time scale of change of the light curve during the observation period is determined to be of the same order of magnitude as the predicted convective turnover time for the star. △ Less

Submitted 22 September, 2015; v1 submitted 7 March, 2012; originally announced March 2012.

Comments: 7 pages, 5 figures, 1 table

Dependence of the large-scale vortex instability on latitude, stratification and domain size

Authors: Maarit J. Mantere, Petri J. Käpylä, Thomas Hackman

Abstract: In an earlier study, we reported on the excitation of large-scale vortices in Cartesian hydrodynamical convection models subject to rapid enough rotation. In that study, the conditions of the onset of the instability were investigated in terms of the Reynolds (Re) and Coriolis (Co) numbers in models located at the stellar North pole. In this study, we extend our investigation to varying domain siz… ▽ More In an earlier study, we reported on the excitation of large-scale vortices in Cartesian hydrodynamical convection models subject to rapid enough rotation. In that study, the conditions of the onset of the instability were investigated in terms of the Reynolds (Re) and Coriolis (Co) numbers in models located at the stellar North pole. In this study, we extend our investigation to varying domain sizes, increasing stratification and place the box at different latitudes. The effect of the increasing box size is to increase the sizes of the generated structures, so that the principal vortex always fills roughly half of the computational domain. The instability becomes stronger in the sense that the temperature anomaly and change in the radial velocity are observed to be enhanced. The model with the smallest box size is found to be stable against the instability, suggesting that a sufficient scale separation between the convective eddies and the scale of the domain is required for the instability to work. The instability can be seen upto the co-latitude of 30 degrees, above which value the flow becomes dominated by other types of mean flows. The instability can also be seen in a model with larger stratification. Unlike the weakly stratified cases, the temperature anomaly caused by the vortex structures is seen to depend on depth. △ Less

Submitted 21 September, 2011; v1 submitted 20 September, 2011; originally announced September 2011.

Comments: 6 pages, 5 figures, submitted to AN (Proc. of the 7th Potsdam Thinkshop)

Report number: NORDITA-2011-79

Doppler images of II Pegasi for 2004-2010

Authors: Thomas Hackman, Maarit J. Mantere, Marjaana Lindborg, Ilya Ilyin, Oleg Kochukhov, Nikolai Piskunov, Ilkka Tuominen

Abstract: Aims. We study the spot activity of II Peg during the years 2004-2010 to determine long- and short-term changes in the magnetic activity. In a previous study, we detected a persistent active longitude, as well as major changes in the spot configuration occurring on a timescale of shorter than a year. The main objective of this study is to determine whether the same phenomena persist in the star du… ▽ More Aims. We study the spot activity of II Peg during the years 2004-2010 to determine long- and short-term changes in the magnetic activity. In a previous study, we detected a persistent active longitude, as well as major changes in the spot configuration occurring on a timescale of shorter than a year. The main objective of this study is to determine whether the same phenomena persist in the star during these six years of spectroscopic monitoring. Methods. The observations were collected with the high-resolution SOFIN spectrograph at the Nordic Optical Telescope. The temperature maps were calculated using a Doppler imaging code based on Tikhonov regularization. Results. We present 12 new temperature maps that show spots distributed mainly over high and intermediate latitudes. In each image, 1-3 main active regions can be identified. The activity level of the star is clearly lower than during our previous study for the years 1994-2002. In contrast to the previous observations, we detect no clear drift of the active regions with respect to the rotation of the star. Conclusions. Having shown a systematic longitudinal drift of the spot-generating mechanism during 1994-2002, the star has clearly switched to a low-activity state for 2004-2010, during which the spot locations appear more random over phase space. It could be that the star is near to a minimum of its activity cycle. △ Less

Submitted 18 January, 2012; v1 submitted 30 June, 2011; originally announced June 2011.

Comments: Accepted for publication in Astron. and Astrophys., 8 pages, 5 figures

Starspots due to large-scale vortices in rotating turbulent convection

Authors: P. J. Käpylä, M. J. Mantere, T. Hackman

Abstract: We study the generation of large-scale vortices in rotating turbulent convection by means of Cartesian direct numerical simulations. We find that for sufficiently rapid rotation, cyclonic structures on a scale large in comparison to that of the convective eddies, emerge, provided that the fluid Reynolds number exceeds a critical value. For slower rotation, cool cyclonic vortices are preferred, whe… ▽ More We study the generation of large-scale vortices in rotating turbulent convection by means of Cartesian direct numerical simulations. We find that for sufficiently rapid rotation, cyclonic structures on a scale large in comparison to that of the convective eddies, emerge, provided that the fluid Reynolds number exceeds a critical value. For slower rotation, cool cyclonic vortices are preferred, whereas for rapid rotation, warm anti-cyclonic vortices are favored. In some runs in the intermediate regime both types of cyclones co-exist for thousands of convective turnover times. The temperature contrast between the vortices and the surrounding atmosphere is of the order of five per cent. We relate the simulation results to observations of rapidly rotating late-type stars that are known to exhibit large high-latitude spots from Doppler imaging. In many cases, cool spots are accompanied with spotted regions with temperatures higher than the average. In this paper, we investigate a scenario according to which the spots observed in the temperature maps could have a non-magnetic origin due to large-scale vortices in the convection zones of the stars. △ Less

Submitted 18 November, 2013; v1 submitted 29 June, 2011; originally announced June 2011.

Comments: 8 pages, 8 figures, version published in ApJ. Version with higher resolution figures can be found at http://www.helsinki.fi/~kapyla/publ.html

Report number: NORDITA-2011-53

Journal ref: Astrophys. J., 742, 34 (2011)

The continuous period search method and its application to the young solar analogue HD 116956

Authors: J. Lehtinen, L. Jetsu, T. Hackman, P. Kajatkari, G. W. Henry

Abstract: We formulate an improved time series analysis method for the analysis of photometry of active stars. This new Continuous Period Search (CPS) method is applied to 12 years of V band photometry of the young solar analogue HD 116956 (NQ UMa). The new method is developed from the previous Three Stage Period Analysis (TSPA) method. Our improvements are the use of a sliding window in choosing the modell… ▽ More We formulate an improved time series analysis method for the analysis of photometry of active stars. This new Continuous Period Search (CPS) method is applied to 12 years of V band photometry of the young solar analogue HD 116956 (NQ UMa). The new method is developed from the previous Three Stage Period Analysis (TSPA) method. Our improvements are the use of a sliding window in choosing the modelled datasets, a criterion applied to select the best model for each dataset and the computation of the time scale of change of the light curve. We test the performance of CPS with simulated and real data. The CPS has a much improved time resolution which allows us to better investigate fast evolution of stellar light curves. We can also separate between the cases when the data is best described by periodic and aperiodic (e.g. constant brightness) models. We find, however, that the performance of the CPS has certain limitations. It does not determine the correct model complexity in all cases, especially when the underlying light curve is constant and the number of observations too small. Also the sensitivity in detecting two close light curve minima is limited and it has a certain amount of intrinsic instability in its period estimation. Using the CPS, we find persistent active longitudes in the star HD 116956 and a "flip-flop" event that occurred during the year 1999. Assuming that the surface differential rotation of the star causes observable period variations in the stellar light curve, we determine the differential rotation coefficient to be |k|>0.11. The mean timescale of change of the light curve during the whole 12 year observing period was T_C=44.1 d, which is of the same order as the predicted convective turnover time of the star. We also investigate the presence of activity cycles on the star, but do not find any conclusive evidence supporting them. △ Less

Submitted 3 March, 2011; v1 submitted 23 July, 2010; originally announced July 2010.

Comments: 14 pages, 11 figures, 3 tables

Dynamical evolution of titanium, strontium, and yttrium spots on the surface of the HgMn star HD 11753

Authors: M. Briquet, H. Korhonen, J. F. Gonzalez, S. Hubrig, T. Hackman

Abstract: Aims. We gathered about 100 high-resolution spectra of three typical HgMn (mercury-manganese) stars, HD 11753, HD 53244, and HD 221507, to search for slowly pulsating B-like pulsations and surface inhomogeneous distribution of various chemical elements. Methods. Classical frequency analysis methods were used to detect line profile variability and to determine the variation period. Doppler imaging… ▽ More Aims. We gathered about 100 high-resolution spectra of three typical HgMn (mercury-manganese) stars, HD 11753, HD 53244, and HD 221507, to search for slowly pulsating B-like pulsations and surface inhomogeneous distribution of various chemical elements. Methods. Classical frequency analysis methods were used to detect line profile variability and to determine the variation period. Doppler imaging reconstruction was performed to obtain abundance maps of chemical elements on the stellar surface. Results. For HD 11753, which is the star with the most pronounced variability, distinct spectral line profile changes were detected for Ti, Sr, Y, Zr, and Hg, whereas for HD 53244 and HD 221507 the most variable line profiles belong to the elements Hg and Y, respectively. We derived rotation periods for all three stars from the variations of radial velocities and equivalent widths of spectral lines belonging to inhomogeneously distributed elements: P_rot (HD 11753)=9.54 d, P_rot (HD 53244)=6.16 d, and P_rot (HD 221507)=1.93 d. For HD 11753 the Doppler imaging technique was applied to derive the distribution of the most variable elements Ti, Sr, and Y using two datasets separated by ~65 days. Results of Doppler imaging reconstruction revealed noticeable changes in the surface distributions of TiII, SrII, and YII between the datasets, indicating the hitherto not well understood physical processes in stars with radiative envelopes that cause a rather fast dynamical chemical spot evolution. △ Less

Submitted 9 March, 2010; originally announced March 2010.

Comments: Accepted for publication in Astronomy and Astrophysics, 6 pages, 9 figures, 2 tables.

Ellipsoidal primary of the RS CVn binary zeta And: Investigation using high-resolution spectroscopy and optical interferometry

Authors: H. Korhonen, M. Wittkowski, Zs. Kovari, Th. Granzer, T. Hackman, K. G. Strassmeier

Abstract: We have obtained high-resolution spectroscopy, optical interferometry, and long-term broad band photometry of the ellipsoidal primary of the RS CVn-type binary system zeta And. Based on the optical interferometry the apparent limb darkened diameter of zeta And is 2.55 +/- 0.09 mas using a uniform disk fit. The Hipparcos distance and the limb-darkened diameter obtained with a uniform disk fit giv… ▽ More We have obtained high-resolution spectroscopy, optical interferometry, and long-term broad band photometry of the ellipsoidal primary of the RS CVn-type binary system zeta And. Based on the optical interferometry the apparent limb darkened diameter of zeta And is 2.55 +/- 0.09 mas using a uniform disk fit. The Hipparcos distance and the limb-darkened diameter obtained with a uniform disk fit give stellar radius of 15.9 +/- 0.8 Rsolar, and combined with bolometric luminosity, it implies an effective temperature of 4665 +/- 140 K. The temperature maps obtained from high resolution spectra using Doppler imaging show a strong belt of equatorial spots and hints of a cool polar cap. The equatorial spots show a concentration around the phase 0.75. This spot configuration is reminiscent of the one seen in the earlier published temperature maps of zeta And. Investigation of the Halpha line reveals both prominences and cool clouds in the chromosphere. Long-term photometry spanning 12 years shows hints of a spot activity cycle, which is also implied by the Doppler images, but the cycle length cannot be reliably determined from the current data. △ Less

Submitted 22 February, 2010; originally announced February 2010.

Comments: 9 pages, 9 figures, accepted for A&A

Surface temperature maps for II Peg during 1999-2002

Authors: M. Lindborg, M. J. Korpi, I. Tuominen, T. Hackman, I. Ilyin, N. Piskunov

Abstract: The active RS CVn star II Peg has been spectroscopically monitored for almost 18 years with the SOFIN spectrograph at NOT, La Palma, Spain. In this paper we present five new surface temperature maps of the object for the years 1999 (two maps), 2001 (one map) and 2002 (two maps). The active RS CVn star II Peg has been spectroscopically monitored for almost 18 years with the SOFIN spectrograph at NOT, La Palma, Spain. In this paper we present five new surface temperature maps of the object for the years 1999 (two maps), 2001 (one map) and 2002 (two maps). △ Less

Submitted 1 November, 2009; originally announced November 2009.

First measurement of the magnetic field on FK Com and its relation to the contemporaneous starspot locations

Authors: H. Korhonen, S. Hubrig, S. V. Berdyugina, Th. Granzer, T. Hackman, M. Schoeller, K. G. Strassmeier, M. Weber

Abstract: In this study we present simultaneous low-resolution longitudinal magnetic field measurements and high-resolution spectroscopic observations of the cool single giant FK Com. The variation of the magnetic field over the rotational period of 2.4 days is compared with the starspot location obtained using Doppler imaging techniques, V-band photometry and V-I colours. The chromospheric activity is st… ▽ More In this study we present simultaneous low-resolution longitudinal magnetic field measurements and high-resolution spectroscopic observations of the cool single giant FK Com. The variation of the magnetic field over the rotational period of 2.4 days is compared with the starspot location obtained using Doppler imaging techniques, V-band photometry and V-I colours. The chromospheric activity is studied simultaneously with the photospheric activity using high resolution observations of the Halpha, Hbeta and Hgamma line profiles. Both the maximum (272 +/- 24 G) and minimum (60 +/- 17 G) in the mean longitudinal magnetic field, <Bz>, are detected close to the phases where cool spots appear on the stellar surface. A possible explanation for such a behaviour is that the active regions at the two longitudes separated by 0.2 in phase have opposite polarities. △ Less

Submitted 2 December, 2008; originally announced December 2008.

Comments: 10 Pages, 11 figures (quality of Figures 7,8 and 10 reduced), accepted for publication in MNRAS

Surface imaging of HD 199178 (V1794 Cygni)

Authors: T. Hackman, L. Jetsu, I. Tuominen

Abstract: We present surface temperature maps for the FK Comae-type star HD 199178 (V1794 Cygni) calculated from high resolution spectra obtained in 1994 and 1995. The spot pattern evolves, but all maps reveal a large cool spot remaining nearly at the same high latitude. The main spot is 1200 -- 1600 K cooler than the mean surface temperature. The observed slightly flat bottomed absorption lines would usu… ▽ More We present surface temperature maps for the FK Comae-type star HD 199178 (V1794 Cygni) calculated from high resolution spectra obtained in 1994 and 1995. The spot pattern evolves, but all maps reveal a large cool spot remaining nearly at the same high latitude. The main spot is 1200 -- 1600 K cooler than the mean surface temperature. The observed slightly flat bottomed absorption lines would usually be interpreted as evidence for a large cool polar spot. We argue that antisolar surface differential rotation offers a better explanation for the box like shape of the line profiles. However, we do not find conclusive evidence for antisolar differential rotation and note that there are still other possible explanations for the slightly flat bottomed line profiles. △ Less

Submitted 6 June, 2001; originally announced June 2001.

Comments: 10 pages, 5 figures, Astron. Astrophys., in press

Study of FK Comae Berenices. II. Spot evolution from 1994 to 1997

Authors: H. Korhonen, S. V. Berdyugina, T. Hackman, K. G. Strassmeier, I. Tuominen

Abstract: We present new surface (Doppler) images of the late-type single giant FK Com for June-July 1996, July-August 1996, April 1997 and June 1997. These images are compared with the previously published images from 1994 and 1995. The consecutive maps are cross-correlated to see the possible migration of the spots and the effects of differential rotation. The cross-correlation confirms an average longi… ▽ More We present new surface (Doppler) images of the late-type single giant FK Com for June-July 1996, July-August 1996, April 1997 and June 1997. These images are compared with the previously published images from 1994 and 1995. The consecutive maps are cross-correlated to see the possible migration of the spots and the effects of differential rotation. The cross-correlation confirms an average longitudinal spot migration of 0.22+/-0.03 in phase within a year. This movement is probably an artifact caused by a difference between the accepted rotation period and the real photometric period for these years. If this is true, then the photometric rotation period for these years is 2.4037+/-0.0005 days. Measurements from these four years and six maps limit the surface differential rotation to alpha=0.0001+/-0.0002, where alpha is the difference between polar and equatorial angular velocities relative to the equatorial angular velocity. △ Less

Submitted 7 July, 2000; originally announced July 2000.

Comments: 10 pages, 7 figures, accepted for publication in A&A