应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【求助】体系添加电荷的描述问题
172/2返回列表上一页12
查看: 833  |  回复: 25
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 positron 的 5 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

positron

木虫 (职业作家)

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

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


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

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

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

» 猜你喜欢
1楼 2009-08-26 06:56:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
引用回帖:
Originally posted by zhxd1028 at 2009-8-27 19:35:


首先我觉得是可以确定H+的位置，而不会是分布几率。
我想大体思路是这样：
1.给定H+的位置坐标以及Na、Cl 原子的坐标，计算结构的电荷分布
2.在给定的电荷分布下计算体系各个原子的Hellmman-Fennyman forc ...

越学习越觉得自己不知道的东西多

有动力，但也头痛，呵呵
赞一下 回复此楼
18楼2009-08-28 15:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
引用回帖:
Originally posted by zhxd1028 at 2009-8-27 19:35:


首先我觉得是可以确定H+的位置，而不会是分布几率。
我想大体思路是这样：
1.给定H+的位置坐标以及Na、Cl 原子的坐标，计算结构的电荷分布
2.在给定的电荷分布下计算体系各个原子的Hellmman-Fennyman forc ...

那氢离子我可以用平面波代替

然后解ＨＦ方程吗？

我看有文献是一个平面波基组，我若假定离子的能量最低（也就是零？），然后就是一个常数可行？

我想问题简化再简化，这种最简单的情况可行不？
若是一个电子呢，我假定能量为零，然后平面波表示，常数可行？
赞一下 回复此楼
19楼2009-08-28 15:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
★ ★
fegg7502(金币+2,VIP+0):鼓励交流 8-30 17:01
引用回帖:
Originally posted by zhxd1028 at 2009-8-28 18:54:

“那氢离子我可以用平面波代替”我觉得这个思路很有意思，对，就把它当做一个在势场中运动的正电荷。然后接HF-方程。就是不知道解出来后的意义是什么，可以理解成质子在势场中的分布纪律吧。你太有才了，这都能 ...

我的意思是假定氢离子能量（动能）最小，也就是是说把它放在了势场中的某个位置（没有动能的引入势场中的某个位置），感觉这样简化是最简单

还有解HF方程时，若动能取常数（平面波），也就是（-h^2/(4*pi^2*2*m）ψ用一个常数表示即可？（这种理解和处理可行不？）

还有在晶体中，势场是三维的，3D偏微分方程解着有点繁，这时可否这样考虑：
其势场是在一系列Z值（晶胞的Z方向）对应的二维势场分布

这样我就可以轻松的用matlab的eig函数求二维矩阵的特征向量和特征值

每一组特征向量归一化---》得到空间分布

还是一定要解3D 偏微分方程，得到空间的特征向量和特征值，得到空间分布？
赞一下(10人) 回复此楼
21楼2009-08-29 01:30:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
★
freshgirl(金币+1,VIP+0):鼓励一下自己纠正问题~ 9-2 11:58
引用回帖:
Originally posted by zhxd1028 at 2009-8-30 18:52:


既然动能取常数，那么常数和0没多大区别。那么好像你的意思是不要动能项了，只考虑在二维势能下H+的分布，是吗？这样处理出来的H+的波函数没意义吧。

而且，H+在势场下平衡后是总能量最小，不是动能最小。 ...

我理解错了，应该是总能量最小的情况下的H+的波函数
赞一下(10人) 回复此楼
24楼2009-08-30 22:34:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
引用回帖:
Originally posted by zhxd1028 at 2009-8-30 18:52:


既然动能取常数，那么常数和0没多大区别。那么好像你的意思是不要动能项了，只考虑在二维势能下H+的分布，是吗？这样处理出来的H+的波函数没意义吧。

而且，H+在势场下平衡后是总能量最小，不是动能最小。 ...

那在简单的考虑简单的，波函数用两个平面波展开：

ψ=a1p1+a2p2

这时得到：
    (a1p1)^2+U*(a1p1)=E1*(a1p1)
{
    (a2p2)^2+U*(a2p2)=E2*(a2p2)

再加上归一化条件的方程，感觉还是未知数多，方程少
赞一下 回复此楼
25楼2009-08-30 23:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)
再细致看了你的建议

我好像还有一个条件没有用上，那就是周期性边界条件问题

这个又要怎么用呢

对于H+在晶体中的运动来说，要施加怎样的周期性边界条件呢？
赞一下 回复此楼
26楼2009-08-30 23:12:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 positron 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭