CO and C_2 Absorption Toward W40 IRS 1a

The H II region W40 harbors a small group of young, hot stars behind roughly 9 magnitudes of visual extinction. We have detected gaseous carbon monoxide (CO) and diatomic carbon (C_2) in absorption toward the star W40 IRS 1a. The 2-0 R0, R1, and R2 lines of 12CO at 2.3 micron were measured using the CSHELL on the NASA IR Telescope Facility (with upper limits placed on R3, R4, and R5) yielding an N_CO of (1.1 +/- 0.2) x 10^18 cm^-2. Excitation analysis indicates T_kin > 7 K. The Phillips system of C_2 transitions near 8775 Ang. was measured using the Kitt Peak 4-m telescope and echelle spectrometer. Radiative pumping models indicate a total C_2 column density of (7.0 +/- 0.4) x 10^14 cm^-2, two excitation temperatures (39 and 126 K), and a total gas density of n ~ 250 cm^-3. The CO ice band at 4.7 micron was not detected, placing an upper limit on the CO depletion of delta < 1 %. We postulate that the sightline has multiple translucent components and is associated with the W40 molecular cloud. Our data for W40 IRS 1a, coupled with other sightlines, shows that the ratio of CO/C_2 increases from diffuse through translucent environs. Finally, we show that the hydrogen to dust ratio seems to remain constant from diffuse to dense environments, while the CO to dust ratio apparently does not.Comment: To appear in The Astrophysical Journal 17 pages total, 5 figures Also available at http://casa.colorado.edu/~shuping/research/w40/w40.htm

Diffusive propagation of UHECR and the propagation theorem

We present a detailed analytical study of the propagation of ultra high energy (UHE) particles in extragalactic magnetic fields. The crucial parameter which affects the diffuse spectrum is the separation between sources. In the case of a uniform distribution of sources with a separation between them much smaller than all characteristic propagation lengths, the diffuse spectrum of UHE particles has a {\em universal} form, independent of the mode of propagation. This statement has a status of theorem. The proof is obtained using the particle number conservation during propagation, and also using the kinetic equation for the propagation of UHE particles. This theorem can be also proved with the help of the diffusion equation. In particular, it is shown numerically, how the diffuse fluxes converge to this universal spectrum, when the separation between sources diminishes. We study also the analytic solution of the diffusion equation in weak and strong magnetic fields with energy losses taken into account. In the case of strong magnetic fields and for a separation between sources large enough, the GZK cutoff can practically disappear, as it has been found early in numerical simulations. In practice, however, the source luminosities required are too large for this possibility.Comment: 16 pages, 13 eps figures, discussion of the absence of the GZK cut-off in strong magnetic field added, a misprint in figure 6 corrected, version accepted for publication in Ap

Diffusion of Cosmic Rays in the Expanding Universe. II. Energy Spectra of Ultra-High Energy Cosmic Rays

We consider the astrophysical implications of the diffusion equation solution for Ultra-High Energy Cosmic Rays (UHECR) in the expanding universe, obtained in paper I (V.Berezinsky & A.Gazizov, ApJ 643 (2006) 8). The UHECR spectra are calculated in a model with sources located in vertices of the cubic grid with a linear constant (source separation) d. The calculations are performed for various magnetic field configurations (B_c,l_c), where l_c is the basic scale of the turbulence and B_c is the coherent magnetic field on this scale. The main purpose of these calculations is to demonstrate the validity of the solution obtained in paper I and to compare this solution with the Syrovatsky solution used in previous works. The Syrovatsky solution must be necessarily embedded in the static cosmological model. The formal comparison of the two solutions with all parameters being fixed identically reveals the appreciable discrepancies between two spectra. These discrepancies are less if in both models the different sets of the best-fit parameters are used.Comment: 17 pages, 10 figure

EGMF Constraints from Simultaneous GeV-TeV Observations of Blazars

Attenuation of the TeV gamma-ray flux from distant blazars through pair production with extragalactic background light leads to the development of electromagnetic cascades and subsequent, lower energy, GeV secondary gamma-ray emission. Due to the deflection of VHE cascade electrons by extragalactic magnetic fields (EGMF), the spectral shape of this arriving cascade gamma-ray emission is dependent on the strength of the EGMF. Thus, the spectral shape of the GeV-TeV emission from blazars has the potential to probe the EGMF strength along the line of sight to the object. We investigate constraints on the EGMF derived from observations of blazars for which TeV observations simultaneous with those by the Fermi telescope were reported. We study the dependence of the EGMF bound on the hidden assumptions it rests upon. We select blazar objects for which simultaneous Fermi/LAT GeV and Veritas, MAGIC or HESS TeV emission have been published. We model the development of electromagnetic cascades along the gamma-ray beams from these sources using Monte Carlo simulations, including the calculation of the temporal delay incurred by cascade photons, relative to the light propagation time of direct gamma-rays from the source. Constraints on EGMF could be derived from the simultaneous GeV-TeV data on the blazars RGB J0710+591, 1ES 0229+200, and 1ES 1218+304. The measured source flux level in the GeV band is lower than the expected cascade component calculated under the assumption of zero EGMF. Assuming that the reason for the suppression of the cascade component is the extended nature of the cascade emission, we find that B>10^{-15} G (assuming EGMF correlation length of ~1 Mpc) is consistent with the data. Alternatively, the assumption that the suppression of the cascade emission is caused by the time delay of the cascade photons the data are consistent with B>10^{-17} G for the same correlation length.Comment: 9 pages, 9 figure

The critical velocity effect as a cause for the H\alpha emission from the Magellanic stream

Observations show significant H\alpha-emissions in the Galactic halo near the edges of cold gas clouds of the Magellanic Stream. The source for the ionization of the cold gas is still a widely open question. In our paper we discuss the critical velocity effect as a possible explanation for the observed H\alpha-emission. The critical velocity effect can yield a fast ionization of cold gas if this neutral gas passes through a magnetized plasma under suitable conditions. We show that for parameters that are typical for the Magellanic Stream the critical velocity effect has to be considered as a possible ionization source of high relevance.Comment: 9 pages, 2 figures. accepted, to appear in The Astrophysical Journa

Simulations of small-scale turbulent dynamo

We report an extensive numerical study of the small-scale turbulent dynamo at large magnetic Prandtl numbers Pm. A Pm scan is given for the model case of low-Reynolds-number turbulence. We concentrate on three topics: magnetic-energy spectra and saturation levels, the structure of the field lines, and the field-strength distribution. The main results are (1) the folded structure (direction reversals at the resistive scale, field lines curved at the scale of the flow) persists from the kinematic to the nonlinear regime; (2) the field distribution is self-similar and appears to be lognormal during the kinematic regime and exponential in the saturated state; and (3) the bulk of the magnetic energy is at the resistive scale in the kinematic regime and remains there after saturation, although the spectrum becomes much shallower. We propose an analytical model of saturation based on the idea of partial two-dimensionalization of the velocity gradients with respect to the local direction of the magnetic folds. The model-predicted spectra are in excellent agreement with numerical results. Comparisons with large-Re, moderate-Pm runs are carried out to confirm the relevance of these results. New features at large Re are elongation of the folds in the nonlinear regime from the viscous scale to the box scale and the presence of an intermediate nonlinear stage of slower-than-exponential magnetic-energy growth accompanied by an increase of the resistive scale and partial suppression of the kinetic-energy spectrum in the inertial range. Numerical results for the saturated state do not support scale-by-scale equipartition between magnetic and kinetic energies, with a definite excess of magnetic energy at small scales. A physical picture of the saturated state is proposed.Comment: aastex using emulateapj; 32 pages, final published version; a pdf file (4Mb) of the paper containing better-quality versions of figs. 5, 8, 12, 15, 17 is available from http://www.damtp.cam.ac.uk/user/as629 or by email upon request

Anti-GZK effect in Ultra High Energy Cosmic Rays diffusive propagation

We discuss the antiGZK effect in the diffusive propagation of ultra high energy protons in intergalactic magnetic fields, which consists in a jump-like increase of the maximum distance from which ultra high energy protons can reach an observer. The position of this jump, $E_j \approx 2\times 10^{18}$ eV, is determined exclusively by energy losses (transition from adiabatic to pair-production energy losses) and it is independent of the diffusion parameters. The diffuse spectrum presents a low-energy steepening approximately at this energy, which is very close to the position of the second knee observed in the cosmic ray spectrum. The dip, seen in the universal spectrum as a signature of the interaction with the cosmic microwave background radiation, is also present in the case of diffusive propagation in magnetic fields.Comment: 19 pages, 8 figures, a typo correcte

Magnetic Reconnection Triggered by the Parker Instability in the Galaxy: Two-Dimensional Numerical Magnetohydrodynamic Simulations and Application to the Origin of X-Ray Gas in the Galactic Halo

We propose the Galactic flare model for the origin of the X-ray gas in the Galactic halo. For this purpose, we examine the magnetic reconnection triggered by Parker instability (magnetic buoyancy instability), by performing the two-dimensional resistive numerical magnetohydrodynamic simulations. As a result of numerical simulations, the system evolves as following phases: Parker instability occurs in the Galactic disk. In the nonlinear phase of Parker instability, the magnetic loop inflates from the Galactic disk into the Galactic halo, and collides with the anti-parallel magnetic field, so that the current sheets are created in the Galactic halo. The tearing instability occurs, and creates the plasmoids (magnetic islands). Just after the plasmoid ejection, further current-sheet thinning occurs in the sheet, and the anomalous resistivity sets in. Petschek reconnection starts, and heats the gas quickly in the Galactic halo. It also creates the slow and fast shock regions in the Galactic halo. The magnetic field ($B\sim 3 \mu$G), for example, can heat the gas ($n\sim 10^{-3}$ cm$^{-3}$) to temperature of $\sim 10^6$ K via the reconnection in the Galactic halo. The gas is accelerated to Alfv\'en velocity ($\sim 300$ km s$^{-1}$). Such high velocity jets are the evidence of the Galactic flare model we present in this paper, if the Doppler shift of the bipolar jet is detected in the Galactic halo. Full size figures are available at http://www.kwasan.kyoto-u.ac.jp/~tanuma/study/ApJ2002/ApJ2002.htmlComment: 13 pages, 12 figures, uses emulateapj.sty, accepted by Ap

The contribution of the Unresolved Extragalactic Radio Sources to the Brightness Temperature of the sky

The contribution of the Unresolved Extragalactic Radio Sources to the diffuse brightness of the sky was evaluated using the source number - flux measurements available in literature. We first optimized the fitting function of the data based on number counts distribution. We then computed the brightness temperature at various frequencies from 151 MHz to 8440 MHz and derived its spectral dependence. As expected the frequency dependence can be described by a power law with a spectral index $\gamma \simeq -2.7$, in agreement with the flux emitted by the {\it steep spectrum} sources. The contribution of {\it flat spectrum} sources becomes relevant at frequencies above several GHz. Using the data available in literature we improved our knowledge of the brightness of the unresolved extragalactic radio sources. The results obtained have general validity and they can be used to disentangle the various contributions of the sky brightness and to evaluate the CMB temperature.Comment: Accepted for publication in the Astrophysical Journa

Diffuse continuum gamma rays from the Galaxy

A new study of the diffuse Galactic gamma-ray continuum radiation is presented, using a cosmic-ray propagation model which includes nucleons, antiprotons, electrons, positrons, and synchrotron radiation. Our treatment of the inverse Compton (IC) scattering includes the effect of anisotropic scattering in the Galactic interstellar radiation field (ISRF) and a new evaluation of the ISRF itself. Models based on locally measured electron and nucleon spectra and synchrotron constraints are consistent with gamma-ray measurements in the 30-500 MeV range, but outside this range excesses are apparent. A harder nucleon spectrum is considered but fitting to gamma rays causes it to violate limits from positrons and antiprotons. A harder interstellar electron spectrum allows the gamma-ray spectrum to be fitted above 1 GeV as well, and this can be further improved when combined with a modified nucleon spectrum which still respects the limits imposed by antiprotons and positrons. A large electron/IC halo is proposed which reproduces well the high-latitude variation of gamma-ray emission. The halo contribution of Galactic emission to the high-latitude gamma-ray intensity is large, with implications for the study of the diffuse extragalactic component and signatures of dark matter. The constraints provided by the radio synchrotron spectral index do not allow all of the <30 MeV gamma-ray emission to be explained in terms of a steep electron spectrum unless this takes the form of a sharp upturn below 200 MeV. This leads us to prefer a source population as the origin of the excess low-energy gamma rays.Comment: Final version accepted for publication in The Astrophysical Journal (vol. 537, July 10, 2000 issue); Many Updates; 20 pages including 49 ps-figures, uses emulateapj.sty. More details can be found at http://www.gamma.mpe-garching.mpg.de/~aws/aws.htm