Temperature quenching in LAB based liquid scintillator

The effect of temperature changes on the light output of LAB based liquid scintillator is investigated in a range from -5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-5$$\end{document} to 30∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$30\,^{\circ }$$\end{document}C with α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-particles and electrons in a small scale setup. Two PMTs observe the scintillator liquid inside a cylindrically shaped aluminum cuvette that is heated or cooled and the temperature dependent PMT sensitivity is monitored and corrected. The α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-emitting isotopes in dissolved radon gas and in natural Samarium (bound to a LAB solution) excite the liquid scintillator mixtures and changes in light output with temperature variation are observed by fitting light output spectra. Furthermore, also changes in light output by compton electrons, which are generated from external calibration γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}-ray sources, is analysed with varying temperature. Assuming a linear behaviour, a combined negative temperature coefficient of (-0.29±0.01)%/∘C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${(-0.29 \pm 0.01)}{\,\%/^{\circ }}\hbox {C}$$\end{document} is found. Considering hints for a particle type dependency, electrons show (-0.17±0.02)%/∘C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${(-0.17 \pm 0.02)}{\,\%/^{\circ }}\hbox {C}$$\end{document}, whereas the temperature dependency seems stronger for α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-particles, with (-0.35±0.03)%/∘C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${(-0.35 \pm 0.03)}{\,\%/^{\circ }}\hbox {C}$$\end{document}. Due to a high sampling rate, a pulse shape analysis can be performed and shows an enhanced slow decay component at lower temperatures, pointing to reduced non-radiative triplet state de-excitations.


Introduction
In neutrino-and astroparticle physics, liquid scintillator detectors have been used in many milestone experiments, and upcoming large-scale liquid scintillator experiments also indicate the importance of this detector material. Linear alkyl benzene (LAB 1 ) as solvent together with PPO 2 as primary solute are the basis for a modern organic liquid scintillator, a e-mail: arnd.soerensen@gmail.com 1 C 6 H 5 C n H 2n+1 , n = 10 − 13. 2 2,5-diphenyloxazole. which is already in use by the Daya-Bay [1] and RENO [2] experiments, will soon be used by SNO+ [3] and is also foreseen for future neutrino detectors like JUNO [4].
While most properties of LAB based liquid scintillators are well studied, the influence of temperature on the light generation is poorly known and only recently interest in this matter emerged again. The SNO+ experiment will be operated at ≈ 12 • C due to the cooled water shielding surrounding the liquid scintillator vessel [3], and JUNO investigated the eventually positive aspects of a significant cooling of the active detector medium on the energy resolution [5].
Previous works On the subject of temperature effects in liquid scintillators, real data and good studies are scarce, especially when focusing on changes in light yield. The very first investigations were carried out by Seliger and Ziegler [6] in the mid 1950s. Unfortunately the temperature of the whole measurement setup was changed, including the PMTs. Therefore that work was critically discussed and commonly disqualified by Laustriat, Coche [7,8] and Birks [9].
The former two did investigate temperature effects themselves with three different toluene solutions. While maintaining the PMTs at constant temperature and exciting the mixture via UV light, negative temperature coefficients in light output for PPO at about 20 % and about 10 % for PBD are reported, both at −10 • C. When 137 Cs was used as a radiation source, the effect is much less expressed with less than 5 % for both solvents, again at −10 • C. This study is difficult to interpret quantitatively since the setup is not described, it is not explained how temperature stability on the PMTs is achieved or how this was controlled, and no uncertainties are reported. Therefore, unfortunately, this can only be seen as an indication that scintillation light output might be quenched at higher temperatures and is enhanced when the liquid is beeing cooled down.
More than two decades later Homma, Murase and Sonehara [10][11][12] revisited this topic and carried out some measurements with αand β-radiation sources dissolved in five different liquid scintillator solvents: o-xylene, m-xylene, ethylbenzene, benzene and cumene. The single fluor solute in all of them was 4 g/l PPO. The radiation sources were 241 Am and 131m Xe, which is, besides γ -rays and X-rays, emitting mono-energetic conversion electrons of E β = 0.164 MeV. Even though pulse-height spectra are shown, again the data is not analysed quantitatively, no uncertainties are given and only a very rough description of the setup is included. It is claimed that the PMT temperature is somehow stabilised for all measurements. In a very coarse 'by hand' approximation of presented energy spectra, a temperature coefficient of about −0.5 %/ • C might be extracted. Nevertheless for all scintillator mixtures and both radiation types a significant increase in light output at lower temperatures is reported.
In 1999 Buontempo et al. [13] investigated a variety of combinations of four solvents and four solutes. Irradiated with a 90 Sr β-source, when a neutral gas was present, a light yield increase of ≈ 10 % was observed at −10 • C. Unfortunately again, temperature effects on the PMT were ignored and a setup scheme shows that the cooled and heated cuvette was directly mounted onto the PMT, thus having direct thermal contact.
Finally, the so far last and best documented work was published by Xia et al. [5]. This group built a detector setup much alike the one developed for this work. A 137 Cs γ -ray source was used in a Compton scattering setup that could be heated and cooled. To account for PMT sensitivity changes, an external LED light source was used to generate pulsed reference light, which was led via an optical fibre onto the PMT. With this technique any 'hysteresis' effects were excluded and changes in the PMT sensitivity were measured and corrected. Unfortunately the exact composition of the tested liquids is not reported, and uncertainties for the quantitative results are not presented. With a linear approximation of the shown LAB data, a temperature coefficient of ≈ −0.30%/ • C can be derived from this study.
Even though most of the previous studies have significant conceptual problems or report only coarse results, all show consistently that organic liquid scintillator solutions yield more light at lower temperatures.
This work The influence of temperature changes in the range of −5 to 30 • C on the light output of LAB based liquid scintillator mixtures are studied, for α-particle and electron signals. A small detector setup was developed (Sect. 2) and special care is taken to control eventual changes of the PMT sensitivity due to varying temperature (Sect. 2.1). An energy calibration for all used LAB mixtures has been carried out by irradiating the setup with γ -ray emitters and fitting the resulting Compton-edge structure in the energy spectrum (1) PMT mounts 3D-printed with ABS (acrylonitrile butadiene styrene) material, (2) UV transparent acrylic light guides, (3) quartz glass windows, (4) cavity of ≈ 110 ml to hold the liquid scintillator, (5) buffer volume & fill and empty port, (6) four pipes, flushed with glycol-water mixture for cooling and heating, (7) optical fibre feed through to point external LED reference light directly onto the PMT cathodes. The aluminium caps holding the temperature sensors are shown in light grey in the upper drawing (Sect. 2.2). As this has also been done at multiple temperatures, the changes of the light output for the Compton scattered electrons can be extracted (Sect. 2.3). After evaluating the detector response with an unloaded liquid scintillator mixture (LAB + 5 g/l PPO + 20 mg/l bis-MSB 3 ) at different temperatures, measurements with a samarium loaded (Sect. 3.1) and a radon loaded (Sect. 3.2) LAB mixtures were carried out in the given temperature range. Beside the studies of the energy spectra, a pulse shape analysis (Sect. 3.3) yields important insights. All liquid scintillator mixtures were bubbled for 20 min with pure nitrogen before they were filled into the detector, which itself was flushed with nitrogen.

Experimental setup
The detector itself, shown in Fig. 1, consists of an aluminium block with a cylindrical cavity to hold the liquid scintillator, and two opposing mounted Hamamatsu R2059 PMTs to detect the scintillation light on either side of the vessel. These are thermally decoupled from the liquid by quartz glass windows and 3 cm long acrylic light guides, since the PMT sensitivity is known to significantly change with tempera-  [7,14,15]. All light guiding parts were optically coupled with silicone grease. The temperature of the active medium is changed by flushing a heated or cooled glycol-water mixture through four pipes in the aluminium block, resulting in temperatures between −5 and 30 • C inside the liquid scintillator. These temperature changes are detected by two independent temperature sensors (DS18BS20), with an accuracy of 0.5 and 0.044 • C precision.
As the thermal decoupling of the PMTs cannot be perfect, PMT sensitivity changes are still expected and addressed by a monitoring/correction technique, see Sect. 2.1. To make sure that all PMT gain changes are noticed and can be corrected, an independent light source is used to monitor the PMT response. A schematic overview of the setup and readout electronics is given in Fig. 2.
To achieve PMT coincidence detection and also to include the signals from the LED, the electronic signals from the PMTs are split and simultaneously lead to a NIM 4 electronics crate for coincidence-and counting logic, and to a 2 GS/s, 10bit FADC 5 (Acquiris DC-282).
In order to extract a reliable information about the detected amount of light from the recorded pulse traces, the integral of the pulse is used. Furthermore, considering isotropically produced light and homogeneous reflectivity on the vessel walls, the use of the geometric mean of the pulse area values eliminates the dependency on the position of the light generation, resulting in a superior energy resolution for this type of detector. In order for this technique to work, the gain of 4 Nuclear instrumentation module. 5 Flash analog-to-digital converter. Comparison of averaged pulse shapes from ≈ 1 million LED pulser events as they appear on one PMT (red dotted line) and the photodiode (blue solid line). With a rise time shorter than 175 ps, the photodiode is able to resolve small 'ripples' of the LED signal, caused by the LED driver. The photodiode signal is linearly amplified the two PMTs are matched by placing a gamma-ray source at the centre (± 0.5 mm) and adjusting the PMT high voltage.

PMT sensitivity control
With a thermal conductivity of λ ≈ 0.2W m −1 K −1 , the acrylic light guides insulate the PMTs well, but not perfectly. Furthermore, the box in which the detector is placed in a nitrogen atmosphere, reaches a thermal equilibrium with the cuvette only after some hours. That is why temperature changes of up to ≈ 10 • C on the PMTs are expected. The manufacturer for the PMTs used in this setup, states a flat 'temperature coefficient' for the anode sensitivity of ≈ −0.4 % / • C [14], which could be confirmed in this work.
As the temperature changes occur on long timescales in an evacuated phototube, they are difficult to quantify for individual and not identical phototubes, and a dedicated control of the PMT sensitivity is necessary. In order to achieve such a PMT sensitivity control, a custom built pulser drives a LED (428 nm) with a variable frequency to generate very short light pulses of ≈ 10 ns (see Fig. 3) that are as similar as possible to the light emitted from the liquid scintillator. This light is then split and lead simultaneously onto the two PMTs and a thermally fully decoupled silicon pin-photodiode (Alphalas UPD-200-PS) outside the N 2 -box, refer to Fig. 2. This independent signal serves as normalisation and control that the amount of emitted light from the LED is constant.
Proof of concept Extensive tests have been undertaken to show that the light reference by the LED pulser is stable, reliable and that it can be actually used to correct for PMT sensitivity changes by ambient conditions. Long time measurements of up to 24 h were made to guarantee thermal equi- librium of all parts and gain insights into time stability. These datasets were divided into 1-h-'blocks' to observe changes on this time base. A pulser rate of 20 Hz guarantees good statistics without interfering with the signals. As the background rates in the basic LAB solution (LAB + 5 g/l PPO + 20 mg/l bis-MSB) are low with about 40 Hz (integral over the full energy spectrum), random coincidences with the pulser are negligible ∼ 10 −5 Hz, considering a coincidence window of 30 ns. The resulting Gaussian shaped peaks in the charge spectra were fitted with a standard Gaussian and the mean and the fit uncertainty on the mean were extracted. Graphs for this analysis with one PMT and the photodiode are shown in Figs. 4 and 5, respectively. In this case, after 24 h at room temperature, the cooling device was used to lower the temperature of the liquid scintillator to 5 • C (left part of the diagrams). After ≈ 4 h, the final temperature is reached, and kept stable by constant cooling. The influence on the PMT sensitivity of the ambient temperature change is clearly visible. After about 5 h, a thermal equilibrium is reached and the differences between the hourly data points turn into a random scatter of the order O(0.1-0.2%). Note that these data points are hour averages and the deviation from the overall mean.
From this observation it was concluded that a time of at least 6 h needs to be in between the set of a new temperature on the cooling device and actual data taking. In practise, the cool down and warm up periods were during the night and the time span before measuring at least 18 h. Furthermore, it could be shown that the LED pulser and / or the photodiode read-out is long time stable to ± 0.2 %. These fluctuations are due to fluctuations of the intensity of the LED pulser source. This allows to measure changes in the amount of light hitting the PMTs of the order O(1 %).
Finally, the application of a correction to real data, that was taken at different temperatures, needed to be shown. Since Comparison of charge spectra from the LED pulser, obtained at 5 • C (blue) and 24 • C (black). The red one is generated from the same data as the one in blue, but with the sensitivity correction applied, event by event. The solid black line shows the Gaussian fit to the low temperature data the pulser trigger signal is also recorded every time the LED is lit, these events can be 'cut' from the data-set and analysed separately. As the area of the pulses are Gaussian distributed around their mean, a standard Gaussian function is used to fit the data and to extract the mean value and fit uncertainty, see Fig. 6. Once this is done for both measurements, the correction factor can be calculated simply as the ratio of the two LED pulser mean values. Since the signal from the photodiode is known to be very stable, any shift from this source can then be subtracted from the correction quotient, as this would indicate a slight change in the light that hits the PMTs. This correction factor can then be used to scale the pulse integral on an event by event basis. As a test, two measurements with a 22 Na calibration source were carried out, one at room temperature (24.20 ± 0.04) • C and one at (4.80 ± 0.04) • C.  22 Na charge spectra. The Compton edge at higher energy is clearly visible. The one at the lower energy is less pronounced and merged with the back-scatter feature. In black, the spectrum obtained at (23.65 ± 0.04) • C, and in blue the one from the measurement at (4.80 ± 0.04) • C. The red spectrum is generated from the same data as the blue one, but with the sensitivity correction applied, event by event. The inset shows a zoom at the upper Compton edge A comparison of the spectra from the 22 Na calibration source can be seen in Fig. 7. A closer look at Fig. 7 reveals a small change of the Compton edge towards higher charge values for the cooled scintillator. With the application of the gain correction, this effect is reduced and the corrected spectra 'moved' towards the room-temperature position. This is the expected behaviour of the sensitivity correction via a separate light source, as it is applied in this setup.
An important test is the effect of temperature changes on the pulse shapes of the PMTs, in this case illuminated by the LED. While the photodiode is perfectly thermally decoupled from the cooled setup, the PMTs undergo temperature changes. With its 2GS/s, the FADC has the capability to record the pulse shape of the PMT voltage pulses. Figure 8 shows an overlay of average pulse shapes of LED pulser events on PMT-1, at five temperatures, normalised to the pulse height at room-temperature. The arithmetic average of the voltage value is calculated at each sample position (0.5 ns bin), as well as the variance of the voltage values. Since the variance is very small and the amount of pulses large (≈ 150 k), the uncertainties are too small to be visible here. It can be noted that, clearly the temperature variation has no effect on the shape of the LED pulses on the PMTs, and therefore that the PMTs do not produce different pulse shapes at different temperatures.

Energy calibration
The energy resolution and calibration is obtained by fitting the upper edge of the Compton continuum in the charge spectra of external gamma-ray sources. Firstly, the ideal Compton spectrum is calculated for the individual isotope energies, making use of the analytic crosssection calculation by Klein and Nishina [16]. This ideal spectrum is generated for each individual gamma-line, and those are then combined into a sum-spectrum with the exact branching ratios for the individual decays [17], see the solid black line spectrum in Fig. 9. The used gamma-ray calibration sources ( 137 Cs, 22 Na, 60 Co, 207 Bi, 88 Y) cover the energy range from 477 keV to 1612 keV with their respective Compton edge energies. Secondly, the actually sharp Compton edge is 'smeared' by the energy resolution of the detector and needs to be taken into account for a fit. For every gamma-line, the Compton edge needs to be fitted locally, since the energy resolution is energy dependent itself. Finally, from each edge-fit the position and energy resolution can be extracted. After a rough 'by hand' approximation, the extracted parameters for the position and resolution are used to start a full likelihood fit, using the Minuit routines of the CERN ROOT software packages. The ideal spectrum convoluted with an energy dependent energy resolution is drawn with a solid red line in Fig. 9, where the following Eq. 1 is used as analytical function for the energy resolution: with p 1 as parameter to describe the finite resolution due to electronic noise, p 2 regarding statistical effects of photon and photo-electron generation and p 3 the influence of the position of the light generation and attenuation. For all three LAB mixtures used in this work, the energy calibration and energy dependent energy resolution obtained from these fits are shown in Figs. 10 and 11, respectively. A detailed description of the used loading technique and properties of these mixtures can be found at the beginning of the Sects. 3.1 and 3.2.
As expected, the LAB mixture with metal loading emits much less light (see Fig. 10), or more precisely, much less light is detected, due to reduced transparency of the scintillator. This effect can also be seen in UV-transmission measurements. It is also the main reason why the energy resolution is worse for samarium loaded LAB. Another factor increasing mainly the uncertainty of the energy resolution for lower energies in this scintillator mixture, is the alpha peak from the 147 Sm, disturbing the Compton edge fitting.
Moreover it can be seen that the energy resolution in the 'pure' scintillator mixture (LAB + 5 g/l PPO + 20 mg/l bis-MSB) is not affected when 222 Rn is added to the liquid only by diffusion, see green and blue data in Fig. 11. Details of the loading process will be discussed in Sect. 3.2. The reason why the green data points are below the blue ones in Fig. 10 is Energy dependent energy resolution for three mixtures: 'pure' LAB + PPO + bis-MSB (blue), samarium loaded LAB (red) and radon loaded LAB (green), at room-temperature (24.5 ± 0.5) • C. The fit function (see Eq. 1) is shown with its 1σ uncertainty band, obtained from the fit the fact that for the radon loaded LAB measurement a lower high voltage on the PMTs had to be used, since they were different PMTs (but same type) and the activity of the radon was high.

Light yield calibration at different temperatures
Previous studies by Homma et al. [10][11][12], Buontempo et al. [13] and especially Xia et al. [5] have shown or at least indicated that an increase in electron light yield with lower temperatures has to be expected. Thus, the calibration technique described in the previous section is applied at different temperatures. This is especially important to identify a possible difference in the temperature quenching between α-particles and electrons.
Four energies from three different gamma-ray emitting isotopes ( 22 Na, 207 Bi and 137 Cs) were used to determine the position of their Compton edges at five temperatures (see Fig. 12) in the 'pure' scintillator mixture (LAB +5 g/l PPO+ 20 mg/l bis-MSB). Combining the individual results for each energy, Fig. 12 shows the change in charge versus temperature. The slope of the linear fit is (−0.173 ± 0.018)%/ • C, indicating a change in light yield in the same order of magnitude as the results reported in [5]. There, only diagrams with data points and splines are presented and no linear approximation is performed. Therefore only a very rough 'by hand' estimate of ≈ −0.32%/ • C can be used for comparison with the data obtained in this work.
A comparison of the linear fit slopes for the individual isotopes and their energies will be presented in final discussion of all results in Sect. 4.

Samarium loaded LAB
At the Brookhaven National Laboratory (BNL), natural samarium was loaded into LAB based liquid scintillators, using the surfactant TMHA 6 in a single-stage solventsolvent extraction procedure, which can be used to form organometallic compounds that can be dissolved in an organic liquid.
Natural samarium contains (15.01±0.03) % of the isotope 147 Sm, which is α-instable with a half-life of T 1/2 = (1.070± 0.009) × 10 11 y and E α = 2.248 MeV [18,19]. The mixture used for the measurements in this work contained about 2 % natural samarium, resulting in a 'constant' activity of about 200 Bq in the 100 ml volume of the detector used, from this decay. Although natural samarium also contains 11.24 % of the α-instable isotope 148 Sm, with its half-life of T 1/2 = 7 × 10 15 y [17], its activity of ≈ 2 × 10 −3 Bq is negligible in this case.
Analysis Figure 13 shows a calibrated light yield spectrum at room temperature. The fit function is a Gaussian, though as the samarium peak showed a high energy tail, other shapes were tried too (Landau-and skewed normal distribution), neither did fit well nor improved the goodness of fit significantly. Thus the model independent, conservative Gaussian shape was chosen, especially as no physics motivation for any other shape was found. The overall activity in the setup was too low for pile-up events, and no database entry of any X-ray emission for 147 Sm was found. 6 3,5,5-trimethyl-hexenoic acid.  Therefore, for all nat Sm-LAB temperatures the following fit routine was applied: First a simple Gaussian fit with a wide range is performed, to roughly identify the peak position and width. In a second step, the fit range is confined to 1σ around the mean of the first fit. The second (actual) fit is in this case restricted to 1σ to avoid any distortion due to the high energy tail. This is incorporated in the final analysis and evaluation by a 'fit range uncertainty', which is calculated as the variance of the fit mean values of Gaussian approximations with 1σ , 2σ and 3σ fit range. For the room-temperature measurements a list of the individual uncertainties considered is given in Table 1.
All temperature varied measurements have been performed at least twice, except the one for 30 • C which has only been done once. The sequence of the measurement was chosen for practical reasons to 24 and then repeated without the two 24 • C steps in between. As this does not allow to check for any systematic effects of the T changes in the detector system (e.g. hysteresis effects), a more systematic approach for the temperature sequence was chosen in the radon loaded LAB studies, presented in the next Sect. 3.2.
As room-temperature is an arbitrary but natural scale for liquid scintillator measurements, all changes observed with temperature are normalised to the room-temperature value.  Figure 14 gives an impression of how the temperature change affects the 147 Sm peak position, after the PMT sensitivity correction is applied to correct for any ambient temperature effects on the PMTs. In Fig. 14 the systematic shift towards a higher light yield with decreasing temperature is clearly visible and a linear fit yields a slope of (−0.26 ± 0.02) %/ • C. The data for the 30 • C measurement seems to be shifted disproportionately towards lower charge values. The visual effect is increased by the observation that the peak is also more narrow. Also in the final analysis of the peak shift, the 30 • C data seems to deviate of the seemingly linear trend of the other data points, see Fig. 15.
In order to incorporate also the results from the radon loaded LAB measurements and pulse shape analysis, a full discussion of the results is given in Sect. 4.

Radon loaded LAB
In order to extend the study of the light yield temperature behaviour with other α-emitting isotopes, preferably at higher energies and with another loading technique, the decay chain of 222 Rn was chosen. This isotope can be obtained from 226 Ra, which is an 100 % α-emitter with a halflife of T 1/2 = 1600y, decaying solely to 222 Rn. As radon is a gas, the isotope escapes the porous ceramic Radium source and can be transferred via gas exchange. This radon gas is used in a nitrogen atmosphere to 'bubble' through a LAB mixture. For better comparison, the same 'pure' LAB + 5 g/l PPO + 20 mg/l bis-MSB solution (without samarium) was used. After circulating the nitrogen -radon atmosphere through the LAB mixture for about 72 h, an activity of about 200 kBq was introduced into the liquid scintillator. This high activity was beyond expectation and the setup had to rest for several weeks to reduce the activity to about 1 kBq. Even though the 222 Rn decay chain ends only at the stable 206 Pb, in the scope of this work already the isotope 210 Pb presents the end of the chain with a half-life of T 1/2 = 22.2y, leading to an activity too low to be visible. This short chain produces three α-particles and has two β − -decays, refer to Table 2.
The first α-decay of 222 Rn has the longest half-life of T 1/2 = 3.82d and thus, as the other isotopes decay with minutes and micro-seconds, the chain is very fast at secular equilibrium. The first α-energy of 5.49 MeV is close to the one from the following 218 Po with 6.00 MeV, and can therefore not be distinguished with this detector setup. Nevertheless, in a combined fit, it is possible to extract the individual light yields. The resulting 214 Pb decays via β − -emission to 214 Bi, with an end-point energy of 1.0 MeV. Subsequently 214 Bi changes to 214 Po, emitting an electron with the maximum energy of 3.3 MeV. Its very short half-life of T 1/2 = 164.3 µs proves useful for tagging this and the following α-decay to the final 210 Pb. This last α-decay has an energy of 7.69 MeV, which separates it well from the preceding α-particles.
This LAB mixture with dissolved radon was used to observe light output changes at five temperatures, in the following sequence: No systematic effect has been observed, repeated measurements scattered around their mean randomly, without a sign preference. A calibrated light yield spectrum is shown in Fig. 16. Analysis As can be seen in Fig. 16, the two β-decays present a significant background for the quenched α-peaks in the energy spectrum. Different strategies can be applied to still extract the relevant peak shifts with varying temperature. One option is to analyse the energy spectrum as it is and either subtract the β-background via fitting it first, or include it in a combined fit with the α-peaks. Another option is to make use of the short 164.3 µs half-life of 214 Po and apply a timecoincidence cut to separate 214 Bi -214 Po (Bi-Po coincidence) event pairs from the decay chain and thus the energy spectra. Even though this second option seems very promising, the first one yields smaller uncertainties and will be presented here. The second analysis can be found in the Appendix A and confirms the absolute numbers of the 'combined fit method' nicely. In a previous study by Krosigk et al. [20], a 1 l LAB + 2 g/l PPO sample was analysed, which had been placed inside the water-filled SNO detector in 2008. The liquid scintillator was also 'contaminated' with dissolved radon gas and the energy spectrum shows very similar features as are visible in Fig. 16. The background, mainly consisting of the two βdecays, was approximated with a 6th order polynomial and then subtracted. Krosigk et al. excluded the α-peak region from their background fit.
In the analysis of this work, the α-peaks and the βbackgrounds were fitted in one combined fit, visible in Fig. 16. The fit function is the sum of three Gaussian functions and a 6th order polynomial. In order for this fit to properly converge, some of the 15 parameters needed to be constrained, mainly those of the two close α-peaks. As they stem from subsequent decays with no other branches, they have the same activity. The strength parameter was therefore the same for both Gaussians, yet it could float freely. Also the meanparameters of the first two Gaussians are constrained, such that they can not take the same value. The mean-parameter [%]  Table 4 Temperature coefficients (slope values) and uncertainties for the three α-isotopes in the radon loaded LAB mixture (+5g/l PPO + 20mg/l bis-MSB), obtained via the linear fits in Fig. 17 Isotope of the first Gaussian has to be in a range below their sum-Gaussian peak mean and vice versa. With these rough constraints, the fits proved to be very stable and well converging throughout the analysis, also for different temperatures. The room-temperature results for the peak positions are listed in Table 3 and are well compatible with the values obtained by Krosigk et al. When analysing the full energy spectrum with the described fit at different temperatures, changes of the individual α-peak positions can be extracted. The resulting slopes are shown in Fig. 17 and the values are listed in Table 4.
The repeated temperatures vary less than 0.5 • C, and for the final analysis their average has been used. Of course, also the PMT sensitivity correction is applied before analysing the spectra.

Pulse shape analysis
Since the Acquiris FADC is able to take a voltage sample from the PMTs every 0.5 ns, it is possible to perform a pulse shape analysis (PSA) with this experimental setup. As the variations between individual pulse shapes are very large, averaged pulses will be used for comparisons and fits. For a given pulse ensemble, at each sample point, the arithmetic average of the voltage values is calculated, as well as the variance between these pulses. With a typical size of one million pulses, the uncertainty (variances) are very small, and mostly not visible in the following graphs.
To illustrate the analysis principle Fig. 18 shows an averaged pulse shape from the 147 Sm α-peak obtained from the samarium loaded LAB mixture, at room-temperature. The typical, very fast rise to the maximum is clearly visible, as well as a slower falling slope, followed by an even less steep slope, which can be attributed to a fast and slow decay component. The slight kink visible at ≈ 123 ns might be attributed to a change in decay constants. Although it was demonstrated that it is not originating from any PMT effects (see PMT pulse shapes in Sect. 2.1, Fig. 8), these deviations from a perfectly smooth curve can not be fully explained. In order to quantify the changes in the slopes, the following composition of two exponential functions was chosen, convoluted with a time response function with a Gaussian resolution and fitted to the average pulse shapes: with U -voltage, t -time, τ f ast , τ slow -decay constants, c 1 ,c 2 -constants and σ -resolution of the time response function. It was found that a two component exponential function fitted the data best, all attempts with a third component lead to fits with a vanishing third part. Also the rising edge was not included here, since it is the same for all temperatures (see e.g. Fig. 19). Once c 1 and c 2 are fixed to their fit value at roomtemperature, the approximations for the other temperatures reveal the change in the decay constants. These are e.g. for 147 Sm at room-temperature: τ f ast = (4.48 ± 0.01) ns and τ slow = (16.90 ± 0.01) ns for the fast and slow component, respectively. Figure 19 shows a comparison of averaged pulse shapes from the 147 Sm α-peak, at different temperatures.

Samarium loaded LAB
It is clearly visible that the shapes change between 115 ns and 140 ns. With lower temperatures, the falling slope of the pulses becomes less steep. In contrast, the rising edge is the same for all temperatures. The changes that these constants undergo is illustrated in Fig. 20. It is quite clear that the second, slow component is stronger expressed with lower temperatures, showing an increase of about 4 % over the range of 35 • C. Radon loaded LAB For the PSA of the radon loaded LAB measurements, it is sensible to use the data-sets with the isolated βand α-events. This is done by applying the same time and energy cuts as in the alternative analysis of the peak shifts, described in the Appendix A. When analysing the averaged pulse shapes from the 214 Bi β-decays, they seemingly resemble the prolonged pulses that can be seen from the α activity in samarium loaded LAB (Fig. 19). Here again, with lower temperatures the pulse get longer, while with rising temperature the slow decay component is to be suppressed (Fig. 21). To quantify this behaviour and compare with the results from the α-decays, the same fit as for the samarium loaded LAB is applied (Fig. 18, Eq. 2). The results for the two decay constants from the fitting are shown in Fig. 22. As it has been observed for the α-decay of 147 Sm, the fast decay component experiences only minor changes with temperature variation of about ≈ 1 % over a range of 35 • C. On the other hand, even though the uncertainties are larger, the slow decay part changes ≈ 5 % over the same range. This is significantly more and close to the observation from 147 Sm, even though here it is electrons instead of α-particles.
This similar behaviour can again be confirmed with yet another α-source, the well isolated 214 Po decay in the radon loaded LAB. Again, applying the energy-and time-cuts described in the Appendix A ensures the very good separation of the polonium decay from unwanted backgrounds. Again, while the fast decay component exhibits only a small change with temperature, the slow decay becomes ≈ 4% stronger at 35 • C lower temperature. This is very similar to changes observed from 147 Sm, and less than that of the 214 Bi β-decays (see Fig. 22).

Discussion and conclusions
If the results from all energies are combined and a linear approximation is assumed, a temperature coefficient of (−0.17 ± 0.02) %/ • C for electrons is obtained in this work via energy calibration at different temperatures in 'pure' LAB based liquid scintillator (Sect. 2.3). This is in the same order of magnitude as the observation by Xia et al. [5], showing diagrams which correspond to coefficients of ≈ −0.30%/ • C. Their setup is well designed and thus their result seems trustworthy and find more or less a confirmation with this independent study. Although there are more studies to mention here, namely those by Seliger and Ziegler [6], Laustriat and Coche [8], Homma et al. [12] and Buontempo [13], all of these have significant problems either in their measurement concept or present their data in way that does not allow for a direct comparison. Nevertheless, all of those publications (except [6]) indicate and show a significant increase in electron generated light output at lower temperatures.
For α-particles previous measurement data is even more scarce, only Homma et al. [10,12] present pulse-height spectra from 241 Am at different temperatures. Unfortunately they do not analyse these spectra quantitatively but show only very coarse diagrams. A very rough calculation from a 'by hand' analysis of the presented spectra, using only two data points for 241 Am in a PPO-benzene mixture, results in ≈ −0.5%/ • C. In this work, the combined average temperature coefficient for all energies, including the values for 222 Rn, 218 Po, 214 Po and 147 Sm, is which is in the same order of magnitude again and confirms the overall trend. It must be pointed out that it is difficult to combine the coefficients originating from two scintillator types. While the diffusion of radon atoms does not involve  Fig. 23 Temperature coefficients for four gamma-ray emitting isotopes (blue, filled data points) and four α-emitters (green hollow points, pink for samarium) in LAB based liquid scintillator. Model testing has been done by fitting a global constant (red dash-dotted line) two separate constants (red solid lines) and a linear function (red dashed line) any significant change in the chemistry of the liquid scintillator, the metal loading technique does. Therefore, a comparison of the independent temperature coefficients is sensible and can be done by analysing these values energy dependent. Figure 23 illustrates the combined data for electrons (generated by gamma-rays in 'pure' scintillator, Sect. 2.3) in blue and for α-particles in radon-and samarium loaded LAB mixtures (green and pink, hollow symbols). The first observation is the already mentioned overall trend of slightly negative temperature coefficients of less than one percent per degree Celsius. Secondly, it is clearly visible that the samarium value (pink data point) deviates from the other α-isotope values, measured in the radon loaded scintillator, and is connected to a much larger uncertainty. This can partly be explained by its position at lower energies, but might also have its origin in the fact that it is measured in another scintillator mixture, which is already expressed in the worse optical properties and worse energy resolution.
Concerning a possible energy dependence, three very simple models have been tested. Fitting a global constant describes the data very badly (p-value < 10 −10 ) but excludes a zero temperature coefficient already by more than 30σ . The fit value for a global constant is (− 0.29 ± 0.01)%/ • C. A better approximation, but still very bad, is possible by assuming a linear function, which yields a p-value of ∼ 10 −6 . Out of the chosen models, a fit with separate constants for electrons (−0.17 ± 0.01) %/ • C and alphas (−0.36 ± 0.01) %/ • C yields the best result with a probability value of 6 × 10 −4 . One reason for the overall bad results in approximating an analytical model is the position and large uncertainty of the 147 Sm data point. With these few data points, a final con-clusion seems difficult and none of the three simple models can describe the data well. Judging purely from the numbers, the two constants model has a slight preference. As it is very difficult to measure α sources at even lower energies (due to the ten-fold quenching in energy), the only realistic option in clarifying this would be measuring higher energy electrons, in the range above 3 MeV. This might be achievable with high energy gamma rays or the measurement of the energy spectrum of the electrons from muon-decay, reaching up to 50 MeV at their upper, sharp end. The data point estimated from Xia et al. [5] for 137 Cs of ≈ −0.30%/ • C would push the first blue data point in Fig. 23 down quite a bit and more towards the regime of the alphas. Beyond these most simple models, it could also be possible that a linear energy dependency for electrons is realised, whereas α-particles might experience a larger shift with a temperature dependent Birks' parameter.
Pulse shape analysis Useful insights for a physics interpretation are also provided by the comparison of the pulse shape analysis results. For this work, the temperature behaviour of averaged pulse shapes from 147 Sm, 214 Po and 214 Bi have been analysed. All of these are averages of a very large amount of shapes, on the order of 1 million each. They have been fitted by a two part exponential function with two separate decay constants (Eq. 2), which turned out to be the only reasonable way to approximate these shapes analytically. All three isotopes show a significant increase in the strength of the later pulse shape with lower temperatures, such that the decay constant characterising the 'slow' decay component increases. An overview of some of these values is given in Table 5. When comparing the diagrams for the relative change of the constants (Figs. 20, 22) the similarities are striking. Especially the trace of the data points for the slow component seem to follow a non-linear curve, maybe quadratic or exponential. Nevertheless, the simplest linear approximation gives the following results for the slope of the slow decay constant versus temperature: 147 Sm (−0.09 ± 0.003) %/ • C, 214 Po (−0.10 ± 0.005) %/ • C and 214 Bi (−0.14 ± 0.007) %/ • C.
As the two α-particles might be actually the same, the bismuth electrons seem to show an even larger increase in the slow decay time with lower temperatures. If this particle dependence is true, a lower scintillator temperature would increase pulse-shape discrimination capabilities, although an improvement on an event-by-event basis needs to be demonstrated since these differences were found in large number averages.
A final physics interpretation needs to take a lot of different information into account. Firstly an overall increase in scintillation light with lower temperatures has been reported Appendix: Alternative analysis: radon loaded LAB This appendix section shall illustrate an alternative approach to analysing the data obtained with the radon loaded LAB mixture (see Sect. 3.2). In the relevant part of the 222 Rn decay chain, the 214 Po decay has very short half-life of 164.3 µs, compared to the other isotopes. This is an opportunity for a time coincidence analysis, which will be used as a first cut in this alternative analysis. A closer look at the times in between two events is insightful, this spectrum from one measurement is shown in Fig. 24. Theoretically, the timedifference distribution between 214 Bi and the subsequent 214 Po event should be an exponential with a decay constant of 4219s −1 . When looking for this structure in the time -difference spectrum in Fig. 24, it can be found at the left end of the spectrum in the interval [0 µs,500 µs]. Events with the less steep slope in centre belong to the other decays with smaller decay constants and therefore longer half-lives. It is clearly visible that the spectrum is continuous and that the overall activity is quite high, which is not beneficial for this kind of analysis. This high activity of about 1 kBq at the beginning of the measurements was chosen due to practical reasons. For each temperature point some hours are needed to perform the measurements and changing the temperature of the liquid scintillator and the detector, so that about one temperature per day could be investigated. Thus, making 10 measurements took about 2 weeks, which is just in the range of the 222 Rn half-life of T 1/2 = 3.82 d.
This relatively high activity mainly disturbs the time coincidence between the polonium and polonium alpha, which expresses itself in a 'contamination' with background events in the energy spectra after a time cut is applied. For this reason  Fig. 25 Charge spectra from radon-loaded LAB, populated with events that are falling into different time windows, ranging from 60 µs to 130 µs. The 214 Po α-peak is clearly pronounced, but also background events from the double-α feature and the β-spectra are still visible to the left. The 70 µs cut has the best peak to background ratio a heuristic approach for the length of the coincidence time window is chosen. In Fig. 25, energy spectra of the second event (supposedly from 214 Po) are plotted, dependent on different coincidence time windows. The 70 µs spectrum shows the best peak to background ratio, and was hereafter chosen as time cut to isolate the 214 Bi-214 Po coincidence. Although the 214 Po α-peak is nicely pronounced and isolated with this first cut, still some events from the β-spectra and the other α events are 'leaking' into these spectra.
To further improve this selection, a second cut is applied, this time on the energy of either the first or the second event in the coincidence, to clean up the spectra of the other event, respectively.
To improve the 214 Bi β-spectra, the following selection of events is performed: 1. events with a time difference of 70 µs 2. first event of those pairs is taken, if the second event lies within the [817,1170] keV energy window (the 214 Po αpeak) The effect of this is demonstrated in Fig. 26. Although the actual signal sacrifice is large, the improvement on the shape and the strong suppression of disturbing α-events make this cut still useful. For β-spectra it is difficult to quantify any small change of light output with temperature change, since they do not represent a steep, peak-like structure. A visual inspection reveals no visible change in the β-spectra with changing temperature. Nevertheless, an analysis of the pulse shapes indicates also a temperature dependence for the decay constants of the betas, see Sect. 3.3.
A similar event selection can be performed to enhance the 214 Po spectra: 1. events with a time difference of 70 µs Comparison of charge spectra obtained from radon loaded LAB, when applying only time cuts (blue) and constraining also the energy of the following event to the energy of 214 Po (red). The peak-like features in the blue spectrum are effectively suppressed when applying the energy cut. Since signal sacrifice is large, the red spectrum is normalised to the area of the blue, for better comparison Again, despite a large signal loss, the effect of reducing background events that disturb the shape and therefore also the position of the α-peak, is worthwhile. Beside the greatly improved signal to background ratio (compare Fig. 25 vs. Fig. 27), also the fit range uncertainty decreased from 1.7 % to 0.3 %. This is already referring to the final analysis of the peak positions, that again is performed in a two step fit procedure. As it has been done for the analysis of the samarium loaded LAB, here again a first standard Gaussian fit determines the rough peak position. Then, in a second step, the mean and sigma from the first fit is extracted to perform the actual Gaussian fit in the 2σ range around this mean. Table 6 gives an exemplary overview of the uncertainty val-   Fig. 27), in blue. And from fitting the double α-feature of 222 Rn and 218 Po, in red. Each temperature point is the average of the repeated individual measurements and these are normalised to the room-temperature value at (24.5 ± 0.04) • C. The shaded areas mark the 1σ uncertainty band from the linear fit ues for the analysis of the first room-temperature measurement.
An illustration of the effect of the temperature changes on the α-peak position can be found in Fig. 27, the PMT sensitivity correction is already applied. Also a graphical representation of the fit results versus temperature is given in Fig. 28. In a last variation of the analysis, a fit to the double α-feature of 222 Rn and 218 Po is performed, mainly for comparison and cross-check. This time, events that fall into the 70 µs coincidence time window of the 214 Bi-214 Po time window, are excluded from the analysis. This suppresses the 214 Po αpeak very effectively as well as the background contribution of the 214 Bi β-spectrum. Thus the shape and the position of the double α-peak is less influenced by the tail of the higher energy αand β-decays. This time, only one single Gaussian is fitted to the double α-feature, though in the same manner as before. That means a first fit to roughly identify peak position and width and then a second more refined approximation. Since there is still one β-spectrum in the background, the uncertainties are much larger for this rough and simple analysis. Nevertheless the results, also shown in Fig. 28, provide a cross-check and comparison to the outcome of the other analyses and are agreeing well within the error margins.