Reactive molecular dynamics simulation