Search

The Euclid mission, designed to map the geometry of the dark Universe, presents an unprecedented opportunity for advancing our understanding of the cosmos through its photometric galaxy cluster survey. This paper focuses on enhancing the precision of halo bias (HB) predictions, which is crucial for deriving cosmological constraints from the clustering of galaxy clusters. Our study is based on the peak-background split (PBS) model linked to the halo mass function (HMF); it extends with a parametric correction to precisely align with results from an extended set of $N$-body simulations carried out with the OpenGADGET3 code. Employing simulations with fixed and paired initial conditions, we meticulously analyze the matter-halo cross-spectrum and model its covariance using a large number of mock catalogs generated with Lagrangian Perturbation Theory simulations with the PINOCCHIO code. This ensures a comprehensive understanding of the uncertainties in our HB calibration. Our findings indicate that the calibrated HB model is remarkably resilient against changes in cosmological parameters including those involving massive neutrinos. The robustness and adaptability of our calibrated HB model provide an important contribution to the cosmological exploitation of the cluster surveys to be provided by the Euclid mission. This study highlights the necessity of continuously refining the calibration of cosmological tools like the HB to match the advancing quality of observational data. As we project the impact of our model on cosmological constraints, we find that, given the sensitivity of the Euclid survey, a miscalibration of the HB could introduce biases in cluster cosmology analyses. Our work fills this critical gap, ensuring the HB calibration matches the expected precision of the Euclid survey. The implementation of our model is publicly available in https://github.com/TiagoBsCastro/CCToolkit., Comment: 20 pages; 12 figures; accepted for publication in A&A; abstract abridged for arXiv submission

We used the Javalambre Photometric Local Universe Survey (J-PLUS) DR2 photometry in twelve optical bands over 2176 deg2 to estimate the fraction of white dwarfs with presence of CaII H+K absorption along the cooling sequence. We compared the J-PLUS photometry against metal-free theoretical models to estimate the equivalent width in the J0395 passband of 10 nm centered at 395 nm (EW_J0395), a proxy to detect calcium absorption. A total of 4399 white dwarfs within 30000 > Teff > 5500 K and mass M > 0.45 Msun were analyzed. Their EW_J0395 distribution was modeled using two populations, corresponding to polluted and non-polluted systems, to estimate the fraction of calcium white dwarfs (f_Ca) as a function of Teff. The probability for each individual white dwarf of presenting calcium absorption, pca, was also computed. The comparison with both the measured Ca/He abundance and the metal pollution from spectroscopy shows that EW_J0395 correlates with the presence of calcium. The fraction of calcium white dwarfs increases from f_Ca = 0 at Teff = 13500 K to f_Ca = 0.15 at Teff = 5500 K. We compare our results with the fractions derived from the 40 pc spectroscopic sample and from SDSS spectra. The trend found in J-PLUS observations is also present in the 40 pc sample, however SDSS shows a deficit of metal-polluted objects at Teff < 12000 K. Finally, we found 39 white dwarfs with pca > 0.99. Twenty of them have spectra presented in previous studies, whereas we observed six additional targets. These 26 objects were all confirmed as metal-polluted systems. The J-PLUS optical data provide a robust statistical measurement for the presence of CaII H+K absorption in white dwarfs. We find a 15 +- 3 % increase in the fraction of calcium white dwarfs from Teff = 13500 K to 5500 K, which reflects their selection function in the optical from the total population of metal-polluted systems., Comment: Submitted to Astronomy and Astrophysics. Comments are welcome

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

With its 12 optical filters, the Javalambre-Photometric Local Universe Survey (J-PLUS) provides an unprecedented multicolor view of the local Universe. The third data release (DR3) covers 3,192 deg$^2$ and contains 47.4 million objects. However, the classification algorithms currently implemented in its pipeline are deterministic and based solely on the sources morphology. Our goal is classify the sources identified in the J-PLUS DR3 images into stars, quasi-stellar objects (QSOs), and galaxies. For this task, we present BANNJOS, a machine learning pipeline that uses Bayesian neural networks to provide the probability distribution function (PDF) of the classification. BANNJOS is trained on photometric, astrometric, and morphological data from J-PLUS DR3, Gaia DR3, and CatWISE2020, using over 1.2 million objects with spectroscopic classification from SDSS DR18, LAMOST DR9, DESI EDR, and Gaia DR3. Results are validated using $1.4 10^5$ objects and cross-checked against theoretical model predictions. BANNJOS outperforms all previous classifiers in terms of accuracy, precision, and completeness across the entire magnitude range. It delivers over 95% accuracy for objects brighter than $r = 21.5$ mag, and ~90% accuracy for those up to $r = 22$ mag, where J-PLUS completeness is < 25%. BANNJOS is also the first object classifier to provide the full probability distribution function (PDF) of the classification, enabling precise object selection for high purity or completeness, and for identifying objects with complex features, like active galactic nuclei with resolved host galaxies. BANNJOS has effectively classified J-PLUS sources into around 20 million galaxies, 1 million QSOs, and 26 million stars, with full PDFs for each, which allow for later refinement of the sample. The upcoming J-PAS survey, with its 56 color bands, will further enhance BANNJOS's ability to detail each source's nature., Comment: Submitted to Astronomy and Astrophysics. 29 pages, 23 figures, 4 tables

We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resolution simulations. These models are benchmarked against comoving snapshots of the Flagship I N-body simulation at $z=(0.9,1.2,1.5,1.8)$, which have been populated with H$\alpha$ galaxies leading to catalogues of millions of objects within a volume of about $58\,h^{-3}\,{\rm Gpc}^3$. Our analysis suggests that both models can be used to provide a robust inference of the parameters $(h, \omega_{\rm c})$ in the redshift range under consideration, with comparable constraining power. We additionally determine the range of validity of the EFTofLSS model in terms of scale cuts and model degrees of freedom. From these tests, it emerges that the standard third-order Eulerian bias expansion can accurately describe the full shape of the real-space galaxy power spectrum up to the maximum wavenumber $k_{\rm max}=0.45\,h\,{\rm Mpc}^{-1}$, even with a measurement precision well below the percent level. In particular, this is true for a configuration with six free nuisance parameters, including local and non-local bias parameters, a matter counterterm, and a correction to the shot-noise contribution. Fixing either tidal bias parameters to physically-motivated relations still leads to unbiased cosmological constraints. We finally repeat our analysis assuming a volume that matches the expected footprint of Euclid, but without considering observational effects, as purity and completeness, showing that we can get consistent cosmological constraints over this range of scales and redshifts., Comment: 38 pages, 19 figures

We employ the corrected Gaia Early Data Release 3 (EDR3) photometric data and spectroscopic data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR7 to assemble a sample of approximately 0.25 million FGK dwarf photometric standard stars for the 12 J-PLUS filters using the Stellar Color Regression (SCR) method. We then independently validated the J-PLUS DR3 photometry, and uncovered significant systematic errors: up to 15 mmag in the results of Stellar Locus (SL) method, and up to 10 mmag mainly caused by magnitude-, color-, and extinction-dependent errors of the Gaia XP spectra with the Gaia BP/RP (XP) Synthetic Photometry (XPSP) method. We have also further developed the XPSP method using the corrected Gaia XP spectra by Huang et al. (2023) and applied it to the J-PLUS DR3 photometry. This resulted in an agreement of 1-5 mmag with the SCR method, and a two-fold improvement in the J-PLUS zero-point precision. Finally, the zero-point calibration for around 91% of the tiles within the LAMOST observation footprint is determined through the SCR method, with the remaining approximately 9% of tiles outside this footprint relying on the improved XPSP method. The re-calibrated J-PLUS DR3 photometric data establishes a solid data foundation for conducting research that depends on high-precision photometric calibration., Comment: ApJS accepted; 21 papes; 20 figures, see main results in Figures 5 and 12

We present the photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) third data release (DR3), comprising 1642 pointings of two square degrees each. We selected nearly 1.5 million main sequence stars with a signal-to-noise ratio larger than ten in the twelve J-PLUS passbands and available low-resolution (R = 20-80) spectrum from the blue and red photometers (BP/RP) in Gaia DR3. We compared the synthetic photometry from BP/RP spectra with the J-PLUS instrumental magnitudes, after correcting for the magnitude and color terms between both systems, to obtain an homogeneous photometric solution for J-PLUS. To circumvent the current limitations in the absolute calibration of the BP/RP spectra, the absolute color scale was derived using the locus of 109 white dwarfs closer than 100 pc with a negligible interstellar extinction. Finally, the absolute flux scale was anchored to the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) photometry in the r band. The precision of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and by comparison with the spectro-photometric standard star GD 153, is ~12 mmag in u, J0378, and J0395; and ~7 mmag in J0410, J0430, g, J0515, r, J0660, i, J0861, and z. The estimated accuracy in the calibration along the surveyed area is better than 1% for all the passbands. The Gaia BP/RP spectra provide a high-quality, homogeneous photometric reference in the optical range across the full-sky, in spite of their current limitations as an absolute reference. The calibration method for J-PLUS DR3 reaches an absolute precision and accuracy of 1% in the twelve optical filters within an area of 3284 square degrees., Comment: Submitted to Astronomy and Astrophysics. 17 pages, 14 figures, 2 appendix. Comments are welcome. arXiv admin note: text overlap with arXiv:2101.12407

The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliability of numerical methods for treating massive neutrinos. The comparison includes eleven full $N$-body implementations (not all of them independent), two $N$-body schemes with approximate time integration, and four additional codes that directly predict or emulate the matter power spectrum. Using a common set of initial data we quantify the relative agreement on the nonlinear power spectrum of cold dark matter and baryons and, for the $N$-body codes, also the relative agreement on the bispectrum, halo mass function, and halo bias. We find that the different numerical implementations produce fully consistent results. We can therefore be confident that we can model the impact of massive neutrinos at the sub-percent level in the most common summary statistics. We also provide a code validation pipeline for future reference., Comment: 44 pages, 17 figures, 2 tables; v2: minor revision, accepted manuscript; published on behalf of the Euclid Consortium; data available at https://doi.org/10.5281/zenodo.7868793

Ultracool dwarfs (UCDs) comprise the lowest mass members of the stellar population and brown dwarfs, from M7 V to cooler objects with L, T, and Y spectral types. Most of them have been discovered using wide-field imaging surveys, for which the Virtual Observatory (VO) has proven to be of great utility. We aim to perform a search for UCDs in the entire Javalambre Photometric Local Universe Survey (J-PLUS) second data release (2176 deg$^2$) following a VO methodology. We also explore the ability to reproduce this search with a purely machine learning (ML)-based methodology that relies solely on J-PLUS photometry. We followed three different approaches based on parallaxes, proper motions, and colours, respectively, using the VOSA tool to estimate the effective temperatures. For the ML methodology, we built a two-step method based on principal component analysis and support vector machine algorithms. We identified a total of 7827 new candidate UCDs, which represents an increase of about 135% in the number of UCDs reported in the sky coverage of the J-PLUS second data release. Among the candidate UCDs, we found 122 possible unresolved binary systems, 78 wide multiple systems, and 48 objects with a high Bayesian probability of belonging to a young association. We also identified four objects with strong excess in the filter corresponding to the Ca II H and K emission lines and four other objects with excess emission in the H$\alpha$ filter. With the ML approach, we obtained a recall score of 92% and 91% in the test and blind test, respectively. We consolidated the proposed search methodology for UCDs, which will be used in deeper and larger upcoming surveys such as J-PAS and Euclid. We concluded that the ML methodology is more efficient in the sense that it allows for a larger number of true negatives to be discarded prior to analysis with VOSA, although it is more photometrically restrictive., Comment: Accepted in A&A

Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be subdominant with respect to other sources of uncertainty in recovering cosmological parameters from Euclid cluster counts. Our model is calibrated against a suite of N-body simulations using a Bayesian approach taking into account systematic errors arising from numerical effects in the simulation. First, we test the convergence of HMF predictions from different N-body codes, by using initial conditions generated with different orders of Lagrangian Perturbation theory, and adopting different simulation box sizes and mass resolution. Then, we quantify the effect of using different halo-finder algorithms, and how the resulting differences propagate to the cosmological constraints. In order to trace the violation of universality in the HMF, we also analyse simulations based on initial conditions characterised by scale-free power spectra with different spectral indexes, assuming both Einstein--de Sitter and standard $\Lambda$CDM expansion histories. Based on these results, we construct a fitting function for the HMF that we demonstrate to be sub-percent accurate in reproducing results from 9 different variants of the $\Lambda$CDM model including massive neutrinos cosmologies. The calibration systematic uncertainty is largely sub-dominant with respect to the expected precision of future mass-observation relations; with the only notable exception of the effect due to the halo finder, that could lead to biased cosmological inference., Comment: 25 pages, 21 figures, 5 tables, 3 appendixes; v2 matches published version

Context. Stellar parameters are among the most important characteristics in studies of stars, which are based on atmosphere models in traditional methods. However, time cost and brightness limits restrain the efficiency of spectral observations. The J-PLUS is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resource for studies of stars. Machine learning provides powerful tools to efficiently analyse large data sets, such as the one from J-PLUS, and enable us to expand the research domain to stellar parameters. Aims. The main goal of this study is to construct a SVR algorithm to estimate stellar parameters of the stars in the first data release of the J-PLUS observational campaign. Methods. The training data for the parameters regressions is featured with 12-waveband photometry from J-PLUS, and is cross-identified with spectrum-based catalogs. These catalogs are from the LAMOST, the APOGEE, and the SEGUE. We then label them with the stellar effective temperature, the surface gravity and the metallicity. Ten percent of the sample is held out to apply a blind test. We develop a new method, a multi-model approach in order to fully take into account the uncertainties of both the magnitudes and stellar parameters. The method utilizes more than two hundred models to apply the uncertainty analysis. Results. We present a catalog of 2,493,424 stars with the Root Mean Square Error of 160K in the effective temperature regression, 0.35 in the surface gravity regression and 0.25 in the metallicity regression. We also discuss the advantages of this multi-model approach and compare it to other machine-learning methods., Comment: 18 pages, 10 figures, accepted by AA

Over the past decades, several studies have discovered a population of galaxies undergoing very strong star formation events, called extreme emission line galaxies (EELGs). In this work, we exploit the capabilities of the Javalambre Photometric Local Universe Survey (J-PLUS), a wide field multifilter survey, with 2000 square degrees observed. We use it to identify EELGs at low redshift by their [OIII]5007 emission line. We intend to provide with a more complete, deep, and less biased sample of local EELGs. We select objects with an excess of flux in the J-PLUS mediumband $J0515$ filter, which covers the [OIII] line at z$<$0.06. We remove contaminants (stars and higher redshift systems) using J-PLUS and WISE infrared data, with SDSS spectra as a benchmark. We perform spectral energy distribution fitting to estimate the properties of the galaxies: line fluxes, equivalent widths (EWs), masses, etc. We identify 466 EELGs at ${\rm z} < 0.06$ with [OIII] EW over 300 \text{\AA} and $r$-band mag. below 20, of which 411 were previously unknown. Most show compact morphologies, low stellar masses ($\log (M_{\star}/M_{\odot}) \sim {8.13}^{+0.61}_{-0.58}$), low dust extinction ($E(B-V)\sim{0.1}^{+0.2}_{-0.1}$), and very young bursts of star formation (${3.0}^{+2.7}_{-2.0}$ Myr). Our method is up to $\sim$ 20 times more efficient detecting EELGs per Mpc$^3$ than broadband surveys, and as complete as magnitude-limited spectroscopic surveys (and reaching fainter objects). The sample is not directly biased against strong H$\alpha$ emitters, in contrast with broadband surveys. We demonstrate the capability of J-PLUS to identify, following a clear selection process, a large sample of previously unknown EELGs showing unique properties. A fraction of them are likely similar to the first galaxies in the Universe, but at a much lower redshift, which makes them ideal targets for follow-up studies., Comment: 22 pages, 22 figures. Submitted to Astronomy & Astrophysics

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), that provides twelve photometric optical passbands over 2176 deg2. We analysed 5926 white dwarfs with r <= 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen (H) and helium (He) dominated atmospheres. We estimated the PDF for effective temperature (Teff), surface gravity, parallax, and spectral type; and the probability of having a H-dominated atmosphere (pH) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as reference, and derived a self-consistent prior for the atmospheric composition as a function of Teff. We described the fraction of He-dominated atmosphere white dwarfs (fHe) with a linear function of Teff at 5000 < Teff < 30000 K. We found fHe = 0.24 +- 0.01 at Teff = 10000 K, a change rate along the cooling sequence of 0.14 +- 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 +- 0.02 at the high-temperature end. We tested the obtained pH by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 Msun, and Teff > 6000 K. The result for H-dominated white dwarfs agrees with previous work, with a dominant M = 0.59 Msun peak and the presence of an excess at M ~ 0.8 Msun. This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provides a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 +- 3 % increase in the fraction of He-dominated white dwarfs from Teff = 20000 K to Teff = 5000 K., Comment: Published in Astronomy and Astrophysics. 18 pages, 15 figures, 4 tables

Context. In modern astronomy, machine learning has proved to be efficient and effective to mine the big data from the newesttelescopes. Spectral surveys enable us to characterize millions of objects, while long exposure time observations and wide surveysconstrain their strides from millions to billions. Aims.In this study, we construct a supervised machine learning algorithm, to classify the objects in the Javalambre Photometric LocalUniverse Survey first data release (J-PLUS DR1). Methods.The sample set is featured with 12-waveband photometry, and magnitudes are labeled with spectrum-based catalogs, in-cluding Sloan Digital Sky Survey spectroscopic data, Large Sky Area Multi-Object Fiber Spectroscopic Telescope, and VERONCAT- Veron Catalog of Quasars & AGN. The performance of the classifier is presented with applications of blind test validations basedon RAdial Velocity Extension, Kepler Input Catalog, 2 MASS Redshift Survey, and the UV-bright Quasar Survey. A new algorithmis applied to constrain the extrapolation that could decrease accuracies for many machine learning classifiers. Results.The accuracies of the classifier are 96.5% in blind test and 97.0% in training cross validation. The F1-scores for each classare presented to show the precision of the classifier. We also discuss different methods to constrain the po, Comment: 26 pages,28 figures, accepted by A&A

We present the updated photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), comprising 1088 pointings of two square degrees, and study the systematic impact of metallicity in the stellar locus technique. The [Fe/H] metallicity from LAMOST DR5 for 146184 high-quality calibration stars, defined with S/N > 10 in J-PLUS passbands and S/N > 3 in Gaia parallax, was used to compute the metallicity-dependent stellar locus (ZSL). The initial homogenization of J-PLUS photometry, performed with a unique stellar locus, was refined by including the metallicity effect in colours via the ZSL. The variation of the average metallicity along the Milky Way produces a systematic offset in J-PLUS calibration. This effect is well above 1% for the bluer passbands and amounts 0.07, 0.07, 0.05, 0.03, and 0.02 mag in u, J0378, J0395, J0410, and J0430, respectively. We modelled this effect with the Milky Way location of the J-PLUS pointing, providing also an updated calibration for those observations without LAMOST information. The estimated accuracy in the calibration after including the metallicity effect is at 1% level for the bluer J-PLUS passbands and below for the rest. We conclude that photometric calibration with the stellar locus technique is prone to significant systematic bias along the Milky Way location for passbands bluer than lambda = 4500 A. The updated calibration method for J-PLUS DR2 reaches 1-2% precision and 1% accuracy for twelve optical filters within an area of 2176 square degrees., Comment: Published in Astronomy & Astrophysics. 18 pages, 12 figures, 3 tables, 2 appendices

Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is currently possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_\chi > 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be correct as other structures, such as filaments and pancakes, exist and make up a significant fraction of the mass in the universe, ranging between 5$\%$ -- 50$\%$ depending on the dark matter model. Using novel fragmentation-free simulations of 1 and 3keV WDM cosmologies we study these "non-halo" structures and estimate their impact on flux-ratio observations. We find that these structures display sharp density gradients with short correlation lengths, and can contribute more to the lensing signal than all haloes up to the half-mode mass combined, thus reducing the differences expected among WDM models. We estimate that non-halo structures can be the dominant cause of line-of-sight flux-ratio anomalies in very warm, but already excluded, $m_x \sim 1 \rm{keV}$ scenarios. For colder cases $m_x \gtrsim 3 \rm{keV}$, we estimate that non-haloes can contribute about $5 - 10\%$ of the total flux-ratio signal., Comment: 14 pages, 10 figures. MNRAS submitted. Comments welcome

Our goal is to estimate the star formation main sequence (SFMS) and the star formation rate density (SFRD) at z <= 0.017 (d < 75 Mpc) using the Javalambre Photometric Local Universe Survey (J-PLUS) first data release, that probes 897.4 deg2 with twelve optical bands. We extract the Halpha emission flux of 805 local galaxies from the J-PLUS filter J0660, being the continuum level estimated with the other eleven J-PLUS bands, and the dust attenuation and nitrogen contamination corrected with empirical relations. Stellar masses (M), Halpha luminosities (L), and star formation rates (SFRs) were estimated by accounting for parameters covariances. Our sample comprises 689 blue galaxies and 67 red galaxies, classified in the (u-g) vs (g-z) color-color diagram, plus 49 AGN. The SFMS is explored at log M > 8 and it is clearly defined by the blue galaxies, with the red galaxies located below them. The SFMS is described as log SFR = 0.83 log M - 8.44. We find a good agreement with previous estimations of the SFMS, especially those based on integral field spectroscopy. The Halpha luminosity function of the AGN-free sample is well described by a Schechter function with log L* = 41.34, log phi* = -2.43, and alpha = -1.25. Our measurements provide a lower characteristic luminosity than several previous studies in the literature. The derived star formation rate density at d < 75 Mpc is log rho_SFR = -2.10 +- 0.11, with red galaxies accounting for 15% of the SFRD. Our value is lower than previous estimations at similar redshift, and provides a local reference for evolutionary studies regarding the star formation history of the Universe., Comment: 26 pages. Accepted for publication in Astronomy and Astrophysics

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged), Comment: The miniJPAS data and associated value added catalogues are publicly accessible via this url: http://archive.cefca.es/catalogues/minijpas-pdr201912

Predicting the spatial distribution of objects as a function of cosmology is an essential ingredient for the exploitation of future galaxy surveys. In this paper we show that a specially-designed suite of gravity-only simulations together with cosmology-rescaling algorithms can provide the clustering of dark matter, haloes, and subhaloes with high precision. Specifically, with only 3 $N$-body simulations we obtain the power spectrum of dark matter at $z=0$ and $z=1$ to better than 3\% precision for essentially all currently viable values of 8 cosmological parameters, including massive neutrinos and dynamical dark energy, and over the whole range of scales explored, 0.03 < $k/h^{-1}Mpc$ < 5. This precision holds at the same level for mass-selected haloes and for subhaloes selected according to their peak maximum circular velocity. As an initial application of these predictions, we successfully constrain $\Omega_{\rm m}$, $\sigma_8$, and the scatter in subhalo-abundance-matching employing the projected correlation function of mock SDSS galaxies., Comment: 12 pages, 8 figures. Submitted. Comments are welcome

From the approximately $\sim$3,500 planetary nebulae (PNe) discovered in our Galaxy, only 14 are known to be members of the Galactic halo. Nevertheless, a systematic search for halo PNe has never been performed. In this study, we present new photometric diagnostic tools to identify compact PNe in the Galactic halo by making use of the novel 12-filter system projects, J-PLUS (Javalambre Photometric Local Universe Survey) and S-PLUS (Southern-Photometric Local Universe Survey). We reconstructed the IPHAS (Isaac Newton Telescope (INT) Photometric H${\alpha}$ Survey of the Northern Galactic Plane) diagnostic diagram and propose four new ones using i) the J-PLUS and S-PLUS synthetic photometry for a grid of photo-ionisation models of halo PNe, ii) several observed halo PNe, as well as iii) a number of other emission-line objects that resemble PNe. All colour-colour diagnostic diagrams are validated using two known halo PNe observed by J-PLUS during the scientific verification phase and the first data release (DR1) of S-PLUS and the DR1 of J-PLUS. By applying our criteria to the DR1s ($\sim$1,190 deg$^2$), we identified one PN candidate. However, optical follow-up spectroscopy proved it to be a H II region belonging to the UGC 5272 galaxy. Here, we also discuss the PN and two H II galaxies recovered by these selection criteria. Finally, the cross-matching with the most updated PNe catalogue (HASH) helped us to highlight the potential of these surveys, since we recover all the known PNe in the observed area. The tools here proposed to identify PNe and separate them from their emission-line contaminants proved to be very efficient thanks to the combination of many colours, even when applied -like in the present work- to an automatic photometric search that is limited to compact PNe., Comment: 12 pages, 10 figures, accepted for publication at A&A

We present the photometric calibration of the twelve optical passbands observed by the Javalambre Photometric Local Universe Survey (J-PLUS). The proposed calibration method has four steps: (i) definition of a high-quality set of calibration stars using Gaia information and available 3D dust maps; (ii) anchoring of the J-PLUS gri passbands to the Pan-STARRS photometric solution, accounting for the variation of the calibration with the position of the sources on the CCD; (iii) homogenization of the photometry in the other nine J-PLUS filters using the dust de-reddened instrumental stellar locus in (X - r) versus (g - i) colours, where X is the filter to calibrate. The zero point variation along the CCD in these filters was estimated with the distance to the stellar locus. Finally, (iv) the absolute colour calibration was obtained with the white dwarf locus. We performed a joint Bayesian modelling of eleven J-PLUS colour-colour diagrams using the theoretical white dwarf locus as reference. This provides the needed offsets to transform instrumental magnitudes to calibrated magnitudes outside the atmosphere. The uncertainty of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and accounting for the absolute colour and flux calibration errors, are ~19 mmag in u, J0378 and J0395, ~11 mmag in J0410 and J0430, and ~8 mmag in g, J0515, r, J0660, i, J0861, and z. We present an optimized calibration method for the large area multi-filter J-PLUS project, reaching 1-2% accuracy within an area of 1 022 square degrees without the need for long observing calibration campaigns or constant atmospheric monitoring. The proposed method will be adapted for the photometric calibration of J-PAS, that will observe several thousand square degrees with 56 narrow optical filters., Comment: Submitted to Astronomy and Astrophysics. 14 figures, 5 tables. Comments are welcome. Extra information about the photometry update can be found at www.j-plus.es/datareleases/dr1_swdl_calibration

In a framework where galaxies mostly migrate on the colour-magnitude diagram from star-forming to quiescent, the green valley is considered a transitional galaxy stage. The details of the processes that drive galaxies from star-forming to passive systems still remain unknown. We developed a method that estimates empirically the star formation quenching times-scales of green valley galaxies, assuming an exponential decay model of the SFH and through a combination of narrow and broad bands from J-PLUS and GALEX. We correlate these quenching time-scales with the presence of bars. We find that the J-PLUS colours F0395-g and F0415-g are sensitive to different SFH, showing, a clear correlation with the Dn(4000) and H-delta,A spectral indices. We find that quenching time-scales obtained with our new approach are in agreement with those determined using spectral indices. We also find that galaxies with high bar probability tend to quench their star formation slowly. We conclude that: 1) J-PLUS filters can be used to measure quenching timescales in nearby green valley galaxies; and 2) the resulting star formation quenching time-scales are longer for barred green valley galaxies. Considering that the presence of a bar indicates that more violent processes (e.g., major mergers) are absent in host galaxies, we conclude that the presence of a bar can be used as a morphological signature for slow star formation quenching., Comment: 14 pages, 17 figures, accepted by Astronomy and Astrophysics

Ultracool dwarfs (UCDs) are objects with spectral types equal or later than M7. Most of them have been discovered using wide-field imaging surveys. The Virtual Observatory (VO) has proven to be of great utility to efficiently exploit these astronomical resources. We aim to validate a VO methodology designed to discover and characterize UCDs in the J-PLUS photometric survey. J-PLUS is a multiband survey carried out with the wide angle T80Cam optical camera mounted on the 0.83-m telescope JAST/T80 in the Observatorio Astrof\'isico de Javalambre. In this work we make use of the Internal Data Release (IDR) covering 528 deg$^2$. We complement J-PLUS photometry with other catalogues in the optical and IR using VOSA, a VO tool that estimates physical parameters from the spectral energy distribution fitting to collections of theoretical models. Objects identified as UCDs are distinguished from background M giants and highly reddened stars using parallaxes and proper motions from Gaia DR2. We identify 559 UCDs, ranging from i=16.2 to 22.4 mag, of which 187 are candidate UCDs not previously reported in the literature. This represents an increase in the number of known UCDs of about 50% in the studied region of the sky, particularly at the faint end of our sensitivity, which is interesting as reference for future wide and deep surveys such as Euclid. Three candidates constitute interesting targets for exoplanet surveys because of their proximity (<40 pc). We also analyze the kinematics of UCDs in our catalogue and find evidence that it is consistent with a Galactic thin-disk population, except for 6 objects that might be members of the thick disk. The results obtained validate the proposed methodology, which will be used in future J-PLUS and J-PAS releases. Considering the region of the sky covered by the IDR used, we foresee to discover 3,000-3,500 new UCDs at the end of the J-PLUS project., Comment: Accepted in A&A

Our goal is to morphologically classify the sources identified in the images of the J-PLUS early data release (EDR) into compact (stars) or extended (galaxies) using a suited Bayesian classifier. J-PLUS sources exhibit two distinct populations in the r-band magnitude vs. concentration plane, corresponding to compact and extended sources. We modelled the two-population distribution with a skewed Gaussian for compact objects and a log-normal function for the extended ones. The derived model and the number density prior based on J-PLUS EDR data were used to estimate the Bayesian probability of a source to be star or galaxy. This procedure was applied pointing-by-pointing to account for varying observing conditions and sky position. Finally, we combined the morphological information from g, r, and i broad bands in order to improve the classification of low signal-to-noise sources. The derived probabilities are used to compute the pointing-by-pointing number counts of stars and galaxies. The former increases as we approach to the Milky Way disk, and the latter are similar across the probed area. The comparison with SDSS in the common regions is satisfactory up to r ~ 21, with consistent numbers of stars and galaxies, and consistent distributions in concentration and (g - i) colour spaces. We implement a morphological star/galaxy classifier based on PDF analysis, providing meaningful probabilities for J-PLUS sources to one magnitude deeper (r ~ 21) than a classical boolean classification. These probabilities are suited for the statistical study of 150k stars and 101k galaxies with 15 < r < 21 present in the 31.7 deg2 of the J-PLUS EDR. In a future version of the classifier, we will include J-PLUS colour information from twelve photometric bands., Comment: Submitted to Astronomy and Astrophysics. 14 pages, 16 figures, 1 tables. Comments are welcome. All extra figures and the number counts files will be available with the paper in press

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag., Comment: Submitted to A&A

The cosmological parameters prefered by the cosmic microwave background (CMB) primary anisotropies predict many more galaxy clusters than those that have been detected via the thermal Sunyaev-Zeldovich (tSZ) effect. This tension has attracted considerable attention since it could be evidence of physics beyond the simplest $\Lambda$CDM model. However, an accurate and robust calibration of the mass-observable relation for clusters is necessary for the comparison, which has been proven difficult to obtain so far. Here, we present new contraints on the mass-pressure relation by combining tSZ and CMB lensing measurements about optically-selected clusters. Consequently, our galaxy cluster sample is independent from the data employed to derive cosmological constrains. We estimate an average hydrostatic mass bias of $b = 0.26 \pm 0.07$, with no significant mass nor redshift evolution. This value greatly reduces the tension between the predictions of $\Lambda$CDM and the observed abundance of tSZ clusters while being in agreement with recent estimations from tSZ clustering. On the other hand, our value for $b$ is higher than the predictions from hydro-dynamical simulations. This suggests the existence of mechanisms driving large departures from hydrostatic equilibrium and that are not included in state-of-the-art simulations, and/or unaccounted systematic errors such as biases in the cluster catalogue due to the optical selection., Comment: 4 pages, 3 figures

Simulations of galaxy evolution aim to capture our current understanding as well as to make predictions for testing by future experiments. Simulations and observations are often compared in an indirect fashion: physical quantities are estimated from the data and compared to models. However, many applications can benefit from a more direct approach, where the observing process is also simulated and the models are seen fully from the observer's perspective. To facilitate this, we have developed the Millennium Run Observatory (MRObs), a theoretical virtual observatory which uses virtual telescopes to `observe' semi-analytic galaxy formation models based on the suite of Millennium Run dark matter simulations. The MRObs produces data that can be processed and analyzed using the standard software packages developed for real observations. At present, we produce images in forty filters from the rest-frame UV to IR for two stellar population synthesis models, three different models of IGM absorption, and two cosmologies (WMAP1/7). Galaxy distributions for a large number of mock lightcones can be `observed' using models of major ground- and space-based telescopes. The data include lightcone catalogues linked to structural properties of galaxies, pre-observation model images, mock telescope images, and Source Extractor products that can all be traced back to the higher level dark matter, semi-analytic galaxy, and lightcone catalogues available in the Millennium database. Here, we describe our methods and announce a first public release of simulated surveys (e.g., SDSS, CFHT-LS, GOODS, GOODS/ERS, CANDELS, and HUDF). The MRObs browser, an online tool, further facilitates exploration of the simulated data. We demonstrate the benefits of a direct approach through a number of example applications (galaxy number counts in CANDELS, clusters, morphologies, and dropout selections)., Comment: MNRAS, in press. Millennium Run Observatory data products, online tools, and more available through http://galformod.mpa-garching.mpg.de/mrobs/

We apply a recently developed scaling technique to the Millennium-XXL, one of the largest cosmological N-body simulations carried out to date 3x10^11 particles within a cube of volume ~70Gpc^3). This allows us to investigate the cosmological parameter dependence of the mass and evolution of haloes in the extreme high-mass tail of the z=6 distribution. We assume these objects to be likely hosts for the population of rare but ultraluminous high-redshift quasars discovered by the Sloan Digital Sky Survey. Haloes with a similar abundance to these quasars have a median mass of 9x10^12 Msun in the currently preferred cosmology, but do not evolve into equally extreme objects at z=0. Rather, their descendants span the full range conventionally assigned to present-day clusters, 6x10^13 to 2.5x10^15 Msun for this same cosmology. The masses both at z=6 and at z=0 shift up or down by factors exceeding two if cosmological parameters are pushed to the boundaries of the range discussed in published interpretations of data from the WMAP satellite. The main factor determining the future growth of a high-mass z=6 halo is the mean overdensity of its environment on scales of 7 to 14 Mpc, and descendant masses can be predicted 6 to 8 times more accurately if this density is known than if it is not. All these features are not unique to extreme high-z haloes, but are generic to hierarchical growth. Finally, we find that extreme haloes at z=6 typically acquired about half of their total mass in the preceding 100 Myr, implying very large recent accretion rates which may be related to the large black hole masses and high luminosities of the SDSS quasars., Comment: 10 pages, 6 figures. Replaced with version published in MNRAS

We present a very large high-resolution cosmological N-body simulation, the Millennium-XXL or MXXL, which uses 303 billion particles to represent the formation of dark matter structures throughout a 4.1Gpc box in a LambdaCDM cosmology. We create sky maps and identify large samples of galaxy clusters using surrogates for four different observables: richness estimated from galaxy surveys, X-ray luminosity, integrated Sunyaev-Zeldovich signal, and lensing mass. The unprecedented combination of volume and resolution allows us to explore in detail how these observables scale with each other and with cluster mass. The scatter correlates between different mass-observable relations because of common sensitivities to the internal structure, orientation and environment of clusters, as well as to line-of-sight superposition of uncorrelated structure. We show that this can account for the apparent discrepancies uncovered recently between the mean thermal SZ signals measured for optically and X-ray selected clusters by stacking data from the Planck satellite. Related systematics can also affect inferences from extreme clusters detected at high redshift. Our results illustrate that cosmological conclusions from galaxy cluster surveys depend critically on proper modelling, not only of the relevant physics, but also of the full distribution of the observables and of the selection biases induced by cluster identification procedures., Comment: 19 pages, 12 figures. Replaced with version published in MNRAS

We propose a new method to linearise cosmological mass density fields using higher order Lagrangian perturbation theory (LPT). We demonstrate that a given density field can be expressed as the sum of a linear and a nonlinear component which are tightly coupled to each other by the tidal field tensor within the LPT framework. The linear component corresponds to the initial density field in Eulerian coordinates, and its mean relation with the total field can be approximated by a logarithm (giving theoretical support to recent attempts to find such component). We also propose to use a combination of the linearisation method and the continuity equation to find the mapping between Eulerian and Lagrangian coordinates. In addition, we note that this method opens the possibility of use directly higher order LPT on nonlinear fields. We test our linearization scheme by applying it to the z~0.5 density field from an N-body simulation. We find that the linearised version of the full density field can be successfully recovered on >~5 h^{-1}Mpc, reducing the skewness and kurtosis of the distribution by about one and two orders of magnitude, respectively. This component can also be successfully traced back in time, converging towards the initial unevolved density field at z~100. We anticipate a number of applications of our results, from predicting velocity fields to estimates of the initial conditions of the universe, passing by improved constraints on cosmological parameters derived from galaxy clustering via reconstruction methods., Comment: 14 pages, 8 figures

We introduce an objective method to assess the probability of finding extreme events in the distribution of cold dark matter such as voids, overdensities or very high mass haloes. Our approach uses an ensemble of N-body simulations of the hierar- chical clustering of dark matter to find extreme structures. The frequency of extreme events, in our case the cell or smoothing volume with the highest count of cluster-mass dark matter haloes, is well described by a Gumbel distribution. This distribution can then be used to forecast the probability of finding even more extreme events, which would otherwise require a much larger ensemble of simulations to quantify. We use our technique to assess the chance of finding concentrations of massive clusters or super- clusters, like the two found in the two-degree field galaxy redshift survey (2dFGRS), using a counts-in-cells analysis. The Gumbel distribution gives an excellent descrip- tion of the distribution of extreme cell counts across two large ensembles of simulations covering different cosmologies, and also when measuring the clustering in both real and redshift space. We find examples of structures like those found in the 2dFGRS in the simulations. The chance of finding such structures in a volume equal to that of the 2dFGRS is around 2%., Comment: 8 pages, 5 figures, modified to match accepted version in MNRAS

We model the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models which feature late $(z<2)$, rapid transitions towards $w=-1$ in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $\Lambda$CDM, despite being very different from $\Lambda$CDM both today and at high redshifts $(z \sim 1000)$. We find that a second class of models which feature non-negligible amounts of dark energy at early times cannot be distinguished from $\Lambda$CDM using measurements of the mass function or the BAO. These results highlight the need to accurately model quintessence dark energy in N-body simulations when testing cosmological probes of dynamical dark energy., Comment: 10 pages, 7 figures, to appear in the Invisible Univers International Conference AIP proceedings series

We present the results of a search for variable stars in the globular cluster NGC 5286, which has recently been suggested to be associated with the Canis Major dwarf spheroidal galaxy. 57 variable stars were detected, only 19 of which had previously been known. Among our detections one finds 52 RR Lyrae (22 RRc and 30 RRab), 4 LPV's, and 1 type II Cepheid of the BL Herculis type. Periods are derived for all of the RR Lyrae as well as the Cepheid, and BV light curves are provided for all the variables. The mean period of the RRab variables is = 0.656 days, and the number fraction of RRc stars is N(c)/N(RR) = 0.42, both consistent with an Oosterhoff II (OoII) type -- thus making NGC 5286 one of the most metal-rich ([Fe/H] = -1.67; Harris 1996) OoII globulars known to date. The minimum period of the \RRab's, namely Pab,min = 0.513 d, while still consistent with an OoII classification, falls towards the short end of the observed Pab,min distribution for OoII globular clusters. As was recently found in the case of the prototypical OoII globular cluster M15 (NGC 7078), the distribution of stars in the Bailey diagram does not strictly conform to the previously reported locus for OoII stars. We provide Fourier decomposition parameters for all of the RR Lyrae stars detected in our survey, and discuss the physical parameters derived therefrom. The values derived for the RRc's are not consistent with those typically found for OoII clusters, which may be due to the cluster's relatively high metallicity -- the latter being confirmed by our Fourier analysis of the ab-type RR Lyrae light curves. We derive for the cluster a revised distance modulus of (m-M)V = 16.04 mag. (ABRIDGED), Comment: 15 pages, 8 figures (plus 57 light curve sets, in an Appendix). The Astronomical Journal, in press

We use the Extended-Press-Schechter (EPS) formalism to study halo assembly histories in a standard $\Lambda$CDM cosmology. A large ensemble of Monte Carlo random walks provides the {\it entire} halo membership histories of a representative set of dark matter particles, which we assume to be neutralinos. The first generation halos of most particles do not have a mass similar to the free-streaming cut-off $M_{f.s.}$ of the neutralino power spectrum, nor do they form at high redshift. Median values are $M_1 = 10^5$ to $10^7M_{f.s.}$ and $z_1 = 13$ to 8 depending on the form of the collapse barrier assumed in the EPS model. For almost a third of all particles the first generation halo has $M_1>10^9M_{f.s.}$. At redshifts beyond 20, most neutralinos are not yet part of any halo but are still diffuse. These numbers apply with little modification to the neutralinos which are today part of halos similar to that of the Milky Way. Up to 10% of the particles in such halos were never part of a smaller object; the typical particle has undergone $\sim 5$ "accretion events' where the halo it was part of falls into a more massive object. Available N-body simulations agree well with the EPS predictions for an "ellipsoidal" collapse barrier, so these may provide a reliable extension of simulation results to smaller scales. The late formation times and large masses of the first generation halos of most neutralinos imply that they will be disrupted with high efficiency during halo assembly., Comment: 7 pages, 7 figures

We use the Millennium Simulation, a large, high resolution N-body simulation of the evolution of structure in a LambdaCDM cosmology, to study the properties and fate of substructures within a large sample of dark matter haloes. We find that the subhalo mass function departs significantly from a power law at the high mass end. We also find that the radial and angular distributions of substructures depend on subhalo mass. In particular, high mass subhaloes tend to be less radially concentrated and to have angular distributions closer to the direction perpendicular to the spin of the host halo than their less massive counterparts. We find that mergers between subhaloes occur. These tend to be between substructures that were already dynamically associated before accretion into the main halo. For subhaloes larger than 0.001 times the mass of the host halo, it is more likely that the subhalo will merge with the central or main subhalo than with another subhalo larger than itself. For lower masses, subhalo-subhalo mergers become equally likely to mergers with the main subhalo. Our results have implications for the variation of galaxy properties with environment and for the treatment of mergers in galaxy formation models., Comment: 14 pages, 11 figures. Submitted to MNRAS

We use an extremely large volume ($2.4h^{-3}{\rm Gpc}^{3}$), high resolution N-body simulation to measure the higher order clustering of dark matter haloes as a function of mass and internal structure. As a result of the large simulation volume and the use of a novel ``cross-moment'' counts-in-cells technique which suppresses discreteness noise, we are able to measure the clustering of haloes corresponding to rarer peaks than was possible in previous studies; the rarest haloes for which we measure the variance are 100 times more clustered than the dark matter. We are able to extract, for the first time, halo bias parameters from linear up to fourth order. For all orders measured, we find that the bias parameters are a strong function of mass for haloes more massive than the characteristic mass $M_{*}$. Currently, no theoretical model is able to reproduce this mass dependence closely. We find that the bias parameters also depend on the internal structure of the halo up to fourth order. For haloes more massive than $M_{*}$, we find that the more concentrated haloes are more weakly clustered than the less concentrated ones. We see no dependence of clustering on concentration for haloes with masses $M