24小时热门版块排行榜    

查看: 833  |  回复: 25
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 positron 的 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

positron

木虫 (职业作家)

[交流] 【求助】体系添加电荷的描述问题

也不知道具体发到哪个版块啦,就借地一问吧


要对一个周期性的系统添加一个带电粒子,然后看其后演化

系统的势能(晶体势能,包括离子势和关联势)已经得到,现在就是编写程序计算演化过程,但是不清楚加入带电粒子薛定谔方程如何描述,也就是初态问题,如何定义t0时刻在体系的某个位置存在一个电荷,然后演化

感觉是平面波,但是没有啥依据,盼看到的朋友能解答。
回复此楼

» 收录本帖的淘帖专辑推荐

good resources

» 猜你喜欢

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

positron

木虫 (职业作家)

想了想,是不是φ=Aexp(-i(k1x+k2y+k3z))的形式呢
2楼2009-08-26 09:42:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
Originally posted by zhxd1028 at 2009-8-26 11:48:
带电粒子的演化?不明白啊,你是想做什么演化呢?电荷移动(就是电子移动吧)还是电荷重新分布?
薛定谔方程描述的是在势场下电子的运动方程啊,描述对象就是电子(也可看成一个电荷),这不就是电荷的演化吗?

应该是这样的

在系统中添加一个电荷(电子),然后看其运动,感觉就是最基本的运动方程啊

但是要给方程一个初始条件才能解,所以才问初始条件

很不懂,所以说不清,原谅
4楼2009-08-26 12:12:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

初始状态可是三维空间自由电子的波函数的形式?

三维空间自由电子波函数的形式是啥我也没有找到,呵呵,手边的书很少
网上也还没找到
5楼2009-08-26 12:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
Originally posted by zhxd1028 at 2009-8-26 22:00:

明白你的意思了,添加一个电子后电子的演化主要看体系的势函数是什么样的,如果势函数变化很缓(最简单的例子就是一个常数),那么添加的电荷就类似于自由电子,可以用平面波展开。这是固体计算中经常用的方法。 ...

那么说是我把问题想简单啦

势函数就是各个离子实的库伦势,然后在加上电子关联势

这个展开看来还是好复杂的。。。

哎,初生牛犊不怕虎,乱写程序啊,没有预先到困难
8楼2009-08-26 22:32:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
Originally posted by zhxd1028 at 2009-8-26 22:00:

明白你的意思了,添加一个电子后电子的演化主要看体系的势函数是什么样的,如果势函数变化很缓(最简单的例子就是一个常数),那么添加的电荷就类似于自由电子,可以用平面波展开。这是固体计算中经常用的方法。 ...

那用平面波展开或者局域化的基组展开,展开系数又该如何确定呢

还有即便得到了波函数的形式,方程又该如何解呢?没有初始条件啊

初始条件的如何描述或者引入,也是我一直搞不明白的

身边都是搞核物理的,我也是闲暇慢慢学习,所以问题很多

盼再指教
9楼2009-08-26 22:37:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
Originally posted by zhxd1028 at 2009-8-27 11:38:

展开系数的确定是解决电子运动方程的核心问题了,现在已有很多成熟的办法。建议看看Kresse写得vasp的理论基础,还有Payne的文章,就是几个软件的理论基础,他们的方法已经成功的应用在软件中了。解决得部分问题 ...

那我就从最简单的开始吧

忍不住还想问一个,就是比如考虑简单金属,例如Na,体心立方结构,计算一个晶胞或者2*2*2的胞,在其中放一个氢离子H+,解HF方程

好想知道这个过程和处理的思路

要是你碰到这种问题,你会如何的处理和考虑呢?

要,要多少平面波展开才行呢?
12楼2009-08-27 13:46:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

等着朋友的回答

呵呵,我就在追加10个金币吧,你也辛苦啦
13楼2009-08-27 13:55:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

★ ★
fegg7502(金币+2,VIP+0):鼓励交流 8-30 17:01
忘了写一句,

就是比如考虑简单金属,例如Na,体心立方结构,计算一个晶胞或者2*2*2的胞,在其中放一个氢离子H+,解HF方程

晶胞内氢离子H+感受到的势能+关联能的分布我通过别的近似已经能得到

之后假如要是你来处理这种HF方程的求解问题,你会如何考虑?

经过你前面的解答,知道了点门道,不过还是劳驾你就这个例子给我缕缕思路
14楼2009-08-27 14:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
Originally posted by zhxd1028 at 2009-8-27 16:12:

那你想知道什么结果呢?是这样一个体系的的电子分布?还是H+(其实就是个质子)的位置,还是体系的总能?或者什么其他的。既然从最简单开始就要只针对一个问题,这样就可以抛弃其他因素的影响了。

H+的位置

这个先考虑吧,看看各处出现的几率吧
16楼2009-08-27 16:54:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 positron 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见