The eROSITA telescope array aboard the Spektrum Roentgen Gamma SRG satellite began surveying the sky in December 2019 with the aim of producing allsky Xray source lists and sky maps of an unprecedented depth. Here we present catalogues of both pointlike and extended sources using the data acquired in the first six months of survey operations eRASS1 completed June 2020 over the half sky whose proprietary data rights lie with the German eROSITA Consortium. We describe the observation process the data analysis pipelines and the characteristics of the Xray sources. With nearly 930 000 entries detected in the most sensitive 0.22.3 keV energy range the eRASS1 main catalogue presented here increases the number of known Xray sources in the published literature by more than 60 and provides a comprehensive inventory of all classes of Xray celestial objects covering a wide range of physical processes. A smaller catalogue of 5466 sources detected in the less sensitive but harder 2.35 keV band is the result of the first true imaging survey of the entire sky above 2 keV. We present methods to identify and flag potential spurious sources in the catalogues which we applied for this work and we tested and validated the astrometric accuracy via crosscomparison with other Xray and multiwavelength catalogues. We show that the number counts of Xray sources in eRASSl are consistent with those derived over narrower fields by past Xray surveys of a similar depth and we explore the number counts variation as a function of the location in the sky. Adopting a uniform allsky flux limit at 50 completeness of FSUB052 keVSUB gt 5  10SUP14SUP erg sSUP1SUP cmSUP2SUP we estimate that the eROSITA allsky survey resolves into individual sources about 20 of the cosmic Xray background in the 12 keV range. The catalogues presented here form part of the first data release DR1 of the SRGeROSITA allsky survey. Beyond the Xray catalogues DR1 contains all detected and calibrated event files source products light curves and spectra and allsky maps. Illustrative examples of these are provided. P The catalogue is available at the CDS via anonymous ftp to A hrefhttpscdsarc.cds.unistra.frcdsarc.cds.unistra.frA ftp130.79.128.5 or via A hrefhttpscdsarc.cds.unistra.frvizbincatJAA682A34httpscdsarc.cds.unistra.frvizbincatJAA682A34A
['catalogs', 'surveys', 'X-rays: general', 'Astrophysics - High Energy Astrophysical Phenomena']
The SRGeROSITA allsky survey. First Xray catalogues and data release of the western Galactic hemisphere
{'First X-ray catalogues and data release of the western Galactic hemisphere ⋆': 'A. Merloni 1 ⋆⋆ , G. Lamer 2 , T. Liu 1 , 25 , 26 , M. E. Ramos-Ceja 1 , H. Brunner 1 , E. Bulbul 1 , K. Dennerl 1 , V. Doroshenko 3 , M. J. Freyberg 1 , S. Friedrich 1 , E. Gatuzz 1 , A. Georgakakis 4 , F. Haberl 1 , Z. Igo 1 , 5 , I. Kreykenbohm 6 , A. Liu 1 , C. Maitra 1 , A. Malyali 1 , M. G. F. Mayer 1 , K. Nandra 1 , P. Predehl 1 , J. Robrade 7 , M. Salvato 1 , 5 , J. S. Sanders 1 , I. Stewart 1 , D. Tubín-Arenas 2 , P. Weber 6 , J. Wilms 6 , R. Arcodia 1 , 8 , E. Artis 1 , J. Aschersleben 9 , A. Avakyan 3 , C. Aydar 1 , 5 , Y. E. Bahar 1 , F. Balzer 1 , W. Becker 1 , 10 , K. Berger 6 , T. Boller 1 , W. Bornemann 1 , M. Brüggen 7 , M. Brusa 11 , 12 , J. Buchner 1 , V. Burwitz 1 , F. Camilloni 1 , N. Clerc 13 , J. Comparat 1 , D. Coutinho 1 , S. Czesla 14 , 7 , S. M. Dannhauer 9 , L. Dauner 6 , T. Dauser 6 , J. Dietl 9 , K. Dolag 15 , T. Dwelly 1 , K. Egg 6 , E. Ehl 1 , S. Freund 1 , 7 , P. Friedrich 1 , R. Gaida 1 , C. Garrel 1 , V. Ghirardini 1 , A. Gokus 6 , 16 , G. Grünwald 1 , S. Grandis 17 , I. Grotova 1 , D. Gruen 15 , A. Gueguen 1 , S. Hämmerich 6 , N. Hamaus 15 , G. Hasinger 18 , 19 , K. Haubner 6 , D. Homan 2 , J. Ider Chitham 1 , W. M. Joseph 3 , A. Joyce 6 , O. König 6 , D. M. Kaltenbrunner 1 , A. Khokhriakova 1 , W. Kink 1 , C. Kirsch 6 , M. Kluge 1 , J. Knies 6 , S. Krippendorf 15 , M. Krumpe 2 , J. Kurpas 2 , P. Li 2 , Z. Liu 1 , N. Locatelli 1 , M. Lorenz 6 , S. Müller 1 , E. Magaudda 3 , C. Mannes 9 , H. McCall 9 , 20 , N. Meidinger 1 , M. Michailidis 3 , K. Migkas 9 , 21 , D. Muñoz-Giraldo 3 , B. Musiimenta 11 , 12 , N. T. Nguyen-Dang 3 , Q. Ni 1 , A. Olechowska 1 , N. Ota 9 , F. Pacaud 9 , T. Pasini 10 , 22 , E. Perinati 3 , A. M. Pires 2 , C. Pommranz 3 , G. Ponti 23 , 1 , K. Poppenhaeger 2 , G. Pühlhofer 3 , A. Rau 1 , M. Reh 6 , T. H. Reiprich 9 , W. Roster 1 , S. Saeedi 6 , A. Santangelo 3 , M. Sasaki 6 , J. Schmitt 7 , P. C. Schneider 7 , T. Schrabback 17 , N. Schuster 15 , A. Schwope 2 , R. Seppi 1 , M. M. Serim 3 , S. Shreeram 1 , E. Sokolova-Lapa 6 , H. Starck 1 , B. Stelzer 3 , J. Stierhof 6 , V. Suleimanov 3 , C. Tenzer 3 , I. Traulsen 2 , J. Trümper 1 , K. Tsuge 6 , T. Urrutia 2 , A. Veronica 9 , S. G. H. Waddell 1 , R. Willer 1 , J. Wolf 1 , 5 , 24 , M. C. H. Yeung 1 , A. Zainab 6 , F. Zangrandi 6 , X. Zhang 1 , Y. Zhang 1 , X. Zheng 1  \n(A ffi liations can be found after the references)  \nReceived 12 June 2023; accepted 26 November 2023', 'ABSTRACT': 'The eROSITA telescope array aboard the Spektrum Roentgen Gamma (SRG) satellite began surveying the sky in December 2019, with the aim of producing all-sky X-ray source lists and sky maps of an unprecedented depth. Here we present catalogues of both point-like and extended sources using the data acquired in the first six months of survey operations (eRASS1; completed June 2020) over the half sky whose proprietary data rights lie with the German eROSITA Consortium. We describe the observation process, the data analysis pipelines, and the characteristics of the X-ray sources. With nearly 930 000 entries detected in the most sensitive 0.2-2.3 keV energy range, the eRASS1 main catalogue presented here increases the number of known X-ray sources in the published literature by more than 60%, and provides a comprehensive inventory of all classes of X-ray celestial objects, covering a wide range of physical processes. A smaller catalogue of 5466 sources detected in the less sensitive but harder 2.3-5 keV band is the result of the first true imaging survey of the entire sky above 2 keV. We present methods to identify and flag potential spurious sources in the catalogues, which we applied for this work, and we tested and validated the astrometric accuracy via cross-comparison with other X-ray and multi-wavelength catalogues. We show that the number counts of X-ray sources in eRASS1 are consistent with those derived over narrower fields by past X-ray surveys of a similar depth, and we explore the number counts variation as a function of the location in the sky. Adopting a uniform all-sky flux limit (at 50% completeness) of F 0 . 5 -2 keV > 5 × 10 -14 erg s -1 cm -2 , we estimate that the eROSITA all-sky survey resolves into individual sources about 20% of the cosmic X-ray background in the 1-2 keV range. The catalogues presented here form part of the first data release (DR1) of the SRG / eROSITA all-sky survey. Beyond the X-ray catalogues, DR1 contains all detected and calibrated event files, source products (light curves and spectra), and all-sky maps. Illustrative examples of these are provided.  \nKey words. catalogues - surveys - X-ray: general', '1. Introduction': "Wide-area multi-wavelength (and multi-messenger) sky surveys play a key role in the development of astrophysics. The design of these surveys is driven both by the desire to open up new discov-  \nery spaces, and by the realisation that charting the structure of the Universe in detail may help solve long-standing open questions of fundamental physics (see e.g. Peebles 1980; Weinberg et al. 2013).  \nSky surveys also occupy a central position in the relatively short history of X-ray astronomy (see e.g. Elvis 2020, for a succinct recap). These surveys were initiated by the SASA Uhuru satellite (1970-1973, Giacconi et al. 1971), a mission designed to conduct a survey of the X-ray sky in the 2-  \n20 keV energy range, resulting in the identification of 339 Xray sources (Forman et al. 1978). Likewise, an X-ray survey was carried out by the cosmic X-ray experiment aboard the HEAO-1 observatory (1977-1979, Rothschild et al. 1979) covering the 0.25-25 keV energy range, enabling the discovery of 842 X-ray sources (Wood et al. 1984). The Einstein (HEAO-2) X-ray Observatory (1978-1981, Giacconi et al. 1979) performed a medium-sensitive survey in the 0.3-3.5 keV energy range, covering roughly one-third of the sky and providing the detection of about 4000 X-ray sources (Harris 1990). However, the first comprehensive all-sky X-ray survey was performed by the ROSAT X-ray observatory in the 0.1-2.4 keV energy range (1990-1999, Truemper 1982).  \nThe ROSAT All-Sky Survey (RASS) and its associated catalogues of X-ray sources, The Bright Source Catalogue (BSC) (Voges et al. 1999), Faint Source Catalogue (FSC) (Voges et al. 2000), and the more comprehensive second catalogue release 2RXS (Boller et al. 2016), marked a new milestone in the quantity and quality of X-ray detections. With more than 10 5 unique sources detected in only six months of survey observations, RASS outnumbered all previous all-sky catalogues by more than two orders of magnitude.  \nSince the turn of the century, XMM-Newton and Chandra , with their large collecting area and high spatial resolution, respectively, have increased the number of known X-ray sources significantly. Chandra and XMM-Newton serendipitous source catalogues, however, only cover a small fraction of the sky, due to their relatively small fields of view and their pointed observation strategy.  \nThese and other X-ray surveys provide a unique view of many celestial phenomena. X-ray emission is a universal signature of accretion of matter onto compact objects. In binary systems, these mark the graveyards of stellar evolution (see e.g. Shapiro & Teukolsky 1983; Warner 1995; Becker & Truemper 1997; Fender et al. 2004; Remillard & McClintock 2006; Haberl 2007; Özel & Freire 2016), where black holes, neutron stars, and white dwarfs accrete matter from a companion. The brightest Xray sources in the sky, including Sco X-1, the first extrasolar one discovered (Giacconi et al. 1962), fall into this category. At cosmological distances, X-rays signpost the growth of the supermassive black holes (SMBHs) that sit at the centres of galaxies and which may strongly influence their formation and subsequent evolution (Brandt & Hasinger 2005; Hopkins et al. 2008; Hickox et al. 2009; Fabian 2012; Alexander & Hickox 2012; Kormendy & Ho 2013; Brandt & Alexander 2015).  \nX-ray emission is seen from gas heated by the shocks generated by rapidly expanding supernova remnants, which are in addition likely responsible for a large fraction of the accelerated particles that di ff use through interstellar space (see e.g. Vink 2012; Blasi 2013, and references therein). Strong X-ray emission is also seen from the coronae of stars (Schmitt 1997; Pizzolato et al. 2003; Wright et al. 2011), which plays an important role in the determination of the physical characteristics of the orbiting planets' atmospheres, with implications for the potential habitability of such planets (Lammer et al. 2003; Sanz-Forcada et al. 2011; Poppenhaeger et al. 2021; Foster et al. 2022).  \nAfurther source of X-ray emission of great wider importance is from hot gas associated with large-scale structures. Most of the baryons in the Universe are indeed predicted to be locked into a warm-hot (X-ray emitting) tenuous phase (Cen & Ostriker 1999; Davé et al. 2001; Nicastro et al. 2018; Tanimura et al. 2020). Direct or indirect detections of these baryons often require sensitive mapping of large volumes in X-rays. Moreover, in the hierarchical distribution of matter characteristic of  \nour observed Universe, the densest knots of the large-scale structure are signposted by the hottest and most massive concentration of di ff use baryons. The clusters of galaxies that identify them are thus extremely sensitive tracers of the underlying geometry and growth history of the cosmos, and therefore prime cosmological tools (Bahcall 1977; Cavaliere & Fusco-Femiano 1978; Sarazin 1986; Rosati et al. 2002; Voit 2005; Arnaud 2005; Norman 2005; Borgani 2008; Vikhlinin et al. 2009; Borgani & Kravtsov 2011; Allen et al. 2011; Reiprich et al. 2013). It is the potential for sensitive cosmological measurements with galaxy clusters that provided the main motivation for the development of eROSITA in the early 2000s, when it became apparent that a significantly larger number of clusters of galaxies compared to those provided by narrow-field instruments would be required to e ff ectively constrain the fundamental parameters of cosmological models (Haiman et al. 2005; Merloni et al. 2012; Pillepich et al. 2012). In particular, these authors argued that sample sizes of order 10 5 clusters were required to provide constraints competitive with other prime cosmological measurement tools.  \neROSITA (extended ROentgen Survey with an Imaging Telescope Array; Predehl et al. 2021), on board the Spektrum Roentgen Gamma (SRG) orbital observatory (Sunyaev et al. 2021), was developed in the period 2007-2019. It is a sensitive, widefield focusing X-ray telescope array, optimised to deliver large e ff ective area and field-of view (hence also grasp and survey speed) in the soft X-ray band. The angular resolution is su ffi -cient to distinguish between the two largest extragalactic source populations, that is, AGN and clusters of galaxies, via measurement of their X-ray extent.  \nThe observing strategy was designed to achieve the needed sensitivity in the soft X-ray band (below 2 keV) to detect at least the requisite 10 5 clusters of galaxies by scanning the entire sky eight times over a period of four years (the eROSITA All-Sky Surveys, hereafter: eRASS).  \nSRG was launched on July 13, 2019 from Baikonur, Kazakhstan, using a Proton-M rocket and a BLOK DM-03 upper stage. On its three months cruise to the second Lagrangian point (L2) of the Earth-Sun system the spacecraft and instruments underwent commissioning and checkout. Since mid-October 2019, SRG was placed in a six-month-periodic halo orbit around L2, with a semi-major axis of about 750 000 km within the ecliptic plane and semi-minor axis of about 400 000 km perpendicular to it (Sunyaev et al. 2021); periodic orbit correction manoeuvres over the intervening years have slightly reduced the size of the orbit in order to satisfy ground segment visibility constraints.  \nFollowing First Light (Maitra et al. 2022; Haberl et al. 2022), a two-months long Calibration and Performance Verification (CalPV) programme was executed between October and December 2019. The large body of publications based on CalPV observations (see Campana et al. 2021, and the associated articles of the A&A special issue) have demonstrated the capabilities of eROSITA, and confirmed the main design characteristics. In particular, the eROSITA Final Equatorial Depth Survey (eFEDS; Brunner et al. 2022), designed to provide uniform exposure over a su ffi ciently large field (140 deg 2 ) about 50% deeper than what is expected for eRASS at the end of the 4 year all-sky survey programme, has shown that large samples of X-ray sources of di ff erent classes can be detected, identified, and characterised making use of the synergy with existing multi-wavelength surveys (see, e.g. Ghirardini et al. 2021; Liu et al. 2022a,b; Salvato et al. 2022; Klein et al. 2022; Schneider et al. 2022; Bulbul et al. 2022; Pasini et al. 2022; Ramos-Ceja et al. 2022, Nandra et al., submitted).  \nThe eROSITA data are shared equally between German and Russian scientists, following an inter-agency agreement signed in 2009. Two hemispheres of the sky have been defined, over which each team has unique scientific data exploitation rights. These data rights are split by Galactic longitude ( l ) and latitude ( b ), with a division marked by the great circle passing through the Galactic poles ( l , b ) = (0 · , + 90 · ); (0 · , -90 · ) and the Galactic Center Sgr A* ( l , b ) = (359 . 9442 · , -0 . 04616 · ): data with -0 . 05576 · < l < 179 . 9442 · (eastern Galactic hemisphere) belong to the Russian consortium, while data with 359 . 9442 · > l > 179 . 9442 · (western Galactic hemisphere) belong to the German eROSITA consortium (eROSITA-DE). Here we only describe and release the data collected in the western Galactic hemisphere.  \nAfter a brief introduction to the salient technical aspects of eROSITA and its calibration (Sect. 2), in Sect. 3 we describe in detail the eROSITA observing procedures in all-sky survey mode. Section 4 is then devoted to a summary of the main data processing stages. Most of the details of the software system used to process eROSITA data have already been presented in Brunner et al. (2022), and we refer the interested reader to that work for more information. The catalogues generated by our standard processing pipeline analysing the first eROSITA allsky survey (eRASS1) are presented in Sect. 5. There we present the preliminary astrometric correction applied to the detected sources as well as our general attempt to isolate and flag potential spurious sources and other artefacts. Section 6 describes further consistency checks on the X-ray catalogue performed by comparing the eRASS1 sources with those of the XMMNewton serendipitous catalogue (photometry) and other multiwavelength Quasi-Stellar Object (QSO) catalogues (astrometry). Section 7 then gives an overview of all the available products, including all-sky maps, and source-specific light curves and X-ray spectra. To conclude, we summarise our work in Sect. 8.", '2. eROSITA technical facts': 'In this section, we provide a compact summary of the main technical characteristics of eROSITA and its calibration; more details, including a summary of the on-ground calibration, can be found in Predehl et al. (2021). Specific technical descriptions of the instrument subsystems can be found in Meidinger et al. (2021) (camera system) and Friedrich et al. (2012) (mirror system), while a thorough description of the on-ground calibration campaign and its results can be found in Dennerl et al. (2020). Indeed, most of the data analysis and pipeline settings adopted for the reduction of the eRASS1 data presented here rely on this extensive on-ground calibration of the instrument. As for all other X-ray space observatories, in-flight calibration is a longterm endeavour; here we present some preliminary results that demonstrate the fidelity of the calibration and the reliability of the data released. An in-depth analysis of the in-flight calibration will be presented elsewhere, and on the website of the first data release (DR1).', '2.1. Instrument characteristics': "eROSITA consists of seven identical and co-aligned X-ray telescopes arranged in a common optical bench. A system of carbon fibre honeycomb panels connects the seven mirror assemblies on the front side with the associated seven camera assemblies on the focal plane side.  \nEach of the mirrors comprises 54 mirror shells in a Wolter-I geometry, with an outer diameter of 360 mm and a common focal  \nlength of 1 600 mm (Friedrich et al. 2008; Arcangeli et al. 2017). The average on-axis resolution of the seven mirrors, as measured during the on-ground calibration, is 16 . 1 '' Half-Energy Width (HEW) at 1.5 keV. The unavoidable o ff -axis blurring typical of Wolter-I optics is compensated by a 0.4 mm shift of the cameras towards the mirrors. This puts each telescope slightly out of focus, leading to a slight degradation of the on-axis performance (18 '' HEW), but improved angular resolution averaged over the field of view. Indeed, in the scanning observational mode adopted for the all-sky survey (see section 3.1), it is the field-of-view average HEW that matters. This is discussed in Section 2.2.1 and Appendix A.  \nEach Mirror Assembly has a CCD camera in its focus (Meidinger et al. 2014). The eROSITA CCDs are advanced versions of the EPIC-pn CCDs on XMM-Newton (Strüder et al. 2001). They have 384 × 384 pixels in an image area of 28 . 8 mm × 28 . 8 mm, yielding a square field of view of 1 · . 03 × 1 · . 03. Each pixel corresponds to a sky area of 9 '' . 6 × 9 '' . 6. The nominal integration time for all eROSITA CCDs is 50 msec. The additional presence of a frame-store area in the CCD reduces substantially the amount of so-called out-of-time events, which are recorded during the CCD read-out, a significant improvement with respect to the EPIC-pn camera on XMM-Newton . For optimal performance during operations, the CCDs are cooled down to about -85 · by means of passive elements (Fürmetz et al. 2008). During survey operation (i.e. in scanning mode of observations), the angle between the scanning direction projected onto the sky and the CCD read-out direction in the focal plane is not the same for the seven cameras 1 . This contributes to averaging out any possible non-circular symmetry of the PSF as well as camera-induced defects. For calibration purposes, each camera has its own filter wheel with a radioactive 55 Fe source and an aluminium and titanium target providing three spectral lines at 5.9 keV (Mn K α ), 4.5 keV (Ti K α ) and 1.5 keV (Al K α ).  \nThe electronics for onboard-processing of the camera data is provided by seven sets of Camera Electronics (CE), each one mounted and interfacing to the Cameras. Each of the CEs provide the proper voltage control and readout timing of the associated camera, and performs the on-board data processing within the time constraints of the camera integration time.", '2.2. eROSITA calibration': "The quantities derived from the eROSITA calibration are stored in a calibration database (CalDB), which is accessed by various processing tasks. Information about the content of the CalDB and how important entries were derived can be found in Brunner et al. (2022), Appendix B. A 'live' online repository can be found on the DR1 website 2 . Here we present a brief assessment of the current status of the in-flight calibration program.", '2.2.1. PSF calibration': "Given the scientific objectives of eROSITA, i.e. imaging the whole sky at soft-to-hard X-ray energies with good sensitivity to low surface brightness di ff use and extended emission, clusters of galaxies in particular, an accurate calibration of the tele-  \nscopes' Point-Spread Function (PSF) is critical. As we describe in greater detail below (§ 4), the catalogues presented here have been constructed using a single photon mode, in which the PSF of each telescope module at the location of each detected event on the CCD is accounted for, using a shapelet-based model calibrated on the extensive dataset accumulated on ground before the launch (Dennerl et al. 2020). A description of the shapeletbased PSF modelling, and its usage in scanning mode can be found in Appendix B.1 of Brunner et al. (2022).  \nWhile work is ongoing to accurately characterize eROSITA's PSF based on the survey data themselves, here we report on a preliminary analysis that confirms the reliability of the onground calibration adopted for the DR1 datasets. Appendix A shows a direct comparison of the average PSF shape (obtained by combining all seven Telescope Modules) from stacking point sources detected in the all-sky survey with the PSF model from the PANTER on-ground calibration and its shapelet representation. The di ff erences between the PSF models are within the ∼ 20% level in the inner 1 ' , although the shapelet PSF drops below the measured PSFs beyond this radius, where it is not used in the source detection process, except for normalisation. The measured HEWs from the source stacking method applied to survey data are 30.0 '' in the 0.2-2.3 keV band and 34.4 '' in the 2.3-5.0 keV band, very close to the pre-flights estimates of 28.3 '' and 36.2 '' , respectively, for the shapelet representation, and 32.0 '' and 38.0 '' for the PANTER ground-based values. The PSF does not appear to be varying across the sky, at least within the limited statistics of our preliminary analysis (see Tab. A.1). In Appendix A we also show an estimate of the PSF azimutal symmetry, and compare positional o ff sets of eRASS1 point sources against Gaia QSOs along equatorial and ecliptic coordinates, demonstrating the goodness of our circular symmetric approximation for the positional uncertainty of the detected sources.", '2.2.2. Energy calibration': 'After launch, the energy calibration obtained on ground was checked by using extensive measurements with the internal 55 Fe calibration sources. These measurements showed that the inflight energy resolution of the detectors was essentially the same as on ground (Tables 1 and 2 in Dennerl et al. 2020). The 55 Fe calibration measurements were also used to derive updated values of the Charge Transfer Ine ffi ciency (CTI) and Gain, in order to minimise their impact on the absolute energy scale. Additional tests of the energy calibration were made with dedicated observations of astrophysical calibration targets, especially of the isolated neutron star RX J1856-3754 and of the supernova remnant 1E0102.2-7219. These demonstrated that the energy calibration is su ffi ciently accurate for the analysis of survey data with their limited photon statistics 3 . Appendix B provides more details on these energy calibration experiments. The currently available energy calibration has already been used successfully for spectroscopic studies (e.g. Camilloni et al. 2023; Mayer et al. 2023; Ponti et al. 2023b; Yeung et al. 2023).', '2.2.3. Flux calibration: Effective area and vignetting': 'The challenge of accurately flux-calibrating space-based X-ray telescopes is as old as X-ray astronomy itself (see e.g. Marshall  \net al. 2021; Madsen et al. 2021, and references therein for recent discussions). While the Quantum E ffi ciency (QE) of the CCD detectors can be accurately measured on ground, but is subject to degradation in the harsh space environment, the main di ffi culty rests with the paucity of suitable stable standard candles and the impossibility to accurately reproduce in the laboratory the observing conditions in space needed to calibrate a telescope e ff ective area and vignetting function (i.e. variation of e ff ective area across the focal plane). Clusters of galaxies, possibly the brightest intrinsically non-varying X-ray sources in the sky, are usually adopted as cross-calibrators among di ff erent Xray observatories, but even after decades, uncertainties remain (Nevalainen & Molendi 2023). For eROSITA, we can take advantage of the all-sky survey nature of the observations to build large samples that can be used to validate the flux calibration in a statistical sense. For point sources, we report here (Section 6.1) a comparison with XMM-Newton in the 0.5-2 keV band, resulting in a possible residual flux systematic uncertainty of just about 6%, similarly to what was found by Maitra et al. (2022) in the early phases of the mission. We have also tested the e ff ect of possible PSF calibration uncertainties on the reconstructed flux, by running an alternative source detection pipeline which makes use of the PSF image preliminary stacks, and found only a 3% source counts systematic o ff set compared to the catalogue described here.  \nThe situation for clusters of galaxies in the survey is still under examination (but see e.g. Whelan et al. 2022; Sanders et al. 2022, for pointed Calibration observations), with Bulbul et al. (2024) reporting a systematic flux deficit with respect to Chandra of about 15%, while Liu et al. (2023a) and Migkas et al. (submitted) found a temperature o ff set with respect to both XMM-Newton and Chandra (with eROSITA measuring lower tempertures than both) that increases with the cluster temperature itself. Work is ongoing to determine to what extent this is induced by calibration uncertainties or by cluster-related astrophysical e ff ects (such as multi-temperature ICM distributions).', '3.1. Scanning strategy': "In order to complete an all-sky survey, the SRG observatory rotates continuously with a period of four hours around an axis pointed near the direction of the Sun. This gives a scan rate of 0.025 deg s -1 . The rotation axis slowly shifts by approximately one degree per day following the motion of the Earth (and of the L2 point) around the Sun. Following this scanning pattern, eROSITA observes the entire sky in about six months, and observes each point in the sky typically six times ('visits') for up to 40 seconds over a day at the ecliptic equator; towards the ecliptic poles, the sources remain observable for more than 24 hours, and are therefore scanned more than six times. Indeed, all great circles (individual scans) intersect in the north and south ecliptic poles in the sky, creating regions of deep exposure and longer visibility periods. In addition, a slight inhomogeneity in the sky coverage is introduced by the elongated halo orbit around L2 and the nonuniform angular movement of the spacecraft rotation axis, which compensate the separations between Sun and Earth, as seen from the spacecraft, in order to maintain the Earth in the cone of the downlink antenna (further details in Predehl et al. 2021; Sunyaev et al. 2021).  \nTable 1. Timeline of the main eROSITA operations and major events during eRASS1.", '3.2. All-sky survey operations': "After the commissioning of the instrument and the CalPV phase, the first all-sky survey started on December 12, 2019, and was completed on June 11, 2020, for a total survey duration of 184 days. During this period, ground contacts with SRG and eROSITA took place every day, without interruption. A timeline of the most significant operations milestones during eRASS1 is presented in Table 1.  \nTable 2. Timeline of extended diagnostic and engineering exposures during eRASS1.  \nNotes. 'Raw Frames' are read-out cycles where each pixel is transmitted (by-passing the event processor), 10 such frames were commanded for each TM consecutively (due to the large amount of telemetry there is some time gap in between). After the seven TMs a second round of 'raw frames' was commanded, giving in total 10 × 7 × 2 = 140 such frames per epoch.  \nAs described in Coutinho et al. (2022) and Predehl et al. (2021), the Mission Control Center, located in Moscow at NPO Lavochkin (NPOL), has mainly used two deep-space antennas for the science downlink (Ussuriysk and Bear Lakes) 4 . In total, approximately 75 GB of telemetry data from eROSITA were dumped during eRASS1, with an average of ∼ 407 MB per day.  \nThe overall observing e ffi ciency of eROSITA during eRASS1 was ∼ 96.5%. This e ffi ciency has been calculated by taking the average of the Good Time Intervals (GTI) with respect to the total observing time for each camera, as they are operating independently. The main observation disruptions responsible for the loss of e ffi ciency come from Camera Electronics (CE) and Interface and Thermal Controller (ITC) anomalies: the CEs and ITC are susceptible to Single Event Upsets (SEU), which can interrupt the functioning of the instrument. Moreover, Telescope Modules (TM) 5 and 7 su ff er from a time-varying light leak, which can lead to loss of data and consequently to gaps in GTI  \n(see Sect. 4.6 and Predehl et al. 2021; Coutinho et al. 2022, for further details). Figure 1 shows the cumulative observing e ffi -ciency as a function of time for each TM during eRASS1.  \nRegular 'Filter Wheel Closed' (FWC) observations were carried out during the survey to monitor the instrumental background, starting in February 2020. The filter wheel of one camera per day was set to the closed position for 30 minutes. To reduce the impact on the exposure of the survey observations, the same camera had an FWC observation every seven days. In total, around 18 short FWC observations for each camera were performed in eRASS1 (exposure fraction 0.3%), providing a precious data set to model and monitor the background (see Yeung et al. 2023, for more details about the instrumental background model). In addition, more extended FWC observations were performed during periods of orbit corrections (see Table 2), along with other diagnostic and engineering exposures. The viewing direction of the telescope was reset after the orbit corrections (and monitoring pointed observations) to ensure survey coverage without gaps.  \nWith the exception of the short-lived SEU-induced malfunctions of the camera electronics and of the ITC (on Feb. 10, 2020), all eROSITA sub-systems were fully functional during eRASS1, and they have not su ff ered any permanent damage, apart from the expected degradation caused by external environmental e ff ects along SRG's L2 orbit. In Coutinho et al. (2022) the interested reader can find more details about the technical performance of eROSITA during the first two years of operations (from eRASS1 to eRASS4).", '3.3. Pre-processing and archiving': "The eROSITA data received at a ground station are forwarded in real time via a socket connection to the SRG operations centre and to IKI (the Space Research Institute of the Russian Academy of Sciences), where they are stored in binary files with a size of approximately 7.5 MB each. These files are transferred to MPE via a data exchange server. The files are then picked up by the pre-processing pipeline, which converts the telemetry into FITS files and forward them to two separate pipelines: the preprocarchiver and the Near Real Time Analysis (NRTA).  \nThe archive is organised into 'erodays', i.e. fixed intervals of 4 hours (also corresponding to the duration of one revolution of SRG in all-sky survey mode). Once an eroday is completed, the data are moved to the regular archive and the post-ingest pipeline is triggered.  \nThe NRTA pipeline has been developed in order to (i) monitor interactively the instrument parameters and check the health of all sub-systems on the shortest possible timescales, and (ii) alert the team and the community of possible interesting timedomain phenomena observed by eROSITA. The NRTA and its functionality is described in Appendix C. At this point, the  \nFig. 1. Cumulative observing e ffi ciency (as a percentage of the elapsed survey time) as a function of time for each eROSITA camera (labeled according to the associated Telescope Module, TM) during eRASS1. With the exception of TM5, badly a ff ected by light leak, all cameras achieved an e ffi ciency of more than 95% by the end of eRASS1, with TM1 reaching close to 99%. The largest loss of net observing time was caused by an ITC malfunction on February 10, 2020.  \n<!-- image -->  \narchived data are ready to be processed by the standard analysis pipeline.", '4. Standard data processing': 'The eRASS1 data were processed with the eROSITA standard data processing pipeline, operated at MPE. The pipeline consists of modules for event processing (TEL processing chain), event file and image creation (EXP processing chain), exposure and background map creation and source detection (DET processing chain), as well as for the creation of source-specific products such as spectra and light curves (SOU processing chain). Each processing chain executes a number of software tasks, which are part of the eROSITA Science Analysis Software System (eSASS).  \nIn comparison with the data processing version (001) from the Early Data Release (EDR) 5 , the event calibration in the processing version (010) of eRASS1 has a stronger telescope module (TM) specific noise suppression of doubles, triples and quadruples 6 , a better computation of the subpixel position, a corrected flagging of pixels next to bad pixels, and improved accuracy of projection. The list of eSASS task and calibration database versions used for this work is provided in Appendix E. A list of standard data products included in the eROSITA DR1 data release is provided in Sect. 7. A detailed description of the main eSASS tasks, associated eROSITA calibration database, and calibrated data products is available in Brunner et al. (2022) (Appendix A-C). In the rest of this section we briefly summarize the organisation of the pipeline and the function of each pipeline chain.  \nFig. 2. Spherical projection of the eRASS1 1B exposure map centred on the western Galactic hemisphere (i.e. the eROSITA-DE area) with overimposed in green the tiling of the sky into equal-area overlapping 3 . 6 · × 3 . 6 · tiles into which data are organised and analysed by the pipeline.  \n<!-- image -->', '4.1. TEL chain': 'The TEL chains perform the functions of event file preparation (tasks evprep , ftfindhotpix ), pattern recombination (task pattern ), energy calibration (task energy ), and attitude calculation (tasks attprep , telatt , and evatt ). The TEL chains are executed separately for each telescope module for data intervals of one eroday.  \nThe sky is divided into 4700 non-overlapping, unique areas ("sky tiles"), organized into 61 equatorial declination zones (see Fig. 2 for a visual representation of the adopted sky tiling). The size of these sky tiles is exactly three degrees in declination and approximately three degrees in right ascension (on the side facing the equator), resulting in an average area of about 8.78 square degrees per tile. Each unique sky tile is embedded in an overlapping, square sky field of size 3 . 6 · × 3 . 6 · centred on the sky tile. The minimum overlap between neighbouring sky fields, introduced in order to avoid edge e ff ects in the source detection, ranges between 15 \' for polar fields and 18 \' for equatorial fields 7 . After the completion of the event calibration, the event data of each TEL chain are sorted into these overlapping sky fields, based on the right ascension and declination of each photon (tasks telgti , telselect , telstage ). Source detection and further source-level analysis are then performed in each sky field.', '4.2. EXP chain': "The event data collected in each sky tile are merged on a 'per TM' basis (task expmerge ) and sky pixel event coordinates with pixel size 0 '' . 05 centred on the sky tile centre are computed (task  \nradec2xy ). Combining these TM-specific files, filtered event files and images that include all TMs (pixel size 4 '' , 3240 × 3240 pixels) are created in a variety of energy bands suitable for source detection (task evtool ). The event filtering excludes invalid patterns (i.e. pixel patterns with a low probability of having been caused by a single X-ray photon), events on or close to bad pixels, as well as events outside a circular detection mask of radius 0.516 · . In addition, good time intervals free of background flares are created, if necessary (task flaregti ) 8 .", '4.3.1. Exposure maps': 'Exposure maps are created on a tile-by-tile basis by the task expmap . For each energy band and for each TM, a map of the active pixels of the CCD (that is, excluding flagged bad pixels, their orthogonal neighbours, as well as out-of-FOV pixels) is first divided by a vignetting map. The latter is generated as a weighted mean of the energy-dependent vignetting function across the respective energy band, the weighting being a power law of spectral index 9 Γ = 1 . 7. The vignetted CCD image is then projected onto the sky (using the same projection as the template image supplied to the task) at a series of time samples (nominally one per second) of the TM attitude, using only samples during GTIs for that TM. The resulting TM-specific sky exposure maps are then combined in a weighted mean for each band, the weights (corresponding to the fraction of area supplied by each TM with respect to the nominal all-TM area) being read from the calibration database. The results are N maps in the N specified energy bands, each giving the spatially varying mean (vignetted) exposure in that given band, in seconds.  \nAs an example, Fig. 3 shows the combined all-sky exposure map (where darker regions mark longer exposures) for the 0.6-2.3 keV band. The main qualitative features of this map are determined by the scanning law and the associated event milestones discussed above (i.e. scan interruptions, overlaps, etc.); these features are common to all exposure maps. Quantitative measures, however, such as the overall normalisation of the exposure, may di ff er for di ff erent energy bands, mainly due to the variations of the vignetting function of the telescope at di ff erent energies.', '4.3.2. Source detection': "The eSASS source detection pipeline was applied to each of the 3 . 6 · × 3 . 6 · tiles with the following steps:  \n- 1. Creation of masks to define image regions with valid data (task ermask );\n- 2. Initial source finding (task erbox );\n- 3. Background determination (task erbackmap );\n- 4. Search for candidate sources (task erbox );\n- 5. Creation of PSF fitting catalogue (task ermldet );\n- 6. Forced PSF photometry in additional energy bands (task ermldet );\n- 7. Aperture photometry in all energy bands and the corresponding sensitivity maps (task apetool ); the aperture is chosen as a fixed one at 75% encircled energy fraction 10 ;  \n8. Creation of sensitivity maps (task ersensmap ).  \nSteps 3. and 4. were iterated twice for an improved separation of sources and background. The source detection setup is very similar to the configuration used for the creation of the eFEDS catalogues (Brunner et al. 2022), where the individual detection tasks are described in more detail.  \nReliable detection and characterisation of X-ray extended sources is crucial for any subsequent application and analysis of galaxy cluster samples. It should be noted that for both point sources and extended sources the rate, count and flux values reported in the eRASS1 catalogues are based on the scaling of the best fit PSF or PSF-folded extent model. By definition, these quantities correspond to values integrated to infinite radii while the model fits are performed on circular sub-images of 1 ' radius. Therefore, the flux related source parameters may be subject to significant systematic uncertainties, in particular for large extended sources. In Bulbul et al. (2024) the interested reader will find a detailed description of the procedure used to define clean clusters samples starting from the extended sources catalogues, and also to derive robust physical quantities from the X-ray data alone.  \nIn each sky tile two catalogues were created: the single-band catalogue with sources detected in the broad 0.2-2.3 keV band (hereafter '1B') for maximal sensitivity, and a 3-bands catalogue (0.2-0.6 keV, 0.6-2.3 keV, 2.3-5.0 keV, hereafter '3B'), for which detection is carried simultaneously and single-band and combined likelihoods for each source are computed.", '4.3.3. Sensitivity maps: ersensmap': 'For both the 1B catalogue and the 3B catalogue sensitivity maps were calculated with the task ersensmap . As described in Brunner et al. (2022), the task ersensmap uses the eSASS exposure maps and background maps to estimate the detection limits in eROSITA observations. For eRASS1, the maps for the master catalogue contain the 0.2-2.3 keV source flux required to reach a typical detection likelihood of 5.0 at the respective pixel position for consistency with the threshold of the PSF fitting catalog created by ermldet . The map values for the 3-band catalogue correspond to the respective 0.2-5 keV flux. For the conversion between the map fluxes and count rates in the detection images the energy conversion factors (ECFs) listed in Table 3 were used. Based on spectral analysis of eFEDS AGN (Liu et al. 2022c), typical eROSITA-detected AGN have a median power-law slope of 2.0. So we adopted this slope and an absorbing column density N H = 3 × 10 20 cm -2 (typical value of the Galactic absorption) to calculate the ECFs (Brunner et al. 2022). The impact of the assumed power-law slope is negligible here.', '4.3.4. Sensitivity maps: apetool': 'The apetool sensitivity maps use aperture photometry to determine the selection function of a sample of X-ray point sources detected in a given spectral band. The selection function is defined as the probability of detecting an X-ray point source with a given count rate or flux in the band of interest across the eROSITA field of view, across all observations of the source throughout the duration of eRASS1. In generating these sensitivity maps, the stacked PSF at a given sky position is the exposuretime weighted superposition of the individual PSFs from all the  \nFig. 3. E ff ective (vignetted) eRASS1 exposure map (Galactic coordinates, Aito ff projection). The values in the map show the exposure time multiplied by the average of the ratio of the (vignetted) e ff ective area to the on-axis e ff ective area in the energy band 0.6-2.3 keV. E ff ective exposure values range from ∼ 100 s at the ecliptic equator to more than 10 000 s close to the ecliptic poles (not visible on the colour scale of this image). The eROSITA-DE (western) Galactic hemisphere is on the right of the central meridian in this map.  \n<!-- image -->  \nTable 3. ECFs used to calculate the 1-band and 3-band sensitivity maps.Notes. For the 3B maps the ECFs refer to the conversion between the individual bands count rates and the total 0.2-5 keV flux. These ECFs are for conversion between an observed count rate and an observed flux, not corrected for absorption, under the spectral model assumption.  \npixels (detector coordinates) that contribute to that sky position as eROSITA slews across the sky. The apetool sensitivity maps can only be used in combination with the aperture photometry of individual X-ray sources provided in the eROSITA catalogues. We refer to the appendices of Brunner et al. (2022) for a full description of the apetool functionality and Georgakakis et al. (2008) for details on the calculation of X-ray sensitivity maps based on aperture photometry.  \nThere are two key parameters that are relevant for the apetool sensitivity maps: the radius of the circular aperture adopted for photometry and the Poisson probability that the observed counts within an aperture are produced by random fluctuations of the background (Poisson false detection probability, or False Alarm Probability, FAP). The lower the latter, the less likely it is that the counts within an aperture are produced by the background, thereby suggesting the presence of an astrophysical source. For a given background level (as specified in the background maps), we use Poisson statistics to estimate the minimum number of photons within the aperture, N min so that the corresponding FAP is below an adopted threshold, thus yielding X-  \nray source detections to a given confidence level. The apetool sensitivity map is an image of N min across the field of view of the eROSITA observations. The sensitivity map can be further combined with the eROSITA exposure and background maps to determine the probability of detecting a source with FAP ≤ P thresh integrated over the eROSITA field of view. This probability as a function of count rate or flux (once an ECF is adopted) is referred to as the sensitivity or area curve and is also provided by apetool . This allows us to study the flux limit at each position, based on the local background level and exposure depth. Two applications of the apetool sensitivity maps (computation of flux limits and number counts) are described in Sects. 5.4 and 5.5.', '4.4. SOU chain': 'Source products (spectra, background spectra, response matrices, ancillary response files and light curves) were created using the srctool task, for the subset of bright sources with a detection likelihood greater than 20 in the single-band 1B catalogue. In terms of net counts, this corresponds to a sample with a median number of ∼ 21, and 90 th (10 th ) percentile of ∼ 79 ( ∼ 12) counts. Further details about the source products are given in Sect. 7.2.', '4.5. Bright sources: pile-up and other losses': "Source parameters derived from the detection pipeline, such as total flux or temporal or spectral variability, carry systematic uncertainties, in particular for very bright X-ray sources. This is mainly due to the pile-up e ff ect (Ballet 1999; Davis 2001; Tamba et al. 2022), which occurs if two or more photons hit the same CCDpixel in the same (50 ms) read-out cycle and the sum of the charges created will enter the event analyser ('energy pile-up'). Pile up also occurs if these photons are recorded in two adjacent pixels where, after recombination of the individual charges, a higher energy value is reconstructed ('pattern pile-up'). It is im-  \nportant to note that for pile-up the total energy band is relevant and also events below the lower event trigger threshold, including optical photons. The probability for these e ff ects mainly depends on the source brightness in photons / cycle / pixel and on the actual shape of the PSF, for example, it is reduced for a deteriorated PSF far o ff -axis. This can lead to apparent spectral and temporal variability within a scan through the FOV, but also between one (e.g., more central) scan and another (more o ff -axis). All this complexity in general requires detailed simulations for an accurate estimate of pile-up e ff ects (see e.g. König et al. 2022). Preliminary studies using SIXTE (Dauser et al. 2019), indicate that for eRASS1 pile-up starts to become important (more than a few per-cent e ff ect) for point sources brighter than ≈ 10 -11 erg s -1 cm -2 in the 0.2-2.3 keV band.  \nPhoton counts can also be lost in eROSITA due to e ff ects in the event analyser and telemetry limitations between camera electronics and ITC. The 'event quota' mechanism for each camera ensures a reasonable telemetry rate, even in the case of, e.g., a CCD column becoming bright between two ground contacts. In the current implementation, the event quota is triggered if in one TM there are more than 50 events for four consecutive readout cycles after onboard rejection of minimum ionising particles (MIP) and bad pixels. In that case, for one minute only frames containing less than 50 events are telemetered. After one minute this is reset and all frames are transmitted until the trigger criteria are fulfilled again. Due to this event quota, complete read-out frames are lost, and all sources within the FOV are a ff ected, differently from pile-up, which only a ff ects the piled-up source itself. This mechanism was designed for instrumental reasons, but very bright point sources, such as Sco X-1 (or even bright optical stars) and extended sources (for example Puppis A) are also (celestial) triggers. Missing read-out cycles are properly handled in the exposure computation, but there remains a bias towards fainter read-out cycles during active event quota triggers.  \nFinally, the source parameters of some of the brightest sources in the catalogues ( ML\\_CTS\\_1 > 1000) su ff er from poor convergence of the PSF fits due to the large number of individual events to be included in the modelling. This may result in larger than expected deviations in both photometry and astrometry.", '4.6. Other known data processing issues': "During the commissioning phase of eROSITA, it was noticed that TM5 and TM7 (the CCDs not equipped with an on-chip Al optical blocking filter) were contaminated by optical light: a small fraction of sunlight reaches the CCD, by-passing the filter wheel. The intensity of this optical contamination depends on the orientation of the telescope with respect to the Sun. This 'light leak', mostly restricted to very low photon energies (typically < 0 . 3 keV), generates a non-negligible amount of telemetry data and decreases the low energy coverage and spectroscopic capabilities for these two cameras. In order to reduce the amount of transmitted data from TM5 and TM7, their primary thresholds are higher (about 125-145 eV) than the ones in the other TMs (about 65-95 eV). The fact that the amount of contamination by optical light is spatially and temporally variable makes it very di ffi cult to derive an accurate energy calibration for these TMs. More details about the light leak in eRASS1 can be found online on the DR1 web portal 11 .  \nThe e ff ect of optical loading, i.e. the appearance of fake Xray sources, or the distortion of the measured X-ray properties  \nof sources associated to (very) bright optical stars, is discussed in § 5.2.4.  \nIn some cases, the source detection algorithm failed to converge during its error estimate procedure, leaving some sources without reliable uncertainty estimates for the position, counts, or extent. Many of the a ff ected sources have low detection likelihood or are related to spurious detections in areas of extended emission. In some cases, however, also significant detections can be a ff ected by this issue. We flag these sources with the labels FLAG\\_NO\\_RADEC\\_ERR , FLAG\\_NO\\_CTS\\_ERR and FLAG\\_NO\\_EXT\\_ERR , respectively (see Tables 5 and 6).  \nFinally, in calculating the number of seconds elapsed between the time reference datum of the mission (00:00 hrs Moscow time on January 1 st , 2000) and a given later date, leap seconds should be added. Five leap seconds occurred between the reference datum and the start of the mission. However, the pipeline software used to create the DR1 data omitted to add these seconds. Therefore, when converting UTC times in the badcamt (bad time intervals) and timeoff (instrumental one second time shifts) calibration components to spacecraft clock a five second shift is introduced, leading to the exclusion of five seconds of good data. In both cases, changes in status (i.e. date entries in the component table) tend to occur when the camera is in an anomalous state, not receiving data (see Brunner et al. 2022, Appendix B). As the cameras are not yet registering photons at these times, not considering leap seconds does not have any other negative consequences (such as inclusion of bad time interval in the data) in this case. These five leap seconds will be properly corrected in the next data release.", '5. The eRASS1 X-ray catalogues': 'Following the approach devised for the Performance Verification eFEDS survey (Brunner et al. 2022), we present here two distinct X-ray catalogues: a catalogue of sources detected in the 0.2-2.3 keV band (selected from the 1B detection process; Main catalogue) and a catalogue of sources detected in the 2.3-5 keV band (selected from the 3B detection process; Hard catalogue).  \nWe generate the single-band 1B catalogue including sources down to a low detection likelihood ( DET\\_LIKE\\_0 ⩾ 5), to maximise completeness. We then make use of the eRASS1 digital twin simulations (Comparat et al. 2019, 2020; Seppi et al. 2022) to estimate the amount of spurious detections as a function of the detection likelihood threshold (see also Liu et al. 2022c). Based on these all-sky survey simulations, we define our Main sample as the one comprising all extended sources and all point sources with DET\\_LIKE\\_0 ⩾ 6. Table 3 of Seppi et al. (2022) indicates that the Main catalogue should contain ≈ 14% spurious detections. This reduces to about 1% for DET\\_LIKE\\_0 > 10. Pointlike sources with 5 ⩽ DET\\_LIKE\\_0 < 6 are released as a (highly contaminated) Supplementary catalogue.  \nTo extract the hard sample from the 3B catalogue, we apply a threshold for the detection likelihood in the 2.3-5 keV band of DET\\_LIKE\\_3 ⩾ 12. The threshold is higher than the one applied for the Main catalogue as the lower sensitivity and higher background in that energy range significantly increases the number of spurious detection at a given detection likelihood (Liu et al. 2022c). Based on the same simulations described in Seppi et al. (2022), we estimate for this threshold a spurious sources fraction of about 8-10% in the hard-band selected sample. We note here that, as shown in Seppi et al. (2022), the expected fraction of spurious sources depends on the exposure, with lower fraction of spurious detections predicted for higher exposures. Here, for simplicity, we have adopted all-sky average estimates.  \nTable 4. Basic eRASS1 catalogue properties.  \n3B detection [2.3-5 keV only]  \nCatalogue  \nDET\\_LIKE\\_3  \nEXT\\_LIKE', 'of Sources': "Hard  \nHard, PS  \nHard, Ext.  \n⩾  \n12  \n⩾  \n12  \n=  \n0  \n5087  \n⩾  \n12  \n>  \n0  \n379  \nNotes. 'PS' indicates point sources (i.e. those with EXT\\_LIKE = 0), and 'Ext.' indicates extent-selected sources with EXT\\_LIKE > 0. We note here that there is a slight di ff erence in EXT\\_LIKE parameter in the 1B and 3B detections, as for the latter all photons in the 0.2-5 keV band are used to evaluate the source extent.  \nTable 4 presents a summary of the catalogues selection criteria and properties. In Fig. 4 we show the distributions of net counts for the Main and the Hard samples, respectively.  \nA summary of the content of each catalogue is presented in the Appendix D (column descriptions, units). Below we describe in greater detail the catalogue creation procedure, the astrometric verification steps and our attempt to flag potential spurious sources.", '5.1. Catalogue creation and preliminary astrometric correction': "The catalogues resulting from the PSF fitting with task ermldet were re-formatted using the eSASS task catprep and then merged into hemisphere catalogues. Since the survey sky tiles overlap each other, the sources outside the nominal, nonoverlapping area were removed from each sky tile catalogue. In order to avoid the loss of significant sources detected by chance just outside the nominal areas in two adjacent tiles, detections in a ± 30 '' strip near the nominal borders were matched, using a matching radius of r = 15 '' for point sources and r = 30 '' for extended sources, and only one detection for each match was kept in the merged catalogue. These matching radii approximately correspond to half of the HEW and one HEW of the survey averaged PSF (see Fig. A.1).  \nAs we discussed above, the eROSITA field of view and survey scanning strategy imply that a source near the ecliptic equator is visited about six times within a time span of one day in each of the eROSITA six-months surveys. The time span and the number of visits increase with the ecliptic latitude of the source. In the catalogues, the epoch of survey coverage was estimated for each source by using the attitude time series for camera TM1 to calculate the times of the observation closest to the optical axis as well as the first and last appearance of the source in the camera field of view. These epochs are listed in the columns MJD , MJD\\_MIN , and MJD\\_MAX .  \nThe positional uncertainty of X-ray sources is an important parameter for their association with multi-wavelength counterparts, especially given the relatively large PSF of eROSITA. The ermldet task provides a statistical estimate of this quantity for individual X-ray sources based on the spatial distribution of their X-ray photons and a PSF model (see Brunner et al. 2022, Appendix A.5, for more details). These measurements may under-  \n⩾  \n0  \n5466  \nestimate the true positional errors because of, e.g., calibration uncertainties and other systematic e ff ects (e.g. Webb et al. 2020). In particular, the astrometric accuracy of the eROSITA all-sky catalogues can be a ff ected by the following systematic factors:  \n- -Errors in the timing between attitude measurements and event arrival;\n- -Boresight calibration errors;\n- -Systematic errors introduced by the PSF fitting.  \nAs for the latter point, eRASS1 data, just like those of the CalPV phase, have been analysed using the PSF model derived from on-ground calibration. The analysis presented in Appendix A demonstrate that this e ff ect is negligible, so we discuss only the former two below.  \nDue to the SRG scanning geometry, a timing mismatch between the attitude measurements and event timing would result in an o ff set along the scanning direction, i.e., ecliptic latitude β . An o ff set in the boresight between the SRG attitude solution and the eROSITA cameras would result in a constant astrometric o ff set in ecliptic coordinates during a 180 · scan between the two ecliptic poles. Assuming that any timing or boresight o ff sets vary only slowly with time, in order to correct for this e ff ect we divided the merged catalogue into stripes of 1 · width in ecliptic longitude λ , corresponding to ∼ 1 day of survey scanning. For each stripe the X-ray positions of point-like detections at ecliptic latitudes ( -60 · to + 60 · ) were matched with mid-infrared counterparts from the AllWISE catalogue (Cutri et al. 2021). After applying a colour cut (0 . 3 mag < W 1 -W 2 < 1 . 7 mag and 2 . 2 mag < W 2 -W 3 < 3 . 6 mag) to select for likely QSO matches, median values of the o ff sets β X -β IR and ( λ X -λ IR) × cos( β ) were calculated. All X-ray positions within each latitude stripe were then corrected using the median o ff sets in ecliptic longitude and latitude ( ∆ λ × cos( β ), ∆ β ). The applied o ff sets range between -4 . 0 '' and + 3 . 9 '' and -1 . 5 '' and + 2 . 5 '' , respectively.  \nThe resulting statistical errors are given in each coordinate as upper and lower bounds (columns RA\\_LOWERR , RA\\_UPERR , DEC\\_LOWERR , DEC\\_UPERR ). An error estimate averaged over both dimensions and directions is given as  \nσ RA = (RA\\_LOWERR + RA\\_UPERR) / 2  \nσ DEC = (DEC\\_LOWERR + DEC\\_UPERR) / 2  \nRADEC\\_ERR = p ( σ RA) 2 + ( σ DEC) 2 , (1)  \nin line with other X-ray catalogues (e.g. Webb et al. 2020).  \nThe upper panels of Figure 5 displays the distribution of RADEC\\_ERR as a function of DET\\_LIKE\\_0 for the point sources in the Main and Hard catalogues, respectively. For sources where the calculation of RADEC\\_ERR failed (see Sect. 4.6), we calculate RADEC\\_ERR using the empirical correlation extracted from Fig. 5. It should be noted here that RADEC\\_ERR does not represent the 68% error radius for two parameters. Under the assumption of a circular error region, the averaged 1-dimensional 68% error as required e.g. for the comparison with a Rayleigh distribution can be derived with σ = RADEC\\_ERR / √ 2. We further elaborate on the astrometric accuracy of the eRASS1 catalogue in Sect. 6.2, where we present a validation method based on a comparison with external catalogues, which reveals the extent of the systematic uncertainty beyond the statistical one described here.  \nThe lower panels of Fig. 5 display the distribution of flux significance as a function of detection likelihood for the point sources in the main and hard catalogues. For the large majority of the sources, the flux measurements have large uncertainties; in order to have at least 3σ flux measurements in the Main catalog, one could adopt an approximate threshold of DET\\_LIKE\\_0 > 20.  \nFig. 4. Distributions of net counts for the Main (0.2-2.3 keV; Left panel ) and Hard (2.3-5 keV; Right panel ) catalogues. Point sources ( EXT = 0) and extended sources ( EXT > 0) are plotted in blue and orange, respectively. Point sources with any SP flag (see Table 6) are displayed as blue shaded histograms. The sources that appear in both the Main and the Hard catalogues are plotted in red in the left panel.  \n<!-- image -->", '5.2. Flagging of problematic sources': 'The imperfect nature of the source detection process inevitably leads to contamination of the eRASS1 catalogue by spurious sources and / or inaccuracies in the derived source properties. The most clearly identifiable examples of spurious detections can be found within the vicinity of extremely bright X-ray point sources, such as Sco X-1, or bright, large extended sources, like supernova remnants, nearby galaxies, or galaxy clusters (see Fig. 6), whereas less clear-cut cases can be found at the lowest detection likelihoods, and their contribution to the eRASS1 catalogue quantified via simulations. Optical loading of the CCDs could also introduce fake X-ray sources in the catalogue.  \nTo warn users of a potential spurious origin for a detection, despite a possibly high-detection likelihood, we have flagged sources that are located within overdensities in the eRASS1 source catalogue associated with systems that could create problems for the automatic background estimation in the detection pipeline (supernova remnants, extremely bright X-ray point sources, Galactic star clusters, local galaxies, and galaxy clusters). We also flag catalogue entries matched to very bright optical stars, as we describe below.', '5.2.1. Identification of overdense regions': "In order to identify those regions on the sky where many potentially spurious sources are clustered, we performed an empirical search for regions with a suspiciously large number of detected sources compared to their surroundings: after performing a uniform cut on the detected flux at F 0 . 2 -2 . 3 keV > 5 × 10 -14 erg s -1 cm -2 , to reduce dependence on the spatially varying exposure, we computed a density map of point-like and extended sources in the single-band catalogue, using a pixelisation of 0 . 25 deg 2 . A 'background' source density map was then created by applying a median filter with a radius of 10 · . By comparing the two maps, we identified all regions with a local source density more than twice the background. The exact shape of the overdensities was extracted by 'zooming in' on each identified region, creating a smoothed histogram of the local source distribution, and selecting all regions with a density larger than three  \ntimes the local background which contribute more than 20 excess sources.  \nWhile this procedure yields many overdensities caused by truly spurious source excesses, some resulting regions are expected to correspond to accumulations of real astrophysical sources. We thus manually identified the correspondence of each of the ∼ 80 localised overdensities to astrophysical sources, using the SIMBAD database (Wenger et al. 2000). Our overdense regions were then classified, and the enclosed sources flagged, according to their correspondence to i) di ff use emission associated with known supernova remnants ( FLAG\\_SP\\_SNR ), ii) excess emission in the vicinity of extremely bright point sources ( FLAG\\_SP\\_BPS ), iii) Galactic star clusters ( FLAG\\_SP\\_SCL ), iv) nearby galaxies ( FLAG\\_SP\\_LGA ).  \nThis classification might be useful, for instance, if one were interested in studying the population of Milky Way point sources, as one would likely want to mask spurious sources caused by mis-classified di ff use emission from supernova remnants, but might not want to mask Galactic stellar clusters.", '5.2.2. Galaxy cluster catalogues': "Another flag is applied for possible spurious sources which are located close to known galaxy clusters. To do that, we make use of published X-ray cluster catalogues, including MCXC (Piffaretti et al. 2011), XXL365 (Adami et al. 2018), XCS (Mehrtens et al. 2012), eFEDS (Liu et al. 2022a), and X-CLASS (Clerc et al. 2012). Sources lying within R = 0 . 5 × R 500 from the cluster center are flagged as FLAG\\_SP\\_GC\\_CONS . When R 500 is not provided in the published catalogue, we adopt a radius of R = 500 kpc. We note that, to avoid over-counting, optical cluster catalogues are not used in this step, because large o ff sets between clusters' X-ray centers and optical centers are frequently observed (see, e.g., Seppi et al. 2022, 2023). Cluster catalogues selected on the basis of the Sunyaev-Zeldovich (SZ) e ff ect are not included either, due to the relatively large positional uncertainties in SZ surveys. Therefore, the catalogue of known galaxy clusters we used in the above approach is a rather conservative and incomplete compilation, and the identified spurious sources should be considered as a supplement for the identifi-  \nFig. 5. Distributions of RADEC\\_ERR (upper panels) and flux measurement significance (in terms of ML\\_RATE\\_1 / ML\\_RATE\\_ERR\\_1 , lower panels) as a function of detection likelihood for point sources in the eRASS1 single-band detected catalogue (left panels) and in the Hard catalogue (right panels). In the left panels, the 40%, 68%, and 95% contours are plotted in orange, while the cyan line indicates an empirical correlation that describes the mode of the distribution, as reported in the label. Sources from the Main catalogue are shown in blue, those from the Supplementary catalogue in grey. The right panels ( Hard catalogue) have only the 68% contour plotted. Note that the lower right panel displays the likelihood and flux significance in the 2.3-5 keV band.  \n<!-- image -->  \nTable 5. Spurious and problematic source flag description.  \nation of overdensities. Further cleaning is performed only for galaxy cluster candidates in the extended source catalogue. An extended source is flagged as a possible spurious detection when it is too close to its neighbour. Visual inspections are also performed on the extended source catalogue to remove any remaining obvious spurious sources and correct for mis-classified cases. We refer the readers to Bulbul et al. (2024) for more details of the cleaning procedure we performed in the extended source catalogue.  \n5.2.3. Summary of spurious sources flagging procedure  \nThe flagging procedure described above identifies cases where high local background levels render the automatic pipeline detection process unreliable. This is illustrated in Fig. 7, which shows the background rate distributions for sources with and without flags: the flagged sources are mostly found in regions with enhanced background rate.  \nFig. 6. Top panel : Aito ff projection of the eRASS1 1B catalogue in Equatorial coordinates (J2000), with each grey point representing a detected source within the catalogue, and the coloured points denoting sources that have been flagged as potentially spurious according to the scheme presented in Table 5. Darker stripes are due to larger sources density due to a higher exposure in those parts of the sky. Bottom left : Zoom in plot of sources detected within the vicinity of Sco X-1, with pink sources flagged as potentially spurious using the FLAG\\_SP\\_BPS column in the catalogue. Bottom right : Similar for sources within the vicinity of the Vela SNR.  \n<!-- image -->  \nA summary of the di ff erent spurious sources flags added is presented in Table 5, while Table 6 reports the number of sources for each of the flag categories. After removing all those sources that are flagged by at least one of the potential spurious categories, the Main catalogue contains 890 036 point sources, with a median sky density of approximately 37 deg -2 .  \nA note of caution is in place here: although steps are taken here to greatly reduce the number of high-detection-likelihood contaminants, it is still recommended that users double-check the relevant science images for their sources of interest before publication (e.g., if their point sources lies within a likely galaxy cluster but is not flagged as spurious here).", '5.2.4. Optical loading': "The eROSITA detectors are prone to optical loading, i.e. the accumulation of low energy (optical / UV) photons within a CCD pixel over the frame integration time of 50 ms, whereas the relevance of the e ff ect is governed by the optical brightness and detector position of the respective source. If the summed energy exceeds the X-ray detection threshold, any optically bright source will start to generate false X-ray events. Due to the na-  \nTable 6. Number of sources flagged as potentially spurious.  \nNotes. 'PS' stands for Point Sources, 'Ext.' for extent-selected, and we separate the Main (' M ') and the Hard (' H ') catalogues. The 'Any SP Flag' row indicates the number of sources that are flagged by any of the five identifiers of potential spurious sources in overdense regions, SP\\_SNR , SP\\_BPS , SP\\_SCL , SP\\_LGA , and SP\\_GC\\_CONS . The last three lines give the number of sources for which statistical error estimate from the pipeline failed (see Sect. 4.6 for more details).  \nFig. 7. Top panel: distribution of local background rate level ( ML\\_BKG\\_1 / ML\\_EXP\\_1 , in units of counts arcmin -2 s -1 ) for sources with and without flags and DET\\_LIKE\\_0 > 10. Bottom panel: the fraction of each corresponding sub-sample among all the sources above any background level.  \n<!-- image -->  \nture of the e ff ect, these events appear predominantly at lower Xray energies. The strength of the optical loading signal strongly increases with increasing source brightness; however, for very bright sources 'saturation-like' e ff ects may occur, when events are removed by the pattern filtering procedure. The optical loading signal further depends on the o ff -axis angle of the source, as the optical PSF is sharper close to FOV center and photons are focussed on a smaller detector area. i.e. fewer pixels. During eRASS data taking, any optically bright source passes several times on di ff erent scan paths through the FOVs and thereby a time dependent optical loading signal may be generated. If a su ffi cient number of X-ray events is created, source detection triggers and the object makes it into the eRASS catalog. Due to the required intrinsic brightness, stars are by far the main contributors to optical loading sources in our data.  \nAmong the practical consequences for astrophysical studies are the presence of fake X-ray sources, pseudo X-ray variability and the distortion of real X-ray sources, as the source may of course be optically bright and also an intrinsic X-ray emitter. If the optically bright source is X-ray dark, or at least faint enough to fall well below the eRASS sensitivity limits, all the registered events would come from optical loading. These fake X-ray sources always have very soft spectra and show pseudo variability. If the optically bright sources are also X-ray bright, these characteristics are likewise present, but the detected signal is a mixture of true X-ray photons plus optical loading events and contamination e ff ects. Any potential disentangling between these e ff ects depends on the individual source properties and on the detailed science objectives, but in general the X-ray properties of sources contaminated by optical loading are uncertain.  \nTo best characterise the optical loading e ff ect, X-ray dark sources are required and two suitable stellar populations are used here: main-sequence and mildly evolved stars with spectral types late-B to mid-A and red giants, i.e. K / M stars with luminosity class III-I. Cross-matching these with the eRASS1 catalogs shows that optical loading is expected to a ff ect DR1 sources if they exceed certain brightness thresholds, and likely a ff ected X-ray sources are flagged in the catalog. The adopted brightness limits are B, V, G ≤ 4 . 5 mag and J ≤ 3 mag; if one or more criteria are fulfilled, then the source is tagged with FLAG\\_OPT . In total, 750 (17, 14) sources are flagged as potentially contaminated by optical loading in the Main ( Supplementary , Hard ) catalogues, respectively. Input for cross-matching are the Tycho-2 (Høg et al. 2000), 2MASS (Cutri et al. 2003), Gaia DR3 (Gaia Collaboration et al. 2023) catalogues plus Simbad database. Where proper motion information is available, optical catalog entries are updated to epoch 2020 positions, and a uniform matching radius of 15 '' is used. A full treatment depends on the individual source properties and is beyond the scope of this work, but users should be aware that specific attention is required when dealing with eROSITA data of brighter stars.", '5.3. Association of soft and hard band selected sources': 'As discussed above, we performed source detection with two di ff erent settings: a single X-ray band (1B, 0.2-2.3 keV) and a three-band detection (3B, 0.2-0.6, 0.6-2.3, and 2.3-5 keV), subsequently down-selected based on the 2.3-5 keV band significance. These produce, respectively, soft and hard X-ray selected source samples. The data used in the 1B and 3B detection are nonetheless largely overlapping. The only photons that are included in the 3B detection but not in the 1B are those in the 2.35 keV band, where the instrumental e ff ective area is relatively low and the background relatively high. For this reason, most of  \nthe Hard catalogue sources are expected to have a matching entry in the 1B catalogues. This was the case also in the eFEDS survey, where 90% of hard band sources were found to have a counterpart in the main, soft X-ray selected catalogue (Nandra et al., in prep.). It is nonetheless of particular interest to identify those sources that are only detected in the hard band, as they signpost objects with extremely hard spectra, most obviously due to heavy obscuration. The association between sources in the two catalogues is not straightforward, however. The complexity arises due to various factors, such as the positional uncertainties, the morphological classification (extent measurement) and blending with nearby objects. In this section, through a specific matching procedure between the 1B and 3B sources, we provide a \'weak\' association that can be used to select sources only detected in the hard band and a \'strong\' association that can be used to select entries in the X-ray catalogues with a high degree of confidence that the X-rays originate from the same astrophysical source.  \nThe weak association is based only on positional information. Since for extended sources the value of the extent parameter EXT (i.e. the best-fit core size of the beta model fitted to an extended source) is very broadly distributed, ranging from ∼ 10 \'\' to the maximum allowed value of 60 \'\' , to search for a counterpart to those sources we adopt a large matching radius of four times the EXT value. For point sources ( EXT = 0), we simply adopt a matching radius of 10 \'\' , which is 99 percentile of the point source RADEC\\_ERR in the main catalog. For each 1B source, we search for 3B sources within its matching radius, and for each 3B source, we search for 1B sources within its matching radius. Then we merge the results, which is equivalent to adopting the larger radius among the two. As a result of such loose association criteria, one source could be associated with multiple sources, some of which are false matches and some of which have di ff erent morphological classifications. To classify these matches, we define P2P, E2E, P2E, and E2P matches to denote pairs of point source (P) and extended sources (E); the former and the latter letters (P or E) indicate the classifications in the 1B and 3B catalogues, respectively. Hard band sources which do not have any counterpart in the soft band catalogues within the match radius are designated as hard band only sources.  \nThe \'strong\' associations require more strict criteria. First, we require that each pair of sources have the same morphological classification, i.e., we adopt only E2E or P2P pairs. With E2E associations, one source can be matched to multiple ones. We adopt only the nearest match in such cases, so that the involved sources are unique. In the cases of P2P matches, one source is never matched to multiple counterparts, but in addition to a maximum separation of 10 \'\' , we further require i) separation / error < 3, where error is the larger of the two RADEC\\_ERR , and ii) the 1B and 3B measured 0.2-2.3 keV (combining band 1 and 2 in the 3B case) source counts cannot di ff er by > 50%. In this way, we enforce not only that the matched sources are su ffi ciently close in sky position, but that also have a consistent broadband brightness, such that sources with the same position but di ff erent extents (caused by de-blending uncertainties in crowded regions) can be excluded. We store the counterpart unique source ID ( UID ) in the UID\\_Hard column of the Main and Supplementary catalogues and the UID\\_1B column of the Hard catalogue as indicators of strongly associated sources.  \nFor the weak associations, we also store the counterpart UID in the UID\\_Hard or UID\\_1B columns but multiplying it by -1. In this way, catalogue columns UID\\_Hard > 0 or UID\\_1B > 0 indicate that a source has a strong association in the other catalogue, while catalogue columns UID\\_Hard < 0 or UID\\_1B < 0  \nFig. 8. Distributions of 0.6-2.3 keV (x-axis) and 2.3-5 keV (y-axis) count rates of the hard-band selected sources. The "Hard-only" sources are marked in red colour. The solid line indicates the ratio (0.084) that corresponds to an unabsorbed power-law spectrum with a slope of 1 . 7. The dashed line indicates a 1:1 ratio.  \n<!-- image -->  \nindicate that a source may have a (weak-association) counterpart in the other catalogue. Finally, UID\\_1B = 0 marks the Hard sources that have no counterpart within their matching radius in the 1B catalogues. Following this procedure, we find 780 \'Hardonly\' sources; of these, 764 are point-like ( EXT = 0) and 16 extended ( EXT > 0). Figure 8 displays the count rate distribution of the hard-band selected sources in the 0.6-2.3 and 2.3-5 keV bands. The Hard-only sources identified above (red points) have significantly higher hardness than the others, by construction.  \nThe matrix of possible identifications of Hard sources in the 1B catalogues, following our association criteria, is given in Table 7. For point sources, about 84% have a strong association, 15% are Hard-only and just 1% have weak associations in the 1B catalogues. These fractions do not change if one considers only sources without any spurious flag. Among extended sources in the Hard catalogue, on the other hand, a larger fraction is flagged as potentially spurious. Of the remaining ones (136 in total), about 80% have a strong association with a 1B catalogue extended source, about 9% have a weak association in the 1B catalogues, and 11% (15 in total) are Hard-only. Such Hard-only extended sources are most probably mis-classified point sources; indeed, only 1 / 15 of these has extent likelihood EXT\\_LIKE > 6 and extent EXT > 30 \'\' .', '5.4. The flux limit and completeness of eRASS1 in the 0.5-2 keV band': "To quantify the flux limit of the eRASS1 survey in the commonly adopted energy band 0.5-2 keV, we use the aperturephotometry-based method based on the apetool task (see Sect. 4.3.4). A local flux limit can be defined based on the local sensitivity curve as the flux where the probability of collecting a su ffi cient number of source photons to constitute a detection at the given threshold reaches a particular value, for example 50%. To calculate a sensitivity curve as a function of flux, we need an ECF to convert count rate to flux. Considering the nonuniform Galactic absorption, we assume a power-law spectral model with Γ = 2 and the total Galactic absorption column density (Willingale et al. 2013) at each source position, and calculate two ver-  \nTable 7. Summary of identification of Hard (3B) sources in the 1B (Main and Supplementary) catalogues.Notes. FLAG = 0 selects sources with all the five spurious flags in the Hard catalogue being 0. Here, P2P, E2P, P2E and E2E refer to possible combinations of matched pairs between the 1B and 3B catalogs, based on the morphology of the sources: point-like (P) or extended (E); the former and the latter letters (P or E) indicate the classifications in the 1B and 3B catalogues, respectively. See text for more details.  \nFig. 9. Bottom panels: Hammer-Aito ff projection maps, in Galactic coordinates, of the logarithm of the 0.5-2 keV flux limit calculated as the flux at 50% sky covering fraction (see text for details). The left and right panels correspond to values uncorrected and corrected for Galactic absorption, respectively. Six particular regions where the flux limits are increased by bright X-ray sources are marked with red circles; from the north (top) to the south (bottom) they are: the Virgo cluster, Sco X-1, the Vela SNR, the Crab pulsar, the LMC, and the SMC. Top panels: histogram of the logarithm of the flux limit in erg s -1 cm -2 . The colour bar on the X-axis illustrates the intensity scale of the corresponding map in the bottom panels.  \n<!-- image -->  \nions of ECF for each source, one correcting the 0.5-2 keV flux for absorption and the other not.  \nForced photometry results at source positions in the 0.5-1 (band P2) and 1-2 keV (band P3) bands are available in the catalogue and can be combined into the 0.5-2 keV values. We sum the total counts ( APE\\_CTS\\_P2 + APE\\_CTS\\_P3 ) and background counts ( APE\\_BKG\\_P2 + APE\\_BKG\\_P3 ) to calculate APE\\_CTS\\_S and APE\\_BKG\\_S , where the su ffi x 'S' (Soft) indicates the 0.52 keV band. To combine the vignetted exposure times in P2 and P3, we measure the weighted average as  \nAPE\\_EXP\\_S = (1 + w ) / 1 APE\\_EXP\\_P2 + w APE\\_EXP\\_P3 ! ,  \nwhere w is the relative weight factor, which depends on the ECF adopted to create the exposure maps in the P2, P3 and S bands. By comparing these three maps, we find a proper factor of w = 0 . 6268. The distribution of sources in the space of background counts and exposure time shows a strong concentration along a linear correlation, which indicates the typical  \nbackground ( APE\\_BKG\\_S / APE\\_EXP\\_S ∼ 0 . 003 counts / s). There is a small fraction ( ∼ 1%) of obvious outliers below 0.002 counts / s, which indicates an underestimated background in low S / N regions. We adopt a minimum value of 0.002 counts / s to correct them. Having computed APE\\_CTS\\_S , APE\\_BKG\\_S , APE\\_EXP\\_S , and the ECF, we can calculate the Poissonian probability ( APE\\_POIS\\_S ) of each source and the sensitivity curve determined by the local background using the Python package 'scipy.special' 12 , which provides the same function as in the apetool task. The values of APE\\_POIS\\_S are reported in the Main and Supplementary catalogues, so that a sample with a fixed FAP threshold can be defined by the user as needed.  \nWe here adopt an aperture size for photometry corresponding to a radius that includes 75% of the PSF photons (i.e. encircled energy fraction of 0.75), while the Poisson false detection probability threshold is set to P thresh = 4 × 10 -6 (following Georgakakis et al. 2008). To project the flux limits to a map,  \nwe subdivide the sky into HEALPix 13 pixels (Zonca et al. 2019; Górski et al. 2005), adopting a HEALPix resolution order of 6. Through Voronoi tessellation, we also pixelise the sky into cells, each of which contains only one source. In each HEALPix pixel, we average the sensitivity curves of all the sources and weigh each source by its Voronoi cell area. Adopting the 50% flux of the averaged sensitivity curve, we obtain a flux limit for each HEALPix pixel and thus a flux limit map of the hemisphere.  \nWe create two distinct maps, as displayed in Fig. 9. The left panel shows the flux limit of a (galactic) X-ray point source without consideration of the galactic absorption. The flux limit map follows the exposure map and the di ff use background map well. The right panel shows the case of an (extragalactic) Xray source obscured by the total Galactic N H. In the absorption corrected case, the high N H in the Galactic plane boosts the flux limit. The hemisphere median flux limits are approximately 5 × 10 -14 erg s -1 cm -2 (before Galactic absorption correction) and 6 × 10 -14 erg s -1 cm -2 (after Galactic absorption correction), respectively. If we were to increase the probability threshold to 80%, the corresponding flux limit would increase by about 50%.  \nWe provide here some basic characteristics for two examples of sub-sample selection with well-defined statistical properties: a flux-limited one, and one above a fixed Poisson false-probability threshold. A selection of point sources obeying F 0 . 5 -2keV > 5 × 10 -14 erg s -1 cm -2 (observed, i.e. not corrected for absorption) returns a sample with 207 439 entries (after removing all those that are flagged by at least one of the potential spurious source categories), with a median sky density of approximately 9.8 deg -2 . On the other hand, a selection of point sources obeying APE\\_POIS\\_S < 4 × 10 -6 returns a sample with 229 266 entries (after removing all those that are flagged by at least one of the potential spurious source categories), with a median sky density of approximately 8.6 deg -2 .  \nFinally, an estimate of the completeness of the eRASS1 catalogue is provided by the analysis of the detailed 'digital twin' eRASS1 simulation presented in Seppi et al. (2022). There it is shown how the point sources' completeness as a function of 0.5-2 keV flux changes with increasing exposure. Based on that work, Figure10 shows the estimated completeness for point sources, expressed as the ratio of detected to simulated AGN from the analysis of the simulations as a function of the net exposure time (for the Main catalog with DET\\_LIKE\\_0 > 6).", '5.5. Number counts of point sources': 'Before eROSITA, deep extra-galactic X-ray surveys 14 were generally confined to particular regions with areas of at most tens of square degrees. The largest contiguous X-ray survey before eRASS1 is the eFEDS (Brunner et al. 2022), which covers an area of 140 deg 2 . In such (relatively) small extra-galactic regions, the point source number density can be considered uniform. Serendipitous surveys, built by combining narrow field exposures (with Chandra , XMM-Newton , Swift ) can provide even larger sky coverage, albeit very patchy. However, in eRASS1 it becomes apparent that the X-ray point source number density is nonuniform across the sky. The first-order reasons for this are obviously the inhomogeneous Galactic absorption (see e.g. Ponti et al. 2023a) and the non-uniform distribution of the Galactic X-ray population itself. Even after excluding these Galactic fea-  \nFig. 10. Estimated completeness of eRASS1 from the simulations of Seppi et al. (2022) as a function of input 0.5-2 keV flux for di ff erent exposure times. The completeness is estimated considering only point sources ( EXT\\_LIKE = 0) with DET\\_LIKE\\_0 > 6.  \n<!-- image -->  \ntures, the number density of distant AGN may not necessarily be uniform either, due to the inhomogeneous large-scale structure and the potential anisotropy (e.g., a dipole structure; see Secrest et al. 2021) of the Universe. More details about the eROSITA point-source number density maps are presented in Liu et al. (in preparation). In this section, we present the point-source number count distributions averaged in a few wide Galactic latitude ranges.  \nWe divide the hemisphere into four Galactic latitude levels, 0-10 · , 10-20 · , 20-40 · , and 40-90 · . As introduced in the previous Sect. 5.4, we calculate the point source number counts using the method described by Georgakakis et al. (2008) based on the products of apetool for the 3B catalogue in the 0.6-2.3 keV band. Through the N H-dependent ECF (see Table 3), we convert the count rate of each source into absorption-corrected flux. The uncertainties of the cumulative X-ray number counts are estimated using a bootstrap re-sampling approach. For a given eROSITA sub-sample selected within a given Galactic latitude interval, the X-ray sources are randomly drawn with replacement to generate 100 new samples with the same size as the original one. The X-ray number counts are then generated for each of the 100 samples following the same approach as with the real data and the uncertainty at fixed flux is then estimated as the 1σ rms scatter of the 100 number count realisations.  \nThe resulting cumulative number count distributions as a function of the absorption-corrected flux and the corresponding sensitivity curves are displayed in Fig. 11. The number counts obtained in sky regions with | b | > 20 · (high Galactic latitudes) are consistent with each other and with the results of Georgakakis et al. (2008), while a ∼ 30% excess is seen in the counts at high fluxes reported in Hasinger et al. (2005). These latter are based on a total of about 200 type-1 AGN selected from the ROSAT Bright Survey (RBS; Schwope et al. 2000). The systematic o ff set with the eRASS1 results may be related to e.g. uncertainties in the sensitivity calculations of the ROSAT sample, or  \n<!-- image -->  \nFig. 11. Left : Cumulative number counts as a function of flux for eRASS1 X-ray point sources selected in the 0.6-2.3 keV band. The apetool sensitivity maps and the eRASS1 aperture photometry for individual sources are used to construct the number counts in di ff erent Galactic latitude intervals, | b | = 0-10 · (brown dashed curves), 10-20 · (yellow solid curves), 20-40 · (purple dashed curves) and 40-90 · (blue solid curves) following methods described in Georgakakis et al. (2008). The shaded region associated with each of these curves corresponds to the 68% uncertainty in the determination of number counts estimated using the bootstrap resampling method described in the text. The curves are plotted up to the flux of 2 × 10 -14 erg s -1 cm -2 where the apetool sensitivity curve corresponds to about 1% of its maximum value. For comparison also shown are the 0.5-2 keV number counts determined by (Hasinger et al. 2005, black solid line) using ROSAT, XMM-Newton and Chandra surveys. The red shaded region and red dashed curve are the 0.5-2 keV number counts estimated from Chandra extra-galactic survey fields (Georgakakis et al. 2008). The extent of the red shaded region at fixed flux corresponds to the 1 σ uncertainty. For fluxes brighter than fX (0.5-2 keV) ≈ 10 -12 erg s -1 cm -2 the red dashed curve is an extrapolation of the best-fit double power determined by Georgakakis et al. (2008). Right : The corresponding sensitivity curve, i.e. area coverage vs. flux limit in the 0.6-2.3 keV band, for di ff erent galactic latitude intervals.  \n<!-- image -->  \nthe inclusion of more extended sources in RBS because of the larger ROSAT PSF.  \nThere is a significant excess for the source population at low Galactic latitude. The absorbed spectral model adopted in the flux calculation is only valid for un-obscured AGN, and the cosmic variance of AGN is of small amplitude and unrelated to Galactic latitude. So the excess must be due to Galactic sources, for which the invalid absorption correction biases high the source fluxes and thus the number counts. Therefore, the low-latitude samples trace the distribution of Galactic X-ray sources (see Freund et al., submitted).', '5.6. The resolved fraction of the Cosmic X-ray Background in eRASS1': 'The cosmic X-ray background (CXB) radiation, discovered by Giacconi and collaborators at the dawn of X-ray astronomy in 1962, can be considered as the ultimate inventory of the energy released by high-energy processes throughout the history of the Universe. As the CXB radiation is dominated by accretion onto black holes, detailed modelling of the CXB over the years has accompanied our deeper understanding of the physical properties of AGN, and of their cosmological evolution. Indeed, deep extra-galactic X-ray surveys have resolved about 80-90% of the CXB using synthesis models of the obscured and un-obscured AGNpopulation (e.g Setti & Woltjer 1989; Comastri et al. 1995; Treister & Urry 2005; Hickox & Markevitch 2006; Gilli et al. 2007; Ueda et al. 2014; Aird et al. 2015; Ananna et al. 2019).  \nHere we provide an estimate of the fraction of the CXB radiation resolved into individual sources by eRASS1 in the soft  \nband (1-2 keV), where the eROSITA sensitivity is around its maximum, and the CXB contribution to the measured background emission dominates over both the foreground Galactic components and over the instrumental background (see e.g. Ponti et al. 2023b). We first adopt a mean value of the CXB intensity in the 1-2 keV band of 12.7 keV cm -2 s -1 sr -1 by fitting the data compilation of Gilli et al. (2007). We then compute the fraction of this intensity in the same energy range contributed by highly reliable eRASS1 sources with DET\\_LIKE\\_0 > 10 (excluding all flagged ones 15 ; see Table 5) above two flux limit thresholds: F 0 . 5 -2 keV > 5 × 10 -14 erg s -1 cm -2 or F 0 . 5 -2 keV > 1 × 10 -14 erg s -1 cm -2 . As discussed in Sect. 5.4 above, the former represents an approximately uniform flux limit achieved over the entire sky, while the latter includes fainter sources detected in the areas of increased exposure near the ecliptic poles.  \nThe two bottom panels of Fig. 12 show a HEALPix map (order 4) of the 1-2 keV CXB intensity resolved fraction in the eROSITA-DE hemisphere. On the right hand side, the fraction computed using the eRASS1 flux limit appears more uniform, with the exception of the low-Galactic latitude regions, where Xray stars in nearby ( < 500 pc) star-forming regions (Sco-Cen association and Gould belt; Schmitt et al. 2022) increase the cumulative average emission. On the left hand side, instead, a significantly higher fraction of the CXB is resolved into fainter sources close to the south ecliptic pole (SEP). Finally, in the top panels of Fig. 12 we show the histogram of the number of HEALPix order 4 pixels at a given resolved fraction in the extra-galactic sky ( | b | > 30 · ) for the two cases. The solid lines are the median  \nvalues: 0.19 and 0.24, respectively. Including only point sources (i.e. removing all those with EXT\\_LIKE > 0) the median value of the resolved fraction is reduced by about 2%, to 0.17 and 0.22, respectively. On the other hand, if we include all (non-flagged) sources with DET\\_LIKE\\_0 > 6 the median value for the CXB resolved fraction remain almost unchanged (0.20) for the uniform flux limit case ( F 0 . 5 -2 keV > 5 × 10 -14 erg s -1 cm -2 ), but increases to 0.33 for F 0 . 5 -2 keV > 1 × 10 -14 erg s -1 cm -2 , due to the large number of fainter sources detected with lower significance in the deeply exposed area near the SEP.', '5.7. Particular regions': 'In the soft X-ray band, a few nearby sources of X-ray emission are very prominent, including the eROSITA bubbles (Predehl et al. 2020), the two Magellanic clouds, the Vela SNR, Virgo cluster, Crab Pulsar, and Sco X-1. They enhance the e ff ective background for the detection of distant sources and thus appear clearly in the flux limit maps displayed in Fig. 9. As discussed in Sect. 5.2.1, our source detection algorithm is not optimized for such regions with bright or extended source contamination. This results in an increased number of spurious sources caused by background fluctuations. There are large uncertainties in the measurement of background and source properties in such regions.  \nIn addition to these, the SEP region also poses particular issues for source detection, because the exposure is deeper than the typical depth across the sky by more than two orders of magnitude. Our detection algorithm is not optimized for such deep exposures either. As discussed by Liu et al. (2022c), one di ffi culty of source detection is de-blending nearby sources (or confusion). This problem is more severe where the fields are crowded, like in the regions near the ecliptic poles. Although the eRASS1 catalogues presented here are complete (i.e. includes also sources in the SEP region), an optimized SEP catalogue will be presented in Liu et al. (in preparation), where a more sophisticated data processing and source detection optimized for crowded regions with deep (and spatially variable) exposure will be introduced.', '6.1. Comparison with XMM source fluxes': "Before eRASS1, the previous best all-sky X-ray catalogue was provided by the ROSAT all-sky survey (2RXS: Boller et al. 2016). In terms of pure source numbers, the largest X-ray catalogue previous to eRASS1 is the XMM-Newton serendipitous catalogue (4XMM: Rosen et al. 2016; Webb et al. 2020). Covering only a small fraction of the sky, 4XMM has a similar resolution to eROSITA and generally deeper exposure than eRASS1. Therefore, we use the 4XMM-DR12 catalogue for a consistency check of our eRASS1 sources.  \nWe selected 1250 . 8k point sources from the eRASS1 1B catalogue and 432 . 9k point sources ( SC\\_SUM\\_FLAG = 0 and SC\\_EXTENT = 0) from 4XMM-DR12 over the entire sky, and performed a simple positional cross-correlation between them, adopting a maximum separation of 10 '' . We found 16 . 5k detection pairs. The majority of the X-ray sources are in the low S / N regime, where their fluxes su ff er large uncertainties and Eddington bias. Taking the ratio of 0.2-2 keV flux to flux uncertainty as a measurement of flux measurement reliability, we selected sources with at least 1 σ reliability in both catalogues and compared their fluxes in the upper panel of Fig. 13. For faint sources, the Eddington bias is still visible in terms of the overestimated  \neROSITA fluxes at the low-flux end. For the brightest sources with 5 σ reliability, the XMM-Newton and eROSITA measured fluxes are consistent, with the mode of the distribution of the flux ratio consistent with unity within 6%. Assuming the XMMNewton fluxes in this band are correct, this can in turn be interpreted as an estimate of the flux calibration uncertainty of eROSITA in the 0.2-2 keV band.  \nWe also compared the hard-band fluxes of the Hard sources with that measured in the 4XMM catalogue, as displayed in the lower panel of Fig. 13. Here we are using slightly di ff erent bands (2.3-5 keV for eROSITA and 2-4.5 keV for XMM-Newton ), but we note here that assuming a power-law spectrum with a slope of Γ = 1 . 7 the fluxes in these two bands are expected to be almost identical, with a di ff erence < 1%. Because of the small sample size (only 131 matches) and the small eROSITA source photon counts in the hard band, this comparison is dominated by selection biases, even though residual e ff ective area calibration issues cannot be ruled out.  \nA small fraction of the sources show significant variability, almost all of which turn brighter in eRASS1. This is because, typically, XMM-Newton observations are much deeper than eRASS1, and thus the sources that become fainter are filtered out by selection bias. The analysis of these variable sources will be presented elsewhere.", '6.2. Astrometric validation': "This section presents a methodology that provides independent constraints on the X-ray positional uncertainty using external multi-wavelength source catalogues with known and accurate positions.  \nWe assume that each source i has a circular, Gaussian positional uncertainty with a standard deviation of σ i . It is then possible to cross-correlate two catalogues and count the number of matches as a function of angular distance. An example of such an angular separation distribution is shown in Fig. 14 and consists of a linear part at large angular distances and a pronounced peak at small separations. The former represents random matches between X-ray sources and the external catalogue, the latter is related to true associations. The number of random matches at a given angular separation depends only on the sky density of the external catalogue and the total number of X-ray sources:  \nN rand( θ ) = NX · 2 π · θ · ρ · d θ, (2)  \nwhere N rand( θ ) is the number of associations for angular separations between θ and θ + d θ and ρ is the sky density of the external catalogue. The peak in Fig. 14 can be expressed as the superposition of NX Rayleigh distributions, each of which depends on the positional uncertainty of the corresponding X-ray source, σ i , so  \nN assoc( θ ) = F NX X i = 1 Pi ( θ | σ i ) , (3)  \nwhere Pi ( θ | σ i ) is the Rayleigh distribution that describes the probability of true associations separated by angular distance θ given the positional uncertainty of the X-ray sources σ i . This parametrisation assumes that the positional errors of the external catalogue are much small than σ i and therefore can be ignored. The factor F , which takes values between zero and one, represents the fraction of X-ray sources expected to have counterparts in the external catalogue. In Fig. 14 the total number of sources at a given angular separation θ can be expressed as the sum of the  \nFig. 12. Resolved CXB fraction in eRASS1. Bottom panels: HEALPix order 4 map of the 1-2 keV CXB intensity resolved fraction (eROSITA-DE hemisphere, Hammer-Aito ff projection, in Galactic coordinates). On the bottom left we include all sources with DET\\_LIKE\\_0 > 10, no flags and X-ray flux F 0 . 5 -2 keV > 1 × 10 -14 erg s -1 cm -2 . On the bottom right we include all sources with DET\\_LIKE\\_0 > 10, no flags and and X-ray flux F 0 . 5 -2 keV > 5 × 10 -14 erg s -1 cm -2 (i.e. the hemisphere median flux limit at 50% completeness). Top panels: distribution of the resolved fraction in the extragalacitc sky ( | b | > 30 · ) in 13.4287 square-degree pixels for the two cases. The solid lines are the median values: 0.24 and 0.19, respectively.  \n<!-- image -->  \nterms given by Equations 2 and 3. Put di ff erently, it is possible to model the observationally determined angular separation distribution of Fig. 14 using the equations above and hence to infer parameters such as F , ρ and the positional uncertainty. The number of matches at a given angular separation bin θ is a Poisson variate with expectation value λ ( θ ) = N rand( θ ) + N assoc( θ ). The likelihood of the observations in Fig. 14 can then be expressed as the product of the Poisson probabilities at each angular separation bin:  \nL = N θ Y j = 1 Pois GLYPH<16> Nj | λ ( θ j ) GLYPH<17> , (4)  \nwhere the index j is for all angular separation bins N θ and Pois( Nj | λ ( θ j )) is the Poisson probability of Nj matches given the expectation value λ ( θ j ). We use the U ltra N est nested sampling algorithm (Buchner 2016, 2019, 2021) to perform Bayesian inference on the likelihood of Equation (4) and determine posteriors for the various model parameters. For this exercise we parameterise the positional uncertainties of individual sources as  \nσ i = POS\\_ERR = q A · σ 2 + σ 2 0 , (5)  \n√  \nwhere σ = RADEC\\_ERR / 2 and RADEC\\_ERR is the catalogued positional uncertainty produced by ermldet (see Eq. 1). In equation 5 above, σ 0 represents systematic uncertainties and A is a multiplicative factor that scales the ermldet uncertainties. Under these assumptions, the four model parameters that are inferred by modelling the distribution of Fig. 14 are ρ , F , A , σ 0. We caution that this parametrisation assumes that the external  \ncatalogue is assumed to have the same sky density, ρ , across the sky. Although the sample selection (see below) minimises such variations, it is inevitable that ρ has an intrinsic scatter that is not accounted for in the current version of the methodology.  \nWe use the catalogue of AGN from Gaia and unWISE Data (Gaia / unWISE; Shu et al. 2019) to cross-match against the eRASS1 X-ray source catalogue. We only consider Gaia / unWISE sources with probability being a QSO PROB\\_RF > 0 . 8 and G -band magnitude < 20 . 5 mag. The latter criterion is adopted to minimise variations in the sky density of QSO candidates because of the variable depth of the GAIA survey as a result of the scanning law of the mission. For this magnitude cut it is empirically found that the sky density of QSO candidates in the extra-galactic sky (Galactic latitudes | b | > 20 · ) is nearly homogeneous. We further limit the eRASS1 catalogue to sources with ermldet detection likelihood DET\\_LIKE\\_0 > 7 (to increase the purity of the sample) that are not spatially extended (parameter EXT = 0) and are not identified as potentially spurious by the algorithms described in Sect. 5.2. We analyse separately eRASS1 sources with low ( < 30) and high ( > 30) number of net counts to accommodate a dependence of the positional uncertainty corrections on source brightness. Figure 14 shows the angular separation distribution from the cross-correlation of the high count eRASS1 sample with the Gaia / unWISE catalogue. The model with parameters set to the median of the corresponding posteriors is shown with the dashed line in that figure. The model represents the observations reasonably well. For the high-counts eRASS1 sub-sample, the parameters that are relevant to the positional error are estimated as A = 1 . 3 ± 0 . 1 and σ 0 = 0 . 9 '' ± 0 . 1 '' .  \nFig. 14. Distribution of the angular separation between eRASS1 X-ray source positions and the Gaia / unWISE QSOs. The red dots correspond to the observed number of associations at a given angular separation bin, δ RX -Opt . We use only eRASS1 X-ray sources with Galactic latitudes in the range | b | = 30 -70 · (i.e. extragalactic sky), DET\\_LKE\\_0 > 7, with more than 30 counts, are not spatially extended (parameter EXT = 0) and are not identified as potentially spurious by the algorithms described in Sect. 5.2. The dashed line is the model described in the text for parameters fixed to the median of the corresponding posteriors.  \n<!-- image -->  \n<!-- image -->  \nFig. 13. Flux comparison between eRASS1 and 4XMM catalogues. The upper panel displays the 0.2-2 keV fluxes of point sources. The brightest sources with at least 5 σ flux measurements are plotted in red points with green 1-sigma uncertainties. The other sources are plotted in blue points without error bars. The black lines indicate 1:1 and deviation from 1:1 by a factor of 2. The lower panel shows a hard-band comparison for the Hard sources.  \n<!-- image -->  \nFor the low-count sub-sample, the 3 σ upper limit of the σ 0 term is about 0.2 '' and consistent with zero within uncertainties. This indicates that the positional error budget for low count sources is dominated by the multiplicative term of Eq. (5), as the typical positional uncertainties of sources at the low count regime are already quite large. For simplicity, we then adopt a single scaling function to determine the positional uncertainty from Eq. 5, which we report as POS\\_ERR in the catalogue  \nPOS\\_ERR = p 1 . 3 · σ 2 + 0 . 9 2 . (6)  \nIndeed, for sources with < 30 counts the di ff erence in positional error between the above more general expression in Eq.(6) and a modified version of this equation where the additive term is assumed to be zero is only ∼ 0.1-0.2 '' on average.  \nAs a further test, we have also tried to calibrate the positional uncertainty of the eRASS1 sources by comparing the catalogue entries with both Chandra and XMM-Newton serendipitous catalogues. The resulting values of the scaling terms A and σ 0 in Eq. 5 are consistent within the uncertainties with those derived above.", '6.3. On the effect of source blending': 'With the relatively large PSF of eROSITA (see Section 2.2.1), source blending could become a issue in crowded regions of the sky. Empirical tests of the e ff ect of blending (or lack thereof) on the catalog would require cross-match against X-ray catalogs where blending is expected to be negligible. Such tests, however, are potentially degenerate with source variability and nonuniformity (of background, exposure, etc.) of the matched catalogues. An alternative method is through simulations, where the full observation and detection process is simulated starting from a realistic sky source population. These have been performed for the eFEDS field (Liu et al. 2022c) and for eRASS1 (Seppi et al. 2022). The biggest impact of source blending is on the detection of extended objects, that is, blending of point sources leads to spurious extended sources in the catalog. Blending also causes incompleteness, since multiple individual sources are considered as only one. At the depth of eFEDS, the incompleteness caused by blending is between 1% and 3%, dependent on source brightness and extent. For eRASS1, at the DET\\_LIKE\\_0 threshold of the main catalog, the fraction of blended sources is less than 1% (see Seppi et al. 2022, Table 3), which is definitely much smaller than the expected level of contamination from spurious detections. Only in the south ecliptic pole region, where the exposure is deeper than average by a factor of more than hundred, this e ff ect becomes more significant. This will be discussed in a separate paper (Liu et al. in prep.).', '7. eROSITA-DE Data Release 1 (DR1)': 'The primary data products of the eROSITA-DE DR1 16 consist of eRASS1 calibrated event files, which contain the information generated by the operating cameras during eRASS1 observations. The data taken during the CalPV phase, which were released as part of the Early Data Release, are not part of DR1.  \nFig. 15. Four examples of half-sky maps released as part of DR1. The maps are colour-coded by the logarithm of the count rate intensity (in cts s -1 arcmin -2 ), displayed in Zenith Equal Area (ZEA) projection in Galactic coordinates, with pixel size of 0.09 deg 2 and centered on ( l , b ) = (270, 0). We show the 0.2-0.6 keV ( top left ), 0.6-2.3 keV ( top right ), 2.3-5 keV ( bottom left ) and 0.2-2.3 ( bottom right ) bands. Please note that each map has a di ff erent (logarithmic) colour scale.  \n<!-- image -->  \nThe eSASS software package to help processing and analysing eROSITA data is also made public, with the name eSASS4DR1.  \nAs discussed above, the eSASS pipeline divides the all-sky observations into 4700 sky tiles for practical purposes of computational tractability (see Fig. 2). eROSITA-DE have proprietary rights on 2248 of these tiles, while 199 of them have shared rights between the German and Russian consortia. eROSITA-DE DR1 finally comprises 2447 sky tiles, of which 199 are only partially filled with eROSITA-DE data. This public eRASS1 data was processed with version 010 of the eSASS pipeline.  \nThe released calibrated event lists of each sky tile contain photons:  \n- -with energies between 0.2-10 keV;\n- -flagged as flag=0xE000F000 , i.e., good events from the nominal field of view, excluding bad pixels;  \n- -with patterns pattern = 15, i.e., including single, double, triple, and quadruple events.  \nBesides calibrated event lists, the DR1 team also releases the following products per sky tile:  \n- -Counts and count rate maps,\n- -Exposure maps,\n- -Background maps,\n- -Tables containing cumulative survey area as a function of limiting flux (for sources with likemin > 5 ),\n- -Sensitivity maps based on aperture photometry. Includes aperture-averaged exposure maps, aperture-integrated background maps, and area curves (survey area sensitive to a given count rate),\n- -Sensitivity maps (for sources with likemin > 5 ),  \n<!-- image -->  \n<!-- image -->  \nFig. 16. Two examples of light curves (top panels) and spectra (bottom panels) for two sources contained in the eRASS1 Main catalogue, generated from the released source products. On the left hand side, we show the light curve and spectrum of the first eROSITA Quasi Periodic Eruptor (eRO-QPE1; Arcodia et al. 2021) (1eRASS J023147.1-102011; DETUID : eb01\\_038099\\_020\\_ML00005\\_002\\_c010). The light curve ( top left ) is binned in individual eroday visits, separated by 4-hours intervals. The spectrum ( bottom left ) is taken from all data, and fitted with an absorbed black body ( zbbody ). On the right hand side, we show one of the brightest AGN in the SEP region (1eRASS J061148.2-662434; DETUID : em01\\_093156\\_020\\_ML00001\\_012\\_c010; Bogensberger et al., in prep.). Given the high ecliptic latitude, the source is visited over two separate long periods of about 14 and 24 days, respectively. The light curve ( top right ) shown has the original binning in 10 seconds intervals. The spectrum ( bottom right ) is fitted with an absorbed power law. Note that both best-fit lines (light blue) are simple physical models aiming to guide the eye and do not intend to encompass all the potential detailed features present in the spectra. The background model is shown as black dotted lines in both figures. For more details on eROSITA spectral modelling and appropriate handling of background spectra, refer to Sect. 3 of Liu et al. (2022b).  \n<!-- image -->', '-Source data products.': 'The data are made available to users using a web-based tool called eRODat. This interface allows one to interactively view the eROSITA sky, using Aladin Lite (Boch & Fernique 2014) to show the eROSITA HiPS (Hierarchical Progressive Surveys) maps and source positions. It also allows the user to identify the sky tiles associated with a given source position or a list of positions, and download the data products for those tiles. In addition, the user can search for sources around a position from the eROSITA X-ray catalogues and view those sources in various surveys, view the details of the catalogue entries, or download the individual source products. Data can be either downloaded by navigating the archive structure to obtain the required data files, following direct links as appropriate, or by adding the products to a virtual basket. The contents of the basket can be download immediately as a single tar file or eRODat can generate a shell script to later download those products.  \nDR1 also includes the publication of the eROSITA upper flux limit server (Tubín-Arenas et al. 2023). The upper limits are calculated using X-ray photometry on the eROSITA standard calibration data products (counts image, background image, and exposure time), following the Bayesian approach described by Kraft et al. (1991). Pre-computed upper limits are  \navailable for every pixel position in the eROSITA-DE sky at a confidence interval of 99.87% (this corresponds to a one-sided 3 σ level). These values are stored using the hierarchical indices (HEALPix) framework to enable a fast search. The products are delivered for the 1B energy band (0.2-2.3 keV), as well as in all sub-bands of the 3B DET run (soft: 0.2-0.6 keV, medium: 0.6-2.3 keV, hard: 2.3-5.0 keV, and the combined total band: 0.2-5.0 keV). The upper flux limit data will be available in two ways: either by downloading the pre-computed data products or through a web tool.', '7.1. Half-sky maps': 'Half-sky counts, count rate and exposure maps are o ff ered (in HiPS format) in seven di ff erent energy bands:  \n- - 0.2-2.3 keV (Main)\n- - 0.2-0.6 keV\n- - 0.6-2.3 keV\n- - 2.3-5.0 keV (Hard)\n- - 0.2-0.5 keV\n- - 0.5-1.0 keV\n- - 1.0-2.0 keV,  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFig. 17. Two more examples of light curves (top panels) and spectra (bottom panels) for two sources contained in the eRASS1 Main catalogue, generated from the released source products. On the left hand side, we show the light curve and spectrum of the bright star AD Leo (1eRASS J101935.6 + 195211; DETUID : em01\\_155069\\_020\\_ML00001\\_011\\_c010). The light curve ( top left ) is shown with the original 10s binning. The spectrum ( bottom left ) is fitted with a multi-temperature thermal collisionally-ionised plasma emission model ( apec+apec+apec , with tied abundances). On the right hand side, we show the light curve and spectrum of SMSS J114447.770-430859.3, the most luminous QSO of the last 9 Gyr (Onken et al. 2022; Kammoun et al. 2023) (1eRASS J114447.6-430858; DETUID : em01\\_176132\\_020\\_ML00001\\_002\\_c010). The light curve ( top right ) is binned in individual eroday visits, while the spectrum ( bottom right ) is fitted with an absorbed power law. Note that both best-fit lines (light blue) are simple physical models aiming to guide the eye and do not intend to encompass all the potential detailed features present in the spectra. The background model is shown as black dotted lines in both figures. For more details on eROSITA spectral modelling and appropriate handling of background spectra, refer to Sect. 3 of Liu et al. (2022b).  \n<!-- image -->  \nwhile the background maps and the sensitivity maps based on aperture photometry are released only for the first four energy bands above. As an example, Fig. 15 shows the half-sky images in Zenith Equal Area projection colour coded by count rate intensity in four bands: 0.2-0.6 keV, 0.6-2.3 keV, 2.3-5 keV and 0.2-2.3 keV, respectively. These images have been created using all TMs (also those a ff ected by light leak), so some artifacts are still present in the lowest-energy map. On the other hand, at the highest energy (2.3-5 keV) the count rate maps are dominated by the (unvignetted) particle background; small temporal fluctuations of such background (including rare Coronal Mass Ejections hitting SRG) are imprinted on those maps as stripes at fixed ecliptic longitudes 17 .  \nZheng et al. 2024 have created custom-made half-sky eRASS1 maps in the standard ROSAT broad bands (R1-R7), with the goal of comparing the eRASS1 and RASS large scale emission in great detail. By only restricting to non-light-leak TMs (1, 2, 3, 4 and 6), and by carefully subtracting the nonvignetted particle background using FWC data, they produce clean maps almost free from artifacts. For more details on top-  \nics such as flux estimation of large-scale di ff use emission for the non-light-leaked TMs, we refer the reader to that work.', '7.2. Source products': "Together with the catalogues of detected sources, eROSITA DR1 also contains a number of source-specific products, which are generated by the eSASS task srctool . In particular, we release source-specific products for 200 217 Main catalogue entries with DET\\_LIKE\\_0 > 20. These products include source and background spectra, their respective ARFs and RMFs, light curves, and source and background event lists. All these source products are o ff ered per telescope module and combined. In this section we describe how these products are generated.  \nsrctool creates spectra, background spectra, response matrices, ancillary response files, and light curves for an input catalogue of sources. The tool chooses a circular source extraction region to optimise the signal to noise ratio of the source spectrum, given the local background surface brightness and the shape of the PSF, clipping the radius to a minimum of 15 '' and maximum of 99% of the PSF encircled energy fraction, assuming a circular PSF, and taking into account any excluded neighbouring contaminating sources. A detailed explanation of the al-  \nthm, including the background definition procedures, can be found online 18 . The final products are produced for each telescope module separately and combined for the telescope module groups TM0 (all seven TMs), TM8 (only those TMs with onchip filter) and TM9 (only those TMs without on-chip filter). The products for each source are placed into separate archive files in order to save disk space usage and reduce the number of files in the archive. At odds with the rest of the pipeline, srctool was configured to use input GTIs calculated from the flaregti task, rather than the standard GTIs. Light curves were created in 10s time bins, discarding those bins when the source is not visible. The light curves were created in energy bands 0.2-0.6, 0.6-2.3 and 2.3-5.0 keV.  \nAs eROSITA scans across the sky during the survey, a source samples di ff erent parts of each telescope and detector over time. The role of srctool is to account for the varying vignetting of the telescope, bad pixels and the e ff ect of photons lost outside the extraction region due to the shape of the PSF. These effects are calculated by stepping through time a set sample points spatially sampling the source. The average corrections are accounted for in the ancillary response matrix for spectral fitting, while the time-variable e ff ects are included in the outputted light curve rates. To avoid double counting these corrections, they are not included in the exposure time of the output files, which instead contain the time eROSITA is looking at the source. Those sources with an extent of zero from the source detection pipeline are assumed to be point sources, while extended sources are assumed to follow a beta model with the core radius given by the extent and an outer power-law slope fixed by the exponent β c = 2 / 3 (see Brunner et al. 2022, Appendix A.5). Table E.1 in the Appendix summarises the changes to the srctool package for DR1.  \nIn Figs. 16 and 17 we present four examples of light curves and spectra of four point sources, generated from the released products.", '8. Summary and outlook': "In this paper we give an overview of the operations, observations, data reduction and analysis of the first eROSITA allsky survey for the hemisphere data (western Galactic) whose proprietary rights lie with the eROSITA-DE Consortium. We present the catalogues of X-ray sources extracted with the standard eSASS pipeline and describe quality control tests we have performed on the catalogued sources (astrometry, photometry, fidelity, etc.). Finally, we summarise the content of the data release associated with the eRASS1 survey so to facilitate interested users from the scientific community at large to interact with the rich data sets generated by SRG / eROSITA.  \nThe scientific and information content of the eRASS1 survey is large, varied and di ffi cult to summarize in a few words. While our focus here is mainly on the eRASS1 X-ray catalogues and on the properties of the X-ray sources that are included there, other works focus on other relevant aspects of the all-sky survey data content. For example, already soon after the completion of eRASS1, Predehl et al. (2020) reported the discovery in the eROSITA all-sky image of large-scale X-ray bubbles in the Milky Way halo (the so-called 'eROSITA bubbles'). Following up on that work, Zheng et al. (2024) have studied maps of the di ff use emission as a function of energy and compared them with the ROSAT ones. Locatelli et al. (2023) have then used  \nTable 8. eRASS1 by (rough) numbers.  \nNotes. We list simplified estimates of the approximate number of calibrated photon counts (in millions) registered by eROSITA in eRASS1 in the 0.2-2.0 keV band in one hemisphere, split into separate physical components: Cosmic X-ray Background (CXB), Milky Way hot Circum-Galactic Medium (MW hot CGM), Instrumental background (FWC), Local Hot Bubble (LHB), Solar Wind Change Exchange (SWCX), Point Sources (PS), Extended Sources (Ext.). For the latter component, we quote in parenthesis the number of net counts from confirmed clusters of galaxies (Bulbul et al., 2024). See text for more details.  \nthose maps to derive structural parameters of the Milky Way Circum-Galactic Medium (CGM) emission. Using dense molecular clouds of known distance as shadowing tools Yeung et al. (2023) have derived constraints on the properties of the Local Hot Bubble from its soft X-ray emission measured by eROSITA.  \nBased on these analysis, here we try to give a simplified answer to the question of what the survey contains by providing a photon (i.e. calibrated events) budget of eRASS1 split into its major physical components. In order to do so, we separate source from background counts on the basis of the catalogues source net counts in di ff erent bands. As for the background and foreground emission, we adopt here the simplifying assumption that the total spectrum in the eFEDS extra-galactic field is representative of the all-sky, and we adopt the best fit model of the di ff use emission spectrum presented in Ponti et al. (2023b) to allocate photons to physical components. The result of this exercise is shown in Table 8. A more detailed inventory, as a function of position in the sky, can be found in Zheng et al. (2024).  \nTime-domain variability analysis of the X-ray sources in eRASS1 will be presented in Boller et al. (in prep.), while Grotova et al. (in prep.) focus on the population of extra-galactic nuclear transients potentially associated to the tidal disruption of stars by SMBHs, which eROSITA is extremely sensitive to (see e.g. Malyali et al. 2021; Sazonov et al. 2021; Malyali et al. 2023; Homan et al. 2023; Liu et al. 2023b).  \nThe details of the cross-identification methodology and of the classification schemes for the various classes of X-ray emitters contained in the eRASS1 catalogues will be described elsewhere. Specifically, the following papers will release the counterpart identification of: clusters of galaxies (Bulbul et al., 2024; Kluge et al., submitted); AGN (Salvato et al., in prep., Waddel et al., submitted); Blazars (Hämmerich et al., in prep.); Cataclysmic Variables (Schwope et al., in prep.); coronal emitting stars (Freund et al., 2024); X-Ray Binaries (Avakyan et al., in prep.); X-ray source population in the Magellanic Clouds system (Maitra et al., in prep., Kaltenbrunner et al., in prep.).  \nWith almost one million entries, most of them never detected before, the eRASS1 half-sky catalogues represent a major advance for our knowledge of the X-ray Universe. Table 9 and Fig. 18 present a comparison of eRASS1 with other catalogues  \nTable 9. Comparison among main catalogs from previous X-ray missions operating, at least partly, in the 'classical' X-ray energy range ( ∼ 0.2-10 keV).  \nNotes. The column N objects lists the approximate number of sources in each catalogue and f Area is the fraction of the sky observed. ∗ The 4XMM-DR12 Hard catalogue (not shown in the figure) is derived from the 4XMM-DR12 by taking all sources for which the 2-5 keV flux is larger than the quoted 2-5 keV flux error.  \nReferences: (1): Forman et al. (1978); (2): Warwick et al. (1981); (3) Wood et al. (1984); (4) heasarc.gsfc.nasa.gov/W3Browse/einstein/ ipc.html , Harris (1990); (5): Boller et al. (2016); (6): heasarc.gsfc.nasa.gov/W3Browse/rosat/wgacat.html ; (7) cxc.cfa.harvard. edu/csc/about2.1.html ; (8) Webb et al. (2020); (9) www.cosmos.esa.int/web/xmm-newton/xmmsl2-ug ; (10) Evans et al. (2020); (11) Brunner et al. (2022).  \nfrom X-ray missions operating in the 'classical' X-ray energy range 0.1-10 keV, and highlight the steady progress in X-ray survey capabilities culminating with the eROSITA catalogues we discuss here. Simply considering the union of all the unique objects catalogued by any previous X-ray mission (without removing possible overlaps), the eRASS1 main catalogue presented here increases the number of known X-ray sources in the published literature by more than 60%.  \nThe sensitive all-sky survey nature of the project implies that data are accumulated for a large variety of astronomical source classes, and for a plethora of possible science applications, well beyond the main mission-design-driving objectives; in other words, eROSITA data are endowed with tremendous legacy value. Indeed, existing all-sky and / or wide-area optical / IR surveys such as Gaia , SDSS, PanSTARRS, DES and Legacy Survey, HSC, VISTA / VHS, WISE already demonstrate the power of combining data sets for a deeper understanding of the high-energy Universe. Such a synergy extends also to longer wavelengths, at a time when the SKA precursors (LOFAR, MWA, ASKAP, MeerKAT) alongside the major observatories (APERTIF, JVLA) are surveying large swaths of the sky at unprecedented depth and speed. Beyond imaging, also large spectroscopic surveys have included among their main programs systematic follow-up of eRASS X-ray sources. Indeed, it is expected that SDSS-V (Kollmeier et al. 2017) and 4MOST (de Jong et al. 2019; Merloni et al. 2019; Finoguenov et al. 2019), by the end of their first survey period, will have accumulated hundreds of thousands of optical spectra of eROSITA X-ray sources.  \nLooking forward, future eROSITA data releases will comprise data from multiple all-sky surveys. Besides providing a deeper view of the X-ray sky, this will open up the systematic study of variability on months / year timescales. Extrapolating from the current data quality, and with the foreseen improvements of the calibration in sight, the next data releases will include high-fidelity catalogues with several million X-ray sources, closely following the expectations laid down in the early phases of the project (Merloni et al. 2012).  \nAcknowledgements. We thank the referee, Prof. F. Bauer, for the careful review and for the many constructive suggestions that significantly improved the final version of the paper. This work is based on data from eROSITA, the soft X-ray instrument aboard SRG, a joint Russian-German science mission supported by the Russian Space Agency (Roskosmos), in the interests of the Russian Academy of Sciences represented by its Space Research Institute (IKI), and the Deutsches Zentrum für Luft- und Raumfahrt (DLR). The SRG spacecraft was built by Lavochkin Association (NPOL) and its subcontractors, and is operated by NPOL with support from the Max-Planck Institute for Extraterrestrial Physics (MPE). The development and construction of the eROSITA X-ray instrument was led by MPE, with contributions from the Dr. Karl Remeis Observatory Bamberg, the University of Hamburg Observatory, the Leibniz Institute for Astrophysics Potsdam (AIP), and the Institute for Astronomy and Astrophysics of the University of Tübingen, with the support of DLR and the Max Planck Society. The Argelander Institute for Astronomy of the University of Bonn and the Ludwig Maximilians Universität Munich also participated in the science preparation for eROSITA. The eROSITA data shown here were processed using the eSASS / NRTA software system developed by the German eROSITA consortium. Part of this work was supported by the German Deutsche Forschungsgemeinschaft, DFG project number Ts 17 / 2-1. E. Bulbul, A. Liu, V. Ghirardini, C. Garrel and X. Zhang acknowledge financial support from the European Research Council (ERC) Consolidator Grant under the European Union's Horizon 2020 research and innovation programme (grant agreement CoG DarkQuest No 101002585). M. Brusa, A. Georgakakis and B. Musiimenta acknowledge funding from the European Union's Horizon 2020 research and innovation program under the Marie SkłodowskaCurie grant agreement No 860744. Z. Igo, C. Aydar, J. Wolf acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy - EXC-2094 - 390783311. G. Ponti acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 865637). P. Li is supported by the Alexander von Humboldt Foundation. D. Tubin-Arenas acknowledges support by DLR grant FKZ 50 OR 2203. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France.", 'References': 'Adami, C., Giles, P., Koulouridis, E., et al. 2018, A&A, 620, A5 Aird, J., Coil, A. L., Georgakakis, A., et al. 2015, MNRAS, 451, 1892 Alexander, D. M. & Hickox, R. C. 2012, New A Rev., 56, 93 Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARA&A, 49, 409 Ananna, T. T., Treister, E., Urry, C. M., et al. 2019, ApJ, 871, 240 Arcangeli, L., Borghi, G., Bräuninger, H., et al. 2017, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 10565, Proc. SPIE, 1056558  \nFig. 18. Graphical comparison of selected X-ray catalogs, based on the data presented in Table 9. For each catalogue we plot the total number of objects vs. the time of the start of the corresponding data-taking period. Each catalogue is represented by a circle, whose radius is proportional to the logarithm of the \'discovery rate" (number of objects divided by the number of months of operations), while the shaded portion of the circle represents the fraction of the sky covered by the catalogue. The colour coding distinguishes telescopes operating mainly in the soft ( < 2 . 5 keV; red), hard ( > 2 keV; green) or broad (both soft and hard; blue) X-ray band. The datapoint corresponding to a putative 4XMM Hard catalog would sit very close to the full-band 4XMM one.  \n<!-- image -->  \nArcodia, R., Merloni, A., Nandra, K., et al. 2021, Nature, 592, 704  \nArnaud, M. 2005, in Background Microwave Radiation and Intracluster Cosmology, ed. F. Melchiorri & Y. Rephaeli (Enrico Fermi International School of Physics Course CLIX), astro-ph / 0508159  \nBahcall, N. A. 1977, ARA&A, 15, 505  \nBallet, J. 1999, A&AS, 135, 371  \nBecker, W. & Truemper, J. 1997, A&A, 326, 682  \nBlasi, P. 2013, A&A Rev., 21, 70  \nBoch, T. & Fernique, P. 2014, in Astronomical Society of the Pacific Conference Series, Vol. 485, Astronomical Data Analysis Software and Systems XXIII, ed. N. Manset & P. Forshay, 277  \nBoller, T., Freyberg, M. J., Trümper, J., et al. 2016, A&A, 588, A103  \nBorgani, S. 2008, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 740, APan-Chromatic View of Clusters of Galaxies and the Large-Scale Structure, ed. M. Plionis, O. López-Cruz, & D. Hughes, 287  \nBorgani, S. & Kravtsov, A. 2011, Advanced Science Letters, 4, 204  \nBrandt, W. N. & Alexander, D. M. 2015, A&A Rev., 23, 1  \nBrandt, W. N. & Hasinger, G. 2005, ARA&A, 43, 827  \nBrunner, H., Liu, T., Lamer, G., et al. 2022, A&A, 661, A1  \nBuchner, J. 2016, Statistics and Computing, 26, 383  \nBuchner, J. 2019, PASP, 131, 108005  \nBuchner, J. 2021, The Journal of Open Source Software, 6, 3001  \nBulbul, E., Liu, A., Pasini, T., et al. 2022, A&A, 661, A10  \nBulbul, E. & Liu, A. e. a. 2024, A&A, in press  \nCamilloni, F., Becker, W., Predehl, P., et al. 2023, A&A, 673, A45  \nCampana, S., Combes, F., Elbaz, D., et al. 2021, A&A, 647, E1  \nCavaliere, A. & Fusco-Femiano, R. 1978, A&A, 70, 677  \nCen, R. & Ostriker, J. P. 1999, ApJ, 514, 1  \n- Churazov, E., Khabibullin, I., Bykov, A. M., Lyskova, N., & Sunyaev, R. 2023, A&A, 670, A156  \nClerc, N., Sadibekova, T., Pierre, M., et al. 2012, MNRAS, 423, 3561 Comastri, A., Setti, G., Zamorani, G., & Hasinger, G. 1995, A&A, 296, 1 Comparat, J., Eckert, D., Finoguenov, A., et al. 2020, The Open Journal of As-  \n- trophysics, 3, 13  \nComparat, J., Merloni, A., Salvato, M., et al. 2019, MNRAS, 487, 2005  \nCoutinho, D., Ramos-Ceja, M. E., Dennerl, K., et al. 2022, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 12181, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. J.-W. A. den Herder, S. Nikzad, & K. Nakazawa, 121811A  \nCutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources.  \n- Cutri, R. M., Wright, E. L., Conrow, T., et al. 2021, VizieR Online Data Catalog, II / 328  \nDauser, T., Falkner, S., Lorenz, M., et al. 2019, A&A, 630, A66  \nDavé, R., Cen, R., Ostriker, J. P., et al. 2001, ApJ, 552, 473  \nDavis, J. E. 2001, ApJ, 562, 575  \nde Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, The Messenger, 175, 3  \nDennerl, K., Andritschke, R., Bräuninger, H., et al. 2020, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 11444, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 114444Q  \nElvis, M. 2020, Nature Astronomy, 4, 23  \nEvans, P. A., Page, K. L., Osborne, J. P., et al. 2020, ApJS, 247, 54  \nFabian, A. C. 2012, ARA&A, 50, 455  \nFender, R. P., Belloni, T. M., & Gallo, E. 2004, MNRAS, 355, 1105  \nFinoguenov, A., Merloni, A., Comparat, J., et al. 2019, The Messenger, 175, 39  \n- Forman, W., Jones, C., Cominsky, L., et al. 1978, ApJS, 38, 357  \nFoster, G., Poppenhaeger, K., Ilic, N., & Schwope, A. 2022, A&A, 661, A23  \nFreund, S. & Czesla, S. e. a. 2024, A&A, in press  \nFriedrich, P., Bräuninger, H., Budau, B., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8443, Space Telescopes and Instrumentation 2012: Ultraviolet to Gamma Ray, ed. T. Takahashi, S. S. Murray, & J.-W. A. den Herder, 84431S  \nFriedrich, P., Bräuninger, H., Budau, B., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7011, Proc. SPIE, 70112T  \nFürmetz, M., Pfe ff ermann, E., Predehl, P., Rohé, C., & Tiedemann, L. 2008, in Space Telescopes and Instrumentation 2008: Ultraviolet to Gamma Ray, Vol. 7011, International Society for Optics and Photonics, 70113Y  \nGaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1 Georgakakis, A., Nandra, K., Laird, E. S., Aird, J., & Trichas, M. 2008, MNRAS, 388, 1205  \nGhirardini, V., Bulbul, E., Hoang, D. N., et al. 2021, A&A, 647, A4 Giacconi, R., Branduardi, G., Briel, U., et al. 1979, ApJ, 230, 540 Giacconi, R., Gursky, H., Paolini, F. R., & Rossi, B. B. 1962, Phys. Rev. Lett., 9,  \n- 439\n- Giacconi, R., Kellogg, E., Gorenstein, P., Gursky, H., & Tananbaum, H. 1971, ApJ, 165, L27\n- Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79\n- Gokus, A., Rau, A., Wilms, J., et al. 2020, The Astronomer\'s Telegram, 13657, 1\n- Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759\n- Haberl, F. 2007, Ap&SS, 308, 181\n- Haberl, F., Maitra, C., Carpano, S., et al. 2022, A&A, 661, A25\n- Haberl, F., Wilms, J., Gokus, A., et al. 2020, The Astronomer\'s Telegram, 13828, 1\n- Haiman, Z., Allen, S., Bahcall, N., et al. 2005, arXiv e-prints, astro  \nHarris, D. E. 1990, The Einstein Observatory Catalog of IPC X-ray Sources  \nHasinger, G., Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417  \nHickox, R. C., Jones, C., Forman, W. R., et al. 2009, ApJ, 696, 891  \nHickox, R. C. & Markevitch, M. 2006, ApJ, 645, 95  \nHøg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27  \nHoman, D., Krumpe, M., Markowitz, A., et al. 2023, A&A, 672, A167  \nHopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, ApJS, 175, 356  \nKammoun, E. S., Igo, Z., Miller, J. M., et al. 2023, MNRAS  \n- Klein, M., Oguri, M., Mohr, J. J., et al. 2022, A&A, 661, A4\n- Koenig, O., Wilms, J., Kreykenbohm, I., et al. 2020, The Astronomer\'s Telegram, 13765, 1\n- Kollmeier, J. A., Zasowski, G., Rix, H.-W., et al. 2017, arXiv e-prints, arXiv:1711.03234\n- König, O., Wilms, J., Arcodia, R., et al. 2022, Nature, 605, 248 Kormendy, J. & Ho, L. C. 2013, ARA&A, 51, 511  \nKraft, R. P., Burrows, D. N., & Nousek, J. A. 1991, ApJ, 374, 344  \nLammer, H., Selsis, F., Ribas, I., et al. 2003, ApJ, 598, L121  \n- Liu, A., Bulbul, E., Ghirardini, V., et al. 2022a, A&A, 661, A2\n- Liu, A., Bulbul, E., Ramos-Ceja, M. E., et al. 2023a, A&A, 670, A96\n- Liu, T., Buchner, J., Nandra, K., et al. 2022b, A&A, 661, A5\n- Liu, T., Merloni, A., Comparat, J., et al. 2022c, A&A, 661, A27\n- Liu, Z., Malyali, A., Krumpe, M., et al. 2023b, A&A, 669, A75  \nLocatelli, N. & Ponti, G. e. a. 2023, A&A, in press [  \narXiv:2310.10715  \n]  \n- Madsen, K. K., Burwitz, V., Forster, K., et al. 2021, arXiv e-prints, arXiv:2111.01613  \nMaitra, C., Haberl, F., Sasaki, M., et al. 2022, A&A, 661, A30  \nMalyali, A., Liu, Z., Rau, A., et al. 2023, MNRAS, 520, 3549  \nMalyali, A., Rau, A., Merloni, A., et al. 2021, A&A, 647, A9  \n- Marshall, H. L., Chen, Y., Drake, J. J., et al. 2021, AJ, 162, 254  \nMayer, M. G. F., Becker, W., Predehl, P., & Sasaki, M. 2023, A&A, 676, A68  \nMehrtens, N., Romer, A. K., Hilton, M., et al. 2012, MNRAS, 423, 1024  \n- Meidinger, N., Andritschke, R., Bornemann, W., et al. 2014, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 9144, Proc. SPIE, 91441W\n- Meidinger, N., Andritschke, R., Dennerl, K., et al. 2021, Journal of Astronomical Telescopes, Instruments, and Systems, 7, 025004  \nMerloni, A., Alexander, D. A., Banerji, M., et al. 2019, The Messenger, 175, 42 Merloni, A., Predehl, P., Becker, W., et al. 2012, arXiv e-prints, arXiv:1209.3114 Nevalainen, J. & Molendi, S. 2023, A&A, 676, A142  \n- Nicastro, F., Kaastra, J., Krongold, Y., et al. 2018, Nature, 558, 406\n- Norman, M. L. 2005, in Background Microwave Radiation and Intracluster Cosmology, ed. F. Melchiorri & Y. Rephaeli (Enrico Fermi International School of Physics Course CLIX), astro-ph / 0511451  \nOnken, C. A., Lai, S., Wolf, C., et al. 2022, PASA, 39, e037  \nÖzel, F. & Freire, P. 2016, ARA&A, 54, 401  \nPasini, T., Brüggen, M., Hoang, D. N., et al. 2022, A&A, 661, A13  \nPeebles, P. J. E. 1980, The large-scale structure of the universe  \n- Pi ff aretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & Melin, J. B. 2011, A&A, 534, A109  \nPillepich, A., Porciani, C., & Reiprich, T. H. 2012, MNRAS, 422, 44 Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., & Ventura, P. 2003, A&A, 397, 147  \nPlucinsky, P. P., Beardmore, A. P., Foster, A., et al. 2017, A&A, 597, A35 Ponti, G., Sanders, J. S., Locatelli, N., et al. 2023a, A&A, 670, A99 Ponti, G., Zheng, X., Locatelli, N., et al. 2023b, A&A, 674, A195 Poppenhaeger, K., Ketzer, L., & Mallonn, M. 2021, MNRAS, 500, 4560 Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&A, 647, A1 Predehl, P., Sunyaev, R. A., Becker, W., et al. 2020, Nature, 588, 227 Ramos-Ceja, M. E., Oguri, M., Miyazaki, S., et al. 2022, A&A, 661, A14 Reiprich, T. H., Basu, K., Ettori, S., et al. 2013, Space Science Reviews, 177, 195  \nRemillard, R. A. & McClintock, J. E. 2006, ARA&A, 44, 49 Rosati, P., Borgani, S., & Norman, C. 2002, ARA&A, 40, 539 Rosen, S. R., Webb, N. A., Watson, M. G., et al. 2016, A&A, 590, A1  \n- Rothschild, R., Boldt, E., Holt, S., et al. 1979, Space Science Instrumentation, 4, 269\n- Rutledge, R. E. 1998, PASP, 110, 754, \\_eprint: astro-ph / 9802256\n- Salvato, M., Buchner, J., Budavári, T., et al. 2018, MNRAS, 473, 4937\n- Salvato, M., Wolf, J., Dwelly, T., et al. 2022, A&A, 661, A3\n- Sanders, J. S., Bi ffi , V., Brüggen, M., et al. 2022, A&A, 661, A36\n- Sanz-Forcada, J., Micela, G., Ribas, I., et al. 2011, A&A, 532, A6 Sarazin, C. L. 1986, Reviews of Modern Physics, 58, 1\n- Sazonov, S., Gilfanov, M., Medvedev, P., et al. 2021, MNRAS, 508, 3820 Scargle, J. D., Norris, J. P., Jackson, B., & Chiang, J. 2013, ApJ, 764, 167 Schmitt, J. H. M. M. 1997, A&A, 318, 215\n- Schmitt, J. H. M. M., Czesla, S., Freund, S., Robrade, J., & Schneider, P. C. 2022, A&A, 661, A40\n- Schneider, P. C., Freund, S., Czesla, S., et al. 2022, A&A, 661, A6\n- Schwope, A., Hasinger, G., Lehmann, I., et al. 2000, Astronomische Nachrichten, 321, 1\n- Secrest, N. J., von Hausegger, S., Rameez, M., et al. 2021, ApJ, 908, L51\n- Seppi, R., Comparat, J., Bulbul, E., et al. 2022, A&A, 665, A78\n- Seppi, R., Comparat, J., Nandra, K., et al. 2023, A&A, 671, A57\n- Setti, G. & Woltjer, L. 1989, A&A, 224, L21\n- Shapiro, S. L. & Teukolsky, S. A. 1983, Black holes, white dwarfs and neutron stars. The physics of compact objects\n- Shu, Y., Koposov, S. E., Evans, N. W., et al. 2019, MNRAS, 489, 4741 Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18\n- Sunyaev, R., Arefiev, V., Babyshkin, V., et al. 2021, A&A, 656, A132\n- Tamba, T., Odaka, H., Bamba, A., et al. 2022, PASJ, 74, 364\n- Tanimura, H., Aghanim, N., Kolodzig, A., Douspis, M., & Malavasi, N. 2020, A&A, 643, L2\n- Treiber, H., Vasilopoulos, G., Bailyn, C. D., et al. 2021, MNRAS, 503, 6187 Treister, E. & Urry, C. M. 2005, ApJ, 630, 115\n- Truemper, J. 1982, Advances in Space Research, 2, 241\n- Tubín-Arenas, D. & Krumpe, M. e. a. 2023, A&A, in press\n- Ueda, Y., Akiyama, M., Hasinger, G., Miyaji, T., & Watson, M. G. 2014, ApJ, 786, 104\n- Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 692, 1060 Vink, J. 2012, A&A Rev., 20, 49\n- Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389\n- Voges, W., Aschenbach, B., Boller, T., et al. 2000, IAU Circ., 7432, 3\n- Voit, G. M. 2005, Reviews of Modern Physics, 77, 207\n- Warner, B. 1995, Cataclysmic variable stars, Vol. 28\n- Warwick, R. S., Marshall, N., Fraser, G. W., et al. 1981, MNRAS, 197, 865\n- Webb, N. A., Coriat, M., Traulsen, I., et al. 2020, A&A, 641, A136\n- Weinberg, D. H., Mortonson, M. J., Eisenstein, D. J., et al. 2013, Phys. Rep., 530, 87\n- Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9\n- Whelan, B., Veronica, A., Pacaud, F., et al. 2022, A&A, 663, A171\n- Willingale, R., Starling, R. L. C., Beardmore, A. P., Tanvir, N. R., & O\'Brien, P. T. 2013, MNRAS, 431, 394\n- Wilms, J., Kreykenbohm, I., Weber, P., et al. 2020, The Astronomer\'s Telegram, 13416, 1\n- Wilms, J., Kreykenbohm, I., Weber, P., et al. 2020, The Astronomer\'s Telegram, 13416, 1\n- Wood, K. S., Meekins, J. F., Yentis, D. J., et al. 1984, ApJS, 56, 507\n- Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48  \nYeung, M. C. H., Freyberg, M. J., Ponti, G., et al. 2023, A&A, 676, A3  \nZheng, X. & Ponti, G. e. a. 2024, A&A, in press  \n- Zonca, A., Singer, L., Lenz, D., et al. 2019, Journal of Open Source Software, 4, 1298\n- 1 Max-Planck-Institut für Extraterrestrische Physik, Gießenbachstraße, D-85748 Garching, Germany\n- 2 Leibniz Institut für Astrophysik Potsdam, An der Sternwarte 16, D14482 Potsdam, Germany\n- 3 Institut für Astronomie und Astrophysik, Universität Tübingen, Sand 1, D-72076 Tübingen, Germany\n- 4 Institute for Astronomy and Astrophysics, National Observatory of Athens, V. Paulou and I. Metaxa, 11532, Greece\n- 5 Exzellenzcluster ORIGINS, Boltzmannstr. 2, 85748, Garching, Germany\n- 6 Dr. Karl Remeis-Sternwarte and Erlangen Centre for Astroparticle Physics, Friedrich-Alexander Universität Erlangen-Nürnberg, Sternwartstraße 7, 96049 Bamberg, Germany\n- 7 Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany\n- 8 MIT Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA  \n- 9 Argelander Institute for Astronomy, University of Bonn, Auf dem Hügel 71, 53121 Bonn, Germany\n- 10 Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69,53121 Bonn, Germany\n- 11 Dipartimento di Fisica e Astronomia "Augusto Righi", Università di Bologna, via Gobetti 93 / 2, 40129 Bologna, Italy\n- 12 INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, via Gobetti 93 / 3, 40129 Bologna, Italy\n- 13 IRAP, Université de Toulouse, CNRS, UPS, CNES, Toulouse, France\n- 14 Tautenburg Landessternwarte, Sternwarte 5, 07778 Tautenburg, Germany\n- 15 University Observatory Munich, Faculty of Physics, LudwigMaximilians-Universität, Scheinerstr. 1, 81679 Munich, Germany\n- 16 Department of Physics & McDonnell Center for the Space Sciences, Washington University in St. Louis, One Brookings Drive, St. Louis, MO63130, USA\n- 17 Universität Innsbruck, Institut für Astro- und Teilchenphysik, Technikerstraße 25 / 8, 6020 Innsbruck, Austria\n- 18 TU Dresden, Institute of Nuclear and Particle Physics, 01062 Dresden, Germany\n- 19 DESY, Notkestraße 85, 22607 Hamburg, Germany\n- 20 Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA\n- 21 Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands\n- 22 Istituto di Radioastronomia IRA-INAF, via Gobetti 101, 40129 Bologna, Italy\n- 23 INAF, Osservatorio Astronomico di Brera, Via E. Bianchi 46, 23807 Merate, (LC), Italy\n- 24 Max-Planck Institut für Astronomie, Königstuhl 17, 69117 Heidelberg\n- 25 CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei 230026, China\n- 26 School of Astronomy and Space Science, University of Science and Technology of China, Hefei 230026, China', 'Appendix A: The eROSITA Survey PSF': "The PSF of the eROSITA telescopes can be obtained from the survey itself to verify the PANTER measurements and compare to its shapelet representation used in the source detection pipeline. We are undertaking a project to obtain the PSF from source stacking and present here preliminary results. 938 sources were taken from the full eROSITA sky, selected from a catalogue of eROSITA point sources with the largest number of counts. To optimise the signal and to avoid bad pixels we used eRASS:4 for TMs 1, 2, 3 and 6, eRASS:3 for TMs 5 and 7, and eRASS:2 for TM4. In this analysis we examined the PSF range out to 4 arcmin radius, to be consistent with the ground based measurements, although there is additional flux beyond this radius, estimated to be about 2.7% from the best fit core PSF (excluding stray light) of Churazov et al. (2023) (their equation A.1). We excluded the very brightest piled-up sources and those in crowded fields or near large extended sources. The position of each source was measured by simultaneously fitting 2D models to two X-ray images, one where the source was inside 15 arcmin of the optical axis, and another when it was outside. We matched against the Gaia G band catalogue to identify those sources potentially a ff ected by optical loading and used the Galactic latitude to flag those potentially a ff ected by scattering halos. We manually inspected the images and profiles of the sources to look for anomalous profiles, such as evidence of extension, or identify neighbouring contaminating sources. For sources which could be contaminated in the outskirts by other sources or by non-flat background, we excluded the source, or manually reduced the outer radius used from 4 arcmin to 1, 2 or 3 arcmin as appropriate. Lightcurves created with srctool were used to quantify pileup, where the inner stacking radius used was chosen as a function of time for a source according to pileup rate thresholds (rates of 15, 44, 100, 155, 250 and 600 s -1 produced inner radii of 20, 30, 40, 60, 90 and 240 arcsec, respectively, based on simulations to estimate the e ff ect of pileup). Source profiles were stacked in each energy band, and for each combination of inner and outer radius. A PSF model plus flat background model was simultaneously fitted to the stacked profiles in all the combinations of radial range, allowing the PSF normalisation and background level to vary in each combination of radial range, to produce a model of the PSF as a function of radius. This process was repeated for each energy band and combinations of TMs.  \nFigure A.1 shows a comparison between the result of our source stacking analysis, ground-based PANTER measurements (Dennerl et al. 2020) and survey-averaged shapelet models of the ground-based data (Brunner et al. 2022). Also shown are comparisons of the PSF for the di ff erent eROSITA TMs. We note here that the survey-averaged shapelet PSF is never fitted to the data by ermldet in photon imaging mode, but rather, a PSF appropriate for each photon is created. Within 1 arcmin radius the deviations of the PSF are within 20%. Much of the deviation in the shapelet PSF is due to 'steps' occurring at intervals of 16 arcsec radius. At larger radius the shapelet PSF in the soft band is lower than the stacked or ground-based PSF. This contributes to the shapelet cumulative fraction being higher in the centre, as the profiles are normalised at 4 arcmin radius. The stacked profiles typically have smaller full-width-halfmaximums (FWHMs) than the PANTER and shapelet PSFs, but comparable HEW. The in-flight eROSITA camera sub-pixel positioning allows a better determination of the PSF core than the 9.6 arcsec pixel would otherwise suggest.  \nFinally, we tested the uniformity of the PSF across the sky. For this, we divided the sample of 938 point sources described  \nabove into four quadrants in ecliptic coordinates (with dividing lines at ecliptic longitude of 180 · and ecliptic latitude of 0 · ) and recompute the HEW of the stacked PSF for TM0 in the 0.22.3 keV band. The values are consistent with the all-sky average HEW(30.0 '' ) within 0.3 '' .  \nTable A.1. Point Spread Function Half-Energy Width (HEW) in the 0.22.3 keV band, computed by stacking point sources in eRASS1 in four di ff erent quadrant across the sky.  \nIn the future we plan to improve our stacking procedure. These improvements could include automated detection of contaminating sources and removal of sources with inconsistent profiles. Future processing versions will improve the accuracy of the boresight calibration as a function of time. By choosing which sources to include in the stacking as a function of radius to optimise the signal to noise, we would also reduce the e ff ect of the sky background. Finally, we will also verify the obtained PSFs by studying di ff erent subsamples of sources.", 'Appendix A.1: On the azimutal symmetry of PSF and positional uncertainties': "A key assumption for the astrometric calibration and positional uncertainty determination of our source detection algorithm (see Section 5.1), and of the srctool extraction tool is that the eROSITA PSF can be approximated with a azimutally symmetric function and that the associated positional errors are also symmetric.  \nHere we validate this assumption in two ways. First of all, we show in Fig. A.2 the profile of the average stacked PSF (using the same 938 sources from the all-sky survey as described in the previous section) along ecliptic longitude and latitude. No significant di ff erence in the two profiles can be seen. Secondly, we show in Fig. A.3 the measured positional o ff set between eRASS1 X-ray source positions and Gaia / unWISE QSOs (see Section 6.2) in both equatorial and ecliptic coordinates directions. We fit these distributions with a Gaussian function, and obtain for the o ff set the following values of mean and Standard Deviation for the two cases: ( ∆ α, ∆ δ )mean = ( -0 . 11 '' , 0 . 08 '' ); ( ∆ α, ∆ δ )SD = (3 . 5 '' , 3 . 5 '' ); ( ∆ λ, ∆ β )mean = ( -0 . 11 '' , 0 . 16 '' ) and ( ∆ λ, ∆ β )SD = (3 . 5 '' , 3 . 5 '' ). From this, we conclude that the symmetry approximation is justified, and that the residual o ff set (of size < 0 . 2 '' ) is small enough when compared to the positional uncertainty to be safely ignored.  \nFig. A.1. Comparison of di ff erent survey-averaged eROSITA PSFs. Shown are the PSFs in the bands 0.2-2.3 keV (left panels), 2.3-5.0 keV (central panels), or for the individual TMs from stacking in 0.2-2.3 keV (right panels). The PSF models shown are the survey-averaged shapelet PSFs, those obtained by stacking sources (see text in Appendix A for more details) and those from ground-based measurements using PANTER. The top row of panels show the PSF surface brightness profiles, normalised within 4 arcmin radius, where the vertical lines plot the FWHM values. The second panels down show the fractional di ff erence of the surface brightness of each PSF from the average stacked profile. The third panels down show the cumulative signal as a function of radius, plotting the HEWs as vertical lines. The lowest panels show the fractional di ff erence of the cumulative profiles to the average stacked profile. The shapelet PSFs are those used for fitting the energy band given, the stacked PSF is weighted by the spectra of the input sources, and the PANTER PSF is obtained at the monochromatic energy specified, which is chosen to be representative of the source photons in the band. The shapelet PSF is plotted as a dashed line outside a radius of 1 arcmin, the maximum used for fitting in the source detection pipeline. In the two leftmost columns, the PSF images were rebinned with 4 arcsec pixels where necessary for a fair comparison. In the rightmost column, we show the best fitting model and uncertainties, rather than stacked, rescaled and rebinned count profiles. The FWHMs quoted are sensitive to the pixelisation used and the inner value, and so should be used with care. HEWs were computed from minimally-binned images.  \n<!-- image -->  \nFig. A.2. Stacked PSF profile from the eRASS1 point sources in the 0.22.3 keV band along ecliptic longitude (red) and latitude (blue), demonstrating the symmetry of the PSF. The measured Half-Energy Widths (HEW) are consistent to within 0.6 '' .  \n<!-- image -->  \n<!-- image -->  \nFig. A.3. Distribution of measured positional o ff set between astrometrically corrected eRASS1 sources with DET\\_LIKE\\_0 > 7 and Gaia / unWISE QSOs (see Section 6.2) along equatorial (left) and ecliptic (right) coordinates in arcseconds. The dashed lines represent the best fit gaussian function to these distributions, the parameters of which are reported in the text.  \n<!-- image -->", 'Appendix B: Testing the energy calibration': "An excellent target for testing the energy calibration is the oxygen-rich SNR 1E 0102-7219 (in the following abbreviated as 1E0102), the brightest SNR in the SMC. It is characterized by strong emission lines of O, Ne, and Mg, exhibits only little 'contaminating' emission from Fe, it is su ffi ciently compact to utilise the high spectral resolution provided by slitless Xray gratings ( XMM-Newton RGS and Chandra HETG), yet extended enough to minimise any problems with pile-up. This object has been adopted as a standard calibration source by the International Astronomical Consortium for High-Energy Calibration (IACHEC), which has developed a standard (purely empirical) model specifically designed for calibration (Plucinsky et al. 2017).  \nIn Fig. B.1 we show spectra of 1E0102 taken in dedicated calibration observations on 2019 Nov 7 - 8 (60 - 61 ks, 16 ks for TM6), and on 2021 Nov 26 - 27 (47 - 49 ks, 21 ks for TM4), thus covering a time span of more than 2 years. The data were taken with the same onboard processing mode which was used in eRASS. A major di ff erence, however, were the CCD temperatures, which ranged during the first observation between -85.5 and -84.6 C, and during the last observation between -77.9 and -77.0 C. When comparing both observations it should also be considered that, although 1E0102 was observed on-axis, its precise location on the CCDs was di ff erent between both observations.  \nSpectra were extracted with the eSASS tasks evtool and srctool for each observation and each TM (except for TM5 and TM7, which are a ff ected by the light leak), separately for a source and a background region. A circular source extraction region was used, with a radius of 1 arcmin centered on the FK5 coordinates α = 16 . 006258 deg , δ = -72 . 032394 deg, and the background was extracted from a circle with a radius of 5 arcmin centered on the FK5 coordinates α = 16 . 169575 deg , δ = -71 . 828828 deg.  \nThe fits are based on the standard IACHEC model for 1E0102, which consists of 52 narrow Gaussian emission lines, superimposed on an absorbed continuum. The emission lines are organised into 4 groups, corresponding to emission from O VII, OVIII, Ne IX, and Ne X 19 .  \nFor an assessment of the quality of the energy reconstruction we performed a combined fit of TMs 12346 from both observations (10 spectra in total), using only single pixel events. We applied the standard IACHEC model and treated the normalisations of the O VII, O VIII, Ne IX, and Ne X line complexes as free, but TMindependent, parameters (4 free parameters). Only the overall normalisation was adjusted individually for each TM (10 free parameters). We allowed for TM specific shifts of the energy scales by XSPEC 'gain fits', with all slopes fixed to 1.0 and the 10 individual o ff sets as additional free parameters. This resulted in a common fit of 10 spectra with 24 free parameters. The fit yields χ 2 = 4055 . 2 for 2290 degrees of freedom, or χ 2 r = 1 . 77 (Fig. B.1), and the mean energy shift is -1 . 3 eV for the first and + 2 . 2 eV for the last observation, with a scatter of ± 1 . 0 eV and ± 3 . 0 eV.  \nThese long pointed observations of a line-rich SNR represent a benchmark test of the energy calibration. Considering that the calibration requirements for eRASS spectra are more relaxed due to the much shorter exposure times, we conclude that the energy calibration is su ffi ciently accurate for the sources detected in the eROSITA all-sky survey.  \nFig. B.1. eROSITA spectra of the SNR 1E 0102-7219, taken with TM12346 in Nov 2019 and 2021.  \n<!-- image -->", 'Appendix C: The NRTA pipeline': "As outlined in Sect. 3.1, the field of view of eROSITA, which has a diameter of roughly one degree, scans the sky continuously with a rotational period of four hours and visits each position in the sky typically for six consecutive scans. Each position is then revisited roughly six months later. This cadence provides a unique opportunity to identify and study transient phenomena in the X-ray sky and o ff ers a compromise between time resolution and sensitivity. The survey schedule of eROSITA does not allow for interruptions in order to perform pointed observations, therefore fast identification of transient events and the communication to other facilities, also in other wavelength bands, is vital. To this end the NRTA pipeline was developed. Its purpose is to analyse the science data on ground as soon as it is ready at MPE and alert the appropriate team of scientists for a given event who can then decide the correct course of action. Beyond this, the NRTA also aids in some technical aspects of the maintenance of the eROSITA instrument.  \nOne of the highest priorities is the identification of bright Xray transients, especially if they are not known. To this end, the NRTA runs an implementation of the Bayesian Blocks algorithm developed by Scargle et al. (2013) on the raw detector count rates of TMs 1, 2, 3, 4, and 6. The algorithm is disabled for TMs 5 and 7 because their optical light leak causes strong fluctuations in the count rates. When a bright source passes the field of view a significant excess of the count rate is expected which can be recognised by the algorithm. This excess is expected to have a duration of roughly 40 s. After additional filtering to exclude such periods of increased count rate caused by artefacts from the CE, each of these time windows is marked for source detection. An exemplary detector light curve of TM 6 for a bright source passing through the field of view and the resulting Bayesian Blocks is shown in Fig. C.1.  \nFor all sources, either detected by source detection or ingested as point sources externally into the pipeline, a variety of products are extracted using the srctool task from the eSASS, most notably the spectrum and the light curve for the pass of the source through the field of view. Based on these data and possibly the information from the source detection, additional custom parameters tailored for specific science cases are calculated, like count rates in di ff erent energy bands, the signal-to-noise ratio,  \nFig. C.1. Example detector light curve (blue) of TM 6 of one telemetry file with bins 5 s in length. The segmentation created by the Bayesian Blocks algorithm is shown in red and the identified intervals of bright sources in the field of view are shaded green. Note that the very short block at t ≈ 1400 s was caused by a remaining corrupt frame. The block of elevated count rate towards the end of the light curve was caused by an unknown transient X-ray source detected on New Year's Eve 2019 (Wilms et al. 2020).  \n<!-- image -->  \nhardness ratios, and the minimal distance of the source to the centre of the field of view. To compare the sources with entries in other external catalogues, the Nway Bayesian algorithm (Salvato et al. 2018) is used to match the positions of the sources against 54 catalogues covering all wavelengths bands.  \nThe notification to the science teams is based on alerts. An alert for a source is generated if a given set of criteria, generically called 'triggers', which can be combined using logical and , or , and not operators, holds true. A variety of these criteria are implemented in a flexible manner and can be used, for example, to test whether a numerical parameter for a source falls within a given range, or equals a specific value, or has a match in a specific catalogue with a specific value. A blacklist is used to suppress alerts from areas which are known to produce large amounts of spurious detections, for example from around Sco X-1 (compare Fig. 6). If a trigger is evaluated as positive for a given source, an e-mail notification with basic information about the source is sent to the science team responsible for the trigger. The alert can be inspected in a web-based front-end which can display more advanced information about the source (images, light curves, and external surveys). A team of two scientists is assigned in weekly rotations to inspect the generated alerts and communicate with the science teams. Some triggers are automatically closed by the system and stored for later bulk analysis. All results of the NRTA are stored on a separate archive and available for later inspection. Typically, it takes less than six hours after the data have arrived for processing until the generation of an alert, which is su ffi cient given that data can be stored on-board for up to 24 hours before being transmitted to the ground station.  \nDuring eRASS1 the NRTA analysed the occurrences of 21.4 million possible sources. Roughly one million of these were found during source detection in regions given externally or identified as containing a bright source, and the remaining ones were ingested as known point sources to monitor. In total, roughly 150 000 alerts were generated. Out of these, we identified several time critical events, which were published mostly as Astronomer's Telegrams (Rutledge 1998). The first such event  \noccurred shortly after the begin of the survey on December 31, 2019 and was caused by a bright, and yet otherwise unidentified, X-ray transient designated SRGt J123822.3-253206 with a flux of 2 × 10 -10 erg s -1 cm -2 in the 0.2-10 keV band (see Fig. C.1, Wilms et al. 2020). Another bright transient, SRGt J071522.1191609, with a flux of 1 . 3 × 10 -11 erg s -1 cm -2 , was detected on April 14, 2020 (Gokus et al. 2020). On May 22, 2020, the NRTA detected a flare from the millisecond pulsar PSR J1023 + 0038 with a flux of 3 . 6 × 10 -11 erg s -1 cm -2 in the 0.3-10 keV band (Koenig et al. 2020). On June 6, 2020, the Be / X-ray binary RX J0529.8-6556 in the LMC was found in outburst by the NRTA for multiple consecutive scans with an average flux of 2 . 5 × 10 -11 erg s -1 cm -2 in the 0.2-8 keV band (Haberl et al. 2020; Treiber et al. 2021).", 'Appendix D: Column descriptions': 'Table D.1 contains a description of the catalogue data model. Our single-band and three-band detections result in the 1B and 3B catalogues, which have di ff erent columns. Each 1B and 3B source has a unique source ID, namely, UID , or DETUID , which are equivalent. The UID is only unique in the 1B or 3B catalogue, and most of the 1B and 3B sources are identical. The 1B3B association results (§ 5.3) are saved in additional columns of UID\\_1B or UID\\_Hard . The 1B catalog is divided into the Main and the Supp catalogues, which have identical columns. The Hard catalogue is a subsample of the 3B catalogue.  \nMost of the source information is calculated in multiple bands and thus a band su ffi x is printed in the column names. Tables D.2 and D.3 provide a dictionary for the energy band definitions in the 1B and 3B catalogues, respectively. A few columns are from the PSF-fitting of source detection using the ermldet task, e.g., DET\\_LIKE , ML\\_RATE , ML\\_CTS , ML\\_FLUX , EXT , and EXT\\_LIKE . In the 3B detection, the band index 1, 2, and 3 indicate 0.2-0.6, 0.6-2.3, and 2.3-5 keV, and band 0 indicates allband summary values. In the 1B detection, only band 1 (0.22.3 keV) is involved, and the all-band summary value (band 0) is identical to that of band 1. For the 1B catalogues, we performed forced PSF-fitting ( ermldet ) and forced aperture photometry ( apetool ) at fixed source positions. They are reported in the bands P1-P9 (Table D.2). Combining band P1 (0.5-1 keV) and P2 (1-2 keV), we calculated the values in band "S" (0.5-2 keV, see § 5.4).  \nTable D.1. eRASS1 catalogs column desciption. All errors are provided as 68% confidence intervals (1σ ).  \nTable D.2. Dictionary of energy band su ffi xes in the eRASS1 1B (Main and Supplementary) catalogsTable D.3. Dictionary of energy band su ffi xes in the eRASS1 Hard catalog', 'Appendix E: Software and calibration versions used in this work': 'The eROSITA data presented in this work were processed in the time period Nov. 2021 to Jan. 2022. The software and calibration versions used (pipeline configuration 010) are therefore di ff erent from those of the earlier eROSITA early data release (EDR, pipeline configuration 001), as well as from later work based on proprietary eROSITA data (pipeline configuration 020). Here we provide a summary of the main di ff erences between the EDR, DR1 and later software (Tab. E.1) and calibration (Tab. E.2) versions. Further details on eSASS software versions are available on the eROSITA DR1 website 20 .  \nTable E.1. EDR / DR1 eSASS task version changes  \nNotes. Only tasks with functional improvements are listed. Task versions for pipeline configuration 020 are as of 2023-01-15. The EDR dataset only provided calibrated event lists; task versions available in the EDR user eSASS package are listed in square brackets. Task srctool has di ff erent but functionally equivalent version numbers in the DR1 pipeline and in the corresponding DR1 user eSASS (user eSASS version marked in square brackets).  \nTable E.2. EDR / DR1 calibration version changes  \nNotes. Listed are the versions numbers which are used in the respective processings. That does not necessarily mean, that they are used for eRASS1 processing, because their application is only necessary at a later date.  \nDi ff erent versions are separated by " / ". In case a calibration file does not exist for a TM it is marked by "-".', 'Appendix F: Essential Dictionary and List of Acronyms': 'Table F.1 presents a list of useful terms and acronyms found in this paper.  \nTable F.1. Essential dictionary and list of acronyms'}