【答案】应助回帖

思想重于技巧，内涵重于表象

11楼2011-09-05 18:45:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

上传到这里了…………麻烦帮忙看一下…………自洽文件里有原子种类和坐标…………
http://u.115.com/file/dn1p8hdm

12楼2011-09-05 21:07:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

★
uuv2010(金币+1): 欢迎参加交流 2011-09-05 21:31:50

引用回帖:

12楼: Originally posted by 改王之王 at 2011-09-05 21:07:01:
上传到这里了…………麻烦帮忙看一下…………自洽文件里有原子种类和坐标…………
http://u.115.com/file/dn1p8hdm

文件拥有者暂时没有权限分享文件。

思想重于技巧，内涵重于表象

13楼2011-09-05 21:19:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

上传完了…………附上了我写的qha文件，和scf、ph文件…………

14楼2011-09-05 22:14:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

★
uuv2010(金币+1): 欢迎参加交流 2011-09-06 22:43:07

引用回帖:

14楼: Originally posted by 改王之王 at 2011-09-05 22:14:23:
上传完了…………附上了我写的qha文件，和scf、ph文件…………

嫩滴matdyn.in文件写错鸟，最后多了个EOF~所以就都悲剧鸟~

思想重于技巧，内涵重于表象

15楼2011-09-06 07:48:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

引用回帖:

15楼: Originally posted by souledge at 2011-09-06 07:48:40:
嫩滴matdyn.in文件写错鸟，最后多了个EOF~所以就都悲剧鸟~

原来是这样…………大感谢呀…………我再去试试…………呵呵

16楼2011-09-06 08:16:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

引用回帖:

15楼: Originally posted by souledge at 2011-09-06 07:48:40:
嫩滴matdyn.in文件写错鸟，最后多了个EOF~所以就都悲剧鸟~

去了EOF，还是出现上面的错，后来我自己又写了遍计算AlAs的文件，就没有问题…………
又重新把ZnO的计算文件删了重写还是错，我写一下步骤：
1.从tetrahedra文件夹中找到相应的ttrinp_hcp文件，该文件名为ttrinp,把里面的'X'字母用fc文件里的数值代替
2..../QHA/bin/tetra.x 命令来产生k点，并把产生的k点复制到matdyn.in文件夹中，其中matdyn.in的文件我是照着example写的。
3.执行.../bin/matdyn.xmatdyn.out命令，重新计算频率
4.执行../QHA/bin/Partial_phonon_DOS.x < phdos1.in，，
执行../QHA/bin/phonon_dos.x 然后就提示上面帖子出现的错误，不知道是哪有错，我把例子里的frequency拷到ZnO的文件夹里计算../QHA/bin/phonon_dos.x 另外我的ZnO文件夹是放在QHA的example文件夹下进行计算的……

17楼2011-09-06 20:05:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

★ ★ ★
贺仪(金币+3): 多谢指导 2011-09-07 20:47:25
改王之王(金币+16): 大感谢………问题终于解决了…………呵呵 2011-09-09 18:12:14

引用回帖:

17楼: Originally posted by 改王之王 at 2011-09-06 20:05:41:
去了EOF，还是出现上面的错，后来我自己又写了遍计算AlAs的文件，就没有问题…………
又重新把ZnO的计算文件删了重写还是错，我写一下步骤：
1.从tetrahedra文件夹中找到相应的ttrinp_hcp文件，该文件名为ttri ...

感觉你做的很不仔细，关键文件是否存在根本没提及。
我说一下详细步骤：
首先是替换ttrinp中的数据，然后生成K点，把这些K点粘贴到matdyn.in的后面（记得把那个多余的EOF去掉），然后matdyn.x生成frequency，同时生成了matdyn.modes。
之后，保证替换后的ttrinp还在当前文件夹中。创建温度输入文件和德拜温度输入文件，然后创建phdos1.in文件。运行Partial_phonon_DOS.x 然后做好phdos.in文件，运行phonon_dos.x 以上所有文件都在同文件夹中，并且文件名和大小写不要轻易更换。最好完全保持与例子中的相同。

CODE:

$ phonon_dos.x 
 natoms==           4

 irec====         168

  0.75000000000000000     

 Ti1 Ti2 

It seems you have imaginary frequences.                     

 Hopefully you know what you are doing

 nstep====         805

  114.5737  114.5743  147.1030

  122.1805  134.7498  138.2352

  125.8343  131.7529  138.9782

  128.8438  128.8446  139.4173

  125.2058  128.6042  153.5996

  130.8430  131.7825  148.8230

  132.6401  133.3911  146.9807

  132.0160  134.5900  146.1146

  133.6986  137.3024  143.4115

  134.1994  140.3732  140.3739

  118.5830  121.3872  162.3758

 E_min=   0.0000000000000000         E_max=   581.27710000000002     

 nstep====         805

****** input tetrahedra for BZ-integration ******

   6   3  12

  0.000000  0.000000  0.000000

  0.000000  0.577350  0.000000

 -0.333330  0.577350  0.000000

  0.000000  0.000000  0.806500

  0.000000  0.577350  0.806500

 -0.333330  0.577350  0.806500

****** input tetrahedra for BZ-integration ******

********************** generate_tetra *********************

  NT0=           3  NTETMX=        1728

    0.0000   0.0000  -0.3333  -0.3333

    0.0000   0.5774   0.5774   0.5774

    0.0000   0.0000   0.0000   0.8065

    1.0000   1.0000   1.0000   1.0000

  volume of tetrahedron =  0.02587

    0.0000   0.0000   0.0000  -0.3333

    0.0000   0.0000   0.5774   0.5774

    0.0000   0.8065   0.8065   0.8065

    1.0000   1.0000   1.0000   1.0000

  volume of tetrahedron =  0.02587

    0.0000   0.0000  -0.3333   0.0000

    0.0000   0.5774   0.5774   0.5774

    0.0000   0.8065   0.8065   0.0000

    1.0000   1.0000   1.0000   1.0000

  volume of tetrahedron =  0.02587

  total volume of BZ is =   0.0776  omg48=  25.7716

****************** end of generate_tetra ******************

 before integration:  E_min=   0.0000000000000000         E_max=   581.27710000000002     

 Freq,Tot_DOS,DOS==   0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000    

 Freq,Tot_DOS,DOS==  0.7500000      0.2083669E-06  0.8334674E-06  0.2084666E-06  0.4792305E-07  0.1302991E-06  0.3024444E-07

 Freq,Tot_DOS,DOS==   1.500000      0.1666935E-05  0.3333870E-05  0.8342646E-06  0.2019754E-06  0.4941277E-06  0.1381614E-06

思想重于技巧，内涵重于表象

18楼2011-09-07 13:58:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

引用回帖:

18楼: Originally posted by souledge at 2011-09-07 13:58:33:
感觉你做的很不仔细，关键文件是否存在根本没提及。
我说一下详细步骤：
首先是替换ttrinp中的数据，然后生成K点，把这些K点粘贴到matdyn.in的后面（记得把那个多余的EOF去掉），然后matdyn.x生成frequency，同 ...

我就是这样算的呀…………不过我产生partial_DOS文件是9.1M，我只要把例子里的别的材料的frequency沾过来就能往下算，可能是我从matdyn.x那步就算错了？我去了EOF文件还是不行…………把asr改成simple也不行…………太纠结了…………

19楼2011-09-07 18:55:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

不会是程序编译错误吧…………我直接双击的compile…………没有显示编译失败呀…………而且例子也能运行…………