24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1194  |  回复: 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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

souledge

专家顾问 (著名写手)

【答案】应助回帖


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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁

引用回帖:
2楼: Originally posted by souledge at 2011-11-01 18:29:12:
如果是3D的,比如选择output_format = 5,是不需要任何矢量的~
否则,就把ei(j)看作是以原点为起点,与欲绘面的边界平行的几个矢量就可以了。不过要求是ei(j)必须正交,也就是ei x ej = 0(叉乘)。
而xi(j)则 ...

非常感谢,还是不行。
&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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖


minmin_0082003(金币+5): 解决鸟,非常感谢 2011-11-02 18:59:52
贺仪(金币+1): 多谢指教! 2011-11-02 21:08:45
引用回帖:
3楼: Originally posted by minmin_0082003 at 2011-11-02 16:47:25:
非常感谢,还是不行。
&INPUTPP
                      prefix = "Test" ,
                      outdir = "/home/gaomin/scr" ,
                     filplot = "charge-te ...

第一部分是&inputpp,第二部分是&plot~不是&input~
思想重于技巧,内涵重于表象
4楼2011-11-02 17:36:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
信息提示
请填处理意见