24小时热门版块排行榜    

查看: 2411  |  回复: 25

小竹子AA

金虫 (正式写手)

[求助] 晶体采用单晶面和双晶面slab模型模拟表面能中的问题 已有3人参与

我模拟的结构是六方晶系的,模拟了晶体的(001)晶面,分别采用了单晶面和双晶面slab模型分别模拟,然后单晶面slab模型模拟的结果约是双晶面slab模型模拟结果的两倍,而这两种模型在表面能公式中正好有分母部分是不是乘2的区别,出现这种情况的原因是什么?有人分别用两个模型算出来是结果一致的么?百思不得其解,先谢谢各位帮助(注:两个模型的真空层厚度是相同的,单晶面slab模型采用了下部原子固定,双晶面slab模型采用中部原子固定;真空层厚度测试表明10埃就够了,在两个模型中均采用了16埃的真空层厚度;单晶面slab模型得到的xtd文件动画中,看到被固定的原子有一点移动,但相对于上部未被固定的原子的移动很小)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

感谢参与,应助指数 +1
蛋蛋小童鞋: 应助指数-1 2016-09-05 21:41:09
单晶面slab模型的表面能是如何计算的?扣掉了底部表面的表面能了吗?
2楼2016-08-23 09:29:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小竹子AA

金虫 (正式写手)

引用回帖:
2楼: Originally posted by c19940408 at 2016-08-23 09:29:42
单晶面slab模型的表面能是如何计算的?扣掉了底部表面的表面能了吗?

就是把下半部分的原子都固定,弛豫过程中这部分原子不能动只有上部原子弛豫,然后得到的就是单晶面slab模型的表面能了,我是这么理解的。怎么扣掉底部的表面能啊?
3楼2016-08-23 09:49:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liupanjin

铜虫 (小有名气)

关注你吧,每次都能提醒,可是我不懂啊

发自小木虫Android客户端
决战2014考研,如今取得胜利!去争取更大胜利
4楼2016-08-23 17:09:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小竹子AA

金虫 (正式写手)

引用回帖:
4楼: Originally posted by liupanjin at 2016-08-23 17:09:24
关注你吧,每次都能提醒,可是我不懂啊

还有人关注我呀真开森那先不要关注我了,我也是半瓶子水,关注我没啥意义撒,不过先谢谢你让我这么开森哈
5楼2016-08-23 20:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liupanjin

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by 小竹子AA at 2016-08-23 20:40:43
还有人关注我呀真开森那先不要关注我了,我也是半瓶子水,关注我没啥意义撒,不过先谢谢你让我这么开森哈...

没事,先关注着吧

发自小木虫Android客户端
决战2014考研,如今取得胜利!去争取更大胜利
6楼2016-08-24 06:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhao8919

银虫 (初入文坛)

【答案】应助回帖

★ ★
蛋蛋小童鞋: 金币+2, 应助指数+1, 谢谢参与,小木虫鼓励积极回帖 2016-09-05 21:41:46
猜想原因:
1、单面slab、双面slab模型,优化后原子位置改变不大,所以能量改变也就不大,这样分母除以二之后出现2倍关系。
2、服务器计算能力吃的消的话,建议增加双面slab的原子层数为单面模型的2倍,再计算表面能。

» 本帖已获得的红花(最新10朵)

7楼2016-08-24 10:07:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小竹子AA

金虫 (正式写手)

引用回帖:
7楼: Originally posted by zhao8919 at 2016-08-24 10:07:13
猜想原因:
1、单面slab、双面slab模型,优化后原子位置改变不大,所以能量改变也就不大,这样分母除以二之后出现2倍关系。
2、服务器计算能力吃的消的话,建议增加双面slab的原子层数为单面模型的2倍,再计算表面 ...

没看懂第一条原因,但觉得第二条原因是有可能的,确实我的单面slab模型中未固定的原子层数,是双面slab模型中的一端未固定原子层数的两倍。还有,我今天看了在模拟的constrain的过程中有固定分数坐标的,有固定笛卡尔坐标的,而普遍认为周期性结构的晶体在模拟过程中要固定分数坐标,而我固定的是笛卡尔坐标,会不会有影响啊?
谢谢您的回复,是不是您没有勾选应助?我这边不能发金币跟您
8楼2016-08-24 10:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小竹子AA

金虫 (正式写手)

送红花一朵
引用回帖:
7楼: Originally posted by zhao8919 at 2016-08-24 10:07:13
猜想原因:
1、单面slab、双面slab模型,优化后原子位置改变不大,所以能量改变也就不大,这样分母除以二之后出现2倍关系。
2、服务器计算能力吃的消的话,建议增加双面slab的原子层数为单面模型的2倍,再计算表面 ...

送你一朵花,谢谢您的回答
9楼2016-08-24 10:22:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hisap-nano

新虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
小竹子AA: 金币+5, ★★★很有帮助, 谢谢 2016-08-24 16:05:01
你应该比较2层和3层的结果,两种情况下都固定表面层以下的原子。第一层和第二层的构型差别很大。
10楼2016-08-24 11:18:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小竹子AA 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见