| 查看: 5227 | 回复: 25 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
frank_zhan金虫 (正式写手)
|
[求助]
怎么利用PWSCF计算差分电荷密度
|
||
|
各位亲虫们， 我想计算一个体系的差分电荷密度，但是不知道过程是怎样的？比如我计算一个SWNT上修饰金原子（Au）的团簇。我想计算这两者之间的电荷转移量。理论上是利用SWNT@Au体系的总电荷密度减掉SWNT和Au团簇的电荷密度。 请求虫友们指导个详细的过程，必有重谢！ 谢谢。 |
回复此楼» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
frank_zhan
金虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 819.6
- 散金: 130
- 红花: 3
- 帖子: 577
- 在线: 136.4小时
- 虫号: 688836
- 注册: 2009-01-06
- 性别: GG
- 专业: 无机非金属类光电信息与功
|
你好，仁兄，最近在做那个查分电荷密度的计算，发现一些问题，请教这位仁兄。我的输入文件如下（省掉了原子的坐标部分）： 第一步，计算3个体系的的scf: SWNT-Au.scf.in &control calculation = 'scf' prefix= 'SWNT-Au' pseudo_dir = '/home/pseudo', outdir= 'mypath/SWNT-Au/tmp' / &system .... / SWNT.scf.in &control calculation = 'scf' prefix= 'SWNT' pseudo_dir = '/home/pseudo', outdir= 'mypath/SWNT-Au/tmp' / &system ... / Au.scf.in &control calculation = 'scf' prefix= 'Au' pseudo_dir = '/home/pseudo', outdir= 'mypath/SWNT-Au/tmp' / &system ... / 第二部. 计算3个体系的电荷密度(using pp.x): SWNT-Au.rho.in &inputpp prefix= 'SWNT-Au' outdir= 'mypath/SWNT-Au/tmp' filplot = 'SWNT-Au' plot_num= 0 spin_component=0 / &plot nfile = 1 filepp(1) = 'SWNT-Au' weight(1) = 1.0 iflag = 3 output_format = 5 fileout = 'SWNT-Au.rho.dat' / SWNT.rho.in &inputpp prefix= 'SWNT' outdir= 'mypath/SWNT-Au/tmp' filplot = 'SWNT' plot_num= 0 spin_component=0 / &plot nfile = 1 filepp(1) = 'SWNT' weight(1) = 1.0 iflag = 3 output_format = 5 fileout = 'SWNT.rho.dat' / Au.rho.in &inputpp prefix='Au' outdir= 'mypath/SWNT-Au/tmp' filplot = 'Au' plot_num= 0 spin_component=0 / &plot nfile = 1 filepp(1) = 'Au' weight(1) = 1.0 iflag = 3 output_format = 5 fileout = 'Au.rho.dat' / 第三步.计算差分电荷密度(using pp.x): SWNT-Au.rho_diff.in &inputpp / &plot nfile = 3 filepp(1) = 'SWNT-Au.rho.dat' weight(1) = 1.0 filepp(2)= 'SWNT.rho.dat' weight(2)= -1 filepp(3)= 'Au.rho.dat' weight(3)= -1 iflag = 3 output_format = 5 fileout = 'SWNT-Au.rho_diff.dat' / 所以的第一步和第二部分的计算都成功了，但是计算最后的差分电荷密度的的时候，结果SWNT-Au.rho_diff.out: Reading header from file SWNT-Au.rho.dat. 也就是根本没有进行计算就退出了。 请仁兄给点建议，金币的话可以追加的，谢谢。 |
22楼2013-07-03 21:15:35
frank_zhan
金虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 819.6
- 散金: 130
- 红花: 3
- 帖子: 577
- 在线: 136.4小时
- 虫号: 688836
- 注册: 2009-01-06
- 性别: GG
- 专业: 无机非金属类光电信息与功
2楼2012-07-09 11:44:23
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
frank_zhan: 金币+10, ★有帮助, 谢谢指导，我先试试，假如可以的话，我会继续给金币的，哈哈。 2012-07-09 14:28:45
uuv2010: 金币+1, 多谢提示 2012-07-11 11:44:59
感谢参与，应助指数 +1
frank_zhan: 金币+10, ★有帮助, 谢谢指导，我先试试，假如可以的话，我会继续给金币的，哈哈。 2012-07-09 14:28:45
uuv2010: 金币+1, 多谢提示 2012-07-11 11:44:59
|
先算SWNT@Au的电荷密度，再算SWNT@Au-Au的电荷密度，再算SWNT@Au-SWNT的电荷密度，然后用前者减去后两者。 过程： scf1.in，->rho1.in->得到ppn1.rho scf2.in, ->rho2.in得到ppn2.rho scf3.in, ->rho3.in得到ppn3.rho 最后算dif.rho.in: / &plot nfile = 3, filepp(1) = 'ppn1.fpl', filepp(2) = 'ppn2.fpl', filepp(3) = 'ppn3.fpl', weight(1) = 1.0, weight(2) = -1.0, weight(3) = -1.0, iflag = 3, output_format = 5, fileout = 'ppn_diff.rho', / 得到差分密度文件ppn_diff.rho和ppn_diff.xsf，其中ppn_diff.xsf可以用xcrysden画差分电荷密度图。 |
3楼2012-07-09 12:50:28
frank_zhan
金虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 819.6
- 散金: 130
- 红花: 3
- 帖子: 577
- 在线: 136.4小时
- 虫号: 688836
- 注册: 2009-01-06
- 性别: GG
- 专业: 无机非金属类光电信息与功
4楼2012-07-09 14:37:15