Source author record

Maciej Bilicki

Maciej Bilicki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2025arXiv

KiDS-Legacy: Cosmological constraints from cosmic shear with the complete Kilo-Degree Survey

We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-infrared, along with an extra $23$ square degrees of KiDS-like calibration observations of deep spectroscopic surveys. Improvements in our redshift distribution estimation methodology, combined with our enhanced calibration data and multi-band image simulations, allow us to extend our lensed sample out to a photometric redshift of $z_{\rm B}\leq2.0$. Compared to previous KiDS analyses, the increased survey area and redshift depth results in a $\sim32\%$ improvement in constraining power in terms of $Σ_8\equivσ_8\left(Ω_{\rm m}/0.3\right)^α= 0.821^{+0.014}_{-0.016}$, where $α= 0.58$ has been optimised to match the revised degeneracy direction of $σ_8$ and $Ω_{\rm m}$. We adopt a new physically motivated intrinsic alignment model that depends jointly on the galaxy sample's halo mass and spectral type distributions, and that is informed by previous direct alignment measurements. We also marginalise over our uncertainty on the impact of baryon feedback on the non-linear matter power spectrum. Comparing to previous KiDS analyses, we conclude that the increase seen in $S_8$ primarily results from our improved redshift distribution estimation and calibration, as well as new survey area and improved image reduction. Our companion paper Stölzner et al. (submitted) presents a full suite of internal and external consistency tests, finding the KiDS-Legacy data set to be the most internally robust sample produced by KiDS to date.

preprint2023arXiv

KiDS-1000 Cosmology: Constraints from density split statistics

Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = σ_8 \sqrt{Ω_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $Ω_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$.

preprint2022arXiv

COWS: A filament finder for Hessian cosmic web identifiers

The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells into an ensemble of individual discreet filaments. Specifically, a thinning algorithm is applied to velocity shear tensor based cosmic web (V-web) to identify the spine of the filaments. This results in a set of filaments with well defined end-point and length. It is confirmed that these sit at local density ridges and align with the appropriate direction defined by the underlying velocity field. The radial density profile of these curved cylindrical filaments, as well as the distribution of their lengths is also examined. The robustness of all results is checked against resolution and the V-web threshold. The code for the COWS method has been made publicly available.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Self-Organizing Map Application on Nearby Galaxies

Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a more complex population structure, represented by between 2 and 6 clusters. Here we use Self Organizing Maps (SOM), an unsupervised learning technique which can be used to visualize similarity in a higher dimensional space using a 2D representation, to map these five-dimensional clusters in the feature space onto two-dimensional projections. To further analyze these clusters, using the SOM information, we agree with previous results that the sub-populations found in the feature space can be reasonably mapped onto three or five clusters. We explore where the "green valley" galaxies are mapped onto the SOM, indicating multiple interstitial populations within the green valley population. Finally, we use the projection of the SOM to verify whether morphological information provided by GalaxyZoo users, for example, if features are visible, can be mapped onto the SOM-generated map. Voting on whether galaxies are smooth, likely ellipticals, or "featured" can reasonably be separated but smaller morphological features (bar, spiral arms) can not. SOMs promise to be a useful tool to map and identify instructive sub-populations in multidimensional galaxy survey feature space, provided they are large enough.

preprint2022arXiv

Joint constraints on cosmology and the impact of baryon feedback: combining KiDS-1000 lensing with the thermal Sunyaev-Zeldovich effect from Planck and ACT

We conduct a pseudo-$C_\ell$ analysis of the tomographic cross-correlation between 1000 deg$^2$ of weak lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev-Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMx, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive constraints on both cosmology and baryon feedback for the first time from these data, marginalising over redshift uncertainties, intrinsic alignment of galaxies, and contamination by the cosmic infrared background (CIB). We find our results to be insensitive to the CIB, while intrinsic alignment provides a small but significant contribution to the lensing--tSZ cross-correlation. The cosmological constraints are consistent with those of other low-redshift probes and prefer strong baryon feedback. The inferred amplitude of the lensing--tSZ cross-correlation signal, which scales as $σ_8(Ω_\mathrm{m}/0.3)^{0.2}$, is low by $\sim 2\,σ$ compared to the primary cosmic microwave background constraints by Planck. The lensing--tSZ measurements are then combined with pseudo-$C_\ell$ measurements of KiDS-1000 cosmic shear into a novel joint analysis, accounting for the full cross-covariance between the probes, providing tight cosmological constraints by breaking parameter degeneracies inherent to both probes. The joint analysis gives an improvement of 40% on the constraint of $S_8=σ_8\sqrt{Ω_\mathrm{m}/0.3}$ over cosmic shear alone, while providing constraints on baryon feedback consistent with hydrodynamical simulations, demonstrating the potential of such joint analyses with baryonic tracers such as the tSZ effect. We discuss remaining modelling challenges that need to be addressed if these baryonic probes are to be included in future precision-cosmology analyses.

preprint2022arXiv

The dark matter halo masses of elliptical galaxies as a function of observationally robust quantities

Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2409 central massive elliptical galaxies from the SDSS spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sersic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame $r-$band luminosity, half-light radius and rest-frame $u-g$ colour, using KiDS weak lensing data and a Bayesian hierarchical approach. For the sake of robustness to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing total luminosity and half-light radius with the luminosity within a 10~kpc aperture, $L_{r,10}$, and the light-weighted surface brightness slope, $Γ_{10}$. Results. We did not detect any correlation between halo mass and either half-light radius or colour, at fixed redshift and luminosity. Conclusions. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than halo mass.

preprint2022arXiv

Universality of the halo mass function in modified gravity cosmologies

We study the halo mass function (HMF) in modified gravity (MG) models using a set of large $N$-body simulations -- the ELEPHANT suite. We consider two popular beyond-general relativity scenarios: the Hu-Sawicki chameleon $f(R)$ model and the normal branch of the Dvali-Gabadadze-Porrati braneworld (nDGP). We show that in MG, analytic formulation based on the Press-Schechter framework offers a grossly inaccurate description of the HMF. We find, however, that once the HMF is expressed in terms of the dimensionless multiplicity function, it approximately assumes a redshift-independent universal character for all the models. Exploiting this property, we propose universal fits for the MG HMF in terms of their fractional departures from the $\mathrm{Λ\text{CDM}}$ case. We find two enclosed formulas, one for $f(R)$ and another for nDGP, that provide a reliable description of the HMF over the mass range covered by the simulations. These are accurate to a few percent with respect to the $N$-body data. We test the extrapolation potential of our fits against separate simulations with a different cosmological background and mass resolution, and find very good accuracy, within $\sim 10\%$. A particularly interesting finding from our analysis is a Gaussian-like shape of the HMF deviation that seems to appear universally across the whole $f(R)$ family, peaking at a mass variance scale characteristic for each $f(R)$ variant. We attribute this behavior to the specific physics of the environmentally-dependent chameleon screening models.

preprint2020arXiv

KiDS-1000 catalogue: Weak gravitational lensing shear measurements

We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our `gold-sample' of galaxies, with well-calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of $6.17$ galaxies per square arcminute. We quantify the accuracy of the spatial, temporal, and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a $0.1σ$ change in the inferred cosmic shear constraints on the clustering cosmological parameter $S_8 = σ_8\sqrt{Ω_{\rm m}/0.3}$. Through a series of two-point null-tests, we validate the shear estimates, finding no evidence for significant non-lensing B-mode distortions in the data. The PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is, however, shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test, we verify the expected redshift scaling of the galaxy-galaxy lensing signal around luminous red galaxies. We conclude that the joint KiDS-1000 shear and photometric redshift calibration is sufficiently robust for combined-probe gravitational lensing and spectroscopic clustering analyses.

preprint2020arXiv

KiDS-1000 Cosmology: Cosmic shear constraints and comparison between two point statistics

We present cosmological constraints from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000), doubling the survey area with nine-band optical and near-infrared photometry with respect to previous KiDS analyses. Adopting a spatially flat $Λ$CDM model, we find $S_8 = σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.759^{+0.024}_{-0.021}$ for our fiducial analysis, which is in $3σ$ tension with the prediction of the Planck Legacy analysis of the cosmic microwave background. We compare our fiducial COSEBIs (Complete Orthogonal Sets of E/B-Integrals) analysis with complementary analyses of the two-point shear correlation function and band power spectra, finding results to be in excellent agreement. We investigate the sensitivity of all three statistics to a number of measurement, astrophysical, and modelling systematics, finding our $S_8$ constraints to be robust and dominated by statistical errors. Our cosmological analysis of different divisions of the data pass the Bayesian internal consistency tests, with the exception of the second tomographic bin. As this bin encompasses low redshift galaxies, carrying insignificant levels of cosmological information, we find that our results are unchanged by the inclusion or exclusion of this sample.

preprint2020arXiv

KiDS-1000 Cosmology: Multi-probe weak gravitational lensing and spectroscopic galaxy clustering constraints

We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of large-scale structure probes breaks the degeneracies between cosmological parameters for individual observables, resulting in a constraint on the structure growth parameter $S_8=σ_8 \sqrt{Ω_{\rm m}/0.3} = 0.766^{+0.020}_{-0.014}$, that has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck. The recovered $S_8$ amplitude is low, however, by $8.3 \pm 2.6$ % relative to Planck. This result builds from a series of KiDS-1000 analyses where we validate our methodology with variable depth mock galaxy surveys, our lensing calibration with image simulations and null-tests, and our optical-to-near-infrared redshift calibration with multi-band mock catalogues and a spectroscopic-photometric clustering analysis. The systematic uncertainties identified by these analyses are folded through as nuisance parameters in our cosmological analysis. Inspecting the offset between the marginalised posterior distributions, we find that the $S_8$-difference with Planck is driven by a tension in the matter fluctuation amplitude parameter, $σ_8$. We quantify the level of agreement between the CMB and our large-scale structure constraints using a series of different metrics, finding differences with a significance ranging between $\sim\! 3\,σ$, when considering the offset in $S_{8}$, and $\sim\! 2\,σ$, when considering the full multi-dimensional parameter space.

preprint2020arXiv

Tomographic measurement of the intergalactic gas pressure through galaxy-tSZ cross-correlations

We cross-correlate maps of the thermal Sunyaev-Zeldovich (tSZ) Compton-$y$ parameter published by Planck with the projected distribution of galaxies in a set of low-redshift tomographic bins. We use the nearly full-sky 2MASS Photometric Redshift and WISE $\times$ SuperCOSMOS public catalogues, covering the redshift range $z\lesssim0.4$. Our measurements allow us to place constraints on the redshift dependence of the mass-observable relation for tSZ cluster count analyses in terms of the so-called 'hydrostatic mass bias' parameter $1-b_{\rm H}$. These results can also be interpreted as measurements of the bias-weighted average gas pressure $\langle bP_e\rangle$ as a function of redshift, a quantity that can be related to the thermodynamics of gas inside haloes and used to constrain energy injection processes. We measure $1-b_{\rm H}$ with $\sim6\%$ precision in 6 equispaced redshift bins, and find no evidence for a redshift-dependent mass bias parameter, in agreement with previous analyses. Our mean value of $1-b_{\rm H} = 0.75\pm0.03$ is also in good agreement with the one estimated by the joint analysis of Planck cluster counts and CMB anisotropies calibrated with CMB lensing. Our measurements of $\langle bP_e\rangle$, at the level of $\sim10\%$ in each bin, are the most stringent constraints on the redshift dependence of this parameter to date, and agree well both with previous measurements and with theoretical expectations from shock-heating models.

preprint2019arXiv

Galaxy classification: A machine learning analysis of GAMA catalogue data

We present a machine learning analysis of five labelled galaxy catalogues from the Galaxy And Mass Assembly (GAMA): The SersicCatVIKING and SersicCatUKIDSS catalogues containing morphological features, the GaussFitSimple catalogue containing spectroscopic features, the MagPhys catalogue including physical parameters for galaxies, and the Lambdar catalogue, which contains photometric measurements. Extending work previously presented at the ESANN 2018 conference - in an analysis based on Generalized Relevance Matrix Learning Vector Quantization and Random Forests - we find that neither the data from the individual catalogues nor a combined dataset based on all 5 catalogues fully supports the visual-inspection-based galaxy classification scheme employed to categorise the galaxies. In particular, only one class, the Little Blue Spheroids, is consistently separable from the other classes. To aid further insight into the nature of the employed visual-based classification scheme with respect to physical and morphological features, we present the galaxy parameters that are discriminative for the achieved class distinctions.

preprint2016arXiv

Dependence of GAMA galaxy halo masses on the cosmic web environment from 100 square degrees of KiDS weak lensing data

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

preprint2016arXiv

Discovery of a supercluster in the ZOA in Vela

We report the discovery of a potentially major supercluster that extends across the Galactic Plane in the constellation of Vela, at a mean recessional velocity of ~18,000 km/s. Recent multi-object spectroscopic observations of this Vela Supercluster (VSCL), using AAOmega+2dF and the Southern African Large Telescope, confirm an extended galaxy overdensity in the Zone of Avoidance (ZOA) located where residual bulk flows predict a considerable mass excess. We present a preliminary analysis of ~4,500 new spectroscopic galaxy redshifts obtained in the ZOA centred on the Vela region (l=272.5+-20 deg, b=0+-10 deg). The presently sparsely-sampled dataset traces an overdensity that covers 25 deg in Galactic longitude on either side of the Plane, suggesting an extent of 25 deg x 20 deg, corresponding to ~115 x 90 $h_{70}$ Mpc at the supercluster redshift. In redshift space, the overdensity appears to consist of two merging wall-like structures, interspersed with clusters and groups. Both the velocity histogram and the morphology of the multi-branching wall structure are consistent with a supercluster classification. $K_s^o$ galaxy counts show an enhancement of ~1.2 over the survey area for galaxies brighter than $M_K^*$ at the VSCL distance, and a galaxy overdensity of $δ=0.50\rm{-}0.77$ within a photometric redshift shell around the VSCL, when compared to various Two-Micron All-Sky Survey samples. Taking account of selection effects, the VSCL is estimated to contribute $v_\rm{LG} \gtrsim 50$ km/s to the motion of the Local Group.

preprint2016arXiv

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

preprint2016arXiv

Towards automatic classification of all WISE sources

The WISE satellite has detected hundreds of millions sources over the entire sky. Classifying them reliably is however a challenging task due to degeneracies in WISE multicolour space and low levels of detection in its two longest-wavelength bandpasses. Here we aim at obtaining comprehensive and reliable star, galaxy and quasar catalogues based on automatic source classification in full-sky WISE data. This means that the final classification will employ only parameters available from WISE itself, in particular those reliably measured for a majority of sources. For the automatic classification we applied the support vector machines (SVM) algorithm, which requires a training sample with relevant classes already identified, and we chose to use the SDSS spectroscopic dataset for that purpose. By calibrating the classifier on the test data drawn from SDSS, we first established that a polynomial kernel is preferred over a radial one for this particular dataset. Next, using three classification parameters (W1 magnitude, W1-W2 colour, and a differential aperture magnitude) we obtained very good classification efficiency in all the tests. At the bright end, the completeness for stars and galaxies reaches ~95%, deteriorating to ~80% at W1=16 mag, while for quasars it stays at a level of ~95% independently of magnitude. Similar numbers are obtained for purity. Application of the classifier to full-sky WISE data, flux-limited to 16 mag (Vega) in the 3.4 μm channel, and appropriate a posteriori cleaning allowed us to obtain reliably-looking catalogues of star and galaxy candidates. However, the sources flagged by the classifier as `quasars' are in fact dominated by dusty galaxies but also exhibit contamination from sources located mainly at low ecliptic latitudes, consistent with Solar System objects. [abridged]

preprint2015arXiv

Automatised classification of WISE sources: first results, future prospects

We present the first results of our dedicated programme of automatised classification of galaxies, stars and quasars in the mid-infrared all-sky data from the WISE survey. We employ the Support Vector Machines (SVM) algorithm, which defines a hyperplane separating different classes of sources in a multidimensional space of arbitrarily chosen parameters. This approach consists of four general steps: 1) selection of the training sample, 2) selection of the optimal parameter space, 3) training of the classifier, 4) application to target data. Here, as the training set, we use sources from a cross-correlation of the WISE catalogue with the SDSS spectroscopic sample. The performance of the SVM classifier was tested as a function of size of the training set, dimension of the parameter space, WISE apparent magnitude and Galactic extinction. We find that our classifier provides promising results already for three classification parameters: magnitude, colour and differential aperture magnitude. Completeness and purity levels as high as 95% are obtained for quasars, while for galaxies and stars they vary between 80-95% depending on the magnitude, deteriorating for fainter sources.

preprint2015arXiv

Cosmology with all-sky surveys

Various aspects of cosmology require comprehensive all-sky mapping of the cosmic web to considerable depths. In order to probe the whole extragalactic sky beyond 100 Mpc, one must draw on multiwavelength datasets and state-of-the-art photometric redshift techniques. Here I summarize our dedicated program that employs the largest photometric all-sky surveys -- 2MASS, WISE and SuperCOSMOS -- to obtain accurate redshift estimates of millions of galaxies. The first outcome of these efforts -- the 2MASS Photometric Redshift catalog (2MPZ) -- was publicly released in 2013 and includes almost 1 million galaxies with a median redshift of z~0.1. I discuss how this catalog was constructed and how it is being used for various cosmological tests. I also present how combining the WISE mid-infrared survey with SuperCOSMOS optical data allowed us to push to depths over 1 Gpc on unprecedented angular scales. These photometric redshift samples, with about 20 million sources in total, provide access to volumes large enough to study observationally the Copernican Principle of universal homogeneity and isotropy, as well as to probe various aspects of dark energy and dark matter through cross-correlations with other data such as the cosmic microwave or gamma-ray backgrounds. Last but not least, they constitute a test-bed for forthcoming wide-angle multi-million galaxy samples expected from such instruments as the SKA, Euclid, or LSST.

preprint2015arXiv

Homogeneity and isotropy in the 2MASS Photometric Redshift catalogue

Using the 2MASS Photometric Redshift catalogue we perform a number of statistical tests aimed at detecting possible departures from statistical homogeneity and isotropy in the large-scale structure of the Universe. Making use of the angular homogeneity index, an observable proposed in a previous publication, as well as studying the scaling of the angular clustering and number counts with magnitude limit, we place constraints on the fractal nature of the galaxy distribution. We find that the statistical properties of our sample are in excellent agreement with the standard cosmological model, and that it reaches the homogeneous regime significantly faster than a class of fractal models with dimensions $D<2.75$. As part of our search for systematic effects, we also study the presence of hemispherical asymmetries in our data, finding no significant deviation beyond those allowed by the concordance model.

preprint2015arXiv

Learning algorithms at the service of WISE survey

We have undertaken a dedicated program of automatic source classification in the WISE database merged with SuperCOSMOS scans, comprehensively identifying galaxies, quasars and stars on most of the unconfused sky. We use the Support Vector Machines classifier for that purpose, trained on SDSS spectroscopic data. The classification has been applied to a photometric dataset based on all-sky WISE 3.4 and 4.6 $μ$m information cross-matched with SuperCOSMOS B and R bands, which provides a reliable sample of $\sim170$ million sources, including galaxies at $z_{\rm med}\sim0.2$, as well as quasars and stars. The results of our classification method show very high purity and completeness (more than 96\%) of the separated sources, and the resultant catalogs can be used for sophisticated analyses, such as generating all-sky photometric redshifts.

preprint2015arXiv

Results from the Wide-field Infrared Survey Explorer (WISE) Future Uses Session at the WISE at 5 Meeting

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

preprint2014arXiv

Mapping the Cosmic Web with the largest all-sky surveys

Our view of the low-redshift Cosmic Web has been revolutionized by galaxy redshift surveys such as 6dFGS, SDSS and 2MRS. However, the trade-off between depth and angular coverage limits a systematic three-dimensional account of the entire sky beyond the Local Volume (z<0.05). In order to reliably map the Universe to cosmologically significant depths over the full celestial sphere, one must draw on multiwavelength datasets and state-of-the-art photometric redshift techniques. We have undertaken a dedicated program of cross-matching the largest photometric all-sky surveys -- 2MASS, WISE and SuperCOSMOS -- to obtain accurate redshift estimates of millions of galaxies. The first outcome of these efforts -- the 2MASS Photometric Redshift catalog (2MPZ, Bilicki et al. 2014a) -- has been publicly released and includes almost 1 million galaxies with a mean redshift of z=0.08. Here we summarize how this catalog was constructed and how using the WISE mid-infrared sample together with SuperCOSMOS optical data allows us to push to redshift shells of z~0.2--0.3 on unprecedented angular scales. Our catalogs, with ~20 million sources in total, provide access to cosmological volumes crucial for studies of local galaxy flows (clustering dipole, bulk flow) and cross-correlations with the cosmic microwave background such as the integrated Sachs-Wolfe effect or lensing studies.

preprint2013arXiv

2MASS Photometric Redshift catalog: a comprehensive three-dimensional census of the whole sky

Key cosmological applications require the three-dimensional galaxy distribution on the entire celestial sphere. These include measuring the gravitational pull on the Local Group, estimating the large-scale bulk flow and testing the Copernican principle. However, the largest all-sky redshift surveys -- the 2MRS and IRAS PSCz -- have median redshifts of only z=0.03 and sample the very local Universe. There exist all-sky galaxy catalogs reaching much deeper -- SuperCOSMOS in the optical, 2MASS in the near-IR and WISE in the mid-IR -- but these lack complete redshift information. At present, the only rapid way towards larger 3D catalogs covering the whole sky is through photometric redshift techniques. In this paper we present the 2MASS Photometric Redshift catalog (2MPZ) containing 1 million galaxies, constructed by cross-matching 2MASS XSC, WISE and SuperCOSMOS all-sky samples and employing the artificial neural network approach (the ANNz algorithm), trained on such redshift surveys as SDSS, 6dFGS and 2dFGRS. The derived photometric redshifts have errors nearly independent of distance, with an all-sky accuracy of sigma_z=0.015 and a very small percentage of outliers. In this way, we obtain redshift estimates with a typical precision of 12% for all the 2MASS XSC galaxies that lack spectroscopy. In addition, we have made an early effort towards probing the entire 3D sky beyond 2MASS, by pairing up WISE with SuperCOSMOS and training the ANNz on GAMA redshift data reaching currently to z_med~0.2. This has yielded photo-z accuracies comparable to those in the 2MPZ. These all-sky photo-z catalogs, with a median z~0.1 for the 2MPZ, and significantly deeper for future WISE-based samples, will be the largest and most complete of their kind for the foreseeable future.

preprint2013arXiv

Hierarchical clustering and the BAO signature

In this contribution we present the preliminary results regarding the non-linear BAO signal in higher-order statistics of the cosmic density field. We use ensembles of N-body simulations to show that the non-linear evolution changes the amplitudes of the BAO signal, but has a negligible effect on the scale of the BAO feature. The latter observation accompanied by the fact that the BAO feature amplitude roughly doubles as one moves to higher orders, suggests that the higher-order correlation amplitudes can be used as probe of the BAO signal.

preprint2013arXiv

Precision cosmology with the 2MASS clustering dipole

Comparison of peculiar velocities of galaxies with their gravitational accelerations (induced by the density field) is one of the methods to constrain the redshift distortion parameter β=(Ω_m^0.55)/b, where Ω_m is the non-relativistic matter density parameter and b is the linear bias. In particular, one can use the motion of the Local Group (LG) for that purpose. Its peculiar velocity is known from the dipole component of the cosmic microwave background, whereas its acceleration can be estimated with the use of an all-sky galaxy catalog, from the so-called clustering dipole. At the moment, the biggest dataset of that kind is the Two Micron All Sky Survey Extended Source Catalog (2MASS XSC) containing almost 1 million galaxies and complete up to ~300 Mpc/h. We applied 2MASS data to measure LG acceleration and used two methods to estimate the beta parameter. Both of them yield β~0.4 with an error of several per cent, which is the most precise determination of this parameter from the clustering dipole to date.

preprint2012arXiv

Motion of the Local Group as a cosmological probe

In this thesis, we use the motion of the Local Group of galaxies (LG) through the Universe to measure the cosmological parameter of non-relativistic matter density, Omega_m. For that purpose, we compare the peculiar velocity of the LG with its gravitational acceleration. The former is known from the dipole of the cosmic microwave background radiation and the latter is estimated here from the clustering dipole of galaxies in the Two Micron All Sky Survey (2MASS) Extended Source Catalog. We start by presenting the general framework of perturbation theory of gravitational instability in the expanding Universe and how it applies to the peculiar motion of the LG. Next, we study a particular effect for the dipole measurement, related to the fact that a nearby Local Void is partially hidden behind our Galaxy. We then describe in detail how we handled the 2MASS extragalactic data for the purpose of our analysis. Finally, we present two methods to estimate the density Omega_m, combined with the linear biasing b into the parameter beta = (Omega_m)^{0.55} / b, from the comparison of the LG velocity and acceleration. The first approach is to study the growth of the 2MASS clustering dipole with increased depth of the sample and compare it with theoretical expectations. The second is to apply the maximum-likelihood method in order to improve the precision of the measurement. With both these methods we find beta=0.4 and Omega_m=0.2, which is consistent with various independent estimates. We also briefly mention some future prospects in the field.

preprint2012arXiv

We do not live in the R_h = c t universe

We analyse the possibility that our Universe could be described by the model recently proposed by Melia & Shevchuk (2012), where the Hubble scale R_h=c/H is at all times equal to the distance ct that light has travelled since the Big Bang. In such a model, the scale factor is proportional to cosmic time and there is neither acceleration nor deceleration of the expansion. We first point out problems with the very foundations of the model and its consequences for the evolution of the Universe. Next, we compare predictions of the model with observational data. As probes of the expansion we use distance data of supernovae type Ia, as well as Hubble rate data obtained from cosmic chronometers and radial baryon acoustic oscillations. We analyse the redshift evolution of the Hubble parameter and its redshift derivatives, together with the so-called O_m diagnostic and the deceleration parameter. To reliably estimate smooth functions and their derivatives from discrete data, we use the recently developed Gaussian Processes in Python package (GaPP). Our general conclusion is that the discussed model is strongly disfavoured by observations, especially at low redshifts (z<0.5). In particular, it predicts specific constant values for the deceleration parameter and for redshift derivatives of the Hubble parameter, which is ruled out by the data.

preprint2011arXiv

Is the 2MASS clustering dipole convergent?

There is a long-standing controversy about the convergence of the dipole moment of the galaxy angular distribution (the so-called clustering dipole). We study the growth of the clustering dipole of galaxies as a function of the limiting flux of the sample from the Two Micron All Sky Survey (2MASS). Contrary to some earlier claims, we find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to 13.5 mag in the near-infrared K_s band (equivalent to an effective distance of 300 Mpc/h). We compare the observed growth of the dipole with the theoretically expected, conditional one (i.e., given the velocity of the Local Group relative to the CMB), for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of its theoretical counterpart once the proper observational window of the 2MASS flux-limited catalog is included. For a contrast, if the adopted window is a top-hat, then the predicted dipole grows significantly faster and converges to its final value for a distance of about 300 Mpc/h. By comparing the observational windows, we show that for a given flux limit and a corresponding distance limit, the 2MASS flux-weighted window passes less large-scale signal than the top-hat one. We conclude that the growth of the 2MASS dipole for effective distances greater than 200 Mpc/h is only apparent. On the other hand, for a distance of 80 Mpc/h (mean depth of the 2MASS Redshift Survey) and the LambdaCDM power spectrum, the true dipole is expected to reach only ~80% of its final value. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.04, which leads to an estimate of the density parameter Omega_m = 0.20+-0.08.

preprint2010arXiv

Influence of the Local Void on measurements of the clustering dipole

In measurements of the clustering dipole from all-sky surveys, an important problem is the lack of information about galaxy distribution in the so-called Zone of Avoidance (ZoA). The existence of the Local Void (LV) has a systematic effect on these measurements. If the ZoA is randomly filled with mock galaxies, then the calculated acceleration of the Local Group of galaxies (LG) has a spurious component, resulting from the lack of real galaxies in the intersection of the LV with the ZoA. This component affects both the misalignment angle between the clustering dipole and the CMB dipole, and the inferred value of mean matter density Omega_m. We calculate the amplitude of the spurious acceleration acting on the LG due to the LV. Its value depends on the geometry and size of the LV, as well as on its density contrast. However, under the simplest assumption of the LV being spherical and completely empty, within the linear theory this amplitude amounts only to about 45 km/s in units of velocity. The resulting change in the misalignment angle is smaller than 1 degree, and the fractional change in the deduced value of Omega_m is about 5%. Accounting for observationally indicated elongation of the LV and maintaining the maximising assumption of a complete lack of galaxies inside increases these numbers only moderately. Specifically, the amplitude of the spurious acceleration rises to about 60 km/s, the misalignment angle remains still smaller than 1 degree, and the fractional change in the deduced value of Omega_m is enhanced to about 7%. Thus, despite the overall importance of the Local Void for the motion of the Local Group, the influence of the intersection of the LV with the ZoA on measurements of the clustering dipole is found to be only a minor systematic effect.

preprint2010arXiv

Is the 2MASS dipole convergent?

We study the growth of the clustering dipole of galaxies from the Two Micron All Sky Survey (2MASS). We find that the dipole does not converge before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to about 300 Mpc/h. We compare the observed growth of the dipole with the theoretically expected, conditional growth for the LambdaCDM power spectrum and cosmological parameters constrained by WMAP. The observed growth turns out to be within 1-sigma confidence level of the theoretical one, once the proper observational window of the 2MASS flux dipole is included. For a contrast, if the adopted window is a top hat, then the predicted dipole grows significantly faster and converges to its final value at a distance of about 200 Mpc/h. We study the difference between the top-hat window and the window for the flux-limited 2MASS survey and we conclude that the growth of the 2MASS dipole at effective distances greater than 200 Mpc/h is only apparent. Eventually, since for the window function of 2MASS the predicted growth is consistent with the observed one, we can compare the two to evaluate beta = (Omega_m)^0.55 / b. The result is beta = 0.38+-0.05, which gives a rough estimate of Omega_m = 0.2+-0.1.

preprint2008arXiv

The velocity-density relation in the spherical model

We study the cosmic velocity-density relation using the spherical collapse model (SCM) as a proxy to non-linear dynamics. Although the dependence of this relation on cosmological parameters is known to be weak, we retain the density parameter Omega_m in SCM equations, in order to study the limit Omega_m -> 0. We show that in this regime the considered relation is strictly linear, for arbitrary values of the density contrast, on the contrary to some claims in the literature. On the other hand, we confirm that for realistic values of Omega_m the exact relation in the SCM is well approximated by the classic formula of Bernardeau (1992), both for voids (delta<0) and for overdensities up to delta ~ 3. Inspired by this fact, we find further analytic approximations to the relation for the whole range delta from -1 to infinity. Our formula for voids accounts for the weak Omega_m-dependence of their maximal rate of expansion, which for Omega_m < 1 is slightly smaller that 3/2. For positive density contrasts, we find a simple relation div v = 3 H_0 (Omega_m)^(0.6) [ (1+delta)^(1/6) - (1+delta)^(1/2) ], that works very well up to the turn-around (i.e. up to delta ~ 13.5 for Omega_m = 0.25 and neglected Omega_Lambda). Having the same second-order expansion as the formula of Bernardeau, it can be regarded as an extension of the latter for higher density contrasts. Moreover, it gives a better fit to results of cosmological numerical simulations.