24小时热门版块排行榜    

查看: 1243  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 一志愿郑大070303,338分,求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见