<div>Dear all,</div><div><br></div><div>To save CPU time, we fixed atoms of the surface slab in the adsorption model during the PHONON calculation.</div><div><br></div><div>the input file are as follows,</div><div><br></div><div><div>WMC-1.0</div><div>&INPUTPH</div><div>   outdir='/fs08/home/js_wjliu/espresso-6.1/qe-6.1/tempdir/WMC-1.0-half-gamma/' ,</div><div>   prefix='WMC-1.0-half-gamma',</div><div>   fildyn='WMC-1.0-half-gamma.dynG',</div><div>   nat_todo=3</div><div>      nogg = .true.,</div><div>alpha_mix(1) = 0.1,</div><div>/</div><div>0.0 0.0 0.0</div><div>49 50 51</div></div><div><br></div><div>So we think only the frequency 145-153 corresponds to atom 49-51 would be reported.</div><div><br></div><div>However, we find the freq 1-9, which are 9 IMAGINARY frequencies  also reported by PHONON. </div><div><br></div><div>Here are some information extracted from output file.</div><div><div><br></div><div><br></div><div>Mode symmetry, C_1 (1)     point group:</div><div><br></div><div><br></div><div>     Atomic displacements:</div><div>     There are 153 irreducible representations</div><div><br></div><div>     Representation     1      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     2      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     3      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     4      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     5      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     6      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     7      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     8      1 modes -A  Not done in this run</div><div><br></div><div>     Representation     9      1 modes -A  Not done in this run</div><div><br></div><div>     Representation    10      1 modes -A  Not done in this run</div></div><div><br></div><div><div>     Representation   144      1 modes -A  Not done in this run</div><div><br></div><div>     Representation   145      1 modes -A  To be done</div><div><br></div><div>     Representation   146      1 modes -A  To be done</div><div><br></div><div>     Representation   147      1 modes -A  To be done</div><div><br></div><div>     Representation   148      1 modes -A  To be done</div><div><br></div><div>     Representation   149      1 modes -A  To be done</div><div><br></div><div>     Representation   150      1 modes -A  To be done</div><div><br></div><div>     Representation   151      1 modes -A  To be done</div><div><br></div><div>     Representation   152      1 modes -A  To be done</div><div><br></div><div>     Representation   153      1 modes -A  To be done</div><div><br></div><div>     Compute atoms:    49,   50,   51,</div></div><div><br></div><div><div>     Alpha used in Ewald sum =   2.8000</div><div><br></div><div>     negative rho (up, down):  9.406E-05 3.627E-03</div><div>     PHONON       :     1h22m CPU       22h 1m WALL</div></div><div><br></div><div><br></div><div><br></div><div><div>     Mode symmetry, C_1 (1)     point group:</div><div><br></div><div>     freq (  1 -  1) =      -1086.5  [cm-1]   --> A               I+R</div><div>     freq (  2 -  2) =       -815.1  [cm-1]   --> A               I+R</div><div>     freq (  3 -  3) =       -764.6  [cm-1]   --> A               I+R</div><div>     freq (  4 -  4) =       -458.8  [cm-1]   --> A               I+R</div><div>     freq (  5 -  5) =       -240.8  [cm-1]   --> A               I+R</div><div>     freq (  6 -  6) =       -206.1  [cm-1]   --> A               I+R</div><div>     freq (  7 -  7) =       -127.7  [cm-1]   --> A               I+R</div><div>     freq (  8 -  8) =        -59.4  [cm-1]   --> A               I+R</div><div>     freq (  9 -  9) =        -21.3  [cm-1]   --> A               I+R</div><div>     freq (145 -145) =        255.2  [cm-1]   --> A               I+R</div><div>     freq (146 -146) =        438.0  [cm-1]   --> A               I+R</div><div>     freq (147 -147) =        507.2  [cm-1]   --> A               I+R</div><div>     freq (148 -148) =        681.1  [cm-1]   --> A               I+R</div><div>     freq (149 -149) =        810.3  [cm-1]   --> A               I+R</div><div>     freq (150 -150) =       1124.5  [cm-1]   --> A               I+R</div><div>     freq (151 -151) =       1568.3  [cm-1]   --> A               I+R</div><div>     freq (152 -152) =       3534.8  [cm-1]   --> A               I+R</div><div>     freq (153 -153) =       3559.3  [cm-1]   --> A               I+R</div></div><div><br></div><div>We also find although we input a unstable adsorption configuration for atom 49-51, there is no imaginary frequency reported in freq 145-153.</div><div><br></div><div>Could anyone help us to deal with this problem. Thanks in advance.</div><div><br></div><div>best wishes</div><div><br></div><div>vega</div><div><div style="color:#909090;font-family:Arial Narrow;font-size:12px">------------------</div><div style="font-size:14px;font-family:Verdana;color:#000;"><div><div style="   ;; color: rgb(144, 144, 144); font-size: 12px;  ">------------------</div><div><div style="line-height: 23px;"><div style="line-height: 23px;">vega</div><div style="line-height: 23px;">Department of physics, Nanjing University, China</div></div></div></div></div></div><div> </div>