Fine structure in swift heavy ion tracks in amorphous SiO2

The Effect of Electrolyte Properties on Ionic Transport through Solid-State Nanopores: Experiment and Simulation

Nanopore membranes enable versatile technologies that are employed in many different applications, ranging from clean energy generation to filtration and sensing. Improving the performance can be achieved by conducting numerical simulations of the system, for example, by studying how the nanopore geometry or surface properties change the ionic transport behavior or fluid dynamics of the system. A widely employed tool for numerical simulations is finite element analysis (FEA) using software, such as COMSOL Multiphysics. We found that the prevalent method of implementing the electrolyte in the FEA can diverge significantly from physically accurate values. It is often assumed that salt molecules fully dissociate, and the effect of the temperature is neglected. Furthermore, values for the diffusion coefficients of the ions, as well as permittivity, density, and viscosity of the fluid, are assumed to be their bulk values at infinite dilution. By performing conductometry experiments with an amorphous SiO2 nanopore membrane with conical pores and simulating the pore system with FEA, it is shown that the common assumptions do not hold for different mono- and divalent chlorides (LiCl, NaCl, KCl, MgCl2, and CaCl2) at concentrations above 100 mM. Instead, a procedure is presented where all parameters are implemented based on the type of salt and concentration. This modification to the common approach improves the accuracy of the numerical simulations and thus provides a more comprehensive insight into ion transport in nanopores that is otherwise lacking.

Direct Fabrication of Sub-10 nm Nanopores in WO3 Nanosheets Using Single Swift Heavy Ions

Extending the fabrication methodology of solid-state nanopores in a wide range of materials is significant in the fields of single molecule detection, nanofluidic devices, and nanofiltration membranes. Here, we demonstrate a new method to directly fabricate size- and density-controllable sub-10 nm nanopores in WO₃ nanosheets using single swift heavy ions (SHIs) without any chemical etching process. By selecting ions of different electronic energy losses (S_e), nanopores with sizes from 1.8 to 7.4 nm can be created in WO₃ nanosheets. The creation efficiency of nanopores achieves ∼100% for S_e > 20 keV/nm, and there exists a critical thickness below which nanopores can be created. Combined with molecular dynamics simulations, we propose that the viscosity and surface tension of the transient molten phase caused by SHIs are the key factors for the formation of nanopores. This method paves a way to fabricate solid-state nanopores in materials with a low viscosity and surface tension.

Insights into nanoparticle shape transformation by energetic ions

Shape modification of embedded nanoparticles can be achieved by means of swift heavy ion irradiation. During irradiation, the particles elongate and align with the direction of the ion beam, presumably due to nanometer-scale phase transitions induced by individual ion impacts. However, the details of this transformation are not fully understood. The shape of metal nanoparticles embedded in dielectric matrices defines the non-linear optical properties of the composite material. Therefore, understanding the transformation process better is beneficial for producing materials with the desired optical properties. We study the elongation mechanism of gold nanoparticles using atomistic simulations. Here we focus on long-timescale processes and adhesion between the nanoparticle and the matrix. Without the necessity of ad-hoc assumptions used earlier, our simulations show that, due to adhesion with the oxide, the nanoparticles can grow in aspect ratio while in the molten state even after silicon dioxide solidifies. Moreover, they demonstrate the active role of the matrix: Only explicit simulations of ion impacts around the embedded nanoparticle provide the mechanism for continuous elongation up to experimental values of aspect ratio. Experimental transmission electron microscopy micrographs of nanoparticles after high-fluence irradiation support the simulations. The elongated nanoparticles in experiments and their interface structures with silica, as characterized by the micrographs, are consistent with the simulations. These findings bring ion beam technology forward as a precise tool for shaping embedded nanostructures for various optical applications.

Examining Different Regimes of Ionization-Induced Damage in GaN Through Atomistic Simulations

The widespread adoption of gGaN in radiation-hard semiconductor devices relies on a comprehensive understanding of its response to strongly ionizing radiation. Despite being widely acclaimed for its high radiation resistance, the exact effects induced by ionization are still hard to predict due to the complex phase-transition diagrams and defect creation-annihilation dynamics associated with group-III nitrides. Here, the Two-Temperature Model, Molecular Dynamics simulations and Transmission Electron Microscopy, are employed to study the interaction of Swift Heavy Ions with GaN at the atomic level. The simulations reveal a high propensity of GaN to recrystallize the region melted by the impinging ion leading to high thresholds for permanent track formation. Although the effect exists in all studied electronic energy loss regimes, its efficiency is reduced with increasing electronic energy loss, in particular when there is dissociation of the material and subsequent formation of N₂ bubbles. The recrystallization is also hampered near the surface where voids and pits are prominent. The exceptional agreement between the simulated and experimental results establishes the applicability of the model to examine the entire electronic energy loss spectrum. Furthermore, the model supports an empirical relation between the interaction cross sections (namely for melting and amorphization) and the electronic energy loss.

Structure Formation and Regulation of Au Nanoparticles in LiTaO3 by Ion Beam and Thermal Annealing Techniques

The size uniformity and spatial dispersion of nanoparticles (NPs) formed by ion implantation must be further improved due to the characteristics of the ion implantation method. Therefore, specific swift heavy ion irradiation and thermal annealing are combined in this work to regulate the size and spatial distributions of embedded Au NPs formed within LiTaO₃ crystals. Experimental results show that small NPs migrate to deeper depths induced by 656 MeV Xe³⁵⁺ ion irradiation. During thermal annealing, the growth of large Au NPs is limited due to the reductions in the number of small Au NPs, and the migrated Au NPs aggregate at deeper depths, resulting in a more uniform size distribution and an increased spatial distribution of Au NPs. The present work presents a novel method to modify the size and spatial distributions of embedded NPs.

SAXS data modelling for the characterisation of ion tracks in polymers

Here, we present new models to fit small angle X-ray scattering (SAXS) data for the characterization of ion tracks in polymers. Ion tracks in polyethylene terephthalate (PET), polycarbonate (PC), polyimide (PI) and polymethyl methacrylate (PMMA) were created by swift heavy ion irradiation using ¹⁹⁷Au and ²³⁸U with energies between 185 MeV and 2.0 GeV. Transmission SAXS measurements were performed at the Australian Synchrotron. SAXS data were analysed using two new models that describe the tracks by a cylindrical structure composed of a highly damaged core with a gradual transition to the undamaged material. First, we investigate the 'Soft Cylinder Model', which assumes a smooth function to describe the transition region by a gradual change in density from a core to a matrix. As a simplified and computational less expensive version of the 'Soft Cylinder Model', the 'Core Transition Model' was developed to enable fast fitting. This model assumes a linear increase in density from the core to the matrix. Both models yield superior fits to the experimental SAXS data compared with the often-used simple 'Hard Cylinder Model' assuming a constant density with an abrupt transition.

Incident Angle Dependent Formation of Ion Tracks in Quartz Crystal with C60+ Ions: Big Ions in Small Channels

Quartz (SiO2) crystals possess intrinsic columnar pores perpendicular to (0001) surfaces, consisting of three- and six-membered ring (3MR and 6MR) structures of Si and O atoms. The diameters of the larger pores, i.e., 6 MRs, are ~0.49 nm, while the diameters of fullerene (C60) ions are 0.7 nm, i.e., larger than either type of the pores. Transmission electron microscopy observation evidenced approximately two times longer ion tracks in the channeling condition, i.e., 0° incidence to (0001) surface, than an off-channeling condition, i.e., 7° incidence in this case, under 6 MeV C60 ion injection. The track length at the 0° incidence decreases more steeply than that at the 7° incidence with decreasing the energy from 6 MeV to 1 MeV. Finally, the track lengths at the 0° and 7° incidences become comparable, i.e., the channeling-like effect disappears at 1 MeV irradiation. This study experimentally indicates that the channeling-like effect of C60 ions is induced in quartz crystals, while the sizes of the channels are smaller than the C60 ions.

Ion track etching of polycarbonate membranes monitored by in situ small angle X-ray scattering

In situ small angle X-ray scattering (SAXS) measurements of ion track etching in polycarbonate foils are used to directly monitor the selective dissolution of ion tracks with high precision, including the early stages of etching. Detailed information about the track etching kinetics and size, shape, and size distribution of an ensemble of nanopores is obtained. Time resolved measurements as a function of temperature and etchant concentration show that the pore radius increases almost linearly with time for all conditions and the etching process can be described by an Arrhenius law. The radial etching shows a power law dependency on the etchant concentration. An increase in the etch rate with increasing temperature or concentration of the etchant reduces the penetration of the etchant into the polymer but does not affect the pore size distribution. The in situ measurements provide an estimate for the track etch rate, which is found to be approximately three orders of magnitude higher than the radial etch rate. The measurement methodology enables new experiments studying membrane fabrication and performance in liquid environments.

C60 ions of 1 MeV are slow but elongate nanoparticles like swift heavy ions of hundreds MeV

This study reports that high fluence fullerene ion (C₆₀⁺) irradiation of 1–6 MeV, which was made possible by a new-type of high-flux ion source, elongates metal nanoparticles (NPs) in amorphous SiO₂ as efficiently as swift heavy ions (SHIs) of 200 MeV Xe¹⁴⁺, i.e., two orders of the magnitude higher energy ions. Comparing the irradiation effects induced by both the beams, the stopping processes of C₆₀ ions in SiO₂ are discussed in this paper. Despite of having almost the same elongation efficiency, the C₆₀⁺ irradiation induced ~10 times more efficient sputtering due to the clustering enhancement and/or the synergy effect. Ion tracks of ~10.4 nm in diameter and 60–80 nm in length were observed in crystalline SiO₂ under 4 MeV C₆₀ irradiation_. While the track diameter was comparable to those by SHIs of the same electronic stopping, much shorter track lengths than those predicted by a rigid C₆₀ molecule model indicates that the fragmentation occurred due to nuclear collisions. The elongation of the metal NPs was induced only down to the depth where the tracks were observed but not beyond.

Etched ion tracks in amorphous SiO2 characterized by small angle x-ray scattering: influence of ion energy and etching conditions

Small angle x-ray scattering was used to study the morphology of conical structures formed in thin films of amorphous SiO₂. Samples were irradiated with 1.1 GeV Au ions at the GSI UNILAC in Darmstadt, Germany, and with 185, 89 and 54 MeV Au ions at the Heavy Ion Accelerator Facility at ANU in Canberra, Australia. The irradiated material was subsequently etched in HF using two different etchant concentrations over a series of etch times to reveal conically shaped etched channels of various sizes. Synchrotron based SAXS measurements were used to characterize both the radial and axial ion track etch rates with unprecedented precision. The results show that the ion energy has a significant effect on the morphology of the etched channels, and that at short etch times resulting in very small cones, the increased etching rate of the damaged region in the radial direction with respect to the ion trajectory is significant.

The Effect of Heavy Ion Irradiation on the Forward Dissolution Rate of Borosilicate Glasses Studied in Situ and Real Time by Fluid-Cell Raman Spectroscopy

Borosilicate glasses are the favored material for immobilization of high-level nuclear waste (HLW) from the reprocessing of spent fuel used in nuclear power plants. To assess the long-term stability of nuclear waste glasses, it is crucial to understand how self-irradiation affects the structural state of the glass and influences its dissolution behavior. In this study, we focus on the effect of heavy ion irradiation on the forward dissolution rate of a non-radioactive ternary borosilicate glass. To create extended radiation defects, the glass was subjected to heavy ion irradiation using 197Au ions that penetrated ~50 µm deep into the glass. The structural damage was characterized by Raman spectroscopy, revealing a significant depolymerization of the silicate and borate network in the irradiated glass and a reduction of the average boron coordination number. Real time, in situ fluid-cell Raman spectroscopic corrosion experiments were performed with the irradiated glass in a silica-undersaturated, 0.5 M NaHCO3 solution at temperatures between 80 and 85 °C (initial pH = 7.1). The time- and space-resolved in situ Raman data revealed a 3.7 ± 0.5 times increased forward dissolution rate for the irradiated glass compared to the non-irradiated glass, demonstrating a significant impact of irradiation-induced structural damage on the dissolution kinetics.

Oriented crystallization of PEG induced by confinement in cylindrical nanopores: structural and thermal properties

Nanoporous ion track-etched polycarbonate is ideally suited for the study of confined polymers via small angle X-ray scattering (SAXS) due to the strictly parallel orientation of the pores as well as their uncorrelated lateral distribution. Nanopores with radii ranging from 17 to 213 nm are prepared and coated with SiO2via atomic layer deposition in order to obtain a well-defined and homogeneous surface. A low molecular weight polyethylene glycol (PEG) homopolymer with a semicrystalline lamellar bulk structure is introduced into the nanopores via melt infiltration. At high temperatures SAXS measurements confirm a uniform filling of the pores with amorphous polymer. Upon cooling below the melting point of PEG, a concentrical structure of semicrystalline lamellae is revealed for large pore radii. We introduce models which successfully describe the combined scattering from nanopores and semicrystalline or amorphous PEG inside. DSC measurements of the confined polymer show a decrease of melting temperature and heat of fusion per gram polymer upon reduction of the pore radius and hint at a change in the lamellar configuration.

Nanoscale Structuring in Confined Geometries using Atomic Layer Deposition: Conformal Coating and Nanocavity Formation

Nanoscale structuring in confined geometries using atomic layer deposition (ALD) is demonstrated for surfaces of nanochannels in track-etched polymer membranes and in mesoporous silica (SBA-15). Suitable process conditions for conformal ALD coating of polymer membranes and SBA-15 with inorganic oxides (SiO2, TiO2, Al2O3) were developed. On the basis of the oxide-coated layers, nanochannels were further structured by a molecular-templated ALD approach, where calixarene macromolecules are covalently attached to the surface and then embedded into an Al2O3 layer. The removal of calixarene by ozone treatment results in 1–2 nm wide surface nanocavities. Surfaces exposed to different process steps are analyzed by small angle X-ray scattering (SAXS) as well as by X-ray photoelectron and infrared spectroscopy. The proposed nanostructuring process increases the overall surface area, allows controlling the hydrophilicity of the channel surface, and is of interest for studying water and ion transport in confinement.

Nanoscale density variations induced by high energy heavy ions in amorphous silicon nitride and silicon dioxide

The cylindrical nanoscale density variations resulting from the interaction of 185 MeV and 2.2 GeV Au ions with 1.0 μm thick amorphous SiN _x :H and SiO _x :H layers are determined using small angle x-ray scattering measurements. The resulting density profiles resembles an under-dense core surrounded by an over-dense shell with a smooth transition between the two regions, consistent with molecular-dynamics simulations. For amorphous SiN _x :H, the density variations show a radius of 4.2 nm with a relative density change three times larger than the value determined for amorphous SiO _x :H, with a radius of 5.5 nm. Complementary infrared spectroscopy measurements exhibit a damage cross-section comparable to the core dimensions. The morphology of the density variations results from freezing in the local viscous flow arising from the non-uniform temperature profile in the radial direction of the ion path. The concomitant drop in viscosity mediated by the thermal conductivity appears to be the main driving force rather than the presence of a density anomaly.

Radiation-induced disorder in compressed lanthanide zirconates

The effects of swift heavy ion irradiation-induced disordering on the behavior of lanthanide zirconate compounds (Ln₂Zr₂O₇ where Ln = Sm, Er, or Nd) at high pressures are investigated. After irradiation with 2.2 GeV ¹⁹⁷Au ions, the initial ordered pyrochlore structure (Fd3[combining macron]m) transformed to a defect-fluorite structure (Fm3[combining macron]m) in Sm₂Zr₂O₇ and Nd₂Zr₂O₇. For irradiated Er₂Zr₂O₇, which has a defect-fluorite structure, ion irradiation induces local disordering by introducing Frenkel defects despite retention of the initial structure. When subjected to high pressures (>29 GPa) in the absence of irradiation, all of these compounds transform to a cotunnite-like (Pnma) phase, followed by sluggish amorphization with further compression. However, if these compounds are irradiated prior to compression, the high pressure cotunnite-like phase is not formed. Rather, they transform directly from their post-irradiation defect-fluorite structure to an amorphous structure upon compression (>25 GPa). Defects and disordering induced by swift heavy ion irradiation alter the transformation pathways by raising the energetic barriers for the transformation to the high pressure cotunnite-like phase, rendering it inaccessible. As a result, the high pressure stability field of the amorphous phase is expanded to lower pressures when irradiation is coupled with compression. The responses of materials in the lanthanide zirconate system to irradiation and compression, both individually and in tandem, are strongly influenced by the specific lanthanide composition, which governs the defect energetics at extreme conditions.

Hillocks created for amorphizable and non-amorphizable ceramics irradiated with swift heavy ions: TEM study

In a previous study, we found that hillocks (i.e. surface ion tracks) can be imaged using transmission electron microscopy (TEM) by irradiating thin CeO₂ samples with swift heavy ions (SHI) at oblique incidence. In the present study, the same TEM method is applied to Y₃Fe₅O₁₂ (YIG) and three fluorides (CaF₂, SrF₂ and BaF₂) for observing hillocks. For YIG, which is one of the amorphizable materials, hillocks are found to have amorphous features consistent with amorphous features of ion tracks. For the fluorides, it is found that the hillocks do not exhibit amorphous features, and they are composed of nanocrystallites. Although hillocks for YIG and CaF₂ exhibit different crystallographic features, hillock diameter agrees with the molten region diameter predicted by the thermal spike model for both materials. It is found that for YIG the hillock diameter is comparable to the ion track diameter, whereas for the fluorides it is always larger than the ion track diameter. The present result shows the existence of the velocity effect for ion track diameter in CaF₂. It is also found that for fluorides both hillock and ion track diameters vary in the order of cation mass (i.e. CaF₂ < SrF₂ < BaF₂). The above results of hillocks and ion tracks for SHI-irradiated fluorides can be consistently interpreted within the framework of the thermal spike model, if melting and successive recrystallization are assumed.

Permanent modifications in silica produced by ion-induced high electronic excitation: experiments and atomistic simulations

The irradiation of silica with ions of specific energy larger than ~0.1 MeV/u produces very high electronic excitations that induce permanent changes in the physical, chemical and structural properties and give rise to defects (colour centres), responsible for the loss of sample transparency at specific bands. This type of irradiation leads to the generation of nanometer-sized tracks around the ion trajectory. In situ optical reflection measurements during systematic irradiation of silica samples allowed us to monitor the irradiation-induced compaction, whereas ex situ optical absorption measurements provide information on colour centre generation. In order to analyse the results, we have developed and validated an atomistic model able to quantitatively explain the experimental results. Thus, we are able to provide a consistent explanation for the size of the nanotracks, the velocity and thresholding effects for track formation, as well as, the colour centre yield per ion and the colour centre saturation density. In this work we will discuss the different processes involved in the permanent modification of silica: collective atomic motion, bond breaking, pressure-driven atom rearrangement and ultra-fast cooling. Despite the sudden lattice energy rise is the triggering and dominant step, all these processes are important for the final atomic configuration.

Dislocation loop formation by swift heavy ion irradiation of metals

A coupled two-temperature, molecular dynamics methodology is used to simulate the structural evolution of bcc metals (Fe and W) and fcc metals (Cu and Ni) following irradiation by swift heavy ions. Electronic temperature dependent electronic specific heat capacities and electron-phonon coupling strengths are used to capture the full effects of the variation in the electronic density of states. Tungsten is found to be significantly more resistant to damage than iron, due both to the higher melting temperature and the higher thermal conductivity. Very interesting defect structures, quite different from defects formed in cascades, are found to be created by swift heavy ion irradiation in the bcc metals. Isolated vacancies form a halo around elongated interstitial dislocation loops that are oriented along the ion path. Such configurations are formed by rapid recrystallization of the molten cylindrical region that is created by the energetic ion. Vacancies are created at the recrystallization front, resulting in excess atoms at the core which form interstitial dislocation loops on completion of crystallization. These unique defect structures could, potentially, be used to create metal films with superior mechanical properties and interesting nanostructures.

An attempt to apply the inelastic thermal spike model to surface modifications of CaF2 induced by highly charged ions: comparison to swift heavy ions effects and extension to some others material

Surface damage appears on materials irradiated by highly charged ions (HCI). Since a direct link has been found between surface damage created by HCI with the one created by swift heavy ions (SHI), the inelastic thermal spike model (i-TS model) developed to explain track creation resulting from the electron excitation induced by SHI can also be applied to describe the response of materials under HCI which transfers its potential energy to electrons of the target. An experimental description of the appearance of the hillock-like nanoscale protrusions induced by SHI at the surface of CaF₂ is presented in comparison with track formation in bulk which shows that the only parameter on which we can be confident is the electronic energy loss threshold. Track size and electronic energy loss threshold resulting from SHI irradiation of CaF₂ is described by the i-TS model in a 2D geometry. Based on this description the i-TS model is extended to three dimensions to describe the potential threshold of appearance of protrusions by HCI in CaF₂ and to other crystalline materials (LiF, crystalline SiO₂, mica, LiNbO₃, SrTiO₃, ZnO, TiO₂, HOPG). The strength of the electron-phonon coupling and the depth in which the potential energy is deposited near the surface combined with the energy necessary to melt the material defines the classification of the material sensitivity. As done for SHI, the band gap of the material may play an important role in the determination of the depth in which the potential energy is deposited. Moreover larger is the initial potential energy and larger is the depth in which it is deposited.

Formation of swift heavy ion tracks on a rutile TiO2 (001) surface

Nanostructuring of surfaces and two-dimensional materials using swift heavy ions offers some unique possibilities owing to the deposition of a large amount of energy localized within a nanoscale volume surrounding the ion trajectory. To fully exploit this feature, the morphology of nanostructures formed after ion impact has to be known in detail. In the present work the response of a rutile TiO2 (001) surface to grazing-incidence swift heavy ion irradiation is investigated. Surface ion tracks with the well known intermittent inner structure were successfully produced using 23 MeV I ions. Samples irradiated with different ion fluences were investigated using atomic force microscopy and grazing-incidence small-angle X-ray scattering. With these two complementary approaches, a detailed description of the swift heavy ion impact sites, i.e. the ion tracks on the surface, can be obtained even for the case of multiple ion track overlap. In addition to the structural investigation of surface ion tracks, the change in stoichiometry of the rutile TiO2 (001) surface during swift heavy ion irradiation was monitored using in situ time-of-flight elastic recoil detection analysis, and a preferential loss of oxygen was found.

Picosecond metrology of laser-driven proton bursts

Tracking primary radiation-induced processes in matter requires ultrafast sources and high precision timing. While compact laser-driven ion accelerators are seeding the development of novel high instantaneous flux applications, combining the ultrashort ion and laser pulse durations with their inherent synchronicity to trace the real-time evolution of initial damage events has yet to be realized. Here we report on the absolute measurement of proton bursts as short as 3.5±0.7 ps from laser solid target interactions for this purpose. Our results verify that laser-driven ion acceleration can deliver interaction times over a factor of hundred shorter than those of state-of-the-art accelerators optimized for high instantaneous flux. Furthermore, these observations draw ion interaction physics into the field of ultrafast science, opening the opportunity for quantitative comparison with both numerical modelling and the adjacent fields of ultrafast electron and photon interactions in matter.

Experimental investigations of the response of matter to ionization would require extremely fast ion pump pulses. Here, the authors explore a different approach observing ionisation dynamics in SiO₂ glass by generating synchronized proton pulses from the interaction of high-power lasers on a solid target.

Tracing temperature in a nanometer size region in a picosecond time period

Irradiation of materials with either swift heavy ions or slow highly charged ions leads to ultrafast heating on a timescale of several picosecond in a region of several nanometer. This ultrafast local heating result in formation of nanostructures, which provide a number of potential applications in nanotechnologies. These nanostructures are believed to be formed when the local temperature rises beyond the melting or boiling point of the material. Conventional techniques, however, are not applicable to measure temperature in such a localized region in a short time period. Here, we propose a novel method for tracing temperature in a nanometer region in a picosecond time period by utilizing desorption of gold nanoparticles around the ion impact position. The feasibility is examined by comparing with the temperature evolution predicted by a theoretical model.

Optical birefringence of Zn nanoparticles embedded in silica induced by swift heavy-ion irradiation

Zn nanoparticles (NPs) embedded in a silica matrix subjected to irradiation with swift heavy ions of 200 MeV Xe¹⁴⁺ have been found to undergo shape elongation from spheres to prolate-spheroids while maintaining the major axes of the NPs in parallel alignment. The directionally-aligned Zn spheroids enable acquisition of optical properties, such as linear dichroism and birefringence. In this paper, the birefringence of the Zn spheroids was evaluated by the crossed-Nicols (XN) transmittance, where a sample was inserted between a pair of optical polarizers that were set in an orthogonal configuration. Linearly-polarized light aligned by the first polarizer was transformed to an elliptic polarization by the birefringence of the Zn spheroids. The existence of the birefringence was confirmed by the non-zero transmittance of the second polarizer in the orthogonal configuration. The sample irradiated with a fluence of 5.0 × 10¹³ ions/cm² exhibited a maximum XN transmittance of 2.1% at a photon energy of ~4 eV. The XN transmission was observed down to a fluence of 1.0 × 10¹² ions/cm², but reduced below the detection limit at a fluence of 1.0 × 10¹¹ ions/cm². The possible application of the elongated Zn NPs as a polarizer with nanometric thickness working in the near- and mid-ultraviolet region is discussed.

Shape elongation of Zn nanoparticles in silica irradiated with swift heavy ions of different species and energies: scaling law and some insights on the elongation mechanism

Zinc nanoparticles (NPs) embedded in silica were irradiated with swift heavy ions (SHIs) of seven different combinations of species and energies. The shape elongation induced by the irradiations was evaluated by optical linear dichroism (OLD) spectroscopy, which is a sensitive tool for determining the change in the mean aspect ratio (AR) of NPs. Although the mean AR change indicated a linear fluence dependence in the low- and medium-fluence regions, it indicated a nonlinear dependence in the high-fluence region. The data reveal that the elongation efficiency of Zn is correlated with the electronic stopping power 'Se in silica' and is not correlated with either the 'Se in Zn' or the nuclear stopping power. The elongation efficiency plotted as a function of the 'Se in silica' revealed a linear relationship, with a threshold value of ∼2 keV nm(-1), which is the same dependence exhibited by the ion-track formation in silica. The log-log plot showed that the elongation efficiency increased linearly with Se above a critical value of ∼3 keV nm(-1) and steeply decreased with Se to the power of 5 below the critical Se. The steep decrease can be ascribed to the discontinuous nature of the ion tracks, which is expected at Se ∼ 2-4 keV nm(-1) in silica. The fluence Φ dependences of AR - 1 under various irradiations are well-normalized with the electronic energy deposition of SHIs, i.e., the product of Se and Φ, with a Se greater than the same critical value of ∼3 keV nm(-1). The normalized data above the critical value fell on a linear relation, AR(Φ) - 1 ∝ SeΦ, for SeΦ < 2 keV nm(-3) and a sublinear relation, AR(Φ) - 1 ∝ (SeΦ)(1/2) for SeΦ > 2 keV nm(-3). On the basis of these experimental results, we discuss some insights into the elongation mechanism.

Temperature dependence of ion track formation in quartz and apatite

Ion tracks were created in natural quartz and fluorapatite from Durango, Mexico, by irradiation with 2.2 GeV Au ions at elevated temperatures of up to 913 K. The track radii were analysed using small-angle X-ray scattering, revealing an increase in the ion track radius of approximately 0.1 nm per 100 K increase in irradiation temperature. Molecular dynamics simulations and thermal spike calculations are in good agreement with these values and indicate that the increase in track radii at elevated irradiation temperatures is due to a lower energy required to reach melting of the material. The post-irradiation annealing behaviour studied for apatite remained unchanged.

Tracks and voids in amorphous Ge induced by swift heavy-ion irradiation

Ion tracks formed in amorphous Ge by swift heavy-ion irradiation have been identified with experiment and modeling to yield unambiguous evidence of tracks in an amorphous semiconductor. Their underdense core and overdense shell result from quenched-in radially outward material flow. Following a solid-to-liquid phase transformation, the volume contraction necessary to accommodate the high-density molten phase produces voids, potentially the precursors to porosity, along the ion direction. Their bow-tie shape, reproduced by simulation, results from radially inward resolidification.

Multi-scale simulation of structural heterogeneity of swift-heavy ion tracks in complex oxides

Tracks formed by swift-heavy ion irradiation, 2.2 GeV Au, of isometric Gd2Ti2O7 pyrochlore and orthorhombic Gd2TiO5 were modeled using the thermal-spike model combined with a molecular-dynamics simulation. The thermal-spike model was used to calculate the energy dissipation over time and space. Using the time, space, and energy profile generated from the thermal-spike model, the molecular-dynamics simulations were performed to model the atomic-scale evolution of the tracks. The advantage of the combination of these two methods, which uses the output from the continuum model as an input for the atomistic model, is that it provides a means of simulating the coupling of the electronic and atomic subsystems and provides simultaneously atomic-scale detail of the track structure and morphology. The simulated internal structure of the track consists of an amorphous core and a shell of disordered, but still periodic, domains. For Gd2Ti2O7, the shell region has a disordered pyrochlore with a defect fluorite structure and is relatively thick and heterogeneous with different degrees of disordering. For Gd2TiO5, the disordered region is relatively small as compared with Gd2Ti2O7. In the simulation, 'facets', which are surfaces with definite crystallographic orientations, are apparent around the amorphous core and more evident in Gd2TiO5 along [010] than [001], suggesting an orientational dependence of the radiation response. These results show that track formation is controlled by the coupling of several complex processes, involving different degrees of amorphization, disordering, and dynamic annealing. Each of the processes depends on the mass and energy of the energetic ion, the properties of the material, and its crystallographic orientation with respect to the incident ion beam.

SAXS investigations of the morphology of swift heavy ion tracks in α-quartz

The morphology of swift heavy ion tracks in crystalline α-quartz was investigated using small angle x-ray scattering (SAXS), molecular dynamics (MD) simulations and transmission electron microscopy. Tracks were generated by irradiation with heavy ions with energies between 27 MeV and 2.2 GeV. The analysis of the SAXS data indicates a density change of the tracks of ~2 ± 1% compared to the surrounding quartz matrix for all irradiation conditions. The track radii only show a weak dependence on the electronic energy loss at values above 17 keV nm(-1), in contrast to values previously reported from Rutherford backscattering spectrometry measurements and expectations from the inelastic thermal spike model. The MD simulations are in good agreement at low energy losses, yet predict larger radii than SAXS at high ion energies. The observed discrepancies are discussed with respect to the formation of a defective halo around an amorphous track core, the existence of high stresses and/or the possible presence of a boiling phase in quartz predicted by the inelastic thermal spike model.

Single ion induced surface nanostructures: a comparison between slow highly charged and swift heavy ions

This topical review focuses on recent advances in the understanding of the formation of surface nanostructures, an intriguing phenomenon in ion-surface interaction due to the impact of individual ions. In many solid targets, swift heavy ions produce narrow cylindrical tracks accompanied by the formation of a surface nanostructure. More recently, a similar nanometric surface effect has been revealed for the impact of individual, very slow but highly charged ions. While swift ions transfer their large kinetic energy to the target via ionization and electronic excitation processes (electronic stopping), slow highly charged ions produce surface structures due to potential energy deposited at the top surface layers. Despite the differences in primary excitation, the similarity between the nanostructures is striking and strongly points to a common mechanism related to the energy transfer from the electronic to the lattice system of the target. A comparison of surface structures induced by swift heavy ions and slow highly charged ions provides a valuable insight to better understand the formation mechanisms.

Role of thermodynamics in the shape transformation of embedded metal nanoparticles induced by swift heavy-ion irradiation

Swift heavy-ion irradiation of elemental metal nanoparticles (NPs) embedded in amorphous SiO(2) induces a spherical to rodlike shape transformation with the direction of NP elongation aligned to that of the incident ion. Large, once-spherical NPs become progressively more rodlike while small NPs below a critical diameter do not elongate but dissolve in the matrix. We examine this shape transformation for ten metals under a common irradiation condition to achieve mechanistic insight into the transformation process. Subtle differences are apparent including the saturation of the elongated NP width at a minimum sustainable, metal-specific value. Elongated NPs of lesser width are unstable and subject to vaporization. Furthermore, we demonstrate the elongation process is governed by the formation of a molten ion-track in amorphous SiO(2) such that upon saturation the elongated NP width never exceeds the molten ion-track diameter.

Atomistic simulation of track formation by energetic recoils in zircon

We have performed classical molecular dynamics simulations of fission track formation in zircon. We simulated the passage of a swift heavy ion through crystalline zircon using cylindrical thermal spikes with energy deposition (dE/dx) of 2.5-12.8 keV nm( - 1) and a radius of 3 nm. At a low dE/dx of 2.55 keV nm( - 1), the structural damage recovered almost completely and a damage track was not produced. At higher values of dE/dx, tracks were observed and the radius of the track increased with increasing dE/dx. Our structural analysis shows amorphization in the core of the track and phase separation into Si-rich regions near the center of the track and Zr-rich regions near the periphery. These simulations establish a threshold dE/dx for fission track formation in zircon that is relevant to thermochronology and nuclear waste immobilization.

Temperature control in molecular dynamic simulations of non-equilibrium processes

Thermostats are often used in various condensed matter problems, e.g. when a biological molecule undergoes a transformation in a solution, a crystal surface is irradiated with energetic particles, a crack propagates in a solid upon applied stress, two surfaces slide with respect to each other, an excited local phonon dissipates its energy into a crystal bulk, and so on. In all of these problems, as well as in many others, there is an energy transfer between different local parts of the entire system kept at a constant temperature. Very often, when modelling such processes using molecular dynamics simulations, thermostatting is done using strictly equilibrium approaches serving to describe the NVT ensemble. In this paper we critically discuss the applicability of such approaches to non-equilibrium problems, including those mentioned above, and stress that the correct temperature control can only be achieved if the method is based on the generalized Langevin equation (GLE). Specifically, we emphasize that a meaningful compromise between computational efficiency and a physically appropriate implementation of the NVT thermostat can be achieved, at least for solid state and surface problems, if the so-called stochastic boundary conditions (SBC), recently derived from the GLE (Kantorovich and Rompotis 2008 Phys. Rev. B 78 094305), are used. For SBC, the Langevin thermostat is only applied to the outer part of the simulated fragment of the entire system which borders the surrounding environment (not considered explicitly) serving as a heat bath. This point is illustrated by comparing the performance of the SBC and some of the equilibrium thermostats in two problems: (i) irradiation of the Si(001) surface with an energetic CaF(2) molecule using an ab initio density functional theory based method, and (ii) the tribology of two amorphous SiO(2) surfaces coated with self-assembled monolayers of methyl-terminated hydrocarbon alkoxylsilane molecules using a classical atomistic force field. We discuss the differences in behaviour of these systems due to applied thermostatting, and show that in some cases a qualitatively different physical behaviour of the simulated system can be obtained if an equilibrium thermostat is used.

Making tracks: electronic excitation roles in forming swift heavy ion tracks

Swift heavy ions cause material modification along their tracks, changes primarily due to their very dense electronic excitation. The available data for threshold stopping powers indicate two main classes of materials. Group I, with threshold stopping powers above about 10 keV nm(-1), includes some metals, crystalline semiconductors and a few insulators. Group II, with lower thresholds, comprises many insulators, amorphous materials and high T(c) oxide superconductors. We show that the systematic differences in behaviour result from different coupling of the dense excited electrons, holes and excitons to atomic (ionic) motions, and the consequent lattice relaxation. The coupling strength of excitons and charge carriers with the lattice is crucial. For group II, the mechanism appears to be the self-trapped exciton model of Itoh and Stoneham (1998 Nucl. Instrum. Methods Phys. Res. B 146 362): the local structural changes occur roughly when the exciton concentration exceeds the number of lattice sites. In materials of group I, excitons are not self-trapped and structural change requires excitation of a substantial fraction of bonding electrons, which induces spontaneous lattice expansion within a few hundred femtoseconds, as recently observed by laser-induced time-resolved x-ray diffraction of semiconductors. Our analysis addresses a number of experimental results, such as track morphology, the efficiency of track registration and the ratios of the threshold stopping power of various materials.

Temperature and pressure spikes in ion-beam cancer therapy

The inelastic thermal spike model is applied to liquid water in relation to high-energy 12C6+ beams (hundreds of MeV/u) used for cancer therapy. The goal of this project is to calculate the heat transfer in the vicinity of the incident-ion track. Thermal spike calculations indicate a very large temperature increase in the vicinity of ion tracks near the Bragg peak during the time interval from 10(-15) to 10(-9) s after the ion's passage and an increase in pressure, as large as tens of MPa, can be induced during that time. These effects suggest a possibility of thermomechanical pathways to disruption of irradiated DNA. An extension of the model for hydrogen, beryllium, argon, krypton, xenon, and uranium ions around the Bragg peak is presented as well.