Athermal repair of nanoscale defects in optical materials using a femtosecond laser

Ion implantation is widely used to fabricate advanced optical and optoelectronic materials and devices. However, nanoscale defects generated during the ion implantation process severely affect the quality and properties of the material and device. Here, combining computational simulations and experiments, we investigate the mechanism for defect repair in fused silica after Cu ion implantation using femtosecond laser irradiation with an energy fluence much lower than the ablation threshold. Atomic force microscopy demonstrates no unexpected ablation. The optical absorption spectra show that various types of defects with formation energies between 1.9 and 6.2 eV can be repaired successfully via an athermal procedure. The Raman spectra imply that the broken chemical bonds reconnect after femtosecond laser irradiation. Our study reveals that low-energy femtosecond laser irradiation can transfer the appropriate energy needed to repair defects; thus it could be useful in fabricating nonlinear optical devices due to its high spatial selectivity and convenience.

With its high optical nonlinearity and ultrafast response time, fused silica, which is doped with metallic particles by ion implantation, is an attractive material for all-optical communication and optoelectronic devices. 1,2 During the ion implantation process, the introduced ions generate trains of defects in lattices. [3][4][5][6][7][8][9][10][11][12][13][14][15][16] As illustrated in SRIM (the Stopping and Range of Ions in Matter) model simulations, [3][4][5] the major defects are vacancies and interstitial atoms. [6][7][8][9][10] These nanoscale defects severely influence the linear and nonlinear optical properties as well as crystalline construction. 3,4 In order to reduce radiation damage, a repair process is followed by or synchronized with ion implantation. Conventionally, thermal annealing is used which heats specimens up to 1300 K. The drawback of this method is that it takes a considerably long time, resulting in a remarkable loss of implanted ions. 11 In addition, when the remaining defects producing special properties are desir-able, such as absorption humps or an enhancement in electrical resistivity, this method cannot be used because of the lack of spatial selectivity. 12,13 Efforts have been made to improve the repair process. Nanosecond laser irradiation after ion implantation has been introduced to decrease the repair time. However, similar to thermal annealing, the method suffers from loss of ions. 14 Low-energy nanosecond laser irradiation is also reported to repair lattice defects in ion implantation. 15 The non-equilibrium state generated by the interaction between the ions and the material facilitates the energy absorption from the lowenergy nanosecond laser, which provides enough energy to shift interstitial atoms to vacancies, avoiding unexpected ablation. However, the setup needs the nanosecond laser system and the ion implantation system to have precise cooperation, which is too complex to be widely used.
For precise and efficient defect repair in the nanoscale, we use a low-energy femtosecond laser to repair lattice defects in Cu ion-implanted fused silica. The process and mechanism of nanodefect repair are investigated. We studied three sets of characterization methods that had shown to prove the effectiveness and selectiveness of this method.

Results and discussion
Characterization using an optical microscope Pure fused silica is green under an optical microscope as shown in Fig. 1a color changed to dark brown (Fig. 1b). When a part of the asimplanted area was irradiated with a femtosecond laser with a fluence of 0.2 J cm −2 , the color of the irradiated area changed remarkably (Fig. 1c), similar to that with pure silica. Furthermore, the micro/nanoscale morphology of the area of interest was scanned by atomic force microscopy (AFM). For the as-implanted specimen, there are a great many nanoscale clusters measured about 200 nm wide and 3 nm high. The nanoscale point defects and smaller clusters are along with nanoscale clusters. However, when this as-implanted specimen has been irradiated with the focused low energy femtosecond laser, the number of clusters is substantially reduced. Besides, there are no craters in the irradiated area, implying no unexpected ablation.
In order to obtain an appropriate laser energy used for irradiation, the laser fluence was varied. When the fluence was above 0.5 J cm −2 , an unwanted nanoscale crater in the irradiated area had been observed, which is shown in Fig. S2 in the ESI. † When it was set below 0.1 J cm −2 , there was no obvious change of color, completely different from Fig. 1c. The optimized fluence is approximately 0.2 J cm −2 , which is much lower than the ablation threshold (2.7 J cm −2 ) of fused silica. 17 Furthermore, to exclude the incubation effects of multiple irradiations which could lead to ablation, the number of laser pulses which stayed at the same point of the sample was reduced to less than 10, and the time interval between two adjacent pulses was set to be longer than 20 milliseconds. So far, for a specimen with an implanted area of 5 × 5 mm 2 , it takes about half an hour to have most of the defects repaired by this method.
We hypothesize that the Cu atoms and the large number of defects approximately 10 nm below the surface might be responsible for the change in color. In addition, the femtosecond laser irradiation with a laser fluence of 0.2 J cm −2 can affect the distribution of the Cu atoms or the defects.

Optical absorption spectra
To understand the mechanism of the femtosecond laser irradiation method, the optical absorption spectra of both the as-implanted and the irradiated samples have been characterized by using an ultraviolet-visual spectrophotometer, as shown in Fig. 2. For the as-implanted sample, an obvious hump from 200 nm to 650 nm was detected, which gives evidence of defects, 18 including the B-band (∼413 nm), nonbridging oxygen holes (∼258 nm), the B2-band (250 nm and 241 nm: oxygen-deficient centers), and E′ centers (214 nm). 31,32 Additionally, a peak near 580 nm was also found, which refers to the surface plasmon resonance (SPR) of Cu nanoparticles. After irradiation with the femtosecond laser, the humps decreased by 25% which indicates a reduction in the defects, and the peak disappeared because the Cu nanoparticles could be splintered and re-shaped by the femtosecond laser irradiation. 19

Raman spectra
To obtain insights into the defect repair mechanism, we studied the Raman spectra of three samples: pure, asimplanted, and irradiated silica (Fig. 3). The spectrum of pure silica has two characteristic peaks which are D1 (492 cm −1 ) and D2 (605 cm −1 ). After implantation, a new peak appears at 1579 cm −1 , which corresponds to the vibration frequency of O 2 . 20 It means that the Si-O bonds were broken and O-center defects appeared during ion implantation. Additionally, D1 was markedly reduced, and D2 almost disappeared, due to the heavy damage of the lattice construction of fused silica. After femtosecond laser irradiation, however, the peak at 1579 cm −1 was significantly decreased, indicating that the O-center defects were rehabilitated. Additionally, D2 reappeared in the irradiated area, which implies the reconstruction of Si-O bonds. Such a reconstruction instead of destruction of Si-O bonds might be attributed to the fact that the repairing fluence of the femtosecond laser is well below the ablation threshold of silica. 34 We might conclude that the low fluence that is smaller than the ablation threshold fluence of the femtosecond laser could be the key for defect repair.
Copper is well known to be extremely mobile within materials, diffusing along interfaces and defects in materials under a host of different environmental factors. As such, the effective control of copper suspensions is critical for defect repair. We found that the copper could exist in two types in our cases: Cu nanoparticles and CuO compounds. Once the Cu ions were implemented with a fluence of 3 × 10 17 cm −2 , there were saturated Cu ions which aggregate and form nanoparticles. This could be verified by the SPR peak shown in Fig. 2. The Raman spectra (Fig. 3) also show the existence of Cu nanoparticles. When irradiated by using the femtoscecond laser, Cu atoms react with SiO 2 to form CuO. In the meantime, some Cu oxide compounds could form Cu atoms which then aggregate into nanoparticles. The concentrations of Cu nanoparticles and compounds depend on the energy density and time duration of the femtosecond laser radiation. It is worth mentioning that we did not observe the loss of Cu in our study as there is no ablation during our low energy fluence femtosecond laser radiation.
To gain more insight into the defect repair, we carried out a simulation of the elemental copper distribution via Stopping and Range of Ions in Matter (SRIM) modeling. The details are in the following simulation subsection. It can be seen that a high concentration of Cu ions is implanted in the target's surface with a depth of approximately 70 to 300 nm ( Fig. 4). At the same time, a train of defects, including vacancies and interstitial atoms, were formed in the material, which can be seen in Fig. 5. The higher the voltage is, the deeper the ions reach and the more defects exist.

Simulations using the SRIM model
To figure out the distribution of the defects induced by the collision between copper ions and fused silica, a SRIM model was built to simulate the process of ion implantation with the acceleration voltages of 30, 60, 100 and 190 kV, in which 60 kV was the most commonly used. The results are shown in Fig. 4. For an acceleration voltage of 30 and 60 kV, Cu ions are implanted in the target's surface with a depth of approximately 300 nm, and the distribution is a Gaussian shape. For 100 and 190 kV, the distributions become irregular and much deeper. Generally, the densities at the surface and the deep end are much smaller, and the density at the middle is obviously Fig. 2 Optical absorption spectra. The optical absorption spectra of the as-implanted specimen (dark square) and the femtosecond laser irradiated area (red circle). The absorption spectra qualitatively show the density of various kinds of defects, such as B-band (413 nm), non-bridging oxygen holes (258 nm), B2-band (250 nm and 241 nm) and E' centers (214 nm). The reduction of the optical absorption in the irradiated area implies repairing of these defects by the femtosecond laser irradiation.   larger. No matter how big the voltage is, the ions are localized at the very shallow surface. It can be seen that a high concentration of Cu ions is implanted in the target's surface with a depth approximately 70 to 300 nm. This localization of Cu ions is a good property for micro/nano-fabrication. Our simulation also shows that at the same time of ion implantation, a train of nanoscale defects, including vacancies and interstitial atoms, are formed in the material, which can be seen in Fig. 5 for the ion implantation with the acceleration voltages of 30, 60, 100 and 190 kV, respectively. Point defects are introduced during collisions between the implanted ions and atoms in the matrix. Avalanche collisions are induced by the recoil atoms. When the kinetic energy is large enough to kick the atom away from the lattice, there will be a vacancy marked in the figure. The higher the voltage is, the deeper the ions reach and the more defects exist. As the voltage grows larger, the defects expand dramatically in a direction parallel to the surface. The depth of the defects is 1.5 times that of the depth of ions due to the recoil atoms. In addition, the higher the voltage is, the deeper the ions reach and the more defects exist. Despite these variations, the processes stay at almost the same level and the evolution after femtosecond laser irradiation is considered to be similar.

Repair process
The evolution of the molecular structures during the repair process is derived from the three sets of aforementioned characterization. Firstly, for the as-implanted specimen, SiO 2 molecules were splintered by the introduced Cu particles and Si-O bonds were broken. O-Center defects appeared, as shown schematically in Fig. 6a. Then it was irradiated with a lowenergy femtosecond laser. In this stage, energy in electromagnetic radiation is transferred to the electrons promoting them from their ground state into some excited state. Due to their much smaller mass, electrons move much faster than ions. Electron-electron collisions are prominent in energy transfer, resulting in the electronic energies adopting a Fermi-Dirac distribution defining an electronic temperature. The maximum electronic temperature will be achieved in approximately 0.1 picoseconds. The energy is then transferred to the ions, via electron-phonon coupling, leading to an increase in the ionic temperature or lattice temperature. The lattice will melt if the ionic temperature is higher than the melting temperature. This process is called thermal melting and it takes on the order of a few picoseconds. For femtosecond laser irradiation, a large fraction (10% or more) of the valence electrons will be excited to conduction electrons 33 and free electrons (Fig. 6c) which strongly modify the interatomic forces, resulting in coherent atomic motions and the irradiated sample will melt on time scales shorter than the thermal melting. The displaced atoms obtain the laser energy in the energy transition (Fig. 6c). Some of them possess enough kinetic energy to move back to some vacant positions. As a result, the chemical bonds are recovered and the defects are repaired (Fig. 6d).
Additionally, the copper suspensions are roughly controlled by the low-energy femtosecond laser irradiation. The copper atoms irradiated by the femtosecond laser begin to accumulate. When the density of the localized copper atoms reaches a certain degree, the Cu nanoparticles were separated out in fused silica, which would play an important role in enhancing the useful nonlinear properties of the material. Furthermore, no copper particles are found on the surface of fused silica from the AFM images, which means that the ejection of large amounts of copper atoms did not happen.
Compared with traditional thermal annealing methods, femtosecond laser irradiation achieves high spatial selectivity because the focused laser can be applied selectively to small areas that require reparation, as shown in Fig. 1c. When remaining defects producing special properties are desirable, such as producing some absorption humps or an enhancement in electrical resistivity, the selectivity of repair could be very useful. In comparison with a nanosecond laser, femtosecond laser irradiation has much less loss of metallic particles. Both experiments [21][22][23][24] and calculations 25,26 have shown that nanosecond laser irradiation can lead to rapid melting of the material and hence annihilate the defects. The advantage of the femtosecond laser over nanosecond laser irradiation is that there is no unexpected ablation, there are no remaining metallic particles in the specimen, and the rapid repair of defects is achieved. The reason for the effective defect repair in our study could be the optimization of two competing factors, laser fluence and pulse duration. On one hand, the laser is strong enough to cause the excitation of a large fraction of valence electrons. On the other hand, the femtosecond laser has a much shorter pulse duration, which is even shorter than the electron-lattice relaxation time to avoid ablation. Therefore, the laser energy absorption is completed during the non-thermal process in significant electron-lattice non-equilibrium states. [27][28][29][30] It is worth noting that this approach of athermal reparation of nanoscale defects via a femtosecond laser is only applicable to optical materials so that the laser can penetrate into the bulk and effectively have the nanoscale defect repaired. In addition, though high spatial selectivity and low efficiency seem to be two sides of the same coin in this study, our approach could be combined with other thermal and ionization processes 35 to further improve the efficiency of defect reparation.

Conclusion
In conclusion, a low energy femtosecond laser has been used to repair the nanoscale defects in fused silica after Cu ion implantation. The optimized laser fluence is approximately 0.2 J cm −2 , which is much lower than the ablation threshold. We investigated the mechanism of defect repair in fused silica after Cu ion implantation using femtosecond laser irradiation with an energy fluence much lower than the ablation threshold. The optical absorption spectra show that the defects can be repaired. The Raman spectra imply that the broken chemical bonds reconnect after the femtosecond laser irradiation. All three sets of characterization showed that interstitial atoms caused by the introduced ions could gain suitable energy to fit vacancies, and then the broken chemical bond recovered. The mechanism of the defect repair in our femtosecond laser irradiation treatment is athermal due to the strongly modified inter-atomic forces caused by the excitation of electrons.

Experimental setup and methods
The fused silica specimens of 10 mm × 10 mm and 1 mm thickness had both surfaces polished, and then were implanted by using a MEVVA ion implanter. The ion beam was Cu with 60 kV acceleration voltage. The dose rate was 2 μA cm −2 and the total dose used for the experiment was 3 × 10 17 ions per cm 2 .
The laser used in the experiment was a commercial chirped-pulse regenerative oscillator-amplifier femtosecond laser system (MaiTai, Spectra-Physics). The laser pulse was linearly polarized, with a central wavelength of 800 nm, a pulse width of 50 fs (FWHM), and a repetition rate of 1 kHz. The fundamental energy is attenuated to the desired level by using a neutral density filter. The pulse was focused on the specimens' surface at normal incidence using a 5× microscope objective lens (NA = 0.13, Olympus), which generated a spot with a diameter of approximately 9 μm. The fused silica substrate was mounted onto a translation stage. To create a suitable irradiation area, the fused silica substrate was translated relative to the laser beam, exposing any given spot on the silicon surface to an average of 100 pulses. The optical properties of the as-implanted silica before and after laser irradiation, as well as pure silica, were characterized by using an optical microscope (BX51, Olympus). The absorption spectra were obtained by using an ultraviolet-visual spectrophotometer (Lambda 35, PE) and the information of lattice defects was characterized by using a Raman spectrometer (inVia-Reflex, Renishaw). The topography of the surface was observed by using an atomic force microscope (Dimension Edge, Bruker).

SRIM simulations details
SRIM, a package of computer programs which calculate the interactions of ions with matter, is based on a Monte Carlo simulation method, the binary collision approximation with a random selection of the impact parameter of the next colliding ion. The SRIM model assumes that when an atom/ion collides with a target atom in the lattice and the energy E trans transferred from the incident atom/ion is larger than the displacement energy E disp , which is the energy used to kick the target atom away from the lattice with a certain distance so that it cannot restore immediately, the target atom will be kicked away from the lattice. Because of the minimum energy loss E latt used for the atom displacement from the lattice, the recoil energy of the displaced atom is derived from E recoil = E trans − E latt . If E recoil is larger than E disp , it will collide with other target atoms and more defects of vacancies and interstitial atoms could exist.
The defect types and their dynamic distributions are important in the process of defect repairs. To shed more light on this, we studied the distribution of the oxygen and silicon atoms through the SRIM simulations for the ion implantation with acceleration voltages of 30, 60, 100 and 190 kV, respectively. The results are shown in Fig. 7. We observed that the number of recoil oxygen atoms is much larger than silicon for the simple reason that oxygen is easier to be knocked away from the lattice by copper ions than silicon due to the lighter mass. The recoil oxygen atoms contribute large numbers of defects. The distribution of recoil atoms is similar to that of the copper ions in Fig. 4, which means that the number of recoil atoms is roughly decided by the number of copper atoms.