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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】PWSCF中Ti的赝势是不是有问题呢？
181/2返回列表12下一页
查看: 2231  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么网上很多人说本周 12号出结果 +4 瞬息宇宙 2026-08-10 4/200 2026-08-11 02:02 by 旁观2020
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +3 Tide man 2026-08-10 3/150 2026-08-10 21:48 by GZUZM
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
关闭
请填处理意见
关闭