Mikiyas Tsegaye wrote: > I thought the solution would be to include epsil = .true. in the ph > input file, but that makes the calculation crash at that step. no, it makes the calculation "stop", not "crash", because you cannot calculate effective charges in a metal P. -- Paolo Giannozzi, Democritos and University of Udine, Italy