No Arabic abstract
We present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15hr and 1hr field data of the GBT observations previously presented in Masui et al (2013) and Switzer et al (2013), covering about 41 square degrees at 0.6<z<1.0, for which cross-correlations may be measured with the galaxy distribution of the WiggleZ Dark Energy Survey. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contamination using an independent component analysis technique (fastica), and develop a Fourier-based optimal estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission through the non-Gaussian nature of their probability distributions. The temperature power spectra of the intensity maps is dominated by instrumental noise on small scales which fastica, as a conservative subtraction technique of non-Gaussian signals, can not mitigate. However, we determine similar GBT-WiggleZ cross-correlation measurements to those obtained by the Singular Value Decomposition (SVD) method, and confirm that foreground subtraction with fastica is robust against 21cm signal loss, as seen by the converged amplitude of these cross-correlation measurements. We conclude that SVD and fastica are complementary methods to investigate the foregrounds and noise systematics present in intensity mapping datasets.
We introduce a new technique to constrain the line-of-sight integrated electron density of our Galactic halo $text{DM}_text{MW,halo}$ through analysis of the observed dispersion measure distributions of pulsars $text{DM}_text{pulsar}$ and fast radio bursts $text{DM}_text{FRB}$. We model these distributions, correcting for the Galactic interstellar medium, with kernel density estimation---well-suited to the small data regime---to find lower/upper bounds to the corrected $text{DM}_text{pulsar}$/$text{DM}_text{FRB}$ distributions: $max[text{DM}_text{pulsar}] approx 7pm2 text{ (stat)} pm 9 text{ (sys) pc cm}^{-3}$ and $min[text{DM}_text{FRB}] approx 63^{+27}_{-21} text{ (stat)} pm 9 text{ (sys) pc cm}^{-3}$. Using bootstrap resampling to estimate uncertainties, we set conservative limits on the Galactic halo dispersion measure $-2 < text{DM}_text{MW,halo} < 123 text{pc cm}^{-3}$ (95% c.l.). The upper limit is especially conservative because it may include a non-negligible contribution from the FRB host galaxies and a non-zero contribution from the cosmic web. It strongly disfavors models where the Galaxy has retained the majority of its baryons with a density profile tracking the presumed dark matter density profile. Last, we perform Monte Carlo simulations of larger FRB samples to validate our technique and assess the sensitivity of ongoing and future surveys. We recover bounds of several tens $text{pc cm}^{-3}$ which may be sufficient to test whether the Galaxy has retained a majority of its baryonic mass. We estimate that a sample of several thousand FRBs will significantly tighten constraints on $text{DM}_text{MW,halo}$ and offer a valuable complement to other analyses.
To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $sim$5~kpc to $sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)
We analyzed the distribution of the RC stars throughout Galactic bulge using 2MASS data. We mapped the position of the red clump in 1 sq.deg. size fields within the area |l|<=8.5deg and $3.5deg<=|b|<=8.5deg, for a total of 170 sq.deg. The red clump seen single in the central area splits into two components at high Galactic longitudes in both hemispheres, produced by two structures at different distances along the same line of sight. The X-shape is clearly visible in the Z-X plane for longitudes close to $l=0 deg axis. Crude measurements of the space densities of RC stars in the bright and faint RC populations are consistent with the adopted RC distances, providing further supporting evidence that the X-structure is real, and that there is approximate front-back symmetry in our bulge fields. We conclude that the Milky Way bulge has an X-shaped structure within $|l|<~2deg, seen almost edge on with respect to the line of sight. Additional deep NIR photometry extending into the innermost bulge regions combined with spectroscopic data is needed in order to discriminate among the different possibilities that can cause the observed X-shaped structure.
Traveltime tomography is a very effective tool to reconstruct acoustic, seismic or electromagnetic wave speed distribution. To infer the velocity image of the medium from the measurements of first arrivals is a typical example of ill-posed problem. In the framework of Tikhonov regularization theory, in order to replace an ill-posed problem by a well-posed one and to get a unique and stable solution, a stabilizing functional (stabilizer) has to be introduced. The stabilizer selects the desired solution from a class of solutions with a specific physical and/or geometrical property; e.g., the existence of sharp boundaries separating media with different petrophysical parameters. Usually stabilizers based on maximum smoothness criteria are used during the inversion process; in these cases the solutions provide smooth images which, in many situations, do not describe the examined objects properly. Recently a new algorithm of direct minimization of the Tikhonov parametric functional with minimum support stabilizer has been introduced; it produces clear and focused images of targets with sharp boundaries. In this research we apply this new technique to real radar tomographic data and we compare the obtained result with the solution generated by the more traditional minimum norm stabilizer.
Antenna layout is an important design consideration for radio interferometers because it determines the quality of the snapshot point spread function (PSF, or array beam). This is particularly true for experiments targeting the 21 cm Epoch of Reionization signal as the quality of the foreground subtraction depends directly on the spatial dynamic range and thus the smoothness of the baseline distribution. Nearly all sites have constraints on where antennas can be placed---even at the remote Australian location of the MWA (Murchison Widefield Array) there are rock outcrops, flood zones, heritages areas, emergency runways and trees. These exclusion areas can introduce spatial structure into the baseline distribution that enhance the PSF sidelobes and reduce the angular dynamic range. In this paper we present a new method of constrained antenna placement that reduces the spatial structure in the baseline distribution. This method not only outperforms random placement algorithms that avoid exclusion zones, but surprisingly outperforms random placement algorithms without constraints to provide what we believe are the smoothest constrained baseline distributions developed to date. We use our new algorithm to determine antenna placements for the originally planned MWA, and present the antenna locations, baseline distribution, and snapshot PSF for this array choice.