| 查看: 1417 | 回复: 3 | ||
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 ... 请高手指点，是什么地方设置的不当，该如何修正，不胜感激 |
回复此楼» 猜你喜欢
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
感觉是下周放榜了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助用MS计算声子谱出错的问题
已经有15人回复
- 求助：有限位移法算声子谱
已经有5人回复
- 【求助】castep计算声子谱，为什么不能进行振动分析
已经有9人回复
- 【求助】abinit声子计算
已经有7人回复
- 【求助】使用vasp进行带心冻结声子计算的几个问题~~~
已经有9人回复
- 【求助】求计算时间参考
已经有9人回复
- 【求助】castep计算声子谱求助
已经有9人回复
- 【求助】请教CASTEP计算声子谱后出现虚频
已经有6人回复
- 【求助】声子谱
已经有5人回复
- 【求助】有人用有限差分方法计算过G点声子谱吗？
已经有10人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
- 【求助】为什么一算声子谱就不出结果了呢？
已经有6人回复
- 【求助】问大家一个基础概念，什么是晶体的简正频率？在声子谱中能读出简正频率吗？
已经有4人回复
- 【求助】castep计算声子谱
已经有3人回复
- 【求助成功】声子能带和声子谱的联系
已经有16人回复
- 【求助】声子谱的计算速度
已经有11人回复
- 【求助】声子谱计算
已经有5人回复
- 【求助】声子谱
已经有8人回复
- 【求助】声子谱
已经有6人回复
- 【求助】用MS中的castep计算声子振动模频率，具体怎么算
已经有10人回复
- 【求助】怎样分析声子谱？
已经有18人回复
fengshiquan
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 5831.8
- 散金: 80
- 红花: 5
- 帖子: 599
- 在线: 1569.4小时
- 虫号: 884725
- 注册: 2009-10-26
- 专业: 原子和分子物理
2楼2013-04-07 10:20:33
3楼2013-04-07 16:08:15
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fengshiquan: 金币+10, ★有帮助, 呵呵，楼上已经提到了，我去改了这个参数，貌似确实是这个问题，还是要谢谢你的帮助，随便问下，改了这个参数出现了较大虚频，需要修改那个收敛精度，上面的in文件中有几个控制收敛精度的参数，不知道调哪个可以消去须频，谢谢 2013-04-08 10:25:14
franch: 金币+2, 谢谢回帖交流，， 2013-04-08 21:13:58
感谢参与，应助指数 +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
25