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.     

Relativistically expanding cylindrical electromagnetic fields

We study relativistically expanding electromagnetic fields of cylindrical geometry. The fields emerge from the side surface of a cylinder and are invariant under translations parallel to the axis of the cylinder. The expansion velocity is in the radial direction and is parametrized by   v = R /( ct )  . We consider force-free magnetic fields by setting the total force the electromagnetic field exerts on the charges and the currents equal to zero. Analytical and semi-analytical separable solutions are found for the relativistic problem. In the non-relativistic limit, the mathematical form of the equations is similar to equations that have already been studied in static systems of the same geometry.     

Toroidal versus poloidal magnetic fields in Sun-like stars: a rotation threshold

From a set of stellar spectropolarimetric observations, we report the detection of surface magnetic fields in a sample of four solar-type stars, namely HD 73350, HD 76151, HD 146233 (18 Sco) and HD 190771. Assuming that the observed variability of polarimetric signal is controlled by stellar rotation, we establish the rotation periods of our targets, with values ranging from 8.8 d (for HD 190771) to 22.7 d (for HD 146233). Apart from rotation, fundamental parameters of the selected objects are very close to the Sun's, making this sample a practical basis to investigate the specific impact of rotation on magnetic properties of Sun-like stars.
We reconstruct the large-scale magnetic geometry of the targets as a low-order  (ℓ < 10)  spherical harmonic expansion of the surface magnetic field. From the set of magnetic maps, we draw two main conclusions. (i) The magnetic energy of the large-scale field increases with rotation rate. The increase in chromospheric emission with the mean magnetic field is flatter than observed in the Sun. Since the chromospheric flux is also sensitive to magnetic elements smaller than those contributing to the polarimetric signal, this observation suggests that a larger fraction of the surface magnetic energy is stored in large scales as rotation increases. (ii) Whereas the magnetic field is mostly poloidal for low rotation rates, more rapid rotators host a large-scale toroidal component in their surface field. From our observations, we infer that a rotation period lower than ≈12 d is necessary for the toroidal magnetic energy to dominate over the poloidal component.     

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.     

A threshold for sequentially self-propagating star formation

When a stellar wind bubble expands into an homogeneous medium, there are two possible outcomes. This is due to the fact that the self-gravity of the swept-up shell acts in two orthogonal directions: tangentially , to promote fragmentation of the shell, and radially , to decelerate expansion of the shell. The outcome depends on whether self-gravity works faster in the tangential or the radial direction.
If the wind power ℒ_o is large and the effective isothermal sound speed a _o in the swept-up gas is small – approximately         – tangential self-gravity works faster. A thin dense shell is swept up and fragments while it is still expanding supersonically. This is the scenario often invoked to explain sequentially self-propagating star formation.
However, if ℒ_o is small and/or a _o is large, radial self-gravity works faster. Expansion of the bubble stalls before the shell can fragment. The expansion speed ceases to be supersonic, the outer shock dissipates, and the shell is neither thin nor dense. Under this circumstance, the shell is unlikely to fragment and star formation will not propagate sequentially.
These conclusions are probably not altered significantly when the medium into which the wind blows is inhomogeneous, provided that the mean density on opposite sides of the bubble does not differ by many orders of magnitude.     

On the magnetic structure and wind parameter profiles of Alfvén wave driven winds in late-type supergiant stars

Cool stars at giant and supergiant evolutionary phases present low-velocity and high-density winds, responsible for the observed high mass-loss rates. Although presenting high luminosities, radiation pressure on dust particles is not sufficient to explain the wind acceleration process. Among the possible solutions to this still unsolved problem, Alfvén waves are, probably, the most interesting for their high efficiency in transfering energy and momentum to the wind. Typically, models of Alfvén wave driven winds result in high-velocity winds if they are not highly damped. In this work, we determine self-consistently the magnetic field geometry and solve the momentum, energy and mass conservation equations, to demonstrate that even a low-damped Alfvén wave flux is able to reproduce the low-velocity wind. We show that the magnetic flux tubes expand with a super-radial factor of S > 30 near the stellar surface, larger than that used in previous semi-empirical models. The rapid expansion results in a strong spatial dilution of the wave flux. We obtained the wind parameter profiles for a typical supergiant star of  16 M_⊙  . The wind is accelerated in a narrow region, coincident with the region of high divergence of the magnetic field lines, up to 100 km s⁻¹. For the temperature, we obtained a slight decrease near the surface for low-damped waves, because the wave heating mechanism is less effective than the radiative losses. The peak temperature occurs at   r ≃ 1.5  r ₀  reaching 6000 K. Propagating outwards, the wind cools down mainly due to adiabatic expansion.     

Magnetic cycle of LQ Hydrae: observational indications and dynamo model

We present a model for the differential rotation and dynamo activity of the young rapidly rotating K0 dwarf LQ Hya ( P _rot=1.6 d). As might be expected from observations of the similar rapid rotator AB Dor, the predicted differential rotation is small (≃0.8 per cent) but extremely efficient in generating magnetic fields. The dynamo, which is of a distributed type, produces a globally axisymmetric field with radial and azimuthal components that are of the same magnitude and display a phase-lag in their evolution of about π/2. This is consistent with the long-term Zeeman–Doppler imaging study by Donati. The latitudinal distribution of flux is, however, a little different from that observed and the cycle period of 3.2 yr is somewhat shorter than suggested by the observations.     

On the quasi-periodic oscillations in magnetars

We study torsional Alfvén oscillations of magnetars, that is neutron stars with a strong magnetic field. We consider the poloidal and toroidal components of the magnetic field and a wide range of equilibrium stellar models. We use a new coordinate system  ( X , Y )  , where     and     and a ₁ is the radial component of the magnetic field. In this coordinate system, the one+two-dimensional evolution equation describing the quasi-periodic oscillations (QPOs), see Sotani et al., is reduced to a one+one-dimensional equation where the perturbations propagate only along the y -axis. We solve the one+one-dimensional equation for different boundary conditions and the open magnetic field lines, that is magnetic field lines that reach the surface and there match up with the exterior dipole magnetic field as well as closed magnetic lines, i.e. magnetic lines that never reach the stellar surface. For the open field lines, we find two families of QPO frequencies: a family of 'lower' QPO frequencies which is located near the x -axis and a family of 'upper' frequencies located near the y -axis. According to Levin, the fundamental frequencies of these two families can be interpreted as the turning point of the continuous spectrum. We find that the upper frequencies are multiples of the lower ones by a constant equalling  2 n + 1  . For the closed lines, the corresponding factor is   n + 1  . By using these relations, we can explain both the lower and the higher observed frequencies in SGR 1806−20 and SGR 1900+14.     

The pulsating hydrogen-deficient star LSS 3184 (BX Cir)

We present extensive photometry and spectroscopy of the extremely hydrogen-deficient star, LSS 3184, recently discovered to be a rapid variable (period ∼0.1066 d) strikingly similar to V652 Her. Over 95 h of photometry confirms the reported variability, which is of rather low amplitude (Δ V ∼0.03 mag), defines the light curve with greater precision and establishes a much more accurate ephemeris (period ∼0.106 578 4 d) to form a basis for detecting period change. Attention is drawn to the usefulness of a period-finding technique that fits harmonic components to the photometric observations. Spectroscopy shows a peak-to-peak variation in radial velocity of about 30 km s⁻¹, which, when combined with the photometric observations, confirms the pulsational nature of the variability and strongly indicates that the pulsations are radial in nature.     

Multiperiodic variations in the last 104-yr light curve of the symbiotic star BF Cyg

We analyse a light curve (LC) of the symbiotic star BF Cyg, covering 114 yr of its photometric history. The star had a major outburst around the year 1894. Since then the mean optical brightness of the system is in steady decline, reaching only in the last few years its pre-outburst value. Superposed on this general decline are some six less intense outbursts of 1–2 mag and duration of 2000–5000 d. We find a cycle of 6376 d, or possibly twice this period, in the occurrence of these outbursts. We suggest that the origin of the system outbursts is in some magnetic cycle in the outer layers of the giant star of the system, akin to the less intense 8000-d magnetic cycle of our Sun. We further find, that in addition to its well-known binary period of 757.3 d, BF Cyg possesses also another photometric period of 798.8 d. This could be the rotation period of the giant star of the system. If it is, the beat period of these two periodicities, 14 580 d, is the rotation period of a tidal wave on the surface of the giant. A fourth period of 4436 d, the beat period of the 14 580-d and the 6376-d cycles is possibly also present in the LC. We predict that BF Cyg will be at the peak of its next outburst around the month of May in the year 2007. The newly discovered 798.8-d period explains the disappearance of the orbital modulation at some epochs in the LC. The 757.3-d oscillations will be damped again around the year 2013.     

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.     

Magnetic deformation of the white dwarf surface structure

The influence of strong, large‐scale magnetic fields on the structure and temperature distribution in white dwarf atmospheres is investigated. Magnetic fields may provide an additional component of pressure support, thus possibly inflating the atmosphere compared to the non‐magnetic case. Since the magnetic forces are not isotropic, atmospheric properties may significantly deviate from spherical symmetry. In this paper the magnetohydrostatic equilibrium is calculated numerically in the radial direction for either for small deviations from different assumptions for the poloidal current distribution. We generally find indication that the scale height of the magnetic white dwarf atmosphere enlarges with magnetic field strength and/or poloidal current strength. This is in qualitative agreement with recent spectropolarimetric observations of Grw+10°8247. Quantitatively, we .nd for e.g. a mean surface poloidal field strength of 100 MG and a toroidal field strength of 2‐10 MG an increase of scale height by a factor of 10. This is indicating that already a small deviation from the initial force‐free dipolar magnetic field may lead to observable effects. We further propose the method of finite elements for the solution of the two‐dimensional magnetohydrostatic equilibrium including radiation transport in the diffusive approximation. We present and discuss preliminary solutions, again indicating on an expansion of the magnetized atmosphere.     

The β Cephei instability strip

We carried out a series of linear stability analyses of the radial and low-degree non-radial p modes for stellar models with initial masses of     . The stellar models were computed by using convective overshoot distance     , 0.25 and 0.40  H _P. Our numerical results show that the β Cephei instability strip forms a horn-shaped region pointing upwards near the main sequence on the Hertzsprung–Russell diagram (HRD). The lower part of the instability strip for the radial modes join the zero-age main-sequence (ZAMS) at     , while the top of the instability strip extends up to     . The instability strip for the non-radial modes is even wider. The overall instability strip is dominated by the radial and non-radial fundamental modes. The first overtone (the radial-order index     is also pulsationally unstable. We have shown that the β Cephei stability is almost independent of the overshoot parameter d _over used for the stellar models, while it depends critically on the metal abundance. With decreasing metal abundance, the instability region shrinks and eventually disappears for     .     

Unstable non-radial modes in radial pulsators: theory and an example

We study the possibility of the excitation of non-radial oscillations in classical pulsating stars. The stability of an RR Lyrae model is examined through non-adiabatic non-radial calculations. We also explore stability in the presence of non-linear coupling between radial and non-radial modes of nearly identical frequency.   In our model, a large number of unstable low-degree (ℓ = 1,2) modes have frequencies in the vicinity of unstable radial mode frequencies. The growth rates of such modes, however, are considerably smaller than those of the radial modes. We also recover an earlier result that at higher degrees (ℓ = 6–12) there are modes trapped in the envelope with growth rates similar to those of radial modes.   Subsequently, monomode radial pulsation of this model is considered. The destabilizing effect of the 1:1 resonance between the radial mode and nearby non-radial modes of low degrees is studied, with the assumption that the excited radial mode saturates the linear instability of all other modes. The instability depends on the radial mode amplitude, the frequency difference, the damping rate of the non-radial mode, and the strength of the non-linear coupling between the modes considered. At the pulsation amplitudes typical for RR Lyrae stars, the instability of the monomode radial pulsation and the concomitant resonant excitation of some non-radial oscillation modes is found to be very likely.     

Discovery and asteroseismological analysis of a new pulsating DB white dwarf star, PG 2246+121

We report 36.6 h of time-resolved CCD photometry of the DB white dwarf star PG 2246+121 and the discovery that it is a new pulsating variable. Analysis of our compact single-site data set allowed the detection of three mode multiplets, two triplets at 256 and 329 s, respectively, and one doublet at 286 s. The frequency splitting within those structures is exactly the same within the length and accuracy of our data set.
We argue that these multiplets are the result of non-radial g-mode pulsations, most probably of spherical degree ℓ=1, which then yields a formal stellar rotation period of 2.00±0.12 d. We suggest that the excited modes are three consecutive radial overtones of order 3–7, most likely k =4,5,6. This discovery's impact on the understanding of pulsating DB white dwarfs is discussed.     

Mode switching in the nearby Mira-like variable R Doradus

We discuss visual observations spanning nearly 70 years of the nearby semiregular variable R Doradus. Using wavelet analysis, we show that the star switches back and forth between two pulsation modes having periods of 332 d and about 175 d, the latter with much smaller amplitude. Comparison with model calculations suggests that the two modes are the first and third radial overtone, with the physical diameter of the star making fundamental-mode pulsation unlikely. The mode changes occur on a time-scale of about 1000 d, which is too rapid to be related to a change in the overall thermal structure of the star and may instead be related to weak chaos.   The Hipparcos distance to R Dor is 62.4 ± 2.8 pc which, taken with its dominant 332-d period, places it exactly on the period–luminosity (P–L) relation of Miras in the Large Magellanic Cloud. Our results imply first-overtone pulsation for all Miras which fall on the P–L relation. We argue that semiregular variables with long periods may largely be a subset of Miras and should be included in studies of Mira behaviour. The semiregulars may contain the immediate evolutionary Mira progenitors, or stars may alternate between periods of semiregular and Mira behaviour.     

Dimensionless measures of turbulent magnetohydrodynamic dissipation rates

The magnetic Reynolds number, R _M, is defined as the product of a characteristic scale and associated flow speed divided by the microphysical magnetic diffusivity. For laminar flows, R _M also approximates the ratio of advective to dissipative terms in the total magnetic energy equation, but for turbulent flows this latter ratio depends on the energy spectra and approaches unity in a steady state. To generalize for flows of arbitrary spectra we define an effective magnetic dissipation number,   R _M,e  , as the ratio of the advection to microphysical dissipation terms in the total magnetic energy equation, incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, and distinct pairs of inner and outer scales for magnetic and kinetic spectra. As expected, for a substantial parameter range   R _M,e∼ O (1) ≪ R _M  . We also distinguish   R _M,e  from     where the latter is an effective magnetic Reynolds number for the mean magnetic field equation when a turbulent diffusivity is explicitly imposed as a closure. That   R _M,e  and     approach unity even if   R _M≫ 1  highlights that, just as in hydrodynamic turbulence, energy dissipation of large-scale structures in turbulent flows via a cascade can be much faster than the dissipation of large-scale structures in laminar flows. This illustrates that the rate of energy dissipation by magnetic reconnection is much faster in turbulent flows, and much less sensitive to microphysical reconnection rates compared to laminar flows.     

A non-adiabatic analysis for axisymmetric pulsations of magnetic stars

This paper presents the results of a non-adiabatic analysis for axisymmetric non-radial pulsations including the effect of a dipole magnetic field. Convection is assumed to be suppressed in the stellar envelope, and the diffusion approximation is used to radiative transport. As in a previous adiabatic analysis, the eigenfunctions are expanded in a series of spherical harmonics. The analysis is applied to a  1.9-M_⊙  , main-sequence model  (log  T _eff= 3.913)  . The presence of a magnetic field always stabilizes low-order acoustic modes. All the low-order modes of the model that are excited by the κ-mechanism in the He  ii ionization zone in the absence of a magnetic field are found to be stabilized if the polar strength of the dipole magnetic field is larger than about 1 kG. For high-order p modes, on the other hand, distorted dipole and quadrupole modes excited by the κ-mechanism in the H ionization zone remain overstable, even in the presence of a strong magnetic field. It is found, however, that all the distorted radial high-order modes are stabilized by the effect of the magnetic field. Thus, our non-adiabatic analysis suggests that distorted dipole modes and distorted quadrupole modes are most likely excited in rapidly oscillating Ap stars. The latitudinal amplitude dependence is found to be in reasonable agreement with the observationally determined one for HR 3831. Finally, the expected amplitude of magnetic perturbations at the surface is found to be very small.     

Relativistic models of magnetars: the twisted torus magnetic field configuration

We find general relativistic solutions of equilibrium magnetic field configurations in magnetars, extending previous results of Colaiuda et al. Our method is based on the solution of the relativistic Grad–Shafranov equation, to which Maxwell's equations can be reduced. We obtain equilibrium solutions with the toroidal magnetic field component confined into a finite region inside the star, and the poloidal component extending to the exterior. These so-called twisted torus configurations have been found to be the final outcome of dynamical simulations in the framework of Newtonian gravity, and appear to be more stable than other configurations. The solutions include higher-order multipoles, which are coupled to the dominant dipolar field. We use arguments of minimal energy to constrain the ratio of the toroidal to the poloidal field.     

Dynamo-generated magnetic fields at the surface of a massive star

Spruit has shown that an astrophysical dynamo can operate in the non-convective material of a differentially rotating star as a result of a particular instability in the magnetic field (the Tayler instability). By assuming that the dynamo operates in a state of marginal instability, Spruit has obtained formulae which predict the equilibrium strengths of azimuthal and radial field components in terms of local physical quantities. Here, we apply Spruit's formulae to our previously published models of rotating massive stars in order to estimate Tayler dynamo field strengths. There are no free parameters in Spruit's formulae. In our models of 10- and  50-M_⊙  stars on the zero-age main sequence, we find internal azimuthal fields of up to 1 MG, and internal radial components of a few kG. Evolved models contain weaker fields. In order to obtain estimates of the field strength at the stellar surface, we examine the conditions under which the Tayler dynamo fields are subject to magnetic buoyancy. We find that conditions for Tayler instability overlap with those for buoyancy at intermediate to high magnetic latitudes. This suggests that fields emerge at the surface of a massive star between magnetic latitudes of about 45° and the poles. We attempt to estimate the strength of the field which emerges at the surface of a massive star. Although these estimates are very rough, we find that the surface field strengths overlap with values which have been reported recently for line-of-sight fields in several O and B stars.