230 likes | 389 Views
Bond-Order Potential for MD Simulation: Relaxation of Semiconductor Nanostructures. tight binding and bond order 4th moment approximation parameterization and fit some examples. Volker Kuhlmann and Kurt Scheerschmidt Max Planck-Institute of Microstructure Physics Halle - Germany.
E N D
Bond-Order Potential for MD Simulation:Relaxation of Semiconductor Nanostructures • tight binding and bond order • 4th moment approximation • parameterization and fit • some examples Volker Kuhlmann and Kurt Scheerschmidt Max Planck-Institute of Microstructure Physics Halle - Germany
large time and length scales accurate atomistic potential quantum mechanics of electrons (slow) empirical potential (fast) pair potential many-body cluster expansion bond order potential density functional theory • - transferable • few parameter • chemical bonds tight binding
Tight Binding exact diagonalisation two-center approximation: Slater-Koster integrals: electronic part (bandstructure) scaling part (elastic constants)
Bond Order Potential Greens function: many atom expansion local density of states moment
2nd moment: contribution negligible normalized moment: angular function: reduced TB parameter:
new contributions to bond terms : torsion angle: on site term :
contribution of largest at constant angle
of most pronounced new angular dependence at constant angle
Potential energy above Si(100) surface BOP2 BOP4 BOP4+ minimum minimum raised maximum
Parametrization and Fit 7 parameter
smooth promotion energy invested energy: promote one electron Gained energy: form new bonds
propose and accept/reject fit via Monte Carlo/ Conjugate gradient fitness of set {r}:
improved 4th moments and promotion energyfor pure carbon systems
simulation of Si(100) waferbonding with rotational twist Scheerschmidt and Kuhlmann, Interface Science 12 (2004)
recursion method and local density of states • solve Gii recursively: • LDOS approximated by moments: moments-theorem • semi-infinite linear chain: ai=a=0 eV bi=b=0.1 eV
adjust parameter to recover properties (Ro,Ucoh,B,C11,…) • s(r) must die out suffic. before cut off via spline • must cut off before 2nd nearest neighbors: • # of paths of length 4 (4th moment) = Nbrs^2 • 256 paths @ 16Nbrs vs. 16 paths @ 4Nbrs • 6th Moment : 64 vs. 4096 • low slopes (n,m) required by elasticity conflict with cutoff -> make a compromise