Introduction

The High Energy Physics (HEP) community is currently discussing the possible upgrades of today’s Large Hadron Collider (LHC), in order to schedule a roadmap for the construction of the next generation HEP experiment. Along with the Compact Linear Collider (CLIC1) and the High Energy Large Hadron Collider (HE-LHC2), the Future Circular Collider (FCC) is the most ambitious of the proposed machines, with its 100 km long tunnel, being designed to reach hadron-hadron collisions at 100 TeV, about 8 times higher than in today’s LHC3. Such new synchrotron ring will allow to expand the current energy and luminosity frontiers, pushing towards novel discoveries in the particle physics field. Along with the civil engineering challenges of building a 100 km long tunnel, several other technological challenges are being considered in order to be able to build and operate such a poweFul machine. Among these challenges, the very harsh radiation environment, evaluated to be several orders of magnitude higher than in the current LHC, will require new dedicated technologies for the detectors and electronics.

The impact of cumulative radiation damage to materials (such as semiconductors) is commonly measured in term of Total Ionizing Dose (TID) and Displacement Damage (DD). The TID is proportional to the ionizing energy released in the material and is expressed in J/Kg or Gray [Gy]. The DD is instead proportional to the number of particles crossing the material and is expressed in particles/cm2 [Ф]. For semiconductor material, the normalized 1 MeV neutron equivalent fluence [n1MeV/cm2] is instead commonly used.

In the LHC, the total Dose and Fluence, cumulated over 10 years of operation, are expected to be about 100 kGy and >1015 particles/cm2 respectively (in the LHC experiments). In order to keep these radiation levels under constant monitoring, the Radiation Monitoring (RADMON) dosimeters are employed in different sections of the tunnel and its experiments (PH-RADMON, targeting the needs of LHC experiments4 and LHC-RADMON for the tunnel area5).

In the case of the FCC, due to an increased collision energy (from 14 TeV to 100 TeV) and rate (from 0.85 GHz to 32.4 GHz)3, radiation levels are expected to be several orders of magnitude higher than today’s LHC, exceeding several tens of MGy (with >1017 particles/cm2) inside the FCC experiments6, and tens of kGy (with >1015 particles/cm2) in certain sections of the FCC tunnel itself7.

Figure 1 shows the integration ranges of the RADMON as well as other types of dosimeters that are today commonly used. Special focus is made on the “MGy range” expected for the FCC experiments, where none of these technologies are capable of properly function, because of saturation problems and/or extensive and irreversible radiation damage.

Figure 1
figure 1

Different active and passive technologies for TID dosimetry available on the market are shown with respect to their monitoring range60. The highlighted and empty area in the Ultra High Dose range indicates the target area of the technology presented in this paper.

These noticeable limitations have triggered an extensive research for new technologies8, and while for the ultra-high fluence monitoring devices based on metallic thin-films are currently under research at CERN and are showing promising results9, no solutions have yet been proposed for an ultra-high dose monitoring technology that could be coupled to the PH-RADMON extending its range of detection across and beyond the MGy range.

In the present work, we report the first proof of principle of an alternative technology for the ultra-high dose monitoring based on the LOF technology10,11,12. The main characteristics, working principle, and target operational range of the work here presented are summarized in Table 1 and compared with the dosimeters found in the PH-RADMON.

Table 1 Operational range of current PH-RADMON devices for the LHC experiment4 as compared with the performance required (after 10 years of operation) from the new Ultra High Fluence and Dose dosimetry for the future HL-LHC upgrade and FCC4,7,60.

This paper is structured as follows: Section “Optical fiber dosimeters” gives an overview of the existing OF dosimeters; Section “Lab on Fiber technology: a promising approach for MGy dosimetry” provides a specific description of the proposed innovative dosimetry technology including the LOF probe design, fabrication and characterization, as well as numerical simulations, based on finite element analysis, predicting the effects of perturbations of different parameters on the LOF spectral response; Section “Irradiation Tests at CERN” deals with the irradiation tests, providing details about the CERN IRRAD facility and the experimental results recorded with LOF probes; Section “Data analysis and discussion” focuses on the processing of the experimental data and comparison of LOF performance versus the numerical simulation and literature; Section “Conclusions and future works” gives our conclusions concerning the first application of LOF technology for dosimetry at ultra-high dose levels and a view on the planned future works to further investigate these innovative dosimetry technique.

Optical fiber dosimeters

In order to overcome the impossibility to monitor ultra-high dose levels with the current technologies, in the last years a great effort has been made to demonstrate the applicability of OFs in this field13,14. The inherent properties of OFs, including the small size, low mass, immunity to electromagnetic interference, possibility of multiplexing and real-time, continuous and remote measurement15, represent unrivaled advantages for OF-based dosimeters16. According to the role played by the fiber within the measurement system, OF dosimeters can be classified, into intrinsic and extrinsic architectures. In intrinsic configurations, where the OF acts both as the sensing and the light-guiding components16, the main effects exploited are the radiation-induced attenuation (RIA)17,18,19 and the radiation-induced changes in the refractive index (RI) of the fiber20, which also led to the application of Fiber Bragg Gratings (FBGs) and Long Period Gratings (LPGs), notoriously sensitive to changes in the fiber RI, as dosimeters21,22,23,24,25,26. Differently, in extrinsic OF dosimeters, the OF serves only as a guide of the light travelling to and from a material sensitive to the radiation. These are the majority and, among their basic exploited principles, they can be recalled the scintillating materials, thermo-luminescence and optically stimulated luminescence27,28,29,30,31,32,33,34,35. It is worth noting that an important contribution to the development of OF dosimeters has been provided by the use of poly(methyl methacrylate) (PMMA) as sensing element13,16.

The performance of OF dosimeters reported so far, including dose sensitivity and operation ranges, is strongly influenced by the fiber chemical composition, radiation environment and manufacturing process36. Indeed, in the case of FBGs written in pure-silica fibers, which exhibit no change of the Bragg wavelength, as well as of the amplitude and the shape of the spectra, at a dose of 100 kGy25, it has been demonstrated that the dose sensitivity and operation range can be significantly modified by doping the fibers, achieving limits of saturation even higher than 200 kGy22,37. Resonance wavelength shifts definitely larger than those measured in FBGs have been observed in turn around points (TAP) LPGs (up to 65 kGy24) and chiral LPGs (up to 100 kGy without complete saturation23). In this context, a clear definition of the limits of operation is not yet possible, as the experiments reported so far are few, as well as the exploited dose rates scenario, which can have a major influence on the response of the sensor and its saturation point. A notable drawback for these dosimeters is the remarkable temperature sensitivity, especially in TAP-LPGs. Moreover, a large number of these approaches, especially those based on the measurement of the optical signal intensity, have the drawback of requiring a highly stable set-up, which is not always available in a hostile radiation environment13,16.

In order to meet the abovementioned needs and to overcome the limitations suffered by the majority of OF-based approaches, especially in terms of operational range and saturation levels, as well as influence of temperature and issues related to intensity-based measurements, in this work, we propose a new extrinsic OF dosimeter for ultra-high dose radiation, based on the so-called LOF technology10,11,12.

The key idea at the basis of LOF technology is to transform a ‘simple’ OF into a miniaturized and multifunctional platform for specific applications through the integration of functionalized materials and components defined at micro or nanoscale10,12,38,39. After devoting a great effort to solve manufacturing issues, our attention has been focused on the development of LOF prototypes for specific applications, mainly (but not limited to) chemical and biological sensing12,39,40.

Here, the proposed LOF platform consists in a metallo-dielectric resonator realized on the OF tip (OFT)38,39 including a patterned PMMA layer covered by a thin gold overlay. Our expectation is to reduce saturation effects taking advantage from the large set of radiation sensitive materials involved in, while wavelength encoding detection allows robust measurements against intensity fluctuations along the whole optical chain.

Preliminary experimental results reported here demonstrated indeed the effectiveness of the proposed OF device as potential dosimeter at MGy dose levels, paving the way for its future application in HEP experiments.

Lab on Fiber technology: a promising approach for MGy dosimetry

Among the different configurations embraced by the LOF concept, the OFT represents an attractive platform for remote sensing applications, thanks to its inherently light coupled micro-sized active area, joined with the faculty of integrating complex micro- and nano-metric structures made of dielectric and metallic materials supporting resonant light trapping effects. The principle of operation of this kind of devices is based on the sensitivity of such resonant phenomena to the optical and geometrical variations induced from the surrounding.

Based on a previous demonstration of LOF device carried out by the Unisannio group38, the present work deals with the realization of a metallo-dielectric resonator on the OFT sensitive to radiation. In particular, the LOF structure can be regarded as the superimposition of two metallic gratings separated by a dielectric layer (Fig. 2).

Figure 2
figure 2

Cross section view of the LOF resonator realized on the OFT.

By illuminating the metal-dielectric structure in ‘out-of-plane’ configuration, localized surface plasmon resonances may be excited by the phase-matching conditions between the scattered components and the modes supported by the hybrid structure38. The spectral response arising from these resonant phenomena is sensitive to variations of both optical properties and size of the platform components. These include real and imaginary parts of the RI of adopted materials, as well as lattice period, holes radius, metal and dielectric thickness. The dielectric layer, in particular, is one of the constituent elements able to play a fundamental role in the design and tailoring of the spectral features of a LOF device for dosimetry application. Among dielectrics, polymers are recognized as highly sensitive to radiations41. Within this class of materials, the attention has been focused here on the PMMA because of the well-known effects induced by proton beams on its optical and geometrical properties42,43,44,45 and since it exhibits a RI close to 1.5, i.e. higher than that of the OF, creating a condition that promotes the excitation of spectral resonances38.

LOF probe design

The LOF probe design is based on the use of the finite element method provided by the Comsol Multiphysics – RF module46. The design parameters, including a holes radius of 240 nm, grating period of 1050 nm, PMMA and gold thickness of 200 nm and 50 nm, respectively (Fig. 3a), have been chosen in order to set the resonance wavelength in the operating range of standard single mode fibers, enabling the use of robust and reliable commercially available optoelectronic interrogation units. As shown elsewhere47,48, the exploitation of the crystal translational symmetries can reduce the computational domain to the unit cell (given by one period of the structure), terminated with Floquet boundary conditions, placed two-by-two in the opposite walls. The resulting structure supports a transverse-electromagnetic wave emulating the normally incident plane-wave. This kind of excitation was preferred to a proper fiber mode with Gaussian profile to simplify the numerical simulations38. The RI data used for modeling gold in the IR region were taken from ref.49, for silica from ref.50, while for PMMA from ref.51. The numerical reflected spectrum of the LOF structure in a normalized scale exhibits a high reflectance value (higher than 0.6) interrupted by a resonance dip centered at λres1566 nm with a Q-factor of 47 (Fig. 3b).

Figure 3
figure 3

(a) Schematic of the LOF nanostructure (unit cell) and electric field distribution in correspondence of the resonance wavelength (λres1566 nm); (b) numerical reflected spectrum.

The exposure to proton beam could involve the modification of one or more parameters of the LOF structure, affecting both its optical properties (variations of the real or imaginary part of the materials) and geometry (such as the thickness of the layers or the dimensions of the holes and their period).

Depending on the class of materials used to build each single component of the structure, the radiation effects on the LOF device can be essentially divided into three main groups: effects occurring on the polymer, on the metal and on the fiber itself (see Table 2).

Table 2 Resume of the main radiation effects reported in literature occurring in the classes of materials used to build the LOF device.

With regards to the polymeric layer, it is well known that proton radiation causes several changes in the PMMA, including an increase of the propagation losses42,44 and phenomena of compaction and densification42,43,44,45. The consequences on the PMMA of these phenomena in terms of optical properties are an increase of the imaginary part of the RI (kPMMA), due to the enhanced losses, and an increase of the real part (nPMMA) due to the densification. Moreover, compaction may modify the geometry of the PMMA layer along both the vertical direction, implying a reduction in thickness (tPMMA), and the horizontal one, causing an enlargement of the holes radius (rholes).

Concerning the metallic layer, experimental results recently reported elsewhere8 showed that thin layers of certain metals, specifically copper, with 500 and 1000 nm thickness, when exposed to proton beam exhibit a positive variation of resistivity with increasing particle fluence due to the accumulated displacement damage.

A variation in resistivity, as it implies a reduction in the concentration of free charge carriers, causes also a change in the real (ngold) and imaginary part (kgold) of the RI of the material52. In particular, an increase of ngold and a decrease of kgold are expected53. Irradiations of 20 µm layer of gold with fast neutron up to 2.2 × 1019 n/cm2 showed an increase of the gold resistance with the irradiation dose, at fixed temperature54. Nevertheless, a deeper understanding of the effects of very high particle fluence, in the range of interest of future HEP experiments, on nanometric metal films including the physical mechanisms leading to an increase of the resistivity is necessary.

Actually, it must be said that the interaction with radiations leads to a number of chemical and/or physical changes also in the fiber glass, including the creation of preferential absorption in a broad spectral range, well-known as RIA36. Indeed, the enhancement of the imaginary part of the silica RI (kfiber) responsible of the absorption is also associated to an increase of the real part of the RI (known as radiation-induced RI variation), according to the Kramers-Kronig relations25. Such phenomena, as stated before, are the basis of intrinsic OF dosimeters.

Once identified the effects that ionizing radiations may induce on the optical and geometric parameters of the LOF device, numerical simulations have been carried out using the aforementioned RF module to evaluate the consequences produced by the variation of each single parameter of the structure on the reflected spectrum, keeping the others fixed.

In particular, the results of the simulations of a perturbation of nPMMA, performed with steps of 10−2, showed a linear dependence of the resonance wavelength from the changes of nPMMA with a positive slope of 310 nm/refractive index unit (RIU) (Fig. 4a). Furthermore, the evolution of the reflected spectrum occurring when nPMMA increases (Fig. 4b) shows a slight decrease of the reflectivity outside the resonant dip (hereafter called baseline).

Figure 4
figure 4

(a) Numerical simulations predicting the effects of a perturbation of nPMMA on the LOF resonance wavelength. Black dashed lines refer to nominal value of nPMMA. The inset box indicates the equation of the linear fitting. (b) Whole spectrum evolution determined by the variation of nPMMA. The black arrow indicates the variation direction of the parameter.

Differently, a perturbation of kPMMA has no consequences on the resonance wavelength until a value of 10–2, beyond which a fast red shift following a quadratic behavior is expected (Fig. 5a). Such large perturbations of kPMMA also cause a strong deformation of the resonant dip, clearly visible in Fig. 5b.

Figure 5
figure 5

(a) Numerical simulations predicting the effects of a perturbation of kPMMA on the resonance wavelength. Logarithmic scale is used for the x-axis. The inset box indicates the equation of the quadratic fitting. (b) Whole spectrum evolution determined by the variation of kPMMA. The black arrow indicates the variation direction of the parameter.

As regards the thickness of the PMMA layer tPMMA, its change leads to an increasing quadratic variation of the resonance wavelength (Fig. 6a), while the baseline is a decreasing tPMMA function, as shown in Fig. 6b.

Figure 6
figure 6

(a) Numerical simulations predicting the effects of a variation of tPMMA(with steps of 10 nm) on the resonance wavelength. Black dashed lines refer to nominal value of tPMMA. The inset box indicates the equation of the quadratic fitting. (b) Whole spectrum evolution determined by the variation of tPMMA. The black arrow indicates the variation direction of the parameter.

The effect of the other parameters of the structure influenced by radiation on the LOF platform spectrum predicted by numerical simulations are reported in detail in the Supplementary Information (Section S1).

Table 3 summarizes the radiation induced changes in the properties of the different materials constituting the LOF platform according to previous literature and includes the sensitivity analysis to these properties expressed in terms of resonant wavelength and reflectance baseline (far from the resonance) provided by our numerical investigation. Both the sign and the range of variation are shown, where present. When the dependence of the spectral observable from a given parameter is not linear, the sensitivity has been calculated by linearizing the characteristic curve around the nominal value of the parameter itself. Since the values reported in the literature refer to different types of radiation, energy levels, doses, fluences, operational environments spectral ranges, it is clear that the reported table is purely indicative.

Table 3 Sensitivity analysis. The signs reported in the second column follow the legend: “ + “ for a positive change, “−“ for a negative change. The reflectance baseline sensitivity is evaluated on the normalized reflectance in the range 0–1.

Probe Fabrication

The realization of the 2D hybrid metallo-dielectric nanostructure consists of three main steps38, as schematically explained in Fig. 7a: i) polymeric overlay deposition onto the OFT, ii) OFT exposure to electron beam and polymer development, iii) gold superstrate deposition.

Figure 7
figure 7

(a) Schematic representation of the technological steps used in the probe realization process: i) polymeric overlay deposition, ii) OFT exposure to electron beam and polymer development, iii) gold superstrate deposition. (b) Optical microscope image of an OFT after the deposition of the PMMA layer, showing a uniform region with a diameter of about 50 μm around the fiber core. (c) SEM image showing the top view of the metallo-dielectric nanostructure realized on the OFT.

The uniform PMMA deposition on the OFT has been achieved by using a customized spin coater involving a rotating plate modified in order to host the OF. For the purpose, a standard single mode OF (Corning SMF-28) was adopted. After the OF housing in the drilled plate, the PMMA solution was spun upon it. As shown in Fig. 7b, the uniformity of the PMMA layer has been guaranteed in the central part of the OFT for a region of about 50 μm, definitely wider than the fiber core diameter. Furthermore, a careful combination between the spin-coater plate velocity (2000 rpm) and the rotation time (1 minute) allowed obtaining a PMMA thickness of about 200 nm, as dictated by the design.

After the PMMA overlay deposition onto the OFT, a periodic pattern has been obtained through Electron Beam Lithography (EBL). Following the design, a periodic structure consisting of a holes array with a period of ~1050 nm and a radius of ~240 nm has been realized in the PMMA layer using the RAITH 150 EBL system. The EBL process involved an electron acceleration voltage of 20 kV, a dose of 120 μC/cm2 and an electron beam aperture of 7.5 μm. After the PMMA electronic impression, the desired 2-dimensional (2D) array has been obtained by developing the polymeric layer.

The last realization step relies on the deposition of a 50 nm-thick gold superstrate onto the polymeric pattern realized by means of a conventional DC magnetron sputtering. To ensure the good uniformity of deposition, a dedicated holder has been designed and realized to keep the OFT perpendicular to the gold target.

Morphological and optical characterization

After manufacturing, two prototypes of the proposed LOF platform (henceforth referred to as prototypes A and B) have been fabricated on two different fibers. A morphological characterization of the devices has been carried out by acquiring several scanning electron microscope (SEM) images of the structures (Fig. 7c).

With the aim of evaluating the spectrum reflected by the LOF probes, an optoelectronic set-up comprising a broadband light source (with bandwidth 1200–1700 nm), a directional coupler and an optical spectrum analyzer (OSA), ANDO AQ6317C, with best wavelength resolution of 10 pm and dynamic range 60 dB, was adopted. The OSA was connected to a computer and controlled by a LabVIEWTM plug-in. In order to compensate for the intensity fluctuations which may occur during the characterization, the LOF reflectance spectrum was calculated by normalizing the reflected spectrum with that reflected by a fiber-optic reference mirror, fabricated by depositing a 160 nm-thick gold film on the tip of a standard single-mode fiber. The experimental reflectance spectra of the LOF prototypes are shown in Fig. 8. Comparing the experimental and simulated spectra (Fig. 3b), slight differences in the exact location of the resonance wavelength and visibility can be noted. These discrepancies can be attributed to the fabrication tolerances influencing both the holes radius and period. Nonetheless, both LOF devices exhibited a resonance dip centered around 1550 nm, enabling the use of a compact optoelectronic interrogation unit to monitor the spectrum during the proton irradiation.

Figure 8
figure 8

Experimental normalized reflectance spectrum of the 2D hybrid metallo-dielectric structure: (a) prototype A; (b) prototype B.

Irradiation Tests at CERN

To test the response of the sensor to ultra-high levels of radiation, as the ones encountered in the future generation of particle accelerators, the first irradiation experiments have been performed at the high-energy proton irradiation facility (IRRAD), located at the East experimental area of the Proton Synchrotron (PS) accelerator at CERN55.

The PS delivers to IRRAD a Gaussian proton beam of 23 GeV with a typical size of 12 × 12 mm2, in spills of ~5 × 1011 p/spill, every 10 seconds on average. With such configuration, samples typically of 10 × 10 mm2 in size can be exposed to a fluence of 1016 p/cm2 in about 2 weeks, while smaller objects as the OF sensors described in this work, over the same period of time, can be irradiated up to higher particle fluence approaching the 1017 p/cm2. Concerning the samples, these are placed on remotely controlled tables that can be moved along the X-, Y-axis in order to precisely align them with the beam as shown in Fig. 9a. The tables are located in three different irradiation zones according to the nature of the samples to be irradiated. Due to the intrinsic “low-mass” of the OF samples, two prototypes of the proposed LOF dosimeter have been installed upstream the IRRAD facility on the table IRRAD3 in Zone 156.

Figure 9
figure 9

(a) OF samples placed on the irradiation table IRRAD3; (b,c) detail of the OF sensor mounted on the sample holders.

The intensity of the extracted proton spills is monitored using a Secondary Emission Chamber (SEC) device provided by PS beam instrumentation, while the profile of the proton beam, with a spill-by-spill resolution, is obtained by a custom-made Beam Profile Monitor (BPM)57. The calibration of this beam instrumentation, as well as the final total number of protons that crossed the sample (normalized per cm2) during the irradiation test, is obtained by means of activation measurements of thin Aluminium (Al) foils accurately inserted in front of the irradiated sample58.

The basic measurement performed on the 23 GeV proton beam is the determination of the proton fluence by evaluating the 24Na and 22Na activity in thin Al foils produced via dedicated nuclear reactions59. With these activation techniques it is possible to obtain fluence measurements with an accuracy of ±7%.

Taking into account that ionization is the main contribution to the energy loss of a charged particle and that its mean value, the stopping power (dE/dx), is given by the Bethe–Bloch law, it is possible to simply convert the fluence into the dose (Gy) deposited in thin samples, using the following formula:

$$D=K\cdot {(dE/dx)}_{m}\cdot \phi $$

where φ is the proton fluence expressed in particles/cm2, K = 1.602·10−10 is the scale factor and (dE/dx)m expressed in MeV·cm2/g is the minimum ionizing energy loss rate. For 23 GeV protons it has value 1.664 MeV·cm2/g for Silicon (Si). Si was chosen as reference material for the absorbed dose conversion from measured particle fluence, since standard dosimetry for electronics is commonly calculated in Si. This choice is also justified by the fact that the energy loss of particles at energies of 23 GeV is about the same for any material60 (e.g.61 Al 1.615 MeV·cm2/g, Li 1.639 MeV·cm2/g).

LOF prototypes irradiation and experimental results

The LOF sensors have been exposed to the 23 GeV proton beam for 9 days, with an average dose-rate of about 200 kGy/day and a total cumulated dose of about 1.8 MGy, corresponding to a total fluence of 0.67 × 1016 protons/cm2 (Fig. 10). The dose was calculated from the protons fluence evaluated through Al-foils activation measurements and normalized to 1 cm2. The Al-foil used for dosimetry, placed upstream the LOF sensor prototypes, is also visible in Fig. 9b,c. During the exposure, relative humidity and temperature have been monitored by using the instrumentation available at IRRAD62. The details on the monitoring of these parameters are shown in the Supplementary Information (Section S2).

Figure 10
figure 10

Daily and cumulative dose absorbed by the LOF devices during the exposure to proton beam.

LOF prototypes were connected to a sm125 Optical Sensing Interrogator (Micron Optics Inc.) covering the wavelength range 1510–1590 nm with an accuracy of 1 pm in order to measure the effects of the proton beam on their spectral response. Spectra reflected by the LOF prototypes were saved by using a time interval of 1 hour. Because of the significant amount of data collected during the 9-days exposure and to better appreciate the spectral variations undergone by the two probes, graphics shown in Fig. 11 report only one spectrum per day, related to the dose accumulated in the previous 24 hours.

Figure 11
figure 11

Spectra reflected by the two LOF prototypes during the exposure to 23 GeV proton beam for 9 days. The legend indicates the dose accumulated at the end of each day; (a) prototype A; (b) prototype B.

The first relevant result arising from observing the spectra shown in Fig. 11 is that the exposure to protons at high fluence levels did not produce a strong degradation of the spectrum, demonstrating the “radiation resistance” feature of the developed device. Moreover, despite the tolerances affecting the realization processes may cause slightly dissimilar spectral features, the two prototypes showed a very similar behavior, consisting of a resonance migration towards lower wavelengths.

Based on this observation, once acquired, LOF spectra were processed by a Matlab® script providing filtering and tracking of the changes occurred in the resonance wavelength. Specifically, in order to take into account also possible spectral deformations and asymmetries (mainly affecting the dips shown in Fig. 11b), an accurate approach based on the determination of the barycenter within an adequate distance (0.3 dBm) from the minimum was used. The average of the wavelength shifts undergone by the two probes with respect to the corresponding starting spectra is reported as a function of the time in Fig. 12a. Furthermore, in the same picture, the cumulative absorbed dose is indicated on a daily basis. The slightly irregular shape of the blue (negative) shift versus time could be explained by considering that the dose-rate was not constant during the exposure, as reported in Fig. 10. The calibration curve reporting the resonance wavelength shift sampled at the end of each day as a function of the cumulative absorbed dose is shown in Fig. 12b.

Figure 12
figure 12

Average of the resonance wavelength shifts calculated for both LOF sensors; (a) wavelength shift reported as a function of the time (1 point per hour); (b) calibration curve displaying the wavelength shift reported as a function of the dose (1 point per day). Error bars refer to the standard deviation calculated on the two prototypes.

The average blue shift in correspondence of a total dose of 1.8 MGy was of 1.4 nm, with a sensitivity of about −0.6 nm/MGy. The slope of the calibration curve suggests that the present device could detect doses even higher than those reported here, also in view of the radiation resistance characteristic of the OF probe. In addition, the very low amplitude of the error bars, referring to the standard deviation between the responses of the two prototypes, demonstrates the high coherence exhibited by the two devices.

In order to estimate the accuracy of the proposed dosimeter, a steady state standard deviation pertaining to the resonant wavelength of 0.010 nm was calculated in time interval of 10 hours in constant environmental conditions and no radiation. By dividing this value for the sensitivity to the dose, a resolution of 16.7 kGy was achieved, corresponding to a resolution normalized to the maximum dose observed of 1.8 MGy equal to 0.9%.

Another effect resulting from the exposure to the proton beam is the progressive rise of the baseline, clearly visible from spectra shown in Fig. 11. In particular, an average variation of 0.16 dBm has been observed as response to a total dose of 1.8 MGy.

Data analysis and discussion

Experimental results reported in the previous section have shown that the main effects of the LOF probes exposure to 23 GeV proton beam on the reflected spectrum are a blue shift of the resonance wavelength of about 1.4 nm (Fig. 12) and a slight increase of the baseline of about 0.16 dBm (Fig. 11). According to the behaviors listed in Table 3, which matches the phenomena reported in the literature with the numerical predictions, there are three parameters able to induce a blue shift as response to radiations: tPMMA, ngold and kfiber. Among these, radiation changes in kfiber and ngold are expected to induce a decrease in the reflectance baseline in contrast with the experimental evidence of a baseline increase (see Table 3). On the contrary, a reduction of tPMMA is consistent with both the trends of the resonance wavelength and of the baseline (Table 3). The only other phenomenon able to justify the experimental increase of the baseline is an enlargement of the rholes. Like the decrease of tPMMA, this too can be considered as a consequence of the compaction of the PMMA layer. However, an extension of rholes is not well-matched with the experimental negative wavelength shift shown in Fig. 12. On this line of arguments, it is possible to conclude that the dominant mechanism induced by the radiation able to explain both the resonance migration to lower wavelengths and at the same time the baseline increase is the reduction of the PMMA layer thickness. Moreover, this behavior is consistent with data reported so far in the current literature42,43. Indeed, the PMMA presents linear polymer chains composed of monofunctional monomers. Therefore, high energy particles traversing a polymer produce bond breaking that creates reactive species, as radicals, promoting chain transfer reaction. Formation of new radicals increases with the deposited dose. As a consequence, it is possible that PMMA becomes saturated with reactive species and the PMMA exposed to ionizing dose exhibits a new morphology, leading to a compaction of the PMMA layer itself63.

Conclusions and future works

The HEP community is currently preparing for the incoming upgrade of the present Large Hadron Collider, into the higher performance High Luminosity LHC, and it also envisages the possible future upgrades such as the HE-LHC and the new machines like the CLIC and the FCC. In all these options the very harsh radiation environment is evaluated to be several order of magnitudes higher than in the current LHC, demanding new dedicated technologies for both the detectors and electronics to be used. These radiation levels are expected to exceed several tens of MGy and currently available technologies for dose monitoring cannot withstand and/or measure up to these radiation levels.

An innovative dosimeter, based on the LOF technology, for ultra-high dose monitoring in the challenging scenario of the FCC at CERN has been the object of the experimental research discussed in the present paper.

In particular, we proposed here a new extrinsic OF dosimeter, based on the so-called LOF technology, consisting in a metallo-dielectric resonator realized on the OF tip including a PMMA layer.

Once identified the effects that ionizing radiation induces on the optical and geometric parameters of the LOF device, numerical simulations have been carried out to evaluate the impact induced by the variation of each single parameter on the reflected spectrum, keeping the others fixed.

Two LOF probes were designed and manufactured in order to set the resonance wavelength in the operating range of standard single mode fibers, enabling the use of robust, reliable, and commercially available optoelectronic interrogation units. Technical details were presented on selection of basic materials, manufacturing techniques and LOF probe final morphological and optical characterization. The two LOF dosimeter prototypes were exposed at the IRRAD Proton Facility at CERN in Geneva with 23 GeV protons for a total deposited dose up to 1.8 MGy in 1 cm2 of Si.

The exposure to protons at high dose levels did not produce any strong degradation of the spectrum, demonstrating the “radiation resistance” feature of the LOF technology. The experimental data from both prototypes show a clear variation of the LOF signal versus the total dose. The main effect consists of a resonance migration towards lower wavelengths with increasing dose level. In particular, a blue shift of the resonance wavelength of about 1.4 nm and a slight increase of the baseline of about 0.16 dBm were recorded. From the simulations, it arises that a possible reduction of the thickness of the PMMA layer of the LOF resonator (tPMMA) is consistent with the observed spectral trends. This hypothesis is also supported by literature42,43 referring to a similar behavior of PMMA exposed at different doses, energies, and radiation sources.

The preliminary experimental results reported in this paper demonstrate the effectiveness of the proposed OF device as potential dosimeter at MGy dose levels, overcoming the limit of 100 kGy achieved by the current techniques. In particular, a maximum dose of 1.8 MGy has been measured with a sensitivity of −0.6 nm/MGy and a resolution of 16.7 kGy. Additional research will follow in order to deeper understand the radiation damage mechanisms occurring in the LOF nanostructure. New LOF prototypes are being currently designed and a second irradiation campaign is being scheduled at the IRRAD Proton Facility targeting higher dose levels, to assess the LOF dosimeters performance at the tens MGy range, and then to study the post-irradiation annealing effects.