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.     

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.     

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.     

A search for Ca ii emission enhancement in stars resulting from nearby giant planets

We present a search for periodicities ( P _chr) in the chromospheric Ca  ii infrared triplet emission of several stars ( τ Boo, 51 Peg, υ And, ρ ¹ Cnc, ρ CrB, 70 Vir and GL 876) which may be directly attributable to interaction with close-in giant planets. Activity enhancements could arise from increased non-radiative heating and dynamo action in planet-induced tidal bulges (with P _chr≈ P _orb/2), or from interactions between the stellar and planetary magnetic fields (with P _chr≈ P _orb). We compare both P _chr and the phase dependence of the activity with the planetary orbital period P _orb, the orbital phase, and models. No significant P _chr or phase dependence attributable to planets can be clearly identified. We place approximate upper limits on the amplitude of any planet-induced activity. We identify a possible stellar rotation period for GL 876, and support previous period determinations for four other stars. We discuss the results and possible directions of future research.     

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.     

UBV photometry, UV spectroscopy and radio observations of the peculiar binary V Sagittae

We have performed high-speed UBV photometric observations on the peculiar binary V Sagittae. Using three new eclipse timings we update the orbital ephemeris and convert it to a dynamical time-scale (TDB). We also searched for quasi-periodic oscillations but did not detect them. Using the Wilson–Devinney algorithm we have modelled the light curve to find the stellar parameters of V Sge. We find that the system is a detached binary but that the primary star is very close to filling its Roche lobe, while the secondary star fills 90 per cent of its Roche lobe volume. We find temperatures of the primary and the secondary star to be T ₁=41 000 K and T ₂=22 000 K. We find i =72° and masses of 0.8 M_⊙ and 3.3 M_⊙ for the primary and secondary stars respectively. De-archived Hubble Space Telescope ( HST ) spectroscopy of V Sge shows evidence of mass loss via a wind or winds. In addition we report radio observations of V Sge during an optical high state at 2 cm, 3.6 cm and 6 cm wavelengths. The 3.6 cm emission is increased by a factor of more than six compared with an earlier detection in a previous optical high state.     

LO Pegasi: an investigation of multiband optical polarization

We present BVR polarimetric study of the cool active star LO Pegasi (LO Peg) for the first time. LO Peg was found to be highly polarized among the cool active stars. Our observations yield average values of polarization in LO Peg:   P_B = 0.387 ± 0.004 per cent, θ_B= 88°± 1°; P_V = 0.351 ± 0.004 per cent, θ_V= 91°± 1°  and   P_R = 0.335 ± 0.003 per cent, θ_R= 91°± 1°  . Both the degree of polarization and the position angle are found to be variable. The semi-amplitude of the polarization variability in B, V and R bands is found to be  0.18 ± 0.02, 0.13 ± 0.01  and  0.10 ± 0.02  per cent, respectively. We suggest that the levels of polarization observed in LO Peg could be the result of scattering of an anisotropic stellar radiation field by an optically thin circumstellar envelope or scattering of the stellar radiation by prominence-like structures.     

Burial of the polar magnetic field of an accreting neutron star – II. Hydromagnetic stability of axisymmetric equilibria

The theory of polar magnetic burial in accreting neutron stars predicts that a mountain of accreted material accumulates at the magnetic poles of the star, and that, as the mountain spreads equatorward, it is confined by, and compresses, the equatorial magnetic field. Here, we extend previous, axisymmetric, Grad–Shafranov calculations of the hydromagnetic structure of a magnetic mountain up to accreted masses as high as   M _a= 6 × 10⁻⁴ M_⊙  , by importing the output from previous calculations (which were limited by numerical problems and the formation of closed bubbles to   M _a < 10⁻⁴ M_⊙  ) into the time-dependent, ideal-magnetohydrodynamic code zeus-3d and loading additional mass on to the star dynamically. The rise of buoyant magnetic bubbles through the accreted layer is observed in these experiments. We also investigate the stability of the resulting hydromagnetic equilibria by perturbing them in zeus-3d . Surprisingly, it is observed that the equilibria are marginally stable for all   M _a≤ 6 × 10⁻⁴ M_⊙  ; the mountain oscillates persistently when perturbed, in a combination of Alfvén and acoustic modes, without appreciable damping or growth, and is therefore not disrupted (apart from a transient Parker instability initially, which expels <1 per cent of the mass and magnetic flux).     

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.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

ROSAT All-Sky Survey X-ray and EUV observations of YY Gem and AU Mic

Simultaneous X-ray and extreme ultraviolet (EUV) ( ROSAT XRT and WFC All-Sky Survey) observations of the highly active dMe flare stars YY Gem and AU Mic show that the two stars displayed an unusual type of flaring behaviour. We detect several X-ray and EUV flares superimposed on an enhanced and smoothly varying quiescent background. The two large impulsive-type X-ray flares on YY Gem reach peak X-ray luminosities of     and we estimate that they had similar integrated luminosities (∼6–8×10³³ erg). AU Mic also produced several X-ray and EUV flares, with one very impulsive flare producing a 10-fold increase in XRT count rate. This flare was even larger than the YY Gem flares (peak L _X of     and integrated L _X of    
The     ratio for both stars is at the 'saturation' limit found in rapidly rotating dwarfs and the most active RS CVn stars. We suggest that the gradually varying components are the result of a period of continuous, unresolved flaring activity. Alternatively, they may be the result of the emergence and subsequent decay of a new magnetic active region on the stellar surface of these stars.     

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.     

Long Secondary Periods in variable red giants

We present a study of a sample of Large Magellanic Cloud red giants exhibiting Long Secondary Periods (LSPs). We use radial velocities obtained from VLT spectral observations and MACHO and OGLE light curves to examine properties of the stars and to evaluate models for the cause of LSPs. This sample is much larger than the combined previous studies of Hinkle et al. and Wood, Olivier & Kawaler.
Binary and pulsation models have enjoyed much support in recent years. Assuming stellar pulsation, we calculate from the velocity curves that the typical fractional radius change over an LSP cycle is greater than 30 per cent. This should lead to large changes in T _eff that are not observed. Also, the small light amplitude of these stars seems inconsistent with the radius amplitude. We conclude that pulsation is not a likely explanation for the LSPs. The main alternative, physical movement of the star – binary motion – also has severe problems. If the velocity variations are due to binary motion, the distribution of the angle of periastron in our large sample of stars has a probability of  1.4 × 10⁻³  that it comes from randomly aligned binary orbits. In addition, we calculate a typical companion mass of  0.09 M_⊙  . Less than 1 per cent of low-mass main-sequence stars have companions near this mass  (0.06–0.12 M_⊙)  whereas ∼25–50 per cent of low-mass red giants end up with LSPs. We are unable to find a suitable model for the LSPs and conclude by listing their known properties.     

VLT spectroscopy of XTE J2123−058 during quiescence: the masses of the two components

We present Very Large Telescope (VLT) low-resolution spectroscopy of the neutron star X-ray transient XTE J2123−058 during its quiescent state. Our data reveal the presence of a K7V companion which contributes 77 per cent to the total flux at λ 6300 and orbits the neutron star at     . Contrary to other soft X-ray transients (SXTs), the H _α emission is almost exactly in antiphase with the velocity curve of the optical companion. Using the light-centre technique we obtain     and hence     This, combined with a previous determination of the inclination angle     yields     and     . M ₂ agrees well with the observed spectral type. Doppler tomography of the H α emission shows a non-symmetric accretion disc distribution mimicking that seen in SW Sex stars. Although we find a large systemic velocity of −     this value is consistent with the galactic rotation velocity at the position of J2123−058, and hence a halo origin. The formation scenario of J2123−058 is still unresolved.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Accretion in stellar clusters and the initial mass function

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses M ∗≳0.5 M_⊙ . The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a ρ  ∝  R ⁻² distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of γ =−3/2 (where Salpeter is γ =−2.35) . Once the stars dominate the potential and are virialized, which occurs first in the cluster core, the accretion radius is the Bondi–Hoyle radius. The resultant mass spectrum has an asymptotic limit of γ =−2 with slightly steeper slopes ( γ ≈−2.5) if the stars are already mass-segregated. Simulations of accretion on to clusters containing 1000 stars show that, as expected, the low-mass stars accumulate the majority of their masses during the gas-dominated phase whereas the high-mass stars accumulate the majority of their masses during the stellar-dominated phase. This results in a mass spectrum with a relatively shallow γ ≈3/2 power law for low-mass stars and a steeper power law for high-mass stars −2.5≲ γ ≤−2 . This competitive accretion model also results in a mass-segregated cluster.     

Self-enrichment by asymptotic giant branch stars in globular clusters: comparison between intermediate and high metallicities

We present theoretical evolutionary sequences of intermediate-mass stars  ( M = 3 − 6.5 M_⊙)  with metallicity   Z = 0.004  . Our goal is to test whether the self-enrichment scenario by massive asymptotic giant branch stars may work for the high-metallicity globular clusters, after previous works by the same group showed that the theoretical yields by this class of objects can reproduce the observed trends among the abundances of some elements, namely the O–Al and O–Na anticorrelations, at intermediate metallicities, i.e.  [Fe/H]=−1.3  . We find that the increase in the metallicity favours only a modest decrease of the luminosity and the temperature at the bottom of the envelope for the same core mass, and also the efficiency of the third dredge-up is scarcely altered. On the contrary, differences are found in the yields, due to the different impact that processes with the same efficiency have on the overall abundance of envelopes with different metallicities. We expect the same qualitative patterns as in the intermediate-metallicity case, but the slopes of some of the relationships among the abundances of some elements are different. We compare the sodium–oxygen anticorrelation for clusters of intermediate metallicity ( Z ≈ 10⁻³) and clusters of metallicity large as in these new models. Although the observational data are still too scarce, the models are consistent with the observed trends, provided that only stars of   M ≳ 5 M_⊙  contribute to self-enrichment.     

Bulk viscosity of mixed nucleon–hyperon–quark matter in neutron stars

We calculate the coefficient of bulk viscosity by considering the non-leptonic weak interactions in the cores of hybrid stars with both hyperons and quarks. We first determine the dependence of the production rate of neutrons on the reaction rate of quarks in the non-leptonic processes, that is,  Γ _n = K _s Γ _s +Γ_Λ+ 2Γ_Σ−  . The conversion rate,   K _s   , in our scenario is a complicated function of baryon number density. We also consider the medium effect of quark matter on bulk viscosity. Using these results, we estimate the limiting rotation of the hybrid stars, which may suppress the r-mode instability more effectively. Hybrid stars should be the candidates for the extremely rapid rotators.