<div dir="ltr"><div><div><div>Dear All,<br></div>    I am trying to calculate IR frequency of NH4Br crystal. In <a href="http://ph.in">ph.in</a> and <a href="http://dynmat.in">dynmat.in</a> I have used asr=.true., and asr='zero-dim' respectively. So, I should get 6 zero frequency due to asr option. But I am getting the following frequencies <br>
<br>-----------------------------------------------------<br># mode   [cm-1]    [THz]      IR<br>    1   -265.41   -7.9566    8.4124<br>    2   -251.33   -7.5346    8.3177<br>    3   -224.98   -6.7446    1.1717<br>    4   -108.27   -3.2458    0.0181<br>
    5     -0.00   -0.0000    0.0000<br>    6     -0.00   -0.0000    0.0000<br>    7     -0.00   -0.0000    0.0000<br>    8      0.00    0.0000    0.0000<br>    9      0.00    0.0000    0.0000<br>   10      0.00    0.0000    0.0000<br>
   11    192.15    5.7604    0.0005<br>   12    238.67    7.1553    0.4008<br>   13    264.00    7.9144    3.5942<br>   14    293.57    8.8012    0.6786<br>   15    329.56    9.8799    1.7381<br>   16    392.07   11.7539    0.0041<br>
   17    409.82   12.2862    0.0005<br>   18    413.05   12.3828    0.0013<br>   19   1088.47   32.6315    1.3513<br>   20   1285.27   38.5313    2.2625<br>   21   1393.04   41.7624    1.9338<br>   22   1433.41   42.9725    3.7008<br>
   23   1476.19   44.2550    1.3273<br>   24   1493.99   44.7888    1.1527<br>   25   1597.18   47.8821    0.6625<br>   26   1635.26   49.0238    0.0000<br>   27   1676.34   50.2554    0.0047<br>   28   1686.18   50.5505    0.0000<br>
   29   2657.92   79.6825   28.6572<br>   30   2931.08   87.8715   14.0421<br>   31   3056.17   91.6215    0.2991<br>   32   3078.70   92.2970   18.7792<br>   33   3094.15   92.7603    7.0421<br>   34   3148.71   94.3958   61.3242<br>
   35   3176.09   95.2167   69.5553<br>   36   3474.05  104.1493   36.0143<br>-----------------------------------------------------------------------------------------------<br><br></div>I am also giving my inputs also  <br>
</div>-----------------------------------<br><a href="http://scf.in">scf.in</a><br><div><div>&control<br>    calculation = 'scf',<br>    title = 'NH4Br'<br>    restart_mode = 'from_scratch'<br>
    tprnfor = .TRUE.<br>    prefix = 'NH4Br',<br>    outdir = '------'<br>    pseudo_dir = '---------'<br>    nstep=1000<br> /<br> &system<br>    ibrav = 0,<br>     nat = 12,<br>    ntyp = 3,<br>
    nbnd = 60<br>    tot_charge = 0<br>   celldm(1) = 1.890359168<br>    ecutwfc = 30.00<br>    ecutrho = 300.00<br>    london      = .true.<br>    london_s6   = 0.75<br>    london_rcut = 200.00<br> /<br> &electrons<br>
 electron_maxstep = 1000<br>  mixing_beta=0.3<br>  conv_thr = 1.0d-8<br> diagonalization='david'<br>  diago_david_ndim = 4<br>/<br> &CELL<br>   cell_dynamics='bfgs'<br>   press = 0.0<br> /<br>ATOMIC_SPECIES<br>
    H    1.00794   H.pbe-van_ak.UPF<br>    Br   79.904   Br.pbe-van_mit.UPF<br>    N    14.0067   N.pbe-van_ak.UPF<br>ATOMIC_POSITIONS angstrom<br>N        0.000000000   0.000000000   0.000000000<br>.<br>.<br>H        0.000000000   4.781723484   3.381166715<br>
K_POINTS automatic<br>6 6 6    0 0 0<br>CELL_PARAMETERS angstrom<br>   5.633079685   0.000000000   0.000000000<br>   0.000000000   5.633079685   0.000000000<br>   0.000000000   0.000000000   3.982678150<br><br><br></div><div>
<a href="http://ph.in">ph.in</a><br><br>&inputph<br>  tr2_ph=4.0d-17,<br>  prefix='---',<br>  outdir='------',<br>  amass(1)=1.d0,<br>  amass(2)=79.9d0,<br>  amass(3)=14.d0,<br>  asr=.true.,<br>  zue=.true.,<br>
  epsil=.true.,<br>  trans=.true.,<br>  fildyn='NH4Br.dyn.xml',<br> /<br>0.0 0.0 0.0<br><br></div><div><br></div><div><a href="http://dynmat.in">dynmat.in</a><br><br> &input<br>  fildyn='NH4Br.dyn.xml',<br>
  asr='zero-dim'<br> /<br><br>------------------------------------------------------------------------------------------<br></div><div>Why I am getting such kind of frequencies and how can I solve it? Any suggestions is highly appreciated. <br>
<br></div><div>Thanking you in advance.<br></div><div>Pratik<br></div><div>Kolkata, India<br></div><div><br></div></div></div>