A generalized solid-state nudged elastic band method.

A generalized solid-state nudged elastic band (G-SSNEB) method is presented for determining reaction pathways of solid-solid transformations involving both atomic and unit-cell degrees of freedom. We combine atomic and cell degrees of freedom into a unified description of the crystal structure so that calculated reaction paths are insensitive to the choice of periodic cell. For the rock-salt to wurtzite transition in CdSe, we demonstrate that the method is robust for mechanisms dominated either by atomic motion or by unit-cell deformation; notably, the lowest-energy transition mechanism found by our G-SSNEB changes with cell size from a concerted transformation of the cell coordinates in small cells to a nucleation event in large cells. The method is efficient and can be applied to systems in which the force and stress tensor are calculated using density functional theory.

[1]  Duane D. Johnson,et al.  bcc-to-hcp transformation pathways for iron versus hydrostatic pressure: Coupled shuffle and shear modes , 2009 .

[2]  H. Jónsson,et al.  Nudged elastic band method for finding minimum energy paths of transitions , 1998 .

[3]  E. Rabani An interatomic pair potential for cadmium selenide , 2002 .

[4]  New mechanism for the alpha to omega martensitic transformation in pure titanium. , 2003, Physical review letters.

[5]  Emily A Carter,et al.  Finding transition states for crystalline solid-solid phase transformations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[6]  G. Henkelman,et al.  Theoretical calculations of CH4 and H2 associative desorption from Ni(111): could subsurface hydrogen play an important role? , 2006, The Journal of chemical physics.

[7]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[8]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[9]  H. Monkhorst,et al.  SPECIAL POINTS FOR BRILLOUIN-ZONE INTEGRATIONS , 1976 .

[10]  J. Polanyi,et al.  Location of Energy Barriers. I. Effect on the Dynamics of Reactions A + BC , 1969 .

[11]  E. Carter,et al.  Nonadiabaticity in the iron bcc to hcp phase transformation. , 2008, The Journal of chemical physics.

[12]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[13]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .