Constraints on the magnetic field geometry of magnetars

We study the effect of the magnetic field geometry on the oscillation spectra of strongly magnetized stars. We construct a configuration of magnetic field where a toroidal component is added to the standard poloidal one. We consider a star with a type I superconductor core so that both components of the magnetic field are expelled from the core and confined in the crust. Our results show that the toroidal contribution does not influence significantly the torsional oscillations of the crust. On the contrary, the confinement of the magnetic field in the crust drastically affects the torsional oscillation spectrum. A comparison with estimations for the magnetic field strength, from observations, excludes the possibility that magnetars will have a magnetic field solely confined in the crust, that is, our results suggest that the magnetic field in whatever geometry has to permeate the whole star.     

The magnetosphere of oscillating neutron stars in general relativity

Just as a rotating magnetized neutron star has material pulled away from its surface to populate a magnetosphere, a similar process can occur as a result of neutron-star pulsations rather than rotation. This is of interest in connection with the overall study of neutron star oscillation modes but with a particular focus on the situation for magnetars. Following a previous Newtonian analysis of the production of a force-free magnetosphere in this way Timokhin et al., we present here a corresponding general-relativistic analysis. We give a derivation of the general relativistic Maxwell equations for small-amplitude arbitrary oscillations of a non-rotating neutron star with a generic magnetic field and show that these can be solved analytically under the assumption of low current density in the magnetosphere. We apply our formalism to toroidal oscillations of a neutron star with a dipole magnetic field and find that the low current density approximation is valid for at least half of the oscillation modes, similarly to the Newtonian case. Using an improved formula for the determination of the last closed field line, we calculate the energy losses resulting from toroidal stellar oscillations for all of the modes for which the size of the polar cap is small. We find that general relativistic effects lead to shrinking of the size of the polar cap and an increase in the energy density of the outflowing plasma. These effects act in opposite directions but the net result is that the energy loss from the neutron star is significantly smaller than suggested by the Newtonian treatment.     

How to drive roAp stars

Rapidly oscillating Ap stars constitute a unique class of pulsators with which to study non-radial oscillations under some — even for stars — unusual physical conditions. These stars are chemically peculiar, they have strong magnetic fields and they often pulsate in several high-order acoustic modes simultaneously. We discuss here an excitation mechanism for short-period oscillation modes based on the classical κ mechanism. We particularly stress the conditions that must be fulfilled for successful driving. Specifically, we discuss the roles of the chemical peculiarity and strong magnetic field on the oscillation modes and what separates these pulsators from δ Scuti and Am-type stars.     

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.     

Equilibrium configurations of strongly magnetized neutron stars with realistic equations of state

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy, FPS, Shen and LS, employing the Tomimura–Eriguchi scheme to construct the equilibrium configurations. We study the basic physical properties of the sequences in the framework of Newtonian gravity. In addition, we take a new step by taking into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g. structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyse this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than  10¹⁵ g cm⁻³  . The maximum baryon mass of the magnetized stars with axis ratio   q ∼ 0.7  increases about up to 20 per cent for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.     

Alfvén polar oscillations of relativistic stars

We study polar Alfvén oscillations of relativistic stars endowed with a strong global poloidal dipole magnetic field. Here, we focus only on the axisymmetric oscillations which are studied by numerically evolving the two-dimensional perturbation equations. Our study shows that the spectrum of the polar Alfvén oscillations is discrete in contrast to the spectrum of axial Alfvén oscillations which is continuous. We also show that the typical fluid modes, such as the f and p modes, are not significantly affected by the presence of the strong magnetic field.     

Atmospheres and spectra of strongly magnetized neutron stars

We construct atmosphere models for strongly magnetized neutron stars with surface fields     and effective temperatures     . The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars, including radio pulsars, soft gamma-ray repeaters, and anomalous X-ray pulsars. In our models, the atmosphere is composed of pure hydrogen or helium and is assumed to be fully ionized. The radiative opacities include free–free absorption and scattering by both electrons and ions computed for the two photon polarization modes in the magnetized electron–ion plasma. Since the radiation emerges from deep layers in the atmosphere with     , plasma effects can significantly modify the photon opacities by changing the properties of the polarization modes. In the case where the magnetic field and the surface normal are parallel, we solve the full, angle-dependent, coupled radiative transfer equations for both polarization modes. We also construct atmosphere models for general field orientations based on the diffusion approximation of the transport equations and compare the results with models based on full radiative transport. In general, the emergent thermal radiation exhibits significant deviation from blackbody, with harder spectra at high energies. The spectra also show a broad feature     around the ion cyclotron resonance     , where Z and A are the atomic charge and atomic mass of the ion, respectively; this feature is particularly pronounced when     . Detection of the resonance feature would provide a direct measurement of the surface magnetic fields on magnetars.     

Crustal oscillations of slowly rotating relativistic stars

We study low-amplitude crustal oscillations of slowly rotating relativistic stars consisting of a central fluid core and an outer thin solid crust. We estimate the effect of rotation on the torsional toroidal modes and on the interfacial and shear spheroidal modes. The results compared against the Newtonian ones for wide range of neutron star models and equations of state.     

On mode conversion and wave reflection in magnetic Ap stars

We investigate the effect of a strong large-scale magnetic field on the reflection of high-frequency acoustic modes in rapidly oscillating Ap stars. To that end, we consider a toy model composed of an isothermal atmosphere matched on to a polytropic interior and determine the numerical solution to the set of ideal magnetohydrodynamic equations in a local plane-parallel approximation with constant gravity. Using the numerical solution in combination with approximate analytical solutions that are valid in the limits where the magnetic and acoustic components are decoupled, we calculate the relative fraction of energy flux that is carried away in each oscillation cycle by running acoustic waves in the atmosphere and running magnetic waves in the interior. For oscillation frequencies above the acoustic cut-off, we show that most energy losses associated with the presence of running waves occur in regions where the magnetic field is close to vertical. Moreover, by considering the depth dependence of the energy associated with the magnetic component of the wave in the atmosphere we show that a fraction of the wave energy is kept in the oscillation every cycle. For frequencies above the acoustic cut-off frequency, such energy is concentrated in regions where the magnetic field is significantly inclined in relation to the local vertical. Even though our calculations were aimed at studying oscillations with frequencies above the acoustic cut-off frequency, based on our results we discuss what results may be expected for oscillations of lower frequency.     

Self-similar pressure oscillations in neutron star envelopes as probes of neutron star structure

We study eigenmodes of acoustic oscillations of high multipolarity l ∼ 100–1000 and high frequency (∼100 kHz), localized in neutron star envelopes. We show that the oscillation problem is self-similar. Once the oscillation spectrum is calculated for a given equation of state (EOS) in the envelope and given stellar mass M and radius R , it can be rescaled to a star with any M and R (but the same EOS in the envelope). For l ≳ 300, the modes can be subdivided into the outer and inner ones. The outer modes are mainly localized in the outer envelope. The inner modes are mostly localized near the neutron drip point, being associated with the softening of the EOS after the neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and accreted matter and show that the spectra of the inner modes are essentially different. A detection and identification of high-frequency pressure modes would allow one to infer M and R and determine also the EOS in the envelope (accreted or ground state) providing a new, potentially powerful method to explore the main parameters and internal structure of neutron stars.     

Seismic signatures of strange stars with crust

We study acoustic oscillations (eigenfrequencies, velocity distributions, damping times) of normal crusts of strange stars. These oscillations are very specific because of huge density jump at the interface between the normal crust and the strange matter core. The oscillation problem is shown to be self-similar. For a low (but non-zero) multipolarity l , the fundamental mode (without radial nodes) has a frequency of ∼300 Hz and mostly horizontal oscillation velocity; other pressure modes have frequencies ≳20 kHz and almost radial oscillation velocities. The latter modes are similar to radial oscillations (having approximately the same frequencies and radial velocity profiles). The oscillation spectrum of strange stars with crust differs from the spectrum of neutron stars. If detected, acoustic oscillations would allow one to discriminate between strange stars with crust and neutron stars and constrain the mass and radius of the star.     

Magnetic fields and rotation in white dwarfs and neutron stars

Recent spectropolarimetric observations of Ap and Bp stars with improved sensitivity have suggested that most Ap and Bp stars are magnetic with dipolar fields of at least a few hundred gauss. These new estimates suggest that the range of magnetic fluxes found for the majority of magnetic white dwarfs is similar to that of main-sequence Ap–Bp stars, thus strengthening the empirical evidence for an evolutionary link between magnetism on the main sequence and magnetism in white dwarfs. We draw parallels between the magnetic white dwarfs and the magnetic neutron stars and argue that the observed range of magnetic fields in isolated neutron stars  ( B_p ∼ 10¹¹–10¹⁵ G)  could also be explained if their mainly O-type progenitors have effective dipolar fields in the range of a few gauss to a few kilogauss, assuming approximate magnetic flux conservation with the upper limit being consistent with the recent measurement of a field of   B_p ∼ 1100 G  for θ Orion C.
In the magnetic field–rotation diagram, the magnetic white dwarfs can be divided into three groups of different origin: a significant group of strongly magnetized slow rotators  ( P _rot∼ 50 –100 yr)  that have originated from single-star evolution, a group of strongly magnetized fast rotators  ( P _rot∼ 700 s)  , typified by EUVE J0317–853, that have originated from a merger, and a group of modest rotators ( P _rot∼ hours–days) of mixed origin (single-star and CV-type binary evolution). We propose that the neutron stars may similarly divide into distinct classes at birth , and suggest that the magnetars may be the counterparts of the slowly rotating high-field magnetic white dwarfs.     

Gaussian process modelling of asteroseismic data

The measured properties of stellar oscillations can provide powerful constraints on the internal structure and composition of stars. To begin this process, oscillation frequencies must be extracted from the observational data, typically time series of the star's brightness or radial velocity. In this paper, a probabilistic model is introduced for inferring the frequencies and amplitudes of stellar oscillation modes from data, assuming that there is some periodic character to the oscillations, but that they may not be exactly sinusoidal. Effectively, we fit damped oscillations to the time series, and hence the mode lifetime is also recovered. While this approach is computationally demanding for large time series (>1500 points), it should at least allow improved analysis of observations of solar-like oscillations in subgiant and red giant stars, as well as sparse observations of semiregular stars, where the number of points in the time series is often low. The method is demonstrated on simulated data and then applied to radial velocity measurements of the red giant star  ξ Hydrae  , yielding a mode lifetime between 0.41 and 2.65 d with 95 per cent posterior probability. The large frequency separation between modes is ambiguous, however we argue that the most plausible value is 6.3 μHz, based on the radial velocity data and the star's position in the Hertzsprung–Russell diagram.     

The influence of Hall drift on the magnetic fields of neutron stars

We consider the evolution of magnetic fields under the influence of Hall drift and Ohmic decay. The governing equation is solved numerically, in a spherical shell with   r _i / r _o = 0.75  . Starting with simple free-decay modes as initial conditions, we then consider the subsequent evolution. The Hall effect induces so-called helicoidal oscillations, in which energy is redistributed among the different modes. We find that the amplitude of these oscillations can be quite substantial, with some of the higher harmonics becoming comparable with the original field. Nevertheless, this transfer of energy to the higher harmonics is not sufficient to accelerate significantly the decay of the original field, at least not at the   R _B = O (100)  parameter values accessible to us, where this Hall parameter   R _B   measures the ratio of the Ohmic time-scale to the Hall time-scale. We do find clear evidence though of increasingly fine structures developing for increasingly large   R _B   , suggesting that perhaps this Hall-induced cascade to ever-shorter length-scales is eventually sufficiently vigorous to enhance the decay of the original field. Finally, the implications for the evolution of neutron star magnetic fields are discussed.     

The sixth frequency of roAp star HR 1217

The frequency of the sixth mode observed in HR 1217 cannot currently be understood in the light of the asymptotic theory for high-frequency acoustic modes, valid for adiabatic pulsations in spherically symmetric stars. Deviations from spherical symmetry will change the frequency of the observed oscillations away from the frequencies predicted by the asymptotic theory. Among other things, the presence of a magnetic field in the star will introduce such deviations from spherical symmetry. In the present work we inspect the effect of the magnetic field on the oscillations of HR 1217. Although we model HR 1217 in a rather simplified manner, we believe, as argued in the text, that the essence of the conclusions would not be modified if a more sophisticated model were to be used instead. The results show that the frequency of the sixth mode observed in HR 1217 might be understood if the effect of the magnetic field on the frequencies of the oscillations is fully taken into account. Moreover, we interpret the sixth frequency as a mode of even degree, and explain the absence from the observations of an intermediate l =1 , m =0 mode as the result of large energy losses associated with the coupling of the latter with Alfvénic waves, in the strongly magnetized outer layers of the star. Finally, we show how, in the light of this explanation, the frequency of the sixth mode might be used in order to put constraints on the parameters characterizing the star.     

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.     

Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests

We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes.     

Multisite campaign on the open cluster M67 – II. Evidence for solar-like oscillations in red giant stars

Measuring solar-like oscillations in an ensemble of stars in a cluster, holds promise for testing stellar structure and evolution more stringently than just fitting parameters to single field stars. The most-ambitious attempt to pursue these prospects was by Gilliland et al. who targeted 11 turn-off stars in the open cluster M67 (NGC 2682), but the oscillation amplitudes were too small (<20 μmag) to obtain unambiguous detections. Like Gilliland et al. we also aim at detecting solar-like oscillations in M67, but we target red giant stars with expected amplitudes in the range 50–  500 μmag  and periods of 1 to 8 h. We analyse our recently published photometry measurements, obtained during a six-week multisite campaign using nine telescopes around the world. The observations are compared with simulations and with estimated properties of the stellar oscillations. Noise levels in the Fourier spectra as low as  27 μmag  are obtained for single sites, while the combined data reach  19 μmag  , making this the best photometric time series of an ensemble of red giant stars. These data enable us to make the first test of the scaling relations (used to estimate frequency and amplitude) with an homogeneous ensemble of stars. The detected excess power is consistent with the expected signal from stellar oscillations, both in terms of its frequency range and amplitude. However, our results are limited by apparent high levels of non-white noise, which cannot be clearly separated from the stellar signal.     

A frequency analysis of the p-mode pulsations of the Ap SrEu(Cr) star HD 119027

We present new high-speed, multisite photometric observations of the rapidly oscillating Ap star HD 119027 acquired over seven nights during 1996. A frequency analysis of these observations reveals the presence of oscillations at 1835, 1875, 1888, 1913, 1940, 1942 and (possibly) 1953 μHz. These frequencies are consistent with a spacing of either 13 or 26 μHz, depending on the reality of the oscillations at 1875 and 1953 μHz. The data in hand do not permit us to discriminate between the two possible spacings. If the smaller value of the spacing is correct, it suggests that HD 119027 is outside the main-sequence band. Two of the frequencies listed above are separated by only 1.95 μHz, suggesting that they are modes of ( n ,ℓ) and ( n  − 1, ℓ + 2), which in roAp stars is a quantity governed by the internal magnetic field.