Al/Ti/Al phonon-mediated KIDs for UV-VIS light detection over large areas
L. Cardani N. Casali A. Cruciani H. le Sueur M. Martinez F. Bellini M. Calvo M. G. Castellano I. Colantoni C. Cosmelli A. D'Addabbo S. Di Domizio J. Goupy L. Minutolo A. Monfardini M. Vignati
physics.ins-det physics.app-ph






^1INFN - Sezione di Roma, Piazzale Aldo Moro 2, 00185, Roma -
Italy
^2CSNSM, Univ. Paris-Sud, CNRS/IN2P3, Universite' Paris-Saclay,
91405 Orsay, France
^3Dipartimento di Fisica - Sapienza Università di Roma,
Piazzale Aldo Moro 2, 00185, Roma - Italy
^4Institut Neel, CNRS/UJF, 25 rue des Martyrs, BP 166, 38042
Grenoble, France
^5Istituto di Fotonica e Nanotecnologie - CNR, Via Cineto Romano
42, 00156, Roma - Italy

^6INFN - Laboratori Nazionali del Gran Sasso, Assergi (L'Aquila)
I-67010 – Italy
^7Dipartimento di Fisica - Università degli Studi di Genova,
Via Dodecaneso 33, 16146, Genova - Italy
^8INFN - Sezione di Genova, Via Dodecaneso 33, 16146, Genova -
Italy
angelo.cruciani@roma1.infn.it


The development of wide-area cryogenic light detectors with baseline energy resolution lower than 20 eV RMS is essential for next generation bolometric experiments searching for rare interactions. Indeed the simultaneous readout of the light and heat signals will enable background suppression through particle identification.

Because of their excellent intrinsic energy resolution, as well as their well-established reproducibility, Kinetic Inductance Detectors (KIDs) are good candidates for the development of next generation light detectors. The CALDER project is investigating the potential of phonon-mediated KIDs.

The first phase of the project allowed to reach a baseline resolution of 80 eV using a single KID made of aluminium on a 2x2 cm2 silicon substrate acting as photon absorber. In this paper we present a new prototype detector implementing a trilayer aluminium-titanium-aluminium KID. Taking advantage of superconducting proximity effect the baseline resolution improves down to 26 eV.


Al/Ti/Al phonon-mediated KIDs for UV-VIS light detection over large areas
    L Cardani^1, N Casali^1, A Cruciani^1, H le Sueur^2, M Martinez^3,1, F Bellini^3,1, M Calvo^4, M G Castellano^5, I Colantoni^4, C Cosmelli^3,1, A D'Addabbo^6, S Di Domizio^7,8, J Goupy^4, L Minutolo^3,1, A Monfardini^4, M Vignati^1
    Received: date / Accepted: date
===========================================================================================================================================================================================================================================





§ INTRODUCTION


In the last decade, Kinetic Inductance Detectors
(KIDs)<cit.> underwent a rapid development, allowing their successful
application to millimeter <cit.> and UV astronomy
<cit.>. Different projects are currently on-going with the aim to
apply this technology to other fields, ranging from the search
of extrasolar planets <cit.> to X-ray spectroscopy
<cit.>.

The CALDER project <cit.> is developing large area phonon-mediated
KIDs with the aim of detecting small amounts of ultraviolet-visible (UV-VIS) light. Such detectors will play a
crucial role for the next-generation bolometric experiments searching for
neutrinoless double-β decay. CUORE <cit.>, the current leading bolometric experiment, is an array of 988 TeO_2 crystals read by NTD sensors.
Its sensitivity to double-β decay is mainly limited by the background induced by the α radioactivity of the material surrounding the detector at cryogenic temperature, mainly copper. The possibility to distinguish different particles would allow to overcome this limit. Particle identification could be accomplished by measuring the Cherenkov light (≈ 100 eV) emitted in a TeO_2 crystal by electrons, that are the signal candidates, but not by α particles <cit.>.
CUPID<cit.>, the planned	 upgrade of CUORE, will therefore require light detectors
able to monitor the whole face of a crystal (≈ 25 cm2) with
a baseline energy resolution good enough to clearly detect this tiny amount of
light (< 20 eV RMS).

One of the limits of KIDs is that they can reach a maximum sensitive area of a few mm^2<cit.>. To increase the sensible area to several cm^2 we opted for a phonon-mediated approach<cit.>: photons are absorbed by the silicon substrate on which the KID is deposited and produce athermal or ballistic phonons able to reach the KID, where they break Cooper pairs and generate a signal.
The use of KIDs as phonon-mediated detectors was proposed for the detection of Dark matter <cit.>, cosmic rays <cit.> and X-rays <cit.>.
The application of such detectors to detect visible light was initially studied by the CALDER project, using an array of 4 KIDs
<cit.>, made of aluminium and lithographed on a 2 × 2cm^2, 300 μm thick silicon substrate. Recently it was demonstrated a baseline resolution of 82 eV using a single KID made of aluminium and combining both phase and amplitude readout, that had respectively a baseline resolution of 105 eV and 115 eV<cit.>. In the following we describe a resolution improvement down to 26 eV, obtained using a trilayer KID made of titanium (Ti) and aluminium (Al).




§ DETECTOR DESCRIPTION


KIDs base their working principle on the kinetic inductance of Cooper pairs,  which can be modified by an energy release able to break them into quasi-particles (QPs). If the superconductor is inserted in a resonant RLC circuit with high quality factor (Q>10^3), the density variation of QPs modifies the transfer function S_21 of the resonator both in phase and in amplitude. 
The phase readout is usually preferred since the phase response is larger (up to a factor 10) and is given by <cit.>:

    d ϕ/dE = ηα S_ϕ (ω, T) Q/N_0 V
      Δ_0^2

where η is the energy to QPs conversion
efficiency, α is the fraction of kinetic inductance with respect to total inductance, Δ_0
is the superconducting gap, N_0 is the single spin density of states
and V the detector volume. S_ϕ (ω, T) is a dimensionless factor given by the Mattis-Bardeen theory, describing the phase responsivity as a function of 
frequency and effective temperature <cit.>. Finally Q is the quality factor of the resonator. The detector described in <cit.>, made of a 60 nm thick Al film, featured η=7.4 %, α=2.5%, Q=1.5 × 10^5, Δ_0=179 μeV, V=2.4 × 10^5 μm^3, S_ϕ=2.85, while for aluminium N_0=1.72 × 10^10 eV^-1μm^-3.


One way to improve the responsivity is to increase the fraction of kinetic inductance α and to reduce the
superconducting gap Δ_0, without affecting significantly the other parameters in Eq.
<ref>. This can be achieved by exploiting the
superconducting proximity effect between aluminium and a
lower gap superconductor, such as titanium. Titanium also features a higher kinetic inductance due to its higher penetration depth <cit.> leading to a significant increase of α /
Δ^2_0. However titanium was never demonstrated to work as
a high quality factor resonator, possibly because of its high sensitivity to
contaminants <cit.>, yielding dissipative conducting states.
This leads us naturally to protect titanium with aluminium, which is known to have an almost
lossless self-protecting oxide layer.  Moreover, to efficiently collect phonons from the substrate into the
superconductor it is important to have an acoustic impedance matching between
the substrate and the KID. In that respect, aluminium provides a better coupling to silicon than
titanium <cit.>. The geometry, we present, is thus an Al-Ti-Al trilayer with the following thicknesses: Al 14 nm / Ti 33 nm / Al 30 nm.

The fabrication of the chip was made according to the following steps. The intrinsic 380 μm thick Si <100> wafer is dipped
for 1 min in dilute 5% HF solution and then rinsed in deionized water water for 1 min. The
wafer is then moved to an electron beam evaporator
dedicated to superconducting materials, where the layers are subsequently
evaporated at rates of 0.5 nm/s in a residual vacuum of 4 × 10^-8 mb. This is important especially for titanium
to avoid trapping impurities in the layer <cit.>.





The film is patterned using standard UV lithography followed by a
two-step wet etching. In a first step, the aluminium is removed
using a commercial etchant (H_3PO_4/HNO_3/CH_3COOH). Then, in the second step, the titanium
and the bottom aluminium layers are etched in diluted HF
(0.05%) at a rate of about 0.2 nm/s.

The sensor implements a LEKID design <cit.> identical to the Al detector
in <cit.>. Its active area consists of an inductive
meander made of 30 connected strips of
62.5μm×2mm. The meander is closed
with a capacitor made of 2 interdigitated fingers of
1.2mm×50μm. The silicon substrate
is diced into a 2x2 cm2 detector with the KID in the center as shown in Figure <ref>.


The chip is assembled in a copper structure and held by 4 PTFE supports with total
contact area of about 3mm2. The copper holder is thermally anchored to the mixing chamber of a 3He/4He dilution refrigerator with base
temperature lower than 10mK. 
 



§ ELECTRICAL AND THERMAL CHARACTERIZATION OF THE DETECTOR


The electrical characterization of the KID is based on the study of the
complex transmission (redS_21) of the resonator around the
resonant frequency (Figure <ref>). The frequency dependence of the
transmission of a resonator is ideally well
described by a single-pole approximation. It is however
necessary to take into account also impedance mismatches, power distortion and
electronics non-linearity to perform a rigorous estimation <cit.>. This measurement allows
redone to extract the coupling quality factor Q_c, the internal
quality factor Q_i redand the total quality factor Q = (Q^-
1_c + Q_i^- 1)^- 1.

The resonant frequency is found to be close to 2.4 GHz. Q_c is 1.48 × 10^5 in very good agreement with the value predicted by the SONNET electromagnetic simulation of about 1.5 × 10^5. Q_i is estimated as 6 × 10^5, lower
than our typical samples made of aluminium <cit.>,
but high enough to allow the detector to operate well.

The critical temperature and the resistance of the trilayer
film have been measured in a 4-terminal geometry. The superconducting transition was
measured at (805 ± 10) mK with a normal state resistance per
square R_sq = 0.58  Ω/sq above the transition. We estimate theoretically the T_C of the trilayer following <cit.>, assuming perfect
interface transparency to phonons across the layers and a superconductor total thickness
comparable or lower to its coherence length. Solving the Usadel equations <cit.> within these
assumptions yields an analytical dependence of T_C on the thickness of
the layers, similar to that found in the Cooper model <cit.>.
To apply straightforwardly this bilayer model to the Al/Ti/Al trilayer, we treat it as a bilayer Al-Ti taking the average thicknesses of the Al
layers <cit.>. The expected critical temperature from this model is T_C = 720 mK, assuming values of T_C for the single layers that take into account the increase of critical temperature in thin films, as described in <cit.>. The discrepancy between the expected T_C and the measured value may come from the limitations of the model: we do not know the precise values of the electron-phonon coupling in all the layers, the actual transparency of
the interfaces and also the actual T_C for each layer. A more complete model was recently proposed in <cit.> and its experimental validation will be the subject of a separate work. 




Using the SONNET electromagnetic simulations
software we estimate the amount of kinetic inductance L_k needed to shift the resonant frequency from the value simulated for a perfect conductor (L_k=0) to the measured value. The simulation is calibrated against a detector with identical geometry but in aluminium.  L_k  of the calibration device was calculated independently using direct measurements of T_C,R_sq, film thickness and an estimation of α from resonant frequency shift with temperature.
 We obtain a value of 1.4 pH/sq, which is significantly higher than the value of 1.0 pH/sq obtained using the BCS form described in <cit.>, which depends on T_C and R_sq. This is a first hint that our case of inhomogenous superconductor cannot be well described by the BCS theory.
The obtained value of kinetic inductance translates into a kinetic inductance fraction α=17%, much larger than that obtained in aluminium
(α=2.5%) using the same design <cit.>.







We have also measured the relative variations of the
resonant frequency (f - f_0) / f_0 as a function of the temperature (Figure
<ref>). These data are compared with a model developed
by Gao et al.<cit.> using the Mattis-Bardeen theory. The model is calculated for our best estimation of the parameters (T_C = 805 mK and α =
17%). Again the BCS form cannot describe properly the case of our
inhomogeneous superconductor and we can indeed see the data departing from
the BCS model at intermediate temperatures. The interpretation of these measurements needs a dedicated model, whose development is beyond the scope of this paper. 



§ DETECTOR PERFORMANCE


Our experimental setup allows to back-illuminate the detector using two different sources: a non-collimated X-ray ^55Fe
source and an optical fiber coupled to a fast room-temperature LED.
The ^55Fe source illuminates in an almost uniform way the whole substrate with a rate of 1 Hz.
The LED emits photons at 400nm in the typical range of Cherenkov light. Depending on the LED
pulse duration and/or amplitude, a variable number of photons can be sent onto the substrate. The whole optical system is intially calibrated with an accuracy of 10 % using a PMT at room temperature and then cross-calibrated with the ^55Fe source. The position and size of the optical system spot on the substrate is shown in Figure <ref>.

The resonator is excited and probed with a monochromatic tone at the resonant frequency f_0. The output signal is fed into a CITLF3 SiGe cryogenic low noise amplifier, downconverted at room temperature using a superheterodyne electronics and then digitized
with an acquisition card at a sampling frequency of 500 kSPS.
Time traces of up to 12ms of the real (I) and imaginary (Q) parts of S_21 are acquired following a software trigger. Finally I and Q variations are converted into changes of phase
and amplitude relative to the center of the resonance loop. A detailed description of the experimental setup of our laboratory at INFN
Rome, including the room-temperature electronics and the acquisition software, can be found in references <cit.>.

As already mentioned, the phase response is typically larger than the amplitude one. However the phase noise exhibits often an excess noise, which sometimes makes it advantageous to consider both phase and amplitude responses in the analysis of the detector sensitivity. When amplitude and phase noises are dominated by the noise of the low noise amplifier, they decrease with microwave power as S_φ∝ P_in^- 1 / 2; however at large powers the
non-linearities of the kinetic inductance <cit.> and heating effects <cit.> become relevant and can suppress the signal. The
recombination time τ_QP, which determines the signal integration time <cit.>, is also reduced by heating effects, lowering the signal-to-noise ratio. For these reasons, it is important to
determine experimentally the bias power leading to the optimal signal-to-noise ratio.


For this purpose, we perform a microwave power scan from
about -95dBm to -65dBm at the
input of the sample, while the monitoring signal and noise levels,
as well as the pulse shape, in both the phase and amplitude readouts. A reference signal is produced by sending light pulses of 15.5 keV. The signals are processed using a matched filter, a software algorithm that allows to improve the signal-to-noise ratio by suppressing the signal frequencies that are mostly affected by noise <cit.>. The height of the filtered pulses is reported in Figure <ref> left. The phase
response is quite flat with increasing power around 19 mrad/keV, a sizeable improvement with respect to the 6
mrad/keV obtained with the Al detector in <cit.>. A similar increase is observed for the amplitude response, around 2.5 mrad/keV with respect to the 0.6 mrad/keV obtained with the Al detector. 

In Figure <ref> right we show the signal-to-noise ratio, obtained by dividing the amplitude of the pulses by the RMS of their baseline (after the matched filter). We determine an optimal working power of - 70 dBm, about 10dB lower than what obtained with
aluminium. The efficiency η at the optimal working power is estimated as in <cit.> and amounts to η=(7.5 ± 1) % consistent with what found in <cit.>.





The trailing edge of the pulse is well described by a single exponential law with a decay time that varies between 800 μs and 100 μs with increasing microwave power (Figure <ref>). As shown in <cit.>, this decay constant can be interpreted as τ_QP. The noise power spectrum at optimal power is reported in Figure <ref> for
both the phase  and amplitude readouts. The flat noise observed in the amplitude readout
and in the high frequency region of the phase readout, is consistent with the
noise temperature of the cold amplifier (T_N ∼ 6  K). This noise
is exceeded by a low-frequency phase noise, whose spectral index is found to be 0.49 ± 0.05, consistent with the presence of a two-level system noise in the detector, observed often in KIDs <cit.>. In order to compare the obtained TLS noise with other works, we convert the noise in fractional frequency noise as in <cit.> and get S_TLS(f=1 KHz)=1 × 10^-21Hz^-1, that is among the best values reported in <cit.>. Both power spectra show also peaks due to  the readout electronics. These peaks are strongly suppressed by the matched filter and the measured energy resolution is always found consistent with the predictions from the noise power spectrum and the matched filter.

  



To obtain a solid estimation of the detector performance we implement a calibration based on the Poisson statistics of the photons absorbed into the substrate. The detector is biased at the optimal microwave power of -70 dBm and is illuminated with optical pulses of energies in the range 0 - 3 keV. For each energy the distribution of the pulse heights in the phase readout is found to be well described by a Gaussian with mean m and standard deviation σ (Figure <ref> left).
The σ vs m data are then fitted (Figure <ref> right) assuming that two uncorrelated components contribute to the detector resolution: the baseline resolution (σ_0) and the Poissonian term from the photon counting that we model as √(m × r)[The number of photons absorbed in the detector is N_ph=m / r. According to the Poisson statistics one has σ_ph=√(N_ph)× r =√(m/r)× r=√(m × r).], where r is the responsivity per photon which is estimated from the fit: σ = √(σ_0^2 + m× r). By dividing r by the energy of a single 400 nm photon we
obtain a responsivity of (17.54 ± 0.58 mrad/keV) in good
agreement with the calibration against ^55Fe presented previously in this Section. Using this energy calibration we find a phase readout baseline resolution σ_ϕ = (25.64 ± 0.85) eV. Applying the same procedure to the amplitude readout, we obtain σ_A = (77.3 ± 2.5) eV.



§ CONCLUSIONS AND DISCUSSIONS


In this paper we presented the results obtained using an Al/Ti/Al trilayer film
to develop a phonon-mediated KID for the detection of UV-VIS light on large surfaces.
The sample features T_C = 805mK, Q_i ≃600 k
and a low-power recombination time τ_QP=800μs. The main features of the detector are compared in Table 1 with the Al detector, presented in <cit.>. The lower critical temperature increases significantly the fraction of kinetic inductance α and thus the responsivity of the detector (factor 3 in the phase readout). The integrated phase noise after the matched filter improves of about 25 %. This improvement can be ascribed likely  to the change in the fabrication process: the new trilayer detector was made by etching while the previous detectors were made by lift-off. The overall result is the improvement of the RMS baseline resolution for phase readout from 105 eV down to 26 eV. The amplitude responsivity also improves (factor 4). Since the amplitude noise is dominated by the LNA noise, the lower bias power implies however a worsening of the integrated noise by a factor 3. Therefore the amplitude baseline resolution improves only from 115 to 77 eV.
The obtained energy baseline resolution approaches the target resolution of the CALDER project of 20 eV. The last phase of the project -
currently in progress- consists in the development of light detectors of 5x5 cm2.






§ ACKNOWLEDGEMENTS


This work was supported by the European Research Council (FP7/2007-2013) under
contract CALDER no. 335359 and by the MIUR under the
FIRB contract no. RBFR1269SL. The authors thank the personnel of INFN Sezione
di Roma for the technical support, in particular M. Iannone, F. Pellegrino, L.
Recchia and D. Ruggeri. HLS, MC, JG and AM acknowledge support from the ANR
grant ELODIS2 no. ANR-16-CE30-0019-02. The authors thank also the anonymous reviewers for useful suggestions and comments.




1
  day Day P K, LeDuc H G, Mazin H G, Vayonakis A and Zmuidzinas J 2003 Nature 425 817 
  	
  monfardini Monfardini A et al 2011 Astrophys.J.Suppl. 
  194
  
  mazinMazin B et al 2013 Pub. of the Astronomical
  Society of the Pacific 125 1348 
  
	mazin2 Cook T et al 2015 J. Astron. Telescopes, Instruments, and Systems 1 id. 044001
	
  quaranta Quaranta O, Cecil T W, Gades L, Mazin B and Miceli A 2013 Sup. Sci. Tech.
  26 105021
  
  faverzaniFaverzani M et al 2014 J. Low. Temp. Phys. 
  176 530
  
  vignatiBattistelli E S et al 2015 EPJ C 
  75 353
  
  cuore Artusa D R et al 2015 Adv. High En. Phys. 
  2015 879871

  tabarelliTabarelli de Fatis T 2010 Eur. Phys. J. C
  65 359
 
  casali2Casali N et al 2015 Eur. Phys. J. C
  75 12

  casali3 Casali N 2017 Astropart. Phys.
  91 44

	Wang:2015raa
  Wang G et al 2015 arXiv:1504.03599 [physics.ins-det]
  
	invhie Artusa D R et al 2014 Eur. Phys. J. C 74 3096  
    
	zmuZmuidzinas J 2012  Annual Review of Condensed Matter Physics  3 169  
	
  golwala Golwala S, Gao J, Moore D, Mazin B, Eckart M, Bumble B , Day P, Le Duc H G, Zmuidzinas J 2008 J. Low Temp. Phys.
  151 550
	
	swensonSwenson L J, Cruciani A, Benoit A, Roesch M, Yung C S, Bideaud A and Monfardini A 2010 Appl.Phys.Lett.
  96 263511
  
  cruciani Cruciani A, Swenson L J, Monfardini A, Boudou N, Calvo M and Roesch M 2012 J. Low Temp. Phys.
  167 311
  
  moore 	Moore D C, Golwala S R , Bumble B, Cornell B, Day P K, LeDuc H G and Zmuidzinas J 2012 Appl. Phys. Lett. 100, 262601,
  
  
  cardani Cardani L et al 2015 Appl.Phys.Lett.
  107 093508
  
  cardani2 Cardani L et al 2017 Appl. Phys. Lett. 110
  033504
  
	mazinphd Mazin B, PhD Thesis, Caltech (2005)
	
	mattis Mattis D C,  Bardeen J , 1958 Physical Review 111 412
	
	
	kittel Tinkham M, 2004 Introduction to Superconductivity 
		
  impu Peruzzi A, Gottardi E, Peroni I, Ponti G and Ventura G 1999 Nucl. Phys. B 78 576
	
	match Kaplan S B 1979 J. Low Temp. Phys. 37 343
	
	 
	
	
	  doyle Doyle S, Naylon J, Mauskopf P D, Porch A and Dunscombe C 2008 J. Low Temp. Phys.
  151 530 
	
	casali Casali N et al 2016 J. Low Temp. Phys.  184 142
	
	Catalano15 Catalano A et al 2015 Astron. Astrophys 580 A15 

	usadel Usadel K D 1970 Phys. Rev. Lett. 25 570
	
  cooper Cooper L N 1961  Phys. Rev. Lett. 6 689

	lesueur Le Suer H, PhD Thesis, Universite Paris 6 (2007)
 
  strongin Strongin M, Kammerer O F and Paskin A 1965 Phys. Rev. Lett. 14 949
  
	zhao Zhao S, Goldie D J , Witington S , Thomas CN , 2018Supercond. Sci. Technol. 31 015007
  
	Gao Gao J 2006 NIMPA 559 585
	
  
  
  
  
  
  bourrionBourrion O , Bideaud A, Benoit A, Cruciani A, Macias-Perez J F, Monfardini A, Roesch M, Swenson L J and Vescovi C 2011 JINST 6 P06012
 
swensonjap Swenson L J, Day P K, Eom B H, Leduc H G, Llombart N, McKenney C M, Noroozian O and Zmuidzinas J 2013 J. Appl. Phys.  113 104501

pippo Goldie D J, Withington S 2013 Supercond. Sci. Technol.  26 1

  cruciani2 Cruciani A et al 2016 J. Low Temp. Phys.
  184 859
 
gatti Gatti E, Manfredi P F 1986 Riv. Nuovo Cimento  9 1
 
  
  
  
	
 
  
 

 gao_TLS Gao J, Daal M, Vajonakis A, Kumar S, Zmuidzinas J, Sadoulet B, Mazin B A, Day Pk, Leduc H G 2008 Appl. Phys. Lett. 92 152505
	



