应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 声子谱计算求助
41/1返回列表
查看: 1152  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fengshiquan

铁杆木虫 (正式写手)

[求助] 声子谱计算求助

想计算SiO2的声子谱，输入文件如下：
   ndtset   10
#Set 1 : ground state self-consistency

  getwfk1   0            # Cancel default
  kptopt1   1            # Automatic generation of k points, taking
                         # into account the symmetry
    nqpt1   0            # Cancel default
  tolvrs1   1.0d-20      # SCF stopping criterion (modify default)
  rfphon1   0            # Cancel default


#Q vectors for all datasets

#Complete set of symmetry-inequivalent qpt chosen to be commensurate
# with kpt mesh so that only one set of GS wave functions is needed.
#Generated automatically by running GS calculation with kptopt=1,
# nshift=0, shiftk=0 0 0 (to include gamma) and taking output kpt set
# file as qpt set. Set nstep=1 so only one iteration runs.

     nqpt   1            # One qpt for each dataset (only 0 or 1 allowed)
                         # This is the default for all datasets and must
                         #  be explicitly turned off for dataset 1.

     qpt2   0.00000000E+00  0.00000000E+00  0.00000000E+00
     qpt3   0.00000000E+00  0.00000000E+00  0.00000000E+00
     qpt4   2.50000000E-01  0.00000000E+00  0.00000000E+00
     qpt5   5.00000000E-01  0.00000000E+00  0.00000000E+00
     qpt6   2.50000000E-01  2.50000000E-01  0.00000000E+00
     qpt7   5.00000000E-01  2.50000000E-01  0.00000000E+00
     qpt8  -2.50000000E-01  2.50000000E-01  0.00000000E+00
     qpt9   5.00000000E-01  5.00000000E-01  0.00000000E+00
     qpt10 -2.50000000E-01  5.00000000E-01  2.50000000E-01

#Set 2 : Response function calculation of d/dk wave function

    iscf2   -3         # Need this non-self-consistent option for d/dk
  kptopt2   2          # Modify default to use time-reversal symmetry
  rfphon2   0          # Cancel default
  rfelfd2   2          # Calculate d/dk wave function only
  tolvrs2   0.0        # Cancel default for d/dk
  tolwfr2   1.0d-22    # Use wave function residual criterion instead

#Set 3 : Response function calculation of Q=0 phonons and electric field pert.

  getddk3   2          # d/dk wave functions from last dataset
  kptopt3   2          # Modify default to use time-reversal symmetry
  rfelfd3   3          # Electric-field perturbation response only

#Sets 4-10 : Finite-wave-vector phonon calculations (defaults for all datasets)

   getwfk   1          # Use GS wave functions from dataset1
   kptopt   3          # Need full k-point set for finite-Q response
   rfphon   1          # Do phonon response
  rfatpol   1 2        # Treat displacements of all atoms
    rfdir   1 1 1      # Do all directions (symmetry will be used)
   tolvrs   1.0d-9     # This default is active for sets 3-10

#######################################################################
#Common input variables

#Definition of the unit cell
    acell   3*7.1600 angstrom         # This is equivalent to   10.61 10.61 10.61
    rprim   0.0  0.5  0.5   # In lessons 1 and 2, these primitive vectors
            0.5  0.0  0.5   # (to be scaled by acell) were 1 0 0  0 1 0  0 0 1
            0.5  0.5  0.0   # that is, the default.

#Definition of the atom types
   ntypat   2         # There are two types of atom
    znucl   14 8     # The keyword "znucl" refers to the atomic number of the
                      # possible type(s) of atom. The pseudopotential(s)
                      # mentioned in the "files" file must correspond
                      # to the type(s) of atom. Here, type 1 is the Aluminum,
                      # type 2 is the Arsenic.

#Definition of the atoms
    natom   6         # There are two atoms
    typat   1 1 2 2 2 2       # The first is of type 1 (Ga), the second is of type 2 (As).
                       
     xred  0.000  0.000  0.000  # Triplet giving the REDUCED coordinate of atom 1.
           0.250  0.250  0.250
           0.125  0.125  0.125  # Triplet giving the REDUCED coordinate of atom 2.
           0.125  0.125  0.625
           0.125  0.625  0.125
           0.625  0.125  0.125  

#Gives the number of band, explicitely (do not take the default)
    nband   30  
      

#Exchange-correlation functional

      ixc   1             # LDA Teter Pade parametrization

#Definition of the planewave basis set

     ecut   25.0           # Maximal kinetic energy cut-off, in Hartree

#Definition of the k-point grid
    ngkpt   4  4  4         
  nshiftk   4              # Use one copy of grid only (default)
   shiftk   0.0 0.0 0.5    # This gives the usual fcc Monkhorst-Pack grid
            0.0 0.5 0.0
            0.5 0.0 0.0
            0.5 0.5 0.5

#Definition of the SCF procedure
     iscf   5          # Self-consistent calculation, using algorithm 5
    nstep   100         # Maximal number of SCF cycles
   diemac   12        # Although this is not mandatory, it is worth to
                       # precondition the SCF cycle. The model dielectric
                       # function used as the standard preconditioner
                       # is described in the "dielng" input variable section.
                       # The dielectric constant of GaAs is smaller that the one of Si (=12).
得到的结果中发下G点的声子频率
-  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
-  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
-  0.000000E+00  0.000000E+00  6.753131E+02  6.753131E+02  8.131370E+02
-  8.376450E+02  8.376450E+02  8.376450E+02
和我在castep中得到的结果差别很大，而且在结果文件中报了如下错误：
chkph3 : WARNING -
  The dynamical matrix was incomplete : phonon frequencies may be wrong ...
请高手指点，是什么地方设置的不当，该如何修正，不胜感激
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-06 14:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
自己顶一下，有没有高手知道问题出在哪儿，还请不吝赐教，谢谢了
赞一下 回复此楼
2楼2013-04-07 10:20:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pariah

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fengshiquan: 金币+10, ★★★★★最佳答案, 谢谢你的帮助，去改一下试试 2013-04-07 16:17:58
franch: 金币+2, 谢谢回帖交流，， 2013-04-08 21:14:05
rfatpol 参数改为 1 natom 对应你的计算也就是 1 6
赞一下 回复此楼
3楼2013-04-07 16:08:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

若博特

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fengshiquan: 金币+10, ★有帮助, 呵呵，楼上已经提到了，我去改了这个参数，貌似确实是这个问题，还是要谢谢你的帮助，随便问下，改了这个参数出现了较大虚频，需要修改那个收敛精度，上面的in文件中有几个控制收敛精度的参数，不知道调哪个可以消去须频，谢谢 2013-04-08 10:25:14
franch: 金币+2, 谢谢回帖交流，， 2013-04-08 21:13:58
rfatpol   1 2        # Treat displacements of all atoms
你的这个参数貌似有问题啊
这个参数让你对那几个原子进行移动你输入的是第一个原子和第二个原子：12
而后边却总共有6个原子这个不对称的。
赞一下 回复此楼
4楼2013-04-08 10:13:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fengshiquan 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭