Abstract

Targets of isotropic graphite and hexagonal boron nitride were exposed to short pulses of uranium ions with ∼1 GeV kinetic energy. The deposited power density of ∼3 MW/cm³ generates thermal stress in the samples leading to pressure waves. The velocity of the respective motion of the target surface was measured by laser Doppler vibrometry. The bending modes are identified as the dominant components in the velocity signal recorded as a function of time. With accumulated radiation damage, the bending mode frequency shifts towards higher values. Based on this shift, Young’s modulus of irradiated isotropic graphite is determined by comparison with ANSYS simulations. The increase of Young’s modulus up to 3 times the pristine value for the highest accumulated fluence of 3 × 1013 ions/cm2 is attributed to the beam-induced microstructural evolution into a disordered structure similar to glassy carbon. Young’s modulus values deduced from microindentation measurements are similar, confirming the validity of the method. Beam-induced stress waves remain in the elastic regime, and no large-scale damage can be observed in graphite. Hexagonal boron nitride shows lower radiation resistance. Circular cracks are generated already at low fluences, risking material failure when applied in high-dose environment.

1. Introduction

With a new generation of particle accelerators like the Facility for Antiproton and Ion Research (FAIR) or the High Luminosity LHC at CERN, unprecedented high beam intensities will be reached. Under conditions of short pulse duration, the impact of intense ion beams on intercepting devices such as production targets [1], collimators [2], or beam dumps might result in failure of these components. Experiments measuring pulsed-beam-induced effects were previously conducted at the HiRadMat beamline at CERN [37], at Brookhaven National Laboratory [8], and at SIS18 at GSI [9]. This work presents results from irradiations with ∼1 GeV uranium ions at the linear accelerator UNILAC at GSI Helmholtz Center for Heavy Ion Research in Darmstadt (Germany). The ion mass is higher, and the energy is lower than beams used in earlier studies [39]. The energy close to the Bragg peak provides maximum energy loss and thus leads to fast dose accumulation [1012]. The lower energy has also the advantage of rather limited sample activation. Such beam conditions allow efficient monitoring of the response of the target materials under high-intensity short-pulse beams and as a function of radiation damage accumulation. Experiments with accelerated dose accumulation are important in order to test the suitability of production target or beam dump materials under long-term operation conditions.

In thin targets, beam-induced rapid heat deposition leads to bending and radial and transversal waves, which we observed by measuring the surface velocities using laser Doppler vibrometry. A critical parameter for the numerical modelling of these effects is Young’s modulus. With increasing radiation damage, Young’s modulus can change [13, 14]. This leads to a shift of the frequency of beam-induced bending modes. Comparing experimental frequency data to results of ANSYS simulations allows in situ monitoring of the evolution of Young’s modulus during irradiation. Independent offline microindentation measurements on a fluence series of irradiated samples were used to validate Young’s modulus values deduced during the irradiation. The investigated materials are isotropic graphite and sintered boron nitride (h-BN). Due to its usage in nuclear reactors, radiation damage induced by neutrons and protons in graphite is well documented [1518] while data of swift heavy ion induced damage are scarcer [13, 14]. The boron nitride grade has a hexagonal structure similar to graphite, but its elastic modulus is larger and its thermal conductivity is lower. In a previous proton irradiation experiment, X-ray diffraction analysis showed a smaller shift and less broadening of the (002) reflection indicating a higher radiation hardness of h-BN compared to graphite [19]. In the past, both materials have been considered as candidates for beam dumps and production targets [1, 20] or as luminescence screens in accelerators. Investigating the response of the two materials contributes to a better understanding of how radiation damage affects their response to beam impact and which online monitoring techniques can be used to detect critical degradation of beam intercepting devices in an early stage and thus avoid catastrophic failures during beamtimes.

2. Materials and Methods

The isotropic graphite grade investigated was SIGRAFINE R6650 [21] produced by SGL Group with an average grain size of 7 μm. The disk-shaped sample (diameter of 20 mm) was cut from a larger piece and lapped on both sides to a thickness of 0.5 mm. The h-BN sample was isotropic HeBoSint D100 grade produced by Henze Boron Nitride Products [22] of the same size as the graphite specimen. Table 1 summarizes the properties of the materials [2124]. Both materials have a similar mass density, Poisson ratio, coefficient of thermal expansion, and heat capacity, but h-BN has a ∼2 times higher Young’s modulus and a ∼4 times lower thermal conductivity than graphite. For the same impacting beam intensity, higher amplitudes of thermal stress waves have to be expected.

The discs of isotropic graphite and h-BN were mounted in a ring-shaped sample holder as shown in Figure 1(a). The outer millimetre of the samples is clamped to the holder, resulting in a free-standing diameter of the disc of 18 mm. The holders were mounted at the M3 beamline of the UNILAC using the setup depicted in Figure 1(b). The targets were exposed to 100 µs-long U-ion pulses with a kinetic energy of 1.14 GeV and a beam repetition rate of 0.5 Hz. The beam was shaped with horizontal and vertical slits, resulting in a homogenous beam spot with a square shape and an area of 10 × 10 mm2 positioned in the centre of the samples. Position, shape, and homogeneity of the beam were controlled with a luminescence screen mounted on the sample holder as shown in Figure 1(a). The energy loss at the surface and the beam penetration depth were calculated with the SRIM-2013 code [25] (Figure 2). In Table 2, they are listed together with the beam parameters. Beam-induced heating of the samples was controlled with a thermal camera and kept below 200°C. This secures that the properties of the target materials remain close to the values at room temperature given in the data sheets. The thermal camera was positioned normal to the beam axis observing at the same time the temperature on the front and back side of the sample with help of two infrared gold mirrors. The surface velocity related to the beam-induced bending and stress waves in the samples was measured with a Polytec OFV-525/-5000-S laser Doppler vibrometer (LDV). The laser is directed to the sample and reflected from the surface. The frequency of the reflected light of a moving object is changed due to the Doppler effect. Comparing the frequency shift of the reflected light to the initial laser signal allowed us to determine the velocity of the sample surface and thus derive the frequency and amplitude of thermoelastic stress waves. The LDV was aligned at the rear side at 45° angle to avoid any beam-induced damage to the LDV and the corresponding view port at the beamline. A tiny piece of adhesive reflective tape with an area of about 1 mm2 was glued at the centre of the sample to increase the reflectivity for the laser beam and to avoid the loss of the reflected laser beam during surface displacement. The LDV was connected to an oscilloscope allowing signal acquisition with a temporal resolution of 1 ns. Before each LDV measurement, the beam was switched off to let the samples cool down to room temperature. Subsequently, the velocity signal generated by the first pulse, after switching on the beam again, was recorded.

To benchmark the results of the LDV measurements, simulations with ANSYS 2019 R2 release were conducted. A quarter of the disc was simulated using the properties given in Table 1. A schematic drawing of geometry and orientation is shown in Figure 3. The energy loss of the ions was simplified by depositing the total kinetic energy across the penetration depth resulting in a constant heat generation of 3 MW/cm³ in the irradiated volume for the pulse duration of 100 μs. To simulate the effects of beam-induced material degradation, simulations were also conducted for isotropic graphite with increased Young’s moduli and reduced thermal conductivities [26] in the irradiated volume.

For indentation experiments, SIGRAFINE R6650 square samples of 10 × 10 mm2 with a thickness of 2 mm were cut and lapped from the top. The samples were irradiated to different fluences with a flux of 1–3 × 109 ions/cm2·s. One half of each sample was masked with a thick aluminium foil which allowed a direct comparison of the irradiated part with the pristine part. Microindentation using a NanoTest Vantage system by Micro Materials with a Berkovich indenter tip was conducted on the irradiated samples up to a maximum indentation depth of 10 μm (contact area larger than 1000 μm2). This is large enough to penetrate several grains and at the same time small enough to remain in the region of constant electronic energy loss (Figure 2) and to avoid the influence of the pristine substrate. The unloading curves were analysed by the Oliver–Pharr method to obtain Young’s modulus [27, 28]. Measurements on at least 25 positions were taken on each sample both for the pristine part and the part exposed to the ion beam.

3. Results and Discussion

3.1. Isotropic Graphite Grade SIGRAFINE R6650

To estimate the maximum temperature rise per pulse, it is assumed that the entire power density of 3 MW/cm³ of a U-ion pulse is transferred into heat. This value was inserted in the ANSYS simulations continuously over the duration of the 100 μs-long ion pulse. We simulated a quarter of a disc with a diameter of 20 mm and a thickness of 0.5 mm with adiabatic boundary conditions (justified by the short time scales of the involved processes). Figure 4(a) shows the simulated temperature evolution of a single beam pulse on an isotropic graphite disc at different depths and radial positions. The temperature of the quasipristine sample increases during the duration of the pulse and reaches a maximum of 156°C at the end of the ion pulse at 100 μs in the centre of the beam spot (y = 0, z = 0). The generated heat is conducted through the sample equilibrating within 1 ms to almost similar temperature values on both surfaces of the disc (y = 0, compare z = 0 with z = 0.5). Over the entire time, the outer rim of the disc (y = 10, z = 0) remains at room temperature, confirming that the use of adiabatic conditions in the investigated time frame is adequate.

To simulate the case of a radiation-damaged sample, the thermal conductivity of the volume exposed to the beam was reduced to 4 W/m·K [26]. The density was kept constant since radiation-induced swelling was shown to be below 5% [29, 30]. The values of specific heat and coefficient of thermal expansion are expected to change towards those of glassy carbon [26, 31]. As the difference of both these properties is small between graphite and glassy carbon, they were kept constant in the simulation. Due to the degraded properties, the temperature at the sample surface (y = 0, z = 0) increases to 240°C after a U-ion pulse of 100 μs duration. A temperature of 255°C is obtained if no heat conduction is included, illustrating the severely limited heat conduction from the front of the sample within the time duration of the pulse. The temperature in a depth where the ions stop (y = 0, z = 0.06) is lower than observed for the pristine sample. For damaged material, the dissipation of the energy deposited in the central beam spot is obviously slowed down, whereas heat conduction into the pristine surrounding is not affected. After 1 ms, the temperature difference between the front and the back surface of the sample, at the central spot of the beam (y = 0, compare z = 0 with z = 0.5), is still larger than 50°C.

The generated heat in the beam spot leads to thermal expansion of the material. Since the heating is initially affecting only one surface of the disc, bending modes are generated. To understand the signal of the LDV, the expected first-order bending frequency fb1 of a fixed and a free circular plate is calculated using equation (1) [32]. The respective radius of the unclamped part (0.9 cm) or the free disc (1 cm) is r, the thickness is t, Young’s modulus is E, the density is ρ, and Poisson’s ratio is ν. λ is a constant depending on the binding condition and Poisson’s ratio and corresponds to 3.196 for a fixed binding condition and 2.958 (graphite) and 2.970 (h-BN) for a free binding condition [32].

Due to the short time of the ion pulse, the thermal inertia of the outer material prevents the expansion of the heated beam spot, causing thermal stress. This stress generates stress waves which travel through the sample and are reflected at the interfaces. The frequency of the transversal elastic wave ftr, which is reflected between the back and front surfaces of the sample, is calculated using the following equation [33]:

The calculated frequencies for the investigated materials are summarized in Table 3. The frequencies expected for the transversal stress waves are three orders of magnitude larger, allowing an easy differentiation from the bending frequencies.

Figures 5(a) and 5(b) show velocity signals and corresponding fast Fourier transform (FFT) spectra measured with the LDV at the rear side (x = y ≈ 0, z = 0.5) of a 500 μm-thick SGL R6650 disc impacted with 1.14 GeV U-ions. For clarity, only data for the quasipristine sample as well as after the accumulation of a fluence of 1 × 1013 and 3.5 × 1013 ions/cm2 are presented. Due to the alignment of the LDV under an angle of 45°, the measured amplitude is significantly reduced and proper transformation from the voltage signal to velocity is difficult. Moreover, the amplitude slightly fluctuates from pulse to pulse due to beam intensity variations. Thus, the raw voltage signal is presented and focus is put on the frequency spectra of the involved processes. In the pristine sample, a harmonic oscillation with a frequency of 5.38 ± 0.02 kHz is generated in excellent agreement with the frequency fb1, free calculated for the free binding condition (Table 3). With increasing U-ion fluence, the frequency signal increases to 6.10 ± 0.02 kHz at the maximum applied fluence of 3.5 × 1013 ions/cm2. We ascribe the increase of the frequency to an increase of Young’s modulus in the irradiated area, assuming that the boundary conditions remained constant. Neither swelling was detected by means of profilometry nor macroscopic failures of the samples such as cracks were observed using optical microscopy. Irradiation of graphite results in the formation of interstitials and vacancies. Agglomeration of these defects creates vacancy clusters [34]. These defects prevent movement of dislocations and sliding of basal planes, causing a significant hardening and increase of Young’s modulus.

ANSYS simulations of the deformation of the graphite disc were performed to infer beam-induced changes of Young’s modulus. In a first step, the thermal input of Figure 4(a) was used on a pristine disc. The velocity and the corresponding FFT spectrum in the centre of the rear side (x = y = 0, z = 0.5) of the disc with free boundary conditions are shown in Figures 5(c) and 5(d). An oscillation with a frequency of 5.32 ± 0.02 kHz is obtained, similar to the calculated value using equation (1) and measured with the LDV on the pristine sample.

In a second step, ANSYS simulations of discs with larger Young’s modulus in the beam spot were conducted. In Figures 5(c) and 5(d), examples of the velocity in the centre of the rear side of discs with 2 and 3 times the pristine Young’s modulus are shown. The properties of the part outside the beam spot were assumed to stay constant since the temperature increase is not high enough to significantly change the properties (as seen in Figure 4). The frequency spectra of the LDV recordings at different fluences can now be compared with the frequency results of the ANSYS simulations as demonstrated in Figure 6(a) for 2 and 3 times the pristine Young’s modulus. The frequencies for those two simulations were measured in the LDV signal at accumulated fluences of 1.5 × 1013 and 3.3 × 1013 ions/cm2, respectively. The same procedure was used to estimate a continuous development of Young’s modulus depending on the fluence as shown in Figure 6(b).

To test the validity of this method, the deduced Young’s modulus values are compared to microindentation measurements on a sample series with different accumulated fluences. From the good agreement in the investigated fluence range, it can be concluded that the LDV online measurement method provides realistic results for beam-induced changes of Young’s modulus for isotropic graphite. This is excellent news because it reduces the need to irradiate samples series with different accumulated fluences and avoids uncertainties due to property fluctuations between several samples.

Both methods show a large increase up to about 3 times the pristine Young’s modulus E0 for irradiation with a fluence of 3 × 1013 U-ions/cm2. The obtained increase is similar to the ratio of Young’s moduli of glassy carbon [31] and SIGRAFINE R6650 [21] of 2.8. This supports the assumption that intensive ion irradiation induces a transition towards a disordered structure similar to glassy carbon [14, 26]. Compared to neutron irradiation, the effects of heavy ion irradiation are much larger. The fluences achieved in this study correspond to less than 0.01 dpa as calculated by SRIM Kinchin and Pease damage calculation [25, 35]. Several tens of dpa are necessary for similar changes induced by neutron irradiation [1517].

Besides the large amplitude caused by the bending of the sample, also elastic or plastic transversal stress waves are expected at a pulse length of 100 μs and a deposited power of approximately 3 MW/cm³ [36]. These transversal waves propagate along the disk axis and, according to the values calculated in Table 3, appear in the MHz frequency regime. Time-resolved frequency analysis of the LDV signal was performed by continuous wavelet transformation with a Morlet wavelet with a wave number of 20. The results for the first dynamic response obtained for a pristine isotropic graphite target are shown in Figure 7. The highest wavelet coefficient appears at a frequency of 2.70 ± 0.02 MHz, very close to the calculated frequency of 2.73 MHz. This signal corresponds to an elastic stress wave as the pressure wave propagates at the speed of sound [37]. The yield strength of graphite is not reached and no plastic deformation occurs. This supports our hypothesis that the sample properties outside the beam spot are not changed as inserted in the ANSYS simulations. The maximum of the wavelet coefficient is obtained at the end of the pulse of 100 μs. At this point, the highest values of temperature and velocity are reached, resulting in the highest thermal stress. After the end of the beam pulse, the temperature difference decreases and the stress wave signal is rapidly damped.

3.2. Hexagonal Boron Nitride

Isotropic hexagonal BN has a similar structure as graphite. However, its polar bonding results in a lower resistance to crack propagation. LDV velocity signals and corresponding FFT spectra measured at the rear side (x = y ≈ 0, z = 0.5) of the sample are shown in Figure 8. At a low fluence of 4 × 1011 ions/cm2, a vibration with a main frequency of 9.53 ± 0.07 kHz was measured. According to the calculation with equation (1), 10.07 kHz corresponds to a bending mode for a fixed sample. ANSYS simulations with a fixed support at the outer 1 mm rim show good agreement of the waveform with the measured velocity signals as presented in Figure 9.

Figure 10 shows the development of the bending frequency as a function of fluence in the frequency range around 10 kHz (Figure 10(a)) and 40 kHz (Figure 10(b)). The low-frequency signal remains rather stable at about 9.5 kHz up to 2 × 1012 ions/cm2 and then increases to 10.6 kHz for 1 × 1013 ions/cm2. For higher fluences, the signal fluctuates irregularly.

In the high-frequency regime, a new double-peak signal appears at a fluence of 2 × 1012 ions/cm2 with centres at 38.4 ± 0.1 and 39.4 ± 0.1 kHz (Figure 8(b)). These frequencies remain stable up to a fluence of 1 × 1013 ions/cm2. With further fluence increase (Figure 10(b)), a fluctuation of high-frequency peaks in the FFT spectra similar to those in the low frequency regime occurs. We ascribe these frequency effects to the formation of cracks (>2 × 1012 ions/cm2) and finally to the complete detachment of the beam-exposed central part of the sample (>1 × 1013 ions/cm2).

Our assumption is supported by the observation during dismounting of the irradiated sample, where we identified a round piece of 8-9 mm in diameter completely detached from the rest of the sample (Figure 11). The fracture piece shows pronounced bending which is probably due to beam-induced swelling. Based on this one sample, quantitative evaluation of this effect was unfortunately not possible.

By using equation (1) for a sample diameter of 8.5 mm, which is the average diameter of the detached central part of the h-BN target, a free bending frequency of 39.0 kHz is calculated, fitting well the measured frequencies. The appearance of two peaks is ascribed to the noncircular shape of the detached piece. ANSYS simulations indicate that a Young's modulus of approximately 4 times the pristine value is required in the beam spot to reproduce the experimentally observed frequency changes at a fluence of 1 × 1013 ions/cm2. For a detached sample piece, this should result in even larger changes of the frequency within the beam spot which is not the case (Figure 10(b)). This indicates that the increase of the low frequency is primarily caused by cracks in the sample and not by an increase of Young’s modulus.

In contrast to isotropic graphite, visible degradation of irradiated h-BN appears at smaller fluences, indicating a lower radiation resistance. While the graphitic samples kept their shape and integrity, h-BN shows a visible deformation which eventually results in cracking and fragmentation of the structure. First failure signs are visible in the velocity and FFT signals at fluences as low as 2 × 1012 ions/cm2. The damage effect becomes more pronounced at higher fluences as shown for the signal at a fluence of 1 × 1013 ions/cm2. While h-BN was demonstrated to be stable under continuous proton irradiation [19], the larger energy loss and the pulsed structure of the U-ion beam result in stress concentration profiles appearing at the interface between irradiated and pristine material leading to crack formation. Due to the brittle nature of h-BN, the irradiated part of the target finally detaches from the sample even though the penetration depth of the ion beam is only about one tenth of the sample thickness. In contrast, graphite is characterized by slow crack growth and self-healing of preexisting microcracks due to growth in the c-axis direction [38, 39], leading to higher critical stress values. Moreover, in irradiated graphite, stress is accommodated by creep deformation preventing catastrophic failures.

4. Conclusions

Discs of isotropic graphite and hexagonal boron nitride were exposed to U-ion pulses with GeV kinetic energy. The dynamic response of the targets was monitored by recording the surface velocity signal at the rear side using laser Doppler vibrometry. Beam-induced changes of the first-order bending frequency of the isotropic graphite disc allow online determination of Young’s modulus of irradiated material by comparison with ANSYS simulations. Young’s modulus rises to approximately 3 times the pristine value for a fluence of 3.5 × 1013 U-ions/cm2. The comparison with microindentation data shows good agreement, suggesting the online method as an interesting alternative to offline measurements of a series of samples that have been irradiated at different fluences. The observation by continuous wavelet transformation enables the detection of transversal elastic stress waves for the investigated pulse intensities of up to 1.6×1010 ions/cm2 per pulse. In the case of h-BN, cracks formed around the beam spot at a fluence of 2 × 1012 ions/cm2 which eventually result in failure of the disc. Our observations provide clear evidence that graphite exhibits a higher radiation hardness making it superior for applications in targets and beam dumps for high power accelerators. To provide safe operation conditions and estimate the life time of beam intercepting devices, the large increase of Young’s modulus in graphite needs to be taken into account when simulating beam conditions and designing components. Our study also demonstrates that laser Doppler vibrometry is a suitable tool to monitor failures of beam-exposed devices by early detection of defect formation in materials, resulting in shifting frequencies or appearance of new frequency peaks, before catastrophic failures happen.

Data Availability

The experimental and simulation data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

The results presented here are based on a UMAT experiment, which was performed at the M-branch of the UNILAC at the GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt (Germany), in the frame of FAIR Phase-0.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was funded by the European Union’s Horizon 2020 Research and Innovation Program under grant agreement no. 730871. PB gratefully acknowledges the support from the HGS-HIRe Graduate School.