vega lew wrote: > I am not sure. But when I specify the 'ion_dynamics = 'bfgs' ,', the > relax calculation could run smoothly. of course it does, but does it enforce constraints? -- Paolo Giannozzi, Democritos and University of Udine, Italy