Separation of instrumental and astrophysical foregrounds for mapping cosmic microwave background anisotropies

For the most sensitive present and future experiments dedicated to cosmic microwave background (CMB) anisotropy observations, the identification and separation of signals coming from different sources is an important step in the data analysis. This problem of the restitution of signals from the observation of their mixture is classically called 'component separation' in CMB mapping. In this paper, we address the general problem of separating, for millimetre-wave sky-mapping applications, components which include not only astrophysical emissions in two-dimensional maps, but also one-dimensional instrumental effects in the data streams. We show that component separation methods can be adapted to separate simultaneously both astrophysical emissions and components coming from time-dependent foreground signals which originate from the instrument itself. Such general methods can be used for the optimal processing of low-redundancy observations where multi-channel observations are a precious tool to remove systematic effects, as may be the case for the Planck mission.     

Radio Cherenkov signals from the Moon: Neutrinos and cosmic rays

Neutrino production of radio Cherenkov signals in the Moon is the object of radio telescope observations. Depending on the energy range and detection parameters, the dominant contribution to the neutrino signal may come from interactions of the neutrino on the Moon facing the telescope, rather than neutrinos that have traversed a portion of the Moon. Using the approximate analytic expression of the effective lunar aperture from a recent paper by Gayley, Mutel and Jaeger, we evaluate the background from cosmic ray interactions in the lunar regolith. We also consider the modifications to the effective lunar aperture from generic non-standard model neutrino interactions. A background to neutrino signals are radio Cherenkov signals from cosmic ray interactions. For cosmogenic neutrino fluxes, neutrino signals will be difficult to observe because of low neutrino flux at the high energy end and large cosmic ray background in the lower energy range considered here. We show that lunar radio detection of neutrino interactions is best suited to constrain or measure neutrinos from astrophysical sources and probe non-standard neutrino-nucleon interactions such as microscopic black hole production.     

Exact reconstruction with directional wavelets on the sphere

A new formalism is derived for the analysis and exact reconstruction of band-limited signals on the sphere with directional wavelets. It represents an evolution of a previously developed wavelet formalism developed by Antoine & Vandergheynst and Wiaux et al. The translations of the wavelets at any point on the sphere and their proper rotations are still defined through the continuous three-dimensional rotations. The dilations of the wavelets are directly defined in harmonic space through a new kernel dilation, which is a modification of an existing harmonic dilation. A family of factorized steerable functions with compact harmonic support which are suitable for this kernel dilation are first identified. A scale-discretized wavelet formalism is then derived, relying on this dilation. The discrete nature of the analysis scales allows the exact reconstruction of band-limited signals. A corresponding exact multi-resolution algorithm is finally described and an implementation is tested. The formalism is of interest notably for the denoising or the deconvolution of signals on the sphere with a sparse expansion in wavelets. In astrophysics, it finds a particular application for the identification of localized directional features in the cosmic microwave background data, such as the imprint of topological defects, in particular, cosmic strings, and for their reconstruction after separation from the other signal components.     

Signals of statistical anisotropy in WMAP foreground-cleaned maps

Recently, a symmetry-based method to test for statistical isotropy of the cosmic microwave background was developed. We apply the method to template-cleaned 3- and 5-years Wilkinson Microwave Anisotropy Probe - Differencing Assembly maps. We examine a wide range of angular multipoles from  2 < l < 300  . The analysis detects statistically significant signals of anisotropy inconsistent with an isotropic cosmic microwave background in some of the foreground-cleaned maps. We are unable to resolve whether the anomalies have a cosmological, local astrophysical or instrumental origin. Assuming the anisotropy arises due to residual foreground contamination, we estimate the residual foreground power in the maps. For the W -band maps, we also find a highly improbable degree of isotropy we cannot explain. We speculate that excess isotropy may be caused by faulty modelling of detector noise.     

Wavelet analysis and the detection of non-Gaussianity in the cosmic microwave background

We investigate the use of wavelet transforms in detecting and characterizing non-Gaussian structure in maps of the cosmic microwave background (CMB). We apply the method to simulated maps of the KaiserStebbins effect resulting from cosmic strings, on to which Gaussian signals of varying amplitudes are superposed. We find that the method significantly outperforms standard techniques based on measuring the moments of the pixel temperature distribution. We also compare the results with those obtained using techniques based on Minkowski functionals, and we again find the wavelet method to be superior. In particular, using the wavelet technique, we find that it is possible to detect non-Gaussianity even in the presence of a superposed Gaussian signal with 3 times the rms amplitude of the original cosmic string map. We also find that the wavelet technique is useful in characterizing the angular scales at which the non-Gaussian signal occurs.     

Neural networks and the separation of cosmic microwave background and astrophysical signals in sky maps

We implement an independent component analysis (ICA) algorithm to separate signals of different origin in sky maps at several frequencies. Owing to its self-organizing capability, it works without prior assumptions on either the frequency dependence or the angular power spectrum of the various signals; rather, it learns directly from the input data how to identify the statistically independent components, on the assumption that all but, at most, one of the components have non-Gaussian distributions.
We have applied the ICA algorithm to simulated patches of the sky at the four frequencies (30, 44, 70 and 100 GHz) used by the Low Frequency Instrument of the European Space Agency's Planck satellite. Simulations include the cosmic microwave background (CMB), the synchrotron and thermal dust emissions, and extragalactic radio sources. The effects of the angular response functions of the detectors and of instrumental noise have been ignored in this first exploratory study. The ICA algorithm reconstructs the spatial distribution of each component with rms errors of about 1 per cent for the CMB, and 10 per cent for the much weaker Galactic components. Radio sources are almost completely recovered down to a flux limit corresponding to ≃0.7 σ _CMB, where σ _CMB is the rms level of the CMB fluctuations. The signal recovered has equal quality on all scales larger than the pixel size. In addition, we show that for the strongest components (CMB and radio sources) the frequency scaling is recovered with per cent precision. Thus, algorithms of the type presented here appear to be very promising tools for component separation. On the other hand, we have been dealing here with a highly idealized situation. Work to include instrumental noise, the effect of different resolving powers at different frequencies and a more complete and realistic characterization of astrophysical foregrounds is in progress.     

Multidetector multicomponent spectral matching and applications for cosmic microwave background data analysis

We present a new method for analysing multidetector maps containing several astrophysical components. Our method, based on matching the data to a model in the spectral domain, permits us to estimate jointly the spatial power spectra of the components and of the noise, as well as their mixing coefficients. It is of particular relevance for analysis of millimetre-wave maps of cosmic microwave background (CMB) anisotropies.     

Constraining Galactic pγ Interactions with Cosmic Ray Electron and Positron Spectra

High energy protons produced by various sources of cosmic rays, e.g., supernova remnants, pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, participate in Pγ and pp interactions. Although pp interactions may be the dominant mechanism in our Galaxy, it is unclear how important pγ process is. We show that the upper bound on the fraction of total number of protons participating in pγ interactions inside all Galactic astrophysical sources of cosmic rays is 10%.     

Imaging the heart of astrophysical objects with optical long-baseline interferometry

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: (1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; (2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long-baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image-reconstruction process was ??blind?? in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image-reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with six to ten telescopes, should be well suited to provide images of complex sceneries.     

A Faraday rotation template for the Galactic sky

Using a set of compilations of measurements for extragalactic radio sources, we construct all-sky maps of the Faraday rotation produced by the Galactic magnetic field. In order to generate the maps, we treat the radio source positions as a kind of 'mask' and construct combinations of spherical harmonic modes that are orthogonal on the masked sky. As long as relatively small multipoles are used, the resulting maps are quite stable to changes in the selection criteria for the sources, and show clearly the structure of the local Galactic magnetic field. We also suggest the use of these maps as templates for cosmic microwave background (CMB) foreground analysis, illustrating the idea with a cross-correlation analysis between the Wilkinson Microwave Anisotropy Probe ( WMAP ) data and our maps. We find a significant cross-correlation, indicating the presence of a significant residual contamination.     

The IAA cosmic dust laboratory: Experimental scattering matrices of clay particles

We present the first results of measurements on solid particles performed at the Instituto de Astrof?´sica de Andaluc?´a (IAA) cosmic dust laboratory located in Granada, Spain. The laboratory apparatus measures the complete scattering matrix as a function of the scattering angle of aerosol particles. The measurements can be performed at a wavelength (λ) of 483, 488, 520, 568, or 647 nm in the scattering angle range from 3° to 177°. Results of special test experiments are presented which show that our experimental results for scattering matrices are not significantly contaminated by multiple scattering and that the sizes/shapes of the particles do not change during the measurements. Moreover, the measured scattering matrix for a sample of green clay particles is compared with measurements previously performed in the Amsterdam light scattering setup for the same sample. New measurements on a white clay sample at 488 and 647 nm are also presented. The apparatus is devoted to experimentally studying the angle dependence of scattering matrices of dust samples of astrophysical interest. Moreover, there is a great interest in similar studies of aerosols that can affect the radiative balance of the atmosphere of the Earth and other planets such as silicates, desert dust, volcanic ashes, and carbon soot particles.     

移动星团视差研究方法的改进

移动星团视差的精确测定对于宇宙距离尺度以及许多天体物理关系的研究起重要的定标作用。因此,在大约一个世纪的时间内,确定移动星团视差的方法一直处于不断改进之中。对此作了简要介绍和评述。     

A multifrequency maximum-entropy joint analysis of COBE and Tenerife data

We use data from the Tenerife 10-, 15- and 33‐GHz beam-switching experiments along with the COBE 53- and 90‐GHz data to separate the cosmic microwave background (CMB) signal from the Galactic signal, and create two maps at high Galactic latitude. The new multi-MEM technique is used to obtain the best reconstruction of the two channels. The two maps are presented, and known features are identified within each. We find that the Galactic contribution to both the 15- and 33-GHz Tenerife data is small enough to be ignored when compared with the errors in the data and the magnitude of the CMB signal.     

Prospects for spectroscopic reflected-light planet searches

High-resolution spectroscopic searches for the starlight reflected from close-in extrasolar giant planets have the capability of determining the optical albedo spectra and scattering properties of these objects. When combined with radial velocity measurements they also yield the true mass of the planet. To date, only two such planets have been targeted for reflected-light signals, yielding upper limits on the optical albedos of the planets. Here we examine the prospects for future searches of this kind. We present Monte Carlo estimates of prior probability distributions for the orbital velocity amplitudes and planet/star flux ratios of six bright stars known to harbour giant planets in orbits with periods of less than 5 d. Using these estimates, we assess the viability of these targets for future reflected-light searches using 4- and 8-m class telescopes.     

Second‐order error propagation in the Mueller matrix of a spectropolarimeter

Accurate measurements of Stokes IQUV in spectral lines is required for precise reconstruction of stellar magnetic field geometries with Zeeman‐Dopper imaging. Spectral Zeeman features are intrinsically weak and subjected to a number of instrumental uncertainties. The aim of this work is to study the details of the instrumental uncertainties in the Stokes IQUV measurements in spectral lines and ways of their reduction. We make a practical comparison of the polarimetric performances of two high‐resolution échelle spectropolarimeters, namely SOFIN at the NOT, and HARPS at ESO. We show the residual spectra for both instruments to characterize the cross‐talk between the observed Stokes parameters. We employ a self‐calibrating least‐squares fit to eliminate some of the polarization uncertainties to derive the full Stokes vector from stellar spectra (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Cosmic‐ray rejection for single spectroscopic CCD images by means of template matching

This article describes a novel algorithm for cosmic‐ray rejection in single spectroscopic CCD images. This algorithm is based on a variation of template matching. It focuses on identifying those pixels belong to spectra, while other conventional algorithms tried to locate the cosmic‐ray hits directly. The main principle is applying template matching to find suspicious blocks, which is followed by surface patching to locate the legitimate pixels accurately. Therefore, the rest pixels are the ones corrupted by cosmic‐ray hits. Meanwhile, almost all the parameters are automatically extracted from the images. Examples of its performance are given for both simulated and observed images. It shows an advantage of significantly low false alarm rate with relatively high detection rate (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

GRBs on probation: Testing the UHE CR paradigm with IceCube

Gamma ray burst (GRB) fireballs provide one of very few astrophysical environments where one can contemplate the acceleration of cosmic rays to energies that exceed 10²⁰ eV. The assumption that GRBs are the sources of the observed cosmic rays generates a calculable flux of neutrinos produced when the protons interact with fireball photons. With data taken during construction IceCube has already reached a sensitivity to observe neutrinos produced in temporal coincidence with individual GRBs provided that they are the sources of the observed extra-galactic cosmic rays. We here point out that the GRB origin of cosmic rays is also challenged by the IceCube upper limit on a possible diffuse flux of cosmic neutrinos which should not be exceeded by the flux produced by all GRB over Hubble time. Our alternative approach has the advantage of directly relating the diffuse flux produced by all GRBs to measurements of the cosmic ray flux. It also generates both the neutrino flux produced by the sources and the associated cosmogenic neutrino flux in a synergetic way.     

Dynamical phasing of Type II Cepheids

In this paper we examine the problems of phasing using light curves and offer an alternate technique using the changes in acceleration to establish the zero point. We give astrophysical justification as to why this technique is useful and apply the technique to a selection of Type II Cepheids. We then examine some limitations of the technique which qualify its use.     

Magnetic-Tower Jet Solution for Launching Astrophysical Jets

In spite of the large number of global three-dimensional (3-D) magnetohydrodynamic (MHD) simulations of accretion disks and astrophysical jets, which have been developed since 2000, the launching mechanisms of jets is somewhat controversial. Previous studies of jets have concentrated on the effect of the large-scale magnetic fields permeating accretion disks. However, the existence of such global magnetic fields is not evident in various astrophysical objects, and their origin is not well understood. Thus, we study the effect of small-scale magnetic fields confined within the accretion disk. We review our recent findings on the formation of jets in dynamo-active accretion disks by using 3-D MHD simulations. In our simulations, we found the emergence of accumulated azimuthal magnetic fields from the inner region of the disk (the so-called magnetic tower) and also the formation of a jet accelerated by the magnetic pressure of the tower. Our results indicate that the magnetic tower jet is one of the most promising mechanisms for launching jets from the magnetized accretion disk in various astrophysical objects. We will discuss the formation of cosmic jets in the context of the magnetic tower model.     

Gamma-ray probe of cosmic ray pressure in galaxy clusters and cosmological implications

Cosmic rays produced in cluster accretion and merger shocks provide pressure to the intracluster medium (ICM) and affect the mass estimates of galaxy clusters. Although direct evidence for cosmic ray ions in the ICM is still lacking, they produce γ-ray emission through the decay of neutral pions produced in their collisions with ICM nucleons. We investigate the capability of the Gamma-ray Large Area Space Telescope ( GLAST ) and imaging atmospheric Čerenkov telescopes (IACTs) for constraining the cosmic ray pressure contribution to the ICM. We show that GLAST can be used to place stringent upper limits, a few per cent for individual nearby rich clusters, on the ratio of pressures of the cosmic rays and thermal gas. We further show that it is possible to place tight (≲10 per cent) constraints for distant  ( z ≲ 0.25)  clusters in the case of hard spectrum, by stacking signals from samples of known clusters. The GLAST limits could be made more precise with the constraint on the cosmic ray spectrum potentially provided by IACTs. Future γ-ray observations of clusters can constrain the evolution of cosmic ray energy density, which would have important implications for cosmological tests with upcoming X-ray and Sunyaev–Zel'dovich effect cluster surveys.