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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】PWSCF中Ti的赝势是不是有问题呢？
181/2返回列表12下一页
查看: 2346  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

identation

金虫 (著名写手)

[交流] 【求助】PWSCF中Ti的赝势是不是有问题呢？ 已有7人参与

Ti的组态应该是3s2 3p6 3d2 4s2，可是pwscf里面库中的赝势的占据态是3s2 3p6 3d1 4s2(下图)，这是怎么回事呢？谢谢大家！
回复此楼

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

pwscf

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

netx_ray

木虫 (小有名气)
★
identation(金币+5): 2010-06-19 11:37:06
hedaors(金币+1):谢谢交流 2010-06-19 22:03:00
一个简单的方法，就是你手动把3d的那个占据数改成2.0，然后计算看看是不是结果一样。

楼主试试看看
赞一下(1人) 回复此楼
2楼2010-06-19 11:33:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by netx_ray at 2010-06-19 11:33:10:
一个简单的方法，就是你手动把3d的那个占据数改成2.0，然后计算看看是不是结果一样。

楼主试试看看

这个可以改吗？
回复此楼
3楼2010-06-19 11:37:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

netx_ray

木虫 (小有名气)
我现在身边没有装pwscf的电脑 你改改试试吧 建个面心立方的Ti计算一下看看呗

据我理解这个地方的占据数 应该没有影响的。

我也是一知半解，试试吧
赞一下 回复此楼
4楼2010-06-19 11:41:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★
identation(金币+1):谢谢参与
identation(金币+20): 2010-06-19 16:23:37
cenwanglai(金币+2):谢谢回复 2010-06-20 20:01:44
引用回帖:
Originally posted by identation at 2010-06-19 08:26:32:
Ti的组态应该是3s2 3p6 3d2 4s2，可是pwscf里面库中的赝势的占据态是3s2 3p6 3d1 4s2(下图)，这是怎么回事呢？谢谢大家！

这样是按Ti的离子组态来构造赝势的，并考虑了半芯态的3s2 3p6，这是可以的。构造赝势时，并不一定需要按中性原子的电子组态来构造。在应用赝势之前，自己一定要测试一下赝势的质量，比如对它的单质体材料或分子，以及最常见的化学物进行计算。
赞一下(2人) 回复此楼
5楼2010-06-19 15:10:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by valenhou001 at 2010-06-19 15:10:09:


这样是按Ti的离子组态来构造赝势的，并考虑了半芯态的3s2 3p6，这是可以的。构造赝势时，并不一定需要按中性原子的电子组态来构造。在应用赝势之前，自己一定要测试一下赝势的质量，比如对它的单质体材料或分 ...

谢谢老师，那布局分析的时候，Ti电荷的计算，是用11减去占据数呢还是用12减去占据数呢？
另外，在scf输出文件中又显示zval=12.0，不是11.0，那Ti布局分析的时候怎么根据占据数求价态呢？
谢谢老师
     PseudoPot. # 2 for Ti read from file Ti.pbe-sp-van_ak.UPF
     Pseudo is Ultrasoft, Zval = 12.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  851 points,  6 beta functions with:
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
                l(5) =   2
                l(6) =   2
     Q(r) pseudized with  8 coefficients,  rinner =    1.000   1.000   1.000
                                                       1.000   1.000
赞一下 回复此楼
6楼2010-06-19 16:25:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)
★ ★ ★
identation(金币+1):谢谢参与
identation(金币+10): 2010-06-19 19:48:02
cenwanglai(金币+2):计算的时候考虑到12个电子,这个程序怎么实现的呢?检查一下总电子树,可能就少一个电子呢. 2010-06-20 20:03:10
这个没错的.计算的时候肯定还是考虑到12个电子.但单原子和构成固体的院子是不一样的.d轨道只有一个也是可以的.还有一个电子在原子间自由运动形成金属键.你可以参考其他的物质,比如氢.有些氢的赝势只有0.5个电子.赝势文件不单单要考虑到单原子情况,也需要考虑到固体/化合物中的真实情况,这样做出来的赝势才能有针对性的符合你的要求.万能的赝势是不存在的.
赞一下(2人) 回复此楼
7楼2010-06-19 18:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by condensed at 2010-06-19 18:26:59:
这个没错的.计算的时候肯定还是考虑到12个电子.但单原子和构成固体的院子是不一样的.d轨道只有一个也是可以的.还有一个电子在原子间自由运动形成金属键.你可以参考其他的物质,比如氢.有些氢的赝势只有0.5个电子.赝 ...

那在布局分析的时候还要用12-占据数得到Ti的有效charge？
赞一下 回复此楼
8楼2010-06-19 19:48:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dapaosheng
★
identation(金币+1):谢谢参与
9楼2010-06-19 20:24:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★
identation(金币+1):谢谢参与
identation(金币+10): 2010-06-19 23:52:53
cenwanglai(金币+2):谢谢回复 2010-06-20 20:03:53
引用回帖:
Originally posted by identation at 2010-06-19 19:48:40:

那在布局分析的时候还要用12-占据数得到Ti的有效charge？

是的，用12.这里的电子占据只是说这个电子没有在考虑的轨道内，而可能在更高能级轨道上，所以没有写出来。所以价电子数目仍然是12.
另外，不知道你用的是哪个Ti的势函数？Ti.pbe-??的那个？或者你可以根据自己的需要，按照这个文件给出的能量截断什么的自己构造一个。
赞一下(2人) 回复此楼
10楼2010-06-19 21:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭