Molecular dynamics study of the finite-size effect in 2D nanoribbon silicene

We present the molecular dynamics (MD) simulation of the finite-size effect in 2D silicene nanoribbons (SiNRs). Five silicene nanoribbon structures were modelled with the same length and buckling, but different width sizes. All models were melted from the perfect honeycomb silicene structure and then cooled from the disordered liquid state. Structural properties such as ring sizes, radial distribution function, coordination numbers, and interatomic distances have been carefully analysed. All models’ critical melting and cooling temperature were determined via total energy per atom. We found that 2D SiNRs transition temperature is lower than the pristine silicene. Particularly, the SiNRs present a noticeable finite-size effect on nanoribbons, resulting in a domain-type structure of 4, 5, and 6-fold rings.