Detection of new point-sources in WMAP Cosmic Microwave Background (CMB) maps at high Galactic latitude. A new technique to extract point sources from CMB maps

In experimental microwave maps, point-sources can strongly affect the estimation of the power-spectrum and/or the test of Gaussianity of the Cosmic Microwave Background (CMB) component. As a consequence, their removal from the sky maps represents a critical step in the analysis of the CMB data. Before removing a source, however, it is necessary to detect it and source extraction consists of a delicate preliminary operation. In the literature, various techniques have been presented to detect point-sources in the sky maps. The most sophisticated ones exploit the multi-frequency nature of the observations that is typical of the CMB experiments. These techniques have "optimal" theoretical properties and, at least in principle, are capable of remarkable performances. Actually, they are rather difficult to use and this deteriorates the quality of the obtainable results. In this paper, we present a new technique, the "weighted matched filter" (WMF), that is quite simple to use and hence more robust in practical applications. Such technique shows particular efficiency in the detection of sources whose spectra have a slope different from zero. We apply this method to three Southern Hemisphere sky regions - each with an area of 400 square degrees - of the seven years Wilkinson Microwave Anisotropy Probe (WMAP) maps and compare the resulting sources with those of the two seven-year WMAP point-sources catalogues. In these selected regions we find seven additional sources not previously listed in WMAP catalogues and discuss their most likely identification and spectral properties.Comment: Astronomy and Astrophysics, 2011, in pres

An approach for the detection of point-sources in very high resolution microwave maps

This paper deals with the detection problem of extragalactic point-sources in multi-frequency, microwave sky maps that will be obtainable in future cosmic microwave background radiation (CMB) experiments with instruments capable of very high spatial resolution. With spatial resolutions that can be of order of 0.1-1.0 arcsec or better, the extragalactic point-sources will appear isolated. The same holds also for the compact structures due to the Sunyaev-Zeldovich (SZ) effect (both thermal and kinetic). This situation is different from the maps obtainable with instruments as WMAP or PLANCK where, because of the smaller spatial resolution (approximately 5-30 arcmin), the point-sources and the compact structures due to the SZ effect form a uniform noisy background (the "confusion noise"). Hence, the point-source detection techniques developed in the past are based on the assumption that all the emissions that contribute to the microwave background can be modeled with homogeneous and isotropic (often Gaussian) random fields and make use of the corresponding spatial power-spectra. In the case of very high resolution observations such an assumption cannot be adopted since it still holds only for the CMB. Here, we propose an approach based on the assumption that the diffuse emissions that contribute to the microwave background can be locally approximated by two-dimensional low order polynomials. In particular, two sets of numerical techniques are presented containing two different algorithms each. The performance of the algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/Submillimeter Array (ALMA).Comment: Accepted for publication on "Astronomy & Astrophysics". arXiv admin note: substantial text overlap with arXiv:1206.4536 Replaced version is the accepted one and published in A&

The Dependence of Galaxy Shape on Luminosity and Surface Brightness Profile

For a sample of 96,951 galaxies from the Sloan Digital Sky Survey Data Release 3, we study the distribution of apparent axis ratios as a function of r-band absolute magnitude and surface brightness profile type. We use the parameter fracDeV to quantify the profile type (fracDeV = 1 for a de Vaucouleurs profile; fracDeV = 0 for an exponential profile). When the apparent axis ratio q_{am} is estimated from the moments of the light distribution, the roundest galaxies are very bright (M_r \sim -23) de Vaucouleurs galaxies and the flattest are modestly bright (M_r \sim -18) exponential galaxies. When the apparent axis ratio q_{25} is estimated from the axis ratio of the 25 mag/arcsec^2 isophote, we find that de Vaucouleurs galaxies are flatter than exponential galaxies of the same absolute magnitude. For a given surface brightness profile type, very bright galaxies are rounder, on average, than fainter galaxies. We deconvolve the distributions of apparent axis ratios to find the distribution of the intrinsic short-to-long axis ratio gamma, assuming constant triaxiality T. For all profile types and luminosities, the distribution of apparent axis ratios is inconsistent with a population of oblate spheroids, but is usually consistent with a population of prolate spheroids. Bright galaxies with a de Vaucouleurs profile (M_r < -21.84, fracDeV > 0.9) have a distribution of q_{am} that is consistent with triaxiality in the range 0.4 < T < 0.8, with mean intrinsic axis ratio 0.66 < gamma < 0.69. The fainter de Vaucouleurs galaxies are best fit with prolate spheroids (T = 1) with mean axis ratio gamma = 0.51.Comment: 32 pages, 12 figures, to appear in Ap

From proper motions to star cluster dynamics: measuring velocity dispersion in deconvolved distribution functions

We investigate the effect that the usually large errors associated with ground-based proper motion (PM) components have on the determination of a star cluster's velocity dispersion (\sv). Rather than histograms, we work with PM distribution functions (PMDFs), taking the $1\sigma$ uncertainties formally into account. In this context, a cluster's intrinsic PMDF is broadened by the error distribution function (eDF) that, given the average error amplitude, has a width usually comparable to the cluster PMDF. Thus, we apply a Richardson-Lucy (RL) deconvolution to the PMDFs of a set of relatively nearby and populous open clusters (OCs), using the eDFs as point spread functions (PSFs). The OCs are NGC\,1039 (M\,34), NGC\,2477, NGC\,2516, NGC\,2682 (M\,67), and NGC\,7762. The deconvolved PMDFs are approximately Gaussian in shape, with dispersions lower than the observed ones by a factor of 4-10. NGC\,1039 and NGC\,2516, the nearest OCs of the sample, have deconvolved \sv\ compatible with those of bound OCs of mass $\sim10^3$\,\ms. NGC\,2477 and NGC\,2682 have deconvolved PMDFs with a secondary bump, shifted towards higher average velocities, which may be an artefact of the RL deconvolution when applied to asymmetric profiles. Alternatively, it may originate from cluster merger, large-scale mass segregation or, least probably, binaries.Comment: Accepted by MNRA

Oscillations of tori in the pseudo-Newtonian potential

Context. The high-frequency quasi-periodic oscillations (HF QPOs) in neutron star and stellar-mass black hole X-ray binaries may be the result of a resonance between the radial and vertical epicyclic oscillations in strong gravity. Aims. In this paper we investigate the resonant coupling between the epicyclic modes in a torus in a strong gravitational field. Methods. We perform numerical simulations of axisymmetric constant angular momentum tori in the pseudo-Newtonian potential. The epicyclic motion is excited by adding a constant radial velocity to the torus. Results. We verify that slender tori perform epicyclic motions at the frequencies of free particles, but the epicyclic frequencies decrease as the tori grow thicker. More importantly, and in contrast to previous numerical studies, we do not find a coupling between the radial and vertical epicyclic motions. The appearance of other modes than the radial epicyclic motion in our simulations is rather due to small numerical deviations from exact equilibrium in the initial state of our torus. Conclusions. We find that there is no pressure coupling between the two axisymmetric epicyclic modes as long as the torus is symmetric with respect to the equatorial plane. However we also find that there are other modes in the disc that may be more attractive for explaining the HF QPOs.Comment: 8 pages, 9 figure

One-Point Probability Distribution Functions of Supersonic Turbulent Flows in Self-Gravitating Media

Turbulence is essential for understanding the structure and dynamics of molecular clouds and star-forming regions. There is a need for adequate tools to describe and characterize the properties of turbulent flows. One-point probability distribution functions (pdf's) of dynamical variables have been suggested as appropriate statistical measures and applied to several observed molecular clouds. However, the interpretation of these data requires comparison with numerical simulations. To address this issue, SPH simulations of driven and decaying, supersonic, turbulent flows with and without self-gravity are presented. In addition, random Gaussian velocity fields are analyzed to estimate the influence of variance effects. To characterize the flow properties, the pdf's of the density, of the line-of-sight velocity centroids, and of the line centroid increments are studied. This is supplemented by a discussion of the dispersion and the kurtosis of the increment pdf's, as well as the spatial distribution of velocity increments for small spatial lags. From the comparison between different models of interstellar turbulence, it follows that the inclusion of self-gravity leads to better agreement with the observed pdf's in molecular clouds. The increment pdf's for small spatial lags become exponential for all considered velocities. However, all the processes considered here lead to non-Gaussian signatures, differences are only gradual, and the analyzed pdf's are in addition projection dependent. It appears therefore very difficult to distinguish between different physical processes on the basis of pdf's only, which limits their applicability for adequately characterizing interstellar turbulence.Comment: 38 pages (incl. 17 figures), accepted for publication in ApJ, also available with full resolution figures at http://www.strw.leidenuniv.nl/~klessen/Preprint

On the Nature of X-ray Variability in Ark 564

We use data from a recent long ASCA observation of the Narrow Line Seyfert 1 Ark 564 to investigate in detail its timing properties. We show that a thorough analysis of the time series, employing techniques not generally applied to AGN light curves, can provide useful information to characterize the engines of these powerful sources.We searched for signs of non-stationarity in the data, but did not find strong evidences for it. We find that the process causing the variability is very likely nonlinear, suggesting that variability models based on many active regions, as the shot noise model, may not be applicable to Ark 564. The complex light curve can be viewed, for a limited range of time scales, as a fractal object with non-trivial fractal dimension and statistical self-similarity. Finally, using a nonlinear statistic based on the scaling index as a tool to discriminate time series, we demonstrate that the high and low count rate states, which are indistinguishable on the basis of their autocorrelation, structure and probability density functions, are intrinsically different, with the high state characterized by higher complexity.Comment: 13 pages, 13 figures, accepted for publication in A&

Velocity Field Statistics in Star-Forming Regions. I. Centroid Velocity Observations

The probability density functions (pdfs) of molecular line centroid velocity fluctuations and fluctuation differences at different spatial lags are estimated for several nearby molecular clouds with active internal star formation. The data consist of over 75,000 $^{13}$CO line profiles divided among twelve spatially and/or kinematically distinct regions. Although three regions (all in Mon R2) appear nearly Gaussian, the others show strong evidence for non-Gaussian, often nearly exponential, centroid velocity pdfs, possibly with power law contributions in the far tails. Evidence for nearly exponential centroid pdfs in the neutral HI component of the ISM is also presented, based on older optical and radio observations. These results are in contrast to pdfs found in isotropic incompressible turbulence experiments and simulations. Furthermore, no evidence is found for the scaling of difference pdf kurtosis with Reynolds number which is seen in incompressible turbulence, and the spatial distribution of high-amplitude velocity differences shows little indication of the filamentary appearance predicted by decay simulations dominated by vortical interactions. The variation with lag of the difference pdf moments is presented as a constraint on future simulations.Comment: LaTeX, 23 pages, with 15 Figures included separately as gif image files. Refereed/revised version accepted to the Astrophysical Journal. A complete (but much larger) postscript version is available from http://ktaadn.gsfc.nasa.gov/~miesc

The Apparent and Intrinsic Shape of the APM Galaxy Clusters

We estimate the distribution of intrinsic shapes of APM galaxy clusters from the distribution of their apparent shapes. We measure the projected cluster ellipticities using two alternative methods. The first method is based on moments of the discrete galaxy distribution while the second is based on moments of the smoothed galaxy distribution. We study the performance of both methods using Monte Carlo cluster simulations covering the range of APM cluster distances and including a random distribution of background galaxies. We find that the first method suffers from severe systematic biases, whereas the second is more reliable. After excluding clusters dominated by substructure and quantifying the systematic biases in our estimated shape parameters, we recover a corrected distribution of projected ellipticities. We use the non-parametric kernel method to estimate the smooth apparent ellipticity distribution, and numerically invert a set of integral equations to recover the corresponding distribution of intrinsic ellipticities under the assumption that the clusters are either oblate or prolate spheroids. The prolate spheroidal model fits the APM cluster data best.Comment: 8 pages, including 7 figures, accepted for publication in MNRA