The irradiation of carbon based nanostructures with ions and electrons has been shown to be an appropriate tool to tailor their properties. The defects induced in the nanostructures during irradiation are able to modify their mechanical and electronic properties. Here we simulate the irradiation of carbon nanotubes with carbon ions using a molecular dynamics code. We use the Tersoff potential joined smoothly to the Universal Ziegler-Biersack-Littmark potential at short distances. We study the number of defects produced after irradiation with a single carbon ion finding a saturation with its energy at ∼ 3 keV. We observe, after continuum irradiation with low energy ions, the formation of bumps in the irradiated region. For larger energy ions we find that the diameter of the nanotube shrinks as shown in previous works.This work has been financially supported by Project FIS2010-17225 by the Spanish Ministerio de Ciencia e Innovación

Denton Zanello, Cristian D.

Heredia-Avalos, Santiago

Moreno Marín, Juan Carlos

Irradiation of carbon nanotubes with carbon projectiles: A molecular dynamics study

Repositorio Institucional de la Universidad de Alicante

  Copyright line will be provided by the publisher pss-Header will be provided by the publisher   1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 Review copy – not for distribution  (pss-logo will be inserted here  by the publisher) Irradiation of carbon nanotubes with carbon projectiles. A molecular  dynamics study Cristian D. Denton*,1, Santiago Heredia-Avalos**,2, and Juan Carlos Moreno-Marín2 1 Departamento de Física Aplicada, Universidad de Alicante, E-03080, Spain 2 Departamento de Física, Ingeniería de Sistemas y Teoría de la Señal, Universidad de Alicante, E-03080, Spain Received ZZZ, revised ZZZ, accepted ZZZ Published online ZZZ                 (Dates will be provided by the publisher.) Keywords (Carbon nanotubes, Ion irradiation, Defect annealing, Computer Simulations)                               * Corresponding author: e-mail denton@ua.es, Phone: +34 965 903 400, Fax: +34 965 909 726 ** e-mail sheredia@ua.es, Phone: +34 965 903 400, Fax: +34 965 909 750   The irradiation of carbon based nanostructures with ions and electrons has been shown to be an appropriate tool to tailor their properties. The defects induced in the nanos-tructures during irradiation are able to modify their me-chanical and electronic properties. Here we simulate the irradiation of carbon nanotubes with carbon ions using a molecular dynamics code. We use the Tersoff potential joined smoothly to the Universal Ziegler-Biersack-Littmark potential at short distances. We study the num-ber of defects produced after irradiation with a single carbon ion finding a saturation with its energy at ~ 3 keV. We observe, after continuum irradiation with low energy ions, the formation of bumps in the irradiated re-gion. For larger energy ions we find that the diameter of the nanotube shrinks as shown in previous works.   Copyright line will be provided by the publisher  1 Introduction Since their discovery in 1991 [1] car-bon nanotubes (CNTs) have attracted a lot of attention due to their outstanding electronic and mechanical properties [2]. It has been shown that CNTs can be metals or narrow band semiconductors depending on their chirality. These properties, in addition to the fact that they have a nanom-eter size, make them ideal candidates for use in nanoelec-tronics [3]. CNTs have also extremely high Young modu-lus and strength, which are expected to be, respectively, about 19 and 56 times that of steel [4]. As a consequence of these excellent mechanical properties CNTs can be used for reinforcements of other materials [2,5]. Many other ex-citing applications for CNTs have emerged in the last years [2].  In order to use CNTs in these potential applications, it is necessary to develop effective methods to tailor their properties. A possible solution is the use of energetic parti-cles such as electrons, protons or ions [6]. Irradiation of CNTs with particles can implant foreign atoms, create de-fects, produce amorphous regions or recristallize the lat-tice. This technique is specially suited for manipulation of CNTs since the graphitic structure has a tendency to recon-struct the point defects by saturating dangling bonds [7]. The irradiation with ion or electron beams can also be used to coalesce [8], weld [9], link [10] or cut [11] the CNTs. This could allow the creation of molecular junctions and multiterminal heterostructures in future carbon-based nanoelectronics.  These technological applications need a thorough knowledge of the irradiation effects in the electrical and mechanical properties of CNT after irradiation. Thus it is important to analyze the mechanisms of production of de-fects under irradiation and also the changes in the structure induced by irradiation due to the saturation of the dangling bonds. Several papers have been devoted, both theoreti-cally and experimentally, to the study of the effects of the ionic irradiation in CNTs and related nanostructures [6-7, 12-16].  In this paper we study the production of defects after irradiation of CNTs with carbon ions, using molecular dy-namics simulations. We analyze the changes in the irradi-ated CNTs for different energies of the incident ion and dose. 2  Author, Author, and Author: Short title  Copyright line will be provided by the publisher 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 2 Simulation details We have developed a code based on classical molecular dynamics (MD) in order to simulate the bombardment of CNTs by carbon projectiles. Using this simulation code we have studied the defects produced on CNTs by carbon projectiles with incident en-ergies from 10 eV to 4 keV. Because of the low energies involved in this work, we have neglected the electronic stopping, being the colli-sional process governed by the nuclear interactions. In order to describe the C-C interactions, we have used the Tersoff potential [17], smoothly linked to the Ziegler-Biersack-Littmark (ZBL) universal potential at short in-teratomic distances [18] to describe realistically the hard collisions between energetic atoms. Therefore, the inter-atomic potential was evaluated according to:  (1) where F(r) is a Fermi function:  (2) being bf = 14 Å–1 and rf = 0.95 Å. This potencial has al-ready been used to study the ionic irradiation of graphite [12], as well as the B and N implantation in CNTs [15]. We have solved the equations of motion of all interact-ing atoms using the velocity Verlet algorithm [19] with a variable time step, which is stated considering that the maximum distance the atoms can move between two con-secutive time steps. This distance was chosen checking that our simulation code satisfies the law of conservation of en-ergy. It is worth to mention that, in order to avoid the dis-placement of the CNT when bombarding with carbon pro-jectiles, we have kept fixed the atoms at both ends of the CNT during the simulation.  On the other hand, we have used the Berendsen ther-mostat to control the temperature of the CNT [20]. In this scheme, we have scaled the velocity of some atoms, placed near the CNT ends, according to the scaling factor:  (3) where Δt is the time step of the simulation, τ is the time constant to control the bath coupling, T0 is the desired tem-perature, and T is the instantaneous temperature. 3 Results and discussion We have considered in our simulations single-walled carbon nanotubes (SWCNT) armchair (10, 10), ~100 Å in length and ~14 Å in diameter; both ends have fixed positions, whereas the velocity of the carbon atoms corresponding to 5% of the total length of the CNT next to its fixed ends are scaled according to the Berendsen thermostat, as described previously. It is worth to mention that the carbon projectiles bom-bard the CNT in the direction perpendicular to the nano-tube axis, although the exact impact coordinates are ran-domly distributed on a central region of 10% of the total length of the CNT, as shown in figure 1.   Figure 1 SWCNT with armchair (10, 10) structure, ~100 Å in length and ~14 Å in diameter. Red atoms at both ends of the CNTs have fixed positions during simulations, whereas Berend-sen thermostat was applied to light blue atoms next to the CNT ends. Carbon projectiles bombard the CNT centre in the direction perpendicular to the CNT axis and along 10% of the CNT length.  We have divided the carbon irradiation of CNTs in three stages. The first one reproduces the impact of the projectile in the direction perpendicular to the CNT axis. We simulate this stage during 20 ps and the temperature of the Berendsen thermostat is stated at 300 K. The second stage is the posterior heating of the CNT during 50 ps, be-ing 3000 K the temperature of the Berendsen thermostat. Finally, at the third stage we cool the CNT to a tempera-ture of  0 K during 15 ps. These stages can be observed on figure 2, where we have depicted the temperature of the CNT as a function of the simulation time after the irradi-ation with a 200 eV carbon projectile at 0.5 ps. Figure 2 CNT temperature as a function of simulation time after the irradiation with a 200 eV carbon projectile at 0.5 ps. The curves represent the CNT temperature during: irradiation stage (red), annealing stage (blue), and cooling stage (brown).  When carbon ions hit the CNTs wall, the main damage mechanism is the transfer of energy that results in the dis-pss-Header will be provided by the publisher 3   Copyright line will be provided by the publisher 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 placement of the atoms of the structure, and hence, the va-cancies generation. Furthermore, depending on the trans-ferred energy, the primary atoms can be backscattered, they can replace or extract atoms of the CNT, or even they can be added to the CNT (adatoms). We have studied two different situations in our simula-tions. Firstly, we have evaluated the damage produced in the CNT due to individual carbon ion bombardment; sec-ondly we have also analyzed the damage caused by consecutive carbon ions bombardment. 3.1 Single projectile irradiation Single-ion bom-bardment of the CNT causes damage, which can be evalu-ated through the formation of defects in the structure. Ac-cording to previous authors [14], we assume that the co-ordination defect number reflects quantitatively the irradi-ation damage produced. This magnitude is defined as the number of carbon atoms with coordination number differ-ent than three.  Figure 3 Number of coordination defects as a function of the in-cident projectile energy. After the collisional stage (black circles) and after the subsequent annealing stage (blue squares). The curves were depicted to guide the eye.  In figure 3 we present the results of the coordination defect number obtained after the irradiation with carbon ions of different incident energies (black circles) and its evolution after the annealing stage of the CNT (blue squares); the curves were depicted to guide the eye. We can see that the coordination defect number increases with incident ion energy up to 3 keV, reaching a saturation value ~15 at higher energies. In addition, figure 3 shows that after the irradiation stage, the heating of the CNT in the annealing stage helps the migration of interstitial atoms and vacancies and their recombination, so that the coordi-nation defect number decreases. The irradiation of low energy carbon ions causes the adhesion of carbon atoms to the CNT wall (single ada-toms), which evolved after annealing to the formation of double-coordinated atoms (bridge defects). Figure 4(a) shows a typical transformation of an adatom to a bridge defect after annealing, when a 10 eV carbon projectile hits on the CNT during the irradiation stage. However, for higher projectile energies, the irradiation also causes the formation of vacancies, resulting in different size holes that can evolve after the annealing stage to non-hexagonal structures, called Stone-Wales defects. We show in figure 4(b) the typical transformation of two vacancies into a Stone-Wales defect after annealing; in this case the projec-tile hits the CNT with 100 eV. In both figures the colours indicate the coordination number: white (1), red (2), blue (3), and green (4).   Figure 4 Typical transformation of defects during annealing af-ter irradiation with: (a) 10 eV and (b) 100 eV carbon ion. Colours mean the coordination number: white (1), red (2), blue (3), green (4).  3.2 Continuum projectile irradiation On the other hand, we have studied the damage produced on the CNT when it is successively irradiated with carbon ions.   Figure 5 Coordination defect number as a function of the num-ber of ions (ion dose) that bombard the CNT with 100 eV energy. 4  Author, Author, and Author: Short title  Copyright line will be provided by the publisher 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 We have depicted in figure 5 the coordination defect number as a function of the number of incident carbon ions; the energy of all carbon projectiles is 100 eV. As ex-pected, we have found in our simulations that the coordina-tion defect number grows proportionally to the incident ionic dose. In addition, we have observed that the coordi-nation defect number increases with the projectile energy.   Figure 6 Appearance of CNT after irradiation with 50 consecu-tive carbon ions with different energies.  The effects of the continuous irradiation of the CNT, depending on the energy of the incident projectile, are shown in figure 6, in which we present the appearance of the CNT after irradiation with 50 projectiles; colours rep-resent the coordination number of the carbon atoms as pre-viously stated. With low energy ions (~10 eV), adatoms are produced in the walls of the CNT, which can evolve into the forma-tion of bumps on the irradiated region, as observed in fig-ure 6. When the projectiles have higher energies (100-1000 eV), vacancies appear during the irradiation stage, which evolve, during the annealing, into polygonal shapes. In ad-dition, a progressive reduction of the CNT diameter is ob-served, as shown in figure 6. Increasing the dose results in a higher damage, which could even break the CNT in the irradiated area. 4 Conclusions We have developed a MD code in order to simulate the irradiation of nanostructures. In this work, we use this code to study the irradiation of CNTs with carbon ions. We have evaluated the damage produced by individual ions in a perfect CNT. We find that for low projectile en-ergies (E < 50 eV) the defects are only adatoms, while va-cancies or cluster of vacancies are also observed for higher energies. These defects are partly healed after annealing through the saturation of dangling bonds. The number of defects generated by a single ion grows with projectile en-ergy but saturates at ~ 3 keV. We have also studied the effect of continuum irradi-ation of the CNT. For low energies the accumulation of adatoms produces a bump in the irradiated region. At intermediate energies the irradiation produces vacancies which are healed through Stone-Wales defects. This gives rise to a shrinking of the CNT diameter in the irradiated region which has been observed before [13]. If the projec-tile energy is high enough the continuum irradiation pro-duces the breaking of the CNT. Acknowledgements This work has been financially sup-ported by Project FIS2010-17225 by the Spanish Ministerio de Ciencia e Innovación. The authors would like to thank M. J. Ca-turla for her help and suggestions. References [1] S. IIjima, Nature 354, 56 (1991). [2] M. S. Dresselhaus. G. Dresselhaus, P. Avouris (Eds.), Car-bon nanotubes, Synthesis, Structure, Properties and Appli-cations (Springer, Berlin, 2001).  [3] P. L. Mc Euen, Nature 393, 15 (1998). [4] R. H. Baughman, A. A. Zakhidov, W. A. de Heer, Science 297, 787 (2002). [5] A. B. Dalton, S. Collins, E. MuÃ±oz, J. M. Razal, V. H. Eb-ron, J. P. Ferraris, J. N. Coleman, B. G. Kim, R. H. Baughman, Nature 703, 423 (2003). [6] A. V. Krasheninnikov and K. Nordlund, J. Appl. Phys. 107, 071301 (2010). [7] A. V. Krasheninnikov, F. Banhart, Nat. Mater. 6, 723 (2007). [8] M. Terrones, H. Terrones, F. Banhart, J.-C. Charlier, P. M. Ajayan, Science 288, 1226 (2000). [9] M. Terrones, F. Banhart, N. Grovert, J.-C. Charlier, H. Ter-rones, P. M. Ajayan, Phys. Rev. Lett. 89, 075505 (2002). [10] H. Stahl. J. Appenzeller, R. Martel, P. Avouris, B. Lenge-ler, Phys. Rev. Lett. 85, 5186 (2000). [11] M. S. Raghuveer, P. G. Ganesan, J. D'Arcy-Gall, G. Ra-manath, M. Marshall, I. Petrov, Appl. Phys. Lett. 84, 4484 (2004). [12] K. Nordlund, J. Keinonen, T. Mattila, Phys. Rev. Lett. 77, 699 (1996). [13] A. V. Krasheninnikov, K. Nordlund and J. Keinonen, Phys. Rev. B 65, 165423 (2002). [14] A. V. Krasheninnikov and K. Nordlund, Nucl. Instrum. Meth. Phys. Res. B 216, 355 (2004). [15] J. Kotakoski, A. V. Krasheninnikov, Yuchen Ma, A. S. Fo-ster, K. Nordlund, R. M. Nieminen, Phys. Rev. B 71, 205408 (2005). [16] A. V. Krasheninnikov, J. Comp. Theor. Nanoscience 5, 1 (2008). [17] J. Tersoff, Phys. Rev. B 39, 5566 (1989). [18] J. F. Ziegler, J. P. Biersack and M. D. Ziegler, SRIM: The stopping and range of ions in matter (Lulu Press, Napa, 2008).  [19] M. P. Allen y D. J. Tildesley, Computer simulation of Li-quids (Oxford Science Publications, Clarendon Press, Ox-ford, 2002). [20] H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak, J. Chem. Phys. 81, 3684 (1984). 

Irradiation of carbon nanotubes with carbon projectiles: A molecular dynamics study

Abstract

Similar works

Full text

Available Versions

RUA

Repositorio Institucional de la Universidad de Alicante