| 查看: 3442 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
justinzhouke新虫 (初入文坛)
|
[交流]
gamma点声子虚频 已有5人参与
|
||
|
我在计算闪锌矿结构的InAs的gamma点频率的时候，结算结果出奇的奇怪。具体如下： Diagonalizing the dynamical matrix q = ( 0.000000000 0.000000000 0.000000000 ) ************************************************************************** omega( 1) = -6.801498 [THz] = -226.873553 [cm-1] ( -0.502968 0.000000 -0.207252 0.000000 -0.106323 0.000000 ) ( 0.755263 0.000000 0.311212 0.000000 0.159655 0.000000 ) omega( 2) = -6.801498 [THz] = -226.873553 [cm-1] ( 0.147751 0.000000 -0.088266 0.000000 -0.526891 0.000000 ) ( -0.221864 0.000000 0.132541 0.000000 0.791185 0.000000 ) omega( 3) = -6.801498 [THz] = -226.873553 [cm-1] ( 0.180077 0.000000 -0.506449 0.000000 0.135339 -0.000000 ) ( -0.270406 0.000000 0.760490 -0.000000 -0.203226 0.000000 ) omega( 4) = -0.136566 [THz] = -4.555354 [cm-1] ( 0.006948 0.000000 -0.160492 0.000000 0.690824 0.000000 ) ( 0.006906 0.000000 -0.159518 0.000000 0.686631 0.000000 ) omega( 5) = -0.136566 [THz] = -4.555354 [cm-1] ( 0.617358 0.000000 0.341420 0.000000 0.073109 0.000000 ) ( 0.613612 0.000000 0.339348 0.000000 0.072665 0.000000 ) omega( 6) = -0.136566 [THz] = -4.555354 [cm-1] ( -0.349090 0.000000 0.600599 -0.000000 0.143042 0.000000 ) ( -0.346972 0.000000 0.596954 -0.000000 0.142174 0.000000 ) ************************************************************************** gamma点虚频居然如此之大。我以及做了收敛性测试和结构优化。 但是为何gamma点虚频居然如此之大？。谁能帮助我一下。  输入文件如下： scf： &control calculation = 'scf', restart_mode= 'from_scratch', pseudo_dir = './', outdir = './tmp', prefix='inas', tstress = .true. tprnfor = .true. / &system ibrav= 0, celldm(1)=11.450120305668, nat= 2, ntyp= 2, ecutwfc =48, / &electrons diagonalization='cg' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 5.0d-12 / ATOMIC_SPECIES In 111.818 In.pbe-hgh.UPF As 74.920 As.pbe-hgh.UPF ATOMIC_POSITIONS crystal In 0.00 0.00 0.00 As 0.25 0.25 0.25 K_POINTS automatic 9 9 9 0 0 0 CELL_PARAMETERS alat 0.000000000 0.508881713 0.508881713 0.508881713 0.000000000 0.508881713 0.508881713 0.508881713 -0.000000000 ph： phonons of InAs &inputph tr2_ph=1.0d-14, prefix='inas', epsil=.true. amass(1) =111.818, amass(2) =74.920, outdir='./tmp', fildyn='inas.dyn', / 0.0 0.0 0.0 自恰输出结果为： 。。。。 ! total energy = -16.65528050 Ry Harris-Foulkes estimate = -16.65528050 Ry estimated scf accuracy < 4.8E-12 Ry The total energy is the sum of the following terms: one-electron contribution = 1.42930748 Ry hartree contribution = 1.95307784 Ry xc contribution = -4.59653835 Ry ewald contribution = -15.44112747 Ry convergence has been achieved in 34 iterations Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 atom 2 type 2 force = 0.00000000 0.00000000 0.00000000 Total force = 0.000000 Total SCF correction = 0.000000 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= 0.00 0.00000001 0.00000000 0.00000000 0.00 0.00 0.00 0.00000000 0.00000001 -0.00000000 0.00 0.00 -0.00 0.00000000 -0.00000000 0.00000001 0.00 -0.00 0.00 。。。。 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
有谁用过 phonon 这个软件，希望交流下使用经验
已经有13人回复
- 声子谱计算在Gamma点出0.2THz虚频，可以忽略吗？
已经有52人回复
- 石墨烯声子谱Gamma点附近虚频，怎么处理
已经有5人回复
- 用VASP-DFPT+phonopy计算石墨烯声子谱求助
已经有37人回复
- 声子谱！！谢谢！！
已经有4人回复
- Gamma点处出现比较大的虚频说明什么
已经有3人回复
- 你们怎么处理Gamma的虚频的？
已经有4人回复
- 声子谱 虚频 求助
已经有10人回复
- degauss的收敛
已经有13人回复
- 哭了，各种虚频！
已经有23人回复
- 用自己生成的赝势算能带效果很好，算声子却一团糟
已经有11人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- 怎么消除虚频，细节
已经有10人回复
- 求助如何分析castep声子计算的结果
已经有7人回复
- 自旋体系的声子计算应该用什么软件呢？
已经有5人回复
7楼2018-02-01 21:14:49
2楼2015-09-04 10:55:04
justinzhouke
新虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 3.5
- 红花: 1
- 帖子: 24
- 在线: 6.5小时
- 虫号: 3513863
- 注册: 2014-11-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2015-09-04 14:31:11
magicmonk
至尊木虫 (著名写手)
- 应助: 172 (高中生)
- 金币: 11451.3
- 红花: 24
- 帖子: 1124
- 在线: 3108.8小时
- 虫号: 1191315
- 注册: 2011-01-17
- 专业: 理论和计算化学
4楼2015-09-05 11:28:12
10