The effect of r-mode instability on the evolution of isolated strange stars

We studied the evolution of isolated strange stars (SSs) synthetically, considering the influence of r -mode instability. Our results show that the cooling of SSs with non-ultrastrong magnetic fields is delayed by heating due to r -mode damping for millions of years, while the spin-down of the stars is dominated by gravitational radiation (GR). Especially for the SSs in a possible existing colour–flavour locked (CFL) phase, the effect of r -mode instability on the evolution of stars becomes extremely important because the viscosity, neutrino emissivity and specific heat involving pairing quarks are blocked. It leads to the cooling of these colour superconducting stars being very slow and the stars can remain at high temperature for millions of years, which differs completely from previous understanding. In this case, an SS in CFL phase can be located at the bottom of its r -mode instability window for a long time, but does not spin-down to a very low frequency for hours.     

Heating and cooling of magnetars with accreted envelopes

We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (   B ≳ 10¹⁴  G) non-accreted and accreted outermost envelopes composed of different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.     

Thermal evolution of rotating hybrid stars

As a neutron star spins down, the nuclear matter is continuously converted into quark matter due to the core density increase, and then latent heat is released. We have investigated the thermal evolution of neutron stars undergoing such deconfinement phase transition. We have taken into account the conversion in the frame of the general theory of relativity. The released energy has been estimated as a function of changed rate of deconfinement baryon number. The numerical solutions to the cooling equation are seen to be very different from those without the heating effect. The results show that neutron stars may be heated to higher temperatures which is well matched with pulsar's data despite the onset of fast cooling in neutron stars with quark matter cores. It is also found that the heating effect has a magnetic field strength dependence. This feature could be particularly interesting for high temperatures of low-field millisecond pulsars at a later stage. The high temperature could fit the observed temperature for PSR J0437−4715.     

Minimal models of cooling neutron stars with accreted envelopes

We study the 'minimal' cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and enhanced cooling is produced by neutrino emission due to the Cooper pairing of neutrons. Extending our recent previous work, we include the effects of surface accreted envelopes of light elements. We employ the phenomenological density-dependent critical temperatures   T _cp(ρ)  and   T _cnt(ρ)  of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature   T _cns(ρ)  of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of our recent previous work and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations.     

Cooling of quark stars in the colour superconductive phase: effect of photons from glueball decay

The cooling history of a quark star in the colour superconductive phase is investigated. Here we specifically focus on the two-flavour colour (2SC) phase where the novel process of photon generation via glueball (GLB) decay has already been investigated. The picture we present here can, in principle, be generalized to quark stars entering a superconductive phase where similar photon generation mechanisms are at play. As much as 10⁴⁵–10⁴⁷ erg of energy is provided by the GLB decay in the 2SC phase. The generated photons slowly diffuse out of the quark star, keeping it hot and radiating as a blackbody (with possibly a Wien spectrum in gamma-rays) for millions of years. We discuss hot radio-quiet isolated neutron stars in our picture (such as RX J185635–3754 and RX J0720.4–3125) and argue that their nearly blackbody spectra (with a few broad features) and their remarkably tiny hydrogen atmosphere are indications that these might be quark stars in the colour superconductive phase where some sort of photon generation mechanism (reminiscent of the GLB decay) has taken place. Fits to observed data of cooling compact stars favour models with superconductive gaps of  Δ_2SC∼ 15–35 MeV  and densities  ρ_2SC= (2.5–3.0) ×ρ_N  (ρ_N being the nuclear matter saturation density) for quark matter in the 2SC phase. If correct, our model combined with more observations of isolated compact stars could provide vital information to studies of quark matter and its exotic phases.     

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.     

Thermal relaxation in young neutron stars

The internal properties of the neutron star crust can be probed by observing the epoch of thermal relaxation. After the supernova explosion, powerful neutrino emission quickly cools the stellar core, while the crust stays hot. The cooling wave then propagates through the crust, as a result of its finite thermal conductivity. When the cooling wave reaches the surface (age 10–100 yr) , the effective temperature drops sharply from 250 eV to 30 or 100 eV, depending on the cooling model. The crust relaxation time is sensitive to the (poorly known) microscopic properties of matter of subnuclear density, such as the heat capacity, thermal conductivity, and superfluidity of free neutrons. We calculate the cooling models with the new values of the electron thermal conductivity in the inner crust, based on a realistic treatment of the shapes of atomic nuclei. Superfluid effects may shorten the relaxation time by a factor of 4. The comparison of theoretical cooling curves with observations provides a potentially powerful method of studying the properties of the neutron superfluid and highly unusual atomic nuclei in the inner crust.     

Almost analytic models of ultramagnetized neutron star envelopes

Recent ROSAT measurements show that the X-ray emission from isolated neutron stars is modulated at the stellar rotation period. To interpret these measurements, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We present nearly analytic models of the thermal structure of the envelopes of ultramagnetized neutron stars. Specifically, we examine the limit in which only the ground Landau level is filled. We use the models to estimate the amplitude of modulation expected from non-uniformities in the surface temperatures of strongly magnetized neutron stars. In addition, we estimate cooling rates for stars with fields B  ∼ 10¹⁵ − 10¹⁶ G, which are relevant to models that invoke 'magnetars' to account for soft γ-ray emission from some repeating sources.     

Crust–core coupling and r-mode damping in neutron stars: a toy model

r-modes in neutron stars with crusts are damped by viscous friction at the crust–core boundary. The magnitude of this damping, evaluated by Bildsten & Ushomirsky (BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency for neutron stars spun up by accretion in low-mass X-ray binaries (LMXBs). In this paper we explore the mechanical coupling between the core r-modes and the elastic crust, using a toy model of a constant-density neutron star having a crust with a constant shear modulus. We find that, at spin frequencies in excess of ≈50 Hz, the r-modes strongly penetrate the crust. This reduces the relative motion (slippage) between the crust and the core compared with the rigid-crust limit. We therefore revise down, by as much as a factor of 10²–10³ , the damping rate computed by BU, significantly reducing the maximal possible spin frequency of neutron stars with solid crusts. The dependence of the crust–core slippage on the spin frequency is complicated, and is very sensitive to the physical thickness of the crust. If the crust is sufficiently thick, the curve of the critical spin frequency for the onset of the r-mode instability becomes multivalued for some temperatures; this is related to avoided crossings between the r-mode and higher-order torsional modes in the crust. The critical frequencies are comparable to the observed spins of neutron stars in LMXBs and millisecond pulsars.     

The Na i D resonance lines in main-sequence late-type stars

We study the sodium D lines (D1: 5895.92 Å; D2: 5889.95 Å) in late-type dwarf stars. The stars have spectral types between F6 and M5.5 ( B − V between 0.457 and 1.807) and metallicity between  [Fe/H]=−0.82  and 0.6. We obtained medium-resolution echelle spectra using the 2.15-m telescope at the Argentinian observatory Complejo Astronómico El Leoncito (CASLEO). The observations have been performed periodically since 1999. The spectra were calibrated in wavelength and in flux. A definition of the pseudo-continuum level is found for all our observations. We also define a continuum level for calibration purposes. The equivalent width of the D lines is computed in detail for all our spectra and related to the colour index ( B − V ) of the stars. When possible, we perform a careful comparison with previous studies. Finally, we construct a spectral index  ( R '_D)  as the ratio between the flux in the D lines and the bolometric flux. We find that, once corrected for the photospheric contribution, this index can be used as a chromospheric activity indicator in stars with a high level of activity. Additionally, we find that combining some of our results, we obtain a method to calibrate in flux stars of unknown colour.     

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment.
The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to  ∼3000 km s⁻¹  ) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.     

Photometric modelling of starspots – I. A Barnes–Evans-like surface brightness–colour relation using (Ic−K)

In the first part of this work, the empirical correlation of stellar surface brightness F _V with ( I _c− K ) broad-band colour is investigated by using a sample of stars cooler than the Sun. A bilinear correlation is found to represent well the brightness of G, K and M giant stars. The change in slope occurs at ( I _c− K )∼2.1 or at about the transition from K to M spectral types. The same relationship is also investigated for dwarf stars and found to be distinctly different from that of the giants. The dwarf star correlation differs by an average of −0.4 in ( I _c− K ) or by a maximum in F _V of ∼−0.1, positioning it below that of the giants, with both trends tending towards convergence for the hotter stars in our sample. The flux distribution derived from the F _V −( I _c− K ) relationship for the giant stars, together with that derived from an F _V −( V − K ) relationship and the blackbody flux distribution, is then utilized to compute synthetic light V and colour ( V − R )_c, ( V − I )_c and ( V − K ) curves of cool spotted stars. We investigate the effects on the amplitudes of the curves by using these F _V –colour relations and by assuming the effective gravity of the spots to be lower than the gravity of the unspotted photosphere. We find that the amplitudes produced by using the F _V −( I _c− K ) relationship are larger than those produced by the other two brightness correlations, meaning smaller and/or warmer spots.     

A young cooling neutron star in the remnant of Supernova 1987A

We describe the cooling theory for isolated neutron stars that are several tens of years old. Their cooling differs greatly from the cooling of older stars that has been well studied in the literature. It is sensitive to the physics of the inner stellar crust and even to the thermal conductivity of the stellar core, which is never important at later cooling stages. The absence of observational evidence for the formation of a neutron star during the explosion of Supernova 1987A is consistent with the fact that the star was actually born there. It may still be hidden in the dense center of the supernova remnant. If, however, the star is not hidden, then it should have a low thermal luminosity (below ～10³⁴ erg s^?1) and a short internal thermal relaxation time (shorter than 13 yr). This requires that the star undergo intense neutrino cooling (e.g., via the direct Urca process) and have a thin crust with strong superfluidity of free neutrons and/or an anomalously high thermal conductivity.     

Cooling of Accretion-Heated Neutron Stars

We present a brief, observational review about the study of the cooling behaviour of accretion-heated neutron stars and the inferences about the neutron-star crust and core that have been obtained from these studies. Accretion of matter during outbursts can heat the crust out of thermal equilibrium with the core and after the accretion episodes are over, the crust will cool down until crust-core equilibrium is restored. We discuss the observed properties of the crust cooling sources and what has been learned about the physics of neutron-star crusts. We also briefly discuss those systems that have been observed long after their outbursts were over, i.e, during times when the crust and core are expected to be in thermal equilibrium. The surface temperature is then a direct probe for the core temperature. By comparing the expected temperatures based on estimates of the accretion history of the targets with the observed ones, the physics of neutron-star cores can be investigated. Finally, we discuss similar studies performed for strongly magnetized neutron stars in which the magnetic field might play an important role in the heating and cooling of the neutron stars.     

Unveiling the thermal and magnetic map of neutron star surfaces though their X-ray emission: method and light-curve analysis

Recent Chandra and XMM–Newton observations of a number of X-ray 'dim' pulsating neutron stars have revealed quite unexpected features in the emission from these sources. Their soft thermal spectrum, believed to originate directly from the star surface, shows evidence for a phase-varying absorption line at some hundred eVs. The pulse modulation is relatively large (pulsed fractions in the range ∼12–35 per cent), the pulse shape is often non-sinusoidal, and the hard X-ray colour appears to be anticorrelated in phase with the total emission. Moreover, the prototype of this class, RX J0720.4−3125, has been found to undergo rather sensible changes in both its spectral and timing properties over a time-scale of a few years. All these new findings seem difficult to reconcile with the standard picture of a cooling neutron star endowed with a purely dipolar magnetic field, at least if surface emission is produced in an atmosphere on top of the crust. In this paper we explore how a dipolar+quadrupolar star-centred field influences the properties of the observed light curves. The phase-resolved spectrum has been evaluated accounting for both radiative transfer in a magnetized atmosphere and general relativistic ray-bending. We computed over 78 000 light curves, varying the quadrupolar components and the viewing geometry. A comparison of the data with our model indicates that higher-order multipoles are required to reproduce the observations.     

NLTE models of line-driven stellar winds – III. Influence of X-ray radiation on wind structure of O stars

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional artificial source of X-rays, assumed to originate in the wind shocks.
We show that the influence of shock X-ray emission on wind mass-loss rate is relatively small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray sources, especially for stars cooler than   T _eff≲ 35 000 K  .
We discuss the origin of the   L _X/ L ∼ 10⁻⁷  relation. For stars with thick wind this relation can be explained assuming that the cooling time depends on wind density. Stars with optically thin winds exhibiting the 'weak wind problem' display enhanced X-ray emission which may be connected with large shock cooling length. We propose that this effect can explain the 'weak wind problem'.
Inclusion of X-rays leads to a better agreement of the model ionization structure with observations. However, we do not find any significant influence of X-rays on P  v ionization fraction implying that the presence of X-rays cannot explain the P  v problem.
We study the implications of modified ionization equilibrium due to shock emission on the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like ions may be affected by the line absorption within the cool wind.     

On the dynamics of superfluid neutron star cores

We discuss the nature of the various modes of pulsation of superfluid neutron stars using comparatively simple Newtonian models and the Cowling approximation. The matter in these stars is described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). In our model, we incorporate the non-dissipative interaction known as the entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. We show that there is no independent set of pulsating g-modes in a non-rotating superfluid neutron star core, even though the linearized superfluid equations contain a well-defined (and real-valued) analogue to the so-called Brunt–Väisälä frequency. Instead, what we find are two sets of spheroidal perturbations whose nature is predominately acoustic. In addition, an analysis of the zero-frequency subspace (i.e. the space of time-independent perturbations) reveals two sets of degenerate spheroidal perturbations, which we interpret to be the missing g-modes, and two sets of toroidal perturbations. We anticipate that the degeneracy of all these zero-frequency modes will be broken by the Coriolis force in the case of rotating stars. To illustrate this we consider the toroidal pulsation modes of a slowly rotating superfluid star. This analysis shows that the superfluid equations support a new class of r-modes, in addition to those familiar from, for example, geophysical fluid dynamics. Finally, the role of the entrainment effect on the superfluid mode frequencies is shown explicitly via solutions to dispersion relations that follow from a 'local' analysis of the linearized superfluid equations.     

Masses and luminosities of O‐ and B‐type stars and red supergiants

Massive stars are of interest as progenitors of supernovae, i.e. neutron stars and black holes, which can be sources of gravitational waves. Recent population synthesis models can predict neutron star and gravitational wave observations but deal with a fixed supernova rate or an assumed initial mass function for the population of massive stars. Here we investigate those massive stars, which are supernova progenitors, i.e. with O‐ and early B‐type stars, and also all supergiants within 3 kpc. We restrict our sample to those massive stars detected both in 2MASS and observed by Hipparcos, i.e. only those stars with parallax and precise photometry. To determine the luminosities we calculated the extinctions from published multi‐colour photometry, spectral types, luminosity class, all corrected for multiplicity and recently revised Hipparcos distances. We use luminosities and temperatures to estimate the masses and ages of these stars using different models from different authors. Having estimated the luminosities of all our stars within 3 kpc, in particular for all O‐ and early B‐type stars, we have determined the median and mean luminosities for all spectral types for luminosity classes I, III, and V. Our luminosity values for supergiants deviate from earlier results: Previous work generally overestimates distances and luminosities compared to our data, this is likely due to Hipparcos parallaxes (generally more accurate and larger than previous ground‐based data) and the fact that many massive stars have recently been resolved into multiples of lower masses and luminosities. From luminosities and effective temperatures we derived masses and ages using mass tracks and isochrones from different authors. From masses and ages we estimated lifetimes and derived a lower limit for the supernova rate of ≈20 events/Myr averaged over the next 10 Myr within 600 pc from the sun. These data are then used to search for areas in the sky with higher likelihood for a supernova or gravitational wave event (like OB associations) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical study of the r-mode instability of rapidly rotating nascent neutron stars

The first results of numerical analysis of classical r-modes of rapidly rotating compressible stellar models are reported. The full set of linear perturbation equations of rotating stars in Newtonian gravity is solved numerically without the slow rotation approximation. A critical curve of gravitational wave emission induced instability, which restricts the rotational frequencies of hot young neutron stars, is obtained. Taking the standard cooling mechanisms of neutron stars into account, we also show the 'evolutionary curves' along which neutron stars are supposed to evolve as cooling and spinning down proceed. Rotational frequencies of 1.4-M_⊙ stars suffering from this instability decrease to around 100 Hz when the standard cooling mechanism of neutron stars is employed. This result confirms the results of other authors, who adopted the slow rotation approximation.     

100-yr mass-loss modulations on the asymptotic giant branch

We analyse the differences in infrared circumstellar dust emission between oxygen-rich Mira and non-Mira stars, and find that they are statistically significant. In particular, we find that these stars segregate in the K–[12] versus [12]–[25] colour–colour diagram, and have distinct properties of the IRAS LRS spectra, including the peak position of the silicate emission feature. We show that the infrared emission from the majority of non-Mira stars cannot be explained within the context of standard steady-state outflow models.
The models can be altered to fit the data for non-Mira stars by postulating non-standard optical properties for silicate grains, or by assuming that the dust temperature at the inner envelope radius is significantly lower (300–400 K) than typical silicate grain condensation temperatures (800–1000 K) . We argue that the latter is more probable and provide detailed model fits to the IRAS LRS spectra for 342 stars. These fits imply that two-thirds of non-Mira stars and one-third of Mira stars do not have hot dust (>500 K) in their envelopes.
The absence of hot dust can be interpreted as a recent (∼100 yr) decrease in the mass-loss rate. The distribution of best-fitting model parameters agrees with this interpretation and strongly suggests that the mass loss resumes on similar time-scales. Such a possibility appears to be supported by a number of spatially resolved observations (e.g. recent Hubble Space Telescope images of the multiple shells in the Egg Nebula) and is consistent with new dynamical models for mass loss on the asymptotic giant branch.