Galaxy clusters appear as extended sources in XMM-Newton images, but not all extended sources are clusters. So, their proper classification requires visual inspection with optical images, which is a slow process with biases that are almost impossible to model. We tackle this problem with a novel approach, using convolutional neural networks (CNNs), a state-of-the-art image classification tool, for automatic classification of galaxy cluster candidates. We train the networks on combined XMM-Newton X-ray observations with their optical counterparts from the all-sky Digitized Sky Survey. Our data set originates from the X-CLASS survey sample of galaxy cluster candidates, selected by a specially developed pipeline, the XAmin, tailored for extended source detection and characterisation. Our data set contains 1 707 galaxy cluster candidates classified by experts. Additionally, we create an official Zooniverse citizen science project, The Hunt for Galaxy Clusters, to probe whether citizen volunteers could help in a challenging task of galaxy cluster visual confirmation. The project contained 1 600 galaxy cluster candidates in total of which 404 overlap with the expert's sample. The networks were trained on expert and Zooniverse data separately. The CNN test sample contains 85 spectroscopically confirmed clusters and 85 non-clusters that appear in both data sets. Our custom network achieved the best performance in the binary classification of clusters and non-clusters, acquiring accuracy of 90 %, averaged after 10 runs. The results of using CNNs on combined X-ray and optical data for galaxy cluster candidate classification are encouraging and there is a lot of potential for future usage and improvements.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogue in its present state, as well as of associated catalogues as super-clusters and fossil groups. We release a sample containing 365 clusters in total. We give the details of the follow-up observations and explain the procedure adopted to validate the cluster spectroscopic redshifts. Considering the whole XXL cluster sample, we have provided two types of selection, both complete in a particular sense: one based on flux-morphology criteria, and an alternative based on the [0.5-2] keV flux within one arcmin of the cluster centre. We have also provided X-ray temperature measurements for 80$\%$ of the clusters having a flux larger than 9$\times$10$^{-15}$$\rm \thinspace erg \, s^{-1} \, cm^{-2}$. Our cluster sample extends from z$\sim$0 to z$\sim$1.2, with one cluster at z$\sim$2. Clusters were identified through a mean number of six spectroscopically confirmed cluster members. Our updated luminosity function and luminosity-temperature relation are compatible with our previous determinations based on the 100 brightest clusters, but show smaller uncertainties. We also present an enlarged list of super-clusters and a sample of 18 possible fossil groups. This intermediate publication is the last before the final release of the complete XXL cluster catalogue when the ongoing C2 cluster spectroscopic follow-up is complete. It provides a unique inventory of medium-mass clusters over a 50~\dd\ area out to z$\sim$1.

We present the version of the point source catalogue of the XXL Survey that was used, in part, in the first series of XXL papers. In this paper we release, in our database in Milan and at CDS: (i) the X-ray source catalogue with 26056 objects in two areas of 25 deg2; (ii) the associated multiwavelength catalogues with candidate counterparts of the X-ray sources in the infrared, near-infrared, optical, and ultraviolet (plus spectroscopic redshift when available); and (iii) a catalogue of spectroscopic redshifts recently obtained in the southern XXL area. We also present the basic properties of the X-ray point sources and their counterparts. Other catalogues described in the second series of XXL papers will be released contextually, and will constitute the second XXL data release.

Context. We present an estimation of cosmological parameters with clusters of galaxies. Aims. We constrain the $\Omega_m$, $\sigma_8$, and $w$ parameters from a stand-alone sample of X-ray clusters detected in the 50 deg$^2$ XMM-XXL survey with a well-defined selection function. Methods. We analyse the redshift distribution of a sample comprising 178 high S/N clusters out to a redshift of unity. The cluster sample scaling relations are determined in a self-consistent manner. Results. In a lambda cold dark matter ($\Lambda$CDM) model, the cosmology favoured by the XXL clusters compares well with results derived from the Planck S-Z clusters for a totally different sample (mass/redshift range, selection biases, and scaling relations). However, with this preliminary sample and current mass calibration uncertainty, we find no inconsistency with the Planck CMB cosmology. If we relax the $w$ parameter, the Planck CMB uncertainties increase by a factor of $\sim$10 and become comparable with those from XXL clusters. Combining the two probes allows us to put constraints on $\Omega_m$=0.316$\pm$0.060, $\sigma_8$=0.814$\pm$0.054, and $w$=-1.02$\pm$0.20. Conclusions. This first self-consistent cosmological analysis of a sample of serendipitous XMM clusters already provides interesting insights into the constraining power of the XXL survey. Subsequent analysis will use a larger sample extending to lower confidence detections and include additional observable information, potentially improving posterior uncertainties by roughly a factor of 3.

A well characterised detection pipeline is an important ingredient for X-ray cluster surveys. We present the final development of the XXL Survey pipeline. The pipeline optimally uses X-ray information by combining many overlapping observations of a source when possible, both for its detection and its characterisation. It can robustly detect and characterise several types of X-ray sources: AGNs (point-like), galaxy clusters (extended), galaxy clusters contaminated by a central AGN, and pairs of AGNs close on the sky. We perform a thorough suite of validation tests via realistic simulations of XMM-Newton images and we introduce new selection criteria for various types of sources that will be detected by the survey. We find that the use of overlapping observations allows new clusters to be securely identified that would be missed or less securely identified by using only one observation at the time. We also find that with the new pipeline we can robustly identify clusters with a central AGN that would otherwise have been missed, and we can flag pairs of AGNs close on the sky that might have been mistaken for a cluster.

We present the results of a study of the AGN density in a homogeneous and well-studied sample of 167 bona fide X-ray galaxy clusters ($0.1<z<0.5$) from the XXL Survey. The results can provide evidence of the physical mechanisms that drive AGN and galaxy evolution within clusters. The XXL cluster sample mostly comprises poor and moderately rich structures ($M=10^{13} - 4\times10^{14} M_{\rm o}$). Our aim is to statistically study the demographics of cluster AGNs as a function of cluster mass and host galaxy position. To investigate the effect of the environment on AGN activity, we computed the fraction of spectroscopically confirmed X-ray AGNs ($L_{\rm X [0.5-10\,keV]}>10^{42}$ erg cm$^{-1}$) in bright cluster galaxies, up to $6r_{500}$ radius. To study the mass dependence and the evolution of the AGN population, we further divided the sample into low- and high-mass clusters and two redshift bins (0.1-0.28 and 0.28-0.5). We detect a significant excess of X-ray AGNs, at the 95% confidence level, in low-mass clusters between $0.5r_{500}$ and 2$r_{500}$, which drops to the field value within the cluster cores ($r<0.5r_{500}$). In contrast, high-mass clusters present a decreasing AGN fraction towards the cluster centres. The high AGN fraction in the outskirts is caused by low-luminosity AGNs. It can be explained by a higher galaxy merging rate in low-mass clusters, where velocity dispersions are not high enough to prevent galaxy interactions and merging. Ram pressure stripping is possible in the cores of all our clusters, but probably stronger in deeper gravitational potentials. Compared with previous studies of massive or high-redshift clusters, we conclude that the AGN fraction in cluster galaxies anti-correlates strongly with cluster mass. The AGN fraction also increases with redshift, but at the same rate with the respective fraction in field galaxies.

Context. The XMM-XXL survey uses observations from XMM-Newton to detect clusters of galaxies over a wide range in mass and redshift. The moderate PSF of XMM-Newton means that point sources within or projected onto a cluster may not be separated from the cluster emission, leading to enhanced luminosities and affecting the selection function of the cluster survey. Aims. We present the results of short Chandra observations of 21 galaxy clusters and cluster candidates at redshifts z $>$ 1 detected in the XMM-XXL survey in X-rays or selected in the optical and infra-red. Methods. With the superior angular resolution of Chandra we investigate whether there are any point sources within the cluster region that were not detected by the XMM-XXL analysis pipeline, and whether any point sources were misclassified as distant clusters. Results. Of the 14 X-ray selected clusters, nine are free from significant point source contamination, either having no previously unresolved sources detected by Chandra, or with less than about 10% of the reported XXL cluster flux being resolved into point sources. Of the other five sources, one is significantly contaminated by previously unresolved AGN, and four appear to be AGN misclassified as clusters. All but one of these cases are in the subset of less secure X-ray selected cluster detections and the false positive rate is consistent with that expected from the XXL selection function modelling. We also considered a further seven optically-selected cluster candidates associated with faint XXL sources that were not classed as clusters. By decontaminating and vetting these distant clusters, we provide a pure sample of clusters at redshift z $>$ 1 for deeper follow-up observations, and demonstrate the utility of using Chandra snapshots to test for AGN in surveys with high sensitivity but poor angular resolution.

Galaxy clusters trace the highest density peaks in the large-scale structure of the Universe. Their clustering provides a powerful probe that can be exploited in combination with cluster mass measurements to strengthen the cosmological constraints provided by cluster number counts. We investigate the spatial properties of a homogeneous sample of X-ray selected galaxy clusters from the XXL survey, the largest programme carried out by the XMM-Newton satellite. The measurements are compared to $\Lambda$-cold dark matter predictions, and used in combination with self-calibrated mass scaling relations to constrain the effective bias of the sample, $b_{eff}$, and the matter density contrast, $\Omega_{\rm M}$. We measured the angle-averaged two-point correlation function of the XXL cluster sample. The analysed catalogue consists of $182$ X-ray selected clusters from the XXL second data release, with median redshift $\langle z \rangle=0.317$ and median mass $\langle M_{500} \rangle\simeq1.3\cdot10^{14} M_\odot$. A Markov chain Monte Carlo analysis is performed to extract cosmological constraints using a likelihood function constructed to be independent of the cluster selection function. Modelling the redshift-space clustering in the scale range $10<r\,[$\Mpch$]<40$, we obtain $\Omega_{\rm M}=0.27_{-0.04}^{+0.06}$ and $b_{eff}=2.73_{-0.20}^{+0.18}$. This is the first time the two-point correlation function of an X-ray selected cluster catalogue at such relatively high redshifts and low masses has been measured. The XXL cluster clustering appears fully consistent with standard cosmological predictions. The analysis presented in this work demonstrates the feasibility of a cosmological exploitation of the XXL cluster clustering, paving the way for a combined analysis of XXL cluster number counts and clustering.

The XMM-XXL Survey spans two fields of $\rm 25$ deg$^2$ each observed for more than 6Ms with XMM, which provided a sample of tens of thousands of point sources with a flux limit of $\sim 2.2 \times 10^{-15}$ and $\sim 1.4 \times 10^{-14}$ erg s$^{-1}$ cm$^{2}$, corresponding to 50% of the area curve, in the soft band and hard band, respectively. In this paper we present the spatial clustering properties of $\sim 3100$ and $\sim 1900$ X-ray active galactic nuclei (AGNs) in the 0.5-2 and 2-10 keV bands, respectively, which have been spectroscopically observed with the AAOmega facility. This sample is 90% redshift complete down to an optical magnitude limit of $r\lesssim 21.8$. The sources span the redshift interval $0<z<5.2$, although in the current analysis we limit our samples to $z\le 3$, with corresponding sample median values of $\bar{z}\simeq 0.96$ and 0.79 for the soft band and hard band, respectively. We employ the projected two-point correlation function to infer the spatial clustering and find a correlation length $r_0=7.0 (\pm 0.34)$ and 6.42$(\pm 0.42)$ $h^{-1}$ Mpc, respectively, for the soft- and hard-band detected sources with a slope for both cases of $\gamma=1.44 (\pm 0.1)$. The power-law clustering was detected within comoving separations of 1 and $\sim$25 $h^{-1}$ Mpc. These results, as well as those derived in two separate redshift ranges, provide bias factors of the corresponding AGN host dark matter halos that are consistent with a halo mass of $\log_{10} [M_h/(h^{-1} M_{\odot})]=13.04\pm 0.06$, confirming the results of most recent studies based on smaller X-ray AGN samples.

Following extensive tests on analytical toy-catalogues (paper III), we present the results of a more realistic study over a 711 deg2 template-based cosmological simulation. Dark matter halos from the Aardvark simulation have been ascribed luminosities, temperatures and core-radii, using local scaling relations and assuming self-similar evolution. Predicted X-ray sky-maps are converted into XMM event lists using an instrumental simulator. The XXL pipeline run on the resulting sky images, produces an 'observed' cluster catalogue over which the tests have been performed. This allowed us to investigate the relative power of various combinations of count-rate, hardness-ratio, apparent-size and redshift information. Two fitting methods were used : a traditional MCMC approach and a simple minimisation procedure (Amoeba) whose mean uncertainties are a posteriori evaluated by means of synthetic catalogues. Results were analysed and compared to the predictions from a Fisher analysis (FA). For this particular catalogue realisation, assuming perfectly known scaling relations, the CR-HR combination gives $\sigma_8$ and $\Omega_m$ at the 10% level, while CR-HR-r$_c$-z improves this to $\leq$ 3%. Adding a second hardness ratio improves results from the CR-HR1-r$_c$ combination, but to a lesser extent than adding redshift information. When all coefficients of the M-T relation (including scatter) are also fitted, cosmological parameters are constrained to within 5-10%, and larger for the M-T coefficients (up to a factor of two for the scatter). Errors returned by the MCMC, by Amoeba, and by FA predictions, are in most cases in excellent agreement and always within a factor of two. We also study the impact of the scatter of the M-Rc relation on the number of detected clusters: for the cluster typical sizes usually assumed, the larger the scatter, the lower the number of detected objects.

Modern cosmological simulations rely heavily on feedback from active galactic nuclei (AGN) in order to stave off overcooling in massive galaxies and galaxy groups and clusters. An important independent test is whether or not the simulations capture the broad demographics of the observed AGN population. Here, we have used the cosmo-OWLS suite of cosmological hydrodynamical simulations to produce realistic synthetic catalogs of X-ray AGN out to $z$=3, with the aim of comparing the catalogs to the observed X-ray AGN population in the XXL survey and other recent surveys. We focused on the unabsorbed X-ray luminosity function (XLF), the Eddington ratio distribution, the black hole mass function, and the projected clustering of X-ray AGN. To compute the unabsorbed XLF of the simulated AGN, we used recent empirically-determined bolometric corrections. We show that the simulated AGN sample accurately reproduces the observed XLF over 3 orders of magnitude in X-ray luminosity in all redshift bins. To compare to the observed Eddington ratio distribution and the clustering of AGN, we produced detailed 'XMM-Newton-detected' catalogs of the simulated AGN. This requires the production of synthetic X-ray images extracted from light cones of the simulations that fold in the relevant instrumental effects of XMM-Newton. We apply a luminosity- and redshift-dependent obscuration function for the AGN and employ the same AGN detection algorithm as used for the real XXL survey. We demonstrate that the detected population of simulated AGN reproduces the observed Eddington ratio distribution and projected clustering from XXL quite well. We conclude that the simulations have a broadly realistic population of AGN and that our synthetic X-ray AGN catalogs should be useful for interpreting additional trends and as a helpful tool for quantifying AGN contamination in galaxy group and cluster X-ray surveys.

Despite compelling theoretical arguments, the use of clusters as cosmological probes is, in practice, frequently questioned because of the many uncertainties impinging on cluster mass estimates. Our aim is to develop a fully self-consistent cosmological approach of X-ray cluster surveys, exclusively based on observable quantities, rather than masses. This procedure is justified given the possibility to directly derive the cluster properties via ab initio modelling, either analytically or by using hydrodynamical simulations. In this third paper, we evaluate the method on cluster toy-catalogues. We model the population of detected clusters in the count-rate -- hardness-ratio -- angular size -- redshift space and compare the corresponding 4-dimensional diagram with theoretical predictions. The best cosmology+physics parameter configuration is determined using a simple minimisation procedure; errors on the parameters are derived by scanning the likelihood hyper-surfaces with a wide range of starting values. The method allows a simultaneous fit of the cosmological parameters, of the cluster evolutionary physics and of the selection effects. When using information from the X-ray survey alone plus redshifts, this approach is shown to be as accurate as the mass function for the cosmological parameters and to perform better for the cluster physics, as modelled in the scaling relations. It enables the identification of degenerate combinations of parameter values. Given the considerably shorter computer times involved for running the minimisation procedure in the observed parameter space, this method appears to clearly outperform traditional mass-based approaches when X-ray survey data alone are available.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

The XMM Cluster Archive Super Survey (X-CLASS) is a serendipitously-detected X-ray-selected sample of 845 galaxy clusters based on 2774 XMM archival observations and covering approximately 90 deg$^2$ spread across the high-Galactic latitude ($|b|>20$ deg) sky. The primary goal of this survey is to produce a well-selected sample of galaxy clusters on which cosmological analyses can be performed. This article presents the photometric redshift followup of a high signal-to-noise subset of 266 of these clusters with declination $\delta<+20$ deg with GROND, a seven channel ($grizJHK$) simultaneous imager on the MPG 2.2m telescope at the ESO La Silla Observatory. We use a newly developed technique based on the red sequence colour-redshift relation, enhanced with information coming from the X-ray detection to provide photometric redshifts for this sample. We determine photometric redshifts for 236 clusters, finding a median redshift of $z=0.39$ with an accuracy of $\Delta z = 0.02 (1+z)$ when compared to a sample of 76 spectroscopically confirmed clusters. We also compute X-ray luminosities for the entire sample and find a median bolometric luminosity of $7.2\times10^{43} \mathrm{erg\ s^{-1}}$ and a median temperature 2.9 keV. We compare our results to the XMM-XCS and XMM-XXL surveys, finding good agreement in both samples. The X-CLASS catalogue is available online at http://xmm-lss.in2p3.fr:8080/l4sdb/.

The growth of brightest cluster galaxies is closely related to the properties of their host cluster. We present evidence for dry mergers as the dominant source of BCG mass growth at $z\lesssim1$ in the XXL 100 brightest cluster sample. We use the global red sequence, H$\alpha$ emission and mean star formation history to show that BCGs in the sample possess star formation levels comparable to field ellipticals of similar stellar mass and redshift. XXL 100 brightest clusters are less massive on average than those in other X-ray selected samples such as LoCuSS or HIFLUGCS. Few clusters in the sample display high central gas concentration, rendering inefficient the growth of BCGs via star formation resulting from the accretion of cool gas. Using measures of the relaxation state of their host clusters, we show that BCGs grow as relaxation proceeds. We find that the BCG stellar mass corresponds to a relatively constant fraction 1\%\ of the total cluster mass in relaxed systems. We also show that, following a cluster scale merger event, the BCG stellar mass lags behind the expected value from the M$_{cluster}$ - M$_{BCG}$ relation but subsequently accretes stellar mass via dry mergers as the BCG and cluster evolve towards a relaxed state.

X-ray extragalactic surveys are ideal laboratories for the study of the evolution and clustering of active galactic nuclei (AGN). The XXL Survey spans two fields of a combined 50 $deg^2$ observed for more than 6Ms with XMM-Newton, occupying the parameter space between deep surveys and very wide area surveys; at the same time it benefits from a wealth of ancillary data. This paper marks the first release of the XXL point source catalogue selected in the 2-10 keV energy band with limiting flux $F_{2-10keV}=4.8\cdot10^{-14}\rm{erg\,s^{-1}\,cm^{-2}}$. We use both public and proprietary data sets to identify the counterparts of the X-ray point-like sources and improved upon the photometric redshift determination for AGN by applying a Random Forest classification trained to identify for each object the optimal photometric redshift model library. We also assign a probability to each source to be a star or an outlier. We model with Bayesian analysis the X-ray spectra assuming a power-law model with the presence of an absorbing medium. We find an average unabsorbed photon index of $\Gamma=1.85$ and average hydrogen column density $\log{N_{H}}=21.07 cm^{-2}$. We find no trend of $\Gamma$ or $N_H$ with redshift and a fraction of 26% absorbed sources ($\log N_{H}>22$). We show that the XXL-1000-AGN number counts extended the number counts of the COSMOS survey to higher fluxes and are fully consistent with the Euclidean expectation. We constrain the intrinsic luminosity function of AGN in the 2-10 keV energy band where the unabsorbed X-ray flux is estimated from the X-ray spectral fit up to z=3. Finally, we demonstrate the presence of a supercluster size structure at redshift 0.14, identified by means of percolation analysis of the XXL-1000-AGN sample. The XXL survey, reaching a medium flux limit and covering a wide area is a stepping stone between current deep fields and planned wide area surveys.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

Context. The XXL Survey is the largest survey carried out by the XMM-Newton satellite and covers a total area of 50 square degrees distributed over two fields. It primarily aims at investigating the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Aims. This article presents the XXL bright cluster sample, a subsample of 100 galaxy clusters selected from the full XXL catalogue by setting a lower limit of $3\times 10^{-14}\,\mathrm{erg \,s^{-1}cm^{-2}}$ on the source flux within a 1$^{\prime}$ aperture. Methods. The selection function was estimated using a mixture of Monte Carlo simulations and analytical recipes that closely reproduce the source selection process. An extensive spectroscopic follow-up provided redshifts for 97 of the 100 clusters. We derived accurate X-ray parameters for all the sources. Scaling relations were self-consistently derived from the same sample in other publications of the series. On this basis, we study the number density, luminosity function, and spatial distribution of the sample. Results. The bright cluster sample consists of systems with masses between $M_{500}=7\times 10^{13}$ and $3\times 10^{14} M_\odot$, mostly located between $z=0.1$ and 0.5. The observed sky density of clusters is slightly below the predictions from the WMAP9 model, and significantly below the predictions from the Planck 2015 cosmology. In general, within the current uncertainties of the cluster mass calibration, models with higher values of $\sigma_8$ and/or $\Omega_m$ appear more difficult to accommodate. We provide tight constraints on the cluster differential luminosity function and find no hint of evolution out to $z\sim1$. We also find strong evidence for the presence of large-scale structures in the XXL bright cluster sample and identify five new superclusters.

Traditionally, galaxy clusters have been expected to retain all the material accreted since their formation epoch. For this reason, their matter content should be representative of the Universe as a whole, and thus their baryon fraction should be close to the Universal baryon fraction. We make use of the sample of the 100 brightest galaxy clusters discovered in the XXL Survey to investigate the fraction of baryons in the form of hot gas and stars in the cluster population. We measure the gas masses of the detected halos and use a mass--temperature relation directly calibrated using weak-lensing measurements for a subset of XXL clusters to estimate the halo mass. We find that the weak-lensing calibrated gas fraction of XXL-100-GC clusters is substantially lower than was found in previous studies using hydrostatic masses. Our best-fit relation between gas fraction and mass reads $f_{\rm gas,500}=0.055_{-0.006}^{+0.007}\left(M_{\rm 500}/10^{14}M_\odot\right)^{0.21_{-0.10}^{+0.11}}$. The baryon budget of galaxy clusters therefore falls short of the Universal baryon fraction by about a factor of two at $r_{\rm 500}$. Our measurements require a hydrostatic bias $1-b=M_X/M_{\rm WL}=0.72_{-0.07}^{+0.08}$ to match the gas fraction obtained using lensing and hydrostatic equilibrium. Comparing our gas fraction measurements with the expectations from numerical simulations, our results favour an extreme feedback scheme in which a significant fraction of the baryons are expelled from the cores of halos. This model is, however, in contrast with the thermodynamical properties of observed halos, which might suggest that weak-lensing masses are overestimated. We note that a mass bias $1-b=0.58$ as required to reconcile Planck CMB and cluster counts should translate into an even lower baryon fraction, which poses a major challenge to our current understanding of galaxy clusters. [Abridged]

NICMOS 2 observations are crucial for constraining distances to most of the existing sample of z > 1 SNe Ia. Unlike the conventional calibration programs, these observations involve long exposure times and low count rates. Reciprocity failure is known to exist in HgCdTe devices and a correction for this effect has already been implemented for high and medium count-rates. However observations at faint count-rates rely on extrapolations. Here instead, we provide a new zeropoint calibration directly applicable to faint sources. This is obtained via inter-calibration of NIC2 F110W/F160W with WFC3 in the low count-rate regime using z ~ 1 elliptical galaxies as tertiary calibrators. These objects have relatively simple near-IR SEDs, uniform colors, and their extended nature gives superior signal-to-noise at the same count rate than would stars. The use of extended objects also allows greater tolerances on PSF profiles. We find ST magnitude zeropoints (after the installation of the NICMOS cooling system, NCS) of 25.296 +- 0.022 for F110W and 25.803 +- 0.023 for F160W, both in agreement with the calibration extrapolated from count-rates 1,000 times larger (25.262 and 25.799). Before the installation of the NCS, we find 24.843 +- 0.025 for F110W and 25.498 +- 0.021 for F160W, also in agreement with the high-count-rate calibration (24.815 and 25.470). We also check the standard bandpasses of WFC3 and NICMOS 2 using a range of stars and galaxies at different colors and find mild tension for WFC3, limiting the accuracy of the zeropoints. To avoid human bias, our cross-calibration was "blinded" in that the fitted zeropoint differences were hidden until the analysis was finalized.

We performed a detailed and, for a large part interactive, analysis of the matching output between the X-CLASS and redMaPPer cluster catalogues. The overlap between the two catalogues has been accurately determined and possible cluster positional errors were manually recovered. The final samples comprise 270 and 355 redMaPPer and X-CLASS clusters respectively. X-ray cluster matching rates were analysed as a function of optical richness. In a second step, the redMaPPer clusters were correlated with the entire X-ray catalogue, containing point and uncharacterised sources (down to a few 10^{-15} erg s^{-1} cm^{-2} in the [0.5-2] keV band). A stacking analysis was performed for the remaining undetected optical clusters. Main results show that neither of the wavebands misses any massive cluster (as coded by X-ray luminosity or optical richness). After correcting for obvious pipeline short-comings (about 10% of the cases both in optical and X-ray), ~50% of the redMaPPer (down to a richness of 20) are found to coincide with an X-CLASS cluster; when considering X-ray sources of any type, this fraction increases to ~ 80%; for the remaining objects, the stacking analysis finds a weak signal within 0.5 Mpc around the cluster optical centers. The fraction of clusters totally dominated by AGN-type emission appears to be of the order of a few percent. Conversely ~ 40% of the X-CLASS clusters are identified with a redMaPPer (down to a richness of 20) - part of the non-matches being due to the fact that the X-CLASS sample extends further out than redMaPPer (z<1 vs z<0.6); extending the correlation down to a richness of 5, raises the matching rate to ~ 65%.

We present a catalog of RR Lyrae stars (RRLs) observed by the Xuyi Schmidt Telescope Photometric Survey (XDSS). The area we consider is located in the North Galactic Cap, covering 376.75 sq deg at RA $\approx$ 150 deg and Dec $\approx$ 27 deg down to a magnitude limit of i $\approx$ 19. Using the variability information afforded by the multi-epoch nature of our XDSS data, combined with colors from the Sloan Digital Sky Survey, we are able to identify candidate RRLs. We find 318 candidates, derive distances to them and estimate the detection efficiency. The majority of our candidates have more than 12 observations and for these we are able to calculate periods. These also allows us to estimate our contamination level, which we predict is between 30% to 40%. Finally we use the sample to probe the halo density profile in the 9-49 kpc range and find that it can be well fitted by a double power law. We find good agreement between this model and the models derived for the South Galactic Cap using the Watkins et al. (2009) and Sesar et al. (2010) RRL data-sets, after accounting for possible contamination in our data-set from Sagittarius stream members. We consider non-spherical double power law models of the halo density profile and again find agreement with literature data-sets, although we have limited power to constrain the flattening due to our small survey area. Much tighter constraints will be placed by current and future wide-area surveys, most notably ESA's astrometric Gaia mission. Our analysis demonstrates that surveys with a limited number of epochs can effectively be mined for RRLs. Our complete sample is provided as accompanying online material.

We present the final release of the multi-wavelength XMM-LSS data set,covering the full survey area of 11.1 square degrees, with X-ray data processed with the latest XMM-LSS pipeline version. The present publication supersedes the Pierre et al.(2007) catalogue pertaining to the initial 5 square degrees. We provide X-ray source lists in the customary energy bands (0.5-2 and 2-10 keV) for a total of 6721 objects in the deep full-exposure catalogue and 5572 in the 10ks-limited one, above a detection likelihood of 15 in at least one band. We also provide a multiwavelength catalogue, cross-correlating our list with IR, NIR, optical and UV catalogues. Customary data products (X-ray FITS images, CFHTLS and SWIRE thumbnail images) are made available together with our interactively queriable database in Milan, while a static snapshot of the catalogues will be supplied to CDS, as soon as final acceptance is completed.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

We study the effect of filter zero-point uncertainties on future supernova dark energy missions. Fitting for calibration parameters using simultaneous analysis of all Type Ia supernova standard candles achieves a significant improvement over more traditional fit methods. This conclusion is robust under diverse experimental configurations (number of observed supernovae, maximum survey redshift, inclusion of additional systematics). This approach to supernova fitting considerably eases otherwise stringent mission calibration requirements. As an example we simulate a space-based mission based on the proposed JDEM satellite; however the method and conclusions are general and valid for any future supernova dark energy mission, ground or space-based.

We report on work to increase the number of well-measured Type Ia supernovae (SNe Ia) at high redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe Ia that were discovered during 2001 are presented. Additionally, for the two SNe with z>1, we present ground-based J-band photometry from Gemini and the VLT. These are among the most distant SNe Ia for which ground based near-IR observations have been obtained. We add these six SNe Ia together with other data sets that have recently become available in the literature to the Union compilation (Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the most important ones being the refitting of all light curves with the SALT2 fitter and an improved handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 compilation. The flat concordance LambdaCDM model remains an excellent fit to the Union2 data with the best fit constant equation of state parameter w=-0.997^{+0.050}_{-0.054} (stat) ^{+0.077}_{-0.082} (stat+sys\ together) for a flat universe, or w=-1.035^{+0.055}_{-0.059} (stat)^{+0.093}_{-0.097} (stat+sys together) with curvature. We also present improved constraints on w(z). While no significant change in w with redshift is detected, there is still considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data.

Dome A, the highest plateau in Antarctica, is being developed as a site for an astronomical observatory. The planned telescopes and instrumentation and the unique site characteristics are conducive toward Type Ia supernova surveys for cosmology. A self-contained search and survey over five years can yield a spectro-photometric time series of ~1000 z<0.08 supernovae. These can serve to anchor the Hubble diagram and quantify the relationship between luminosities and heterogeneities within the Type Ia supernova class, reducing systematics. Larger aperture (>4-m) telescopes are capable of discovering supernovae shortly after explosion out to z~3. These can be fed to space telescopes, and can isolate systematics and extend the redshift range over which we measure the expansion history of the universe.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

We present data from the Supernova Cosmology Project for five high redshift Type Ia supernovae (SNe Ia) that were obtained using the NICMOS infrared camera on the Hubble Space Telescope. We add two SNe from this sample to a rest-frame I-band Hubble diagram, doubling the number of high redshift supernovae on this diagram. This I-band Hubble diagram is consistent with a flat universe (Omega_Matter, Omega_Lambda= 0.29, 0.71). A homogeneous distribution of large grain dust in the intergalactic medium (replenishing dust) is incompatible with the data and is excluded at the 5 sigma confidence level, if the SN host galaxy reddening is corrected assuming R_V=1.75. We use both optical and infrared observations to compare photometric properties of distant SNe Ia with those of nearby objects. We find generally good agreement with the expected color evolution for all SNe except the highest redshift SN in our sample (SN 1997ek at z=0.863) which shows a peculiar color behavior. We also present spectra obtained from ground based telescopes for type identification and determination of redshift.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

We present observations of SCP 06F6, an unusual optical transient discovered during the Hubble Space Telescope Cluster Supernova Survey. The transient brightened over a period of ~100 days, reached a peak magnitude of ~21.0 in both i_775 and z_850, and then declined over a similar timescale. There is no host galaxy or progenitor star detected at the location of the transient to a 3 sigma upper limit of i_775 = 26.4 and z_850 = 26.1, giving a corresponding lower limit on the flux increase of a factor of ~120. Multiple spectra show five broad absorption bands between 4100 AA and 6500 AA and a mostly featureless continuum longward of 6500 AA. The shape of the lightcurve is inconsistent with microlensing. The transient's spectrum, in addition to being inconsistent with all known supernova types, is not matched to any spectrum in the Sloan Digital Sky Survey (SDSS) database. We suggest that the transient may be one of a new class.

The recent robust and homogeneous analysis of the world's supernova distance-redshift data, together with cosmic microwave background and baryon acoustic oscillation data, provides a powerful tool for constraining cosmological models. Here we examine particular classes of scalar field, modified gravity, and phenomenological models to assess whether they are consistent with observations even when their behavior deviates from the cosmological constant Lambda. Some models have tension with the data, while others survive only by approaching the cosmological constant, and a couple are statistically favored over LCDM. Dark energy described by two equation of state parameters has considerable phase space to avoid Lambda and next generation data will be required to constrain such physics.

We present results from a study of the orbits of eclipsing binary stars (EBs) in the Magellanic Clouds. The samples comprise 4510 EBs found in the Large Magellanic Cloud (LMC) by the MACHO project, 2474 LMC EBs found by the OGLE-II project (of which 1182 are also in the MACHO sample), 1380 in the Small Magellanic Cloud (SMC) found by the MACHO project, and 1317 SMC EBs found by the OGLE-II project (of which 677 are also in the MACHO sample); we also consider the EROS sample of 79 EBs in the bar of the LMC. Statistics of the phase differences between primary and secondary minima allow us to infer the statistics of orbital eccentricities within these samples. We confirm the well-known absence of eccentric orbit in close binary stars. We also find evidence for rapid circularization in longer period systems when one member evolves beyond the main sequence, as also found by previous studies.

We present a new sample of 4634 eclipsing binary stars in the Large Magellanic Cloud (LMC), expanding on a previous sample of 611 objects and a new sample of 1509 eclipsing binary stars in the Small Magellanic Cloud (SMC), that were identified in the light curve database of the MACHO project. We perform a cross correlation with the OGLE-II LMC sample, finding 1236 matches. A cross correlation with the OGLE-II SMC sample finds 698 matches. We then compare the LMC subsamples corresponding to center and the periphery of the LMC and find only minor differences between the two populations. These samples are sufficiently large and complete that statistical studies of the binary star populations are possible.

The light curve of a blue variable in the MACHO LMC database (FTS ID 78.5979.72) appeared nearly unvarying for ~4 years (quasi-flat segment) but then rapidly changed to become periodic with noisy minima for the remaining 4 years (periodic segment); there are no antecedent indications of a gradual approach to this change. Lomb Periodogram analyses indicate the presence of two distinct periods of ~61 days and 8 days in both the quasi-flat and the periodic segments. Minima of the periodic segment cover at least 50% of the orbital period and contain spikes of light with the 8-day period; maxima do not show this short period. The system typically shows maxima to be redder than minima. The most recent OGLE-III light curve shows only a 30-day periodicity. The variable's V and R magnitudes and color are those of a Be star, and recent sets of near infrared spectra four days apart, secured during the time of the OGLE-III data, show H-alpha emission near and at a maximum, confirming its Be star characteristics. The model that best fits the photometric behavior consists of a thin ring-like circumstellar disk of low mass with four obscuring sectors orbiting the central B star in unison at the 61-day period. The central star peers through the three equi- spaced separations between the four sectors producing the 8-day period. The remainder of the disk contains hydrogen in emission making maxima appear redder. A companion star of lower mass in an inclined and highly eccentric orbit produces an impulsive perturbation near its periastron to change the disk's orientation, changing eclipses from partial to complete within ~ 10 days.

We present a methodology to discover outliers in catalogs of periodic light-curves. We use cross-correlation as measure of ``similarity'' between two individual light-curves and then classify light-curves with lowest average ``similarity'' as outliers. We performed the analysis on catalogs of variable stars of known type from the MACHO and OGLE projects and established that our method correctly identifies light-curves that do not belong to those catalogs as outliers. We show how our method can scale to large datasets that will be available in the near future such as those anticipated from Pan-STARRS and LSST.