Measurement of the proton light response of various LAB based scintillators and its implication for supernova neutrino detection via neutrino–proton scattering

The proton light output function in electron-equivalent energy of various scintillators based on linear alkylbenzene (LAB) has been measured in the energy range from 1 MeV to 17.15 MeV for the first time. The measurement was performed at the Physikalisch–Technische Bundesanstalt (PTB) using a neutron beam with continuous energy distribution. The proton light output data is extracted from proton recoil spectra originating from neutron–proton scattering in the scintillator. The functional behavior of the proton light output is described successfully by Birks’ law with a Birks constant kB between (0.0094±0.0002) cm MeV−1 and (0.0098±0.0003) cm MeV−1 for the different LAB solutions. The constant C, parameterizing the quadratic term in the generalized Birks law, is consistent with zero for all investigated scintillators with an upper limit (95 % CL) of about 10−7 cm2 MeV−2. The resulting quenching factors are especially important for future planned supernova neutrino detection based on the elastic scattering of neutrinos on protons. The impact of proton quenching on the supernova event yield from neutrino–proton scattering is discussed.


Introduction
In the last decades enormous progress has occurred in the field of neutrino physics. A new field of astrophysics has been opened with the observation of solar and supernova neutrinos [1][2][3][4]. First indications of ultra-high energy neutrino detection exist [5]. From all sources the appearance of supernova neutrinos is the least predictable but a lot can be learned if such an event occurs in the Milky Way. In the final stages of a core-collapse of massive stars nearly all a e-mail: belina.von_krosigk@tu-dresden.de the binding energy of the star is expected to be released in the form of neutrinos and antineutrinos of all flavors. The expected spectrum consists basically of two components: The deleptonization burst consisting of ν e released within a few ms and the emission of all kinds of neutrinos from the Helmholtz-Kelvin cooling phase of the protoneutron star lasting several seconds. For recent reviews see [6][7][8].
The intense neutrino pulse can be observed on the Earth by neutrino detectors with a sufficiently low energy threshold of about 1 MeV. For the SN 1987A this was achieved by Kamiokande II, IMB and the Baksan Scintillator Telescope [2][3][4]. The signal is expected to be dominated byν e because of the large cross section for the respective reaction. However, it is of great interest to measure the flux of ν μ , ν τ , ν μ andν τ , collectively called ν x , and ν e as well. The detailed shape of the overall spectrum depends strongly on the model simulation. Recent progress has been made in using 3-dimensional simulations taking into account large scale convection. Furthermore magnetic fields, general relativity and other effects were included resulting in different spectral details. As the average neutrino energy is below about 30 MeV, charged current interactions are only possible for electron type neutrinos, leaving only the far more challenging neutral current reactions as possible experimental detection channel for supernova ν x . In [9,10], neutrino-proton elastic scattering, ν + p → ν + p, in liquid scintillator detectors is proposed for the measurement of ν x since it is the only neutral current channel providing spectral information. Though the total cross section of this process is about a factor of four smaller than the cross section of inverse-beta decay, the yield above a realistic threshold of about 200 keV is of the same order since this reaction is possible for all six neutrino types [9].
This detection channel plays an essential role for large scale deep underground liquid scintillator detectors like KamLAND, Borexino and SNO+, in the pending measurement of the total energy and temperature of ν μ , ν τ ,ν μ andν τ . It should be noted, however, that only a fraction of the deposited kinetic energy of highly ionizing particles like protons is visible. Hence to properly reconstruct the true kinetic proton energy E p from the visible energy E vis p , the light response of the liquid scintillator for these particles has to be known. The ionization quenching of the light yield is well-described by Birks' law [11] and can be quantified in terms of the parameter product kB known as Birks' constant and a second parameter C (Eq. (7)). While the proton light response function of the KamLAND scintillator, which consists mainly of normal-paraffine (80 %) and pseudocumene (20 %), is already measured up to 10 MeV [12], the response to protons of scintillators based on linear alkylbenzene (LAB) has not been measured before. LAB is employed for novel liquid scintillators which are used or planned as the neutrino target in recently commissioned, upcoming and potential future neutrino experiments like Daya Bay, RENO, SNO+ and LENA. Also further potential experiments like LENS, HANOHANO and Daya Bay II are considering LAB as possible solvent for their scintillator.
This article presents the measurement of the light response function L(E) for protons relative to the electron response function L e (E), carried out at the accelerator facility PIAF of the Physikalisch-Technische Bundesanstalt (PTB) for different LAB based scintillators. Recoil protons are produced in the liquid scintillator by neutron-proton scattering events using a neutron beam with continuous energy distribution. The maximum energy of the recoil protons is determined from the energy of the incident neutrons, which is measured using the time-of-flight (TOF) method. Hence, the proton light output function can be determined from the position of the recoil proton edge in the pulse-height spectra produced by mono-energetic neutrons. The presented method was already used successfully at the PTB for the characterization of NE213, 1 BC501 2 and BC501A scintillation detectors [13][14][15][16].

Measurements
In the present experiment, the light output function of binary and ternary LAB based scintillation systems is determined. In these solutions, 2,5-diphenyloxazole (PPO) is employed as primary fluor and p-bis-(o-methylstyryl)-benzene (bis-MSB) as secondary fluor, which acts as wavelength shifter. Two concentrations of PPO in the LAB solvent, 2 g/l and 1 Nuclear Enterprise Ltd. 2 Saint-Gobain Ceramics & Plastics, Inc. 3 g/l, are studied, each with and without bis-MSB, for which a concentration of 15 mg/l is chosen.
After dissolving the solute in LAB, the solutions were purified from oxygen by bubbling with argon for 30 min. Without further contact to air, the solutions are subsequently filled in a scintillation detector consisting of a cylindrical dural cell with one port covered by a window made of ground and polished fused silica for higher UV transparency. The diameter and inner height of the cell is 50.8 mm. The cell is coupled to a XP2020Q 3 photomultiplier tube (PMT), also equipped with a fused silica window, by a conical UV transparent poly(methyl methacrylate) (PMMA) light guide. 4 The inner walls of the cell and the light guide are coated with BC-622A 5 reflective paint to increase the light collection efficiency. The detector is equipped with a LED gain stabilization system which regulates the high voltages such that the maximum drift of the gain is smaller than 1 %. Standard NIM modules are used to obtain a pulse-height (PH) signal from the ninth dynode out of twelve, as well as a pulse-shape (PS) signal and a time-of-flight (TOF) signal from the anode of the photomultiplier. The ninth dynode is chosen for obtaining the PH signal to avoid a nonlinear gain of the PMT. The PH, PS and TOF signals are registered by a multi-parameter data acquisition system and stored event-by-event. The PS discrimination module, as described in detail in [17], uses the zero crossing method to derive a timing marker from the trailing edge of the anode pulse. The time difference of this PS timing marker to the TOF timing marker derived from the fast leading edge of the anode pulse, using a constant fraction discriminator (CFD), is measured with a time-to-amplitude converter (TAC). It is related to the decay time of the scintillation light and therefore to the ionization density of the charged particle producing the event. The non-linearity of the TAC converter is measured using a time calibration module which produces a set of pulses at multiple intervals of t = 40 ns. The variations of the peak positions from the nominal values quantify the non-linearity and are determined to be smaller than ±0.23 ns. Figure 1 shows the registered events sorted into a PS versus PH matrix independent of the energy of the incident neutron. The structures due to neutron-induced events ('n'), like recoil protons or other light charged particles, are partially separated from those caused by Compton electrons ('γ '). The electron events are due to ambient γ -radiation and photons resulting from the inelastic neutron scattering on carbon nuclei in the scintillator and on aluminum nuclei of the walls. A suppression of events with a PS signal below the indicated line strongly reduces the influence of photon-induced events with respect to the PH spectra of neutron-induced events. However, for LAB based scintillators a fraction of misidentified photon-induced events remains in the overlap region, as indicated by Fig. 1. This is mainly effected by the small number of photoelectrons produced per scintillation event which is primarily due to the wavelength mismatch between the emission spectrum and the sensitivity of the photocathode.
The energy calibration of the PH spectra and determination of the light response to electrons for each scintillator is obtained from the position of the Compton edges produced by a set of standard γ -ray sources ( 137 Cs, 22 Na, 207 Bi). Ambient background is measured separately and subtracted from each source measurement. The analysis of the calibration data is described in detail in Sect. 3.
To obtain a proton light output function over an energy range as large as possible in a single experiment, a neutron beam with a continuous spectral distribution is used, which is produced at the PTB TOF facility by bombarding a stopping-length beryllium target with a 19 MeV proton beam from the isochronous CV28 cyclotron [18]. At an emission angle of 0 • , the maximum energy of the neutrons is 17.15 MeV. The center of the cell was positioned at a distance of 12.11 m from the Be target under 0 • with respect to the proton beam. The distance between target and detector cell is measured with an accuracy of ±2 mm. To consider the target thickness of 3 mm and include the uncertainty of the effective center of the cell caused by the fluence attenuation, ±6 mm are conservatively assumed. With a beam current of about 50 nA, the neutron rate at the detector front face was about 3 × 10 4 s −1 . The repetition frequency of the proton beam was set to 481.3 kHz. The trigger rate of the CFD module was about 3 × 10 3 s −1 for a threshold set close to the electronic noise. The probability for multiple neutron interactions within one beam pulse is thus smaller than 1 % and does not cause a deterioration of the shape of the pulse-height spectra. In the given configuration the neutron energy threshold for time-frame overlap is 0.15 MeV, which corresponds to an electron-equivalent energy for the given scintillators of less than 20 keV. Hence, with a hardware threshold of about 200 keV, time-frame overlap is avoided.
Proton bunches that are not entirely deflected by the internal beam pulse selector system cause additional pulses, so-called satellite pulses, on the target. The spectral distributions of outgoing satellite neutrons are assumed to be the same as the ones from the main proton bunches and are visible in the PH vs. TOF matrix as faint distributions which are shifted in TOF relative to the distribution resulting from the main beam pulse. Hence, satellites cause spurious events above the recoil proton edge in the pulse-height spectra which could affect the determination of the position of the recoil proton edge. A subtraction of the shifted TOF spectra, though, showed no significant influence on the position determination of the proton recoil edge.
To improve the dynamic range of the pulse-height measurement for neutron energies below 5 MeV, each measurement is repeated with an increased amplification. This amplification mode is referred to as high gain (HG) mode, the other one as low gain (LG) mode. For the HG measurement, n/γ -separation by pulse-shape discrimination (PSD) is not possible. However, at these energies, inelastic scattering on carbon does not contribute to the neutron induced response because of the energy of the first excited state in 12 C which is at 4.439 MeV. The contribution from prompt γ -radiation from the target and from ambient γ -radiation is subtracted experimentally.

Calibration
The amount of detected scintillation light and the pulseheight resolution strongly depend on the specific composition of the scintillator. Therefore the pulse-height calibration of the scintillation detector is repeated for each filling and amplification mode using three γ -ray sources with a total of six γ -lines, namely 137 Cs (E γ = 0.662 MeV), 22 Na The position of each Compton edge, with energy where all energies are given in MeV, is defined by comparing the measured pulse-height spectrum with a distribution simulated using the Monte Carlo photon transport code GRESP [19]. To adjust the theoretical shape to the measurement, the simulated spectrum is folded using a Gaussian resolution function with a pulse-height dependent width and fitted to the measured pulse-height spectra (Fig. 2). In this way the measured pulse-height signal PH is calibrated in units of the maximum kinetic energy E e of the Compton electrons. The calibrated pulse-height (light output) L is given by

Fig. 2
Pulse-height spectrum of 137 Cs given in pulse-height channels. The calculated Compton electron spectrum before and after folding with the Gaussian resolution function is compared to the experimental spectrum for LAB with 2 g/l PPO and 15 mg/l bis-MSB, measured in HG mode where n denotes the calibration factor. The relative resolution L at pulse-height L is parameterized by [20] The fit parameters α, β, γ can be attributed to the individual contributions from spatial dependence of the light collection efficiency, statistical variation of the number of photoelectrons and electronic noise, respectively. The values of the resolution parameters for each scintillator filling are obtained by fitting Eq. (3) to the locally determined full width at half maximum L at each Compton edge position. Figure 3 shows the resulting resolution functions for every sample as a function of the light output. The lower resolution for the scintillators without wavelength shifter is mainly attributed to detector properties. The reflectivity of the BC-622A coating inside the cell is reduced below 420 nm and the XP2020Q PMT has its peak sensitivity at 420 nm, while the emission peak of PPO lies around 360 nm. Thus the spatial dependency and statistical variation are increased. However, the properties are more fortunate for the scintillators with bis-MSB which has a maximum emission around 420 nm. For this reason, besides the two binary scintillators, the two systems with additional wavelength shifter are investigated, though no difference in the ionization quenching properties is to be expected.
According to [11], the light output function L(E) is related to the stopping power dE/dx of a charged particle of kinetic energy E in the traversed scintillator via Fig. 3 Relative pulse-height resolution L/L for electrons as function of the light output L in electron-equivalent energy. The respective solute admixture to LAB is given in the legend. Shown are the data points acquired from source measurements in HG mode and the respective adapted resolution function Eq. (3). In the shown total uncertainties, the single contributions are added quadratically if the particle is stopped in the scintillator. This relation is known as Birks' formula with the product kB being the Birks constant. The light output L is the total light emitted when a charged particle loses all of its energy E within the scintillator. For electrons with energies ≥125 keV, the stopping power becomes very small [11] and Eq. (4) reduces to the linear expression The subscript 'e' is added here to underline that this description of the light output only holds for electrons. The parameter S denotes the chosen scaling between electron energy and light output and is set to S = 1, which yields the light output in electron-equivalent energy, if the calibration procedure described above is used. The constant E 0 accounts for the fact that the response of the scintillator to electrons is only approximately linear. The non-linearity at small energies leads to an energy offset, which is assumed to be E 0 = 5 keV for NE213 and BC501. Recent measurements with LAB and EJ-301 scintillators [21], though, revealed a non-linear behavior already starting at 400 keV, resulting in an offset of about 50 keV. Within the present work, the potential non-linearity of the electron light output can not be discriminated from the effect of the remaining electronic offset of the amplifiers which was corrected for during the analysis of the pulse-height calibration, i.e. E 0 is set to zero.
A comparison of Eq. (2) and Eq. (5) gives the relation between PH and electron energy where the second equality only holds for the case described here. The calibration factor n −1 can now be derived from a linear fit to the positions of all Compton edges in the PH spectrum as performed in Fig. 4 and is independent of the total offset of the pulse-height scale. The values obtained in the present measurements are summarized in Table 1.

Proton light output function
For the determination of the proton light output function L(E), basically the same technique is used as for the analysis of the photon-induced events in the calibration measurements described in Sect. 3.1. The events measured with the white neutron beam are first sorted into a TOF vs. PH matrix. As mentioned above, PSD is used to select only events produced by light charged particles for the spectra obtained with the LG setting of the amplifier, while in the HG mode n/γ -separation by PSD is not feasible. The TOF range covers the full neutron energy range from about 1 MeV to 17.15 MeV. The TOF scale is established using the arrival time difference between prompt gammas and neutrons of  maximum energy and a time determined using an electronic time calibration module. The prompt gamma line is broadened due to the finite time resolution of the accelerator and the uncertainty of the centroid position of the resulting peak is calculated to be ±0.1 ns. Calibration gives a TOF channel width of 0.7764 ns and the relative uncertainty of the time calibration is 0.05 %. Pulse-height spectra induced by quasi mono-energetic neutrons are then extracted by collecting all events contained in small windows around the TOF value corresponding to an individual neutron energy. The size of the TOF window is always made smaller than the pulse-height resolution L for protons of the respective energy. A second window of the same size is placed on the random background at TOF values smaller than those of neutrons with the maximum energy and the extracted background spectrum is subtracted from the signal spectrum. Pulse-height spectra for different incoming neutron energies are shown in Fig. 5. The energy of the recoil protons, pro-ducing the edge in the pulse-height spectra, equals that of the incident neutrons. Before extracting the position of the recoil edge, the dependence of the TOF measurement on the pulse-height is investigated. The so-called time-walk is caused by imperfections in the constant fraction timing technique employed in the CFD module which show up in particular at small pulse-height. It leads to a distortion of the PH vs. TOF matrices and finally in the extracted proton recoil spectra at low pulse-heights. A correction of the PH vs. TOF matrices of LAB with 2 g/l PPO and 15 mg/l bis-MSB revealed a deviation of about only one channel in the TOF spectrum. For comparison, the smallest window width, which is chosen at high neutron energies, has a width of five channels. A spectral distortion in the region of small pulse-heights, however, has the biggest influence on the definition of the proton recoil edge at low neutron energies, where the TOF window width has to be increased because of a decreased resolution, as described above. The maximum width of the TOF windows amounts to 25 channels, which is much larger than the observed deviation in the measured TOF. Hence, time-walk has a minor effect. As a correction of the PH vs. TOF matrices, to account for time-walk, can bring in new artefacts due to the binned structure of the data, no correction is performed within this analysis.  5 Calculated (line) and measured (histogram) proton recoil spectra at different incoming neutron energies. The measured spectra are extracted from the data sets for LAB with 2 g/l PPO and 15 mg/l bis-MSB. For the spectrum at 3 MeV neutron energy, HG data is used, for the rest LG data. The calculation is performed using the NRESP7 code and does not consider a threshold (leading to counts starting at channel 1, where there is no data available) and photon induced background (leading to less counts at low pulse-heights above threshold compared to data) To calculate the detector response to monoenergetic neutrons, the Monte Carlo code NRESP7 [22] is used. NRESP7 models all reactions of neutrons in the scintillator, in the detector housing and in the Lucite light guide, including the production of secondary charged particles, but it does not simulate the interaction of deexcitation photons (e.g. from the first excited state in 12 C at 4.439 MeV) resulting from inelastic neutron scattering. Instead, it is assumed that such events are suppressed by PSD. If the incoming neutrons have enough energy to excite 12 C, signals from deexcitation photons with pulse-heights up to about 180 channels can add up to the proton recoil spectra. As this background is not fully suppressed by PSD, it leads to more counts above threshold in the measured pulse-height spectra than in the simulated ones (see Fig. 5) for neutron energies higher than ∼4.5 MeV. However, this mainly affects the proton recoil edge determination, if the recoil edge is in the same region as the photon induced background. This is only the case for a maximum proton energy (i.e. incoming neutron energy) up to ∼7 MeV. For the fit of calculated to the measured spectrum, a pulse-height interval is chosen, which starts below the proton recoil edge and ends above. For neutron energies above about 7 MeV, this interval is chosen such, that photon induced background is excluded. The pulse-height given in channels is finally translated into light output in electron-equivalent energy L using the factor n −1 obtained in the calibration with photon sources. The resulting relative proton light output resolution functions are plotted in Fig. 6.
The simulated pulse-height spectra are adapted to the experimental ones in an iterative process. As an initial step, an existing reference light output function for NE213 detectors is used to calculate pulse-height spectra using NRESP7. The simulated spectra are compressed or stretched and adjusted in height to fit the experimental spectra in the pulseheight region around the recoil edge and the respective compression or stretching factors are applied to correct the light output function. The corrected function is inserted as input for the next iteration step. The behavior of the applied light output function follows a straight line for E ≥ 8 MeV, where the parameters of the line are gained from a fit to the extracted data points. After about four steps, the procedure converges. Figure 7 compares the initial light response function with the final one for LAB with 2 g/l PPO and 15 mg/l bis-MSB. The data points extracted in this way are presented in Fig. 8 together with the electron data points gained from calibration. The results for the remaining scintillators are shown in Figs. 9-11. The contributions to the total uncertainty in L and E are summarized in Table 2.
The most popular analytical description for the energy dependency of the light output is given by Birks' formula Eq. (4), assuming ionization quenching to be the reason for the reduced light output. This expression, though widely used for all kinds of scintillators, was originally developed for anthracene crystals. However, the underlying physical processes in liquid scintillators are different and this description should be considered as being semi-empirical. For large values of dE/dx, such as for incident ions, Eq. (4) is often encountered in its generalized form including a quadratic correction term parameterized by C, as proposed in [23]. Within the present work, electron and proton light output functions of a particular scintillator are measured under identical conditions so that the scaling factor S for proton data is the same as the one for electron data, i.e. S = 1. The Birks constant kB and the parameter C for LAB based scintillators are obtained by means of a χ 2 fit of the theoretically expected light output L theo in Eq. (7) to the experimentally determined values L exp presented in Figs. [8][9][10][11]. Extracted data points, instead of the final input function in NRESP7, are used for this fit to be independent from the mentioned straight line fit for E ≥ 8 MeV. The stopping powers dE/dx for protons are calculated using the code SRIM [24]. The χ 2 calculation is carried out following the pull approach [25], where . The data is taken with LAB, 2 g/l PPO and 15 mg/l bis-MSB. In the shown total uncertainties, the single contributions (listed in Table 2) are added quadratically. For some values of L(E), the resulting error bars are smaller than the marker size Within this calculation K = 8 systematic uncertainties, listed in Table 2, are included as nuisance parameters. f k n describes the fractional change of the n-th value of L theo if the k-th source of systematics is varied by 1 σ k and ξ k is a standard normal deviate. The normalization condition for the ξ k 's is realized through quadratic penalties, summed over the K sources of systematics. The statistical uncertainty, denoted by u, is also indicated in Table 2. For each combination of kB and C, Eq. (8) is minimized with respect to all ξ k 's using MINUIT and the minimal value, χ 2 pull , is stored in a (kB, C) map. The best fit values are finally extracted at the minimal χ 2 pull in the full parameter space. The uncertainties are extracted from the χ 2 pull projections on the  Table 1) Statistical uncertainty 1σ stat.
recoil edge position <0.02 % kB resp. C axis. The best fit value of kB has a two-sided limit and the kB values at χ 2 pull − χ 2 pull,min = 1 give the 1σ uncertainty. C has a one-sided limit and the C value at χ 2 pull − χ 2 pull,min = 1.645 refers to an upper limit at 95 % confidence level. Besides the experimental uncertainties, a further uncertainty has to be considered, namely the uncertainty of the stopping power calculation entering Eq. (7). With the SRIM code, a determination of the stopping power for elemental materials is possible with an accuracy of a few percent [26], thus the described fitting procedure is repeated for dE/dx + σ dE/dx and dE/dx − σ dE/dx . Assuming a shift of 2 % of the stopping power already leads to a 2 % difference in the resulting quenching parameters. This uncertainty is added quadratically to the experimental one and is the dominant contribution to the resulting total uncertainty. The fit results are summarized in Table 3. C is consistent with zero for all investigated scintillator samples with an upper limit (95 % CL) of about 10 −7 cm 2 MeV −2 and the values for kB range between (0.0094 ± 0.0002) cm MeV −1 and (0.0098 ± 0.0003) cm MeV −1 .
According to current knowledge, all ionization quenching processes are primary processes in the scintillator, i.e. processes that transfer ionization energy to excitation energy of the solvent. Quenching processes compete with the excitation of the solvent molecules into π -electron singlet states [11]. This is the primary process which is essential for the final scintillation emission in unitary systems as well as in binary and ternary systems. Consequently, for scintillators with the same solvent, the magnitude of ionization quenching should be the same independent of possible so- lutes. Within the presented measurement no significant deviation from these considerations could be observed. For practical use, the proton quenching factor Q p , the ratio of proton and electron light output functions for each examined LAB scintillator is presented in Fig. 12. Above a kinetic energy of 8 MeV, the proton light output is about half the electron light output. Below that energy the quenching factor strongly decreases, i.e. the proton light output is further reduced down to <20 % of the electron light output for E ≈ 1 MeV.

Application of the results
The detection of supernova neutrinos of all flavors using neutrino-proton elastic scattering is discussed in detail in [10]. Within that reference, different experiments capable of the respective detection are compared, namely Borexino, KamLAND, SNO+ and LENA. For the LAB scintillator detector SNO+, a Birks' constant of kB = 0.0073 cm MeV −1 is used. However at that time no measurement of proton Fig. 12 Proton quenching factor Q p as function of kinetic energy E for LAB based scintillators. The respective solute admixture is given in the legend quenching in LAB was published. The measurements presented here for LAB with 2 g/l PPO (as will be used in SNO+) reveals with kB = 0.0098 cm MeV −1 a stronger proton quenching. To show the impact of a stronger quenching on the event yield induced by a supernova, the same detector (i.e. SNO+) and assumptions for the supernova are used as in [10]. In the same reference, the prospect of the larger detector LENA is discussed, assuming kB = 0.01 cm MeV −1 and 3.3 × 10 33 free protons in the liquid scintillator. As LENA favors LAB as solvent with 3 g/l PPO and 20 mg/l bis-MSB, 6 this detector is also discussed in the present article with the obtained result of kB = 0.0098 cm MeV −1 for LAB with 3 g/l PPO and 15 mg/l bis-MSB. For both detectors, C is considered to be zero in agreement with the results summarized in Table 3.
The supernova is assumed to occur at a distance d = 10 kpc, releasing over the duration of the burst a total energy of ε = 3 × 10 53 erg which is equipartitioned among all flavors. The energy is distributed according to where dF α /dE ν is the neutrino fluence per flavor. The energies are given in MeV, d in 10 kpc and ε α in 10 52 erg. The mean energy E ν α of the individual flavors differs and is exemplarily set to 12 MeV for ν e , 15 MeV forν e and 18 MeV for ν x , following the approach in [10]. The respective fluence in ν e ,ν e and the sum of all ν x as well as the total fluence are shown in Fig. 13. The fluences are of the order of 10 10 cm −2 MeV −1 and the neutrino energy E ν ranges up to about 60 MeV. For these low neutrino energies, the differential neutrinoproton elastic scattering cross section [27] can be simplified to (considering only the lowest order in E ν /m p ) and the cross section for neutrinos and anti-neutrinos can be regarded as identical [9]. For the last step, a proton mass of m p = 938 MeV is used and the neutral-current coupling constants are c V = 0.04 and c A = 0.635 [9]. E ν and the proton recoil energy E p are given in MeV. Using Eq. (10) and Eq. (11), the true recoil energy spectrum of recoil protons produced by supernova neutrinos is given by where N p is the number of free protons in the scintillator target. The minimum neutrino energy needed to accelerate a proton to energy E p is ν,min ≈ m p E p /2 [9]. Figure 14 shows the expected proton recoil spectra for the various flavors in a detector with N p = 5.9 × 10 31 (corresponding to the SNO+ target mass of about 0.8 kt [28]). Despite the small cross section, the intense neutrino fluence apparently causes a significant neutrino-proton elastic scattering yield of the order of 10 2 events per kt LAB. The proton energies range from 0 MeV to about 5 MeV. However, proton quenching markedly changes the shape of proton energy spectra and shifts the true proton recoil energies to lower observed energies. Even though there is no energy threshold to produce scintillation light, below 0.2 MeV organic liquid scintillators are strongly affected by β-decays of 14 C, which is intrinsic to organic solvents. The high decay rate dominates this region and thus sets an effective threshold of 0.2 MeV (further backgrounds are not taken into account within this article as they are described in detail in [9] and found to be negligible for the presented supernova signal channel). Hence, the observed event yield finally depends on the quenching strength. The expected event yield above threshold will be discussed in the following for SNO+ and LENA. The respective detector specifications are given in Table 4. The main difference between both detectors is the amount of target protons, which is about 5.9 × 10 31 in SNO+ and 3.25 × 10 33 in LENA, given a proton density in LAB of 6.31 × 10 22 cm −3 [29].  Table 4 Detector properties assuming LAB with 2 g/l PPO as target in SNO+ and LAB with 3 g/l PPO and 15 mg/l bis-MSB in LENA. The SNO+ target mass and energy resolution is the same as employed in [10] for better comparison of the results and in agreement with [28]. The LENA target mass and resolution is taken from [30], where the expected resolution is quoted to be ≥200 p.e./MeV. The number of free protons N p is determined using a proton density in LAB of 6.31 × 10 22 cm −3 [29]. C is assumed to be zero To convert the true proton energy spectra in Fig. 14 into spectra of visible energy E vis p in electron equivalent, the quenching factor Q p (E p ), as presented in Fig. 12, is applied to Eq. (12). The resulting spectra are folded with the energy resolution and scaled according to the respective number of free protons as given in Table 4. Figures 15 and 16 finally show the expected supernova neutrino signal spectra from ν-p elastic scattering in SNO+ and LENA under the given assumptions for the supernova scenario and the detectors. The event yield above a threshold of 0.2 MeV, using kB = 0.0098 cm MeV −1 and C equal to zero for LAB, is about 98 events in SNO+ and 5403 events in LENA. The yield scales linearly with the target size. For LAB scintillators it is in the range of ∼123 events/kt, which makes this detection channel very promising for large scale detectors. Applying the listed energy resolution only changes both yields by less than 0.5 % compared to perfect resolution and seems to have a negligible effect. The dependence on energy resolution thus is not further studied here, however, it needs to be investigated when a profound sensitivity study to supernova neutrino fluence parameters is performed. As mentioned earlier, the determination of kB is most sensitive to the calculation of the stopping power dE/dx and vice versa the calculation of Q p for a fixed Birks' constant. In [10] dE/dx tables from PSTAR [31] are used while for the presented analysis SRIM tables are employed. Though the stopping power data is similar, the expected yield for SNO+ is recalculated with kB = 0.0073 cm MeV −1 and C considered to be zero, using SRIM tables, to show the impact of a higher kB value without bias. Figure 17 compares the true proton energy sum spectrum and the observed sum spectra assuming kB = 0.0073 cm MeV −1 and kB = 0.0098 cm MeV −1 respectively. The higher quenching reduces the expected event yield by about 16 % from ∼116 events 7 to ∼98 events. This figure demonstrates how proton quenching pushes supernova ν-p elastic scattering signal events to lower energies and that the yield above threshold strongly depends on the strength of the quenching. Therefore it is of great importance for liquid scintillator detectors, sensitive to supernova neutrinos, to properly measure the proton quenching parameters of the used scintillator to extract valuable information about supernova neutrino fluences of all flavors.

Summary and outlook
This article presents a measurement of the light output function L(E) of protons relative to the one of electrons for four different LAB based scintillators. A measurement of the relative proton light output reveals the magnitude of light output reduction due to ionization quenching. Fitting the functional description of L(E) proposed by J.B. Birks, parameterized by kB and C, to the experimental data yielded values for kB ranging between (0.0094 ± 0.0002) cm MeV −1 and (0.0098 ± 0.0003) cm MeV −1 and an upper limit of about 10 −7 cm 2 MeV −2 (95 % CL) for the parameter C for the investigated scintillators.
It should be noted, however, that the parameters within Birks' law, kB and C, alone are not conclusive for the comparison of different scintillators (i.e. scintillators with different solvents) as they are correlated with the stopping power of the scintillator material. Consequently, to decide between two different scintillators, always the full light output functions need to be compared. Different solutes in the same solvent, on the other hand, do not show a difference in the calculated stopping powers due to their small concentration. For the latter case, kB and C are thus sensitive to a potential influence on the strength of ionization quenching by different solutes. The presented measurement does not reveal a significant disparity of the resulting parameters of all LAB based scintillators and is in accordance with the expectation that all ionization quenching processes are primary processes.
A further remark is that the values for kB extracted from measurements of the non-linearity in the light output of lowenergy electrons is not necessarily identical to the one of heavy ionizing particles like ions [32]. Whereas it is assumed that the kB value is the same for different ions, like protons and alphas, once the conditions of the measurement are fixed, as discussed in detail in [33]. This hypothesis can be tested with the taken data by determining the α light output function for the investigated scintillators. 12 C(n, α) 9 Be and 12 C(n, n 3α) reactions, induced by neutrons from the beam, lead to characteristic structures in the pulse-height response of the scintillator which can be reproduced with the NRESP7 code. A proper knowledge of the scintillator response to α's is important for the development of a sophisticated background model of a liquid scintillator detector, including α-induced backgrounds. This analysis and discussion though is beyond the scope of this article and an outlook to future possibilities.
Using the measured proton light output functions for LAB solutions, the potential of LAB based scintillator detectors to observe supernova neutrinos of all flavors via ν-p elastic scattering is discussed and it is shown how a different Birks constant affects the event yield. Due to this impact on the observation, well-known quenching parameters are important for liquid scintillator detectors, capable of supernova neutrino detection.