Latent Tracks in Ion-Irradiated LiTaO3 Crystals: Damage Morphology Characterization and Thermal Spike Analysis

Systematic research on the response of crystal materials to the deposition of irradiation energy to electrons and atomic nuclei has attracted considerable attention since it is fundamental to understanding the behavior of various materials in natural and manmade radiation environments. This work examines and compares track formation in LiTaO3 induced by separate and combined effects of electronic excitation and nuclear collision. Under 0.71–6.17 MeV/u ion irradiation with electronic energy loss ranging from 6.0 to 13.8 keV/nm, the track damage morphologies evolve from discontinuous to continuous cylindrical zone. Based on the irradiation energy deposited via electronic energy loss, the subsequently induced energy exchange and temperature evolution processes in electron and lattice subsystems are calculated through the inelastic thermal spike model, demonstrating the formation of track damage and relevant thresholds of lattice energy and temperature. Combined with a disorder accumulation model, the damage accumulation in LiTaO3 produced by nuclear energy loss is also experimentally determined. The damage characterizations and inelastic thermal spike calculations further demonstrate that compared to damage-free LiTaO3, nuclear-collision-damaged LiTaO3 presents a more intense thermal spike response to electronic energy loss owing to the decrease in thermal conductivity and increase in electron–phonon coupling, which further enhance track damage.


Introduction
Recent research has demonstrated that perovskite-type LiTaO 3 crystals, as representative functional oxides, have extensive potential applications in various radiation environments, such as tritium breeding in fusion reactors [1,2], nuclear waste immobilization [3], and radiation detection [4]. In the research field of ion-solid interactions and irradiation effects, systematic study of the response of crystal materials to the deposition of irradiation energy to electrons and atomic nuclei has always attracted considerable attention. As a research focus, understanding and clarifying the relevant physical mechanisms is not only fundamental to predicting and evaluating the damage behaviors of materials in natural and man-made radiation environments [5], but it is also key to the applications of ion and electron beam techniques in atomic-level defect manipulation, material modification, and micro/nanofabrication [6][7][8][9][10][11][12][13][14][15][16].
In this work, through appropriate choices of irradiating ion energies, velocities, and energy loss intensities, the role of irradiation energy deposited to electron and lattice subsystems and the induced lattice-defect production and structure evolution in LiTaO 3 crystals are studied comprehensively. By combining the microstructure characterizations with inelastic thermal spike model calculations, the formation of latent tracks with different damage morphologies under the action of electronic energy loss and the corresponding threshold conditions are analyzed. Based on a disorder accumulation model [17,18], the damage accumulation in LiTaO 3 produced by the nuclear energy loss process is also determined. Then, the coupled effect between the electronic excitation and nuclear collision on the thermal spike response and the induced enhanced track damage behavior are further experimentally characterized and theoretically demonstrated.

Energetic Ion Irradiation and Damage Characterization
Optically polished single-crystal LiTaO 3 samples with a (006) surface normal zone axis direction were used in this work. Irradiation experiments with 1 MeV (0.005 MeV/u) 197 Au + and 20 MeV (0.71 MeV/u) 28 Si 3+ at 300 K were performed utilizing a 1.7 MV tandem acceleratior (NEC, Middleton, USA) and a 6 MV tandem accelerator (HVEE, Amersfoort, The Netherlands), respectively, within the State Key Laboratory of Nuclear Physics and Technology, Peking University. Irradiation experiments with 247 MeV (6.17 MeV/u) 40 Ar 12+ and 358 MeV (6.17 MeV/u) 58 Ni 19+ at 300 K were performed at the Heavy Ion Research Facility in Lanzhou (HIRFL), Institute of Modern Physics, Chinese Academy of Sciences. Relatively low ion-beam current density and ion flux (8.3 × 10 10 cm −2 ·s −1 for Au + , 1.9 × 10 10 cm −2 ·s −1 for Si 3+ , 6.5 × 10 10 cm −2 ·s −1 for Ar 12+ , and 8.1 × 10 10 cm −2 ·s −1 for Ni 19+ ) were maintained during the irradiation process to avoid undesired heating and charge accumulation effects on the LiTaO 3 samples. Rutherford backscattering (RBS)/channeling analysis was carried out with the abovementioned 1.7 MV tandem accelerator. A 2.0 MeV He + beam was extracted, and a Si detector located at a scattering angle of 160 • relative to the He + beam was used to collect the backscattered He + signal and assess the irradiation damage. Cross-sectional LiTaO 3 samples for transmission electron microscopy (TEM) observations were prepared through polishing, dimpling, and Ar-ion milling with a Gatan 695 precision polishing system (Gatan Inc., Pleasanton, USA), and the corresponding irradiation damage regions were characterized utilizing an FEI Tecnai G2 F20 transmission electron microscope (FEI Co., Hillsboro, OR, USA).

Simulations and Calculations of Irradiation Damage Models
As indicated in Table 1, the electronic/nuclear energy losses (E ele /E nucl ) and displacements per atom (dpa) in LiTaO 3 induced by energetic Au + , Si 3+ , Ar 12+ , and Ni 19+ irradiation were determined through the stopping and range of ions in matter (SRIM) 2013 simulation code with full-damage cascade mode [19,20], in which the density and default displacement energies of Li, Ta, and O atoms were set to 7.45 g·cm −3 , 25 eV, 25 eV, and 28 eV, respectively [21,22]. Based on the irradiation energy deposited by electronic energy loss, the induced energy exchange and diffusion and temperature evolution processes in the electron and lattice subsystems of the LiTaO 3 crystal were numerically calculated through two heat-conduction differential equations described in the inelastic thermal spike (iTS) model [23][24][25]. In the iTS model, the energy exchange between the electronic and lattice subsystems is described by a set of thermal diffusion equations, one for the electronic subsystem (Equation (1)) and one for the atomic subsystem (Equation (2)) Crystals 2020, 10, 877 where the thermal conductivity K a and specific heat coefficient C a of the lattice subsystem versus lattice temperature were extracted from [26]; the thermal conductivity K e and specific heat coefficient C e of the electron subsystem were fixed to 100 W m −1 ·K −1 and 1.0 J cm −3 ·K −1 , respectively, owing to the constant thermodynamic parameters in the electron subsystem [23,24]; the electron-phonon mean free path λ was deduced from the band gap (4.7 eV) and set to 4.2 nm [27]; the electron-phonon coupling parameter g was set to 5.7 × 10 18 W·m −3 ·K −1 [26]; the electron-phonon mean free time τ was 1.75 × 10 −13 s owing to τ = C e /g [28]; and the latent heat of fusion was set to 150 J/g. The term A(r, t) describes the spatial and temporal energy deposition from the incident ion to the electrons within 10 −17 -10 −15 s, and the electrons would deposit their energy close to the ion path within 10 −15 s [24,29]. The damage accumulation in LiTaO 3 produced by nuclear energy loss originating from low-energy ion irradiation was also fitted through a disorder accumulation model.

Latent Ion Tracks Induced by Electronic Energy Loss
As mentioned above, SRIM 2013 simulation code is primarily used to determine the profile of electronic energy loss as a function of ion penetration path, and could not provide the radial distribution of deposited irradiation energy via electronic energy loss process at a certain penetration depth, which actually is essential for the iTS calculations so as to obtain the radial evolutions of electron and lattice temperatures and further describe the lattice melting behavior. During the ion slowing process, the interaction cross section between the irradiating ion and target electrons and the subsequently induced electron cascade in the target directly depend on the ion velocity, which would further determine the radial distribution of the energy deposited in the electron subsystem via the electronic energy loss process. Considering a cylindrical volume in which 0.66 of the total electronic energy is deposited, the corresponding cylindrical radius (radial distance from the ion trajectory) is defined as the absorption radius α e to describe the energy distribution profile [24]. The value of α e would increase with increasing ion velocity (velocity effect), resulting in a large energy spread region and locally low energy deposition density [24]. In this work, as shown in Figure 1a, the α e of the LiTaO 3 crystal used for electronic energy deposition increases from 2.1 to 4.7 nm as the ion velocity increases from 0.71 to 6.17 MeV/u. Under ion irradiation with velocities of 0.71, 1.15, 2.15, and 6.17 MeV/u, the fraction of irradiation energy deposited on the electron subsystem derived from the Katz model and MC calculations [24,29,30] and the corresponding normalized energy deposition distribution F(r) versus the radial distance from the ion trajectory are shown in Figure 1b. Then, considering the SRIM-simulated electronic energy loss (E ele ) relative to the ion velocity and ensuring that the integral of the energy deposition density over the radial area equals E ele [24], the determined distribution profile of the energy deposition density in the electron subsystem is shown in Figure 1c, illustrating that during the ion irradiation process, for a similar E ele (6.0-6.2 keV/nm), the deposition density of electronic energy would accordingly increase with decreasing ion velocity from 6.17 to 0.71 MeV/u (orange and wine curves in Figure 1c). For the same ion velocity (6.17 MeV/u), the deposition density of electronic energy would accordingly increase with increasing E ele from 6.0 to 12.0 keV/nm (orange and blue curves in Figure 1c). For the ion (1.15 MeV/u) with an E ele of 8.4 keV/nm and the ion (6.17 MeV/u) with an E ele of 12.0 keV/nm, the deposition density of electronic energy in the cylindrical core area are similar (pink and blue curves in Figure 1c). For the ion (2.15 MeV/u) with an E ele of 13.8 keV/nm, the deposition density of electronic energy reaches the maximum (green curve in Figure 1c).
Crystals 2020, 10, x FOR PEER REVIEW 4 of 12 (orange and blue curves in Figure 1c). For the ion (1.15 MeV/u) with an Eele of 8.4 keV/nm and the ion (6.17 MeV/u) with an Eele of 12.0 keV/nm, the deposition density of electronic energy in the cylindrical core area are similar (pink and blue curves in Figure 1c). For the ion (2.15 MeV/u) with an Eele of 13.8 keV/nm, the deposition density of electronic energy reaches the maximum (green curve in Figure 1c). Under the action of electronic energy deposition, the induced energy exchange and diffusion and temperature evolution in electron and lattice subsystems are discussed utilizing the inelastic thermal spike (iTS) model [23][24][25]. Figure 1d illustrates the iTS-calculated energy evolution in a LiTaO3 crystal under 0.71 MeV/u ion irradiation with an Eele of 6.2 keV/nm. As shown, the kinetic energy of the incident ion is gradually deposited into the electronic subsystem via electronic energy loss, reaching a peak (6.2 keV/nm) at approximately 10 −15 s. The radial distribution of the energy deposition density in the electron subsystem at this time is described by the wine curve in Figure 1c. Then, as shown in the inset of Figure 1d, owing to the much higher temperature in the electron subsystem than in the lattice subsystem, the deposited electronic energy is gradually transferred from the hot electronic subsystem to the cooler lattice through electron-phonon coupling, leading to significant evolution of the electronic and lattice temperatures; once the energy transferred and deposited to lattice atoms exceeds a certain threshold, local melting occurs, and lattice damage along the ion trajectory is formed via a subsequent rapid quenching process. In this energy exchange stage, Under the action of electronic energy deposition, the induced energy exchange and diffusion and temperature evolution in electron and lattice subsystems are discussed utilizing the inelastic thermal spike (iTS) model [23][24][25]. Figure 1d illustrates the iTS-calculated energy evolution in a LiTaO 3 crystal under 0.71 MeV/u ion irradiation with an E ele of 6.2 keV/nm. As shown, the kinetic energy of the incident ion is gradually deposited into the electronic subsystem via electronic energy loss, reaching a peak (6.2 keV/nm) at approximately 10 −15 s. The radial distribution of the energy deposition density in the electron subsystem at this time is described by the wine curve in Figure 1c. Then, as shown in the inset of Figure 1d, owing to the much higher temperature in the electron subsystem than in the lattice subsystem, the deposited electronic energy is gradually transferred from the hot electronic subsystem to the cooler lattice through electron-phonon coupling, leading to significant evolution of the electronic and lattice temperatures; once the energy transferred and deposited to lattice atoms exceeds a certain threshold, local melting occurs, and lattice damage along the ion trajectory is formed via a subsequent rapid quenching process. In this energy exchange stage, the sum of irradiation energy deposited into the electron and lattice subsystems remains constant (6.2 keV/nm). Finally, due to the decrease in electron temperature and increase in lattice temperature, the energy exchange driven by the temperature difference ceases to occur since the temperatures between the electron and lattice subsystems reach equilibrium within approximately 10 −10 s.
As shown in Figure 2, the temporal evolution profiles of the energy deposition to atoms (Figure 2a,d,g,j,m) and the temporal evolution of the atomic temperature (Figure 2b,e,h,k,n) are shown for the different irradiation conditions at different radii (distance from the ion path), obtained by the iTS model [23][24][25]. The related cross-sectional TEM images in LiTaO 3 under ion irradiation with different ion velocities and electronic energy losses are also presented (Figure 2c,f,i,l,o). For 0.71 MeV/u Si 3+ irradiation with an E ele of 6.2 keV/nm, the energy deposited to the lattice atoms within a 1 nm radius exceeds the criterion E m for melting phase formation, and the corresponding lattice temperature exceeds the melting temperature T m . It should be pointed out that, according to the distribution of specific heat C a as a function of the temperature T a in LiTaO 3 crystal [26], the E m for melting phase formation is calculated to be 0.54 eV/atom (1113 J/g) via E m = T m T 0 C a (T a )dT a + E l , where T 0 is set to 300 K, and the latent heat of fusion E l is set to 150 J/g [31]. This results in the formation of discontinuous tracks with a radius of~1.0 nm, as shown in Figure 2c. For 6.17 MeV/u Ar 12+ irradiation with an E ele of 6.0 keV/nm, the energy deposited to the lattice atoms along the ion path reaches a peak at 0.47 eV/atom within~10 −13 s, and the corresponding atomic temperature increases to peak at 1860 K, which is still lower than the melting temperature T m (1923 K). Therefore, as shown in Figure 2f, no melting phase or lattice damage would occur [32]. For 1.15 MeV/u Ar 12+ irradiation (E ele of 8.4 keV/nm) and 6.17 MeV/u Ni 19+ irradiation (E ele of 12.0 keV/nm), accompanied by increasing energy deposition on the lattice atoms, the radii of the produced cylindrical melting zones increase to~1.5 nm and~2.0 nm, respectively. As shown in Figure 2i,l, this results in the formation of discontinuous tracks of larger size, with~1.5 nm and~2.0 nm diameter, respectively. In the case of 2.15 MeV/u Ni 19+ irradiation with an E ele of 13.8 keV/nm, the energy deposited to the lattice atoms in the cylinder center reaches 1.20 eV/atom, significantly exceeding the melting criterion, and lattice melting occurs in a cylindrical region with a radius of~3.0 nm. Owing to the much higher distributions of atomic energy and temperature, continuous track damage with a radius of~2.5 nm shown in Figure 2o is formed during the subsequent rapid quenching process [26].

Damage Accumulation Induced by Nuclear Energy Loss
Utilizing 1 MeV Au + irradiation, the damage production and accumulation of LiTaO3 crystals due to nuclear energy loss only, originating from the elastic collisions between irradiating ions and target atomic nuclei, were determined. The SRIM-simulated nuclear energy loss (Enucl) profile and displacements per atom (dpa) profile corresponding to the Au + fluence of 1 × 10 14 cm −2 are shown in Figure 3a, illustrating that the dpa peak at 0.71 is located at a depth of ~100 nm. The measured RBS/channeling spectra of Au + -irradiated LiTaO3 under different fluences are shown in Figure 3b,

Damage Accumulation Induced by Nuclear Energy Loss
Utilizing 1 MeV Au + irradiation, the damage production and accumulation of LiTaO 3 crystals due to nuclear energy loss only, originating from the elastic collisions between irradiating ions and target atomic nuclei, were determined. The SRIM-simulated nuclear energy loss (E nucl ) profile and displacements per atom (dpa) profile corresponding to the Au + fluence of 1 × 10 14 cm −2 are shown in Figure 3a, illustrating that the dpa peak at 0.71 is located at a depth of~100 nm. The measured RBS/channeling spectra of Au + -irradiated LiTaO 3 under different fluences are shown in Figure 3b, and then, based on an iterative procedure [33], the profiles of Ta sublattice disorder as a function of depth are calculated and shown in Figure 3c, indicating the long-range disorder in LiTaO 3 under Au + irradiation to the local dpa of 0.25. The Ta sublattice disorder in the damage peak region versus the local irradiation dose (dpa) is further plotted in Figure 3d, and the obtained damage accumulation curve is described by a disorder accumulation model [17,18], where the irradiation-induced disorder S (black curve) measured by the RBS/channeling technique would arise from the sum of two contributions: amorphous fraction f a (blue curve) and interstitials and small interstitial clusters in residual crystalline regions S d (red curve). The related fitting parameters-including the amorphization cross section σ a , the effective cross section σ s for defect-stimulated amorphization, the saturation value S * d for the defect-induced disorder, and the coefficient B proportional to an effective defect recombination volume-are also indicated.
Crystals 2020, 10, x FOR PEER REVIEW 7 of 12 and then, based on an iterative procedure [33], the profiles of Ta sublattice disorder as a function of depth are calculated and shown in Figure 3c, indicating the long-range disorder in LiTaO3 under Au + irradiation to the local dpa of 0.25. The Ta sublattice disorder in the damage peak region versus the local irradiation dose (dpa) is further plotted in Figure 3d, and the obtained damage accumulation curve is described by a disorder accumulation model [17,18], where the irradiation-induced disorder S (black curve) measured by the RBS/channeling technique would arise from the sum of two contributions: amorphous fraction (blue curve) and interstitials and small interstitial clusters in residual crystalline regions (red curve). The related fitting parameters-including the amorphization cross section , the effective cross section for defect-stimulated amorphization, the saturation value * for the defect-induced disorder, and the coefficient proportional to an effective defect recombination volume-are also indicated.

Enhanced Track Damage Induced by Coupled Nuclear and Electronic Energy Losses
To study the coupled effects of electronic and nuclear energy losses, the Au + -irradiated LiTaO3 sample with a fluence of 2 × 10 13 cm −2 was subsequently irradiated with 20 MeV (0.71 MeV/u) Si 3+ to a fluence of 1 × 10 13 cm −2 . The measured RBS/channeling spectra, corresponding Ta disorder profiles, cross-sectional TEM images, and FFT patterns of LiTaO3 crystals under different irradiation conditions are shown in Figure 4, and they are used to discuss the Eele-induced track damage response in pristine and Enucl-damaged LiTaO3. According to the RBS/channeling spectra shown in Figure 4a, the determined disorder profiles of the Ta sublattice in ion-irradiated LiTaO3 samples shown in Figure 4b,e,h quantitatively indicated that in contrast to pristine LiTaO3, the same electronic energy loss induced by Si 3+ irradiation in Enucl-damaged LiTaO3 leads to complete disorder in the region at a depth of ~100 nm, demonstrating the additive effect of electronic and nuclear energy losses on

Enhanced Track Damage Induced by Coupled Nuclear and Electronic Energy Losses
To study the coupled effects of electronic and nuclear energy losses, the Au + -irradiated LiTaO 3 sample with a fluence of 2 × 10 13 cm −2 was subsequently irradiated with 20 MeV (0.71 MeV/u) Si 3+ to a fluence of 1 × 10 13 cm −2 . The measured RBS/channeling spectra, corresponding Ta disorder profiles, cross-sectional TEM images, and FFT patterns of LiTaO 3 crystals under different irradiation conditions are shown in Figure 4, and they are used to discuss the E ele -induced track damage response in pristine and E nucl -damaged LiTaO 3 . According to the RBS/channeling spectra shown in Figure 4a, the determined disorder profiles of the Ta sublattice in ion-irradiated LiTaO 3 samples shown in Figure 4b,e,h quantitatively indicated that in contrast to pristine LiTaO 3 , the same electronic energy loss induced by Si 3+ irradiation in E nucl -damaged LiTaO 3 leads to complete disorder in the region at a Crystals 2020, 10, 877 8 of 12 depth of~100 nm, demonstrating the additive effect of electronic and nuclear energy losses on damage production. Compared to pristine LiTaO 3 , the enhanced track damage response in the predamaged LiTaO 3 upon electronic energy loss is also confirmed by TEM observations. As shown in Figure 4b,c, the lattice defects and defect clusters were first produced by nuclear energy loss via Au + irradiation; then, upon subsequent 0.71 MeV/u Si 3+ irradiation, the introduced E ele approximatively locates at the Bragg peak position, and remains constant (6.2 keV/nm) in the near surface region; finally, the induced track damage in the corresponding E nucl -damaged region is shown in Figure 4h,i. This damage is much more significant than the E ele -induced damage in pristine LiTaO 3 (Figure 4e,f).
Crystals 2020, 10, x FOR PEER REVIEW 8 of 12 damage production. Compared to pristine LiTaO3, the enhanced track damage response in the predamaged LiTaO3 upon electronic energy loss is also confirmed by TEM observations. As shown in Figure 4b,c, the lattice defects and defect clusters were first produced by nuclear energy loss via Au + irradiation; then, upon subsequent 0.71 MeV/u Si 3+ irradiation, the introduced Eele approximatively locates at the Bragg peak position, and remains constant (6.2 keV/nm) in the near surface region; finally, the induced track damage in the corresponding Enucl-damaged region is shown in Figure 4h,i. This damage is much more significant than the Eele-induced damage in pristine LiTaO3 (Figure 4e,f). To understand the damage response of Au + -irradiated (Enucl-damaged) LiTaO3 upon subsequent Si 3+ irradiation with intense electronic energy loss, the corresponding energy deposited to atoms and the consequent atomic temperature evolutions were further calculated using the iTS model [23][24][25]. In contrast to pristine LiTaO3, the lattice defects and defect clusters existing in the Au + -irradiated sample scatter electrons and phonons, decreasing thermal conductivities ( and ) of the electron and lattice subsystems and the electron-phonon mean free path , thus leading to an increasing electron-phonon coupling parameter . Based on the previous research work [5,23,[34][35][36], during the iTS calculation process for Au + -irradiated LiTaO3 crystals, and are assumed to reduce to 75% and is assumed to be 40% larger compared to the original values of pristine LiTaO3. To understand the damage response of Au + -irradiated (E nucl -damaged) LiTaO 3 upon subsequent Si 3+ irradiation with intense electronic energy loss, the corresponding energy deposited to atoms and the consequent atomic temperature evolutions were further calculated using the iTS model [23][24][25]. In contrast to pristine LiTaO 3 , the lattice defects and defect clusters existing in the Au + -irradiated sample scatter electrons and phonons, decreasing thermal conductivities (K e and K a ) of the electron and lattice subsystems and the electron-phonon mean free path λ, thus leading to an increasing electron-phonon coupling parameter g. Based on the previous research work [5,23,[34][35][36], during the iTS calculation process for Au + -irradiated LiTaO 3 crystals, K e and K a are assumed to reduce to 75% and g is assumed to be 40% larger compared to the original values of pristine LiTaO 3 . Comparing Figure 5a,c to Figure 5b,d, under 0.71 MeV/u Si 3+ irradiation with an E ele of 6.2 keV/nm, the energy deposited to the lattice along the ion path increases from 0.68 to 1.07 eV/atom within~10 −13 s, and the corresponding atomic temperature increases from 2400 to 3860 K, demonstrating that the introduction of E nucl -induced lattice damage would effectively enhance the response intensity of the E ele -induced thermal spike. In addition, the melting temperature of the crystal decreases owing to the lattice defects and defect clusters embedded in the crystalline system, further promoting the formation of the melting phase and track damage via the subsequent rapid quenching process. Therefore, the experimentally observed additive impacts of the coupled effects of nuclear and electronic energy losses upon track damage production is theoretically demonstrated based on the evolution of energy deposition and atomic temperature in the lattice subsystem.
Crystals 2020, 10, x FOR PEER REVIEW 9 of 12 electronic energy losses upon track damage production is theoretically demonstrated based on the evolution of energy deposition and atomic temperature in the lattice subsystem.

Conclusions
In the research field of ion-solid interactions and irradiation effects, a key challenge is to clarify the deposition of irradiation energy to electron and lattice subsystems and further understand the induced damage production of various materials in irradiation environments. In this work, under ion irradiation with 0.71-6.17 MeV/u with electronic energy loss ranging from 6.0 to 13.8 keV/nm, the experimentally observed latent tracks with different damage morphologies were theoretically discussed based on inelastic thermal spike analysis. The effect of electronic energy loss and absorption radius relative to ion velocity on irradiation energy deposition and exchange processes between electron and lattice subsystems and the further induced evolutions of electronic and atomic temperatures are determined, describing the formation of the lattice melting phase and latent track damage in LiTaO3 crystals. In addition, the lattice defects and defect clusters produced via nuclear energy loss (i) decrease the thermal conductivity and increase the electron-phonon coupling, leading to the enhanced thermal spike response under electronic energy loss and (ii) accordingly reduce the melting temperature of the crystalline system, similarly to [34,[36][37][38], demonstrating the coupled mechanism by which electronic and nuclear energy losses promote track damage.

Conclusions
In the research field of ion-solid interactions and irradiation effects, a key challenge is to clarify the deposition of irradiation energy to electron and lattice subsystems and further understand the induced damage production of various materials in irradiation environments. In this work, under ion irradiation with 0.71-6.17 MeV/u with electronic energy loss ranging from 6.0 to 13.8 keV/nm, the experimentally observed latent tracks with different damage morphologies were theoretically discussed based on inelastic thermal spike analysis. The effect of electronic energy loss and absorption radius relative to ion velocity on irradiation energy deposition and exchange processes between electron and lattice subsystems and the further induced evolutions of electronic and atomic temperatures are determined, describing the formation of the lattice melting phase and latent track damage in LiTaO 3 crystals. In addition, the lattice defects and defect clusters produced via nuclear energy loss (i) decrease the thermal conductivity and increase the electron-phonon coupling, leading to the enhanced thermal spike response under electronic energy loss and (ii) accordingly reduce the melting temperature of the crystalline system, similarly to [34,[36][37][38], demonstrating the coupled mechanism by which electronic and nuclear energy losses promote track damage.