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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf的AIMD参数
51/1返回列表
查看: 3350  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

c_chi

木虫 (正式写手)

[求助] pwscf的AIMD参数 已有1人参与

我刚接触pwscf，需要用来做AIMD的计算，参数如下，但总是报错
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine  read_namelists (17):
      reading namelist control
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
我知道是&system中参数有问题，但仔细检查和调整仍然不行，有没有人火眼金睛帮我检查一下，谢谢了！！

&control
    calculation='md'
    restart_mode='from_scratch',
    pseudo_dir = ./,
    outdir= ./,
    prefix='Gra'
    tprnfor = .true.
    tstress = .true.
    etot_conv_thr = 1.0D-4    !Default: 1.0D-4  (a.u)
    forc_conv_thr = 1.0D-3    !Default: 1.0D-3  (a.u)
    dt=20,        !time step for molecular dynamics, in Rydberg atomic units;1 a.u.=4.8378 * 10^-17 s
    nstep = 100   !number of ionic + electronic steps
    verbosity = 'high'
/
&system
    ibrav = 0,
    nat = 2, ntyp=1,
    ecutwfc = 30.0,              ! kinetic energy cutoff (Ry) for wavefunctions
    ecutrho = 150.0,             ! kinetic energy cutoff (Ry) for charge density and potential;Default:4 * ecutwfc
    occupations='smearing',      ! 'tetrahedra' especially suited for calculation of DOS,Not suitable for force/optimization/dynamics calculations
    smearing='mp', degauss=0.03  ! 'fermi-dirac'
    nspin=1,                     !1:non-polarized calculation ;2:spin-polarized calculation, LSDA;4,spin-polarized calculation, noncollinear
    tot_charge=0                 !tot_charge=+1 means one electron missing from the system
/
&electrons
    electron_maxstep=100
    diagonalization = 'david'  ! cg
    mixing_mode     = 'plain'  ! TF
    mixing_beta     = 0.5
    conv_thr        = 1.D-6
/
&ions
    ion_dynamics = 'verlet'
    pot_extrapolation = 'second_order'
    wfc_extrapolation = 'second_order'
    ion_temperature = 'rescale-v'
    tempw = 293.0             !Starting temperature (Kelvin)
    nraise=1
    tolp=100                  !Tolerance for velocity rescaling
   #delta_t = 1
/
ATOMIC_SPECIES
C   12.0107   C.pbe-rrkjus.UPF
ATOMIC_POSITIONS crystal    !alat | bohr | angstrom | crystal ; 0, fixed
C    0.00000   0.00000   0.5   0 0 0
C    0.33333   0.66667   0.5   1 1 1
K_POINTS automatic
8 8 1 0 0 0
CELL_PARAMETERS   angstrom    ! alat | bohr | angstrom ; needed only if ibrav = 0 is specified, ignored otherwise
   2.13042   -1.23 0
   0.000    2.46  0
  0   0  18
回复此楼

» 猜你喜欢

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

c_chi

木虫 (正式写手)
引用回帖:
8楼: Originally posted by goldenfisher at 2014-04-21 19:59:45
以ang为单位肯定是可以的。如果你设置了celldm（1），那么就是以celldm（1）为单位，而celldm一般默认就是以bohr为单位。我就经常用ang单位。
celldm默认是bohr为单位，A，B，C默认是ang为单位，说明文件里有...

恩，谢谢你的耐心解释
赞一下 回复此楼
9楼2014-04-21 20:49:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
c_chi: 金币+5, ★★★很有帮助 2014-04-19 14:00:59
看提示应该是Control中的问题。试试pseudo和out的目录打上引号
赞一下 回复此楼
2楼2014-04-18 21:03:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)
引用回帖:
2楼: Originally posted by goldenfisher at 2014-04-18 21:03:32
看提示应该是Control中的问题。试试pseudo和out的目录打上引号

谢谢，确实是这个问题！
但是又遇到一个新问题：计算的每次MD的前两部步scf都会出现很多
     c_bands:  4 eigenvalues not converged
最后在4~5个MD后就出现下面的错误提示跳出计算。我尝试修改晶格参数，截断能都不能解决。我想问的是下面的错误是我的什么参数有问题才会使计算进行不下去吗？谢谢了。。

在论坛上也看到别人提到“c_bands:  4 eigenvalues not converged
”的警告可以忽略，只要在随后scf收敛就行。这个警告是和模型有关吗？是不是想其他人说的这个警告不能通过修改参数解决？

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine c_bands (1):
     too many bands are not converged
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
赞一下 回复此楼
3楼2014-04-19 14:00:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)
引用回帖:
2楼: Originally posted by goldenfisher at 2014-04-18 21:03:32
看提示应该是Control中的问题。试试pseudo和out的目录打上引号

我后来测试发现没有定义celldm(1)，于是我添加celldm(1)=1.8897269（bhor和angstrom的换算），上面的error没有了，但是有新的error：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine cell_base_init (1):
     lattice vectors in Bohr or in a0 units?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
我找到源程序cell_base.f90看了一下，但是没读懂，不知道这是因为什么？能不能指点一下，谢谢了！
      SELECT CASE ( TRIM( cell_units ) )
        CASE ( 'bohr' )
          IF( celldm( 1 ) /= 0.0_DP .OR. a /= 0.0_dp ) CALL errore &
              ('cell_base_init','lattice vectors in Bohr or in a0 units?',1)
          units = 1.0_DP
        CASE ( 'angstrom' )
          IF( celldm( 1 ) /= 0.0_DP .OR. a /= 0.0_dp ) CALL errore &
              ('cell_base_init','lattice vectors in A or in a0 units?',2)
          units = 1.0_DP / bohr_radius_angs
        CASE DEFAULT
          IF( celldm( 1 ) /= 0.0_DP ) THEN
             units = celldm( 1 )
          ELSE IF ( a /= 0.0_dp ) THEN
             units = a / bohr_radius_angs
          ELSE
             units = 1.0_DP
          END IF
     END SELECT
赞一下 回复此楼
4楼2014-04-19 20:35:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭