Density profiles of dark matter halos with anisotropic velocity tensors

We present density profiles, that are solutions of the spherical Jeans equation, derived under the following two assumptions: (i) the coarse grained phase-density follows a power-law of radius, ρ/σ³ ∝ r⁻, and (ii) the velocity anisotropy parameter is given by the relation β_a(r)=β₁+2β₂ (r/r_*)/[1+(r/r_*)²] where β₁, β₂ are parameters and r_* equals twice the virial radius, r_vir, of the system. These assumptions are well motivated by the results of N-body simulations. Density profiles have increasing logarithmic slopes γ, defined by γ=−d ln ρ/d ln r. The values of γ at r=10^−2.5r_vir, a distance where the systems could be resolved by large N-body simulations, lie in the range 1.0–1.6. These inner values of γ increase for increasing β₁ and for increasing concentration of the system. On the other hand, slopes at r=r_vir lie in the range 2.42–3.82. A model density profile that fits well the results at radial distances between 10⁻³r_vir and r_vir and connects kinematic and structural characteristics of spherical systems is described.     

Photoelectric monitoring of a shell star in perseus

We report on the 3-colour photoelectric monitoring of the shell star in Perseus, discovered in 1978, [1]. On the night of 1983 Feb 4–5, this star still showed obvious light variations, though the amplitude was less than in 1979. On the next night, the light variation was close to the observing error. We surmise that the shell activity is now at a stable phase for a few years.     

A theoretical study of the direction and polarization of solar radio millisecond spike radiation

We investigated the angular direction and polarization of the solar radio millisecond spike emission in the model in which the spike emission is due to the second harmonic instability modes driven by electron cyclotron maser of loss cone distributed electrons during the propagation of a nonlinear plasma density wave near the magnetic mirror. We found that, when the angle θ between the wave vector and the magnetic field is > 60 °, the emission is in 100% X-mode polarization; when 40 ° < θ<60 °, the emission is in 100% O-mode polarization provided the amplitude of the density wave is below a certain limit; above that limit, the polarization will fall from 100% O-mode to even the X-mode. We also found that only 0.1% of the free energy of energy carrying electrons in the source region is converted into radiation wave energy.     

The cosmological constant and the clusters of galaxies

The dependence of the formation redshift of clusters of galaxies on the cosmological constant is discussed on the basis of recent measurements of gravitational masses and radii from X-ray astronomy techniques. For a flat CDM model with several different values of the cosmological constant, the cluster number density as a function of the redshift is calculated. Comparing the calculated number density at z − 0 with the observed value limits the spectrum parameter Γ to the range 0.15–0.32. We also show that if Γ and/or Ω_Λ were larger, many more clusters of galaxies would have been formed at high redshifts.     

LS I +61 303 as a potential neutrino source on the light of magic results

Very high energy γ-rays have recently been detected from the microquasar LS I +61 303 using the MAGIC telescope. A phenomenological study on the concomitant neutrinos that would be radiated if the γ-ray emission is hadronic in origin is herein presented. Neutrino oscillations are considered, and the expected number of events in a km-scale detector such as ICECUBE is computed under different assumptions including orbital periodicity and modulation, as well as different precision in the modeling of the detector. We argue that the upper limits already imposed on the neutrino emission of LS I +61 303 using AMANDA-II and the forthcoming measurements by ICECUBE may significantly constrain – in an independent and unbiased way – the γ-ray to neutrino flux ratio, and thus the possibility of a hadronic origin of the γ-rays. The viability of hadronic models based on wind–jet interactions in the LS +61 303 system after MAGIC measurements is discussed.     

Gamma-ray and synchrotron emission from neutralino annihilation in the Large Magellanic Cloud

We calculate the expected flux of γ-ray and radio emission from the LMC due to neutralino annihilation. Using rotation curve data to probe the density profile and assuming a minimum disk, we describe the dark matter halo of the LMC using models predicted by N-body simulations. We consider a range of density profiles including the NFW profile, a modified NFW profile proposed by Hayashi et al. (2003) to account for the effects of tidal stripping, and an isothermal sphere with a core. We find that the γ-ray flux expected from these models may be detectable by GLAST for a significant part of the neutralino parameter space. The prospects for existing and upcoming Atmospheric Cherenkov Telescopes (ACTs) are less optimistic, as unrealistically long exposures are required for detection. However, the effects of adiabatic compression due to the baryonic component may improve the chances for detection by ACTs. The maximum flux we predict is well below EGRET's measurements and thus EGRET does not constrain the parameter space. The expected synchrotron emission generally lies below the observed radio emission from the LMC in the frequency range of 19.7–8550 MHz. As long as σv<2×10⁻²⁶ cm³ s⁻¹ for a neutralino mass of 50 GeV, the observed radio emission is not primarily due to neutralinos and is consistent with the assumption that the main source is cosmic rays. We find that the predicted fluxes, obtained by integrating over the entire LMC, are not very strongly dependent on the inner slope of the halo profile, varying by less than an order of magnitude for the range of profiles we considered.     

Spectrum of the relic neutrino background from past supernovae and cosmological models

It is greatly expected that the relic neutrino background from past supernovae will be detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, Λ, on the spectrum, since some recent cosmological observations strongly suggest the existence of Λ. We find following results. (1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimates (6–10 MeV). (2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae are dominant, is about 25h₅₀²(R_SN/0.1 yr⁻¹)(n_Gh₅₀⁻³/0.02 Mpc⁻³) events per year, where R_SN is the supernova rate in a galaxy, n_G is the number density of galaxies, and h₅₀ = H₀/(50 km/s Mpc), where H₀ is the Hubble constant. (3) The event rate is almost insensitive to Λ. The flux increases in the low energy side (< 10 MeV) with increasing Λ, but decreases in the high energy side (> 10 MeV) in models in which the integrated number of supernovae in one galaxy is fixed.     

Probing pulsar winds using inverse Compton scattering

We investigate the effects of inverse Compton scattering by electrons and positrons in the unshocked winds of rotationally-powered binary pulsars. This process can scatter low energy target photons to produce gamma rays with energies from MeV to TeV. The binary radio pulsars PSR B1259−63 and PSR J0045−73 are both in close eccentric orbits around bright main sequence stars which provide a huge density of low energy target photons. The inverse Compton scattering process transfers momentum from the pulsar wind to the scattered photons, and therefore provides a drag which tends to decelerate the pulsar wind. We present detailed calculations of the dynamics of a pulsar wind which is undergoing inverse Compton scattering, showing that the deceleration of the wind of PSR B1259−63 due to ‘inverse Compton drag' is small, but that this process may confine the wind of PSR J0045−73 before it attains pressure balance with the outflow of its companion star. We calculate the spectra and light curves of the resulting inverse Compton emission from PSR B1259−63 and show that if the size of the pulsar wind nebula is comparable to the binary separation, then the γ-ray emission from the unshocked wind may be detectable by atmospheric Cherenkov detectors or by the new generation of satellite-borne γ-ray detectors such as INTEGRAL and GLAST. This mechanism may therefore provide a direct probe of the freely-expanding regions of pulsar winds, previously thought to be invisible.     

The potential of the ground based arrays of imaging atmospheric Cherenkov telescopes. II. Gamma ray flux sensitivities

We discuss the capability of ‘100 GeV’ class imaging atmospheric Cherenkov telescope (IACT) arrays as future powerful instruments of ground-based gamma-ray astronomy. It is assumed that the array is gathered from individually triggered quadrangular 4-IACT ‘cells’ with a linear size of about 100 m. The multi-cell concept allows coverage of large detection areas economically, and at the same time the effective exploitation of the stereoscopic approach of determination of the shower parameters using information obtained by several IACTs simultaneously. Determination of arrival directions of γ-ray primaries on an event-by-event basis with accuracy δθ ≤ 0.1° combined with high suppression efficiency (at both the hardware and software levels) of the background hadronic showers by a factor of ≈ 10³, and large, up to 1 km² collection areas, can provide minimum detectable energy fluxes of ≥ 100 GeV γ-rays from point sources down to 10⁻¹³ erg/cm² s which is about 3 orders of magnitude lower than the current sensitivities achieved by the satellite-borne detectors at MeV and GeV energies. High sensitivities of multi-IACT arrays would partially compensate the limited efficiency of the technique for all-sky surveys, as well as allow study of moderately extended (≤ 1°) γ-ray sources. IACT arrays with minimum detectable fluence of ≥ 100 GeV γ-rays S_γ < 10⁻⁸ erg/cm² are well suited for effective exploration of highly sporadic nonthermal phenomena from different classes of astrophysical objects on time-scales from ≤ 1 s to several minutes.     

The pionic contribution to diffuse γ-rays: upper limits

Diffuse γ-rays probe the highest-energy processes at the largest scales. Here we derive model-independent constraints on the hadronic contribution to the Galactic and extragalactic γ-ray spectra at in the energy range 50 MeVE_γ10 GeV. The hadronic component is dominated by emission from neutral pions, with a characteristic spectrum symmetric about m_π⁰/2. We exploit the well-defined properties of the pion decay spectrum to quantify the maximum pionic fraction of the observed γ-ray intensity. We find that the Galactic spectrum above 30 MeV can be atmost about 50% pionic. The maximum pionic contribution to the extragalactic spectrum is energy dependent; it also depends on the redshift range over which the sources are distributed, ranging from as low as about 20% for pions generated very recently, to as much as 90% if the pions are generated around redshift 10. The implications of these constraints for models of γ-ray and neutrino emission are briefly discussed.     

Constraining the power spectrum using clusters

We analyze an extended redshift sample of Abell/ACO clusters and compare the results with those coming from numerical simulations of the cluster distribution, based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, so that we generate a set of 48 independent mock Abell/ACO cluster samples per model, on which we estimate cosmic variance effects. Other than the standard CDM model, we consider (a) Ω₀ = 1 CDM models based on lowering the Hubble parameter and/or on tilting the primordial spectrum; (b) Ω₀ = 1 Cold + Hot DM models with 0.1 ≤Ω_ν ≤0.5; (c) low-density flat ΛCDM models with 0.3 ≤Ω₀ ≤0.5. We compare real and simulated cluster distributions by analysing correlation statistics, the probability density function, and supercluster properties from percolation analysis. We introduce a generalized definition of the spectrum shape parameter Γ in terms of σ₂₅/σ₈, where σ_ris the rms fluctuation amplitude within a sphere of radius r. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 Γ 0.25 and an effective spectral index at the 20 h⁻¹ Mpc scale −1.1 n_eff −0.9 are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press-Schechter approach, whose reliability is explicitly tested against N-body simulations. By combining results from the analysis of the distribution and the abundance of clusters we conclude that, of the cosmological models considered here, the only viable models are either Cold + Hot DM ones with 0.2 Ω_ν 0.3, better if shared between two massive ν species, and ΛCDM ones with 0.3 Ω₀0.5.     

Separation of gamma and hadron initiated air showers with energies between 20 and 500 TeV

The discrimination between air showers initiated by γ rays and by hadrons is one of the fundamental problems in experimental cosmic-ray physics. The physics of this ‘γ/hadron separation’ is discussed in this paper. We restrict ourselves to the energy range from about 20 to 500 TeV, and take only the information contained in the lateral Čerenkov light distribution and the number of electrons at the detector level into consideration. An understanding of the differences between air showers generated by γ rays and those due to hadrons leads us to formulate suitable observables for the separation process. Angle integrating Čerenkov arrays (AICA) offer a promising new approach to ground-based γ-ray astronomy in the energy region from about 20 to 500 TeV. In order to establish this technique, an efficient suppression of the overwhelming hadronic background radiation is required. As an example for our general discussion, we present one method for γ/hadron separation in AICAs called ‘LES’. It is based on the simultaneous determination of the shower size and some characteristic parameters of the lateral distribution of the Čerenkov light. The potential inherent within this technique is demonstrated in quantitative detail for the existing ‘AIROBICC’ AICA. We also propose an objective measure of the intrinsic sensitivity of a detection scheme in ground-based γ-ray astronomy, the ‘reduced quality factor’. It is shown that AICAs may reach a sensitivity to γ-ray point sources in the high VHE range similar to that of the Čerenkov-telescope imaging technique in the low VHE region.     

The self-gravity of Uranian rings

We derive an approximate, analytic expression for the self-gravity of Uranian elliptical rings at any point in the N-streamline formalism (N → ∞), and numerically evaluate by a Monte Carlo calculation, the correction factor γ that corrects for the approximations used in the derivation. We confirm that the 2-streamline (N = 2) model overestimates the mass of the ring: for the ε, , ßrings, the mass given by the N-streamline model (N → ∞) is only 2/5 that by the 2-streamline model.     

Report on the physical characteristics of Vaidya-Tikekar's exact relativistic model for a superdense star

We report here the results of our examination of the physical properties of Vaidya-Tikekar's model for a relativistic star. Full details will be published elsewhere. The analysis yields a strong indication that the model is stable with respect to infinitesimal radial oscillations. We find that the adiabatic speed of sound is smaller than the speed of light everywhere inside the fluid sphere if the radius of the sphere is larger than 1.46 times its Schwarzschild radius. We also find that the fluid must necessarily be supraluminal somewhere if the speed of sound is decreasing outwards close to the center. We further find that the strong energy condition is fulfilled everywhere if it is fulfilled at the center. Since the ratio of the pressure p and the density ⋅ is decreasing outwards, this indicates that the temperature gradient is negative. We also find that the relativistic adiabatic index is larger than two. Demanding the fluid to be causal, and taking the pressure and the density to be somewhere given by 7.4 ⋅ 10³³ dynes/cm³ and 5.1 ⋅ 10¹⁴ g/cm³, we calculate the maximum mass of the fluid sphere to be 3 solar masses.     

Gamma–hadron separation methods for the VERITAS array of four imaging atmospheric Cherenkov telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive γ-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 10⁵ m², combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic-ray initiated air showers that resemble γ-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of σ_63% = 0.1° (0.2°) and σ_68% = 15% (22%), respectively. The main emphasis of the paper is the discussion of γ–hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a γ–hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.     

General relativistic effects in the γ-ray radiation from pulsars

General relativistic and quantum-mechanical effects in the γ-ray radiation from neutron stars are studied in this paper. For three values of the magnetic dipole moment, we calculated the upper limit of the escaping energy of the γ-photon as a function of the angle between the line of sight and the magnetic axis, with and without these effects. We conclude that for low magnetic moments, the quantum-mechanical effects are unimportant but the general relativistic effects are to reduce significantly the upper limit. For high magnetic moments, the two effects nearly cancel out and the net effect is small. We discuss the implications of these results on the γ-ray spectrum.     

Strong correlation between the peak power density in the time domain and the luminosity of long γ-ray bursts

This work presents a possible luminosity estimator for long γ-ray bursts (GRBs) based on their light curves. We use the method of variability analysis in the time domain to calculate the power density spectrum (PDS) for each of the 12 GRBs with known redshifts observed by CGRO/BATSE. The peak of the power density spectrum P is a measure of the intensity of variability of the given light curve and a strong correlation is found between P and the isotropic peak luminosity L of the GRB. It is a successor to the lag-luminosity relation of Norris et al. (2000) and the variability-luminosity relation of Reichart et al. (2001).     

Full-wave VLF modes in a cylindrically symmetric enhancement of plasma density

VLF waves in the magnetosphere may be guided by ducts, i.e. enhancements of plasma density aligned with the geomagnetic field. It is required to account for the good transmission properties of these ducts as some leakage of energy is expected, with consequent attenuation. We consider a cylindrical duct, whose axis is parallel to a uniform magnetic field. The electron density of a cold plasma is a function only of radial distance r from the axis, taking constant values in inner and outer regions r < a and r > b, and varying smoothly in the duct wall a < r < b. We compute full-wave ‘trapped’ VLF modes. Results are presented for a range of values of the parameters of the model. In general we find that attenuation of leading modes is very low, agreeing with observation. Specific features are that leakage (arising from mode conversion in the duct wall) increases as the wall is made thinner, that leakage for the leading modes is less for weaker ducts, that many modes may propagate, and that, because of interference effects, some higher order modes also suffer little attenuation.     

Implications of modified AGN spectra due to absorption by infrared photons

Recent observations of high energy photons from active galactic nuclei (AGN) have lead to a renewed interest in the effect of gamma-ray cascading in background radiation fields. In recent years, numerous authors have explored possible modifications to the observed high energy photon spectra from AGN due to propagation through such fields. This paper will re-examine a number of these issues. The major conclusions of this paper are: (1) The reaction γ + γ → e⁺ + e⁻ provides a highly sensitive probe of background infrared (IR) fields. However little can be inferred about the specific nature of the IR background or the effect on spectra from more distant AGN from a study of nearby sources alone. Currently, only upper limits to the background IR density can be established. (2) The contribution of secondary photons from pair-cascading off of background microwave and infrared radiation is most likely unobservable in the regime of 1 TeV unless the strength of extragalactic magnetic fields are much less than 10⁻¹³ G and the inherent source spectra continue to much higher energies. The possible contribution from this process may be conclusively ruled out through further spectral measurements of AGN in the TeV regime during high and low states. (3) There is little hope of unambiguously extracting the values of either the Hubble constant or Ω from γ-γ attenuation measurements. (4) The sensitivity of the attenuation to the density of the background IR produced prior to the epoch of the observed source suggests a future possibility to probe directly the evolution of the background IR radiation through multiple source observations.     

Improvement of gamma-hadron discrimination at TeV energies using a new parameter, image Surface Brightness

The recent advances in TeV γ-ray astronomy using imaging Cerenkov telescopes is mainly due to the good separation between the properties of γ-ray and hadronic showers. A new method for this discrimination is described. The method is mainly based on the change of the image's Size and Surface Brightness (Size/Area) with energy and the differences in these parameters between those initiated by γ-rays and hadrons. The application of the method to the Whipple's 10 m telescope database from the Crab Nebula and Markarian 421 demonstrates its sensitivity to gamma-hadron separation. This study shows the importance of the image Surface Brightness, in addition to the previously used Shape and Orientation parameters.