GeV emission from neutron-rich internal shocks of some long γ-ray bursts

In the neutron-rich internal shocks model for γ-ray bursts (GRBs), the Lorentz factors (LFs) of ion shells are variable, and so are the LFs of accompanying neutron shells. For slow neutron shells with a typical LF of approximate tens, the typical β-decay radius is  ∼10¹⁴–10¹⁵ cm  . As GRBs last long enough  [ T ₉₀ > 14(1 + z ) s]  , one earlier but slower ejected neutron shell will be swept successively by later ejected ion shells in the range  ∼10¹³–10¹⁵ cm  , where slow neutrons have decayed significantly. Part of the thermal energy released in the interaction will be given to the electrons. These accelerated electrons will mainly be cooled by the prompt soft γ-rays and give rise to GeV emission. This kind of GeV emission is particularly important for some very long GRBs and is detectable for the upcoming satellite Gamma-Ray Large Area Space Telescope (GLAST).     

Jet and accretion power in the most powerful Fermi blazars

Among the blazars detected by the Fermi satellite, we have selected the 23 blazars that in the 3 months of survey had an average γ-ray luminosity above 10⁴⁸ erg s⁻¹. For 17 out of the 23 sources we found and analysed X-ray and optical–ultraviolet data taken by the Swift satellite. With these data, implemented by archival and not simultaneous data, we construct the spectral energy distributions, and interpreted them with a simple one-zone, leptonic, synchrotron and inverse Compton model. When possible, we also compare different high-energy states of single sources, like 0528+134 and 3C 454.3, for which multiple good sets of multiwavelength data are available. In our powerful blazars the high energy emission always dominates the electromagnetic output, and the relatively low level of the synchrotron radiation often does not hide the accretion disc emission. We can then constrain the black hole mass and the disc luminosity. Both are large (i.e. masses equal or greater than  10⁹ M _⊙  and disc luminosities above 10 per cent of Eddington). By modelling the non-thermal continuum we derive the power that the jet carries in the form of bulk motion of particles and fields. On average, the jet power is found to be slightly larger than the disc luminosity, and proportional to the mass accretion rate.     

On the electric field screening by electron–positron pairs in a pulsar magnetosphere

We present a steady one-dimensional model for a pulsar polar cap accelerator, where the field-aligned electric field and flow are solved self-consistently with a given current density. It is assumed that no particles return to the star. It is known that the space-charge-limited flow is accelerated to energies high enough to create electron–positron pairs if the assumed current density is high enough. We find that when pairs are created in such a space-charge-limited flow, the accelerating electric field is screened out within a short distance after pair creation, if the pair particle flux is larger than a critical value. We also find that a space charge density wave is excited in the screening region.
We find that a pair flux larger than the critical value M _c=10³–10⁵ must be reached in a layer with thickness equal to the braking distance for the decelerating component. Therefore, the required multiplicity – the number of pairs created by one primary particle – is too large to be realized in the actual pulsar magnetosphere. We suggest that in order to obtain a localized potential drop along the polar cap magnetic flux, one needs to take into account additional effects such as wave–particle interaction or quasi-periodic pair creation.     

Hubble Space Telescope ultraviolet spectroscopy of blazars: emission-line properties and black hole masses

The ultraviolet (UV) spectra of 16 blazars  (〈 z 〉≃ 1)  from the archives of the Hubble Space Telescope Faint Object Spectrograph have been analysed in order to study in a systematic way the properties of their broad UV emission lines. We find that the luminosities of the most prominent and intense lines, Lyα and C  iv λ1549, are similar to those of normal radio-loud quasars at comparable redshifts. However, the equivalent widths of blazar lines are significantly smaller than those of radio-loud quasars. Therefore, while the intrinsic broad-line region luminosity of blazars appears to be indistinguishable from that of radio-loud quasars, their continuum must be comparatively higher, most probably due to relativistic beaming. We have combined the UV luminosities of the debeamed continuum with the emitting gas velocity to derive estimates of the masses of the central supermassive black holes. The size of the broad-line region was computed in two ways: (1) via an empirical relationship between UV continuum luminosity and broad-line region size, and (2) through the external photon density required by blazar models to reproduce the inverse Compton components observed at γ-rays. The second method yields significantly different results from the first method, suggesting that it provides only a very rough estimate or a lower limit on the size of the broad-line region. We find that the average mass of the central black holes in blazars is  ∼2.8 × 10⁸ M_⊙  , with a large dispersion, comparable to those computed for other radio-loud active galactic nuclei.     

The Deep X-Ray Radio Blazar Survey (DXRBS) – II. New identifications

We have searched the archived, pointed ROSAT Position Sensitive Proportional Counter data for blazars by correlating the WGACAT X-ray data base with several publicly available radio catalogues, restricting our candidate list to serendipitous X-ray sources with a flat radio spectrum ( α _r≤0.70, where S _ν ∝ ν ^{− α} ). This makes up the Deep X-ray Radio Blazar Survey (DXRBS). Here we present new identifications and spectra for 106 sources, including 86 radio-loud quasars, 11 BL Lacertae objects, and nine narrow-line radio galaxies. Together with our previously published objects and already-known sources, our sample now contains 298 identified objects: 234 radio-loud quasars [181 flat-spectrum quasars: FSRQ ( α _r≤0.50) and 53 steep-spectrum quasars: SSRQ], 36 BL Lacs and 28 narrow-line radio galaxies. Redshift information is available for 96 per cent of these. Thus our selection technique is ∼90 per cent efficient at finding radio-loud quasars and BL Lacs. Reaching 5-GHz radio fluxes ∼50 mJy and 0.1–2.0 keV X-ray fluxes a few ×10⁻¹⁴ erg cm⁻² s⁻¹, DXRBS is the faintest and largest flat-spectrum radio sample with nearly complete (∼85 per cent) identification. We review the properties of the DXRBS blazar sample, including redshift distribution and coverage of the X-ray-radio–power plane for quasars and BL Lacs. Additionally, we touch upon the expanded multiwavelength view of blazars provided by DXRBS. By sampling for the first time the faint end of the radio and X-ray luminosity functions, this sample will allow us to investigate the blazar phenomenon and the validity of unified schemes down to relatively low powers.     

Coulomb corrections to the equation of state of nuclear statistical equilibrium matter: implications for SNIa nucleosynthesis and the accretion-induced collapse of white dwarfs

Coulomb corrections to the equation of state of degenerate matter are usually neglected in high-temperature regimes, owing to the inverse dependence of the plasma coupling constant, Γ, on temperature. However, nuclear statistical equilibrium matter is characterized by a large abundance by mass of large- Z (iron group) nuclei. It is found that Coulomb corrections to the ion ideal gas equation of state of matter in nuclear statistical equilibrium are important at temperatures T ≲5–10×10⁹ K and densities ρ ≳10⁸ g cm⁻³. At a temperature T =8.5×10⁹ K and a density ρ =8×10⁹ g cm⁻³, the neutronization rate is larger by ≳28 per cent when Coulomb corrections are included. However, the conductive velocity of a thermonuclear deflagration wave in C–O drops by ∼16 per cent when Coulomb corrections to the heat capacity are taken into account. The implications for SNIa models and nucleosynthesis, and also for the accretion-induced collapse of white dwarfs, are discussed. Particularly relevant is the result that the minimum density for collapse of a white dwarf to a neutron star is shifted down to 5.5–6×10⁹ g cm⁻³, a value substantially lower than previously thought.     

Scattered emission from a relativistic outflow and its application to gamma-ray bursts

We investigate a scenario of photon scattering by electrons within a relativistic outflow. The outflow is composed of discrete shells with different speeds. One shell emits radiation for a short duration. Some of this radiation is scattered by the shell(s) behind. We calculate in a simple two-shell model the observed scattered flux density as a function of the observed primary flux density, the normalized arrival time delay between the two emission components, the Lorentz factor ratio of the two shells and the scattering shell's optical depth. Thomson scattering in a cold shell and inverse Compton scattering in a hot shell are both considered. The results of our calculations are applied to the gamma-ray bursts and the afterglows. We find that the scattered flux from a cold slower shell is small and likely to be detected only for those bursts with very weak afterglows. A hot scattering shell could give rise to a scattered emission as bright as the X-ray shallow decay component detected in many bursts, on a condition that the isotropically equivalent total energy carried by the hot electrons is large, ∼10⁵²–10⁵⁶ erg. The scattered emission from a faster shell could appear as a late short γ-ray/MeV flash or become part of the prompt emission depending on the delay of the ejection of the shell.     

Monopole gravitational waves from relativistic fireballs driving gamma-ray bursts

Einstein's general relativity predicts that pressure, in general stresses, plays a similar role to energy density,  ε=ρ c ²  (with ρ being the corresponding mass density), in generating gravity. The source of gravitational field, the active gravitational mass density, sometimes referred to as Whittaker's mass density, is  ρ_grav=ρ+ 3 p / c ²  , where p is pressure in the case of an ideal fluid. Whittaker's mass is not conserved, hence its changes can propagate as monopole gravitational waves. Such waves can be generated only by astrophysical sources with varying gravitational mass. Here we show that relativistic fireballs, considered in modelling gamma-ray burst phenomena, are likely to radiate monopole gravitational waves from high-pressure plasma with varying Whittaker's mass. Also, ejection of a significant amount of initial mass-energy of the progenitor contributes to the monopole gravitational radiation. We identify monopole waves with   h ¹¹+ h ²²  waves of Eddington's classification which propagate (in the z -direction) together with the energy carried by massless fields. We show that the monopole waves satisfy Einstein's equations, with a common stress-energy tensor for massless fields. The polarization mode of monopole waves is  Φ₂₂  , i.e. these are perpendicular waves which induce changes of the radius of a circle of test particles only (breathing mode). The astrophysical importance of monopole gravitational waves is discussed.     

Calculated spectra for HeH+ and its effect on the opacity of cool metal-poor stars

The wavelength and Einstein A coefficient are calculated for all rotation–vibration transitions of  ⁴He¹H⁺, ³He¹H⁺, ⁴He²H⁺  and  ³ He²H⁺  , giving a complete line list and the partition function for  ⁴HeH⁺  and its isotopologues. This opacity is included in the calculation of the total opacity of low-metallicity stars and its effect is analysed for different conditions of temperature, density and hydrogen number fraction. For a low helium number fraction (as in the Sun), it is found that HeH⁺ has a visible but small effect for very low densities  (ρ≤ 10⁻¹⁰ g cm⁻³)  , at temperatures around 3500 K. However, for high helium number fraction, the effect of HeH⁺ becomes important for higher densities  (ρ≤ 10⁻⁶ g cm⁻³)  , its effect being most important for a temperature around 3500 K. Synthetic spectra for a variety of different conditions are presented.     

The neutrino emission due to plasmon decay and neutrino luminosity of white dwarfs

One of the effective mechanisms of neutrino energy losses in red giants, pre-supernovae and in the cores of white dwarfs is the emission of neutrino–antineutrino pairs in the process of plasmon decay. In this paper, we numerically calculate the emissivity due to plasmon decay in a wide range of temperatures 10⁷–10¹¹ K and densities (2 × 10²–10¹⁴) g cm⁻³. Numerical results are approximated by convenient analytical expressions. We also calculate and approximate by analytical expressions the neutrino luminosity of white dwarfs due to plasmon decay, as a function of their mass and internal temperature. This neutrino luminosity depends on the chemical composition of white dwarfs only through the parameter μ_e (the net number of baryons per electron) and is the dominant neutrino luminosity in all white dwarfs at the neutrino cooling stage.     

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.     

Wisps and knots in the central Crab nebula

The fast-spinning Crab pulsar (∼30 turn s⁻¹), which powers the massive expansion and synchrotron emission of the entire Crab nebula, is surrounded by quasi-stationary features such as fibrous arc-like wisps and bright polar knots in the radial range of 2×10¹⁶≲ r ≲2×10¹⁷ cm, as revealed by high-resolution (∼0.1 arcsec) images from the Wide Field and Planetary Camera 2 (WFPC2) on board the Hubble Space Telescope ( HST ). The spin-down energy flux (∼5×10³⁸ erg s⁻¹) from the pulsar to the luminous outer nebula, which occupies the radial range 0.1≲ r ≲2 pc, is generally believed to be transported by a magnetized relativistic outflow of an electron–positron e^± pair plasma. It is then puzzling that mysterious structures like wisps and knots, although intrinsically dynamic in synchrotron emission, remain quasi-stationary on time-scales of a few days to a week in the relativistic pulsar wind. Here we demonstrate that, as a result of slightly inhomogeneous wind streams emanating from the rotating pulsar, fast magnetohydrodynamic (MHD) shock waves are expected to appear in the pulsar wind at relevant radial distances in the forms of wisps and knots. While forward fast MHD shocks move outward with a speed close to the speed of light c , reverse fast MHD shocks may appear quasi-stationary in space under appropriate conditions. In addition, Alfvénic fluctuations in the shocked magnetized pulsar wind can effectively scatter synchrotron beams from gyrating relativistic electrons and positrons.     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the gamma-ray burst (GRB) prompt emission site from the explosion centre R , by determining the location of the electron's self-absorption frequency in the GRB prompt optical-to-X/γ-ray spectral energy distribution, assuming that the optical and the γ-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely the bulk Lorentz factor of the emitting source Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 300  . For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R -constraints for this non-detection sample and those for the detection sample.     

Can a slowly rotating neutron star be a radio pulsar?

It is shown that the radius of curvature of magnetic field lines in the polar region of a rotating magnetized neutron star can be significantly less than the usual radius of curvature of the dipole magnetic field. The magnetic field in the polar cap is distorted by toroidal electric currents flowing in the neutron star crust. These currents close up the magnetospheric currents driven by the electron–positron plasma generation process in the pulsar magnetosphere. Owing to the decrease in the radius of curvature, electron–positron plasma generation becomes possible even for slowly rotating neutron stars, with   PB ^−2/3₁₂ < 10 s  , where P is the period of star rotation and   B ₁₂= B /10¹² G  is the magnitude of the magnetic field on the star surface.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .     

Atomic hydrogen shells in the Messier 82 starburst

We report on a search for atomic hydrogen holes and shells in the nearby starburst galaxy M82, using high angular resolution (∼1.3 arcsec) VLA H  i absorption observations. From this study, we have detected four H  i shells in the central kiloparsec of M82. The sizes of these shells (∼30–50 pc) are smaller than those of the majority of shells observed in the Large Magellanic Cloud, although the M82 shells have higher expansion velocities (∼30 km s⁻¹) and typical kinetic energies of  10⁵¹–10⁵² erg  . Because our observations were made in absorption, strong selection effects are present which hinder the detection of shells that could be present outside, or behind, the extended radio continuum associated with the starburst. Nevertheless, our detection of four shells in M82 actually represents a higher density of shells per unit area compared with the Large Magellanic Cloud.
We also discuss the gas dynamics in the central kiloparsec of M82, and discuss the velocity structure of gas in a barred potential and in wind-driven shells. We conclude that in M82 the observed gas dynamics are most likely a superposition of both effects.     

Fe xiii emission lines in active region spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent fully relativistic calculations of radiative rates and electron impact excitation cross-sections for Fe  xiii are used to generate emission-line ratios involving 3s²3p²–3s3p³ and 3s²3p²–3s²3p3d transitions in the 170–225 and 235–450 Å wavelength ranges covered by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS). A comparison of these line ratios with SERTS active region observations from rocket flights in 1989 and 1995 reveals generally very good agreement between theory and experiment. Several new Fe  xiii emission features are identified, at wavelengths of 203.79, 259.94, 288.56 and 290.81 Å. However, major discrepancies between theory and observation remain for several Fe  xiii transitions, as previously found by Landi and others, which cannot be explained by blending. Errors in the adopted atomic data appear to be the most likely explanation, in particular for transitions which have 3s²3p3d ¹D₂ as their upper level. The most useful Fe  xiii electron-density diagnostics in the SERTS spectral regions are assessed, in terms of the line pairs involved being (i) apparently free of atomic physics problems and blends, (ii) close in wavelength to reduce the effects of possible errors in the instrumental intensity calibration, and (iii) very sensitive to changes in N _e over the range  10⁸–10¹¹ cm⁻³  . It is concluded that the ratios which best satisfy these conditions are 200.03/202.04 and 203.17/202.04 for the 170–225 Å wavelength region, and 348.18/320.80, 348.18/368.16, 359.64/348.18 and 359.83/368.16 for 235–450 Å.     

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.     

The role of non-thermal electrons in the optical continuum of stellar flares

The continuum emission of stellar flares in UV and visible bands can be enhanced by two or even three orders of magnitude relative to the quiescent level and is usually characterized by a blue colour. It is difficult for thermal atmospheric models to reproduce all these spectral features. If the flaring process involves the acceleration of energetic electrons which then precipitate downwards to heat the lower atmosphere, collisional excitation and ionization of ambient hydrogen atoms by these non-thermal electrons could be important in powering the continuum emission. To explore such a possibility, we compute the continuum spectra from an atmospheric model for a dMe star, AD Leo, at its quiescent state, when considering the non-thermal effects by precipitating electron beams. The results show that if the electron beam has an energy flux large enough (for example, ℱ₁∼10¹² erg cm⁻² s⁻¹), the U -band brightening and, in particular, the U − B colour are roughly comparable with observed values for a typical large flare. Moreover, for electron beams with a moderate energy flux ℱ₁≲10¹¹ erg cm⁻² s⁻¹, a decrease of the emission at the Paschen continuum appears. This can explain at least partly the continuum dimming observed in some stellar flares. Adopting an atmospheric model for the flaring state can further raise the continuum flux, but it yields a spectral colour incomparable with observations. This implies that the non-thermal effects may play the chief role in powering the continuum emission in some stellar flares.     

Faraday rotation in the VLBI core of BL Lacertae

External Faraday rotation has been detected in both the core and the parsec-scale jet of BL Lac in a four-frequency very long baseline interferometry (VLBI) experiment. This unexpected result indicates the presence of significant amounts of thermal gas close to the nucleus of this object. The rotation measure (RM) in the jet components is constant, and differs from the currently accepted Galactic RM, indicating that this value (−205 rad m⁻²) is not applicable to the components in the parsec-scale jet. The similarity of the RM in these jet components leads us to suspect that the jet RM is caused by a foreground screen in our Galaxy, although we cannot rule out a combination of Galactic RM and RM local to the jet. If the jet RM is due solely to the Galaxy, this would mean that the currently accepted value of the foreground RM (−205 rad m⁻² ) is not correct, either because the value changed between 1982 and 1997, or because the assumption of no intrinsic source rotation was incorrect, as it was at our later epoch of observation. Our observations suggest a value of     .
After correcting for the foreground RM, the core value is −427 rad m⁻², which is unexpected since, owing to the weakness of their line emission, BL Lac objects are often assumed to be depleted in gas. The core RM appears to be variable, probably because of the presence of at least two polarized components close to the core the relative contributions of which vary with time.