Cool gas is the most abundant phase of matter in the interstellar medium of galaxies and in the Earth’s atmosphere. Hence, precise knowledge of the spatial variation and dynamics of the physical and chemical properties of cool gas is of crucial importance for astrophysics,1 as well as weather forecasting and global climate modeling.2,3 These properties can be diagnosed using passive terahertz remote sensing, because many molecules and atoms in cool gas make energy level transitions that leave distinct fingerprints in the terahertz emission spectrum.4 However, most of the currently used coherent receivers5 can observe only about 1 to 10 GHz of bandwidth at a time (exceptionally up to 6,7). This makes it very time consuming to cover a significant fraction of the 0.1- to 1-THz band (referred to as “terahertz radiation” in the remainder of this article), where the atmosphere is partially transparent. Indeed, it is often essential to probe multiple energy transitions of multiple chemical species to make a meaningful diagnosis, calling for much wider simultaneous bandwidths.2,8
The on-chip filterbank spectrometer9–12 is a recently proposed concept that aims to enable wideband terahertz remote-sensing spectroscopy with a threefold advance. The first step is to use an array of incoherent detectors8 to measure the signal, after dispersion, in one detector for each frequency bin. This decouples the observation bandwidth from the detector readout bandwidth, allowing very efficient back-ends for medium resolution spectroscopy. The second step is to use superconducting millimeter wave electronics to: (i) couple the signal to the chip, (ii) disperse it, and (iii) measure it. This integrates the functionality of a classical dispersion spectrometer8 with a typical size of onto a chip of a few , which makes the spectrometer much more scalable toward larger bandwidths, longer wavelengths, and a larger number of spectrometers to form a spectral-imaging array. The emergence of such a passive spectroscopic imager could also reform the landscape of terahertz applications in our daily lives, because it combines the advantages of a broad input frequency band, ideally suited for pressure broadened lines, fast imaging without scanning, real-time material diagnosis even through clothing, and not irradiating the subject.13–16 The last advance is the fact that incoherent detectors allow for extremely sensitive spectrometers, especially for low radiation environments, because they are not subject to quantum noise.8
Here we demonstrate the design, performance, and operation of a terahertz on-chip filterbank spectrometer system, consisting of a spectrometer chip, cryogenic system, and readout electronics. The spectrometer chip, as shown in Fig. 1, covers a frequency band of 332 to 377 GHz with 49 spectral channels with a constant spectral resolving power of . It is based on incoherent NbTiN-Al hybrid microwave kinetic inductance detector (MKID) technology17 that can be scaled to near the NbTiN gap frequency of . The sensitivity we demonstrate here is set by the photon noise inherent to the incoming optical signal, and by the instrument radiation coupling efficiency. Using large bandwidth antennas such as the leaky wave antenna18 and the sinuous antenna,19 input bandwidths in excess of one octave are possible (e.g., 300 to 900 GHz), allowing truly wideband spectroscopy with up to channels per spectrometer.
Filterbank Spectrometer Chip
The chip is fabricated from a 100-nm-thick NbTiN film, which is deposited on the c-plane sapphire substrate using reactive magnetron sputtering.20 This film has a critical temperature , which implies that the material has low losses up to near the gap frequency ( is the Boltzmann constant, and is the Planck constant). The chip couples to linearly polarized radiation using a double slot antenna, patterned using UV contact lithography and plasma etching in a plasma, followed by an in situ cleaning.21 The antenna is placed in the focus of a extended-hemispherical Si lens,22 which is antireflection coated with a -thick layer of Parylene-C.23 The received signal is guided to the filterbank through the terahertz line: a coplanar waveguide (CPW) patterned in the NbTiN ground plane. The filterbank sorts the wideband signal into 49 sub-bands by means of narrow band pass filters. These filters are coplanar, interdigitated resonators, as shown in Fig. 1(f), which are coupled to the terahertz line on one side, and to a NbTiN-Al hybrid MKID24 on the other side. At the resonant frequency of the filter, maximum signal power is transferred from the terahertz line to the MKID. Each MKID is a CPW quarter wavelength resonator, with a resonant frequency in the order of 5 GHz. Its open end is coupled to the readout line, and its shorted end is coupled to the terahertz filter. Near the terahertz filter the MKID CPW has a central line made from 40-nm-thick aluminum with a resistivity of and . For the signal frequencies of 332 to 377 GHz, the aluminum acts as a radiation absorber, because is smaller than the signal frequency band of the spectrometer. The radiation power absorbed creates a proportional shift in the resonant frequency of the MKID, which is read out as a change in the transmission phase of a microwave tone in the readout line that connects to all MKIDs. The ground planes of the readout line and the terahertz line are balanced using aluminum bridges with dielectric bricks fabricated from spin-coated polyimide LTC9505 from Fujifilm.21 As a reference for the terahertz signal power at the input of the filterbank, three “wideband MKIDs” are placed before the filterbank; these are MKIDs that are weakly coupled to the terahertz line directly without a filter, with a nearly constant power coupling of over the relevant frequency range. Similarly, another three wideband MKIDs are placed after the filterbank for a reference of the power that runs through the filterbank without being drawn out by the filter channels. Furthermore, there are four “dark MKIDs” that are placed away from the terahertz line as a reference for signal power coupled to the MKIDs by way of stray light or surface waves.25 The chip is equipped with a mesh of 40-nm-thick -phase Ta on its backside to reduce the propagation of stray radiation inside the chip. The terahertz signal line is terminated by a CPW with an aluminum center line, which absorbs the remaining terahertz signal after the filterbank, to prevent reflection of power that is not drawn out by the filter channels.
The filterbank spectrometer chip is integrated into a system, which includes the detector readout using microwave frequency division multiplexing for simultaneously reading out all detectors, optics to couple the terahertz signal to the chip, and a refrigerator for cooling them. The spectrometer chip is mounted in a housing as shown in Fig. 2(a). The optical entrance of the chip housing is closed off above the lens with a 1-THz low-pass filter, and a 350-GHz bandpass filter (BPF) with a transmission as shown in Fig. 4(e). The chip housing is mounted inside a light-tight chamber equipped with coaxial cable filter feedthroughs.26 The chamber is cooled via copper thermal straps to 120 mK using a two-stage adiabatic demagnetization refrigerator, connected to the 4-K stage provided by a two-stage pulse tube cooler. The chamber itself is isolated from the 4-K environment by a thermal-mechanical structure based on Vespel tubes. Both the chamber and chip housing operated at 120 mK are coated on the inside by a radiation absorber that consists of 3% by weight carbon powder mixed with Epotek 92 epoxy, in which we embed 1- or 0.5-mm SiC grains as diffusive scattering elements. This recipe is adapted from Ref. 27 to yield better absorption at long wavelengths and better adhesion on large surfaces. A microwave analog/digital readout system28 and a 4-K low noise amplifier are used to read out the phase response of all of the MKIDs simultaneously. This is achieved by creating a time-dependent signal that is the reverse Fourier transform of a set of single frequency tones, one for each MKID, with the addition of a few “blind” tones used for system calibration and monitoring. These tones are sent to the chip, where radiation absorbed in an MKID modifies its associated single frequency tone. After passing through the chip the signal is amplified at 4 K using the low noise amplifier with a noise temperature of 5 K, after which it is further amplified and analyzed in the readout system. Here, the change in complex transmission at each tone is obtained using an on-board fast Fourier transform (FFT) engine, and it is converted to a resonance frequency response in the postanalysis to increase the linearity.29 Radiation is coupled to the detector chip using a 4-K light-tight box, equipped with an optical relay of two parabolic mirrors, baffling structures coated with the same radiation absorber as the 120-mK box and sample holder, infrared filters, and a polarizing filter mounted in copolarization with the antenna, as shown in Fig. 2(b). The relay creates a pupil where the beam is tapered at the level, limiting large angular stray radiation. Even in a single beam system, the pupil has a specific meaning in optics, being the conjugate plane of the focal plane that limits the angular extent of the beam. Complete baffling to limit the optical throughput incident on the chip in both area extent and angle therefore requires limiting apertures (stops) in both conjugate planes (focal and pupil planes). Here aperture baffles at the window, on the array and the lens itself act as an effective focal plane aperture stop. While at the pupil, there is an aperture limiting the angle to an ratio of 3.9 (7.7 deg). This entire optical chain is cooled to 4 K. In addition, the chip is surrounded by two layers of magnetic shielding (Nb and Cryophy30), also cooled to 4 K. The combined shields attenuate the magnetic field perpendicular to the chip by a factor .
The frequency response of the spectrometer is measured with a photomixing terahertz continuous wave (CW) source (Toptica Terabeam 1550) in a setup as shown in Fig. 3(a). We couple the beam from the CW source to the cryostat by reflecting it from a -thick sheet of Mylar to attenuate the signal power by to . We then measure the response of each MKID while sweeping the frequency of the CW source from 320 to 380 GHz. In all results, we apply a frequency shift of to the CW source frequency as a result of a frequency calibration of the CW source using a methanol emission line as will be discussed in Sec. 6. Figure 3(b) shows the response of one filter channel MKID near the filterbank entrance, which has less influence from neighboring channels and transmission losses compared to filters further downstream. The calculation of the response takes the following steps: (i) divide the frequency-dependent response of the filter-MKID with the response of a wideband MKID at the entrance of the filterbank (see Fig. 1); (ii) normalize the peak of the response to unity to calculate the instrument optical efficiency (), see Fig. 4(e); (iii) return to the frequency-dependent response and adjust it so that the peak height is equal to divided by the losses between the entrance of the filterbank to the cryostat window. In this way, the only result from the photomixing-source scan is the relative frequency response, and the peak height of correctly represents the magnitude of the scattering parameter according to the definition in the inset of Fig. 3(b). A Lorentzian fit to the measured transmission peak yields a quality factor of . From the measured peak and , we can infer an internal and a coupling of , assuming the same coupling strength to the feed line and to the MKID.10 The simulation using CST Microwave Studio shows that the filter has a measured transmission to the MKID of on resonance, which is consistent with the measured value. The peak transmission is smaller than the theoretical maximum transmission of 50% for a single BPF.31 According to the power budget around the resonance frequency simulated with CST Microwave Studio32 shown in Fig. 3(b), this can be understood as a result of the following contributions: of the power passing through the terahertz signal line; being reflected back toward the input of the filterbank; 13% being lost as radiation from the filter into the substrate; and 10% being absorbed by the aluminum bridges; leaving to couple to the MKID. Note that the simulated coupling numbers do not add up to 100% but to 97%, due to the total power budget from the CST simulation. This uncertainty is much smaller than the error in the absolute efficiency measurement.
The measured transmission from the filterbank entrance to the MKIDs for all 49 channels is plotted in Fig. 3(c). With being the frequency spacing between adjacent channels, the peak frequencies of the filters follow well a geometric sequence of , as shown in Fig. 3(d). The convex trend seen in the small residual as presented in Fig. 3(e) is likely due to a gradient in sheet inductance of the NbTiN film along the filterbank,20 and a scatter on the order of 0.1% of the resonance frequency as seen here is typical for the readout frequency of MKID arrays.22,33 In Fig. 3(f), we present the quality factors of all channels, showing an average . The frequency spacing of is an intentional oversampling compared to the mean filter of 300, to conservatively avoid gaps in the spectral coverage (this resolution is well suited for the detection of redshifted emission lines from terahertz-bright galaxies34). As a consequence of the power being partially shared among neighboring channels, the oversampling causes the peak transmission of many channels to be slightly lower (mean of 8%) than the highest-frequency channels close to the entrance of the filterbank, such as the one shown in Fig. 3(b). Collectively, the sufficiently small scatter in both filter frequency and with a channel yield of 100% makes it possible to cover the entire 332- to 377-GHz band without significant gaps, as shown in Fig. 3(c).
Because the MKID detectors used in this study are photon-noise limited in the relevant range of terahertz loading power,24 we can measure the sensitivity of the spectrometer using the stream of photons as an absolute calibration source. A pair of black body radiators at room temperature and liquid temperature are placed at reflection and transmission positions of a wire grid seen from the spectrometer window, respectively, as shown in Fig. 4(a). Looking into the wire grid, the spectrometer sees an effective load brightness temperature of , where is the reflectance in copolarization to the cryogenic wire grid and on-chip antenna.
Figure 4(b) shows the relative frequency response ( is the MKID resonant frequency for ) of a representative filter-channel MKID, as a function of . From a square-root-law fit (expected if the responsivity is scaling as , where is the absorbed power,35 and the total optical loading is dominated by the variable thermal load outside of the cryostat), we determine the temperature response . To be able to measure the spectrometer sensitivity expressed in the noise equivalent power (NEP), we convert the load temperature to radiation power , which is defined as the single-mode, single polarization radiation power outside of the cryostat window that can couple to a single filter channel:
Here, is the transmission of the filter channel, as shown in Fig. 3(c), with the peak transmission normalized to unity. Note that we have limited the integral bounds of Eq. (1) to the bandwidth of the quasioptical filter to avoid integrating noise as signal. Furthermore, is the speed of light in vacuum, and is the Planck brilliance as a function of and frequency . The NEP can now be determined experimentally by combining the responsivity with MKID photon noise level usingFig. 4(c). The amplitude noise shown in Fig. 4(c), at a level of compared to the phase noise, corresponds to the noise in the readout system,28 proving that the phase noise level is not affected by the readout noise. It has to be noted that both and depend on background load . We plot in Fig. 4(d) the experimental NEP obtained using Eq. (2) at . In the range of 345 to 365 GHz, where the quasioptical filter stack has its highest transmission of , the instrument has an optical NEP of .
Here, is the Planck constant, is the Boltzmann occupation number of the radiation outside of the cryostat window, and is the pair-breaking efficiency.38 We show the optical efficiency in Fig. 4(e), together with the band-pass characteristics of the quasioptical filter stack. The figure shows that the instrument optical efficiency in the passband of the filter stack is . This is close to the product of the following transmissions: (1) quasioptical filter stack with a transmission of , (2) coupling of the cryogenic optics of , simulated with physical optics using GRASP,39 (3) radiation efficiency of the lens-antenna of , simulated in CST Microwave Studio, (4) ohmic loss of the 30 bridges across the terahertz line in between the antenna and the filterbank of , estimated from an independent measurement of the loss of the bridges, and (5) on-chip filter has a mean peak transmission of as shown in Fig. 3(c).
It is informative to calculate the optical efficiency for the reference MKIDs, by approximating with the full passband of the quasioptical filter stack normalized to unity at the maximum. The results are shown in the right panel of Fig. 4(e). The optical efficiency of the wideband MKIDs at the input of the filterbank is consistent with the designed value of when multiplied with . The wideband MKIDs at the output of the filterbank receive less power than those at the input because of the fraction of power taken out by the filterbank, including radiation losses, ohmic losses, and reflections. The dark MKIDs receive 30 to 40 dB less power than the filter channel MKIDs, showing that the amount of stray light that couples to the MKID detectors is sufficiently small for constructing a filterbank with channels.
Detection of Methanol Gas Emission Spectrum
To demonstrate that the system can spectroscopically observe emission or absorption lines from molecular gas, we have measured the emission spectrum of methanol gas at 17 mbar using the experimental configuration as shown in Fig. 5(a). We couple the system beam to a gas cell that can be filled with methanol gas at 293 K so that the spectrometer looks through the gas cell into a black body cooled to 77 K with liquid . Initially, the gas cell is pumped to a near vacuum of 0.014 mbar. While the spectrometer is continuously observing with all channels at a sampling rate of 160 Hz, we increase the methanol pressure to 17 mbar. The relative frequency shift , where is the initial MKID frequency with the gas cell at vacuum, is converted to an effective brightness temperature by using the response curve as shown in Fig. 4(b) for each MKID. To convert the individual filter channel response to a spectral brightness, we use the measured spectral response of the filterbank, obtained using the photomixing CW source, as shown in Fig. 3(c). The spectral brightness obtained with our spectrometer, is shown in Fig. 5(b). We can compare this response to a simulation of the expected emission spectrum of methanol gas: we present the opacity of methanol gas at 17 mbar by the green dashed line in Fig. 5(b). This spectrum is calculated using the expected line frequencies and intensities from the JPL line catalog40 and taking different line broadening mechanisms into consideration in order to comply with the length (270 mm), pressure, and temperature of the gas cell.41 By further taking into account the losses at the gas cell window and convoluting the intensity spectrum with the bandpass characteristics of the filters as presented in Fig. 3(c), we can calculate the expected response of each channel of the filterbank, as shown in Fig. 5(b). Note that here we have applied to all of the spectral channels a common multiplication factor and an offset, to compensate for the uncertainty in the beam spillover at the windows of the gas cell and at the cold source. Interestingly, we find that we have to use an overall frequency shift of to get a maximum correspondence between the measured filterbank response and the methanol spectrum. This is consistent with the 2-GHz absolute accuracy quoted by the supplier of the CW source and shows that it is possible to use a gas cell as a method to improve the absolute frequency calibration of the spectrometer.
From the good correspondence between measurement and simulation as shown in Fig. 5(b), we can conclude that the spectral shape observed by the on-chip spectrometer reproduces very well the peak frequencies and their relative strengths expected from the database.
Measurement of the System Beam Pattern
The on-chip filterbank spectrometer uses direct detectors, MKIDs, but the radiation coupling via the filterbank and lens-antenna is phase coherent. The coupling between the lens-antenna and finally the telescope is sensitive to the exact phase- and amplitude distribution of the beam, details of which are easily overlooked when doing a conventional scalar beam pattern measurement using thermal sources. To avoid this, and to be able to predict the coupling to the telescope quantitatively, we have measured the phase- and amplitude pattern, using a method pioneered in Ref. 42 and explained in Fig. 6(a). We use two coherent sources, each driven from a signal generator and a frequency multiplier. The two synthesizers have a small frequency difference and a source offset frequency of , well within the bandwidth of the detectors and readout system, which modulates the detector response in the time domain. The complex field parameters can be obtained as a function of the position of RF1 by a complex FFT of the time-domain data. However, MKIDs are not phase-sensitive detectors and the phase information is not conserved in the detection process. To overcome this, we create a stable phase reference by adding a small amplitude modulation to all the readout tones coming from the cryostat at a frequency using the “tone modulator.” The exact technique is described in great detail in Ref. 43.
We show, in Fig. 6(b), the final result of this measurement, which is the far-field pattern calculated from the measured data. We observe a reasonably clean pattern, which can be fitted well to a Gaussian beam, yielding a Gaussicity of 0.82, angular offsets in and given by and , and Gaussian beam far-field divergence angle and , corresponding to the value in power of the best fit Gaussian beam. The reference plane of this direction is the front side of the mounting structure of the cryostat. The cause of the relatively large angular offset is not fully understood but likely a misalignment between the field lens and the beam impinging on it.
Summary and Outlook
We have demonstrated the operation of a full on-chip filterbank spectrometer with 49 spectral channels with a resolution of that observes over an instantaneous band of 332 to 377 GHz. The photon-noise limited sensitivity of , together with an actual detection of methanol gas, proves that the spectrometer can already be applied in atmospheric and astronomical sciences. Indeed, the spectrometer presented here meets all the interface requirements for immediate installation on the 10-m ASTE telescope44 as the first generation of the DESHIMA spectrometer45 for astronomy: with the excellent atmospheric transmission of the Atacama Desert at around 350 GHz, the on-chip spectrometer will more than quadruple the simultaneous bandwidth of existing heterodyne spectrometers at this type of facility,46 matching well the full spectral coverage of direct-detection cameras.
Yet the system demonstrated in this article should be regarded as a narrow-band, single-pixel demonstrator, when compared to the vast potential scalability of the on-chip filterbank spectrometer concept. By using quarter-wavelength BPFs17 and wideband antennas,18,19 the filterbank can be naturally extended to a bandwidth of 1:3 (1.5 octave).9 For a 1:3 broad band system, the out-of-band coupling of the individual filters must be below with respect to the peak power coupling to limit out-of-band power loading the detectors. The stray radiation coupling to “blind MKIDs” in our design is already good enough, but the current filter design has spurious resonances and will not fulfill this requirement. Better instrument sensitivity can be reached by increasing the system coupling efficiency, by improving the coupling efficiency of the chip and the optics. The on-chip filter transmission can be improved from to values approaching unity by a combination of two methods: first, the single filter efficiency is now limited by radiation loss and ohmic losses in the aluminum bridges. This can be mitigated by using microstrip filters and a microstrip THz line, which will both eliminate radiation loss and the need for bridges. This will bring the single filter efficiency close to the theoretical limit of 50%.31 The second step requires a more advanced filter design, an example would be to combine several oversampled filters to absorb more power31 and to incoherently couple that power into one MKID.47 The instrument optical efficiency can be improved by using isotropic substrates (e.g., Si) for the antenna, and by a very careful choice of the infrared and quasioptical filters. These improvements on the chip and optics should collectively bring the instrument optical efficiency up to , corresponding to an instrument NEP of . As we show in Sec. 10, an incoherent spectrometer with this NEP value will have the same per-channel sensitivity as a coherent receiver with a single-sideband receiver noise temperature of : this is very close to the standard quantum limit8 of for coherent detection and a factor of better than the state-of-the-art superconductor-insulator-superconductor receivers for astronomy in this frequency range.46 With the dual advantage in bandwidth and sensitivity, near-future incoherent spectrometers will offer a substantial sensitivity margin over coherent spectrometers for medium resolution spectroscopy ( a few thousands), opening up new observational parameter space in fields such as very wide-band spectroscopy and blind spectroscopic surveys.
Furthermore, the compactness of the spectrometer chip allows small spectrometer units that can be combined into a focal plane array of spectrometer pixels—what is often referred to as a hyperspectral imager or imaging spectrometer. Indeed, the on-chip spectrometer is regarded as the most viable path toward a multiobject spectrometer, which is expected to substantially improve the galaxy-surveying capabilities of existing ground-based observatories for terahertz astronomy48,49 and enable future satellite missions for climatology and meteorology applications.2
Appendix A: Proof of Photon Noise Limited Radiation Detection
The noise spectrum given in Fig. 4(c) is white for both phase and amplitude, which is indicative of photon noise but not a hard proof. Referring to Ref. 50, photon noise limited performance is proven if (i) the noise is white, (ii) has a roll off determined by the quasiparticle lifetime, which depends on the absorbed power. The DESHIMA system as described does not allow to perform such a test due to the limited bandwidth of the readout system. We therefore remeasured the same chip-holder combination in the exact same setup as in Ref. 50. and measured the noise spectra using the cryogenic black body calibrator at temperature of 3 and 30 K. A set of band-pass and low-pass filters defining a 320- to 380-GHz band coupling between calibration load and the DESHIMA holder, the power to an individual MKID, taking into account the DESHIMA filter bandwidth, is 0.015 fW at 3 K and 15 fW at 30 K, the power absorbed in DESHIMA is 20 fW at a 77 K load. The resulted spectrum from one of the resonators is shown in Fig. 7. We clearly observe a white noise spectrum, with a roll-off depending on the absorbed power. At the 3-K black body temperature, the noise is still not quite flat, indicating that two-level-system noise still has a small contribution. At 30 K, the noise is white in the 10- to 100-Hz frequency range used in the article to determine the coupling efficiency. Also, we clearly see a power dependent roll off in the spectrum. This shows that the device is intrinsically photon noise limited at all power levels relevant for the experiment. Note that the amplitude noise shows a similar behavior, but the photon signature exceeds the amplifier noise, given by the level at the highest frequency, only marginally, implying a significant system noise contribution even at low frequencies.
Appendix B: Equivalent Coherent Receiver Noise Temperature of an MKID-Based Incoherent Spectrometer
The calculation is similar to the case in which the incoherent system is completely photon-noise limited,8 except that here we include the relatively small effect of the quasiparticle recombination noise of MKIDs. We will consider a case representative for a spectrometer system installed on a ground-based astronomical observatory, where the signal to be detected has a power that is much smaller compared to a background temperature of typical for the 0.3- to 1.0-THz range.8
The signal-to-noise ratio (SNR) of a coherent receiver is given by the radiometer equation:
Similarly, the SNR for an incoherent spectrometer is given as
Equating using and , we obtain
Note that is the independent of for a photon-noise limited MKID, because the NEP scales with .
As discussed in the article, an on-chip filterbank spectrometer with an instrument optical efficiency of yields an NEP of . Here, we have taken and , where the factor is the area under a Lorentzian curve whose peak amplitude and full-width-half-maximum are both unity. According to Eq. (6), this NEP would correspond to a of 22 K, which is close to the standard quantum limit for coherent detection.8
We thank Peter Hargrave for suggesting promising applications of the on-chip spectrometer to atmospheric sciences and Klaas Keizer for the precise mechanical work on the cryostat. This research was supported by the Netherlands Organization for Scientific Research NWO (Vidi Grant No. 639.042.423, NWO Medium Investment Grant No. 614.061.611 DESHIMA), the European Research Counsel ERC (Consolidator Grant No. 648135 MOSAIC), and the Japan Society for the Promotion of Science JSPS (KAKENHI Grant Nos. JP25247019 and JP17H06130). P.J. de V. was supported by the NWO (Veni Grant No. 639.041.750). T.M.K. was supported by the ERC Advanced Grant No. 339306 (METIQUM) and the Russian Science Foundation (Grant No. 17-72-30036). N.L. was supported by ERC (Starting Grant No. 639749). J.S. and M.N. were supported by the JSPS Program for Advancing Strategic International Networks to Accelerate the Circulation of Talented Researchers (Program No. R2804). B.M. was supported by the European Union’s Horizon 2020 research and innovation program under Grant Agreement No. 730562 (RadioNet). The authors have no relevant financial interests in the manuscript and no other potential conflicts of interest to disclose.
Akira Endo is an assistant professor in the Terahertz Sensing Group of TU Delft. He is interested in studying galaxies at high redshift using superconducting detectors. He obtained his PhD in astronomy from the University of Tokyo in 2009. He has been focusing on the DESHIMA experiment since 2009, when he began a postdoc at Kavli Institute of Nanoscience Delft. In his current position (since 2014), he also enjoys teaching astronomical instrumentation and systems engineering.
Kenichi Karatsu is a postdoctoral researcher at SRON. He received his PhD in science in 2011 from Kyoto University with study of proton spin structure at the RHIC-PHENIX experiment. In 2015, he joined the DESHIMA project at TU Delft/SRON. His main role is to lead laboratory evaluation and telescope campaign of the instrument. His research interest is to develop an experimental instrument for revealing mysteries of the universe.
Alejandro Pascual Laguna received his BSc degree in telecommunications engineering from the Universidad Pontificia Comillas, ICAI School of Engineering, Madrid, Spain, in 2014. In 2016, he received his MSc degree (cum laude) in electrical engineering from Delft University of Technology (TUD), Delft, The Netherlands, where he is currently working toward his PhD with the Terahertz Sensing Group. His research interests include on-chip solutions for efficient broadband sub-mm wave spectrometric and imaging systems based on kinetic inductance detectors.
Behnam Mirzaei received his PhD from the Quantum Nanoscience Department of TUD in 2018. His focus was on the application of quantum cascade lasers as local oscillators (LO) in super-THz heterodyne receivers. He is currently contributing to the development of the receiver unit of GUSTO (Galactic/extragalactic Ultra long duration balloon Spectroscopic Stratospheric THz Observatory), as a postdoc in TU Delft together with a team in SRON (the Netherlands Institute for Space Research).
Robert Huiting is a mechanical design engineer. He received his BS degree precision engineering in 2001 and technology management in 2003 from the Hogeschool van Utrecht. He has been working as a design engineer for SRON Netherlands Institute for Space Research since 2007. He has worked as a research engineer at the FOM institute for plasma physics for the XUV Optics Group the years before. His experience from the last 6 years lies in designing instrumentation for ground-based space research.
David J. Thoen is a cleanroom engineer working at TUD. He received his BS degree in applied physics from Fontys University of Applied Sciences, Eindhoven, in 2008 while he worked as a microwave engineer at FOM Institute Rijnhuizen since 2007. Since 2010, he has been working at TUD on the development of Microwave-kinetic-inductance-detectors far-infrared detectors for astronomy. He has extensive experience in cleanroom processing, process development, vacuum, and cryogenic technology.
Vignesh Murugesan received his MSc degree in microsystem integration technology from Chalmers University of Technology in 2007. He worked as a process integration engineer from 2007 to 2008 for Infineon Technologies, Germany. From 2010 to 2013, he worked as an MEMS process engineer for Thermo Fisher Scientific, Netherlands. Since 2013, he has been working as a process engineer at SRON Netherlands for MKIDS Group. He is currently responsible for the fabrication and process development of MKID chips.
Stephen J. C. Yates has been an instrument scientist at SRON, the Netherlands Institute for Space Research since 2006. He received his PhD from the University of Bristol in 2003, followed by work at the CNRS-CRTBT (now Institut Neél) Grenoble on experimental low temperature techniques in 2003 to 2006. His current interests concentrate on MKID development for sub-mm astronomy applications but also include a wider interest in device physics and superconductivity, optics, and full end to end instrument characterization and performance.
Juan Bueno graduated in physics from the University of Cantabria in 2003 and received his PhD from the University of Leiden in 2007. In 2008, he was awarded with a NASA postdoctoral position, joining the Jet Propulsion Laboratory, where he pioneered a new type of pair-breaking radiation detector, the quantum capacitance detector. He became an instrument scientist at SRON in 2012 working on the development of far-IR and sub-mm wave kinetic inductance detectors.
Nuri van Marrewijk received his MSc degree in applied physics from the Kavli Institute of Nanoscience, TUD, in 2014. He participated in the DESHIMA project as a PhD student at the Faculty of Electrical Engineering, Mathematics and Computer Science, TUD.
Sjoerd Bosma received his BSc and MSc degrees in 2015 and 2017 (cum laude) in electrical engineering from the TUD, Delft, The Netherlands (TU Delft). He is currently a PhD candidate at the Terahertz Sensing Group at the TU Delft, where he works on leaky-wave lens antenna arrays for submillimeter-wave spectrometers. From September 2018 to February 2019, he participated in the JPL Visiting Student Researcher Program at the Jet Propulsion Laboratory, California, USA.
Ozan Yurduseven is a senior antenna engineer at Huawei Technologies Dusseldorf GmbH, Germany. He received his BSc and MSc (Hons.) degrees in electronics and communications engineering from Yildiz Technical University, Istanbul, Turkey, in 2009 and 2011, respectively. He received his PhD in the Microelectronics Department from the TUD, Delft, The Netherlands, in 2016. His research interests include wideband integrated lens antennas for THz imaging, numerical techniques in electromagnetics and metamaterials.
Nuria Llombart is a professor at the THz Sensing Group of TU Delft, where she has been working since 2012. Her current research interests include the analysis and design of planar antennas, periodic structures, reflector antennas, lens antennas, and waveguide structures, with emphasis in the terahertz range. In 2015, she was the recipient of a European Research Council Starting Grant. She serves as a board member of the IRMMW-THz International Society.
Junya Suzuki is a postdoctoral fellow in the Institute of Particle and Nuclear Studies (IPNS) at KEK. He received his BS degree in physics from the University of Tokyo in 2011 and his MS and PhD degrees in in physics from the University of Tokyo in 2013 and 2016, respectively. He joined IPNS in April 2016 and is involved in the development of mm-wave telescopes.
Masato Naruse is an assistant professor at the Graduate School of Science and Engineering at Saitama University. He received his BS degree in physics from Kyoto University in 2007 and his MS and PhD degrees in science from the University of Tokyo in 2009 and 2012. He joined the Saitama University Faculty in 2012 and develops high-sensitivity superconducting detectors for astronomy and industry.
Pieter J. de Visser is an instrument scientist at SRON, Netherlands Institute for Space Research. He received his BS and MS degrees in applied physics from TUD in 2007 and 2009 with a specialization in astronomy and instrumentation at Leiden University. He received his PhD in 2014 from TUD (cum laude). After a postdoc at the University of Geneva, Switzerland, he now focuses at SRON on energy-resolving detectors for visible/near-infrared light.
Paul P. van der Werf is a professor of extragalactic astrophysics at Leiden Observatory, Leiden University, The Netherlands. His interests are in the interstellar medium of galaxies at low and high redshift, in the cosmic evolution of the galaxy population, and in the fundamentals of science. His observational work focuses on spectroscopy in wavelength regions from the near-infrared to the radio regime.
Teunis M. Klapwijk is at the Kavli Institute of Nanoscience, mostly known for his fundamental contributions to superconductivity. He obtained his PhD in 1977, continued his work at Delft until 1985, with a break at Harvard in 1979 to 1980, followed by a professorship at Groningen. He returned to Delft in 1999, where he became a research-only professor in 2013. He has also been a key-contributor to research on superconductors for astronomical detection.
Jochem J. A. Baselmans is a senior instrument scientist at the SRON Netherlands Institute for Space Research and professor at TUD. He received his PhD (summa cum laude) from the University of Groningen in 2002, titled “Controllable Josephson Junctions.” He leads the Dutch effort on the development of microwave kinetic inductance detectors, where his main interests are ultrasensitive MKIDs for THz radiation detection and advanced on-chip imaging spectrometers for sub-THz imaging spectroscopy.