| 查看: 1169 | 回复: 3 | |||
WDD880227木虫 (正式写手)
|
[求助]
pwscf如何设置nelec 已有2人参与
|
|
我想用pwscf计算n型石墨烯，例如一个n型石墨烯原胞的电子数比本征石墨烯多0.01， 那么我设置nelec=8.01，然后进行relax计算，但是输出文件中错误信息“Error in routine read_namelists (19): reading namelist system” 另外，如果不设置nelec=8.01进行relax，能够顺利优化完成。 那么问题出在nelec的设置上，nelect该怎么设置呢？ 下面是我的in 文件： &control prefix = 'graphene' calculation = 'relax' restart_mode = 'from_scratch' tstress = .true. tprnfor = .true. outdir = './' wfcdir = './' pseudo_dir = './' nstep = 1000 etot_conv_thr = 3.0E-9 forc_conv_thr = 1.0D-5 / &system ibrav = 0 celldm(1) = 4.6487281193142 celldm(3) = 4.8780487804878 nat = 2 ntyp = 1 nbnd = 10 nelec = 8.02 ecutwfc = 60.0 / &electrons electron_maxstep = 100 conv_thr = 1.0d-8 / CELL_PARAMETERS 0.86602540378443864 -0.5 0.0 0.0 1.0 0.0 0.0 0.0 1.0 ATOMIC_SPECIES C 12.011 C.UPF ATOMIC_POSITIONS crystal C 0.0000000000000000 0.0000000000000000 0.5000000000000000 C 0.3333333333333333 0.6666666666666667 0.5000000000000000 K_POINTS automatic 36 36 1 1 1 1 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2014-07-31 20:44:22
3楼2014-08-01 17:59:19
4楼2014-08-02 01:11:12
