<div>Dear All,<br><br>Steve and Jibiao thank you for your help.<br><br>I performed the calculations of two OH <br>groups on the surface using the nat_todo <br>setting implemented in ph.x. <br><br>However, there is a problem with the dynmat.x <br>program, namely, when I used asr='zero-dim',  <br>the program stucks on the ASR calculations.<br>....</div>
<div><br>Reading Dynamical Matrix from file vac_2OH.dyn<br>...Force constants read<br>...epsilon and Z* not read (not found on file)<br>Acoustic Sum Rule: || Z*(ASR) - Z*(orig)|| =    0.000000E+00<br><br>This is my input file for ph.x:<br><br>&inputph<br>  prefix = 'vac_2OH'<br>  nat_todo = 4<br>  amass(1) = 12.01<br>  amass(2) = 14.0067<br>  amass(3) = 1.0079<br>  amass(4) = 15.9994<br>  outdir = './'<br>  fildyn = 'vac_2OH.dyn'<br>  epsil = .false.<br>/<br>0.0 0.0 0.0<br>252, 253, 254, 255<br><br>and for dynmat.x<br><br>&input<br> fildyn = 'vac_2OH.dyn'<br> asr = 'zero-dim'<br>/</div>
<div> </div>
<div><br>Have you had a similar problem?</div>
<div>Do you know how to solve it?<br><br>Best,</div>
<div>Tom</div>