Improved orbital parameters of accreting millisecond pulsar SAX J1808.4-3658

We analyze the three outbursts of the X-ray millisecond pulsar SAX J1808.4-3658 that occurred in 1998, 2000, and 2002 observed with RXTE. With a technique based on epoch folding search we find an unique orbital solution valid over the five years of high temporal resolution data available. We revise the estimate of the orbital period, P_orb = 7249.1569(1) s and of its error, which we decrease by one order of magnitude. We also give the first constraint on the orbital period derivative, . We find that in 2002 the pulse profile shape is clearly asymmetric, showing a secondary peak at about 145° from the main pulse, which is different from the sinusoidal shape reported at the beginning of the 1998 outburst.     

Force-free magnetic field in a cylindrical flux rope without a constant alpha

It is generally assumed that magnetic fields inside interplanetary magnetic clouds and flux ropes in the solar photosphere are force-free. In order to model such fields, the solution of rot B = B is commonly used where  = const. But comparisons of this solutions with observations show significant difference. To treat this problem,we examine the solutions with .     

From solar nanoflares to stellar giant flares: Scaling laws and non-implications for coronal heating

In this study we explore physical scaling laws applied to solar nanoflares, microflares, and large flares, as well as to stellar giant flares. Solar flare phenomena exhibit a fractal volume scaling, V(L)  L^1.9, with L being the flare loop length scale, which explains the observed correlation between the total emission measure EM_p and flare peak temperature T_p in both solar and stellar flares. However, the detected stellar flares have higher emission measures EM_p than solar flares at the same flare peak temperature T_p, which can be explained by a higher electron density that is caused by shorter heating scale height ratios s_H/L ≈ 0.04–0.1. Using these scaling laws we calculate the total radiated flare energies E_X and thermal flare energies E_T and find that the total counts C are a good proxy for both parameters. Comparing the energies of solar and stellar flares we find that even the smallest observed stellar flares exceed the largest solar flares, and thus their observed frequency distributions are hypothetically affected by an upper cutoff caused by the maximum active region size limit. The powerlaw slopes fitted near the upper cutoff can then not reliably be extrapolated to the microflare regime to evaluate their contribution to coronal heating.     

AL index dependence on the solar wind input revisited

Estimates of the geomagnetic indices made with real-time solar wind measurements form the basis of many space weather forecast techniques. We analyze 20 years of hourly AL and OMNI solar wind data to determine geomagnetic importance of various solar wind and IMF parameters. Besides the solar wind driver of primary importance (VB_s), the first-order contributions, significantly increasing the quality of the model are: solar wind velocity, 2 h of solar wind prehistory and 1 h of AL history. The factors of secondary importance, marginally contributing to overall statistical quality, are IMF B_y, solar wind density, and IMF fluctuations. The dynamic pressure is geomagnetically effective only if the pressure is lower than the average. Modelling of the same data set with an artificial neural network (ANN) confirmed our selection of important factors. Statistically the ANN model was just marginally better than our analytic expression . The AU index dependence is principally different from AL in several respects; therefore modelling of the AE composite index is physically misleading.     

X-ray spectroscopy and variability of AGN detected in the 2 Ms Chandra Deep Field-North Survey

We investigate the nature of the faint X-ray source population through X-ray spectroscopy and variability analyses of 136 active galactic nuclei (AGN) detected in the 2 Ms Chandra Deep Field-North Survey with >200 background-subtracted 0.5–8.0 keV counts [F_{0.5–8.0 keV} = (1.4−200) × 10⁻¹⁵ erg cm⁻² s⁻¹]. Our preliminary spectral analyses yield median spectral parameters of Γ = 1.61 and intrinsic N_H = 6.2 × 10²¹ cm⁻² (z = 1 assumed when no redshift available) when the AGN spectra are fitted with a simple absorbed power-law model. However, considerable spectral complexity is apparent (e.g., reflection, partial covering) and must be taken into account to model the data accurately. Moreover, the choice of spectral model (i.e., free vs. fixed photon index) has a pronounced effect on the derived JVH distribution and, to a lesser extent, the X-ray luminosity distribution. We also find that among the 136 AGN, 10 (≈7%) show significant Fe K emission-line features with equivalent widths in the range 0.1–1.3 keV. Two of these emission-line AGN could potentially be Compton thick (i.e., Γ < 1.0 and large Fe K equivalent width). Finally, we find that 81 (≈60%) of the 136 AGN show signs of variability, and that this fraction increases significantly (≈80–90%) when better photon statistics are available.     

Simultaneous determination of aerosol- and surface characteristics from top-of-atmosphere reflectance using MERIS on board of ENVISAT

The determination of aerosol optical thickness (AOT) from nadir scanning multi-spectral radiometers, like SeaWiFS, MERIS or MODIS, requires the separation of spectral atmospheric and surface properties. Since SeaWiFS and MERIS do not provide information at 2.1 μm, like MODIS, the estimation of the surface reflectance cannot be made by the cross correlation approach described by Kaufman et al., 1997. The BAER approach (Bremen AErosol Retrieval), von Hoyningen-Huene et al., 2003, uses a linear mixing model of spectra for ‘green vegetation’ and ‘bare soil’, tuned by the NDVI, determining an apparent surface to enable this separation of aerosol and surface properties from VIS and NIR channels. Thus AOT can be derived over a wide range of land surfaces for wavelengths <0.67 μm. Using MERIS L1 data over Europe, the AOT retrieved is comparable with ground-based observations, provided by AERONET. Regional variation of AOT can be observed, showing the atmospheric variability for clear sky conditions by: large scale variation of aerosol turbidity, regional pollution, urban regions, effects of contrails and cases of aerosol-cloud interaction. Simultaneously with the spectral AOT also spectral surface reflectance is obtained, where all atmospheric influences have been considered (molecules, aerosols and absorbing gases (O₃)) for channels with wavelengths <0.67 μm. The AOT is extrapolated by Angström power law to NIR channels and the atmospheric correction for land surface properties is performed, enabling the further investigation of land use and spectral land properties.     

XMM-Newton observations of low luminosity Be/X-ray candidates

A small number of early Be stars exhibit X-ray luminosities intermediate between those typical of early type stars and those radiated by Be/X-ray binaries in the quiescent state. We report on XMM-Newton observations of two such Be stars, HD 161103 and SAO 49725 which were originally discovered in a systematic cross-correlation between the ROSAT all-sky survey and SIMBAD. The new observations confirm the X-ray luminosity detected by ROSAT (L_X  10³² erg s⁻¹) and the hardness of their X-ray spectra (thin thermal with kT  8–10 keV or power law with photon index of 1.7) which are both unusual for normal early type stars. We discuss the possible origin of this excess X-ray emission in the light of the models proposed for γ-Cas, magnetic disc-star interaction or accretion onto a compact companion object, neutron star or white dwarf, and compare the properties of these two sources with those of the new massive systems discovered in the XMM- Newton/SSC survey of the Galactic plane.     

Multiband echo tomography of Sco X-1

We present preliminary results of a simultaneous X-ray/optical campaign of the prototypical LMXB Sco X-1 at 1–10 Hz time resolution. Lightcurves of the high excitation Bowen/HeII emission lines and a red continuum at λ_c  6000 Å were obtained through narrow interference filters with ULTRACAM, and these were cross-correlated with simultaneous RXTE X-ray lightcurves. We find evidence for correlated variability, in particular when Sco X-1 enters the Flaring branch. The Bowen/HeII lightcurves lag the X-ray lightcurves with a light travel time which is consistent with reprocessing in the companion star while the continuum lightcurves have shorter delays consistent with reprocessing in the accretion disc.     

Storm-time formation of a relativistic electron belt and some relevant phenomena in other magnetospheric plasma domains

An empirical formula relating the strength of a storm given by its |Dst|_max with the L-coordinate of the peak of storm-injected relativistic electrons is one of a few well-confirmed quantitative relations found in the magneto-spheric physics. We successively extended a dataset of the formula’s basic storms with several events of high Dst-amplitude up to the highest observed |Dst|_max = 600 nT. Possible applying of the formula to the predicting of the ring-current plasma-pressure distribution and the lowest westward electrojet position for a storm are discussed. We have also analyzed the 2000–2001 years’ data on relativistic electrons from our instruments installed on EXPRESS-A (geosynchronous orbit; E_e = 0.8–6 MeV), Molniya-3 (h = 500 × 40 000 km, i = 63°; E_e = 0.8–5.5 MeV) and GLONASS (h = 20 000 km, i = 64°; E_e  l MeV) along with other correlated measurements: GOES series (E_e > 2 MeV), geomagnetic indices (Dst, AE, AL) and interplanetary parameters (solar wind, IMF). The goal is to investigate which outer conditions are most responsible for the high/low output of the storm-injected relativistic electrons. For the geosynchronous orbit, two factors are found as the necessary condition of the highest electron output: high and long-lasting substorm activity on a storm recovery phase and high velocity of solar wind. On the contrary, extremely low substorm activity surely observed during whole the storm recovery phase constitutes a sufficient condition of the non-increased after-storm electron intensity. For the first time found cases of the increased after-storm electron intensity observed at the inner L-shells with no simultaneously seen increase in the geosynchronous distances are presented.     

X-ray measurements of the dark matter distribution in clusters of galaxies with Chandra

We investigate the dark matter distributions in the central region of two clusters of galaxies (A1835 and MKW3S) using Chandra data. N-body simulations in the standard cold dark matter (CDM) model predict the dark matter distribution shows a cuspy dark matter profile: ρ(r) ∝ r, with in the range 1–2, while observations of dwarf and low surface brightness galaxies seem to favor the presence of a relatively flat core: 0 <  < 1. To investigate the dark matter distributions in the central region of clusters of galaxies, we analyze the Chandra data of A1835 and MKW3S with a deprojection method. We derive the mass profiles without the assumption of analytical models. We examine the inner slope of derived mass profiles assuming the dark matter profile is described with a power-law expression. The values of the slope are 0.95 ± 0.10 for A1835 and 1.33 ± 0.12 for MKW3S within the radius of 200 kpc. These are consistent with the result of the CDM simulations. However, within the radius of 100 kpc, the value of is less than unity for A1835 (0.47 ± 0.31). Our result implies that the central dark matter profile of some clusters cannot be described by CDM halos.     

Observations of diffuse galactic positron annihilation radiation

The Galactic plane was scanned nearly three times by the UCSD/MIT Hard X-Ray and Low Energy Gamma-Ray Experiment on HEAO-1 from August 1977 through September 1978. Its Medium Energy Detectors were of the NaI/CsI phoswich type and operated over the 100 keV to 2 MeV range, with a 17° FWHM field of view and a 9% energy resolution at 511 keV. Sky maps for each epoch of observation were constructed in several energy bands. After subtraction of known point sources, a component associated with the galactic plane remains, whose spectrum is consistent with a power law and a positron annihilation spectrum. In the 333 to 583 keV energy band the flux is concentrated within ±35° of the galactic center, and the ratio of flux/radian (anticenter) is high, with a 2σ lower limit of 13. The parameters of the galactic center region's annihilation spectrum are positronium fraction of 0.9±0.1 and 511 keV flux of (2.0±0.7)×10⁻³ photons/cm²-sec-rad.     

The dark matter halo structure of “fossil” groups candidates

We investigated properties of four isolated giant elliptical galaxies with extended X-ray halo using ASCA data. The derived size of X-ray halo, X-ray luminosity, and gravitational mass of the dark halo are unusually large those of X-ray halo of a single galaxy, but are typical for X-ray halos of groups and poor clusters of galaxies. The measured temperatures and abundances of the X-ray halo gas in these galaxies are also similar to those of the groups and poor clusters. Based on these results we identified these galaxies as “isolated X-ray overluminous elliptical galaxy” (IOLEG). The radial profiles of dark halo in these objects were derived from X-ray data. It is found that some are similar to those of compact groups while others are the same as those of normal ellipticals. The dark halos of lOLEGs are thus indistinguishable from those of groups (and poor clusters), which appears to be consistent with a widely believed idea that lOLEGs are a product of dynamical evolution of a compact group. However, mass-to-light ratios of IOLEGs (M₂₀₀/L_B  100–1000) are far greater than those of Hickson compact groups M₂₀₀/L_B  40–60). Since it is hard to consider that total optical luminosity of a compact group decreases by an order of magnitude in the course of dynamical evolution, such difference in the observed mass-to-light ratio between IOLEGs and Hickson compact groups strongly suggests that most IOLEGs have not evolved from compact groups which are observed at present.     

Intensity profiles of the multiple scattering light near the cometary nucleus

The multiple scattering of solar radiation in the cometary atmosphere is treated with the method of successive scattering. Referring to in situ measurements of comet Halley about the size and spatial distributions of dust, the optical thickness τ₁ of dust has been estimated, i.e. τ₁=0.03 at wavelength λ=0.62μm in a quiet time, but τ₁=0.3 when the outbursts/jets occur. In the derivation of τ₁, optical properties of dust including a mixing ratio of absorbing to silicate grains, are determined based on the polarimetry of P/Halley at λ=0.62μm observed during the phase angles over Nov. 1985 to May 1986 at the Dodaira Station of Tokyo Astronomical Observatory.

It is found that a temporary enhancement of τ₁ leads an increase of the upward reflected intensity when the surface albedo A of the nucleus is less than 0.04, but the reverse is true when A>0.04. On the other hand, the intensity of the downward radiation at the surface of the nucleus always decreases as an increase of τ₁.     

Hard X-ray spectral properties and discovery of narrow annihilation line in the spectrum of Nova Muscae

The observations of X-ray Nova in Musca (GRS1124-684) by two coded mask telescopes on board GRANAT observatory provided spectral data in broad 3 – 1300 keV band. During these observations, spanned over a year, the Nova was detected in a three apparently different spectral states, corresponding to different epochs of the soft X-ray light curve: (1) A spectrum with two distinct components (soft, below 8 keV and hard power law tail with slope 2.5, detected up to 300 keV). The soft emission changed gradually with characteristic decay time around 30 days, while power law component exhibited strong variability on the time scales of several hours and decreased much more slowly. (2) A soft spectrum (without hard power law tail), observed during the “kick” of the soft X-ray light curve. (3) A hard power law spectrum with slope 2.2. Thus, while the 3 – 300 keV luminosity decreased by more than order of magnitude, the source passed through all spectral states known for galactic black hole candidates (Cyg X-1, GX339-4, 1E1740.7-2942, GRS1758-258 etc.).

On January 20–21 1991, the SIGMA telescope aboard GRANAT detected a relatively narrow variable emission line near 500 keV (Fig.1,2) with net flux ≈ 6 · 10⁻³ phot/s/cm², most probably related with electron-positron annihilation processes, occurring in the source /1–4/. Additional excess above power law continuum, centered around 200 keV, was found during this observation.     

Fractal analysis of short and intermediate GRB and the fireball model

In the present work the possibility of the fractal analysis application for GRB temporal profiles was studied. We have analysed the 4B revised BATSE catalog: temporal profiles of GRB with t₉₀ < 3 s (287 short and 100 intermediate) were studied on TTE data, a sample of 278 intermediate GRB with t₉₀  3 s were studied on DISCSC data. An analysis of the background fractal dimension distributions obtained using TTE and DISCSC data (143 and 110 background regions, respectively), indicates that for both datasets background fractal dimensions D_bgr = 1.5 that the fractal dimension distributions obtained by using these data can be processed simultaneously. The change of the fractal index D_bgr for Poisson statistics – dominated sets with different coefficients of error in counting (up to 10) was studied and D_bgr = 1.5. The ranges of fractal dimension (0.80  D  2.25 for short and 0.85  D  2.01 for intermediate GRB) are shifted over range for theoretical fractal curve (1 < D < 2) due to the finite detector time resolution. There are four subgroups in fractal dimension distribution for short GRB (D = 1.05 ± 0.03, D = 1.31 ± 0.05, D = 1.51 ± 0.04, D = 1.90 ± 0.03) and six subgroups for intermediate one (D = 1.05 ± 0.09, D = 1.24 ± 0.08, D = 1.44 ± 0.07, D = 1.51 ± 0.08, D = 1.64 ± 0.07, D = 1.91 ± 0.1). Time profiles with fractal dimension smaller then background can be obtained by using models with many short chaotic processes in sources, for example, fireball model with shock waves. The range of fractal dimensions for the modelled temporal profiles is 1.213  D  1.400, which can correspond to subgroups of short and intermediate GRB with D = 1.31 and D = 1.24; moreover, the fractal dimension of a simulated indented event and GRB990208 are equal within the error limits for some model parameters and it is possible to obtain smooth temporal profiles with D = D_bgr.     

The GPS measured SITEC caused by the very intense solar flare on July 14, 2000

This work studies the sudden increases in total electron content of the ionosphere caused by the very intense solar flare on July 14, 2000. Total electron content (TEC) data observed from a Global Positioning System (GPS) network are used to calculate the flare-induced TEC increment, δTEC_f, and variation rate, dTEC_f/dt. It is found that both dTEC_f/dt and δTEC_f are closely related with the solar zenith angles. To explain the observation results, we derived a simple relationship between the partial derivative of the flare-induced TEC, ∂TEC_f/∂t, which is a good approximation for dTEC_f/dt, and the solar zenith angle χ, as well as the effective flare radiation flux I_f, according to the well-known Chapman theory of ionization. The derived formula predicted that ∂TEC_f/∂t is proportional to I_f and inverse proportional to Chapman function ch(χ). This theoretical prediction not only explains the correlation of dTEC_f/dt and δTEC_f with χ as shown in our TEC observation, but also gives a way to deduce I_f from TEC observation of GPS network. Thus, the present work shows that GPS observation is a powerful tool in the observation and investigation of solar flare effects on the ionosphere, i.e., the sudden ionospheric disturbances, which is a significant phenomenon of space weather.     

Study of the luminosity-dependent cyclotron resonance energies of the binary X-ray pulsar 4U 0115+63 with RXTE

We report on the luminosity dependent change of the cyclotron resonance energy obtained from a transient X-ray pulsar, 4U 0115+63. Using RXTE data observed on 1999 March, we found that the fundamental resonance energy stayed constant (11 keV) when the source luminosity was above 5 × 10³⁷ erg s⁻¹. As the luminosity decreased below 5 × 10³⁷ erg s⁻¹, the fundamental resonance energy gradually increased up to 16 keV at 0.16 × 10³⁷ erg s⁻¹. The luminosity dependence of the resonance energy can be understood by the change of the accretion column height.     

LMXB luminosity function and connection to globular cluster in early type galaxies

We derive bias-corrected X-ray luminosity functions (XLFs) of LMXBs detected in 14 E and S0 galaxies observed with Chandra. After correcting for incompleteness, the individual XLFs are statistically consistent with a single power-law. A break at or near L_X,Eddington , as previously reported, is not required in any individual case. The combined XLF with a reduced error, however, suggests a possible break at L_X = 5 × 10³⁸ erg s⁻¹, which may be consistent with the Eddington luminosity of neutron stars with the largest possible mass (3 M), or of He-enriched neutron star binaries. We confirm that the total X-ray luminosity of LMXBs is correlated with the the near-IR luminosities, but the scatter exceeds that expected from measurement errors. The scatter in L_X(LMXB)/L_K appears to be correlated with the specific frequency of globular clusters, as reported earlier.

We cross-correlate X-ray binaries with globular clusters determined by ground-based optical and HST observations in 6 giant elliptical galaxies. With the largest sample reported so far (300 GC LMXBs with a 5:2 ratio between red and blue GCs), we compare their X-ray properties, such as X-ray hardness, XLF and L_X/L_B and find no statistically significance difference between different groups of LMXBs. Regardless of their association with GCs, both GC and field LMXBs appear to follow the radial profile of the optical halo light, rather than that of more extended GCs. This suggests that while metallicity is a primary factor in the formation of LMXBs in GCs, there may be a secondary factor (e.g., encounter rate) playing a non-negligible role.     

Compact X-ray binaries in and out of core collapsed globulars

We review new Chandra and HST observations of the core collapsed cluster NGC 6397 as a guide to understanding the compact binary (CB) populations in core collapse globulars. New cataclysmic variables (CVs) and main sequence chromospherically active binaries (ABs) have been identified, enabling a larger sample for comparison of the L_x, F_x/F_V and X-ray vs. optical color distributions. Comparison of the numbers of CBs with L_x  10³¹ erg s⁻¹ in 4 core collapse vs. 12 King model clusters reveals that the specific frequency S_X (number of CBs per unit cluster mass) is enhanced in core collapse clusters, even when normalized for their stellar encounter rate. Although core collapse is halted by the dynamical heating due to stellar (and binary) interaction with CBs in the core, we conclude that production of the hardest CBs – especially CVs – is enhanced during core collapse. NGC 6397 has its most luminous CVs nearest the cluster center, with two newly discovered very low luminosity (old, quiescent) CVs far from the core. The active binaries as well as neutron star systems (MSP and qLMXB) surround the central core. The overall CB population appears to be asymmetric about the cluster center, as in several other core collapse clusters observed with Chandra, suggesting still poorly understood scattering processes.     

Geminga pulsar observations with gamma-telescope Gamma-1

The Geminga light curve obtained with the “Gamma-1” telescope features two peaks separated by 0.5 ± 0.03 period. The light curve is pronounced for γ-quanta energies higher than 400 MeV. The pulsed flux upper limit (1σ) in the energy interval 50 – 300 MeV is 6·10⁻⁷ cm⁻²sec⁻¹. For energies >300 MeV the pulsed component power law spectrum has an exponent 1.1 _−0.3^+1.1 and an integral flux (1.1±0.3)·10⁻⁶ cm⁻²sec⁻¹.