The first magnetic maps of a pre-main-sequence binary star system – HD 155555

We present the first maps of the surface magnetic fields of a pre-main-sequence binary system. Spectropolarimetric observations of the young, 18 Myr, HD 155555 (V824 Ara, G5IV+K0IV) system were obtained at the Anglo-Australian Telescope in 2004 and 2007. Both data sets are analysed using a new binary Zeeman–Doppler imaging (ZDI) code. This allows us to simultaneously model the contribution of each component to the observed circularly polarized spectra. Stellar brightness maps are also produced for HD 155555 and compared to previous Doppler images.
Our radial magnetic maps reveal a complex surface magnetic topology with mixed polarities at all latitudes. We find rings of azimuthal field on both stars, most of which are found to be non-axisymmetric with the stellar rotational axis. We also examine the field strength and the relative fraction of magnetic energy stored in the radial and azimuthal field components at both epochs. A marked weakening of the field strength of the secondary star is observed between the 2004 and 2007 epochs. This is accompanied by an apparent shift in the location of magnetic energy from the azimuthal to radial field. We suggest that this could be indicative of a magnetic activity cycle. We use the radial magnetic maps to extrapolate the coronal field (by assuming a potential field) for each star individually – at present ignoring any possible interaction. The secondary star is found to exhibit an extreme tilt (≈75°) of its large-scale magnetic field to that of its rotation axis for both epochs. The field complexity that is apparent in the surface maps persists out to a significant fraction of the binary separation. Any interaction between the fields of the two stars is therefore likely to be complex also. Modelling this would require a full binary field extrapolation.     

Stellar surface imaging: mapping brightness and magnetic fields

Rapidly rotating late‐type stars typically display signs of magnetic activity that exceed those seen on the Sun by over two orders of magnitude. The techniques of Doppler imaging and Zeeman Doppler imaging have been instrumental in unveiling magnetic activity patterns at the photospheres of these active stars. Essentially, these techniques work by inverting time‐series of high resolution spectra to produce temperature, brightness and/or magnetic field maps at the surfaces of stars. I will describe how these techniques work and review what they have taught us about the nature of magnetic activity in rapid rotators over the last 20 years. Finally, I will conclude by outlining the capabilities of these techniques in light of new instrumentation that is now becoming available. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spots on FK Com: active longitudes and “flip‐flops”

We have earlier investigated the surface structures of a late‐type, single, giant FK Com for the years 1994–1998 using Doppler imaging. These surface temperature maps revealed long‐lived active regions at high latitudes. Long‐term photometric observations also show that these active regions tend to occur at two permanent active longitudes which are 180 degrees apart from each other, and that the activity switches the longitude with an average period of about 3 years (the “flip‐flop” phenomenon). In this work we present new Doppler maps of FK Com obtained 1998‐2003 and light‐curve maps obtained 2002–2003. These new maps are investigated together with the earlier temperature maps and light‐curve maps, with an aim of further studying the active longitudes, “flip‐flop” phenomenon and surface differential rotation on FK Com. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The symbiotic star CH Cygni – I. Non-thermal bipolar jets

Very Large Array surface brightness and spectral index maps of the evolving extended emission of the triple symbiotic star CH Cygni are presented. These are derived from observations at 4.8, 8.4 and 14 GHz between 1985 and 1999. The maps are dominated by thermal emission around the central bright peak of the nebula, but we also find unambiguous non-thermal emission associated with the extended regions. Our observations confirm that this is a jet. The central region has been associated with the stellar components through Hubble Space Telescope imaging. If the jets are the result of ejection events at outburst, expansion velocities are consistent with those from other measurement methods. We propose that the non-thermal emission is caused by material ejected in the bipolar jets interacting with the circumstellar wind envelope. The resulting shocks lead to local enhancements in the magnetic field from the compact component of the order of 3 mG.     

Surface magnetic fields on two accreting T Tauri stars: CV Cha and CR Cha

We have produced brightness and magnetic field maps of the surfaces of CV Cha and CR Cha: two actively accreting G- and K-type T Tauri stars in the Chamaeleon I star-forming cloud with ages of 3–5 Myr. Our magnetic field maps show evidence for strong, complex multipolar fields similar to those obtained for young rapidly rotating main-sequence stars. Brightness maps indicate the presence of dark polar caps and low-latitude spots – these brightness maps are very similar to those obtained for other pre-main-sequence and rapidly rotating main-sequence stars.
Only two other classical T Tauri stars have been studied using similar techniques so far: V2129 Oph and BP Tau. CV Cha and CR Cha show magnetic field patterns that are significantly more complex than those recovered for BP Tau, a fully convective T Tauri star.
We discuss possible reasons for this difference and suggest that the complexity of the stellar magnetic field is related to the convection zone; with more complex fields being found in T Tauri stars with radiative cores (V2129 Oph, CV Cha and CR Cha). However, it is clearly necessary to conduct magnetic field studies of T Tauri star systems, exploring a wide range of stellar parameters in order to establish how they affect magnetic field generation, and thus how these magnetic fields are likely to affect the evolution of T Tauri star systems as they approach the main sequence.     

Doppler images and chromospheric variability of TWA 17

We present Doppler imaging and a Balmer line analysis of the weak-line T Tauri star TWA 17. Spectra were taken in 2006 with the University College London Echelle Spectrograph on the Anglo-Australian Telescope. Using least-squares deconvolution to improve the effective signal-to-noise ratio, we produced a Doppler map of the surface spot distribution. This shows similar features to maps of other rapidly rotating T Tauri stars, i.e. a polar spot with more spots extending out of it down to the equator.
In addition to the photospheric variability, the chromospheric variability was studied using the Balmer emission. The mean Hα profile has a narrow component consistent with rotational broadening and a broad component extending out to 220 km s⁻¹. The variability in Hα suggests that the chromosphere has at least one slingshot prominence  3 R _*  above the surface.     

The global magnetic topology of AB Doradus

We have used Zeeman–Doppler maps of the surface field of the young, rapid rotator AB Dor  ( P _rot=0.514 d)  to extrapolate the coronal field, assuming it to be potential. We find that the topology of the large-scale field is very similar in all three years for which we have images. The corona divides cleanly into regions of open and closed field. The open field originates in two mid-latitude regions of opposite polarity separated by about 180° of longitude. The closed field region forms a torus extending almost over each pole, with an axis that runs through these two longitudes. We have investigated the effect on the global topology of different forms of flux in the unobservable hemisphere and in the dark polar spot where the Zeeman signal is suppressed. The flux distribution in the unobservable hemisphere affects only the low-latitude topology, whereas the imposition of a unidirectional polar field forces the polar cap to be open. This contradicts observations that suggest that the closed field corona extends to high latitudes and leads us to propose that the polar cap may be composed of multipolar regions.     

Indirect imaging of the accretion stream in eclipsing polars – III. HU Aquarii low state

We apply our technique for indirect imaging of the accretion stream to the polar HU Aqr, using eclipse profiles observed when the system was in a low-accretion state. The eclipse profile is different from that in the high state, and more variable from cycle to cycle. We find that the stream maps are brightest near the white dwarf and there is no significant brightening in the threading region. In the low state the stream threads on to the magnetic field closer to the L₁ point than in the high state, with a footpoint of the accreting field line at high latitude. We then produce maps of the accretion region from polarimetry using Stokes imaging. These show that the majority of the accretion occurs near the equator. The difference between the maps may be explained if most of the stream material is not emitting significantly in the low state. If so, neither the stream eclipse mapping nor Doppler tomography techniques will trace the bulk of the accretion flow between the two stars.     

The chromospherically active binary star EI Eridani: II. Long‐term Doppler imaging

Data from 11 years of continuous spectroscopic observations of the active RS CVn‐type binary star EI Eridani – gained at NSO/McMath‐Pierce, KPNO/Coudé Feed and during the MUSICOS 98 campaign – were used to obtain 34 Doppler maps in three spectroscopic lines for 32 epochs, 28 of which are independent of each other. Various parameters are extracted from our Doppler maps: average temperature, fractional spottedness, and longitudinal and latitudinal spot‐occurrence functions. We find that none of these parameters show a distinct variation nor a correlation with the proposed activity cycle as seen from photometric long‐term observations. This suggests that the photometric brightness cycle may not necessarily be due to just a cool spot cycle. The general morphology of the spot pattern remains persistent over the whole period of 11 years. A large cap‐like polar spot was recovered from all our images. A high degree of variable activity was noticed near latitudes of ≈60–70° where the appendages of the polar spot emerged and dissolved (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hubble Space Telescope observations of SV Cam – I. The importance of unresolved star-spot distributions in light-curve fitting

We have used maximum entropy eclipse-mapping to recover images of the visual surface brightness distribution of the primary component of the RS CVn eclipsing binary SV Cam, using high-precision photometry data obtained during three primary eclipses with Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope ( HST ). These were augmented by contemporaneous ground-based photometry secured around the rest of the orbit. The goal of these observations was to determine the filling factor and size distribution of star-spots too small to be resolved by Doppler imaging. The information content of the final image and the fit to the data were optimized with respect to various system parameters using the χ² landscape method, using an eclipse-mapping code that solves for large-scale spot coverage. It is only with the unprecedented photometric precision of the HST data (0.000 15 mag) that it is possible to see strong discontinuities at the four contact points in the residuals of the fit to the light curve. These features can only be removed from the residual light curve by the reduction of the photospheric temperature, to synthesize high unresolvable spot coverage, and the inclusion of a polar spot. We show that this spottedness of the stellar surface can have a significant impact on the determination of the stellar binary parameters and the fit to the light curve by reducing the secondary radius from  0.794 ± 0.009  to  0.727 ± 0.009 R_⊙  . This new technique can also be applied to other binary systems with high-precision spectrophotometric observations.     

Hubble Space Telescope observations of SV Cam – II. First derivative light-curve modelling using phoenix and atlas model atmospheres

The variation of the specific intensity across the stellar disc is an essential input parameter in surface brightness reconstruction techniques such as Doppler imaging, where the relative intensity contributions of different surface elements are important in detecting star-spots. We use phoenix and atlas model atmospheres to model light curves derived from high precision (signal-to-noise ratio ≃ 5000) Hubble Space Telescope ( HST ) data of the eclipsing binary SV Cam (F9V+K4V), where the variation of specific intensity across the stellar disc will determine the contact points of the binary system light curve. For the first time, we use χ² comparison fits to the first derivative profiles to determine the best-fitting model atmosphere. We show the wavelength dependence of the limb darkening and that the first derivative profile is sensitive to the limb-darkening profile very close to the limb of the primary star. It is concluded that there is only a marginal difference (<1σ) between the χ² comparison fits of the two model atmospheres to the HST light curve at all wavelengths. The usefulness of the second derivative of the light curve for measuring the sharpness of the primary's limb is investigated, but we find that the data are too noisy to permit a quantitative analysis.     

Light‐curve inversions with truncated least‐squares principal components: Tests and application to HD 291095 = V1355 Orionis

We present a new inversion code that reconstructs the stellar surface spot configuration from the light curve of a rotating star. Our code employs a method that uses the truncated least‐squares estimation of the inverse problem's objects principal components. We use spot filling factors as the unknown objects. Various test cases that represent a rapidly‐rotating K subgiant are used for the forward problem. Tests are then performed to recover the artificial input map and include data errors and input‐parameter errors. We demonstrate the robustness of the solution to false input parameters like photospheric temperature, spot temperature, gravity, inclination, unspotted brightness and different spot distributions and we also demonstrate the insensitivity of the solution to spot latitude. Tests with spots peppered over the entire stellar surface or with phase gaps do not produce fake active longitudes. The code is then applied to ten years of V and I ‐band light curve data of the spotted sub‐giant HD291095. A total of 22 light curves is presented. We find that for most of the time its spots were grouped around two active longitudes separated on average by 180°. Switches of the dominant active region between these two longitudes likely occurred about every 3.15±0.23 years while the amplitude modulation of the brightness occurred with a possible period of 3.0±0.15 years. For the first time, we found evidence that the times of the activity flips coincide with times of minimum light as well as minimum photometric amplitude, i.e. maximum spottedness. From a comparison with simultaneous Doppler images we conclude that the activity flips likely take place near the rotational pole of the star. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Near-infrared and optical studies of the fast nova V4643 Sgr (Nova Sagittarii 2001)

It is generally accepted that the presence of a hot magnetic corona provides the source of X-ray emission in cool stars. With this connection one could expect to see the variation of magnetic flux in the activity cycle of a star mirrored by a similar variation in the stars X-ray emission. Using magnetic maps produced from flux emergence and transport simulations and assuming a potential field for the corona, we can extrapolate the coronal magnetic field and hence calculate the variation of the X-ray emission. We consider three types of activity cycle that successfully reproduce the pattern of intermingled magnetic flux at high latitudes, a feature observed with Zeeman–Doppler imaging. The three different cycles take the form of (1) an enhanced butterfly pattern where flux emergence is extended to a latitude of 70°, (2) an extended emergence profile as before but with an overlap of 4 yr in the butterfly diagram and (3) where no butterfly diagram is used. The cyclic variation in the X-ray emission is around two orders of magnitude for cases (1) and (3), but less than one order of magnitude for case (2). For all three cases, the rotational modulation of the X-ray emission is greatest at cycle minimum, but the emission measure weighted density varies little over the cycle. For cases (1) and (2) the fraction of the total flux that is open (along which a wind can escape) varies little over the cycle, but for case (3) this is three times larger at cycle minimum than at maximum. Our results clearly show that although magnetic cycles may exist for stars they are not necessarily observable in the X-ray emission.     

Long‐term photometry of the ultra‐fast rotator LO Peg

LO Peg is a young main‐sequence star of spectral type K3. With its equatorial rotation velocity of 65 km s^–1 it is amongst the ultra‐fast rotators. Its high equatorial rotation velocity and rapidly changing surface activity features make it an important object in terms of both stellar activity and the evolution of stellar rotation and angular momentum. Since its discovery as a variable star, it has mostly been subject to spectral surface mapping studies such as Doppler Imaging, while there have been very few photometric studies on it. This paper aims to present the first long‐term photometric observations and its results covering the years between 2003 and 2009. The UBVR Johnson wide band photometric data showed that the surface activity structures of LO Peg vary in timescales changing between days and months, and parallel to this, the mean, maximum and minimum brightness and amplitudes change dramatically between years and sometimes even within the same observation season. Long‐term changes in system brightness and colours, both characteristic features of active stars, were also seen in this ultra‐fast young star. The active longitudes, which has a life time of ∼1.3 years and an activity cycle period of ∼4.8 years for LO Peg were estimated (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling surface magnetic field evolution on AB Doradûs due to diffusion and surface differential rotation

From Zeeman–Doppler images of the young, rapidly-rotating K0 dwarf AB Doradûs, we have created a potential approximation to the observed radial magnetic field and have evolved it over 30 d subject to the observed surface differential rotation , meridional flow and various diffusion rates. Assuming that the dark polar cap seen in Doppler images of this star is caused by the presence of a unipolar field, we have shown that the observed differential rotation will shear this field to produce the observed high-latitude band of unidirectional azimuthal field. By cross-correlating the evolved fields with the initial field each day we have followed the decay with time of the cross-correlation function. Over 30 d it decays by only 10 per cent. This contrasts with the results of Barnes et al. , who show that on this time-scale the spot distribution of He699 is uncorrelated. We propose that this is due to the effects of flux emergence changing the spot distributions.     

On the time dependence of differential rotation in young late-type stars

A model for the angular momentum transfer within the convection zone of a rapidly rotating star is introduced and applied to the analysis of recent observations of temporal fluctuations of the differential rotation on the young late-type stars AB Doradus (AB Dor) and LQ Hydrae (LQ Hya). Under the hypothesis that the mean magnetic field produced by the stellar dynamo rules the angular momentum exchanges and that the angular velocity depends only on the distance s from the rotation axis and the time, the minimum azimuthal Maxwell stress  | B_sB _φ|  , averaged over the convection zone, is found to range from ∼0.04 to  ∼0.14 T²  . If the poloidal mean magnetic field   B _s   is of the order of 0.01 T, as indicated by the Zeeman–Doppler imaging maps of those stars, then the azimuthal mean field   B _φ  can reach an intensity of several teslas, which significantly exceeds equipartition with the turbulent kinetic energy. Such strong fields can account also for the orbital period modulation observed in cataclysmic variables and RS Canum Venaticorum systems with a main-sequence secondary component. Moreover, the model allows us to compute the kinetic energy dissipation rate during the maintenance of the differential rotation. Only in the case of the largest surface shear observed on LQ Hya may the dissipated power exceed the stellar luminosity, but the lack of a sufficient statistic on the occurrence of such episodes of large shear does not allow us to estimate their impact on the energy budget of the convection zone.     

类太阳活动恒星表面黑子结构的多普勒成像(Ⅰ)：技术

回顾了恒星表面结构多普勒成像研究工作的发展历史,重点评述了该项技术在类太阳活动恒星研究领域的应用。介绍了该项技术在国际上的发展现状,给出了恒星表面黑子结构形态的多普勒成像技术的基本原理,讨论了该项技术对观测仪器、观测过程和观测对象的要求,并对未来这一领域的研究工作进行了展望。     

Deciphering the pulsations of G 29-38 with optical time series spectroscopy

We present optical time series spectroscopy of the pulsating white dwarf star G 29-38 taken at the Very Large Telescope (VLT). By measuring the variations in brightness, Doppler shift and line shape of each spectrum, we explore the physics of pulsation and measure the spherical degree (ℓ) of each stellar pulsation mode. We measure the physical motion of the g modes correlated with the brightness variations for three of the eight pulsation modes in this data set. The varying line shape reveals the spherical degree of the pulsations, an important quantity for properly modelling the interior of the star with asteroseismology. Performing fits to the Hβ, Hγ and Hδ lines, we quantify the changing shape of the line and compare them to models and previous time series spectroscopy of G 29-38. These VLT data confirm several ℓ identifications and add four new values, including an additional ℓ= 2 and a possible ℓ= 4. In total, from both sets of spectroscopy of G 29-38, eleven modes now have known spherical degrees.     

Differential rotation of giant stars

From a set of high-resolution spectral observations of late type giant stars we used Doppler imaging to derive time-series temperature maps of the stellar surfaces. Using these temperature maps, it is possible to track the temporal changes of the spot features and derive estimates of the strength and sign of the differential surface rotation of these stars. Looking into the latitudinal changes of the surface maps, it is also possible to derive meridional flows on these stars. But due to the lower accuracy of the latitudes of the reconstructed spot features, the data requirements are higher than for the detection of differential rotation. Nevertheless, a correlation between the differential rotation and meridional flow estimates is suggested. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)