| 查看: 3314 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
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 |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有7人回复
-
申请2026年博士
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
关于pwscf的ibrav=0和不为0时的问题
已经有12人回复- 新手求助pwscf的进行磁性运算和LDA+U运算
已经有35人回复
- 编译pwscf提示找不到F77的库
已经有10人回复
- Pwscf编译出错
已经有4人回复
- pwscf是如何计算debye温度的
已经有7人回复
- pwscf并行计算时总是隔一段时间重新开始
已经有5人回复
- 如何让Pwscf找出全部对称性
已经有9人回复
- 请问用pwscf可以计算光学性质吗？
已经有8人回复
- 求PWSCF使用说明书
已经有3人回复
- PWscf计算homo lumo或者valence conduction band，问题已解决
已经有14人回复
- PWscf中CP的温度单位
已经有5人回复
- 【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
已经有14人回复
- 【求助】求助pwscf一些参数的物理意义和设定经验
已经有13人回复
- 【原创】PWSCF中的nbnd和degauss
已经有11人回复
- 【其他】pwscf中的原胞体积计算！
已经有4人回复
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1551
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
6楼2014-04-19 23:54:09
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1551
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
2楼2014-04-18 21:03:32
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2448.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
谢谢，确实是这个问题！ 但是又遇到一个新问题：计算的每次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
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2448.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
我后来测试发现没有定义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