Probing magnetic turbulence by synchrotron polarimetry: statistics and structure of magnetic fields from Stokes correlators

We describe a novel technique for probing the statistical properties of cosmic magnetic fields based on radio polarimetry data. Second-order magnetic field statistics like the power spectrum cannot always distinguish between magnetic fields with essentially different spatial structure. Synchrotron polarimetry naturally allows certain fourth-order magnetic field statistics to be inferred from observational data, which lifts this degeneracy and can thereby help us gain a better picture of the structure of the cosmic fields and test theoretical scenarios describing magnetic turbulence. In this work we show that a fourth-order correlator of specific physical interest, the tension force spectrum, can be recovered from the polarized synchrotron emission data. We develop an estimator for this quantity based on polarized emission observations in the Faraday rotation free frequency regime. We consider two cases: a statistically isotropic field distribution, and a statistically isotropic field superimposed on a weak mean field. In both cases the tension force power spectrum is measurable; in the latter case, the magnetic power spectrum may also be obtainable. The method is exact in the idealized case of a homogeneous relativistic electron distribution that has a power-law energy spectrum with a spectral index of   p = 3  , and assumes statistical isotropy of the turbulent field. We carry out numerical tests of our method using synthetic polarized emission data generated from numerically simulated magnetic fields. We show that the method is valid, that it is not prohibitively sensitive to the value of the electron spectral index p , and that the observed tension force spectrum allows one to distinguish between e.g. a randomly tangled magnetic field (a default assumption in many studies) and a field organized in folded flux sheets or filaments.     

A new moment method for continuum radiative transfer in cosmological re-ionization

We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation, closed by an Eddington tensor computed using a long characteristics (LC) method. We show that traditional short characteristics and the optically thin approximation are inappropriate for computing Eddington factors for the problem of cosmological re-ionization. We evolve the non-equilibrium ionization field via an efficient and accurate (errors <1 per cent) technique that switches between fully implicit or explicit finite differencing depending on whether the local time-scales are long or short compared to the time-step. We tailor our code for the problem of cosmological re-ionization. In tests, the code conserves photons, accurately treats cosmological effects and reproduces analytic Strömgren sphere solutions. Its chief weakness is that the computation time for the LC calculation scales relatively poorly compared to other techniques  ( t _LC∝ N ^∼1.5_cells)  ; however, we mitigate this by only recomputing the Eddington tensor when the radiation field changes substantially. Our technique makes almost no physical approximations, so it provides a way to benchmark faster but more approximate techniques. It can readily be extended to evolve multiple frequencies, though we do not do so here. Finally, we note that our method is generally applicable to any problem involving the transfer of continuum radiation through a periodic volume.     

Interpreting the transmission windows of distant quasars

We propose the apparent shrinking criterion (ASC) to interpret the spatial extent, R _w, of transmitted flux windows in the absorption spectra of high- z quasars. The ASC can discriminate between the two regimes in which R _w corresponds either to the physical size,   R _{H  ii}   , of the quasar H  ii region or to the distance,   R ^max_w  , at which the transmitted flux drops to  =0.1  and a Gunn–Peterson (GP) trough appears. In the first case [H  ii region (HR) regime], one can determine the intergalactic medium mean H  i fraction,   x _H I  ; in the second [proximity region (PR) regime], the value of R _w allows one to measure the local photoionization rate and the local enhancement of the photoionization rate,  Γ_G  , due to nearby/intervening galaxies. The ASC has been tested against radiative transfer+smoothed particle hydrodynamics numerical simulations, and applied to a sample of 15 high-   z ( z > 5.8  ) quasar spectra. All sample quasars are found to be in the PR regime; hence, their observed spectral properties (inner flux profile, extent of transmission window) cannot reliably constrain the value of   x _{H  i}   . Four sample quasars show evidence for a local enhancement (up to 50 per cent) in the local photoionization rate possibly produced by a galaxy overdensity. We discuss the possible interpretations and uncertainties of this result.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

Software timing calibration of the ARGO-YBJ detector

The ARGO-YBJ experiment is mainly devoted to search for astronomical gamma sources. The arrival direction of air showers is reconstructed thanks to the times measured by the pixels of the detector. Therefore, the timing calibration of the detector pixels is crucial in order to get the best angular resolution and pointing accuracy. Because of the large number of pixels a hardware timing calibration is practically impossible. Therefore an off-line software calibration has been adopted. Here, the details of the procedure and the results are presented.     

Interstellar extinction toward the ultradeep Galactic field of the Chandra observatory from 2MASS data

We have measured the interstellar extinction in the region of ultradeep Galactic-field observations by the Chandra telescope (l ^II, b ^II) ≈ 0.1–1.42 using photometric data from the 2MASS infrared allsky survey. The angular resolution of our interstellar extinction map is 1′.8. We show that the interstellar extinction has a minimum, A _V ～ 3.4, near the center of the Chandra field of view and increases to A _V ～ 5.8–6 at the edge of the field of view. In addition, we show that the bulk of the extinction is gained in the Galactic disk and is approximately the same for all bulge stars. Our results will be subsequently used to process the Chandra data and to estimate the properties of the stellar population in this region.