1 Introduction

Orbiting X-ray observatories are subjected to the impact of charged particles of different origins, e.g. cosmic rays, solar wind and solar events, electrons and protons trapped in the Earth magnetosphere, etc. Amongst them, protons with typical energies ranging from a few tens to a few hundreds of kiloelectronvolts, the so-called ‘soft protons’ (SP), are scattered and funneled towards the focal plane when impinging at grazing incidence on X-ray optics. Once they reach the instruments at the focal plane, they significantly contribute to the level of non-X-ray background (NXB), and potentially damage the detectors in the most severe cases. The effects of SP were first observed on the Chandra X-ray Observatory [12] and on XMM-Newton [6], both launched in 1999. To this day SP affect the operability of these X-ray missions by significantly reducing their useful observing time and their duty cycles – for instance, the observing time of XMM-Newton is reduced by \(\sim \)30–40 % [5].

SP constitute an even bigger issue for Athena (Advanced Telescope for High Energy Astrophysics, [9]), an L-class ESA mission, planned to be launched in the early 2030s. The telescope will be equipped with Silicon Pore Optics (SPO) [13], wafers of silicon coated with high-Z materials, stacked one upon the other and arranged to cover the whole pupil of the telescope.

Athena will orbit around the second Lagrangian point L2Footnote 1 of the Earth-Sun system, at a distance of about 1.5 million kilometer from Earth in the anti-Sun direction. This orbit lies in the tail of the Earth magnetosphere, where the particle environment is extremely variable [7]. At present, the science driven requirement for SP is that their flux at the focal plane has to be < 5 × 10− 4 cts s− 1 cm− 2 keV− 1 (corresponding to 10 % of the total NXB), in the 2–10 keV energy band, for 90 % of the observing timeFootnote 2. In order to have a thorough estimate of the SP flux expected at the focal plane, knowing the scattering efficiency of soft protons at grazing incidence from Athena’s SPO is fundamental. To achieve this goal, we performed experimental scattering efficiency measurements with a SPO sample.

Similar experimental measurements have already been performed for XMM-Newton [10] and eROSITA [2, 3] mirror samples. Especially the measurements on eROSITA showed that the scattering is always non-elastic, with energy losses of the order of a few tens of kiloelectronvolts, depending on the energy of the incident beam. The scattering efficiency always peaks close to the angle of specular reflection and is higher for lower incident angles, while it is almost independent of the energy of the impinging protons.

A model for the scattering efficiency of ions at grazing incidence from the surface of a material was developed by [11], in non-elastic approximation. However, the model is generally valid for any surface and contains some approximations, therefore a simple application of the Remizovich formula to the eROSITA data did not exactly reproduce the efficiencies measured in the laboratory [1]. A fit of the experimental data with the Remizovich formula led to the evaluation of the parameter enclosing the micro-physics of the scattering, so that a new analytical semi-empirical model that better reproduces the data from the eROSITA mirror sample was derived [1]. This model can be used to assess the SP flux expected at the instrumental focal plane of the satellite. The model can also be applied to Athena, provided that the proton scattering properties of the SPO are experimentally determined.

In this publication, we present the first measurements of scattering efficiencies of low-energy protons off a SPO sample. These experimental activities were conducted in the framework of the EXACRAD (Experimental Evaluation of Athena Charged Particle Background from Secondary Radiation and Scattering in Optics) project funded by ESA. The paper is structured as follows: we first describe the laboratory set-up and the elements along the beam line in Section 2; in Section 3 we illustrate how the scattering efficiency is derived from the raw data; the new data are presented in Section 4 and are compared to the data from eROSITA as well as to the semi-empirical model mentioned above in Section 9; finally, we draw our conclusions in Section 5.

2 Experimental setup

The experiment was conducted at the 2.5 MV Van de Graaff accelerator at the Goethe University (Riedberg Campus) in Frankfurt am Main. The setup of the beamline, similar to that of [2, 3], is given in Figs. 12, and 3.

Fig. 1
figure 1

Schematic drawing (not in scale) of the beamline setup. The proton beam enters the setup from the left-hand side. It encounters the pinhole aperture (1 mm in diameter), the Al degrader foil (0.002 mm thick) and the collimator. Inside the target chamber, the normalisation detector can be lowered down to intercept the beam for the normalisation measurements. If the normalisation detector is not in the line of the beam, then protons are reflected from the SPO sample (in yellow) towards the detector chamber, where they hit the central and lateral detectors. The incident angle 𝜃0 between the line of the beam and the mirror varies with the inclination of the target plate, while the scattering angle 𝜃 between the mirror and the detectors in the detector chamber varies with the their height h. The distance d between the target plate and the vertical ax of the central detectors is fixed to 942 mm

Fig. 2
figure 2

A CAD model of the beamline (same as[2]). The proton beam enters the setup from the right and moves towards the left. The SPO sample is located in the target chamber, while the detector is placed in the chamber at the end of the beamline (detector chamber). A second detector (not shown in the picture) was placed next to the central one, with an angular distance of \(\sim 2^{\circ }\)

Fig. 3
figure 3

Picture of the beamline at the Van der Graaff accelerator facility. The beam direction is from the right to the left. The position of the degrader foils and the collimator inside the beamline are pointed out, as well as the detector and target chambers

2.1 Beamline setup

Protons enter the beamline through a copper pinhole aperture of the diameter of 1 mm, which reduces the size of the incoming beam to prevent pile-up and to maintain reasonable rates on the detectors. Successively, the beam goes through a 0.002 mm thick aluminium foil, which degrades the incoming beam energy below the lower limit of the accelerator. The degraded beam enters, at this point, a 78 cm long collimator, which directs part of the widened beam directly to the target. Two further apertures are positioned at the entrance and at the exit of the collimator, both with a diameter of 1 mm. This diameter limits the smallest possible incident angle to \(\sim 0.5^{\circ }\). The apertures are supported in their position by 2 mm aluminum plates, which absorb any proton of the degraded beam not entering the apertures and being scattered by the inner walls of the collimator and of the beamline itself.

The SPO target (Fig. 4), provided by cosineFootnote 3, consists of a 110 mm long single silicon wafer, 0.775 mm thick, grooved in the bottom, and coated on top with a 10 nm of iridium and 7 nm of silicon carbideFootnote 4. It is located in an apposite chamber (hereafter called target chamber) and mounted on a tiltable plate. The height of the target can be adjusted by a set of screws underneath the plate. A linear manipulator is used to change the inclination of the plate, i.e., of the incident angle (𝜃0). The pivoting point is is several centimeters below the line of the beam, so that the target can be completely removed from the beam, allowing for a determination of the primary beam position on the detector plane. The manipulator is set below the target chamber and hence can be easily accessed when the system is under vacuum.

Fig. 4
figure 4

Pictures of the SPO sample, with the reflecting surface up (left) and down (right)

Between the exit of the collimator and the target plate, a Passivated Implanted Planar Silicon (PIPS) detectorFootnote 5 is mounted on a push-pull manipulator, at the same height of the beamline. This detector is used to register the amount of flux of the incident beam impinging on the target, useful to have normalisation measurements. This detector will be called hereafter ‘normalisation detector’. The push-pull manipulator permits a fast removal of the detector, guaranteeing a measure of the impinging proton flux (Φinc, cfr. (2)) for each measure of the scattered beam (see Section 3 for the need of having frequent normalisation measurements). An aluminium blind with an aperture of 3 mm is set on top of the normalisation detector to avoid saturation. Lastly, downstream of the target chamber, a thick aluminum sheet, with a slit of 3 cm height and 1 cm width, is installed a few centimeters after the target plate. This window lets pass only the protons on the line of the beam, while the sheet absorbs all the ones that have been scattered by the inner walls or by the other elements in the target chamber.

At the end of the beamline, a second chamber (hereafter detector chamber) hosts two more PIPS detectors, called ‘central detector’ and ‘lateral detector’, respectively, used to register the on-axis and off-axis fluxes (Φscat(𝜃0,𝜃,ϕ), cfr. (2)) of the beam scattered by the target. They are mounted on a second linear manipulator, which allows a spatially resolved sampling of the scattered beam. The distance between the center of the target plate and the detection plane is 942 mm. The central detector is aligned with the beam, while the lateral detector is set on the left. This configuration allows for a coverage of the scattered beam on the the incident direction and at an azimuthal angle ϕ of 1.97± 0.13. On top of each detector there is a blind with an aperture of a diameter of 1 mm for the central detectors and of 3 mm for the lateral detector, respectively. They reduce the solid angle of the detectors with respect to the mirror center to about 8 × 10− 7 sr and 2 × 10− 5 sr for the central and lateral detector, respectively.

2.2 Data acquisition chain

The pulse signal produced by the PIPS is amplified and digitalised trough several analogical/digital electronic components. A flow chart is given in Fig. 5.

Fig. 5
figure 5

Data flow of the electronic chain for the acquisition of the experimental data. The analogical signal from the PIPS detectors first goes through a pre-amplifier and an amplifier, then it is converted into a digital signal by the ADC, and finally it is stored in the histogramming memory. Contemporary, a pulse genarator sends signal at a fixed frequency to the ADC and to a scaler. The digitised signals are read out by a CAMAC controller unit, which transmits them to a computer once the measurement is finished

The PIPS detectors produce a pulse with an amplitude proportional to the energy of the incident particle. The pulse signal from each PIPS goes through its own pre-amplifier and amplifier and it is then digitalised by the Analog to Digital Converter (ADC). The ADC receive the continuous signal (from 0 to \(\sim \)10 V) from the three channels – one for each detector – and convert them into discrete signals, distributing it into 8192 bins, with a resolution of 1.22 mV. The digitised signals are then passed to the histogramming memory, which produces an histogram for each channel. Once the measurement is done, the histograms are read out by the CAMAC module and are transferred to a computer, which acquires and stores them as raw data files.

The process of digitalisation of the data within the ADC takes a certain time (fractions of seconds), so that if a new signal comes within that time, it is not registered. To account for this dead-time, a pulse generator, which generates pulses at a fixed frequency, is connected to the ADC and to a scaler, which counts the number of pulses produced by the pulse generator during the acquisition time. The scaler is also fed to the CAMAC control module. The difference between the reading of the counts from the ADC and that from the scaler gives the dead-time correction factor (cfr. (4)). The pulse generator fed to the ADC constitutes another channel, so that the whole acquisition systems consists of four channels, all working simultaneously, and the scaler.

2.3 Alignment and angular calibration

The alignment of all the movable elements on the beamline, i.e., the pinhole aperture, the slits, the normalisation detector, and the central detector, is done by using a telescope previously aligned with the exit of the accelerator.

A 520 nm laser is employed to perform the angular calibration. The laser is set right after the pinhole aperture and goes through all the slits. When the target plate is down, the laser reaches the central detector in the detector chamber. In this way, the zero of the beamline, corresponding to 𝜃 = 0, can be established. This measurement gives also the vertical offset on the linear manipulator of the central/lateral detectors.

To calibrate the incident and scattering angles, we use the property of the mirror target to reflect optical light. Hence, we raise the target plate, using its own manipulator, until the light is blocked. Then, we raise the central detector till the laser beam is detected again. Assuming a specular reflection, the angle subtended by the height h of the manipulator will be ζ = 𝜃 + 𝜃0 = 2𝜃0, so that the incident angle can be computed as:

$$ \theta_{0}=\frac{\zeta}{2} $$
(1)

This operation is repeated several time, so that we end up with different angles corresponding to different readings on the linear manipulator of the target plate. The incident angle can then be determined with a simple linear interpolation.

3 Efficiency definition and normalisation measurement

In the laboratory system of reference, the scattering efficiency per unit solid angle can be defined as:

$$ \eta(\theta_{0},\theta,\phi)=\frac{1}{{{\varOmega}}(\theta,\phi)}\frac{{{\varPhi}}_{\text{scat}}(\theta_{0},\theta,\phi)}{{{\varPhi}}_{\text{inc}}} $$
(2)

where 𝜃0 is the incident angle, 𝜃 and ϕ are the polar and azimuthal scattering angles (see Fig. 6), Φscat and Φinc are the scattered and incident proton count rates, and Ω(𝜃,ϕ) is the solid angle subtended by the detector.

Fig. 6
figure 6

Geometrical scheme of the incident plane. The proton beam hits the mirror sample, in the plane xz with an incident angle 𝜃0 and is pseudo-reflected with a polar scattering angle 𝜃 and an azimuthal scattering angle ϕ

The count rate of the scattered particles is given by the number of protons Nscat scattered by the SPO sample reaching the detectors in the detector chamber divided by the integration time Δtscat. In a similar way, the count rate of the incident particles is given by the number of particles Ninc intercepted by the normalisation detector in front of the mirror chamber divided by the integration time Δtinc. The number of counts of incident and scattered protons, Ninc and Nscat, is obtained by integrating the ADC histograms. This number must be corrected for the dead-time of the ADC (cfr. Section 2.2), so that the effective count rates can be expressed as:

$$ {{\varPhi}}_{\text{scat}}(\theta_{0},\theta,\phi)=\alpha \frac{N_{\text{scat}}(\theta_{0},\theta,\phi)}{{{\varDelta}} t_{\text{scat}}} \quad,\quad\quad {{\varPhi}}_{\text{inc}}=\alpha \frac{N_{\text{inc}}}{{{\varDelta}} t_{\text{inc}}} $$
(3)

with the correction factor α given by:

$$ \alpha=\frac{N_{\text{scaler}}}{(N_{\text{pulser}})_{\text{ADC}}} $$
(4)

where Nscaler is the number of counts from the pulse generator as read out from the scaler fed to the CAMAC controller module and (Npulser)ADC is the number of pulses from the pulse generator as read out from the ADC (see Fig. 5). For an ideal incoming proton beam, the number of incident particles Ninc is constant in time. However, the beam exiting the Van de Graaff accelerator was not stable, with fluctuations in the direction of the beamline varying in a time range from a few to several tens of minutes. This made necessary to take normalisation measurements before and after each scattering measurement and average them for each scattering data point, so that:

$$ \frac{N_{\text{inc}}}{{{\varDelta}} t_{\text{inc}}}=\frac{1}{2}\left( \frac{N_{\mathrm{inc,1}}}{{{\varDelta}} t_{\mathrm{inc,1}}}+\frac{N_{\mathrm{inc,2}}}{{{\varDelta}} t_{\mathrm{inc,2}}}\right) $$
(5)

where Ninc,1 and Ninc,2 are the counts in two consecutive normalisation measurements with integration times Δtinc,1 and Δtinc,2, respectively.

Concerning the uncertainties, the one on the scattering angle is given mainly by the errors on the angular calibration, the detector aperture, and the indeterminate position of the impact point of the beam on the mirror surface. The uncertainty on the incident angle 𝜃0 is dominated by the dimension of the aperture on the central detector and by the length of the target. It resulted in \(\sim {0.1^{\circ }}\) for all the chosen scattering angles. Lastly, the uncertainty on the scattering efficiency is mainly given by the intrinsic fluctuation of the proton beam. Minor contributions are due to the count statistics and to the error on the solid angle Ω(𝜃,ϕ). The sum of this contributions results in statistical fluctuations of ± 20% on the scattering efficiencies.

4 Results and discussion

We measured the scattering efficiency at two different energies (high- and low-energy data sets, hereafter) and at four different incident angles: 0.6, 0.8, 1.0, and 1.2, both on-axis and off-axis, the latter at an angle ϕ of about 2. Each data set consists of on-axis and off-axis scattering efficiencies. Results are shown in Figs. 7 and 8, where the scattering efficiencies have been multiplied by the square of the incident angle (as in [1]) and are displayed as a function of the scattering angle divided by the incident one, i.e., Ψ = 𝜃/𝜃0.

Fig. 7
figure 7

Scattering efficiencies of the low-energy data set as a function of the scattering angle Ψ = 𝜃/𝜃0, for the incident angels of 0.6, 0.8, 1.0, and 1.2, for the on-axis (top panel) and off-axis (bottom panel) configurations. Energy beam of 172 ± 30 keV

Fig. 8
figure 8

Scattering efficiencies of the high-energy data set as a function of the scattering angle Ψ, for the incident angels of 0.6, 0.8, 1.0, and 1.2, for the on-axis (top panel) and off-axis (bottom panel) configurations. Energy beam of 471 ± 25 keV

For the high-energy data set (Fig. 8) we used a beam at \(\sim \)590 keV from the accelerator, which was degraded by the Al foil down to 471 ± 25 keV. This energy value was chosen mainly for purposes of comparison with the previous measurements on the eROSITA mirror sample ([2, 3] cfr. Fig. 9). The rationale behind the low-energy value can be found in the work of [7], who showed that the highest transmission efficiency of soft protons is observed for those protons impacting the mirrors with 40-60 keV, for both instruments on board of Athena. Hence, it is crucial to investigate the scattering of soft protons at energies around and below 100 keV. Unfortunately, the present setup could not reach such low energies, limiting us to a proton beam with an energy of \(\sim \)340 keV at the exit of the accelerator, degraded to 172 ± 30 keV by the Al foil. In both cases, the values of the incident energies were determined by simulations with the software TRIMFootnote 6 (TRansport of Ions in Matter, [14]), already validated in [2].

Fig. 9
figure 9

Comparison of the eROSITA scattering efficiencies (blue dots) with the SPO ones (green dots for the low-energy set and red dots for the high-energy set), for the on-axis (top panels) and off-axis (bottom panels) data

As expected, the on-axis scattering efficiencies peak at the specular angle (\(\psi \simeq 1\)) and are consistent with each other within the uncertainties. However, a higher spread is observed for the high-energy on-axis data set (Fig. 8, top panel), with efficiencies ranging from 0.03 to 0.07 at the peak of the distribution. Also the off-axis data show a significant spread, which is expected in this case.

Overall, the maximum scattering efficiency values are \(\sim \)0.07 and \(\sim \)0.02 for the on-axis and off-axis configurations, respectively, with the low-energy data set showing slightly smaller efficiencies than the high-energy one.

4.1 Comparison with the eROSITA measurements

Figure 9 shows the eROSITA measurements [2, 3] overlapped to the SPO data, for both the energies and the ox-axis and off-axis configurations. Though the SPO efficiencies are systematically higher than the eROSITA data, they are consistent within the error bars.

Due to this consistency, we applied the semi-empirical model proposed in [1] to the scattering efficiencies of SPO, as shown in Fig. 10.

Fig. 10
figure 10

Comparison between the experimental scattering efficiency of SPO (points) and the semi-empirical model developed from eROSITA data (solid line), for the low-energy (green) and high-energy (red) data sets and for the on-axis (top panels) and off-axis (bottom panels) measurements

Overall, the model well reproduces the scattering efficiency of the low-energy data set, but overestimates the efficiency of the high-energy data set by a factor of \(\sim \)1.5 times. Nonetheless, it has to be borne in mind that at this stage we simply overlapped the semi-empirical model developed for the eROSITA to the new SPO experimental data. A more accurate model, specific for SPO, can be obtained by fitting the data with the formula of [11] in non-elastic approximation, as in [1], provided that energy loss measurements are retrieved from the raw data.

Lastly, we group the efficiency values from the two data sets by the incident angle, irrespective of the energy of the incident beam. Figure 11 shows that the data are perfectly consistent with each other and with the old eROSITA measurements when grouped by the incident angle, without accounting for the energy. Once again, the semi-empirical model derived for the eROSITA data is overlapped with the data, resulting in a general acceptable agreement.

Fig. 11
figure 11

Scattering efficiencies from eROSITA mirror and the SPO samples, grouped by different incident angles, irrespective of the energy of the impinging beam. The solid line stand for the semi-empirical scattering model, with the error on the resulting efficiency given by the green area

5 Conclusions

Within the EXACRAD project, we measured for the first time the scattering efficiency of a single wafer of SPO hit by low-energy protons at grazing incidence. Measurements were performed at two different energies, of about 470 keV and 170 keV, and at four different incident angles, 0.6, 0.8, 1.0, and 1.2, both on-axis and at an off-axis angle of about 2.

Hereafter some major remarks:

  • the scattering efficiencies show the trend expected from [11] and from the experimental data on eROSITA [2, 3]. The on-axis data peak close to the specular reflection, while the off-axis data show a peak shifted to higher Ψ; the off-axis data reach lower efficiencies than the on-axis ones; higher incident angles resulted in higher scattering efficiency;

  • the SPO data are generally consisted with the eROSITA data, though the high-energy data set show a higher spread in efficiency;

  • as for the eROSITA data, the scattering efficiency very weakly depends on the energy of the incident beam;

  • the semi-empirical model developed from eROSITA experimental data is able to acceptably reproduce the low-energy data set, while it results in higher efficiencies for the high-energy data set. The same model can be improved specifically for SPO, with a direct fit of the experimental data, provided that energy loss measurements are retrieved from the raw data.

As stated above, the experimental configuration used for this experiment is not suited for measurements at energies below 100 keV, which are those expected to contribute the most to the background of Athena [8]. Therefore, some possible changes to the setup could be explored in the future to enable measurements at lower energy ranges.

The work here presented is only the first step towards a thorough estimation of the SP flux expected at the focal plane of Athena. Indeed, after being scattered by the optics, SP cross all the other elements along their path, and, in particular, the filters located in front of each detector. As showed in [8], the filters not only reduce the energy of the protons, but also act as degraders, altering their original trajectories. Hence, a thorough prediction of the overall transmission energy of SP can be achieved only combining the processes of scattering from the optics, crossing of the filters and energy release within the detectors. The latter two phenomena can be investigated, for instance, by simulations, as done, e.g., by [4]. If the estimated soft proton flux should still result higher than the scientific requirement of 5 × 10− 4 cts s− 1 cm− 2 keV− 1, in the 2–10 keV energy band, for 90 % of the observing time (cfr. Section 1), then appropriate solution should be adopted. Currently, the hypothesis of a magnetic diverter specific for protons is under discussion and possible designs are under investigation.