Title: Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Journal Article · · Materials Research Express (Online)
ORCiD logo; ORCiD logo

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Research Organization:
Univ. of Illinois at Urbana-Champaign, IL (United States)
Sponsoring Organization:
USDOE; USDOE Office of Science (SC); USDOE Office of Science (SC), Fusion Energy Sciences (FES)
Grant/Contract Number:
SC0018141
OSTI ID:
2251547
Journal Information:
Materials Research Express (Online), Journal Name: Materials Research Express (Online) Journal Issue: 12 Vol. 10; ISSN 2053-1591
Publisher:
IOP PublishingCopyright Statement
Country of Publication:
United Kingdom
Language:
English

References (27)

Chebfun: A New Kind of Numerical Computing book January 2010
Tridyn-binary collision simulation of atomic collisions and dynamic composition changes in solids journal November 1988
The optimization of convergence for chebyshev polynomial methods in an unbounded domain journal January 1982
A modification of the delves-lyness method for locating the zeros of analytic functions journal July 1985
A Monte Carlo computer program for the transport of energetic ions in amorphous targets journal August 1980
Tridyn — A TRIM simulation code including dynamic composition changes journal March 1984
The temporal development of collision cascades in the binary-collision approximation journal March 1990
Computing real roots of a polynomial in Chebyshev series form through subdivision journal August 2006
Exponentially-convergent strategies for defeating the Runge Phenomenon for the approximation of non-periodic functions, part two: Multi-interval polynomial schemes and multidomain Chebyshev interpolation journal April 2011
Numerical experiments on the accuracy of the Chebyshev–Frobenius companion matrix method for finding the zeros of a truncated series of Chebyshev polynomials journal August 2007
F-TRIDYN simulations of tungsten self-sputtering and applications to coupling plasma and material codes journal June 2018
General numerical algorithm for classical collision integral calculation journal June 2008
F-TRIDYN: A Binary Collision Approximation code for simulating ion interactions with rough surfaces journal October 2017
Modelling radiation effects using the ab-initio based tungsten and vanadium potentials
  • Björkas, C.; Nordlund, K.; Dudarev, S.
  • Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, Vol. 267, Issue 18 https://doi.org/10.1016/j.nimb.2009.06.123
journal September 2009
SRIM – The stopping and range of ions in matter (2010)
  • Ziegler, James F.; Ziegler, M. D.; Biersack, J. P.
  • Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, Vol. 268, Issue 11-12 https://doi.org/10.1016/j.nimb.2010.02.091
journal June 2010
Assessment of surface potential models by molecular dynamics simulations of atom ejection from (100)-Si surfaces journal May 2013
The Lennard-Jones potential: when (not) to use it journal January 2020
Attractive interaction potentials and the binary collision approximation journal June 1997
Interatomic potentials for modelling radiation defects and dislocations in tungsten journal September 2013
Embedded-atom method potential for modeling hydrogen and hydrogen-defect interaction in tungsten journal September 2017
A numerical method for locating the zeros of an analytic function journal January 1967
Polynomial roots from companion matrix eigenvalues journal May 1995
Multiscale modeling of crowdion and vacancy defects in body-centered-cubic transition metals journal August 2007
Computer simulation of atomic-displacement cascades in solids in the binary-collision approximation journal June 1974
Solving Transcendental Equations book January 2014
Finding the Zeros of a Univariate Equation: Proxy Rootfinders, Chebyshev Interpolation, and the Companion Matrix journal January 2013
Transport Properties of the Simple Lennard-Jones/Spline Fluid I: Binary Scattering and High-Accuracy Low-Density Transport Coefficients journal July 2020

Similar Records

Related Subjects

Materials Science