| 查看: 1414 | 回复: 3 | ||||
minmin_0082003金虫 (正式写手)
莫愁莫愁
|
[求助]
计算3D电子密度，如何确定相关矢量呢？
|
|
e1(i), e2(i), e3(i), i=1,3 的值 （3D vectors which determine the plotting parallelepiped） 如何确定的？ 数学问题，忘光了。求助 算了一个2D的，含有3个原子，假设 ABC 三个原子，以B为起点，A减去B为矢量1,C减去B为矢量2 但是总是会出现错误“ Writing data to file charge-test Message from routine chdens: namelist plot not found or invalid, exiting” 输入文件 “&INPUTPP prefix = "Test" , outdir = "/home/scr" , filplot = "charge-test" , plot_num = 0, spin_component = 0 , / &PLOT nfile = 1, filepp(1) = "charge-test" iflag = 2 output_format = 3 fileout = "output.XSF" e1(1) = 0.4, e2(1) = 1.03, e3(1) = 0.55, e1(2) = -4.5, e2(2) = -2.74, e3(2) = 0.5, nx = 50, ny =50 x0(1)= 11.373864755 x0(2)=4.507957875 x0(3)= 6.846860808 ” [ Last edited by minmin_0082003 on 2011-11-1 at 16:47 ] |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请看这个电子密度图是如何做的呀？
已经有6人回复
- 如何做电子密度差分图
已经有16人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有106人回复
- 【求助】CASTEP计算的电子密度的单位？
已经有17人回复
- 【求助】利用Gaussian计算Nd2Fe14B单胞的电子密度、单点能用什么方法和基组，谢谢！
已经有3人回复
- 【求助】电子态密度计算dos图中feimi能级的确定
已经有7人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】Dmol3计算outmol分析电子密度，谢谢！
已经有11人回复
- 【求助】计算电子密度图【已解决】
已经有14人回复
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★
uuv2010(金币+1): 多谢提示 2011-11-01 18:31:20
uuv2010: 专家考核存档 2011-11-01 18:31:26
minmin_0082003(金币+10): 非常感谢，但是我试了3d， 是一样的错误，太奇怪了 2011-11-02 16:45:23
uuv2010(金币+1): 多谢提示 2011-11-01 18:31:20
uuv2010: 专家考核存档 2011-11-01 18:31:26
minmin_0082003(金币+10): 非常感谢，但是我试了3d， 是一样的错误，太奇怪了 2011-11-02 16:45:23
|
如果是3D的，比如选择output_format = 5，是不需要任何矢量的～ 否则，就把ei(j)看作是以原点为起点，与欲绘面的边界平行的几个矢量就可以了。不过要求是ei(j)必须正交，也就是ei x ej = 0（叉乘）。 而xi(j)则表示欲绘面上几个“边”的交点坐标。 比如，这是交点在(0, 0)，两个边的矢量分别是(1, 0, 0)和(0, 1, 0)的TiC的(100)面Fermi面附近电荷密度的图。当然，这个晶胞进行了剪切应变，所以电荷密度看起来比较奇怪～  |
2楼2011-11-01 18:29:12
minmin_0082003
金虫 (正式写手)
莫愁莫愁
- 1ST强帖: 1
- 应助: 2 (幼儿园)
- 贵宾: 0.08
- 金币: 921.3
- 散金: 149
- 帖子: 910
- 在线: 218.4小时
- 虫号: 366989
- 注册: 2007-05-10
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
|
非常感谢，还是不行。 &INPUTPP prefix = "Test" , outdir = "/home/gaomin/scr" , filplot = "charge-test" , plot_num = 0, spin_component = 0 , / &INPUT nfile = 1 , filepp(1) = "charge-test", weight(1) = 1.0, fileout = output , iflag = 3 , output_format = 5 , / filplot 是出来了，但是还是提示“ Calling punch_plot, plot_num = 0 Writing data to file charge-test Message from routine chdens: namelist plot not found or invalid, exiting ” |
3楼2011-11-02 16:47:25
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
4楼2011-11-02 17:36:31
