| 查看: 1411 | 回复: 8 | ||||
785524519金虫 (正式写手)
|
[求助]
求助！关于Au将5d轨道电子考虑成芯电子
|
|
小木虫的各位虫有麻烦过来看一下！比较重要的问题！！！！ 请问各位虫友，计算Au时siesta的SZP默认选5d10 6s1 6p轨道，如果将d考虑为芯电子，计算只算s p轨道 将要如何在fdf写自己的轨道？并且修改Au的赝势psf文件？ 如果有这样的赝势能否不吝提供下。对于修改赝势和自己写价电子轨道小弟并没有写过，比较迷茫。 感谢大家！！！ |
回复此楼» 猜你喜欢
326求调剂
已经有8人回复
-
333求调剂
已经有7人回复
-
一志愿东华大学控制学硕320求调剂
已经有3人回复
-
考研化学学硕调剂，一志愿985
已经有7人回复
-
化学工程321分求调剂
已经有21人回复
-
0703化学调剂
已经有15人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有14人回复
-
326求调剂
已经有4人回复
-
311求调剂
已经有13人回复
-
311求调剂
已经有7人回复
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3511.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
785524519(蛋蛋小童鞋代发): 金币+30, 谢谢参与，小木虫鼓励积极回答问题者 2017-04-13 18:29:30
蛋蛋小童鞋: 应助指数+1 2017-04-13 18:29:39
785524519(蛋蛋小童鞋代发): 金币+30, 谢谢参与，小木虫鼓励积极回答问题者 2017-04-13 18:29:30
蛋蛋小童鞋: 应助指数+1 2017-04-13 18:29:39
|
修改芯电子的定义并非在fdf，而是生成psf的时候给定的参数，而非直接简单修改psf。生成psf你可以使用ATOM（https://departments.icmab.es/lee ... Code/downloads.html）或者APE。 赝势参数的变化会导致赝原子轨道形状发生变化，原来的赝原子基组不再是最优选择，原则上应该重新优化基组，在fdf有专门的block来自定义基组，但具体细节我不是很清楚。我个人目测不考虑d半芯态可能结果会不好。 |
7855245192楼2017-04-12 19:55:03
785524519
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2240.2
- 红花: 2
- 帖子: 329
- 在线: 157小时
- 虫号: 4078654
- 注册: 2015-09-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2017-04-12 22:08:56
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3511.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2017-04-12 22:18:04
785524519
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2240.2
- 红花: 2
- 帖子: 329
- 在线: 157小时
- 虫号: 4078654
- 注册: 2015-09-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2017-04-12 22:28:51
785524519
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2240.2
- 红花: 2
- 帖子: 329
- 在线: 157小时
- 虫号: 4078654
- 注册: 2015-09-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2017-04-13 21:08:09
785524519
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2240.2
- 红花: 2
- 帖子: 329
- 在线: 157小时
- 虫号: 4078654
- 注册: 2015-09-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2017-04-13 21:10:23
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3511.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
8楼2017-04-13 21:43:21
785524519
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2240.2
- 红花: 2
- 帖子: 329
- 在线: 157小时
- 虫号: 4078654
- 注册: 2015-09-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
您好这边又有问题要请教您。是关于siesta得到的电荷密度分布的，目前只是定性分析所以建了一个21个原子，其中20个Au,Ni放在中间位置。然后计算得到的哈密顿量H画出电荷分布图。但是问题来了，计算的结果和siesta中mulliken电荷分部，Hirsh(手册中提到的分布)不一样！最主要的是左右两边的电荷密度不同。siesta中得到的左右两边原子带正电荷，但是按道理来说，左右两边是真空层没有东西，电荷会往最左右两边电荷移动.。想请问下siesta中得到的Mulliken分布难道不是按照体系的哈密顿量得出的吗？请问您有什么好的建议吗。十分感谢！下面是算的结果  Au20Ni1self.png  Au20Ni1Mulliken.jpg |
9楼2017-04-16 21:39:39
