A numerical study of residual stress induced in machined silicon surfaces by molecular dynamics simulation