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.     

Starspots and relativity: Applied Doppler imaging for the Gravity Probe B mission

We present Doppler images and surface differential rotation measurements for the primary of the RS CVn binary IM Pegasi, the guide star for the Gravity Probe B experiment. The data used is a subset of that taken during optical support of the mission and was obtained almost nightly over a near three year period from the Automatic Spectroscopic Telescope operated by Tennessee State University. Using the technique of least-squares deconvolution to increase the signal-to-noise ratio of the data, we have reconstructed 31 maximum entropy Doppler images of the star. The images show that the spot features are relatively stable for over a year (and possibly longer) with both a polar spot and lower latitude features. The most intense features are located on the side facing the secondary. In addition, we have incorporated a solar-like differential rotation law into the imaging process to determine the level of surface differential rotation for IM Peg for 22 epochs. A weighted least-squares average of the measurements gives a surface shear of 0.0142 ± 0.0007 rad/d, meaning that the equator takes ∼440 ± 20 days to lap the poles. Although the level of surface differential rotation was shown to vary over the period of the observations, this may indicate an underestimate in the errors of the method rather than any temporal evolution in the differential rotation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spot patterns and differential rotation in the eclipsing pre-cataclysmic variable binary, V471 Tau

We present surface spot maps of the K2V primary star in the pre-cataclysmic variable binary system, V471 Tau. The spot maps show the presence of large high-latitude spots located at the sub-white dwarf longitude region. By tracking the relative movement of spot groups over the course of four nights (eight rotation cycles), we measure the surface differential rotation rate of the system. Our results reveal that the star is rotating rigidly with a surface shear rate,  dΩ= 1.6 ± 6 mrad d⁻¹  . The single active star AB Dor has a similar spectral type, rotation period and activity level as the K star in V471 Tau, but displays much stronger surface shear  (46 < dΩ < 58 mrad d⁻¹)  . Our results suggest that tidal locking may inhibit differential rotation; this reduced shear, however, does not affect the overall magnetic activity levels in active K dwarfs.     

The rapidly oscillating Ap star HD 99563 and its distorted dipole pulsation mode

We undertook a time-series photometric multisite campaign for the rapidly oscillating Ap (roAp) star HD 99563 and also acquired mean light observations over four seasons. The pulsations of the star, which show flatter light maxima than minima, can be described with a frequency quintuplet centred on 1557.653 μHz and some first harmonics of it. The amplitude of the pulsation is modulated with the rotation period of the star that we determine with 2.91179 ± 0.00007 d from the analysis of the stellar pulsation spectrum and of the mean light data. We break up the distorted oscillation mode into its pure spherical harmonic components and find it is dominated by the ℓ= 1 pulsation, and also has a notable ℓ= 3 contribution, with weak ℓ= 0 and 2 components. The geometrical configuration of the star allows us to see both pulsation poles for about the same amount of time; HD 99563 is only the fourth roAp star for which both pulsation poles are seen and only the third where the distortion of the pulsation modes has been modelled. We point out that HD 99563 is very similar to the well-studied roAp star HR 3831. Finally, we note that the visual companion of HD 99563 is located in the δ Scuti instability strip and may thus show pulsation. We show that if the companion was physical, the roAp star would be a 2.03-M_⊙, object, seen at a rotational inclination of 44°, which then predicts a magnetic obliquity     .     

Optical, X-ray and radio observations of HD 61396: a probable new RS CVn-type binary

We have carried out BVR photometric and H spectroscopic observations of the star HD 61396 during 1998 March 20 to 1999 April 3. We have discovered regular optical photometric variability from this star, with an inferred period of 31.95±0.10 d, and an amplitude of 0.18 mag. A possible period of 35.34±0.12 d, as determined with Hipparcos , cannot be completely ruled out, however. Modelling of its photometric light curve with two circular spots indicates that 521 per cent of the stellar surface is covered by dark starspots which are 830 K cooler than the surrounding photosphere, and produce the observed rotational modulation of the optical flux. Optical spectroscopy reveals a variable H emission feature, indicating that it is an unusually active star.
In addition, we have analysed archival X-ray data of HD 61396, obtained from serendipitous observations with the ROSAT X-ray observatory, and we also discuss the radio properties of this star, based on both published Green Bank and unpublished VLA observations. The strong photometric variability and H emission, the relatively hard X-ray spectrum, and the high X-ray and radio luminosities imply that HD 61396 is most likely to be a member of the RS CVn class of evolved active binary stars. Its X-ray and radio luminosities place it among the five most luminous active binaries detected so far.     

Eight new δ Scuti stars

HD 23194, a member of the Pleiades, was found to pulsate with a period of about 30 min. The literature on the star is reviewed, and it is concluded that it may be a marginal Am star in a binary system. HD 95321 is an evolved Am ( ρ Puppis) star with a 5.1-h periodicity. Mode identification of its pulsation, based on multicolour photometry, suggests that the oscillation is probably non-radial with ℓ=2. We also report on the discovery of six other new δ Scuti stars, some of which may be pulsating in gravity modes.     

Stability of magnetic flux tubes in close binary stars

Photometric and Doppler imaging observations of active binaries indicate the existence of starspots at preferred longitudes (position angles with respect to the companion star). We investigate the stability of magnetic flux tubes in the convection zone of close, fast‐rotating binary stars and explore whether the observed preferred longitudes could be caused by tidal forces and the deformation of the active star. We assume a synchronized binary system with spin axes perpendicular to the orbital plane and a rotation period of a few days. The tidal force and the deviation from spherical structure are considered in lowest‐order perturbation theory. The magnetic field is in the form of toroidal magnetic flux rings, which are stored in mechanical equilibrium within the stably stratified overshoot region beneath the convection zone until the field has grown sufficiently strong for the undulatory instability to initiate the formation of rising loops. Frequencies and geometry of stable as well as growth rates of unstable eigenmodes are determined by linear stability analysis. Particular consideration is given to the question whether the effects of tidal forces and perturbations of the stellar structure can force a rising flux loop to enter the convection zone at specific longitudes.     

Towards an understanding of the Of?p star HD 191612: optical spectroscopy

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 (O6.5f?pe–O8fp). The Balmer and He  i lines show strongly variable emission which is highly reproducible on a well-determined 538-d period. He  ii absorptions and metal lines (including many selective emission lines but excluding He  ii λ4686 Å emission) are essentially constant in line strength, but are variable in velocity, establishing a double-lined binary orbit with   P _orb= 1542 d, e = 0.45  . We conduct a model-atmosphere analysis of the spectrum, and find that the system is consistent with a ∼O8 giant with a ∼B1 main-sequence secondary. Since the periodic 538-d changes are unrelated to orbital motion, rotational modulation of a magnetically constrained plasma is strongly favoured as the most likely underlying 'clock'. An upper limit on the equatorial rotation is consistent with this hypothesis, but is too weak to provide a strong constraint.     

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.     

Massive stars

We describe the present state of massive star research seen from the viewpoint of stellar evolution, with special emphasis on close binaries. Statistics of massive close binaries are reasonably complete for the Solar neighbourhood. We defend the thesis that within our knowledge, many scientific results where the effects of binaries are not included, have an academic value, but may be far from reality. In chapter I, we summarize general observations of massive stars where we focus on the HR diagram, stellar wind mass loss rates, the stellar surface chemistry, rotation, circumstellar environments, supernovae. Close binaries can not be studied separately from single stars and vice versa. First, the evolution of single stars is discussed (chapter I). We refer to new calculations with updated stellar wind mass loss rate formalisms and conclusions are proposed resulting from a comparison with representative observations. Massive binaries are considered in chapter II. Basic processes are briefly described, i.e. the Roche lobe overflow and mass transfer, the common envelope process, the spiral-in process in binaries with extreme mass ratio, the effects of mass accretion and the merging process, the implications of the (asymmetric) supernova explosion of one of the components on the orbital parameters of the binary. Evolutionary computations of interacting close binaries are discussed and general conclusions are drawn. The enormous amount of observational data of massive binaries is summarized. We separately consider the non-evolved and evolved systems. The latter class includes the semi-detached and contact binaries, the WR binaries, the X-ray binaries, the runaways, the single and binary pulsars. A general comparison between theoretical evolution and observations is combined with a discussion of specially interesting binaries: the evolved binaries HD 163181, HD 12323, HD 14633, HD 193516, HD 25638, HD 209481, Per and Sgr; the WR+OB binary V444 Cyg; the high mass X-ray binaries Vela X-1, Wray 977, Cyg X-1; the low mass X-ray binaries Her X-1 and those with a black hole candidate; the runaway Pup, the WR+compact companion candidates Cyg X-3, HD 50896 and HD 197406. We finally propose an overall evolutionary model of massive close binaries as a function of primary mass, mass ratio and orbital period. Chapter III deals with massive star population synthesis with a realistic population of binaries. We discuss the massive close binary frequency, mass ratio and period distribution, the observations that allow to constrain possible asymmetries during the supernova explosion of a massive star. We focuss on the comparison between observed star numbers (as a function of metallicity) and theoretically predicted numbers of stellar populations in regions of continuous star formation and in starburst regions. Special attention is given to the O-type star/WR star/red supergiant star population, the pulsar and binary pulsar population, the supernova rates. Received 17 July 1998     

Modelling the differential rotation of F stars

We model stellar differential rotation based on the mean-field theory of fluid dynamics. DR is mainly driven by Reynolds stress, which is anisotropic and has a non-diffusive component because the Coriolis force affects the convection pattern. Likewise, the convective heat transport is not strictly radial but slightly tilted towards the rotation axis, causing the polar caps to be slightly warmer than the equator. This drives a flow opposite to that caused by differential rotation and so allows the system to avoid the Taylor-Proudman state. Our model reproduces the rotation pattern in the solar convection zone and allows predictions for other stars with outer convection zones. The surface shear turns out to depend mainly on the spectral type and only weakly on the rotation rate. We present results for stars of spectral type F which show signs of very strong differential rotation in some cases. Stars just below the mass limit for outer convection zones have shallow convection zones with short convective turnover times. We find solar-type rotation and meridional flow patterns at much shorter rotation periods and horizontal shear much larger than on the solar surface, in agreement with recent observations. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Activity cycle of the giant star of Z Andromedae and its spin period

We have re-analysed the long-term optical light curve (LC) of the symbiotic star Z Andromedae, covering 112 yr of mostly visual observations. Two strictly periodic cycles and one quasi-periodic cycle can be identified in this LC. A   P 1 = 7550  d quasi-periodicity characterizes the repetition time of the outburst episodes of this symbiotic star. Six such events have been recorded so far. During quiescence states of the system, that is, in time-intervals between outbursts, the LC is clearly modulated by a stable coherent period of   P 2 = 759.1  d. This is the well-known orbital period of the Z Andromedae binary system that has been measured also spectroscopically. A third coherent period of   P 3 = 658.4  d is modulating the intense fluctuations in the optical brightness of the system during outbursts. We attribute the trigger of the outburst phenomenon and the clock that drives it, to a solar-type magnetic dynamo cycle that operates in the convection and the outer layers of the giant star of the system. We suggest that the intense surface activity of the giant star during maximum phases of its magnetic cycle is especially enhanced in one or two antipode regions, fixed in the atmosphere of the star and rotating with it. Such spots could be active regions around the North Pole and the South Pole of a general magnetic dipole field of the star. The P3 periodicity is half the beat of the binary orbital period of the system and the spin period of the giant. The latter is then either 482 or 1790 d. If only one pole is active on the surface of the giant, P3 is the beat period itself, and the spin period is 352 d. It could also be 5000 d if the giant is rotating in a retrograde direction. We briefly compare these findings in the LC of Z Andromedae to similar modulations that were identified in the LC of two other prototype symbiotics, BF Cyg and YY Her.     

The massive double-lined O-type binary HD 165052

We present a new optical spectroscopic study of the O-type binary HD 165052 based on high- and intermediate-resolution CCD observations. We re-investigated the spectral classification of the binary components, obtaining spectral types of O6.5 V and O7.5 V for the primary and secondary, respectively, finding that both stars display weak C  iii λ 5696 emission in their spectra. We also determined a radial-velocity orbit for HD 165052 with a period of  2.95510±0.00001 d  , and semi-amplitudes of 94.8 and  104.7±0.5 km s^-1  , resulting in a mass ratio   Q =0.9  . From a comparison with previous radial-velocity determinations, we found evidence of apsidal motion in the system. Several signatures of wind–wind collision, such as phase-locked variability of the X-ray flux and the Struve–Sahade effect, are also considered. It was also found that the reddening in the region should be normal, in contrast with previous determinations.     

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.     

Composite spectraPaper 9: HR 5983

HR 5983 is a sixth-magnitude composite-spectrum binary system. It has a circular orbit and a period of 108 d. We separate the component spectra and show that they are of types ∼G6 IIIa and A2.5 IV, with a mass ratio (giant:dwarf) of ∼1.04. We find that the secondary has a mass close to 2.6 M_⊙ and has already evolved somewhat from the main sequence, and that the primary appears to have reached the helium-burning 'blue loop' stage and is in synchronous rotation. The system resembles Capella in several respects, and may represent an earlier phase in the evolution of the latter.     

The nearest star of spectral type O3: a component of the multiple system HD 150136

From radial velocities determined in high signal-to-noise digital spectra, we report the discovery that the brightest component of the binary system HD 150136 is of spectral type O3. We also present the first double-lined orbital solution for this binary. Our radial velocities confirm the previously published spectroscopic orbital period of 2.6 d. He  ii absorptions appear double at quadratures, but single lines of N  v and N  iv visible in our spectra define a radial velocity orbit of higher semi-amplitude for the primary component than do the He  ii lines. From our orbital analysis, we obtain minimum masses for the binary components of 27 and  18 M_⊙  . The neutral He absorptions apparently do not follow the orbital motion of any of the binary components, thus they most probably arise in a third star in the system.     

The spectral variations of the O-type runaway supergiant HD 188209

We report spectral time series of the late O-type runaway supergiant HD 188209. Radial velocity variations of photospheric absorption lines with a possible quasi-period of ∼6.4 d have been detected in high-resolution echelle spectra. Night-to-night variations in the position and strength of the central emission reversal of the H α profile occurring over ill-defined time-scales have been observed. The fundamental parameters of the star are derived using state-of-the-art plane-parallel and unified non-LTE model atmospheres, the latter including the mass-loss rate. The derived helium abundance is moderately enhanced with respect to solar, and the stellar masses are lower than those predicted by the evolutionary models. The binary nature of this star is not suggested either from Hipparcos photometry or from radial velocity curves.     

WR 7a: a V Sagittae or a qWR star?

The star WR 7a, also known as SPH 2, has a spectrum that resembles that of V Sagittae stars although no O  vi emission has been reported. The Temporal Variance Spectrum – TVS – analysis of our data shows weak but strongly variable emission of O  vi lines which is below the noise level in the intensity spectrum.
Contrary to what is seen in V Sagittae stars, optical photometric monitoring shows very little, if any, flickering. We found evidence of periodic variability. The most likely photometric period is   P _phot= 0.227(±14) d  , while radial velocities suggest a period of   P _spec= 0.204(±13) d  . One-day aliases of these periods can not be ruled out. We call attention to similarities with HD 45166 and DI Cru (= WR 46), where multiple periods are present. They may be associated to the binary motion or to non-radial oscillations.
In contrast to a previous conclusion by Pereira et al., we show that WR 7a contains hydrogen. The spectrum of the primary star seems to be detectable as the N  v 4604 Å  absorption line is visible. If so, it means that the wind is optically thin in the continuum and that it is likely to be a helium main sequence star.
Given the similarity to HD 45166, we suggests that WR 7a may be a qWR – quasi Wolf–Rayet – star. Its classification is WN4h/CE in the Smith, Shara & Moffat three-dimensional classification system.     

High-resolution optical and infrared spectroscopic observations of Cir X-1

We present new optical and infrared (IR) observations of Cir X-1 taken near apastron. Both sets of spectra show asymmetric emission lines. Archival optical observations show that an asymmetric H α emission line has been in evidence for the past 20 years, although the shape of the line has changed significantly. We present an eccentric ( e ∼0.7–0.9) low-mass binary model, where the system consists of a neutron star orbiting around a (sub)giant companion star of 3–5 M_⊙. We suggest that the broad components of the emission lines arise in a high-velocity, optically thick flow near the neutron star, while the narrow components of the optical and the IR lines arise near the companion star and a heated ejecta shell surrounding the binary respectively. In this model, the velocity of the narrow component reflects the space velocity of the binary; the implied radial velocity (+430 km s⁻¹ after correcting for Galactic rotation) is the highest velocity known for an X-ray binary.     

Differential rotation of cool active stars: the case of intermediate rotators

In this paper, we present a new method for measuring the surface differential rotation of cool stars with rotation periods of a few days, for which the sparse phase coverage achievable from single-site observations generally prevents the use of more conventional techniques. The basic idea underlying this new analysis is to obtain the surface differential rotation pattern that minimizes the information content of the reconstructed Doppler image through a simultaneous fit of all available data.
Simulations demonstrate that the performance of this new method in the case of cool stars is satisfactory for a variety of observing strategies. Differential rotation parameters can be recovered reliably as long as the total data set spans at least 4 per cent of the time for the equator to lap the pole by approximately one complete cycle. We find in particular that these results hold for potentially complex spot distributions (as long as they include a mixture of low- and high-latitude features), and for various stellar inclination angles and rotation velocities. Such measurements can be obtained from either unpolarized or polarized data sets, provided their signal-to-noise ratio is larger than approximately 500 and 5000 per 2 km s⁻¹ spectral bin, respectively.
This method should therefore be very useful for investigating differential rotation in a much larger sample of objects than what has been possible up to now, and should hence give us the opportunity of studying how differential rotation reacts to various phenomena operating in stellar convective zones, such as tidal effects or dynamo magnetic field generation.