Discovery of pulsations in the X-ray transient 4U 1901+03

We describe observations of the 2003 outburst of the hard-spectrum X-ray transient 4U 1901+03 with the Rossi X-ray Timing Explorer. The outburst was first detected in 2003 February by the All-Sky Monitor, and reached a peak 2.5-25 keV flux of 8x10^-9 ergs/cm^2/s (around 240 mCrab). The only other known outburst occurred 32.2 yr earlier, likely the longest presently known recurrence time for any X-ray transient. Proportional Counter Array (PCA) observations over the 5-month duration of the 2003 outburst revealed a 2.763 s pulsar in a 22.58 d orbit. The detection of pulsations down to a flux of 3x10^-11 ergs/cm^2/s (2.5-25 keV), along with the inferred long-term accretion rate of 8.1x10^-11 M_sun/yr (assuming a distance of 10 kpc) suggests that the surface magnetic field strength is below ~5x10^11 G. The corresponding cyclotron energy is thus below 4 keV, consistent with the non-detection of resonance features at high energies. Although we could not unambiguously identify the optical counterpart, the lack of a bright IR candidate within the 1' RXTE error circle rules out a supergiant mass donor. The neutron star in 4U 1901+03 probably accretes from the wind of a main-sequence O-B star, like most other high-mass binary X-ray pulsars. The almost circular orbit e=0.036 confirms the system's membership in a growing class of wide, low-eccentricity systems in which the neutron stars may have received much smaller kicks as a result of their natal supernova explosions.Comment: 7 pages, 6 figures, accepted by ApJ. Very minor addition in response to referee's comment; updated author affiliatio

Parameterization of the Angular Distribution of Gamma Rays Produced by p-p Interaction in Astronomical Environment

We present the angular distribution of gamma rays produced by proton-proton interactions in parameterized formulae to facilitate calculations in astrophysical environments. The parameterization is derived from Monte Carlo simulations of the up-to-date proton-proton interaction model by Kamae et al. (2005) and its extension by Kamae et al. (2006). This model includes the logarithmically rising inelastic cross section, the diffraction dissociation process and Feynman scaling violation. The extension adds two baryon resonance contributions: one representing the Delta(1232) and the other representing multiple resonances around 1600 MeV/c^2. We demonstrate the use of the formulae by calculating the predicted gamma-ray spectrum for two different cases: the first is a pencil beam of protons following a power law and the second is a fanned proton jet with a Gaussian intensity profile impinging on the surrounding material. In both cases we find that the predicted gamma-ray spectrum to be dependent on the viewing angle.Comment: 8 pages, 7 figures, figure 7 updated, accepted for publication in ApJ, text updated to match changes by the editor, two refs updated from preprints to full journal

The Plasma Structure of the Cygnus Loop from the Northeastern Rim to the Southwestern Rim

The Cygnus Loop was observed from the northeast to the southwest with XMM-Newton. We divided the observed region into two parts, the north path and the south path, and studied the X-ray spectra along two paths. The spectra can be well fitted either by a one-component non-equilibrium ionization (NEI) model or by a two-component NEI model. The rim regions can be well fitted by a one-component model with relatively low \kTe whose metal abundances are sub-solar (0.1--0.2). The major part of the paths requires a two-component model. Due to projection effects, we concluded that the low kTe (about 0.2 keV) component surrounds the high kTe (about 0.6 keV) component, with the latter having relatively high metal abundances (about 5 times solar). Since the Cygnus Loop is thought to originate in a cavity explosion, the low-kTe component originates from the cavity wall while the high-kTe component originates from the ejecta. The flux of the cavity wall component shows a large variation along our path. We found it to be very thin in the south-west region, suggesting a blowout along our line of sight. The metal distribution inside the ejecta shows non-uniformity, depending on the element. O, Ne and Mg are relatively more abundant in the outer region while Si, S and Fe are concentrated in the inner region, with all metals showing strong asymmetry. This observational evidence implies an asymmetric explosion of the progenitor star. The abundance of the ejecta also indicates the progenitor star to be about 15 M_sun.Comment: 24 pages, 9 figures, Astrophysical Journal in pres

Recovery of the orbital parameters and pulse evolution of V0332+53 during a huge outburst

The high mass X-ray binary (HMXB) V0332+53 became active at the end of 2004 and the outburst was observed at hard X-rays by RXTE and INTEGRAL. Based on these hard X-ray observations, the orbital parameters are measured through fitting the Doppler-shifted spin periods. The derived orbital period and eccentricity are consistent with those of Stella et al. (1985) obtained from EXOSAT observations, whereas the projected semimajor axis and the periastron longitude are found to have changed from 48$\pm$4 to 86$^{+6}_{-10}$ lt-s and from 313$^{\circ}$$\pm$10 to 283$^{\circ}$$\pm$14, respectively. This would indicate an angular speed of $\geq$ 1.5$^{\circ}$$\pm$0.8 yr$^{-1}$ for rotation of the orbit over the past 21 years. The periastron passage time of MJD 53367$\pm$1 is just around the time when the intensity reached maximum and an orbital period earlier is the time when the outburst started. This correlation resembles the behavior of a Type I outburst. During outburst the source spun up with a rate of 8.01$^{+1.00}_{-1.14}$$\times10^{-6}$ s day$^{-1}$. The evolution of pulse profile is highly intensity dependent. The separation of double pulses remained almost constant ($\sim$ 0.47) when the source was bright, and dropped to 0.37 within $\leq$ 3 days as the source became weaker. The pulse evolution of V0332+53 may correlate to the change in dominance of the emission between fan-beam and pencil-beam mechanisms.Comment: 13 pages, 3 figures, accepted for publication in ApJ

Chandra view of Kes 79: a nearly isothermal SNR with rich spatial structure

A 30 ks \chandra ACIS-I observation of Kes 79 reveals rich spatial structures, including many filaments, three partial shells, a loop and a ``protrusion''. Most of them have corresponding radio features. Regardless of the different results from two non-equilibrium ionization (NEI) codes, temperatures of different parts of the remnant are all around 0.7 keV, which is surprisingly constant for a remnant with such rich structure. If thermal conduction is responsible for smoothing the temperature gradient, a lower limit on the thermal conductivity of $\sim$ 1/10 of the Spitzer value can be derived. Thus, thermal conduction may play an important role in the evolution of at least some SNRs. No spectral signature of the ejecta is found, which suggests the ejecta material has been well mixed with the ambient medium. From the morphology and the spectral properties, we suggest the bright inner shell is a wind-driven shell (WDS) overtaken by the blast wave (the outer shell) and estimate the age of the remnant to be $\sim$ 6 kyr for the assumed dynamics. Projection is also required to explain the complicated morphology of Kes 79.Comment: 12 pages, 6 figures (3 in color), ApJ, in press, April 20, 200

Asymmetric Ejecta Distribution of the Cygnus Loop revealed with Suzaku

We observed a linearly sliced area of the Cygnus Loop from the north-east to the south-west with Suzaku in seven pointings. After dividing the entire fields of view (FOV) into 119 cells, we extracted spectra from all of the cells and performed spectral analysis for them. We then applied both one- and two-component non-equilibrium ionization (NEI) models for all of the spectra, finding that almost all were significantly better fitted by the two-component NEI model rather than the one-component NEI model. Judging from the abundances, the high-kT_e component must be the ejecta component, while the low-kT_e component comes from the swept-up matter. Therefore, the ejecta turn out to be distributed inside a large area (at least our FOV) of the Cygnus Loop. We divided the entire FOV into northern and southern parts, and found that the ejecta distributions were asymmetric to the geometric center: the ejecta of Si, S, and Fe seem to be distributed more in the south than in the north of the Cygnus Loop by a factor of about 2. The degree of ejecta-asymmetry is consistent with that expected by recent supernova explosion models.Comment: 14 pages, 4 figure

Abundance Inhomogeneity in the Northeastern Rim of the Cygnus Loop Revealed by Suzaku Observatory

We present the results of a spatially resolved spectral analysis from four Suzaku observations covering the northeastern rim of the Cygnus Loop. A two-kT_e non-ionization equilibrium (NEI) model fairly well represents our data, which confirms the NEI condition of the plasma there. The metal abundances are depleted relative to the solar values almost everywhere in our field of view. We find abundance inhomogeneities across the field: the northernmost region (Region A) has enhanced absolute abundances compared with other regions. In addition, the relative abundances of Mg/O and Fe/O in Region A are lower than the solar values, while those in the other regions are twice higher than the solar values. As far as we are concerned, neither a circumstellar medium, fragments of ejecta, nor abundance inhomogeneities of the local interstellar medium around the Cygnus Loop can explain the relatively enhanced abundance in Region A. This point is left as an open question for future work.Comment: 13 pages, 6 figure

Changes in the long term intensity variations in Cyg X-2 and LMC X-3

We report the detection of changes in the long-term intensity variations in two X-ray binaries, Cyg X-2 and LMC X-3. In this work, we have used the long-term light curves obtained with the All-Sky Monitors (ASMs) of the Rossi X-Ray Timing Explorer (RXTE), Ginga, Ariel 5, and Vela 5B and the scanning modulation collimator of HEAO 1. It is found that in the light curves of both the sources, obtained with these instruments at various times over the last 30 years, more than one periodic or quasi-periodic component is always present. The multiple prominent peaks in the periodograms have frequencies unrelated to each other. In Cyg X-2, RXTE-ASM data show strong peaks at 40.4 and 68.8 days, and Ginga-ASM data show strong peaks at 53.7 and 61.3 days. Multiple peaks are also observed in LMC X-3. The various strong peaks in the periodograms of LMC X-3 appear at 104, 169, and 216 days (observed with RXTE-ASM) and 105, 214, and 328 days (observed with Ginga-ASM). The present results, when compared with the earlier observations of periodicities in these two systems, demonstrate the absence of any stable long period. The 78 day periodicity detected earlier in Cyg X-2 was probably due to the short time base in the RXTE data that were used, and the periodicity of 198 days in LMC X-3 was due to a relatively short duration of observation with HEAO 1.Comment: 11 pages, 7 postscript figures include