---
author:
- 'The CONCERTO collaboration: P. Ade , M. Aravena, E. Barria , A. Beelen , A. Benoit , M. Béthermin , J. Bounmy , O. Bourrion , G. Bres , C. De Breuck , M. Calvo , Y. Cao , A. Catalano , F.-X. Désert , C.A Durán , A. Fasano , T. Fenouillet , J. Garcia , G. Garde , J. Goupy , C. Groppi , C. Hoarau , G. Lagache , J.-C. Lambert , J.-P. Leggeri , F. Levy-Bertrand , J. Macias-Perez , H. Mani, J. Marpaud , P. Mauskopf , A. Monfardini , G. Pisano , N. Ponthieu , L. Prieur , S. Roni , S. Roudier , D. Tourres , C. Tucker'
bibliography:
- 'conc.bib'
title: 'A wide field-of-view low-resolution spectrometer at APEX: instrument design and science forecast'
---

Introduction {#sec1}
============

Modern imaging and polarimetry cameras, at millimetre and sub-millimetre wavelengths, are currently operating on large (e.g. D $>$ 10 meters) single-dish telescopes. The main goal of such instruments is to map, at relatively high angular resolution (e.g. 5-30 arcseconds) large portions of the sky (e.g. several deg$^2$) with high sensitivity (e.g. RMS$_{\mathrm{MAP}} \lesssim$ 1mJy). Polarised emissions are also measured with similar specifications.

In this context, the dual-band NIKA2 camera represents the first kilo-pixels instrument operating at these wavelengths based on the Kinetic Inductance Detectors (KID) technology [@nika2_1; @nika2_2]. The particular flavour of KID used for NIKA2 (and CONCERTO) are front-illuminated Lumped Elements KID (LEKID) [@Doyle2010]. They consist in inductor-capacitor (LC) superconducting planar resonators made by a long meandered inductor (wire) terminated at both ends by an inter-digitated capacitor. NIKA2 supersedes previous cameras based on Transition Edge Sensors (referred to as TES bolometers) in the frequency range 150-360GHz, such as MAMBO2 at IRAM [@kreysa1998], LABOCA at APEX [@Siringo2009] and SCUBA-2 at the JCMT [@holland2013].

In order to extend the capabilities of the existing instruments, and open new observational windows of the millimetre sky, the spectral dimension has to be added, without sacrificing the instantaneous field-of-view. The large field-of-view and the mapping speed are, actually, the main asset of single dish telescopes when compared to variable baseline interferometers like ALMA[^1] or NOEMA [@NOEMA2020]. For this reason, we are developing a millimetre-wave low spectral resolution (R=$\nu / \Delta \nu \leq 300$) spectrometer with an instantaneous field-of-view of 20arcminutes. In order to preserve the angular resolution, at frequencies of around 300GHz and assuming a 10-meters class telescope, a focal-plane containing around 2,000 spatial pixels is needed. To achieve these figures, we adopt a room-temperature Martin-Puplett Interferometer (MpI) [@MPI70] coupled to a large field millimetre-wave camera. The instrument, named CONCERTO, has been designed to interface with the Atacama Pathfinder EXperiment (APEX) 12-meters telescope [@Gusten2006]. A pathfinder instrument, named KISS and based on the same concept, has been built by our collaboration and deployed in November, 2018, at the Teide Observatory [@KISS2020].

One important science driver that has motivated our developments is the study of \[CII\] emission line at high redshift. \[CII\] is among the brightest lines originating from star-forming galaxies and a reliable tracer of star formation on global scales. With CONCERTO at APEX, we will map in three dimensions the fluctuations of the \[CII\] line intensity in the reionisation and post-reionisation epoch ($z\gtrsim5$). This technique, known as “intensity mapping”, will allow to answer the questions of whether dusty star-formation contributes to early galaxy evolution, and whether \[CII\]-emitters play an important role in shaping cosmic reionisation. The dedicated \[CII\] survey will provide a (spatial-spectral) data cube in which intensity is mapped as a function of sky position and frequency. The 3-D fluctuations are then studied in Fourier space with the power spectrum. The \[CII\] survey will also be sensitive to the CO intensity fluctuations arising from $0.3<z<2$ galaxies, giving the spatial distribution and abundance of molecular gas over a broad range of cosmic time. \[CII\] intensity mapping is also one of the main goals of CCAT-prime [@choi2020] and TIME [@crites2014], two experiments based on different technologies than CONCERTO: gratings and TES (Transition Edge Sensors) bolometers for TIME, KID and Fabry-Perot interferometers for CCAT-p. In addition to the main \[CII\] survey, we expect CONCERTO to bring a significant contribution in a number of areas, including the study of galaxy clusters (via the thermal and kinetic SZ effect), the observation of local and intermediate-redshift galaxies, and the study of Galactic star-forming clouds. In this paper, we detail the main goals of the \[CII\] intensity mapping and galaxy clusters surveys.

The paper is organised as follow. We present the instrumental concept, design and preliminary results in Sect.\[sec2\]. The discussion includes KID detectors, cryogenics and optics. In Sect.\[sec3\], we describe the first laboratory tests (still on-going for a full characterisation). In Sect.\[sec4\], we present the sensitivity estimates, while Sect.\[sec5\] is dedicated to the \[CII\] intensity mapping and SZ surveys.

CONCERTO instrument {#sec2}
===================

CONCERTO has been specifically designed to fit into the Cassegrain cabin (C-cabin) of the APEX telescope. It is composed of two main components: the so-called “chassis” and the “optics box”. The “chassis” includes the camera (cryostat), the MpI interferometer, the readout and control electronics. The “optics box” includes a number of mirrors and polarisers and a cold reference for the MpI. In Table \[tab1\] we summarise the main instrument characteristics. The location of CONCERTO sub-systems are shown in Fig.\[whole\].

  ----------------------------------------------------- ------------------------
         Telescope primary mirror diameter \[m\]                   12
            Field-of-view diameter \[arcmin\]                      20
          Absolute spectral resolution \[GHz\]                  $\geq$ 1
          Relative spectral resolution R \[\#\]                  1-300
          Frequency range HF $\mid$ LF \[GHz\]           195-310 $\mid$ 130-270
            Pixels on Sky HF $\mid$ LF \[\#\]               2152 $\mid$ 2152
       Angular resolution HF $\mid$ LF \[arcsec\]          20-32 $\mid$ 23-45
   Average angular resolution HF $\mid$ LF \[arcsec\]         26 $\mid$ 34
      Instrument geometrical throughput \[srm$^2$\]       2.5$\times$10$^{-3}$
   Single Pixel geometrical throughput \[sr$\,$m$^2$\]   1.16$\times$10$^{-6}$
                Data rate \[MBytes/sec\]                          128
  ----------------------------------------------------- ------------------------

  : Main characteristics of CONCERTO. For comparison, the NIKA2 instrument (6.5’ field-of-view, 30-meters telescope) geometrical throughput (A$\times\Omega$), characterising the camera collecting power, is 1.7$\times$10$^{-3}$ sr$\,$m$^2$. Concerning the overall optical transmission of CONCERTO compared to NIKA2, we refer the reader to Sect.\[sec4.1\].[]{data-label="tab1"}

In this section, we describe in detail the camera and its content (Sect.\[subsec2-1\]), the MpI (Sect.\[subsec2-2\]), and the chassis and related electronics (Sect.\[subsec2-3\]). The room temperature optics, including the cold reference source, is presented in Sect.\[subsec2-4\]. Section\[subsec2-5\] is devoted to a brief description of the CONCERTO hardware components located elsewhere than in the telescope tower. The installation at the telescope is the purpose of Sect.\[subsec2-6\].

![Location of the CONCERTO sub-systems: the chassis and optics in the C-cabin, the Data AcQuisition computers (DAQ) in the middle container (“instrumentation container”) and the gas handling system and the pulse-tube compressor in the bottom container (“compressors room”). Real Time Analysis (RTA) computers and hard disks are not in the telescope tower.[]{data-label="whole"}](whole.pdf){width="9cm"}

The camera {#subsec2-1}
----------

CONCERTO camera is based on a cryogenic-liquid-free custom dilution cryostat. The dilution insert and the pulse-tube orientation, in particular, have been specifically designed to allow the rotation of the cryostat axis following telescope movements. The cryostat is optimised for the range of telescope elevations (EL) comprised between 30 and 90 degrees. The best working point is achieved for EL=60 degrees.

![CONCERTO camera cross-section and 3-D view (inset). The positions of the three HDPE lenses are shown (L1, L2, L3), together with the Image Stop (IS), the Cold Pupil (CP), the cold polariser (P3) and the two focal planes (FP). The diameter of L1 is around 250 mm.[]{data-label="cryostat"}](Cryostat.pdf){width="9cm"}

The main camera optical features, shown in Fig. \[cryostat\], are the Image Stop (IS) at a temperature of 4K and the Cold Pupil (CP) at the base temperature of around 0.1K. Three HDPE (High Density PolyEthilene) lenses are used in the camera: L1 (room temperature), L2 (4K) and L3 (0.1K). In order to analyse the polarised signal, the last polariser of the MpI (P3) is placed just in front of the LEKID arrays (FP) at base temperature. P3 is a custom wire-grid polariser. It is realised on a 12 $\mu$m-thick Polyimide membrane, and with Copper wires with a pitch of 50$\mu$m. A number of IR-blocking (thermal) and metallic multi-meshes filters [@Pisano] are mounted at different stages. In particular, we have thermal filters on the warmest stages (room temperature to 50K), low-pass multi-mesh filters at the intermediate temperatures (50K to 1K) and band-defining filters at base temperature, i.e. just in front of each focal-plane array. A specially blackened baffle is installed at 4K, between IS and L2, in order to suppress the stray-light.

Since LEKID are sensitive to variations of magnetic fields, they have to be protected by a multi-stage B field screen. Four concentric high-permittivity alloys (i.e. mu-metal and cryogenics variations) cylinders are installed at 300K, 50K (double screen) and 4K. An additional superconducting screen will be wrapped around the focal planes section.

The focal plane arrays are microstrip-coupled LEKID similar to those used for NIKA2 [@nika2_1]. Six excitation/readout lines (feed-lines) are needed to readout each of the 2152 pixels array. A total of twelve pairs of coaxial cables are thus running into the cryostat. The pixels design itself is derived from NIKA [@nika1]. The LEKID details have been optimised to meet the CONCERTO specifications. In particular, the shape of the meander and its coupling quality factor have been adjusted to the target range of frequencies and expected background. The coupling quality factor is designed to be Q$_c\approx\Delta f_{-3dB} / f_0 \approx$ 2.5$\times$10$^4$, with $\Delta f_{-3dB}$ the typical width of the resonance under dark conditions. The thickness of the dielectric substrate is calibrated in order to maximise the quantum efficiency, and is in the range 100-120$\mu$m.

![Picture of the first “100mK block” including the cold pupil, the L3 lens, the P3 polariser and the two arrays (HF and LF) containing 2152 pixels each. The arrays holders are in this case realised in Aluminium. A version of the block with Copper holders is also available.[]{data-label="array"}](Blocco100mK.pdf){width="7.5cm"}

The fabrication process [@dicing2016] looks straightforward when compared to competing detectors having similar performance. The substrate, a 100-millimetres high-purity mono-crystalline Silicon wafer[^2], is prepared in the deposition chamber by a soft ion milling. The superconducting film deposition, Aluminium with thickness of 20nm, is achieved by e-beam evaporation and under a residual chamber pressure of 5$\times$10$^{-8}$ mbars. The deposition rates is fixed at 0.1nm/sec. The UV photo-lithography step is based on a positive resist, and is followed by wet etching. The etching is done using a standard Aluminium etching solution based on phosphoric acid. The diced detectors arrays are packaged in custom holders and bonded, via 17$\mu$m Aluminium wires, to the 50-Ohms micro-strip launchers. Those are then tin-soldered to the inner pin of the SMA[^3] feed-throughs.

The front-end electronics stage is installed in the cryostat, at a temperature of 4K. It is made of a series of twelve low-noise amplifiers (LNA) operating at the resonance frequencies[^4], i.e. in the range 1.5-2.5GHz. A second stage of cryogenic amplification, i.e. twelve commercial LNA, has been added on the 50K cryogenic stage to simplify the room-temperature electronics and reduce its power consumption. The connections between the cold electronics stages, the arrays and the SMA vacuum feed-throughs plate (see Fig.\[cryostat\]), are ensured by commercial semi-rigid cryogenic coaxial cables. In particular, we adopt NbTi superconducting coaxial cables for the portion connecting the output of the LEKID arrays to the input of the front-end amplifiers. Fixed attenuators are mounted on each input line, at the 4K stage. The overall electrical gain of each radio-frequency line to/from the room temperature electronics has been measured and is about +25dB.

The Martin-Puplett Interferometer {#subsec2-2}
---------------------------------

The MpI is a particular kind of Fourier Transform Spectrometer (FTS). It is capable of measuring the differential spectrum of a source, with respect to a given reference. The key elements are three polarisers (P1-beam divider, P2-splitter and P3-analyser), two (fixed/variable) arms, and two (fixed/moving) rooftop mirrors (see Fig.\[MPIreference\]). This technique is widely adopted in the millimetre and sub-millimetre domains, mostly for laboratory characterisations, but also, in the past, for narrow[^5] field-of-view observations from space [@firas99; @Griffin2010]. The first examples of wide field instruments making use of an MpI to obtain spectral information are the stratospheric balloon OLIMPO [@OLIMPO2014] and the already mentioned KISS ground-based spectro-photometre [@Fasano2020].

![Schematics of the MpI concept. Two options are shown for the reference source: a) a de-focused image of the instantaneous field-of-view; b) a cold reference. The polariser P1 provides the needed polarised input to the MpI. P2 is the beam splitter defining the two arms while P3, in the cryostat, dispatches the two projections of the polarised signal to the focal-plane arrays. The incoming beam is represented by a spectral distribution S($\nu$).[]{data-label="MPIreference"}](MpIreference2.pdf){width="9cm"}

To achieve a spectral resolution better than 1GHz, the maximum range $\Delta$l$_{max}$ of the moving rooftop mirror has to larger than 75mm. In CONCERTO, the motors can move by up to 90mm. The range $\Delta$l spanned by the interferogram can be adjusted, on a scan-by-scan basis and depending on the science target, from zero to the maximum. The spectral resolution will then be $$\Delta \nu = \frac{c}{4 \times \Delta l}\,.$$ The optical path difference is thus OPD=$2\Delta l$. Leaving the rooftop mirror stopped at the zero-path difference position results in using CONCERTO as a broad- and dual-band large field-of-view imager.

The distinctive feature of the CONCERTO (and KISS) MpI is the combination of the speed of movement of the rooftop mirror and its size/mass. In order to avoid atmospheric drifts during a single interferogram, the mechanical frequency of the motors is set to around 4Hz, i.e. 8 full interferograms (and spectra) per second are produced by each of the pixels. The lateral size of the mirror to be moved exceeds 0.5meters, for a mass exceeding 3kg. In order to counter-balance the linear momentum associated to such a moving mass, a second motor, with an equivalent mass, oscillates with an opposite instantaneous velocity. With a maximum force around 1000N, the acceleration that can be imposed to the moving mass, including the motor piston, exceeds 100m/s$^2$. The theoretical curve that is commanded to the motor is a square wave, i.e. constant speed for both ways and maximum acceleration at the turn-backs. The real curve is of course smoothed out by the finite acceleration near the extremes. We present a picture of the system that has been built for CONCERTO in Fig.\[MPI3D\].

![CONCERTO double-motors MpI. Two identical linear motors, developing a force of $\geq$1000N each, are acted on in counter-phase to null the total momentum.[]{data-label="MPI3D"}](MpI_picture.pdf){width="9cm"}

The “chassis” and electronics {#subsec2-3}
-----------------------------

The chassis is a single, compact support structure to which many of the core components of CONCERTO are attached. These include the camera itself, the MpI motors and moving mirror, and the electronics boards, along with a large number of modules devoted to the monitoring and control of the instrument. The chassis has been designed and fabricated to match the constraints related to the limited space available in the C-cabin. It allows installing multiple sub-systems of CONCERTO as a single element inside the APEX C-cabin. The chassis is laboratory pre-mounted and can slide through the C-cabin door.

Five microTCA[^6] racks mounted on the side of the chassis host the 12 Advanced Mezzanine Cards (AMC) used to readout the two arrays. The cards have an architecture similar to those used for NIKA2 [@bourrion2012], but they have been improved to be able to generate up to 400 excitation tones spanning 1GHz bandwidth. The data acquisition rate has been increased from less than 100Hz up to 4kHz, in order to properly sample the interferograms generated by the MpI. The calibration strategy has been inherited from the NIKA and NIKA2 instruments [@calvo], but the continuous frequency modulation used there is no longer viable because of the high-sampling rate. As a consequence, in CONCERTO, the position and shape of each resonance circle, which is used to calibrate the data, is reconstructed by sampling three points around the resonance: $f_0$, $f_0+\delta f$, $f_0-\delta f$. The $\delta f$ is much smaller than the resonance width, and of the order of a few kHz. The calibration step is performed at the beginning of each interferogram, while the MpI rooftop mirror is changing the direction of its motion [@KISS2020]. This approach represents in our opinion the best trade-off between optimal calibration and observing efficiency, i.e. the fraction of time that is devoted to science data stream.

The moving elements inside the chassis, in particular the MpI motors and the gas flowing in the Pulse Tube head, generate vibrations which could affect the detector performances. In order to suppress their propagation to the focal plane, the camera is fixed to the chassis only via a series of soft rubber pneumatic actuators, which strongly dampen the vibrations. Furthermore, the pressure inside the actuators (eight in total, with different orientations, as shown in Fig.\[rotation\]) is constantly adjusted by a dedicated software, so that the cryostat position and axis, monitored by means of linear position transducers, are kept constant independently of the telescope elevation. The cryostat position adjustment, requiring a few seconds to complete, is done automatically after each re-pointing and, upon request, between two subsequent observing blocks (scans). The position is, on the other hand, monitored in real time even during scans.

![Rotation of the camera (and chassis) following the telescope elevation (EL). The position of the eight soft rubber pneumatic actuators is shown. Two of them are dedicated to the pulse-tube head. For the six remaining, we indicate (with black arrows) those at action for three representative elevation cases.[]{data-label="rotation"}](Rotation.pdf){width="9cm"}

The cabin optics and the cold reference (optics box) {#subsec2-4}
----------------------------------------------------

The first CONCERTO element along the optical axis, after the telescope mirrors (M1 and M2), is the M3 foldable mirror mounted on the chassis. With a diameter of 900mm, it is the largest among the CONCERTO mirrors. The next mirror M4, attached to C-cabin upper ring, reflects the beam toward the C-cabin floor and directly into the so-called “optics box”. A first virtual image is generated, by the combination M1-M2-M3-M4, before M5. This virtual image plane will be used for some of the CONCERTO qualification tests.

The “optics box” includes a large number of mirrors (M5 to M11), the two polarisers P1 and P2 (Fig.\[MPIreference\]) and the part of the optics providing the cold reference for the MpI. It also includes the fixed rooftop mirror of the interferometer. A general 3-D view is shown in Fig.\[optics1\].

The mirrors are held at the C-cabin temperature, that is regulated at 11$^{\mathrm o}$C=284K. We expect an emissivity of the order of 1% per mirror, equivalent to an additional background of about 3K per surface, so not smaller than 30K in total. The stability of the temperature in the cabin is $\pm$1K. This means that, for an emissivity of 1%, the effective background temperature variation induced per mirror is around 10mK. Considering roughly 10 mirrors, this translates to a total effective background temperature variation of 0.1K. The KID detectors, with a NET of the order of $\approx$mK$\cdot\sqrt{s}$ per pixel, are sensitive to this drift that will produce a correlated signal on all the pixels. However, these instabilities are harmless since they are slower and smaller than the atmospheric fluctuations sitting on top.

![3-D view of the CONCERTO optics, evidencing the M3 mirror interfaced to the APEX sub-reflector (M2). M4 is attached to the ceiling of the C-cabin and represents the only reflective optics component, with M3 and the MpI rooftop moving mirror, outside of the “optics box”. Completing the optical chain: a large number of mirrors (M5 to M11), the two polarisers P1 and P2, and the cold reference optics.[]{data-label="optics1"}](Optics_v16.pdf){width="9cm"}

We performed a trade-off between the requirements related to the image quality (and the interferometry efficiency) and the lateral size of the MpI. In order to obtain a diffraction-limited combined beam for each position of the movable roof mirror (in the range 0–90mm), we imposed the criterion, for each field on the sky, of producing a quasi-parallel beam inside the MpI. According to the geometrical throughput conservation rule, the field-to-field divergence is thus fixed by the diameter of the beam. For the 20arcmin field-of-view, and considering a 12m primary mirror, we obtained an overall beam diameter of about 420mm inside the interferometer. A consequence of this method is that, imposing that the combined beam does not “walk” in the focal plane (see Fig.\[walk\]), we must accept a jitter on the entrance pupil of the optical system, i.e. the “active” portion of the primary mirror. This is, in the end, the main reason why we have decided to under-sample the size of the illuminated primary mirror, to about 11 metres.\

![Simulated focal plane image for the full, 20arcmin, field-of-view. Each spot is shown for the two extreme positions of the roof mirror, i.e. green (0mm) and blue (90mm). The black circles represent the Airy disks. It is clear that the “walking” of the beam is well contained in the diffraction disk.[]{data-label="walk"}](walk.png){width="9cm"}

A remotely controllable three-positions mirror is inserted in the optics chain in order to select the type of reference input for the MpI. The three options are:

-   *Sky:* a de-focused image of the full 20arcminutes instantaneous field-of-view. The distinctive advantage is that the atmospheric common-mode spectrum is optically subtracted, providing a differential measurement of the astrophysics source spectrum with respect to the atmosphere along the line-of-sight. This means that, when targeting a field populated by weak sources, and at the first order, we will obtain a null interferogram. This configuration is ideally suited for compact object, i.e. angular extension smaller than 20arcminutes.

-   *External cold black-body:* a highly-emissive ($\epsilon \geq $0.98) cold disk cooled down by an independent pulse-tube cryostat (T$_{BB}\approx$8K). This configuration is mostly adapted to extended emission observations, i.e. when the spectral and photometric gradients extend on average more than 20arcmin.

-   *Cold cryostat:* a de-focused image of the CONCERTO cryostat cold (inner) parts. In other words CONCERTO “looks” into itself to find a cold auto-reference. Considering the number of optical elements (six mirrors and three lenses) lying between the three-positions mirror and the coldest stage of the cryostat, we expect an equivalent effective temperature of the order of 20K. This is in any case lower than the loading of the sky plus the whole optics train between M1 and the focal plane (which is not lower than 50K). This innovative configuration will be investigated as a simpler alternative to the external cold black body.

We stress the fact that the external cold black-body, or its alternative “cold cryostat”, are not used as spectral calibrators. They simply represent cold, i.e. colder that the combined thermal emission of the atmosphere and the optical chain, references. In contrast to the “Sky” reference case, the interferogram is not expected to be null at first order. We expect in this case potential systematic effects to be minimised by the fact that both the reference and the dominant target (mostly the atmosphere in the common case of weak astrophysical sources) exhibit thermal black-body spectra.

The best choice between the three references will depend on the particular science target, the observing conditions and the still unknown systematic effects affecting this new kind of large field-of-view spectro-photometre. A crucial phase of the on-sky commissioning will be dedicated to investigating this item. We will report in further publications the results of this study, as well as a more detailed description of the CONCERTO MpI spectral reference system that is beyond the scope of the present paper.

CONCERTO hardware outside the C-cabin {#subsec2-5}
-------------------------------------

On top of the elements described above and located in the C-cabin, CONCERTO is made also by modules elsewhere in the telescope tower and beyond, in particular, the commercial pulse-tube compressor (Cryomech CPA289C) and, more interestingly, the dilution cryostat Gas Handling System (GHS) and the Data AcQuisition (DAQ) and Real Time Analysis (RTA) computers.

The GHS is composed of: a) a series of pumps and compressors for the circulation of the $^3$He-$^4$He mixture and for providing compressed air to CONCERTO, and, b) an electronic cabinet hosting a National Instrument CompactRIO[^7] real-time controller and multiple analogue and digital input-output modules. A dedicated Labview-based software is loaded on the CompactRIO. The software continuously monitors the state of the cryostat and controls all the pumps, compressors, valves and actuators of the dilution circuit. It can perform many tasks automatically, such as pre-cooling the system or putting it in a safe mode if the security thresholds are exceeded. It also acts as a server-side program for the client Graphical User Interface (GUI) that is used on remote computers. The GUI allow to easily see the state of the system and control its components (Fig.\[map\]).

![The GUI used to control the dilution refrigerator circuit.[]{data-label="map"}](ConcCryo.png){width="9cm"}

All the pipes and cables that are needed to interconnect the elements of CONCERTO situated in different rooms (C-Cabin, instrumentation and compressors containers) are routed through flexible hoses. The hoses protect the CONCERTO cabling from the environment and guarantee the flexibility required by the movement of the telescope. Inside the C-cabin, all the connections are centralized on a dedicated panel located on the front side of the chassis, thus easing the procedures of plugging and unplugging.

The DAQ system is located in the middle container and is connected to CONCERTO in the C-cabin, producing 128MBytes per second, through five dedicated Ethernet cables. The DAQ consists in two commercial computers with 48GB of RAM and 24 cores each. The disk storage (432TB) and the RTA systems are installed in the so-called “servers room”, located a few tens of meters away from the telescope tower. The RTA computer has 32 cores and 512GB of RAM. The network connection between the DAQ and the disk/RTA is ensured by two 10 Gigabit switches and underground cables.

Installation at the APEX telescope {#subsec2-6}
----------------------------------

The Atacama Pathfinder EXperiment (APEX) telescope is a modified prototype ALMA antenna with a primary mirror diameter of 12meters and a usable field-of-view of about 20arcminutes. The location at around 5100meters a.s.l. on the Chajnantor plateau ensures optimal observing conditions (see Fig.\[atmosphere\]). In particular, the fraction of time showing a precipitable water vapour (PWV) column lower than 2mm is of the order of 70% or more.

![Atmosphere transmission at APEX under 1mm PWV (blue) and 2mm PWV (orange) conditions. The CONCERTO bands for the HF and LF arrays are shown.[]{data-label="atmosphere"}](atmosphere.png){width="9cm"}

The structure of the Cassegrain cabin, and in general the telescope infrastructure, had been designed to host large field-of-view instrumentation. The primary mirror surface has been recently refurbished, and achieves in some conditions a precision of the order of 10$\mu m$ RMS. APEX is thus, to date, a state-of-the-art installation for millimetre and sub-millimetre Astronomy. The telescope has hosted since 2007, in the same place that will be occupied by CONCERTO, the Large APEX BOlometer CAmera (LABOCA) operating at 360GHz [@Siringo2009].

The CONCERTO optics box and chassis are slided separately through the C-cabin door and then fixed to the floor by a sufficient number of 12mm metric screws. Since the beam will bounce between the floor and the top of the cabin (M4), we have measured the deformations of the C-cabin itself under typical APEX observing conditions. This measurement was achieved using two linear wire sensors with range of 3meters and a single measurement precision of 0.1mm. As shown in Fig.\[deformations\], the RMS of both sensors, over 48hours of data taken during standard APEX observations, is smaller than 0.2mm and thus negligible for our purposes.

![Deformations of a wire sensor running between the floor and the top of the APEX Cassegrain cabin recorded during 48hours. The average absolute length of the wire is 2675.3mm. Inset: picture of the wire sensor.[]{data-label="deformations"}](Deformation.pdf){width="9cm"}

The alignment of the mirrors in the optics box will be achieved in laboratory. The position of each mirror is adjusted with three micro-metric screws. No tuning will be possible at the telescope. The alignment of the optics box with respect to the chassis is ensured by the mechanical fixations. A set of specific lasers will be mounted to achieve in laboratory the internal alignments. At the telescope, we will use these lasers to align the optics box with respect to M4, M4 with respect to M3 and M3 with respect to M2. The alignment procedures will represent a critical step of the installation.

Detectors laboratory characterisation {#sec3}
=====================================

We describe in this section the first tests done on CONCERTO detectors. Some of the electrical tests on the resonances occurred in the CONCERTO cryostat itself. On the other hand, the optical characterisation of the detectors has been achieved in the former NIKA2 test-bench. The so-called NIKA1.5 camera is an easily re-configurable optical dilution cryostat with a base temperature of 60mK. It has been recently modified to host one CONCERTO array at a time. In particular, the optical filters can be easily replaced, and NIKA1.5 can be interfaced to a custom MpI for spectral characterisation, or alternatively to a Sky Simulator (described in detail in @nika1) for sensitivity and beams geometry measurements.

![Frequency sweep (transmission of the feed-line between port 1 and 2, i.e. mod(S21)) of five blocks of resonances out of the baseline HF array. They have been acquired under dark conditions, in the CONCERTO cryostat, at T=70mK.[]{data-label="S21"}](S21.pdf){width="9cm"}

Since the twelve readout lines of the HF and LF arrays share a common Local Oscillator (LO, frequency reference for the readout electronics), it is of vital importance to accommodate all the resonances block in a common $\leq$1 GHz band. This is nicely achieved, for example, in the case of the HF array shown in Fig.\[S21\]. The spread between blocks of resonances belonging to the same array, and between different arrays, is mainly due to inhomogeneities and uncertainties in the thickness of the Aluminium film.

![Quality factors distribution for one representative block of resonances in the CONCERTO HF array. For this particular block Q$_c$=23k+/-12k, in line with all the other blocks and with the designed Q$_c$=25k. These quality factors have been measured in the CONCERTO cryostat at T=70mK.[]{data-label="Q"}](Qc_good.pdf){width="9cm"}

Another important electrical parameter to be studied for large arrays of LEKID is the coupling quality factor Q$_c$. The micro-strip configuration that has been chosen has the advantage of guaranteeing a relatively well-peaked distribution of Q$_c$. This is achieved without requiring complicated and risky additional technology steps like cross-the-line micro-bondings or suspended micro-bridges. Figure \[Q\] shows an example of the statistics obtained for one readout line of the CONCERTO HF array. All the lines, as well as the LF arrays tested so far, exhibit similar behaviours. The quality factor distribution is, as designed, peaking around 25k.

![Spectral response of two CONCERTO HF arrays (solid lines) and one NIKA2 260GHz detector (dashed line). Red: HR Silicon substrate thickness of 110$\pm$5$\mu$m; blue: thickness of 100$\pm$5$\mu$m. Spectral responses have been measured in the NIKA1.5 cryostat with low-pass filters defining an open band up to 300GHz.[]{data-label="Spectra"}](SpectraAll.png){width="9cm"}

The spectral response of two HF arrays, with slightly different substrate thicknesses, has been measured. The results are reported in Fig. \[Spectra\].

![Sky simulator trace. A fake “planet” (point-like source) is crossing the field-of-view of the considered pixel. A raster scan with sub-scans at fixed elevation is simulated. The “elevation” steps are of 4mm each. This measurement has been obtained using the NIKA1.5 cryostat at T=120mK and under a background temperature around 50K.[]{data-label="Planet"}](Planet.pdf){width="9.5cm"}

The sensitivity has been measured in terms of NET (Noise Equivalent Temperature), specific for the NIKA1.5 optics system. The average NET per pixel of the CONCERTO arrays, in NIKA1.5, is around 2mK/$\sqrt{Hz}$. This results in an NET of about 45$\mu$K/$\sqrt{Hz}$ per array (polarisation), or 32 $\mu$K/$\sqrt{Hz}$[^8] when combining both polarisations. Since the sensitivities per pixel are in accordance to what had been measured for the very similar NIKA2 detectors [@nika2_1], we will base our sensitivity estimate in Sect.\[sec4\] on NIKA2 values measured on-sky. We believe actually that the sensitivities measured on the maps on the sky for similar detectors are a more realistic prediction compared to somewhat ideal values estimated in laboratory.

The good imaging characteristics of the CONCERTO arrays are demonstrated by the Sky Simulator tests. An example is shown in Fig.\[Planet\]. A deeper geometrical characterisation of the thousands beams is beyond the scope of this paper, and will be only performed on the final arrays.

The typical response time of the LEKID used for CONCERTO ranges between 30 $\mu$s and 100$\mu$s, depending on the background. Even at the chosen sampling rate of 4kHz, a cosmic-ray hit will thus represent a single-point glitch in the CONCERTO raw-time traces. The order of magnitude of expected rate is 0.1Hz per pixel.

Sensitivity estimates {#sec4}
=====================

Due to the similarities between the NIKA2 and CONCERTO detectors, we use the NIKA2 sensitivity measured on sky, and on reduced maps, as a base to estimate the sensitivity for CONCERTO. As already advocated, we think that this approach, coupled to our NET laboratory measurements, provides quite realistic predictions.

CONCERTO as an photometer: dual-band sensitivity \[sec4.1\]
-----------------------------------------------------------

We first compute the sensitivity for CONCERTO as if it was a dual-band imager (LF and HF). For that, we rely on NIKA2 sensitivity measurements on the IRAM 30-meter telescope. Average NEFDs for NIKA2 (NEFD$_{NIKA2}$) are equal to 9.8 and 36.1 mJysec$^{1/2}$, at 150 and 260GHz, respectively, for pwv=2 and an elevation of 60 degrees [@nika2_2]. These numbers already suffer from the transmission of the whole experiment.\
We observe a large difference between the NIKA2 260 and 150GHz channel performances. A combination of known effects explains the gap in sensitivity. Indeed, at 260GHz, i) the beam efficiency of the 30-meters telescope is about 55%, ii) the sensitivity is strongly affected (by 35%) by a known defect of the NIKA2 dichroïc, and iii) the contribution of residual sky noise to the average NEFDs is important. Therefore, as a realistic starting point for CONCERTO, we assume a sensitivity for the LF array equals to that of NIKA2 at 150GHz. For the HF array, we assume a sensitivity better than that of the NIKA2 260GHz channel, thanks to a gain in beam efficiency at APEX (which is at the order of 80% at the CONCERTO wavelengths), the lack of dichroïc in CONCERTO, and better atmospheric conditions. The values are thus (for pwv=2 and an elevation of 60 degrees): $$NEFD_{NIKA2}^{LF} = 10\,[7.5-15] \,\mathrm{mJy\,sec^{1/2} \,and,}$$ $$NEFD_{NIKA2}^{HF} = 15\,[10-20] \,\mathrm{mJy\,sec^{1/2}.}$$ The numbers in bracket give the uncertainties on our assumption. As it was the case for NIKA2, we make the hypothesis that no excess noise will appear in CONCERTO at APEX compared to CONCERTO in laboratory. Of course, this cannot be verified until installation.

Compared to NIKA2, for CONCERTO, we have to scale the sensitivities to match the APEX telescope size, and we add two polarisers in the optical path (P1 and P2, see Fig.\[MPIreference\]). Sensitivity loss is only due to P1, by a factor between $\sqrt{2}$ (if photon noise dominated) and 2. To be conservative, we consider a factor 2. Thus, for a single array of CONCERTO, the NEFD becomes: $$NEFD^{LF, HF} = NEFD_{NIKA2}^{LF, HF} \times 2 \times \left( \frac{27.5}{11} \right)^2$$ where 27.5 and 11 meters are the IRAM and APEX telescopes effective sizes, respectively, i.e. the portion of the primary mirrors that are optically conjugated to the cold pupils of the instruments (aperture stops).

Then we assume a frequency window of $\Delta\nu$=115GHz, making the assumption that the two arrays cover the frequency range 195-310GHz (HF) and 130-270GHz (LF), with a notch filter removing 25GHz of the low-frequency bandpass (around 183GHz). Finally, we also have to take into account the decrease in transmission due to increase optics complexity of CONCERTO compared to NIKA2 (in particular FTS optics will have some transmission loss and additional loading), which we estimate to be $T=0.8$ for an unpolarised source (this is only an additional loss of transmission compared to NIKA2+P1 and not the overall transmission).

The sensitivity of CONCERTO as a dual-band photometer (set when the optical path difference in the FTS is null) is thus: $$NEFD_{\mathrm{phot}}^{LF, HF} = NEFD^{LF, HF} \times \sqrt{\frac{\Delta\nu_{NIKA2}^{LF, HF}}{\Delta\nu}} \times \frac{1}{\sqrt{T}}$$ We have $\Delta\nu_{NIKA2}^{HF}$=48GHz and $\Delta\nu_{NIKA2}^{LH}$=39.2GHz [@nika2_2] and thus:\
$NEFD_{\mathrm{phot}}^{LF}$ = 81.6 \[61.2 - 122.4\] mJysec$^{1/2}$ and\
$NEFD_{\mathrm{phot}}^{HF}$ = 135.4 \[90.3 - 180.6\]mJysec$^{1/2}$.\
Note that due to the FTS in front of the cryostat, CONCERTO is a non-optimal instrument for imaging. However, the option of removing the first polariser for purely photometric campaigns could be studied. In that case the sensitivity is expected to be 2 times better for each individual array (as we took a factor 2 of penalty for P1).

Sensitivity in spectroscopy
---------------------------

For the spectroscopic mode, we consider a fix value for spectral resolution $\delta\nu$=1.5GHz. The number of spectral elements in the frequency range is $N_{se}$= $\Delta\nu$/$\delta\nu$.

The sensitivity per spectral element (in mJy s$^{1/2}$) for a single spectrometer (note that in our case, with an FTS, the number of pixel equals the number of spectrometers) is given by $$NEFD_{FTS}= NEFD_{\mathrm{phot}} \times N_{se} \,.
\label{Eq_1}$$

The beam area is computed assuming a Gaussian beam, $$\begin{aligned}
\Omega_{\mathrm{beam}} &=& 2\pi \left( \frac{\theta_{\mathrm{beam}}}{2\sqrt{2 \log{2}}} \right)^2\,,\end{aligned}$$ with a FWHM determined by the Rayleigh criterion for a D=11m antenna (our illumination of the APEX 12m antenna) at a given frequency (corresponding to a given redshift for the \[CII\] line), $$\begin{aligned}
\theta_{\mathrm{beam}} &=& 1.22\lambda_{\mathrm{obs}} / D\,.\end{aligned}$$

We can then convert the sensitivity per spectral element from point source (Eq.\[Eq\_1\]) to diffuse emission (in MJysr$^{-1}$s$^{1/2}$) following $$NEI_{FTS} = NEFD_{FTS} \times 10^{-9} / \Omega_{\mathrm{beam}}
\label{Eq_2}$$ This is the Noise Equivalent Intensity, on sky, per KIDS, per spectral bin (taken as $\delta\nu$=1.5GHz).

We can finally compute the mapping speed MS (per spectral element) following, $$\mathrm{MS} = \mathrm{FOV} / NEI_{FTS}^2$$ where FOV is the field of view area (with a diameter of 20arcminutes). Numbers are given in Table\[Tbl\_sens\]. Note that we ignored the frequency overlap between the two arrays (and thus a gain of $\sim \sqrt{2}$ on the sensitivity in the frequency overlap region). We considered $NEFD_{DB}^{LF}$ for $\nu \le$ 150GHz, $NEFD_{DB}^{HF}$ for $\nu \ge$ 260GHz, and a linear interpolation between the two NEFDs for $150<\nu<260$GHz.\
We also give in Table\[Tbl\_sens\] the sensitivity for the whole array, per spectral element, which is: $$\sigma_{\mathrm{array}} = NEI_{FTS} / \sqrt{N_{\mathrm{KIDS}}}\,,
\label{sigma_array}$$ where N$_{\mathrm{KIDS}}$ is the number of pixels (KIDS) of each array (we use 1720 KIDS, which correspond to 80% of valid KIDS in each array). This would be the sensitivity of each voxel[^9] of large maps, assuming a RA-DEC (or AZ-EL) raster scan-like scanning strategy (and assuming pixel sizes of the map equal to beam sizes). Each voxel of the observed map would then be observed by each KIDS. We checked these numbers using a scanning strategy similar to NIKA2 raster scans, with 3 interferograms per beam.

  $\nu$ \[GHz\]                                                 131                  156                211                 238                 272                 302
  ----------------------------------------------------- -------------------- ------------------- ------------------ ------------------- ------------------- --------------------
  Redshift of the \[CII\] line                                  13.5                11.2                8.0                 7.0                 6.0                 5.3
  Beam size \[arcsec\]                                          52.4                44.0                32.5               28.8                25.2                 22.7
  Beam solid angle \[$\times$10$^{-8}$ sr\]                     7.30                5.15                2.81               2.21                1.70                 1.37
  Mapping Speed                                          42.8 \[19.0-76.0\]   19.8 \[9.0-35.9\]   3.4 \[1.8-7.1\]     1.7 \[0.9-3.7\]     0.8 \[0.5-1.9\]     0.5 \[0.3-1.2\]
  \[$\times$10$^{-3}$ deg$^2$ / (MJy/sr)$^2$ / hour\]                                                                                                       
  On sky map sensitivity $\sigma_{array}$                 2.1 \[1.5-3.1\]      3.0 \[2.3-4.5\]    7.3 \[5.1-10.2\]   10.4 \[7.1-14.1\]   14.8 \[9.9-19.7\]   18.2 \[12.2-24.3\]
  \[(MJy/sr) $\mathrm{sec}^{1/2}$\]                                                                                                                         

Low spectral-resolution spectroscopic surveys {#sec5}
=============================================

CONCERTO will offer a generic access to a large FoV and low-frequency resolution spectroscopic instrument. This opening of 3D large-scale surveys is the next step after the broad-band photometric experiments, either from the ground (e.g. LABOCA, SCUBA2, NIKA2) or from space (e.g. Herschel and Planck). The first scientific aim of CONCERTO is to map in three dimensions the fluctuations of the \[CII\] line intensity in the reionisation and post-reionisation epoch (z$\ge$5.3). This technique, known as “intensity mapping”, will measure the clustering of \[CII\] emissivity and allow answering questions on how and when galaxies and quasars formed, and the history and topology of reionisation. Even if \[CII\] intensity mapping has been the basis of instrument definition, we extended the instrument capabilities to make CONCERTO a multi-purpose instrument (e.g., extending the frequency range down to 130GHz for galaxies clusters observations). Thus we expect CONCERTO to bring a significant contribution in a number of areas, including the study of galaxy clusters (via the thermal and kinetic SZ effect), the follow-up of cosmological deep surveys, the observation of local and intermediate-redshift galaxies, and the study of Galactic star-forming clouds. In this section we give some forecast on the expected signal to noise ratio that can be obtained on the \[CII\]-emission power spectrum (Sect.\[CII\_IM\]). In addition, we give some predictions for observing the SZ signal of galaxy clusters (Sect.\[SZ\]).

\[CII\_IM\] \[CII\] intensity mapping with CONCERTO
---------------------------------------------------

\[CII\] is one of the brightest emission lines in the spectra of galaxies. It is an excellent coolant for neutral gas in photo-dominated regions and an extinction-free tracer of star formation at high $z$. Being redshifted into the sub-millimetre and millimetre atmospheric windows for $z>4.5$, it has become one of the most popular line at high $z$. Pointing on known objects, with e.g. ALMA, NOEMA, APEX/FLASH, \[CII\] is now detected in a large number of galaxies at high $z$ ($>$150 at $z>4.5$, with a large contribution from the ALMA ALPINE survey, e.g. @bethermin2020). Such observations are a tremendous step forward but we also need to look at the overall population, i.e. observe large volume and unbiased surveys. First observations with ground-based interferometers e.g., with ALMA (“ALMA Spectroscopic Survey in the Hubble Ultra Deep Field” large program - ASPECS, @walter2016) or JVLA (CO Luminosity Density at High Redshift survey - COLDz, @riechers2019) offer a three-dimensional view of the molecular gas content of galaxies. The covered areas are about 5-60 square arcminutes.\
Intensity mapping complements these efforts beautifully [e.g @kovetz2019] providing an unbiased view of the distribution of the gas that is difficult to assemble from targeted measurements of individual galaxies, and probing cosmological volumes, with maps on several-degree scale and large frequency (and thus redshift) coverage. Intensity mapping exploits the confusion-limited regime and measures the integrated light emission from all sources, including unresolved faint galaxies.\
We will conduct with CONCERTO a major survey of about one square degree with 1,200hours of APEX telescope time. The survey will provide a data cube in which intensity is mapped as a function of sky position and redshift. Our main target is the \[CII\] line emission at $z\ge5.3$. But CONCERTO will also observe the CO intensity fluctuations arising from $0.3<z<2$ galaxies, giving the spatial distribution and abundance of molecular gas over a broad range of cosmic time. The 3-D fluctuations will be studied in Fourier space with the power spectrum.\
To compute the expected SNR on the \[CII\] power spectrum at high $z$, we used the \[CII\] model presented in [@Serra2016]. Using measurements of the Cosmic Infrared Background (CIB) angular power spectra from Herschel/SPIRE together with star formation rate density (SFRD) measurements, they constrain the galaxy FIR luminosity as a function of dark-matter halo mass at all relevant redshifts, in the halo model framework. By using scaling relations from [@Spinoglio2012] to link the intensity of emission lines to the galaxy infrared luminosity, they compute 3D emission line power spectra for all relevant lines, including \[CII\]. They compute the expected SNR of cross-power spectra between \[CII\] and other emission lines, that will constrain the mean amplitude of each signal, and thereby gain insight into the mean properties of the ISM of high-$z$ galaxies. Note that in their paper, they use for CONCERTO a constant sensitivity for all redshifts of $\sigma_{array}$= 155mJysec$^{1/2}$, while we have here $\sigma_{array}$=\[156, 206, 230, 250, 250\] mJy s$^{1/2}$ at z=\[11.2, 8.0, 7.0, 6.0, 5.3\].\
We follow [@gong2012] to compute uncertainties on the power spectra. The observing time per map voxel (considering one pixel equals one beam) is given by $$\begin{aligned}
t_{\mathrm{voxel}} &=& t_{\mathrm{survey}} \frac{\Omega_{\mathrm{beam}} N_{\mathrm{KIDS}}}{A}\,,\end{aligned}$$ with A the survey area, $t_{\mathrm{survey}}$ the on-sky survey time (i.e. 1200$\times$0.7=840hours in our case, considering 30% of overheads), $\Omega_{\mathrm{beam}}$ the solid angle of the beam (Table\[Tbl\_sens\]) and $N_{\mathrm{KIDS}}$, the number of pixels (we consider 80% of valid KIDS, thus $N_{\mathrm{KIDS}}$=1720).

Assuming a spherically averaged power spectrum measurement, and a directionally independent on sky sensitivity, the variance of the power spectrum is: $$\begin{aligned}
\mathrm{var}[\bar{P}_{CII}(k)] = \frac{[P_{CII}(k)+
\bar{P}_{CII}^{\mathrm{N}}(k)]^2}{N_{\mathrm{m}}(k,z)},
\label{eqn:avg_spectrum}\end{aligned}$$ where $\mathrm{N_m(k,z)}$ is the number of modes that leads to the power spectrum measurement at each $k$ and $$\begin{aligned}
\bar{P}_{CII}^{\mathrm{N}}(k) = V_{\mathrm{voxel}} \frac{\sigma_{voxel}^2}{t_{\mathrm{voxel}}}\,,
\label{PCII}\end{aligned}$$ with $V_{\mathrm{voxel}}$ the volume surveyed by each voxel. In the case of CONCERTO, each KIDS gets a spectrum, and considering one KID per beam, we have $\sigma_{voxel}$ = $\sigma_{array}$.

The number of modes at each $k$ is given by $$\begin{aligned}
N_m(k,z) = 2\pi\,k^2\Delta\,k\frac{V_s(z)}{(2\pi)^3}\,,\end{aligned}$$ where $\mathrm{\Delta\,k}$ is the Fourier bin size, and V$_\mathrm{s}$(z) the survey volume, expressed as $$\begin{aligned}
V_s(z) = \chi(z)^2 y_{CII} B_{\nu} A\,,
\label{Vs}\end{aligned}$$ with $$\begin{aligned}
y_{CII}(z) &=& \lambda_{CII}(1+z)^2/H(z),\end{aligned}$$ being the factor to convert the frequency intervals to the comoving distance at the wavelength $\lambda_{CII}$ (rest frame wavelength of the \[CII\] line). In Eq.\[Vs\], $B_{\nu}$ is the considered bandwidth of the measurement. The volume surveyed by each voxel ($V_{\mathrm{voxel}}$ in Eq.\[PCII\]) is $$\begin{aligned}
V_{\mathrm{voxel}}&=&\chi(z)^2 y_{CII}(z)\Omega_{\mathrm{beam}}\delta_\nu \,.\end{aligned}$$

We consider measurements spanning a redshift range $\mathrm{\Delta\,z\sim0.6}$ which corresponds to a frequency range of $\mathrm{B_{\nu}}\sim20$ GHz at $z=6.1$ for the \[CII\] line. Note that [@gong2012] give useful relations for computing $V_s$ and $V_{\mathrm{voxel}}$ for the \[CII\] line.\
Table\[SNR\_PK\_CII\] gives the numbers derived from the above computations and Fig.\[fig\_pk\] shows the predicted \[CII\] power spectrum with its error bars. Because such predictions are very uncertain, we assumed two extreme models, giving respectively low and high SFRD at z$>$3 (see Fig.3 of @Lagache2018_proc). Low SFRD is the pessimistic prediction as it corresponds to the lowest UV-driven SFRD; high SFRD is the optimistic prediction as it corresponds to the CIB-driven SFRD derived from the halo modelling of Planck CIB measurements [@PlanckXXX].\

![Predicted \[CII\] power spectrum at z=6. Three cases are shown, corresponding to three scenarios of SFRD at high $z$ [see @Lagache2018_proc]: High SFRD (dark blue), Low SFRD (light blue) and the geometrical mean of the two (red). Only points with SNR$>$1 are shown. SNRs have been computed considering A=1.4 square degrees, $t_{survey}$=840 hours (which corresponds to a total observation time of 1200 hours taking into account the overheads) and sensitivities as estimated in Sect.\[sec4\]. \[CII\] power spectra have been derived from the modelling of CIB power spectra [@Serra2016], using a conversion from SFR to \[CII\] that conservatively underestimates the \[CII\] luminosity by a factor 6 at z=5 compared to recent semi-analytical models [e.g., @Lagache2018] or ALMA ALPINE measurements [@Schaerer2020]. Our estimates of \[CII\] power spectra are thus likely to be underestimated. Also shown are the predicted \[CII\] power spectra from @yue2019 (using the local SFR-\[CII\] relation, black line) and @chung2020 (dashed black line).[]{data-label="fig_pk"}](Pk_CII_z6.eps){width="50.00000%"}

  Redshift $z$                     5.5               6.2                 7                   8
  ------------------------- ----------------- ----------------- ------------------- --------------------
  SNR P$^{CII}$ mean SFRD     23 \[14-44\]      12 \[6.9-24\]     5.7 \[3.1-12\]      2.0 \[1.0-4.1\]
  SNR P$^{CII}$ low SFRD     4.5 \[2.6-9.8\]   1.9 \[1.1-4.3\]   0.78 \[0.48-1.7\]   0.23 \[0.12-0.48\]
  SNR P$^{CII}$ high SFRD     79 \[57-112\]     55 \[36-87\]       34 \[21-60\]        16 \[8.5-29\]

![image](SZ_spec_sensitivity.jpeg){width="19cm"}

Covering an area of 1.4 square degrees, our survey will provide the first measurements of the \[CII\] power spectrum up to z$\sim$7, considering the mean SFRD and average sensitivity estimate, and up to z$\sim$8 in the best case of sensitivity estimate. Note that the low SFRD case is unlikely as it gives low shot-noise levels for the CIB, not compatible with the Planck and Herschel measurements, and not compatible with current SFRD measurements at high-z based on \[CII\] or far-infrared measurements. Moreover, on top of the exact level of the SFRD at high $z$, the relation used to convert SFR to \[CII\] luminosity is another capital ingredient of such model. [@Serra2016] used the relation of [@Spinoglio2012], which gives, for a given SFR, a \[CII\] luminosity 6 times lower than that obtained with the \[CII\]-SFR relations of [@Lagache2018] or observed with the ALMA ALPINE survey [@Schaerer2020]. Our power spectra are thus very likely to be underestimated. In terms of point source sensitivities, our survey will reach 1$\sigma$=\[14.4, 14.9, 19.7, 22.0, 23.9, 23.9\]mJy, for a spectral element $\delta \nu$=1.5GHz, at $\nu$=\[131.0, 156.0, 211.0, 238.0, 272.0, 302\]GHz, respectively.\
Finally, one major difficulty in intensity mapping surveys is the problem of foregrounds. For \[CII\], the main foregrounds will be the contamination from emission lines from lower redshifts, in particular emission from CO rotation transitions . [@silva2015] and [@breysse2015] showed that this contamination can be partially removed by masking out the brightest pixels in the survey, or the low-redshift galaxies selected from other surveys. For that, CONCERTO will highly benefit from the extensive visible-IR photometry and spectroscopy galaxy survey data that are available in the chosen field (i.e. the COSMOS field). In addition, one of the strengths of CONCERTO is in its wide frequency range: several CO lines are simultaneously observed at the same redshift for all redshifts $z>0.35$. Cross-correlation between these lines will be a very powerful method to remove the contamination. This will be specifically addressed for CONCERTO in a future paper.

\[SZ\] Observing galaxy clusters with CONCERTO
----------------------------------------------

Cluster of galaxies are the largest gravitationally bound objects in the Universe and so they are key for understanding the hierarchical large scale structure [@Kravtsov2012]. The study of galaxy clusters and in particular of their number as a function of mass and redshift allows us to constrain cosmological parameters [@Allen2011]. In the frequency range covered by CONCERTO, cluster of galaxies will be mainly detected via the thermal and kinetic Sunyaev-Zel’dovitch (SZ) effect [@Sunyaev1972; @Sunyaev1980]. The thermal SZ (tSZ) effect [@Sunyaev1972] refers to the interaction of the hot electrons in clusters with the CMB photons, and results in a distortion in the CMB spectrum at the position of the cluster. In the case of hot clusters, it will be also affected by relativistic corrections named rSZ . The kinetic SZ (kSZ) [@Sunyaev1980] is a Doppler shift of the CMB photons induced by the proper motion of clusters of galaxies along the line-of-sight.\
Observations of the tSZ effect have been successfully performed at high angular resolution using continuum cameras based on KIDs, such as NIKA and NIKA2 at the IRAM 30-m telescope [@Adam2014; @Adam2015; @Adam2016a; @Adam2017; @Ruppin2016; @Ruppin2018]. For the kSZ effect, also observed with e.g. NIKA2 [@Adam2016b], a multi-wavelength spectrometer as CONCERTO would be an unique tool to separate the tSZ and kSZ as well as the different foreground components (CIB, CMB and the Galactic emission) and extract information about the cluster physics. Indeed, with sufficiently precise spectroscopy measurements [@Birkinshaw1999], we can measure the cluster mass (from the tSZ effect), proper motion along the line-of-sight (from the kSZ), and temperature (from the relativistic corrections to the tSZ).\
Particularly, considering the angular resolution and the mapping speed of CONCERTO, we expect to perform a precise estimate of the shape of the SZ spectrum for clusters of galaxies for redshifts between 0.2 and 0.8. As an illustration we present in Fig.\[concerto\_sz\] a simulation of a typical cluster that could be observed with CONCERTO. In this simulation the mass of the cluster is equal to 10$^{15}$M$_{\sun}$ and the cluster is located at a redshift z = 0.4. We assume an universal pressure profile model from to compute the cluster Compton parameter map. From left to right we present the cluster SED as expected to be measured by CONCERTO at different radial distances from the center of the cluster. We have considered 6 bands in frequency with typical bandwidths of 10 to 25GHz. We explored both the sensitivity to the relativistic SZ effect (blue points) assuming a cluster temperature of 20 keV, and to the kinetic SZ (red points) assuming the cluster is moving towards the observer with a velocity of 1000km/h. Uncertainties were computed from the sensitivity estimates given in Sect. \[sec4\] assuming a mapping area of 310arcmin$^2$ and a total integration time of 30 hours. We also show in the Fig.\[concerto\_sz\] the individual tSZ, kSZ and rSZ effect contributions.\
We find that for reasonable observation times (tens of hours) CONCERTO would provide first spectral 2D mapping of the intracluster medium of high redshift clusters and should be able to measure cluster velocities via the kSZ effect. CONCERTO should also be able to detect the relativistic SZ effect and measure cluster temperature. A more detailed mapping of the cluster temperature would require observation times of about hundreds of hours.

Conclusions {#sec6}
===========

We have presented the design of the CONCERTO instrument, a novel spectrometer to be installed on the APEX telescope. CONCERTO is based on the development of new arrays in the millimetre using Kinetic Inductance Detectors. It will contain two arrays of 2152 KIDS, mounted in a dilution cryostat that has a base temperature of 0.1K. Spectra are obtained by a fast Martin-Puplett interferometer located in front of the cryostat. Frequency resolution can be up to $\delta_{\nu}$=1GHz. The technological choices leading to the final instrument design have been explained in detail. Promising detectors characterisation have been obtained. Estimates of expected sensitivity are given, based mostly on the NIKA2 experience, i.e. an instrument on sky subject to similar constraints as CONCERTO. The expected sensitivity, combined with the large field-of-view (20arcminute diameter) will provide an unprecedented mapping speed for such an instrument.

CONCERTO will cover the frequency range 130-310GHz, allowing in particular the observation of the \[CII\]-emission line at high redshift and the SZ signal from galaxy clusters. These two main science drivers for CONCERTO have been detailed in the paper. CONCERTO is particularly well suited for these two objectives since i) the \[CII\] line intensity mapping requires a combination of sensitivity, large field-of-view and spectral resolution R$\geq$100 and ii) the study of galaxies clusters via the SZ effect, requires a multi-frequency analysis (to separate the different SZ signals, as well as the SZ signals from CIB, CMB and Galactic dust), the angular resolution of the 12-meters telescope and R$\approx$10.\
CONCERTO is at present in advanced stage of fabrication. The installation and technical commissioning at the APEX telescope is scheduled for the first semester of 2021. The commissioning, science verification and observations are foreseen by the end of 2022.

Besides the authors, the technicians and engineers more involved in the experimental setup development have been Maurice Grollier, Olivier Exshaw, Anne Gerardin, Gilles Pont, Guillaume Donnier-Valentin, Philippe Jeantet, Mathilde Heigeas, Christophe Vescovi, and Marc Marton. We acknowledge the crucial contributions of the whole Cryogenics and Electronics groups at Institut Néel and LPSC. The arrays described in this paper have been produced at the PTA Grenoble microfabrication facility. We warmly thank the support from the APEX staff for their help in CONCERTO pre-installations and design. The flexible pipes, in particular, have been routed under the competent coordination of Jorge Santana and Marcelo Navarro. We acknowledge support from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (project CONCERTO, grant agreement No 788212) and from the Excellence Initiative of Aix-Marseille University-A\*Midex, a French “Investissements d’Avenir” programme. GL warmly thanks Matt Bradford, Jamie Bock and Tzu-Ching Chang for insightful discussions on CONCERTO sensitivity computation and J.-G. Cuby for his help and support for the ERC proposal. We are grateful to our administrative staff in Grenoble and Marseille, in particular Patricia Poirier, Mathilde Berard, Lilia Todorov and Valérie Favre, and the Protisvalor team. We acknowledge the crucial help of the Institut Néel and MCBT Heads (Etienne Bustarret, Klaus Hasselbach, Thierry Fournier, Laurence Magaud) during the COVID-19 restriction period.

[^1]: https://almascience.eso.org

[^2]: https://www.sil-tronix-st.com/en/

[^3]: SubMiniature version A

[^4]: http://thz.asu.edu/products.html

[^5]: In this context, the field-of-view has to be expressed in number of beams. The number of beams is around 2,000 for CONCERTO and 300 for KISS, to be compared to few tens at maximum for the previous instruments.

[^6]: https://en.wikipedia.org/wiki/MicroTCA

[^7]: https://en.wikipedia.org/wiki/CompactRIO

[^8]: 32 $\mu$K/$\sqrt{Hz}$ is equivalent to 22.6 $\mu$K$\, \mathrm{sec}^{1/2}$.

[^9]: A voxel represents a value on a regular grid in three-dimensional space.
