In condensed-matter physics, a collision cascade (also known as a displacement cascade or a displacement spike) is a set of nearby adjacent energetic (much higher than ordinary thermal energies) collisions of atoms induced by an energetic particle in a solid or liquid. [1] [2]
If the maximum atom or ion energies in a collision cascade are higher than the threshold displacement energy of the material (tens of eVs or more), the collisions can permanently displace atoms from their lattice sites and produce defects. The initial energetic atom can be, e.g., an ion from a particle accelerator, an atomic recoil produced by a passing high-energy neutron, electron or photon, or be produced when a radioactive nucleus decays and gives the atom a recoil energy.
The nature of collision cascades can vary strongly depending on the energy and mass of the recoil/incoming ion and density of the material (stopping power).
When the initial recoil/ion mass is low, and the material where the cascade occurs has a low density (i.e. the recoil-material combination has a low stopping power), the collisions between the initial recoil and sample atoms occur rarely, and can be understood well as a sequence of independent binary collisions between atoms. This kind of a cascade can be theoretically well treated using the binary collision approximation (BCA) simulation approach. For instance, H and He ions with energies below 10 keV can be expected to lead to purely linear cascades in all materials.
The most commonly used BCA code SRIM [3] can be used to simulate linear collision cascades in disordered materials for all ion in all materials up to ion energies of 1 GeV. Note, however, that SRIM does not treat effects such as damage due to electronic energy deposition or damage produced by excited electrons. The nuclear and electronic stopping powers used are averaging fits to experiments, and are thus not perfectly accurate either. The electronic stopping power can be readily included in binary collision approximation [4] or molecular dynamics (MD) simulations. In MD simulations they can be included either as a frictional force [5] [6] [7] [8] [9] [10] [11] [12] or in a more advanced manner by also following the heating of the electronic systems and coupling the electronic and atomic degrees of freedom. [13] [14] [15] However, uncertainties remain on what is the appropriate low-energy limit of electronic stopping power or electron-phonon coupling is. [12] [16]
In linear cascades the set of recoils produced in the sample can be described as a sequence of recoil generations depending on how many collision steps have passed since the original collision: primary knock-on atoms (PKA), secondary knock-on atoms (SKA), tertiary knock-on atoms (TKA), etc. Since it is extremely unlikely that all energy would be transferred to a knock-on atom, each generation of recoil atoms has on average less energy than the previous, and eventually the knock-on atom energies go below the threshold displacement energy for damage production, at which point no more damage can be produced.
When the ion is heavy and energetic enough, and the material is dense, the collisions between the ions may occur so near to each other that they can not be considered independent of each other. In this case the process becomes a complicated process of many-body interactions between hundreds and tens of thousands of atoms, which can not be treated with the BCA, but can be modelled using molecular dynamics methods. [1] [17]
Typically, a heat spike is characterized by the formation of a transient underdense region in the center of the cascade, and an overdense region around it. [1] [18] After the cascade, the overdense region becomes interstitial defects, and the underdense region typically becomes a region of vacancies.
If the kinetic energy of the atoms in the region of dense collisions is recalculated into temperature (using the basic equation E = 3/2·N·kBT), one finds that the kinetic energy in units of temperature is initially of the order of 10,000 K. Because of this, the region can be considered to be very hot, and is therefore called a heat spike or thermal spike (the two terms are usually considered to be equivalent). The heat spike cools down to the ambient temperature in 1–100 ps, so the "temperature" here does not correspond to thermodynamic equilibrium temperature. However, it has been shown that after about 3 lattice vibrations, the kinetic energy distribution of the atoms in a heat spike has the Maxwell–Boltzmann distribution, [19] making the use of the concept of temperature somewhat justified. Moreover, experiments have shown that a heat spike can induce a phase transition which is known to require a very high temperature, [20] showing that the concept of a (non-equilibrium) temperature is indeed useful in describing collision cascades.
In many cases, the same irradiation condition is a combination of linear cascades and heat spikes. For example, 10 MeV Cu ions bombarding Cu would initially move in the lattice in a linear cascade regime, since the nuclear stopping power is low. But once the Cu ion would slow down enough, the nuclear stopping power would increase and a heat spike would be produced. Moreover, many of the primary and secondary recoils of the incoming ions would likely have energies in the keV range and thus produce a heat spike.
For instance, for copper irradiation of copper, recoil energies of around 5–20 keV are almost guaranteed to produce heat spikes. [21] [22] At lower energies, the cascade energy is too low to produce a liquid-like zone. At much higher energies, the Cu ions would most likely lead initially to a linear cascade, but the recoils could lead to heat spikes, as would the initial ion once it has slowed down enough. The concept subcascade breakdown threshold energy signifies the energy above which a recoil in a material is likely to produce several isolated heat spikes rather than a single dense one.
Computer simulation-based animations of collision cascades in the heat spike regime are available on YouTube. [23]
Swift heavy ions, i.e. MeV and GeV heavy ions which produce damage by a very strong electronic stopping, can also be considered to produce thermal spikes [24] [25] in the sense that they lead to strong lattice heating and a transient disordered atom zone. However, at least the initial stage of the damage might be better understood in terms of a Coulomb explosion mechanism. [26] Regardless of what the heating mechanism is, it is well established that swift heavy ions in insulators typically produce ion tracks forming long cylindrical damage zones [24] [27] of reduced density. [28] [29]
To understand the nature of collision cascade, it is very important to know the associated time scale. The ballistic phase of the cascade, when the initial ion/recoil and its primary and lower-order recoils have energies well above the threshold displacement energy, typically lasts 0.1–0.5 ps. If a heat spike is formed, it can live for some 1–100 ps until the spike temperature has cooled down essentially to the ambient temperature. [30] The cooling down of the cascade occurs via lattice heat conductivity and by electronic heat conductivity after the hot ionic subsystem has heated up the electronic one via electron–phonon coupling. Unfortunately the rate of electron-phonon coupling from the hot and disordered ionic system is not well known, as it can not be treated equally to the fairly well known process of transfer of heat from hot electrons to an intact crystal structure. [31] Finally, the relaxation phase of the cascade, when the defects formed possibly recombine and migrate, can last from a few ps to infinite times, depending on the material, its defect migration and recombination properties, and the ambient temperature.
Since the kinetic energies in a cascade can be very high, it can drive the material locally far outside thermodynamic equilibrium. Typically this results in defect production. The defects can be, e.g., point defects such as Frenkel pairs, ordered or disordered dislocation loops, stacking faults, [32] or amorphous zones. [33] Prolonged irradiation of many materials can lead to their full amorphization, an effect which occurs regularly during the ion implantation doping of silicon chips. [34]
The defects production can be harmful, such as in nuclear fission and fusion reactors where the neutrons slowly degrade the mechanical properties of the materials, or a useful and desired materials modification effect, e.g., when ions are introduced into semiconductor quantum well structures to speed up the operation of a laser. [35] or to strengthen carbon nanotubes. [36]
A curious feature of collision cascades is that the final amount of damage produced may be much less than the number of atoms initially affected by the heat spikes. Especially in pure metals, the final damage production after the heat spike phase can be orders of magnitude smaller than the number of atoms displaced in the spike. [1] On the other hand, in semiconductors and other covalently bonded materials the damage production is usually similar to the number of displaced atoms. [1] [22] Ionic materials can behave like either metals or semiconductors with respect to the fraction of damage recombined. [37]
Collision cascades in the vicinity of a surface often lead to sputtering, both in the linear spike and heat spike regimes. [21] Heat spikes near surfaces also frequently lead to crater formation. [38] [39] This cratering is caused by liquid flow of atoms, [40] but if the projectile size above roughly 100,000 atoms, the crater production mechanism switches to the same mechanism as that of macroscopic craters produced by bullets or asteroids. [41]
The fact that many atoms are displaced by a cascade means that ions can be used to deliberately mix materials, even for materials that are normally thermodynamically immiscible. This effect is known as ion beam mixing. [42]
The non-equilibrium nature of irradiation can also be used to drive materials out of thermodynamic equilibrium, and thus form new kinds of alloys. [43]
A crystallographic defect is an interruption of the regular patterns of arrangement of atoms or molecules in crystalline solids. The positions and orientations of particles, which are repeating at fixed distances determined by the unit cell parameters in crystals, exhibit a periodic crystal structure, but this is usually imperfect. Several types of defects are often characterized: point defects, line defects, planar defects, bulk defects. Topological homotopy establishes a mathematical method of characterization.
In physics, sputtering is a phenomenon in which microscopic particles of a solid material are ejected from its surface, after the material is itself bombarded by energetic particles of a plasma or gas. It occurs naturally in outer space, and can be an unwelcome source of wear in precision components. However, the fact that it can be made to act on extremely fine layers of material is utilised in science and industry—there, it is used to perform precise etching, carry out analytical techniques, and deposit thin film layers in the manufacture of optical coatings, semiconductor devices and nanotechnology products. It is a physical vapor deposition technique.
Ionization is the process by which an atom or a molecule acquires a negative or positive charge by gaining or losing electrons, often in conjunction with other chemical changes. The resulting electrically charged atom or molecule is called an ion. Ionization can result from the loss of an electron after collisions with subatomic particles, collisions with other atoms, molecules and ions, or through the interaction with electromagnetic radiation. Heterolytic bond cleavage and heterolytic substitution reactions can result in the formation of ion pairs. Ionization can occur through radioactive decay by the internal conversion process, in which an excited nucleus transfers its energy to one of the inner-shell electrons causing it to be ejected.
Neutron radiation is a form of ionizing radiation that presents as free neutrons. Typical phenomena are nuclear fission or nuclear fusion causing the release of free neutrons, which then react with nuclei of other atoms to form new nuclides—which, in turn, may trigger further neutron radiation. Free neutrons are unstable, decaying into a proton, an electron, plus an electron antineutrino. Free neutrons have a mean lifetime of 887 seconds.
A Rydberg atom is an excited atom with one or more electrons that have a very high principal quantum number, n. The higher the value of n, the farther the electron is from the nucleus, on average. Rydberg atoms have a number of peculiar properties including an exaggerated response to electric and magnetic fields, long decay periods and electron wavefunctions that approximate, under some conditions, classical orbits of electrons about the nuclei. The core electrons shield the outer electron from the electric field of the nucleus such that, from a distance, the electric potential looks identical to that experienced by the electron in a hydrogen atom.
Imperfections in the crystal lattice of diamond are common. Such defects may be the result of lattice irregularities or extrinsic substitutional or interstitial impurities, introduced during or after the diamond growth. The defects affect the material properties of diamond and determine to which type a diamond is assigned; the most dramatic effects are on the diamond color and electrical conductivity, as explained by the electronic band structure.
In condensed-matter physics, channelling (or channeling) is the process that constrains the path of a charged particle in a crystalline solid.
In nuclear and materials physics, stopping power is the retarding force acting on charged particles, typically alpha and beta particles, due to interaction with matter, resulting in loss of particle kinetic energy. Stopping power is also interpreted as the rate at which a material absorbs the kinetic energy of a charged particle. Its application is important in a wide range of thermodynamic areas such as radiation protection, ion implantation and nuclear medicine.
Car–Parrinello molecular dynamics or CPMD refers to either a method used in molecular dynamics or the computational chemistry software package used to implement this method.
In materials science, an interstitial defect is a type of point crystallographic defect where an atom of the same or of a different type, occupies an interstitial site in the crystal structure. When the atom is of the same type as those already present they are known as a self-interstitial defect. Alternatively, small atoms in some crystals may occupy interstitial sites, such as hydrogen in palladium. Interstitials can be produced by bombarding a crystal with elementary particles having energy above the displacement threshold for that crystal, but they may also exist in small concentrations in thermodynamic equilibrium. The presence of interstitial defects can modify the physical and chemical properties of a material.
Ion beam mixing is the atomic intermixing and alloying that can occur at the interface separating two different materials during ion irradiation. It is applied as a process for adhering two multilayers, especially a substrate and deposited surface layer. The process involves bombarding layered samples with doses of ion radiation in order to promote mixing at the interface, and generally serves as a means of preparing electrical junctions, especially between non-equilibrium or metastable alloys and intermetallic compounds. Ion implantation equipment can be used to achieve ion beam mixing.
In materials science, the threshold displacement energy is the minimum kinetic energy that an atom in a solid needs to be permanently displaced from its site in the lattice to a defect position. It is also known as "displacement threshold energy" or just "displacement energy". In a crystal, a separate threshold displacement energy exists for each crystallographic direction. Then one should distinguish between the minimum and average over all lattice directions' threshold displacement energies. In amorphous solids, it may be possible to define an effective displacement energy to describe some other average quantity of interest. Threshold displacement energies in typical solids are of the order of 10-50 eV.
Stopping and Range of Ions in Matter (SRIM) is a group of computer programs which calculate interactions between ions and matter; the core of SRIM is a program called Transport of Ions in Matter (TRIM). SRIM is popular in the ion implantation research and technology community, and also used widely in other branches of radiation material science.
Swift heavy ions are the components of a type of particle beam with high enough energy that electronic stopping dominates over nuclear stopping. They are accelerated in particle accelerators to very high energies, typically in the MeV or GeV range and have sufficient energy and mass to penetrate solids on a straight line. In many solids swift heavy ions release sufficient energy to induce permanently modified cylindrical zones, so-called ion tracks. If the irradiation is carried out in an initially crystalline material, ion tracks consist of an amorphous cylinder. Ion tracks can be produced in many amorphizing materials, but not in pure metals, where the high electronic heat conductivity dissipates away the electronic heating before the ion track has time to form.
In condensed-matter physics, the binary collision approximation (BCA) is a heuristic used to more efficiently simulate the penetration depth and defect production by energetic ions in solids. In the method, the ion is approximated to travel through a material by experiencing a sequence of independent binary collisions with sample atoms (nuclei). Between the collisions, the ion is assumed to travel in a straight path, experiencing electronic stopping power, but losing no energy in collisions with nuclei.
Radiation materials science is a subfield of materials science which studies the interaction of radiation with matter: a broad subject covering many forms of irradiation and of matter.
In condensed-matter physics, a primary knock-on atom (PKA) is an atom that is displaced from its lattice site by irradiation; it is, by definition, the first atom that an incident particle encounters in the target. After it is displaced from its initial lattice site, the PKA can induce the subsequent lattice site displacements of other atoms if it possesses sufficient energy, or come to rest in the lattice at an interstitial site if it does not.
Double ionization is a process of formation of doubly charged ions when laser radiation is exerted on neutral atoms or molecules. Double ionization is usually less probable than single-electron ionization. Two types of double ionization are distinguished: sequential and non-sequential.
Interatomic potentials are mathematical functions to calculate the potential energy of a system of atoms with given positions in space. Interatomic potentials are widely used as the physical basis of molecular mechanics and molecular dynamics simulations in computational chemistry, computational physics and computational materials science to explain and predict materials properties. Examples of quantitative properties and qualitative phenomena that are explored with interatomic potentials include lattice parameters, surface energies, interfacial energies, adsorption, cohesion, thermal expansion, and elastic and plastic material behavior, as well as chemical reactions.