Boundary value approaches to molecular dynamics simulation