Pith. sign in

REVIEW 126 references

Repulsive interatomic potentials calculated at three levels of theory

Not yet reviewed by Pith; the record is open.

This paper has not been read by Pith yet. Machine review is queued; the pith claim, tier, and objections will appear here once it completes.

SPECIMEN: schema-true, not a live event

T0 review · schema-true

One-sentence machine reading of the paper's core claim.

pith:XXXXXXXX · record.json · timestamp

arxiv 2501.06617 v3 pith:YKSLSZAK submitted 2025-01-11 physics.comp-ph

classification physics.comp-ph
keywords potentialscalculationsinteratomicrepulsivetheorydatadensitymodel
verification ladder T0 review T1 audit T2 compute T3 formal

Signed reviews

No signed human review yet.

0 comments
abstract

The high-energy repulsive interaction between nuclei at distances much smaller than the equilibrium bond length is the key quantity determining the nuclear stopping power and atom scattering in keV and MeV radiation events. This interaction is traditionally modeled within orbital-free density functional theory with frozen atomic electron densities, following the Ziegler-Biersack-Littmark (ZBL) model. In this work, we calculate atom pair specific repulsive interatomic potentials with the ZBL model, and compare them to two kinds of quantum chemical calculations - second-order M{\o}ller-Plesset perturbation theory in flexible Gaussian basis sets as well as density functional theory with numerical atomic orbital basis sets - which go well beyond the limitations in the ZBL model, allowing the density to relax in the calculations. We show that the repulsive interatomic potentials predicted by the two quantum chemical models agree within $\sim$ 1% for potential energies above 30 eV, while the ZBL pair-specific potentials and universal ZBL potentials differ much more from either of these calculations. We provide new pair-specific fits of the screening functions in terms of 3 exponentials to the calculations for all pairs $Z_1$-$Z_2$ for $1 \leq Z_i \leq 92$, and show that they agree within $\sim 2$% with the raw data. We use the new potentials to simulate ion implantation depth profiles in single crystalline Si and show very good agreement with experiment. However, we also show that under channeling conditions, the attractive part of the potential can affect the depth profiles. The full data sets of all the calculated interatomic potentials as well as analytic fits to the data are shared as open access.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

126 extracted references · 79 canonical work pages

  1. [1]

    Often the thickness of the surface oxide layers is not known, and in many ex- periments the implantation tilt and twist angles are not accurately controlled or reported

    Simulation setup Direct comparison of range profiles with experiments is complicated by several factors. Often the thickness of the surface oxide layers is not known, and in many ex- periments the implantation tilt and twist angles are not accurately controlled or reported. The implantation pro- cess may also amorphize semiconductor samples. More- over, i...

  2. [2]

    The results of this comparison are given in fig

    DMol reference We also carried out a similar comparison by using the DMol data for all element pairs up to Z = 92 as reference, as MP2 data is not available for the heavier elements. The results of this comparison are given in fig. 10 and table II. Similarly to the differences to MP2 discussed for lower- Z, both the ZBL universal and ZBL pair-specific pot...

  3. [3]

    muffin-tin

    Comparison of results with experiments The results of the above approach are shown in fig. 18 for As-Si implantation, in fig. 19 for B-Si implantation, and in fig. 20 for P-Si implantation. To put the results into context, a range profile calculated with the com- monly used SRIM code [101, 102] that does not include channeling are also shown in the figure...

  4. [4]

    1609 ± 2 1620 ± 2 1617 ± 1 1630 ± 2 1710 ± 2 H Si 10 Non-channeled 1753 ± 2 1751 ± 2 1758 ± 2 1750 ± 2 1750 ± 2

  5. [5]

    3123 ± 2 3212 ± 2 3310 ± 2 3307 ± 2 3319 ± 2 Fe Fe 100 Non-channeled – 314 ± 1 312 ± 1 292 ± 1 313 ± 1

  6. [6]

    Nordlund, S

    K. Nordlund, S. J. Zinkle, A. E. Sand, F. Granberg, R. S. Averback, R. Stoller, T. Suzudo, L. Malerba, F. Banhart, W. J. Weber, F. Willaime, S. Dudarev, and D. Simeone, Primary radiation damage: a review of current understanding and models, J. Nucl. Mater. 512, 450 (2018)

  7. [7]

    18: Range profiles for As ions implanted into Si compared with experimental data [76]

    Effect of oxide thickness To examine how sensitive the results are to the oxide layer thickness, we ran the case of 15 keV As implanta- tion into the 100 channel and 100 keV As implantation into the non-channeling thickness with only 1 nm oxide 10 16 10 17 10 18 Concentration (ions/cm 3 ) 0 500 1000 1500 2000 2500 3000 3500 Depth (Å) MD, DMol MD, NLH MD, ...

  8. [8]

    Bohr, The penetration of atomic particles through matter, Mat

    N. Bohr, The penetration of atomic particles through matter, Mat. Fys. Medd. Dan. Vid. Selsk. 18, 1 (1948)

Show all 126 references
  1. [9]

    A. R. Leach, Molecular modelling: principles and appli- cations, 2nd ed. (Pearson Education, Harlow, England, 2001)

  2. [10]

    Calculate the screening function ϕ on an equidis- tant grid of R values with a spacing of 0.16 ˚A, start- ing from R = 0, so that the last value of R does not exceed Rmax−0.04 ˚A, where Rmax = r1max + r2max

  3. [11]

    Insert the new data point into the list of ϕ(R) values

    Bisect each interval, calculate ϕ, and evaluate the interpolation error when evaluating the spline through all previously calculated ϕ(R) values at the new point. Insert the new data point into the list of ϕ(R) values. If the interpolation error is below the tolerance, mark th...

  4. [12]

    If there are more than two intervals of equal length in a row, form pairs of adjacent intervals and consider only points in 24 the centers of the pairs

    Starting from the smallest spacing, iterate over all spacings ∆R, and determine points between inter- vals of equal length ∆ R. If there are more than two intervals of equal length in a row, form pairs of adjacent intervals and consider only points in 24 the centers of the pai...

  5. [13]

    G. S. Was, Fundamentals of Radiation Materials Science (Springer, Berlin, 2012)

  6. [14]

    Bernas, ed., Materials Science with Ion Beams , Top- ics in Applied Physics, Vol

    H. Bernas, ed., Materials Science with Ion Beams , Top- ics in Applied Physics, Vol. 116 (Springer, 2010)

  7. [15]

    J. S. Williams, Ion implantation of semiconductors, Mater. Sci. Eng. A 253, 8 (1998)

  8. [16]

    Seeger, The Nature of Radiation Damage in Met- als, in Radiation Damage in Solids, Vol

    A. Seeger, The Nature of Radiation Damage in Met- als, in Radiation Damage in Solids, Vol. 1 (International Atomic Energy Agency, Vienna, 1962) pp. 101–127

  9. [17]

    Lindhard, M

    J. Lindhard, M. Scharff, and H. E. Schiøtt, Range con- cepts and heavy ion ranges, Kgl. Danske Vid. Selskab, Mat.- Fys. Medd. 33 (14), 1 (1963)

  10. [18]

    P. M. Echenique, R. M. Nieminen, and R. H. Ritchie, Density functional calculation of stopping power of an electron gas for slow ions, Solid State Comm. 37, 779 (1981)

  11. [19]

    Nordlund, Historical review of computer simulation of radiation effects in materials, J

    K. Nordlund, Historical review of computer simulation of radiation effects in materials, J. Nucl. Mater.520, 273 (2019), Invited review in Diamond Anniversary issue

  12. [20]

    Kittel, Introduction to Solid State Physics , 3rd ed

    C. Kittel, Introduction to Solid State Physics , 3rd ed. (John Wiley & Sons, New York, 1968)

  13. [21]

    Nordlund, D

    K. Nordlund, D. Sundholm, P. Pyykk¨ o, D. Martinez Zambrano, and F. Djurabekova, Nuclear stopping power of antiprotons, Phys. Rev. A 96, 042717 (2017)

  14. [22]

    M. T. Robinson and I. M. Torrens, Computer Simu- lation of atomic-displacement cascades in solids in the binary-collision approximation, Phys. Rev. B 9, 5008 (1974)

  15. [23]

    Moli` ere, Theorie der Streuung schneller geladener Teilchen .I

    G. Moli` ere, Theorie der Streuung schneller geladener Teilchen .I. Einzelstreuung am abgeschirmten Coulomb- Feld, Z. Naturforsch. A 2, 133 (1947)

  16. [24]

    Eckstein, In Ref

    W. Eckstein, In Ref. [114], eq. (4.4.3) on p. 55

  17. [25]

    J. F. Ziegler, J. P. Biersack, and U. Littmark, The Stopping and Range of Ions in Matter (Pergamon, New York, 1985)

  18. [26]

    A. N. Zinoviev, Electron screening of the Coulomb po- tential at small internuclear distances, Nuclear Instru- ments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 354, 308 (2015)

  19. [27]

    Born and R

    M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln, Ann. Phys. (Berlin) 389, 457 (1927)

  20. [28]

    Fermi and R

    E. Fermi and R. D. Richtmyer, Note on census-taking in Monte Carlo calculations, A declassified report by Enrico Fermi. From the Los Alamos Archive. (1948)

  21. [29]

    L. M. Kishinevskii, Cross sections for inelastic atomic collisions, Bull. Acad. Sci. USSR, Phys. Ser. 26, 1433 (1962)

  22. [30]

    Bloch, Bemerkung zur Elektronentheorie des Fer- romagnetismus und der elektrischen Leitf¨ ahigkeit, Z

    F. Bloch, Bemerkung zur Elektronentheorie des Fer- romagnetismus und der elektrischen Leitf¨ ahigkeit, Z. Phys. 57, 545 (1929)

  23. [31]

    Schiwietz and P

    G. Schiwietz and P. L. Grande, Introducing electron capture into the unitary-convolution-approximation energy-loss theory at low velocities, Phys. Rev. A 84, 052703 (2011)

  24. [32]

    W. Mi, K. Luo, S. B. Trickey, and M. Pavanello, Orbital- free density functional theory: An attractive electronic structure method for large-scale first-principles simula- tions, Chem. Rev. 123, 12039 (2023)

  25. [33]

    Alml¨ of, K

    J. Alml¨ of, K. Faegri, and K. Korsell, Principles for a di- rect SCF approach to LCAO-MO ab-initio calculations, J. Comput. Chem. 3, 385 (1982)

  26. [34]

    W. D. Wilson, L. G. Haggmark, and J. P. Biersack, Calculations of nuclear stopping, ranges, and straggling in the low-energy region, Phys. Rev. B 15, 2458 (1977), publisher: American Physical Society

  27. [35]

    Jensen, Die Ladungsverteilung in Ionen und die Git- terkonstante des Rubidiumbromids nach der statistis- chen Methode, Z

    H. Jensen, Die Ladungsverteilung in Ionen und die Git- terkonstante des Rubidiumbromids nach der statistis- chen Methode, Z. Physik 77, 722 (1932)

  28. [36]

    Sommerfeld, Asymptotische Integration der Differen- tialgleichung des Thomas-Fermischen Atoms, Z

    A. Sommerfeld, Asymptotische Integration der Differen- tialgleichung des Thomas-Fermischen Atoms, Z. Physik 78, 283 (1932)

  29. [37]

    J. P. Biersack and J. F. Ziegler, Refined universal po- tentials in atomic collisions, Nucl. Instr. Meth. 194, 93 (1982)

  30. [38]

    S. T. Nakagawa and Y. Yamamura, Interatomic poten- tial in solids and its application to range calculations, 25 Radiat. Eff. 105, 239 (1988)

  31. [39]

    A. N. Zinoviev, Interaction potentials for modeling of ion-surface scattering, Nucl. Instrum. Meth. Phys. Res. B 269, 829 (2011)

  32. [40]

    L. H. Thomas, The calculation of atomic fields, Math. Proc. Cambridge Philos. Soc. 23, 542 (1927)

  33. [41]

    Fermi, Un metodo statistico per la determinazione di alcune propriet´ a dell’atomo., Rendiconti dell’Accademia Nazionale dei Lincei 6, 602 (1927)

    E. Fermi, Un metodo statistico per la determinazione di alcune propriet´ a dell’atomo., Rendiconti dell’Accademia Nazionale dei Lincei 6, 602 (1927)

  34. [42]

    D. Hait, N. M. Tubman, D. S. Levine, K. B. Whaley, and M. Head-Gordon, What Levels of Coupled Cluster Theory Are Appropriate for Transition Metal Systems? A Study Using Near-Exact Quantum Chemical Values for 3d Transition Metal Binary Compounds, J. Chem. Theory Comput. 15, 5370 (2019)

  35. [43]

    P. A. M. Dirac, Note on exchange phenomena in the Thomas atom, Math. Proc. Cambridge Philos. Soc. 26, 376 (1930)

  36. [44]

    Froese Fischer, G

    C. Froese Fischer, G. Gaigalas, P. J¨ onsson, and J. Biero´ n, GRASP2018—A Fortran 95 version of the General Relativistic Atomic Structure Package, Com- puter Physics Communications 237, 184 (2019)

  37. [45]

    Schiffmann, J

    S. Schiffmann, J. G. Li, J. Ekman, G. Gaigalas, M. Godefroid, P. J¨ onsson, and J. Biero´ n, Relativistic radial electron density functions and natural orbitals from GRASP2018, Computer Physics Communications 278, 108403 (2022)

  38. [46]

    J. H. Van Lenthe, R. Zwaans, H. J. J. Van Dam, and M. F. Guest, Starting SCF calculations by superposition of atomic densities., J. Comput. Chem. 27, 926 (2006)

  39. [47]

    Lehtola, Assessment of initial guesses for self- consistent field calculations

    S. Lehtola, Assessment of initial guesses for self- consistent field calculations. superposition of atomic po- tentials: Simple yet efficient, J. Chem. Theory Comput. 15, 1593 (2019), arXiv:1810.11659

  40. [48]

    N. H. Sabelli, R. Benedek, and T. L. Gilbert, Ground- state potential curves for Al2 and Al2 6+ in the repulsive region, Phys. Rev. A 20, 677 (1979)

  41. [49]

    Lehtola, Accurate reproduction of strongly repul- sive interatomic potentials, Phys

    S. Lehtola, Accurate reproduction of strongly repul- sive interatomic potentials, Phys. Rev. A 101, 032504 (2020), arXiv:1912.12624

  42. [50]

    R. K. Pathak and A. J. Thakkar, Very short-range in- teratomic potentials, J. Chem. Phys. 87, 2186 (1987)

  43. [51]

    N. H. Sabelli, M. Kantor, R. Benedek, and T. L. Gilbert, SCF potential curves for AlH and AlH+ in the attractive and repulsive regions, J. Chem. Phys. 68, 2767 (1978)

  44. [52]

    R. S. Mulliken, The interpretation of band spectra part III. electron quantum numbers and states of molecules and their atoms, Rev. Mod. Phys. 4, 1 (1932)

  45. [53]

    Hirao, Correlation diagram approach as a tool for in- terpreting chemistry: an introductory overview, WIREs Comput

    H. Hirao, Correlation diagram approach as a tool for in- terpreting chemistry: an introductory overview, WIREs Comput. Mol. Sci. 1, 337 (2011)

  46. [54]

    Tarus, K

    J. Tarus, K. Nordlund, J. Sillanp¨ a¨ a, and J. Keinonen, Heat spike and ballistic contributions to mixing in Si, Nucl. Instr. Meth. Phys. Res. B 153, 378 (1999)

  47. [55]

    Nordlund, N

    K. Nordlund, N. Runeberg, and D. Sundholm, Repul- sive interatomic potentials calculated using Hartree- Fock and density-functional theory methods, Nucl. In- str. Meth. Phys. Res. B 132, 45 (1997)

  48. [56]

    K. O. E. Henriksson, K. Nordlund, A. Krasheninnikov, and J. Keinonen, Differences in hydrogen and helium cluster formation, Appl. Phys. Lett. 87, 163113 (2005)

  49. [57]

    Kotakoski, A

    J. Kotakoski, A. V. Krasheninnikov, and K. Nordlund, A molecular dynamics study of the clustering of ion ir- radiation induced potassium in multiwalled carbon nan- otubes, Nucl. Instr. Meth. Phys. Res. B240, 810 (2005)

  50. [58]

    S. F. Boys and F. Bernardi, The calculation of small molecular interactions by the differences of separate to- tal energies. Some procedures with reduced errors, Mol. Phys. 19, 553 (1970)

  51. [59]

    Jensen, Polarization consistent basis sets: Principles, J

    F. Jensen, Polarization consistent basis sets: Principles, J. Chem. Phys. 115, 9113 (2001)

  52. [60]

    Lehtola, Curing basis set overcompleteness with piv- oted Cholesky decompositions, J

    S. Lehtola, Curing basis set overcompleteness with piv- oted Cholesky decompositions, J. Chem. Phys. 151, 241102 (2019), arXiv:1911.10372

  53. [61]

    Lehtola, A review on non-relativistic, fully numerical electronic structure calculations on atoms and diatomic molecules, Int

    S. Lehtola, A review on non-relativistic, fully numerical electronic structure calculations on atoms and diatomic molecules, Int. J. Quantum Chem. 119, e25968 (2019), arXiv:1902.01431

  54. [62]

    Lehtola, Fully numerical Hartree–Fock and density functional calculations

    S. Lehtola, Fully numerical Hartree–Fock and density functional calculations. II. Diatomic molecules, Int. J. Quantum Chem. 119, e25944 (2019), arXiv:1810.11653

  55. [63]

    Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cui, J. J. Eriksen, Y. Gao, S. Guo, J. Hermann, M. R. Hermes, K. Koh, P. Koval, S. Lehtola, Z. Li, J. Liu, N. Mardirossian, J. D. McClain, M. Motta, B. Mussard, H. Q. Pha...

  56. [64]

    Seeger and J

    R. Seeger and J. A. Pople, Self-consistent molecular orbital methods. XVIII. Constraints and stability in Hartree–Fock theory, J. Chem. Phys. 66, 3045 (1977)

  57. [65]

    Møller and M

    C. Møller and M. S. M. Plesset, Note on an approxima- tion treatment for many-electron systems, Phys. Rev. 46, 618 (1934)

  58. [66]

    S. H. Vosko, L. Wilk, and M. Nusair, Accurate spin- dependent electron liquid correlation energies for local spin density calculations: a critical analysis, Can. J. Phys. 58, 1200 (1980)

  59. [67]

    Peltola, K

    J. Peltola, K. Nordlund, and J. Keinonen, Molecular dynamics study on stopping powers of channeled He and Li ions in Si, Nucl. Instr. Meth. Phys. Res. B 217, 25 (2003)

  60. [68]

    Hobler, The role of lattice vibrations and inter- atomic potentials in the simulation of ion implantation, in NASECODE IX , edited by J

    G. Hobler, The role of lattice vibrations and inter- atomic potentials in the simulation of ion implantation, in NASECODE IX , edited by J. J. H. Miller (Front Range Press, Boulder, Colorado, 1993) pp. 55–56

  61. [69]

    Hobler, Monte Carlo simulation of two-dimensional implanted dopant distributions at mask edges, Nucl

    G. Hobler, Monte Carlo simulation of two-dimensional implanted dopant distributions at mask edges, Nucl. In- strum. Meth. Phys. Res. B 96, 155 (1995)

  62. [70]

    Juslin and K

    N. Juslin and K. Nordlund, Pair potential for Fe-He, J. Nucl. Mater. 382, 143 (2008)

  63. [71]

    Nordlund, F

    K. Nordlund, F. Djurabekova, and G. Hobler, Large fraction of crystal directions leads to ion channeling, Phys. Rev. B 94, 214109 (2016)

  64. [72]

    Delley, J

    J. Delley, J. Chem. Phys. 92, 508 (1990)

  65. [73]

    DMol is a trademark of AccelRys. Inc

  66. [74]

    2.3.5 ed

    DMol User Guide, Biosym Technologies Inc., San Diego, California, v. 2.3.5 ed. (1993)

  67. [75]

    Hohenberg and W

    P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Phys. Rev. 136, B864 (1964)

  68. [76]

    Kohn and L

    W. Kohn and L. J. Sham, Self-consistent equations in- cluding exchange and correlation effects, Phys. Rev. 26 140, A1133 (1965)

  69. [77]

    Keinonen, A

    J. Keinonen, A. Kuronen, K. Nordlund, R. M. Niemi- nen, and A. P. Seitsonen, First-Principles Simulation of Collision Cascades in Si to Test Pair-Potential for Si-Si Interaction at 10 eV – 5 keV, Nucl. Instr. Meth. Phys. Res. B 88, 382 (1994)

  70. [78]

    J. F. Ziegler, (1995), TRIM-95 computer code, private communication

  71. [79]

    C. Park, K. M. Klein, A. F. Tasch, and J. F. Ziegler, Critical angles for channeling of boron ions implanted into single-crystal silicon, J. Electrochem. Soc. 138, 2107 (1991)

  72. [80]

    Hobler and G

    G. Hobler and G. Betz, On the useful range of applica- tion of molecular dynamics simulations in the recoil in- teraction approximation, Nucl. Instr. Meth. Phys. Res. B 180, 203 (2001)

  73. [81]

    Nordlund, Molecular dynamics simulation of ion ranges in the 1 – 100 keV energy range, Comput

    K. Nordlund, Molecular dynamics simulation of ion ranges in the 1 – 100 keV energy range, Comput. Mater. Sci. 3, 448 (1995)

  74. [82]

    P. T. Wedepohl, Influence of electron distribution on atomic interaction potentials, Proc. Phys. Soc. 92, 10.1088/0370-1328/92/1/313 (1967)

  75. [83]

    Teller, On the stability of molecules in the Thomas– Fermi theory, Rev

    E. Teller, On the stability of molecules in the Thomas– Fermi theory, Rev. Mod. Phys. 34, 627 (1962)

  76. [84]

    W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C; The Art of Scientific Computing , 2nd ed. (Cambridge University Press, New York, 1995)

  77. [85]

    Eriksson and J

    L. Eriksson and J. A. Davies and P. Jespersgaard, Range Measurements in Oriented Tungsten Single Crys- tals (0.1-1.0 MeV). I. Electronic and Nuclear Stopping Powers, Phys. Rev. 161, 219 (1967)

  78. [86]

    Torri, J

    P. Torri, J. Keinonen, and K. Nordlund, A low-level detection system for hydrogen analysis with the reaction 1H(15N, αγ)12C, Nucl. Instr. Meth. Phys. Res. B84, 105 (1994)

  79. [87]

    Haussalo, K

    P. Haussalo, K. Nordlund, and J. Keinonen, Stopping of 5 – 100 keV helium in tantalum, niobium, tungsten, and AISI 316L steel, Nucl. Instr. Meth. Phys. Res. B 111, 1 (1996)

  80. [88]

    D. Cai, N. Grønbech-Jensen, C. M. Snell, and K. M. Beardmore, Phenomenological electronic stopping- power model for molecular dynamics and Monte Carlo simulation of ion implantation into silicon, Phys. Rev. B 54, 17147 (1996)

  81. [89]

    D. J. Larson, A. K. Petford-Long, A. Cerezo, G. D. W. Smith, D. T. Foord, and T. C. Anthony, Three- dimensional atom probe field-ion microscopy observa- tion of Cu/Co multilayer film structures, Appl. Phys. Lett. 73, 1125 (1998)

  82. [90]

    Buschorn, E

    G. Buschorn, E. Diedrich, W. Kufner, H. G. M.Rzepka, P. Hoffmann-Stascheck, and A. Richter, , Phys. Rev. B 55, 6196 (1997)

  83. [91]

    Nordlund and R

    K. Nordlund and R. S. Averback, Handbook of Ma- terials Modeling (Kluwer Academic, Dordrecht, The Netherlands, 2005) Chap. 6.2. Point defects in metals

  84. [92]

    See Supplemental Material at [URL will be inserted by publisher] for the ZBL pair-specific, DMol and MP2 data sets and the NLH potential coefficients, which are available in the compressed archive package file nlh potentials opendata.tar.gz

  85. [93]

    G. C. Abell, Empirical chemical pseudopotential theory of molecular and metallic bonding, Phys. Rev. B 31, 6184 (1985)

  86. [94]

    A presentation of the MDRANGE computer code is available at http://beam.mv.helsinki.fi/∼knordlun/mdh/mdh program.html

    Open source code available at https://gitlab.com/acclab/MDRANGE. A presentation of the MDRANGE computer code is available at http://beam.mv.helsinki.fi/∼knordlun/mdh/mdh program.html

  87. [95]

    Sillanp¨ a¨ a, J

    J. Sillanp¨ a¨ a, J. Peltola, K. Nordlund, J. Keinonen, and M. J. Puska, Electronic stopping calculated using explicit phase shift factors, Phys. Rev. B 63, 134113 (2001)

  88. [96]

    S. M. Hogg, B. Pipeleers, A. Vantomme, and M. Swartz, Channeling of low energy heavy ions: Er in Si ¡111¿, Appl. Phys. Lett. 80, 4363 (2016)

  89. [97]

    Peltola, K

    J. Peltola, K. Nordlund, and J. Keinonen, Electronic stopping power calculation method for molecular dy- namics simulations using local Firsov and free electron- gas models, Rad. Eff. & Def. in Sol. 161, 511 (2006)

  90. [98]

    D. S. Gemmell, Channeling and related effects in the motion of charged particles through crystals, Rev. Mod. Phys. 46, 129 (1974)

  91. [99]

    Nordlund and G

    K. Nordlund and G. Hobler, Dependence of ion chan- neling on relative atomic number in compounds, Nucl. Instr. Meth. Phys. Res. B 435, 61 (2018)

  92. [100]

    direction (fig. 20b). The data labeled ”MD” are simulated with the MD-RIA approach with the MDRANGE code, and ”SRIM” with the commonly used BCA code SRIM. Experimental data [99] shown for comparison. that are not compared to experiments (subsection VI A) have much stronger cha...

  93. [101]

    I. A. Stegun, Handbook of Mathematical Functions (Na- tional Bureau of Standards, Washington, DC, USA,

  94. [102]

    N. W. Ashcroft and N. D. Mermin, Solid State Physics (Saunders College, Philadelphia, 1976)

  95. [103]

    P. R. Bevington, Data reduction and error analysis for the physical sciences (McGraw-Hill, New York, 1992)

  96. [104]

    M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, Oxford, England, 1989)

  97. [105]

    Nordlund, J

    K. Nordlund, J. Keinonen, and T. Mattila, Formation of ion irradiation-induced small-scale defects on graphite surfaces, Phys. Rev. Lett. 77, 699 (1996)

  98. [106]

    K. Albe, K. Nordlund, and R. S. Averback, Model- ing metal-semiconductor interaction: Analytical bond- order potential for platinum-carbon, Phys. Rev. B 65, 195124 (2002)

  99. [107]

    Tersoff, New Empirical approach for the structure and energy of covalent systems, Phys

    J. Tersoff, New Empirical approach for the structure and energy of covalent systems, Phys. Rev. B 37, 6991 (1988)

  100. [108]

    D. W. Brenner, The art and science of an analytical potential, physica status solidi (b) 217, 23 (2000)

  101. [109]

    Behler, Perspective: Machine learning potentials for atomistic simulations, J

    J. Behler, Perspective: Machine learning potentials for atomistic simulations, J. Chem. Phys. 145, 170901 (2016)

  102. [110]

    986 ± 1 987 ± 1 986 ± 1 982 ± 1 1001 ± 2 Ar Al 30 Non-channeled 295 ± 1 295 ± 1 294 ± 1 291 ± 1 311 ± 1

  103. [111]

    Non-channeled means implantation in the non-channeling direction θ = ϕ = 20◦

    – 2180 ± 4 2180 ± 2 1997 ± 4 2072 ± 4 TABLE IV: Mean ranges ¯R in ˚A for the universal ZBL potential, the pair-specific ZBL potential, the DMol potential, and the MP2 potential for the studied test cases. Non-channeled means implantation in the non-channeling direction θ = ϕ =...

  104. [112]

    Byggm¨ astar, A

    J. Byggm¨ astar, A. Hamedani, K. Nordlund, and F. Djurabekova, Machine-learning interatomic poten- tial for radiation damage and defects in tungsten, Phys. Rev. B 100, 144105 (2019)

  105. [113]

    D. Cai, C. M. Snell, K. M. Beardmore, and N. Grønbech-Jensen, Simulation of phosphorus implan- tation into silicon with a single parameter electronic stopping power model, International J. Modern Physics C 9, 459 (1998)

  106. [114]

    https://gitlab.com/pn2200/ g3data

    G3data—grab graph data, a program for extract- ing data from graphs. https://gitlab.com/pn2200/ g3data. 27

  107. [115]

    J. F. Ziegler, SRIM-2013 software package, available on- line at http://www.srim.org

  108. [116]

    J. F. Ziegler, J. P. Biersack, and M. D. Ziegler, SRIM - The Stopping and Range of Ions in Matter (SRIM Co., Chester, Maryland, USA, 2008)

  109. [117]

    Sillanp¨ a¨ a, K

    J. Sillanp¨ a¨ a, K. Nordlund, and J. Keinonen, Electronic stopping of Silicon from a 3D Charge Distribution, Phys. Rev. B 62, 3109 (2000)

  110. [118]

    Juslin, P

    N. Juslin, P. Erhart, P. Tr¨ askelin, J. Nord, K. O. E. Hen- riksson, K. Nordlund, E. Salonen, and K. Albe, Analyt- ical interatomic potential for modelling non-equilibrium processes in the W-C-H system, J. Appl. Phys. 98, 123520 (2005)

  111. [119]

    Ghaly and R

    M. Ghaly and R. S. Averback, Effect of Viscous Flow on Ion damage near Solid Surfaces, Phys. Rev. Lett. 72, 364 (1994)

  112. [120]

    Nordlund, L

    K. Nordlund, L. Wei, Y. Zhong, and R. S. Averback, Role of electron-phonon coupling on collision cascade development in Ni, Pd and Pt, Phys. Rev. B (Rapid Comm.) 57, 13965 (1998)

  113. [121]

    Bj¨ orkas and K

    C. Bj¨ orkas and K. Nordlund, Comparative study of cas- cade damage in Fe simulated with recent potentials, Nucl. Instr. Meth. Phys. Res. B 259, 853 (2007), notes that ’The cutoff radius in the cluster analysis was sec- ond nearest neighbor and third nearest neighbor for the v...

  114. [122]

    A. E. Sand, J. Dequeker, C. S. Becquart, C. Domain, and K. Nordlund, Non-equilibrium properties of inter- atomic potentials in cascade simulations in tungsten, J. Nucl. Mater. 470, 119 (2016)

  115. [123]

    Belko, M

    V. Belko, M. Posselt, and E. Chagarov, Improvement of the repulsive part of the classical interatomic potential for SiC, Nucl. Instr. Meth. Phys. Res. B 202, 18 (2002)

  116. [124]

    https://www.scipy.org/scipylib/

  117. [125]

    Gombas, Die Statistische Theorie des Atoms und ihre Anwendungen (Springer-Verlag, Vienna, Austria, 1949)

    P. Gombas, Die Statistische Theorie des Atoms und ihre Anwendungen (Springer-Verlag, Vienna, Austria, 1949)

  118. [126]

    Eckstein, Computer Simulations of Ion-Solid Inter- actions (Springer, Berlin, 1991) p

    W. Eckstein, Computer Simulations of Ion-Solid Inter- actions (Springer, Berlin, 1991) p. 40

Pith tools