An improved QM/MM approach for metals

We present an improved quantum mechanical (QM) and molecular mechanical (MM) coupling method for the study of metallic systems. The improved method is based on the earlier work of Choly et al (2005 Phys. Rev. B 71 094101). In this approach, quantum mechanical treatment is spatially confined to a small region, surrounded by a larger molecular mechanical region. This approach is particularly useful for systems where quantum mechanical interactions in a small region, such as lattice defects or chemical impurities, can affect the macroscopic properties of a material. We discuss how the coupling across the different scales can be accomplished efficiently and accurately for metals. The method is tested by performing a multiscale simulation of bulk aluminium (Al) where the coupling errors can be easily analysed. We then apply the method to study the core structure and Peierls stress of an edge dislocation in Al.

[1]  M. Baskes,et al.  Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals , 1984 .

[2]  G. Robert,et al.  Computation of the Peierls stress in tantalum with an extended-range modified embedded atom method potential , 2005 .

[3]  W. Benoit,et al.  The vacancies, lubrication agents of dislocation motion in aluminium , 1989 .

[4]  Y. H. Zhao,et al.  Compact and dissociated dislocations in aluminum: implications for deformation. , 2005, Physical review letters.

[5]  Robert E. Rudd,et al.  COARSE-GRAINED MOLECULAR DYNAMICS AND THE ATOMIC LIMIT OF FINITE ELEMENTS , 1998 .

[6]  Noam Bernstein,et al.  Mixed finite element and atomistic formulation for complex crystals , 1999 .

[7]  A. Nakano,et al.  Coupling length scales for multiscale atomistics-continuum simulations: atomistically induced stress distributions in Si/Si3N4 nanopixels. , 2001, Physical review letters.

[8]  M. Ortiz,et al.  Quasicontinuum analysis of defects in solids , 1996 .

[9]  E Weinan,et al.  Multiscale simulations in simple metals: A density-functional-based methodology , 2004, cond-mat/0404414.

[10]  W. Benoit,et al.  Internal friction and microdeformation due to the intrinsic properties of dislocations: The Bordoni relaxation , 1982 .

[11]  C. Woodward,et al.  Flexible Ab initio boundary conditions: simulating isolated dislocations in bcc Mo and Ta. , 2002, Physical review letters.

[12]  G. Gremaud,et al.  Lubrication agents of dislocation motion at very low temperature in cold-worked aluminium , 1990 .

[13]  William A. Curtin,et al.  A coupled atomistic/continuum model of defects in solids , 2002 .

[14]  M. Ortiz,et al.  An adaptive finite element approach to atomic-scale mechanics—the quasicontinuum method , 1997, cond-mat/9710027.

[15]  Efthimios Kaxiras,et al.  From Electrons to Finite Elements: A Concurrent Multiscale Approach for Metals , 2005, cond-mat/0506006.

[16]  David L. Olmsted,et al.  Lattice resistance and Peierls stress in finite size atomistic dislocation simulations , 2000, cond-mat/0010503.

[17]  F. Nabarro.  Theoretical and experimental estimates of the Peierls stress , 1997 .

[18]  Jiali Gao,et al.  Methods and Applications of Combined Quantum Mechanical and Molecular Mechanical Potentials , 2007 .

[19]  E. Kaxiras,et al.  Generalized-stacking-fault energy surface and dislocation properties of aluminum , 1999, cond-mat/9903440.

[20]  N. Govind,et al.  Electronic-structure calculations by first-principles density-based embedding of explicitly correlated systems , 1999 .

[21]  E. Kaxiras,et al.  Polarization switching in PbTiO3: an ab initio finite element simulation , 2002 .

[22]  Jacobsen,et al.  Interatomic interactions in the effective-medium theory. , 1987, Physical review. B, Condensed matter.

[23]  C. Kölmel,et al.  Combining ab initio techniques with analytical potential functions for structure predictions of large systems: Method and application to crystalline silica polymorphs , 1997 .

[24]  E Weinan,et al.  The Heterognous Multiscale Methods , 2003 .

[25]  Ru-wu Wang,et al.  Atomistic simulation of the atomic structure and diffusion within the core region of an edge dislocation in aluminum , 2000 .

[26]  W. E,et al.  Matching conditions in atomistic-continuum modeling of materials. , 2001, Physical review letters.

[27]  J. Q. Broughton,et al.  Concurrent coupling of length scales: Methodology and application , 1999 .

[28]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[29]  Michael J. Mills,et al.  A study of the structure of Lomer and 60° dislocations in aluminium using high-resolution transmission electron microscopy , 1989 .

[30]  Frederick H. Streitz,et al.  Quantum-based atomistic simulation of materials properties in transition metals , 2002 .

[31]  Emily A. Carter,et al.  Prediction of dislocation nucleation during nanoindentation of Al3Mg by the orbital-free density functional theory local quasicontinuum method , 2006 .

[32]  Rajiv K. Kalia,et al.  Coupling Length Scales for Multiscale Atomistics-Continuum Simulations , 2001 .

[33]  G. Gremaud,et al.  Study of the dislocation mechanism responsible for the Bordoni relaxation in aluminum by the two-wave acoustic coupling method , 1987 .