Optical Follow-up of Planck Cluster Candidates with ... - IOPscience

0 downloads 0 Views 1MB Size Report
Aug 9, 2018 - We report on the search for optical counterparts of Planck Sunyaev-Zel'dovich ... large-scale structure of universe – galaxies: clusters: general – (cosmology:) .... Although the primary mirror is of 0.62 m diameter, the holding.

Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October


© 2018. The Astronomical Society of the Pacific. Printed in the U.S.A.

Optical Follow-up of Planck Cluster Candidates with Small Instruments Vincent Boucher1, Simon de Visscher2, and Christophe Ringeval3 1

CREEF asbl, 3 rue d’Ostin, 5080 La Bruyère, Belgium 4 rue du collège, 1407 Mézery-près-Donneloye, Switzerland 3 Centre for Cosmology, Particle Physics and Phenomenology, Institute of Mathematics and Physics, Louvain University, 2 chemin du cyclotron, 1348 Louvain-la-Neuve, Belgium; [email protected] Received 2018 June 4; accepted 2018 July 23; published 2018 August 9 2

Abstract We report on the search for optical counterparts of Planck Sunyaev-Zel’dovich (SZ) cluster candidates using a 0.6 m non-professional telescope. Among the observed sources, an unconfirmed candidate, PSZ2 G156.24 +22.32, is found to be associated with a region of more than 100 galaxies within a 3 arcminute radius around the Sunyaev-Zel’dovich maximum signal coordinates. Using 14 hours of cumulated exposure over the Sloan color filters g′, r′, i′, and z′, we estimate the photometric redshift of these galaxies at zphot=0.29±0.08. Using the red+0.03 sequence galaxy method gives a photometric redshift of 0.300.05 . Combined with the Planck SZ proxy mass 14 function, this would favor a cluster of 4.4×10 solar masses. This result suggests that a dedicated pool of observatories equipped with such instruments could collectively contribute to optical follow-up programs of massive cluster candidates at moderate redshifts. Key words: telescopes – (cosmology:) large-scale structure of universe – galaxies: clusters: general – (cosmology:) cosmic background radiation – catalogs Online material: color figures PSZ2 catalog in Planck Collaboration et al. (2016a), containing more than 1600 SZ sources. Among them, 1200 objects have been confirmed as clusters by the Planck collaboration through their cross identification in the Meta-Catalog of X-ray detected clusters (MCXC; Piffaretti et al. 2011), with optical counterparts in the Sloan Digital Sky Survey (SDSS; York et al. 2000), in the redMaPPer catalog(Rozo et al. 2015; Rykoff et al. 2016), in the Nasa/IPAC Extragalactic Database (NED), and with infrared galaxy overdensities in the Wide-field Infrared Survey (WISE; Wright et al. 2010). More than 400 SZ sources were still unconfirmed at the time of publication of the catalog, and not all confirmed sources have redshift information. Cluster counts using the PSZ2 catalog have been used for Cosmology in Planck Collaboration et al. (2016b). They allow CMB-independent estimations of the amplitude of the matter power spectrum σ8 and of the matter density parameter Ωm, which are currently in mild tension with the best fit ΛCDM model obtained from CMB primary anistropies alone(Planck Collaboration et al. 2016c). Doing cosmology with cluster counts requires the determination of a scaling relation between the total integrated Compton parameter Y5R500, the cluster angular size, θs, and its total mass, M500 (defined over a radius enclosing 500 times the critical density at redshift z). As explained in Planck Collaboration et al. (2014b), SZ measurements give information on the relation between Y5R500 and θs while breaking the degeneracy between Y5R500 and M500, at a given redshift z, relies on a scaling relation extracted from X-ray observations and assuming hydrostatic

1. Introduction Among the foregrounds imprinting temperature anisotropies and spectral distortions on the cosmic microwave background (CMB) radiation, the thermal Sunyaev-Zel’dovich effect is of particular interest for astrophysics and cosmology (Sunyaev & Zeldovich 1970, 1972, 1980; Sazonov & Sunyaev 1999; Challinor et al. 2000; Bunn 2006). Inverse Compton scattering of the CMB photons by the hot electrons present within galaxy clusters results in a shift and distortion of their blackbody spectrum towards higher frequencies. The intracluster plasma cools down the CMB radiation at frequencies typically lower than 217 GHz, while warming it up at higher frequencies. Such a signature is unique among foregrounds and has been intensively used by the Planck satellite collaboration and other CMB ground telescopes, such as the Atacama Cosmology Telescope (ACT) and the South Pole Telescope, to provide unprecedented catalogs of Sunyaev-Zel’dovich (SZ) sources (Hasselfield et al. 2013; Planck Collaboration et al. 2014a; Bleem et al. 2015; Planck Collaboration et al. 2015a). The SZ effect does not depend on the redshift, and therefore provides a new high-redshift observable of galaxy clusters. The full sky coverage of the Planck satellite has allowed the release of the Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October equilibrium of the intracluster gas (Arnaud et al. 2010, 2007). For each of the SZ sources, the PSZ2 catalog provides the most probable hydrostatic mass MSZ(z) (and its standard deviation) assuming the scaling relation to hold. Follow-up programs of SZ sources are therefore of immediate interest in evaluating MSZ for a given cluster by the determination of its redshift z. Dedicated optical follow-up programs of unconfirmed Planck clusters have been carried on using professional telescopes since the publication of the PSZ2 catalog. The Pan-STARRS 1.8 m telescope survey(Liu et al. 2015) has provided 60 confirmations combined with spectroscopic redshift measurements. Another 16 high-redshift clusters have been confirmed by the Canada-France-Hawaï telescope (CFHT) with photometric redshifts together with richness mass estimates(van der Burg et al. 2016). Spectroscopic redshifts and confirmation of 13 more clusters have been provided by the 1.5 m Russian-Turkish telescope, the 2.2 m Calar Alto Observatory telescope, and the 6 m Bolshoi (BTA) telescope in Vorobyev et al. (2016; see also Planck Collaboration et al. 2015b). The Planck collaboration has used a month of observations at the Canary Islands Observatories for the confirmation of 73 more candidates(Planck Collaboration et al. 2016d). Although these numbers could suggest that unconfirmed SZ sources would soon be exhausted, artificial neural networks have recently been used on the Planck CMB measurements to improve the detection threshold of SZ sources. A new catalog discussed by Hurier et al. (2017) contains almost 4000 galaxy cluster candidates. In addition, ground-based telescopes are still contributing to new SZ detection, thereby calling for an increase of the telescope time needed for follow-up programs(Hilton et al. 2017). In this paper, we explore and report on the possibility to use sub-meter, non-professional telescopes (of typically 0.5 m– 0.6 m diameter) to carry on optical follow-up searches of unconfirmed SZ sources. In spite of obvious technical challenges, some scientific measurements can be achieved such as galaxy counts and photometric redshift estimates. As a test case, we report on a cluster candidate, PSZ2 G156.24 +22.32, which has not been confirmed in any of the aforementioned catalogs, nor in the follow-up programs, and which is not included in the SDSS sky coverage(Ahn et al. 2012; Wen et al. 2012). The paper is organized as follows. Section 2 describes the observations and the method used for data reduction. Some details are provided concerning the observatory, the hardware used, as well as the technical difficulties overcome for achieving sufficiently good tracking and imaging for photometry. In Section 3, we present our main results, namely galaxy counts, photometric redshift, and optical richness estimates for PSZ2 G156.24+22.32. In Section 4, various tests to assess the validity of our measurements are presented. In the conclusion, we briefly discuss a proof of concept on

Boucher, de Visscher, & Ringeval

what small telescopes can and cannot achieve, would they massively contribute to the search of optical counterparts of SZ sources.

2. Observations and Data Reduction 2.1. Observatory and Instrumentation The observatory of Saint-Véran is located in France, at 44° 42′03″ N, 6°52′06″ E, at an altitude of 2930 m a.s.l. (IAU code 615). It has been managed by the Astroqueyras Association4 for roughly 30 years, and is intended to allow amateur astronomers to observe under better conditions compared to that of typical lower altitude suburban skies. The observatory comprises three main instruments: a Cassegrain telescope with a diameter of 0.62 m and 9 m of focal length (T62), and two Ritchey-Chrétien telescopes, each having a diameter of 0.5 m and 4 m of focal length. The instrument mainly used in the context of this work is the T62. Although the primary mirror is of 0.62 m diameter, the holding baffle reduces the effective aperture to 0.60 m. The T62 is mounted on a German-type REOSC equatorial mount, equipped electronically for semi-automatic target pointing. The imaging device used is an APOGEE U16M CCD camera equipped with a KAF-16803 sensor, cooled at −20° C. The sensor has 4096×4096 pixels of area 9×9 μm2, which implies a genuine angular sampling of 0 21×0 21. The field of view is a square of 14′ side. The King (1989) polar drift method has been applied to align the mount right ascension (R.A.) and Earth rotation axes, and therefore improves the pointing and tracking accuracy. This operation must be done regularly, due to local seismic activity and a slow rotation of the mountain where the observatory is installed. In the context of this work, the residual angle between the right ascension axis of the mount and the Earth rotation axis was reduced to a value smaller than roughly 10″. The native periodical error of the RA movement is estimated to be around 5″ at vanishing declination (DEC), with smooth amplitude variations during tracking. The tracking accuracy is yet improved by auto-guiding the T62, i.e., by correcting any deviations between the sky and mount movements. This is achieved using a second camera installed on a 0.2 m SchmidtCassegrain telescope. It is mounted in parallel with the T62, and has a focal length of 2 m. The camera has a EXVIEW sensor made of 752×580 pixels of area 8.2×8.4 μm2. This allows for exposures up to 10 to 15 minutes long, depending on the telescope orientation and the target declination. The resulting guiding error is estimated from the average point-spread function of the stars and found to be less or equal to than the seeing. The measured FWHM varies between 1.2″ to 2″ across the exposure series. These FWHM values led us to optimize the sampling for the APOGEE U16M camera. A binning grouping 4



Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October

Boucher, de Visscher, & Ringeval

hemisphere), ended up showing a significant object count. We therefore started data acquisition in the color filter g′, r′, i′, and z′ from 2017 February 21 until the end of the mission on February 26. According to the above-mentioned criteria, PSZ2 G156.24+22.32 was of sufficiently high altitude in the sky only during the second half of the night. The first part of the night was therefore dedicated in the color imaging of our reference photometric SDSS DR9 galaxy field, located nearby M63, which is roughly at the same altitude in the sky than the one occupied by the cluster later on. Around 1.5 hours of cumulated exposure per filter were taken for the reference field. In total, around 14 hours of color exposure have been used to derive the results presented below.

four pixels into one has been used, yielding an angular sampling of 0 42×0 42. The search for optical counterparts of the SZ targets has been performed in luminance, i.e., without a specific spectral filtering of the light input. Once the galaxies have been located, color imaging of the suspected cluster is done using interferential filters. For this work, we have used a set of ASTRODON GEN 2 SLOAN photometric filters corresponding to the SDSS photometric standards(Fukugita et al. 1996; Smith et al. 2002). The low quantum efficiency of the KAF-16803 sensor in the ultraviolet wavelength range led us to consider only the g′, r′, i′ and z′ filters.

2.2. Data Acquisition 2.3. Reduction

The search of optical counterparts of the Planck SZ candidates took place during two one-week long observation missions, one from 2015 November 7 to November 15, and the other between 2017 February 18 to February 26. The goal of the first mission was to assess the feasibility of the project and was dedicated to the observations of confirmed Planck clusters having known redshifts. Individual galaxy identifications for five clusters having redshifts ranging from 0.1 to 0.6 were successfully performed. In particular, color imaging of two of them, PSZ2 G138.32−39.82 and PSZ2 G164.18−38.88, located at redshift 0.280 and 0.08, respectively, allowed us to test the instrumentation for photometric redshift estimation. A low-signal test of our apparatus was achieved by successfully resolving individual galaxies in luminance within a distant cluster PSZ2 G144.83 +25.11, having a redshift of 0.584 (also known as MACS J0647.7+7015). These results are discussed in more details in Section 4. The 2017 winter mission was dedicated to the search of nonconfirmed Planck SZ sources optical counterparts. Five objects from the PSZ2 catalog, having no X-ray nor optical counterparts, have been pointed. The choice of the targets has been made taking into account several constraints. First, we have kept the airmass reasonably close to unity and only objects well above the horizon by at least 40° have been considered. Second, the targets should not be too close from the Milky Way to minimize possible confusion between faint stars and distant galaxies. A technical difficulty associated with the use of a small telescope is the exposure time needed to collect enough photons for galaxy detection. Based on empirical results from the first mission, around 1 hour of exposure has been accumulated per target, made from 10 minutes luminance slices. After performing some minimal data reduction, a visual inspection allowed us to check for overdensity within the 14′ field of view around the SZ signal coordinates. Among the observed targets, PSZ2 G156.24+22.32 located in the Lynx constellation (northern celestial

To produce photometric usable data, image reduction has been performed in two steps. In a first pass, we have used the Image Reduction and Analysis Facility (IRAF) to correct each image from bias, dark currents and nonlinear response pixels (Tody 1986, 1993). Basic astrometric registration and airmass calculations are also performed for each color images to provide a world coordinate system (WCS) compatible header for each of them. The second pass consists in an accurate astrometry and photometry solving starting from the WCS single-epoch images produced by the first pass. For this purpose, we have used the ASTROMATIC software suite (Bertin & Arnouts 1996; Bertin 2006; Bertin et al. 2002). Sources are extracted in all images using SEXTRACTOR (version 2.19.5) with a detection and analysis threshold set at 1.5σ with respect to the local background noise. An astrometric and photometric solution is then obtained by using the SCAMP code. It has been run in multi-instruments mode, one astrometric and photometric instrument for each color filter, plus luminance, and we have allowed for astrometric distortions within a given instrument. The reference astrometry catalog chosen is the USNO-B1 (Monet et al. 2003). The reported internal astrometric errors are small, 0.05″ (FWHM), and remain negligible with respect to the ones of the reference catalog, evaluated at 0.5″ (FWHM). However, a shear contraction amplitude of typically 0.15% has been corrected for each instrument, far above what one could have expected for an airmass never exceeding 1.4 (0.05%). The distortion map of the pixel scale computed by SCAMP indeed reveals an almost constant hyperbolic pattern over all filters and luminance. Although the reason for this distortion is unclear, different sources are possible, among which a stress on the primary mirror, or still a minor mirror alignment issue leading to astigmatism appearing at moderate defocus. Time constraint did not allow for a deeper exploration and resolution of this issue. Concerning the photometric solution, zero-point changes across the exposure series have been corrected within each filter. The standard deviations of these corrections are found in


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October

Boucher, de Visscher, & Ringeval

the following ranges: Dg0¢ Î [0.1, 0.4], Dr0¢ Î [0.05, 0.3], Di0¢ Î [0.05, 0.4], and Dz 0¢ Î [0.05, 0.5]. The first (and lowest) value is for high signal-over-noise matched sources only while the second one includes all matched objects between the singleepoch images and the reference catalog. Although corrected, these numbers show some degradation of our photometric stability up to half a magnitude for faint objects in the z′ filter, and some significant scatter in the g′ band as well. This could be attributed to the parasitic presence of a V=7.2 magnitude double star in the field of view, as a much lower scatter is obtained for the SDSS reference field (see below and Figure 1). In the last step, coaddition of the corrected single-epoch images for each filter (and luminance) has been delegated to the SWARP code, which takes care of background estimation, as well as resampling, for each image. Highest contrasts have been obtained by using a median weighted stacking, the weights being estimated from the background noise(see Bertin et al. 2002).

2.4. Calibration The absolute zero-point magnitudes for the four coadded images in the g′, r′, i′, and z′ filter have been determined from our SDSS DR9 reference galaxy field (see Section 2.2). The reference images have been reduced and flattened, for each color filter, exactly as described above. Then, we have used the SEXTRACTOR code over the coadded images to create a catalog of photometric reference sources. The magnitude measurement method chosen is mag_auto and only high signal-over-noise detections have been kept, 5σ above the estimated background(Bertin & Arnouts 1996). Our reference catalogs, in g′, r′, i′, and z′, are then matched together and to the SDSS DR9 catalog by cross identifying objects of same position at less than 1″. For this purpose, the imaging software DS9 has been used to create a set of roughly 30 galaxies for which one has both the non-calibrated mag_auto magnitudes and the modelMag magnitudes as provided by the SDSS DR9 catalog(Ahn et al. 2012). The zero-point values for each filter, g′, r′, i′, and z′, have been determined with IRAF using a photometric weighting based on the magnitude error estimates made by SEXTRACTOR. In the following, unless specified otherwise, our calibrated mag_auto magnitudes will be referred to as griz since they are indistinguishable of the modelMag SDSS DR9 magnitudes within the uncertainties of our apparatus and methodology.

Figure 1. Upper panel: circular regions of radius 2′ and 3′ centered at the Planck SZ maximum signal location (marked by a cross) used for galaxy counts and photometric redshift estimates. Objects with a small circle are matched sources across the color filters (stars excluded; see the text). The lower panel is a zoom within the 3′ region. (A color version of this figure is available in the online journal.)

suspected to be galaxies within the 14′ field of view around the SZ coordinates associated with PSZ2 G156.24+22.32. Again, the SEXTRACTOR code has been used for this purpose, in single-image mode and with a detection threshold at 2σ for the g′, r′, and i′ filter. As discussed in the previous section, the z′

3. Results 3.1. Color Catalog of Matched Objects From the photometric usable images obtained as described in Section 2, we have constructed a griz catalog of all objects


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October

Boucher, de Visscher, & Ringeval

from the USNO-URAT1 catalog is unexpectedly doing a good job in removing bright and close galaxies visible in the picture. A closer look to Figure 1 also reveals that many small objects close to the bright star, and which are certainly galaxies belonging to the cluster, are not present in our griz catalog. This is due to the fact that these sources have been excluded by SEXTRACTOR as being too contaminated, in at least one color, by the parasitic glow of the star. As discussed below, five more objects close to the star end up being obvious outliers and will be removed from the subsequent analysis (drawn with a small white box). From the Planck catalog MMF3, we have extracted the twodimensional probability distribution in the plane (Y5R500, θs) associated with PSZ2 G156.24+22.32 (Planck Collaboration et al. 2016a). By marginalizing over the integrated Comptonization parameter Y5R500, the most probable angular size of the SZ region is found to be qs = 2.5¢  0.5¢

(1 )

at 68% of confidence. As a result, we have defined a circular region around the SZ maximum having a radius of 3′ to encompass the whole SZ region. It is represented as the yellow circle in Figure 1. To assess the sensitivity of the results with respect to the chosen cluster angular size, another circular region has been defined, centered at the SZ location, with a radius of 2′ (dashed inner circle in Figure 1). The intersection of our catalog with these two regions will be referred to as griz3 and griz2 in the following; they contain 114 and 58 objects, respectively (including the five outliers; see Figure 1). Finally, let us mention that PSZ2 G156.24+22.32 is associated in the PSZ2 catalog with a value of Q_NEURAL= 0.96, which indicates a high probability of being a genuine SZ source. Moreover, as reported by van der Burg et al. (2016), all the SZ candidates validated by optical counterparts in this work have a Q_NEURAL>0.8, thereby rendering our observation of galaxy overdensity around the SZ location not unlikely.

Figure 2. Photometric redshift estimates of the artificial neural network over an independent testing set of 250,000 galaxies from the SDSS DR9 catalog. The upper panel shows the neural net estimated redshift zphot as a function of photometric SDSS redshift (the plot has been truncated to 8000 random samples for illustration convenience). The lower panel shows the distribution of the residuals Δzphot, and has a standard deviation of sDzphot = 0.027. (A color version of this figure is available in the online journal.)

image showing higher noise than the others, SEXTRACTOR has been used in two-images mode, 2σ detection coming from the i′ image, photometry coming from the z′ one with an analysis threshold set at 1.5σ. Separation between stars and galaxies is made using the neural network classifier of SEXTRACTOR by discarding all objects having a CLASS_STAR greater than 0.98. The remaining objects over all filters are then matched together according to their position, required to be the same at 1″, and then rejected if they correspond to an optical counterpart in the USNO-URAT1 catalog(Zacharias et al. 2015). In Figure 1, the SZ signal center, as reported in the PSZ2 catalog, is represented by a cross, and small circles have been drawn around all the remaining objects of our griz catalog. A visual inspection of their shape suggests that many of them should be galaxies, while some overdensity is visible left from the SZ center. Notice that having chosen to exclude objects

3.2. Photometric Redshift and Cluster Mass To evaluate the photometric redshift of the objects belonging to our color catalogs, we have used an artificial neural network (ANN) regression of the SDSS DR9 color-redshift relationship, but only for the reduced set of griz colors. For this purpose, the publicly available SKYNET code has been run to train a simple feed-forward artificial neural network (Graff et al. 2012, 2014; Hobson et al. 2014). The ANN input layer takes as argument the griz modelMag magnitudes, and the output layer returns the photometric redshift. All input and hidden layer nodes of the ANN are evaluating the hyperbolic tangent of the biased weighted addition of the input signals. For faster training, the output layer has, however, been chosen to only perform a biased linear combination of the last hidden layer nodes. The


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October training set chosen consists of 250,000 photometric clean objects from the SDSS DR9 data release, selected to have nearest neighbors in the SDSS color space, photometric redshift errors not exceeding 0.03, and color magnitudes in the range [12, 29]. Another group of 250,000 objects, selected with the same criteria, has been used as a testing set to check the quality of the ANN regression. We have found that a good enough learning of the ANN over the training set requires four hidden layers of 90 nodes each, in addition to the input layer of four nodes and the single-node output layer. In total, SKYNET has been used to find the optimal values of more than 30,000 weights and biases as described in Graff et al. (2014). The upper panel of Figure 2 shows the ANN output values obtained by inputting the griz color magnitudes of the SDSS DR9 testing set as a function of the photometric redshift estimated by the SDSS DR9 collaboration. The lower panel shows the distribution of the residuals. The ANN reproduces very well the SDSS estimated photometric redshifts for almost all values from zphot;0 to zphot;0.8. The standard deviation of the residuals, over a quarter million of objects, reads sDzphot = 0.027. Because it is of the same order as the SDSS photometric redshift errors of the training catalog, the ANN regression is nearly optimal. Applying the trained ANN over the color magnitudes of our griz3 catalog gives the distribution of photometric redshifts of all galaxies located at less than 3′ from the SZ center associated with PSZ2 G156.24+22.32. In total, griz3 contains 114 objects, five of them end up having a negative photometric redshift. A closer examination of these outliers show that they have unusually large values for the estimated error on mag_auto (Δm0.2), and they are all located close to the bright star. They have therefore been removed from the analysis and are represented as white boxes in Figure 1. The photometric redshift distribution of all the 109 remaining objects has been plotted in the lower panel of Figure 3, while the upper panel shows their color–magnitude distribution. A clear reddening of the objects can be observed between the g and r magnitude, and to a lesser extend between r and i as well, supporting the claim that the sources are redshifted galaxies. The photometric redshift distribution indeed shows a peak around a redshift of 0.3 and no object are found with a redshift less than 0.08 (up to the five outliers). To estimate the mode of the redshift distribution, the “ash” package of the R-PROJECT software suite has been used to compute a polynomial density estimate and extract its maximum. One finds zphot = 0.29  0.08.

Boucher, de Visscher, & Ringeval

Figure 3. Upper panel: magnitude distribution in griz for the 109 galaxies observed within the 3′ radius region centered at the SZ coordinates. The lower panel is the corresponding photometric redshift distribution as inferred by our neural network trained over the SDSS DR9 catalog (five outliers having negative redshifts removed). (A color version of this figure is available in the online journal.)

the intrinsic ANN residuals plotted in Figure 2 due to the color–magnitude errors of our measurements. The most probable hydrostatic mass, MSZ(z), as well as its standard deviation, as a function of the cluster redshift z is given by the Planck collaboration in the MMF3 catalog(Planck Collaboration et al. 2016a). Combined with Equation (2), it yields MSZ = (4.4  1.3) ´ 1014 M☉.

(2 )

(3 )

The quoted errors are the intersect of the median absolute deviations associated with the photometric redshift and the standard deviation on MSZ(z). To test the robustness of Equation (2) with respect to the chosen condition under which galaxies belong to the cluster, we have also considered a region of 2′ radius around the SZ coordinates. This corresponds to the lower 1σ limit of θs in

To minimize sensitivity to potential systematics, the quoted error stands for the median absolute deviation around the mode, normalized such that it would match the usual standard deviation for a Gaussian distribution. Changing the smoothing kernel does not affect the estimation of the mode by more than 0.01. As expected, the median absolute deviation is larger than


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October

Boucher, de Visscher, & Ringeval

Equation (1), and to the objects constituting the griz2 color catalog. After the removal of the same five outliers as in the griz3 catalog, it remains 53 galaxies for which one finds zphot=0.29±0.11. As before, the first value is the mode and the error stands for the median absolute deviation around the mode.

3.3. Optical Richness A visual inspection of the red luminous objects contiguous to the ones contained within the 2′- and 3′-radius regions suggests that a galaxy population is also present further southeast (SE), i.e., toward the lower left corner of Figure 1. Let us notice that the NED database reports one X-ray source in this region, 1RXS J064506.9+592603(Voges et al. 1999). To quantify which objects are clustered, we have estimated the optical richness, λ, as discussed in Rozo et al. (2009). In practice, λ has been computed according to the method presented in Rykoff et al. (2012) and Rykoff et al. (2014). To use the same luminosity tracer as in these references, we have recalibrated our i′ filter over the SDSS DR9 reference field by using the i CmodelMag magnitude (instead of the modelMag). The color has been estimated as before, namely our g−r comes from g′ and r′ filter calibrated from modelMag magnitudes. This color should be an accurate tracer of the red-sequence galaxies up to redshift 0.35(Gladders & Yee 2000; Koester et al. 2007). All of the magnitudes have been corrected for Galactic extinction and reddening(Schlegel et al. 1998; Schlafly & Finkbeiner 2011), and the cosmological parameters have been chosen according to the Planck 2015 favored values(Planck Collaboration et al. 2016c). The presence of the parasiting star and the other uncertainties associated with our small instrumentation are expected to introduce various systematics in the determination of the photometric zero-points. As reported in Rykoff et al. (2012), the determination of λ is expected to be accurate for zero-point shifts up to ±0.05 magnitudes, and we might have systematic uncertainties as large as eight times this value for some objects (see Section 2.3). Adding an additional scatter of 0.1 magnitudes to luminosity and colors is indeed found to increase the value of λ by a factor of two. As a result, there is a possibility that the following evaluations of λ are underestimated, which is inherent to the experimental limitations of our setup. Notice, however, that the dependence of λ with respect to the coordinates of the clustering center and to the redshift, remains mostly insensitive to the uncertainties in the photometric zero-points. This is what we focus on in the following. To minimize the propagation of these potential systematics, we have chosen the radial filter to be ⎛ l ⎞b ⎟ , R c (l ) = R 0 ⎜ ⎝ 100 ⎠

Figure 4. Upper panel: map of the optical richness λ over the field of view of Figure 1 (the red-sequence galaxies are assumed to be at z = 0.29). The lower panel shows the dependence of λ(z) for the two locations exhibiting maximal richness in the upper panel. Only the one located East to the SZ center is well bounded in redshift, and can be associated with a galaxy cluster (red points). (A color version of this figure is available in the online journal.)

with R0=0.95h−1 Mpc, and β=0.1 instead of R0= 1h−1 Mpc, and β=0.2, as done in Rykoff et al. (2012). Our choice is still in the region of minimal intrinsic scatter for redsequence galaxies, but has the advantage of lowering the impact of underestimated λ values. The background galaxies and the red-sequence fit in g−r have been taken as in Rykoff et al. (2012), while the luminosity filter is the one of Rykoff et al. (2014). The upper panel of Figure 4 shows a map of λ over the field of view of Figure 1, obtained by postulating that the redsequence galaxies are at the redshift z=0.29. The two regions exhibiting maximal richness encompass the one visually expected east from the SZ center (left), around (α, δ)= (101.31, 59.48), and the one associated with the SE population at (α, δ)=(101.27, 59.36) (lower left). In the lower panel of

(4 )


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October

paper(Liu et al. 2015). However, the recent data release of the Pan-STARRS1 (PS1) survey(Chambers et al. 2016) covers the location of PSZ2 G156.24+22.32. As can be checked on the PS1 public archive5, the grizy stack images encompassing the SZ location show significant background subtraction and saturation artifacts coming from the central bright star. The depth of the griz channels also appears to be smaller than the one we obtained from the 14 hours of exposure (see Figure 3). This is not very surprising, as the small aperture of the 0.6 m telescope renders imaging less sensitive to saturation.6 Concerning the limiting magnitudes, as explained in Chambers et al. (2016), although the Pan-STARRS telescope collects 10 times more light than the T62, the exposure time on each patch is relatively small. For the plates encompassing PSZ2 G156.24+22.32, this is about 10 minutes for the grz channels, and 20 minutes for i. Nevertheless, the objects we have identified as galaxies are also present in the PS1 redder colors, while being hardly visible in g only. More interestingly, and although we have not attempted a quantitative comparison, the other southeast galaxy overdensity is present and more easily seen owing to the larger field of view accessible in the PS1.

Figure 5. Color–magnitude diagram for the all objects located within a Rc(λmax);1.1 Mpc radius around the cluster center (see the text and Figure 4). Red crosses are for galaxies identified as cluster members, while the black dots are for all the other objects. (A color version of this figure is available in the online journal.)

Figure 4, we have plotted the optical richness λ(z), as a function of the redshift, for these two locations. Only the first one ends up being well localized in redshift space, and this allows us to improve the cluster’s redshift estimate to +0.03 zl = 0.300.05.

Boucher, de Visscher, & Ringeval

4.2. Other Redshift Estimates We have also attempted to estimate the photometric redshift associated with the griz3 catalog by using a simple linear color fit over the g, r, i, and z magnitudes. The best fit has been determined over 100,000 galaxies of the SDSS DR9 catalog by iterative rejections of all objects remaining at more than 3σ from the best fit. This rejection is necessary as a brute force linear color fit over all SDSS DR9 galaxies ends up being reasonably tracking the photometric redshift only within the limited range 0.1zphot0.35. Applying the best fit to the griz3 catalog yields a photometric redshift distribution, again centered around a redshift 0.3, but with a much larger spreading than the one represented in Figure 3. We find zphot=0.35±0.16, the first value being again the mode and the second one stands for the absolute median deviation. Finally, we have tested an alternative redshift estimate by ignoring data taken with the z′ filter, as it shows some higher photometric errors than in the g′, r′, and i′ bands for faint objects (see Section 2.3). Following in all points the method presented in Section 3, but for the gri magnitudes only, one finds, within the 3′-radius region, zphot=0.27±0.06 for the mode and absolute median deviation. This is again compatible with the photometric estimate in Equation (2).

(5 )

The quoted errors are the intercepts for the FWHM. As can be seen in Figure 4, for the SE population of galaxies, λ(z) steadily increases up to z=0.42 to vanish at z=0.7 (not represented). For these redshifts, the choice of the color filter g−r is no longer accurately tracing red-sequence galaxies, but we can, however, conclude that most of these objects do not belong to the population east from the SZ center. Finally, in Figure 5, we have represented the color–magnitude diagram for all objects found within a Rc(λmax) radius around the location (101.31, 59.48), corresponding to the east population of Figure 4. The red crosses represent red-sequence galaxies, which have been identified as belonging to the cluster with a probability greater than 0.5. The scatter of half a magnitude in g−r for bright objects (i < 17) might be the result of the aforementioned systematics in the zero-points. As already discussed, they may bias λ to lower values by pulling galaxies off the red sequence. For fainter objects, the dispersion in g−r shows a clear degradation, which is expected due to the loss of signal, especially toward the limiting magnitude i20.

4. Discussion 4.1. Comparison with Pan-STARRS1


https://confluence.stsci.edu/display/PANSTARRS Two-color fits available http://cp3.irmp.ucl.ac.be/∼ringeval/upload/aape5/ twocolors

The SZ source under scrutiny belongs to the second Planck SZ catalog, and it does not appear in the Pan-STARRS follow-up



Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October

Boucher, de Visscher, & Ringeval

PSZ2 G138.32−39.82has a X-ray counterpart known as RX J0142.0+2131, which is located at a redshift zspec=0.280 (Barr et al. 2005). This cluster has been observed with the 0.5 m telescope and four non-standard interferential color filters available at that time, centered over the red, green, blue and near-infrared wavelengths, with about 2.5 hours of exposure per filter. Calibration over SDSS galaxies was used to determine a transformation matrix between these colors and the SDSS ugri modelMag, and a linear color-redshift fit gives zphot=0.30±0.1 (mode and absolute median deviation). PSZ2 G164.18−38.88,a close cluster belonging to the Abell (1958) catalog and known as Abel 401, at a redshift of z=0.0739. It has been observed in luminance with the T62 and for a total exposure time of 4.5 hours. The shape of galaxies is well resolved, and their internal structure, such as spiral arms, can be observed for many of them. However, no obvious galaxy overdensity can be inferred, as the 14′ field of view is too small for encompassing all the cluster’s galaxies. We have also observed the same field of view with the 0.5 m telescope in three non-standard colors (blue, green, and near-infrared) up to 1.5 hours per filter. A linear color-redshift fit gives zphot=0.09±0.08 for the mode and the absolute median deviation.

Figure 6. Low-signal instrumental test for the 0.6 m telescope during the 2015 mission. The image corresponds to a total of 1.5 hours of luminance exposure. Many galaxies of MACS J0647.7+7015, located at redshift z=0.584, are individually resolved, but the elliptic shape of only a few can be guessed.

4.3. Instrumental Tests on Confirmed Clusters

5. Conclusions

As mentioned in the introduction, instrumental tests have been undertaken during the 2015 mission to assess the feasibility of resolving individual galaxies in distant clusters with small telescopes. Successful identification of individual galaxies has been obtained for all the observed clusters, and we briefly report below our results for each of them.

This article presents the results of optical follow-up searches of Planck SZ sources with non-professional 0.5 m and 0.6 m telescopes. Our main result is the confirmation of the cluster candidate PSZ2 G156.24+22.32 by observing more than 100 galaxies within a 3′-radius region around the SZ coordinates and by estimating their photometric redshift at zphot=0.29±0.08. We have also estimated the optical richness, λ, around this location. Although systematics in the photometric zero-points are expected to bias λ to lower values, its dependence in redshift allows us to narrow down the cluster’s redshift to +0.03 zl = 0.300.05 . Imaging also reveals a contiguous population of galaxies extending a few arcminutes southeast from the SZ region. However, their optical richness suggests that they do not belong to the cluster, and are loosely concentrated around z;0.42. The robustness of these results have been tested against the choice of the region in which galaxies are counted, the colorredshift relationship, and by removing the z′ color data. Moreover, photometric redshift estimates of two already confirmed clusters of the PSZ2 catalog, made with the same telescopes, are compatible with their actual redshifts. The potential of individual sub-meter instruments in the optical confirmation of SZ sources and the measurement of their redshift is obviously limited with respect to what can be achieved with much larger telescopes. However, compared

PSZ2 G144.83+25.11,referred to as MACS J0647.7 +7015, is a massive and distant cluster located at z=0.584. It has been used as a low-signal test of the T62, and imaged in luminance for a total exposure time of 1.5 hours. As can be seen in Figure 6, individual galaxies are resolved, but the elliptical shape of only a few can be guessed. Overdensity in galaxy counts is obvious, but occupies a relatively small angular size, about 2′. PSZ2 G149.75+34.68,also known as RXC J0830.9 +6551, is reported in the Planck catalog with a X-ray redshift of zx=0.1818 (Mullis et al. 2003). The elliptic shape of the most massive galaxies is resolved after less than 1 hour of luminance exposure. PSZ2 G157.32−26.77,also MACS J0308.9+2645, at a redshift of z=0.356 (Ebeling et al. 2001). Galaxies are individually resolved after 1.6 hours of luminance exposure, but elliptic shape can only be inferred for the largest ones.


Publications of the Astronomical Society of the Pacific, 130:104001 (10pp), 2018 October with automated acquisition campaigns, their advantage is the possibility to acquire sensibly deeper data by accumulating longer exposure time over one target (see Section 4.1). The present results suggest that observing galaxy overdensity and estimating redshift are achievable in a redshift range 0.1 to 0.6, albeit with relatively large uncertainties for the redshift, up to 40% for the closest clusters. In addition, there are potentially hundreds of sub-meter class telescopes available around the world, mostly within the amateur astronomer community. With even a fraction of them, dedicated searches of SZ sources could be collectively envisaged, in a way similar to what the Galaxy Zoo project is providing(Lintott et al. 2008). A case-by-case tuning of the acquisition parameters (unitary exposure time in particular) can be operated to account for the presence of bright objects in the vicinity of the faint galaxies to be detected, which is more complicated to achieve in the case of an automatic survey. Such a collective campaign would certainly bring complementarity to automated surveys. At the very least, such a project could be used to provide a first information on the redshift before starting more accurate investigations with professional telescopes.

Boucher, de Visscher, & Ringeval

Bertin, E. 2006, in ASP Conf. Ser, Vol. 351, Astronomical Data Analysis Software and Systems XV, ed. C. Gabriel et al. (San Francisco, CA: ASP), 112 Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 Bertin, E., Mellier, Y., Radovich, M., et al. 2002, in ASP Conf. Ser., Vol. 281, Astronomical Data Analysis Software and Systems XI, ed. D. A. Bohlender, D. Durand, & T. H. Handley (San Francisco, CA: ASP), 228 Bleem, L. E., Stalder, B., de Haan, T., et al. 2015, ApJS, 216, 27 Bunn, E. F. 2006, PhRvD, 73, 123517 Challinor, A. D., Ford, M. T., & Lasenby, A. N. 2000, MNRAS, 312, 159 Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, arXiv:1612.05560 Ebeling, H., Edge, A. C., & Henry, J. P. 2001, ApJ, 553, 668 Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748 Gladders, M. D., & Yee, H. K. C. 2000, AJ, 120, 2148 Graff, P., Feroz, F., Hobson, M. P., & Lasenby, A. 2012, MNRAS, 421, 169 Graff, P., Feroz, F., Hobson, M. P., & Lasenby, A. 2014, MNRAS, 441, 1741 Hasselfield, M., Hilton, M., Marriage, T. A., et al. 2013, JCAP, 7, 008 Hilton, M., Hasselfield, M., Sifón, C., et al. 2017, arXiv:1709.05600 Hobson, M., Graff, P., Feroz, F., & Lasenby, A. 2014, in IAU Symp., Vol. 306, Statistical Challenges in 21st Century Cosmology, ed. A. Heavens, J-L Starck, & A. Krone-Martins (Cambridge: Cambridge Univ. Press), 279 Hurier, G., Aghanim, N., & Douspis, M. 2017, arXiv:1702.00075 King, E. S. 1989, S&T, 77, 387 Koester, B. P., McKay, T. A., Annis, J., et al. 2007, ApJ, 660, 221 Lintott, C. J., Schawinski, K., Slosar, A., et al. 2008, MNRAS, 389, 1179 Liu, J., Hennig, C., Desai, S., et al. 2015, MNRAS, 449, 3370 Monet, D. G., Levine, S. E., Casian, B., et al. 2003, AJ, 125, 984 Mullis, C. R., McNamara, B. R., Quintana, H., et al. 2003, ApJ, 594, 154 Piffaretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & Melin, J.-B. 2011, A&A, 534, A109 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014a, A&A, 571, A29 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014b, A&A, 571, A20 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2015a, A&A, 581, A14 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2015b, A&A, 582, A29 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016a, A&A, 594, A27 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016b, A&A, 594, A24 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016c, A&A, 594, A13 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016d, A&A, 586, A139 Rozo, E., Rykoff, E. S., Bartlett, J. G., & Melin, J.-B. 2015, MNRAS, 450, 592 Rozo, E., Rykoff, E. S., Koester, B. P., et al. 2009, ApJ, 703, 601 Rykoff, E. S., Koester, B. P., Rozo, E., et al. 2012, ApJ, 746, 178 Rykoff, E. S., Rozo, E., Busha, M. T., et al. 2014, ApJ, 785, 104 Rykoff, E. S., Rozo, E., Hollowood, D., et al. 2016, ApJS, 224, 1 Sazonov, S. Y., & Sunyaev, R. A. 1999, MNRAS, 310, 765 Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 Smith, J. A., Tucker, D. L., Kent, S., et al. 2002, AJ, 123, 2121 Sunyaev, R. A., & Zeldovich, I. B. 1980, MNRAS, 190, 413 Sunyaev, R. A., & Zeldovich, Y. B. 1970, Ap&SS, 7, 3 Sunyaev, R. A., & Zeldovich, Y. B. 1972, CoASP, 4, 173 Tody, D. 1986, Proc. SPIE, 627, 733 Tody, D. 1993, in ASP Conf. Ser., Vol.52, Astronomical Data Analysis Software and Systems II, ed. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes (San Francisco, CA: ASP), 173 van der Burg, R. F. J., Aussel, H., Pratt, G. W., et al. 2016, A&A, 587, A23 Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389 Vorobyev, V. S., Burenin, R. A., Bikmaev, I. F., et al. 2016, AstL, 42, 63 Wen, Z. L., Han, J. L., & Liu, F. S. 2012, ApJS, 199, 34 Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579 Zacharias, N., Finch, C., Subasavage, J., et al. 2015, AJ, 150, 101

We warmly thank the Astroqueyras Association for having allocated the time and resources needed to carry on this project. In particular, the acquisition of the set of Sloan filters was greatly appreciated. It is a pleasure to thank E.Bertin, H.J.McCracken, and the two anonymous referees for precious advice. We also acknowledge usage of the SDSS-III data, whose funding and participating institutions can be found athttp:// www.sdss3.org/. The optical richness code has been developed based on the one publicly provided by R.H.Wechsler athttp://risa.stanford.edu/redmapper. We have made use of the NED database, operated by the Jet Propulsion Laboratory, the California Institute of Technology, under NASA contract, and the SIMBAD database, operated at the Centre de Données Astronomique de Strasbourg, France (Wenger et al. 2000). Preprocessing of raw images has been made with IRAF, distributed by the National Optical Astronomy Observatories. Catalog manipulations have been made usingSAOImage DS9.

References Abell, G. O. 1958, ApJS, 3, 211 Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS, 203, 21 Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, L37 Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92 Barr, J., Davies, R., Jørgensen, I., Bergmann, M., & Crampton, D. 2005, AJ, 130, 445


Suggest Documents