Pycnonuclear burning of 34Ne in accreting neutron stars

We study the pycnonuclear burning of ³⁴Ne in the inner crust of an accreting neutron star. We show that the associated energy production rate can be calculated analytically for any arbitrary temporal variability of the mass accretion rate. We argue that the theoretical time-scale for ³⁴Ne burning is currently very uncertain and ranges from a fraction of a millisecond to a few years. The fastest allowable burning may change the composition of the accreted crust while the slowest burning leads to a time-independent nuclear energy generation rate for a variable accretion. The results are important for constructing self-consistent models of the accreted crust and deep crustal heating in neutron stars which enter soft X-ray transients.     

Diamagnetic screening of the magnetic field in accreting neutron stars

A possible mechanism for screening of the surface magnetic field of an accreting neutron star, by the accreted material, is investigated. We model the material flow in the surface layers of the star by an assumed two-dimensional velocity field satisfying all the physical requirements. Using this model velocity we find that, in the absence of magnetic buoyancy, the surface field is screened (i.e. there is submergence of the field by advection) within the time-scale of material flow of the top layers. On the other hand, if magnetic buoyancy is present, the screening happens over a time-scale that is characteristic of the slower flow of the deeper (and hence, denser) layers. For accreting neutron stars, this longer time-scale turns out to be about 10⁵ yr, which is of a similar order of magnitude to the accretion time-scale of most massive X-ray binaries.     

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.     

Spectral evolution of magnetic flares and time lags in accreting black hole sources

We present a model for the short time-scale spectral variability of accreting black holes. It describes the time-averaged spectra well, and also temporal characteristics such as the power-density spectrum, time/phase lags, and coherence function of Cygnus X-1. We assume that X/ γ -rays are produced in compact magnetic flares at radii ≲100 GM c ² from the central black hole. The tendency for magnetic loops to inflate and detach from the underlying accretion disc causes the spectrum of a flare to evolve from soft to hard because of the decrease of the feedback from the cold disc, so causing time delays between hard and soft photons. We identify the observed time lags with the evolution time-scales of the flares, which are of the order of the Keplerian time-scale. We model the overall temporal variability using a pulse avalanche model in which each flare has a certain probability of triggering a neighbouring flare, thus occasionally producing long avalanches. The duration of the avalanches determines the Fourier frequencies at which most of the power emerges.     

A Chandra X-ray observation of the globular cluster Terzan 1

We present a ∼19-ks Chandra Advanced CCD Imaging Spectrometer (ACIS)-S observation of the globular cluster Terzan 1. 14 sources are detected within 1.4 arcmin of the cluster centre with two of these sources predicted to be not associated with the cluster (background active galactic nuclei or foreground objects). The neutron star X-ray transient, X1732−304, has previously been observed in outburst within this globular cluster with the outburst seen to last for at least 12 yr. Here, we find four sources that are consistent with the ROSAT position for this transient, but none of the sources are fully consistent with the position of a radio source detected with the Very Large Array that is likely associated with the transient. The most likely candidate for the quiescent counterpart of the transient has a relatively soft spectrum and an unabsorbed 0.5–10 keV luminosity of  2.6 × 10³² erg s⁻¹  , quite typical of other quiescent neutron stars. Assuming standard core cooling, from the quiescent flux of this source we predict long (>400 yr) quiescent episodes to allow the neutron star to cool. Alternatively, enhanced core cooling processes are needed to cool down the core. However, if we do not detect the quiescent counterpart of the transient this gives an unabsorbed 0.5–10 keV luminosity upper limit of  8 × 10³¹ erg s⁻¹  . We also discuss other X-ray sources within the globular cluster. From the estimated stellar encounter rate of this cluster we find that the number of sources we detect is significantly higher than expected by the relationship of Pooley et al.     

Emergence of magnetic field due to spin-polarized baryon matter in neutron stars

A model of the ferromagnetic origin of magnetic fields of neutron stars is considered. In this model, the magnetic phase transition occurs inside the core of neutron stars soon after formation. However, owing to the high electrical conductivity the core magnetic field is initially fully screened. We study how this magnetic field emerges for an outside observer. After some time, the induced field that screens the ferromagnetic field decays enough to uncover a detectable fraction of the ferromagnetic field. We calculate the time-scale of decay of the screening field and study how it depends on the size of the ferromagnetic core. We find that the same fractional decay of the screening field occurs earlier for larger cores. We conjecture that weak fields of millisecond pulsars, B ∼10⁸–10⁹ G, could be identified with ferromagnetic fields of unshielded fraction ε ∼10⁻⁴–10⁻³ resulting from the decay of screening fields by a factor 1− ε in ∼10⁸ yr since their birth.     

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.     

Repetitive structure in the stellar wind of HD 93843: a normal O-type star

We present the results from a 28-day IUE time-series campaign monitoring the stellar wind of the O5-type giant HD 93843. The principal aim was to study variability in the wind of a star with a normal projected rotation velocity. Systematic changes are identified, amidst continuous line-profile variability, in the absorption troughs of the Si  iv and N  v resonance lines. The patterns observed have characteristic time-scales of several days and are mimicked by fluctuations (of several 100 km s⁻¹) in the blue wings of the saturated C  iv P Cygni profile.   Fourier analysis provides support for the repeatability of wind structures in HD 93843 on a 7.1-d 'period'. Power at this frequency is evident only at intermediate and high velocities (i.e., above ∼0.3 of the terminal velocity). The long modulation time-scale suggests that changes in the star itself probably provide the physical source for triggering the onset of wind structure. Unfortunately the rotational, photometric, pulsational and magnetic properties of HD 93843 are too poorly constrained or known to permit a more detailed interpretation of the 7.1-d wind modulation in terms of potential inhomogeneities at the stellar surface. Nevertheless, our study demonstrates that the incidence of cyclic, possibly regular, stellar-wind variability is not restricted to rapid rotators. Comparisons with other OB stars which have exhibited repetitive wind changes on 'periods' of several days suggest that the time-dependent UV properties of HD 93843 are more akin to those of the O4-type supergiant ζ Puppis.     

Flare variability in the close binary FL Vir

The results of U -filter flare monitoring of the binary flare star FL Vir = Wolf 424 is presented. 57 flares with energies between  2 × 10²⁸  and  2 × 10³¹ erg  were recorded in 20 h of observation. The properties of flare occurrence and flare time-scales are analysed, and the flare activity level in 1980 April is determined to be   L _f( U ) = 8.0 × 10²⁶ erg s⁻¹  . This is larger than previously published results and may indicate a variation in the flare activity level on a time-scale of years. An analysis of existing data indicates that the flare activity level correlates with the relative orbital positions of the stars.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

On the origin of high-velocity runaway stars

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary  50–100 M_⊙  star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of  ≳200–400 km s⁻¹  (typical of pulsars), while  3–4 M_⊙  stars can attain velocities of  ≳300–400 km s⁻¹  (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.     

The physical parameters of Markarian 501 during flaring activity

We analyse an N -body simulation of the Small Magellanic Cloud (SMC), that of Gardiner & Noguchi, to determine its microlensing statistics. We find that the optical depth owing to self-lensing in the simulation is low, 0.4×10⁻⁷, but still consistent (at the 90 per cent level) with that observed by the EROS and MACHO collaborations. This low optical depth is due to the relatively small line-of-sight thickness of the SMC produced in the simulation. The proper motions and time-scales of the simulation are consistent with those observed assuming a standard mass function for stars in the SMC. The time-scale distribution from the standard mass function generates a significant fraction of short time-scale events: future self-lensing events towards the SMC may have the same time-scales as events observed towards the Large Magellanic Cloud (LMC). Although some debris was stripped from the SMC during its collision with the LMC about 2×10⁸ yr ago, the optical depth of the LMC owing to this debris is low, a few ×10⁻⁹, and thus cannot explain the measured optical depth towards the LMC.     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

Heavy-element abundances in seven SC stars and several related stars

We employ spectra of resolution 20–35000 of seven SC stars, four S stars, two Ba stars and two K–M stars to derive abundances of a variety of elements from Sr to Eu relative to iron. Special attention is paid to Rb and Tc, and to the ratio of the heavy s-process species to the light s-process elements. Abundances are derived in LTE, both by using model atmospheres in which the carbon and oxygen abundances are nearly equal and by using curves of growth. Spectrum synthesis is used for critical lines such as the 5924-Å line of Tc and the 7800-Å line of Rb. For most of the heavy-element stars the enhancement of the s-process elements is about a factor of 10. The ratio of the heavy to light s-process species is not far from solar, except for RR Her for which the same ratio is +0.45 dex. For Tc the blending by other lines is severe. While we have probably detected the 5924-Å line, we can only present abundances in the less-than-or-equal-to category. For Rb, whose abundance is sensitive to the ⁸⁵Rb/⁸⁷Rb ratio and hence to the neutron density during s-process production, we find a considerable range of abundances, indicating a neutron density from 10⁶ to ≳10⁸ cm⁻³ for the SC stars. For the four S stars the range is from 10⁷ to ≳10⁸ cm⁻³. Recent calculations by Gallino et al. show that neutron densities near 10⁷ cm⁻³ favour the ¹³C source for neutrons, while densities greater than 10⁸ cm⁻³ may be associated with neutrons from the ²²Ne source.     

Fourier-resolved energy spectra of the Narrow-Line Seyfert 1 Mkn 766

We compute Fourier-resolved X-ray spectra of the Seyfert 1 Markarian 766 to study the shape of the variable components contributing to the 0.3–10 keV energy spectrum and their time-scale dependence. The fractional variability spectra peak at 1–3 keV, as in other Seyfert 1 galaxies, consistent with either a constant contribution from a soft excess component below 1 keV and Compton reflection component above 2 keV or variable warm absorption enhancing the variability in the 1–3 keV range. The rms spectra, which show the shape of the variable components only, are well described by a single power law with an absorption feature around 0.7 keV, which gives it an apparent soft excess. This spectral shape can be produced by a power law varying in normalization, affected by an approximately constant (within each orbit) warm absorber, with parameters similar to those found by Turner et al. for the warm-absorber layer covering all spectral components in their scattering scenario  [ N _H∼ 3 × 10²¹ cm⁻², log(ξ) ∼ 1]  . The total soft excess in the average spectrum can therefore be produced by a combination of constant warm absorption on the power-law plus an additional less variable component. On shorter time-scales, the rms spectrum hardens and this evolution is well described by a change in power-law slope, while the absorption parameters remain the same. The frequency dependence of the rms spectra can be interpreted as variability arising from propagating fluctuations through an extended emitting region, whose emitted spectrum is a power law that hardens towards the centre. This scenario reduces the short time-scale variability of lower energy bands making the variable spectrum harder on shorter time-scales and at the same time explains the hard lags found in these data by Markowitz et al.     

Dust dynamics in dense molecular cores

We show that in a quiescent, dense pre-stellar core, exposed to the average interstellar radiation field, radiation pressure can cause the dust to migrate inwards, relative to the gas, on a time-scale of a few megayears – and faster if the radiation field is stronger than average. This has two potentially important effects.
First, there is an increase in the abundance of dust relative to gas in the inner parts of the core, and hence also in the efficiency of gas-cooling by dust. The increased cooling efficiency predisposes these regions to dynamical collapse and star formation. Additionally, it predisposes them to fragmentation, particularly if – as seems likely – the dust enhancements are stochastic and inhomogeneous, due to anisotropy of the incident radiation field and/or to directing of the migration by the local magnetic field. It also increases the metallicities of the resulting stars, and hence presumably the likelihood of planet formation in their accretion discs.
Secondly, there is a steepening of the optical-depth profile, especially at those impact parameters b where the visual optical depth through the core   τ _t∼1  . Since the observational evidence for steep optical-depth profiles in the outer envelopes of some pre-stellar cores (specifically   τ _t∝ b ^{- β}   , with   β ≳2)  constrains only the dust column density, this leaves open the possibility that the gas has a shallower column-density profile.     

Tidal friction in triple stars

Tidal friction in close binaries, with periods of a few days, is expected to circularize the orbit on a time-scale long compared with human observation but shorter than, or comparable to, the lifetimes of main-sequence stars. In a hierarchical triple star, however, the perturbing effect of the distant third star may decircularize the inner orbit significantly on a time-scale of the order of days (as in λ Tau) or centuries (as in β Per). If the inner pair is observed to be semidetached, however, it is plausible to assume that the eccentricity is small. This may be because tidal friction is operating on a comparably short time-scale, and so it is in principle amenable to observation. We attempt to determine a lower limit to the strength of tidal friction in λ Tau and β Per, on the basis of this consideration. Tidal friction will also lead to a secular transfer of angular momentum from the inner orbit to the outer orbit. Too rapid a transfer may lead to orbital shrinkage that is fast compared with the nuclear time-scales of the inner systems, and this can also be ruled out on observational grounds. Thus we may be able to set an upper as well as a lower limit to the strength of tidal friction, on the basis of observations. In a young hierarchical triple, provided that the orbits are fairly nearly orthogonal, tidal friction can serve to reduce the inner orbital period from months to days within a fairly short period of time, of order P ²_out/ P _in. This may be a significant mechanism for producing young short-period binaries.     

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.     

The multiperiodic δ Scuti star 4 Canum Venaticorum: amplitude variability

Recent multisite campaigns of the Delta Scuti Network have revealed 34 frequencies of pulsation for the star 4 CVn. Our present knowledge of the frequencies makes it possible to reanalyse the shorter data sets in the literature, photometric observations from 1966 to 1997.
4 CVn shows strong amplitude variability with time-scales of ten years or longer, although for neighbouring years the amplitudes usually are similar. Seven of the eight dominant modes show annual variability of ∼12 per cent. The variability increases to ∼40 per cent over a decade. The formally derived time-scale of variation of 30 years can only be a rough estimate, since this is also the length of the available data span. The variability is compared with that of FG Vir, which shows lower amplitude variability.
The cyclic behaviour of the amplitude variations excludes an evolutionary origin. There exists some evidence that a mode at 6.12 d⁻¹, which appeared during 1996 and 1997, may have been present with small amplitudes in the 1976–1978 time period.
The pulsation mode at 7.375 d⁻¹ exhibited the most rapid decrease found so far: the V amplitude dropped from the highest known value of 15 mmag in 1974 to 4 mmag in 1976 and 1 mmag in 1977. After that the mode has been increasing in amplitude. There exists a phase jump between 1976 and 1977, suggesting the growth of a new mode. It is interesting to note that this mode also has the strongest coupling with other modes with combination frequencies, f _i ± f _j . The amplitudes of these combination frequencies are also strongly variable from year to year. We speculate that power is transferred between the modes through mode-coupling.     

Correlated spectral variability in brown dwarfs

Models of brown dwarf atmospheres suggest they exhibit complex physical behaviour. Observations have shown that they are indeed dynamic, displaying small photometric variations over time-scales of hours. Here, I report results of infrared (0.95–1.64 μm) spectrophotometric monitoring of four field L and T dwarfs spanning time-scales of 0.1–5.5 h, the goal being to learn more about the physical nature of this variability. Spectra are analysed differentially with respect to a simultaneously observed reference source in order to remove Earth-atmospheric variations. The variability amplitude detected is typically 2–10 per cent, depending on the source and wavelength. I analyse the data for correlated variations between spectral indices. This approach is more robust than single band or  χ²  analyses, because it does not assume an amplitude for the (often uncertain) noise level (although the significance test still assumes a shape for the noise power spectrum). Three of the four targets show significant evidence for correlated variability. Some of this can be associated with specific features including Fe, FeH, VO and K  i , and there is good evidence for intrinsic variability in H₂O and possibly also CH₄. Yet some of this variability covers a broader spectral range which would be consistent with dust opacity variations. The underlying common cause is plausibly localized temperature or composition fluctuations caused by convection. Looking at the high signal-to-noise ratio stacked spectra, we see many previously identified spectral features of L and T dwarfs, such as K  i , Na  i , FeH, H₂O and CH₄. In particular, we may have detected methane absorption at 1.3–1.4 μm in the L5 dwarf SDSS 0539−0059.