use c++11 threads to sum charge density
output format fix
start form saved potential