| 查看: 3774 | 回复: 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 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有82人回复
-
感觉是下周放榜了
已经有8人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
2楼2014-04-18 21:03:32
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.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 (小学生)
- 金币: 2449.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
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
' ) notconv5楼2014-04-19 20:43:18
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
6楼2014-04-19 23:54:09
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
我是仔细研究了一下说明文件的，里面关于cell parameter的参数描述只有alat | bohr | angstrom 。 我后来又测试了一下，发现是不用定义celdm，但必须cell parameter里面的晶格矢量必须是a.u.，而我这里是想以Ang为单位，于是设置参数为angstrom，但这样似乎不可行，输出的alat大小始终是a矢量的模，单位是a.u.；即使设置成CELL_PARAMETERS (alat=4.64873），但输出仍然是alat=2.46 a.u.。貌似加了参数Hexagonal 也不行。（另外我不知道出了说明文档以外的参数描述在哪里可以找到，这些参数我只是在其他人的例子里面看到过） 所以，解决办法有三个： 1. ibrav=0， 然后CELL_PARAMETERS 里面的矢量转换成a.u.单位 2. ibrav=4，a=2.46，c=14，就不用设置CELL_PARAMETERS 3.ibrav = 4, celldm(1)=4.64873, celldm(3)=6 ，也不用设置CELL_PARAMETERS 但是，这里有个问题，a和celldm(1)不应该是等价的吗？单位都是a.u.；但第二种里面的单位显然是Ang，这里有点迷糊。 还有就是，我如果想使用ibrav=0，并且使得CELL_PARAMETERS 里面的矢量以Ang为单位是不是不可行？ |
7楼2014-04-21 16:27:52
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
8楼2014-04-21 19:59:45
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
9楼2014-04-21 20:49:58
