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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助成功】castep Si高压相变分析
411/1返回列表
查看: 5304  |  回复: 40
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

susie8642

铜虫 (小有名气)

[交流] 【求助成功】castep Si高压相变分析

使用MS中的castep模块进行Si从c到beta的相变分析，通过加静水压，比较输出的final enthalpy进行相变分析，发现得到结果为20多GPa，与文献中的12.5GPa相差甚远，已经把截止能和k点选的相当高了，不知道问题出在哪里，请各位专家帮我解答一下吧！万分感激！是需要设constraint吗？还是输出并不是final enthalpy呢？困惑中！

[ Last edited by susie8642 on 2011-4-1 at 17:19 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-03-28 20:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xk6891

至尊木虫 (著名写手)
★ ★
susie8642(金币+1):谢谢参与
susie8642(金币+1): 2011-03-31 20:49:14
susie8642(金币+1): 非常感谢！ 2011-04-01 08:12:36
zzy870720z(金币+1): 谢谢建议 2011-04-01 08:25:35
个人感觉这个值还是较为正常，加静水压的算法准确度不高感觉，不如改变晶格常数，试一试
赞一下(3人) 回复此楼
13楼2011-03-31 17:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
引用回帖:
Originally posted by susie8642 at 2011-04-08 13:36:53:
那么请问我现在方法是不是已经正确了呢？我又把点取得密了一些，发现比较接近正确值了。

我实在不敢妄加指导，只是个人经验
赞一下(1人) 回复此楼
33楼2011-04-08 14:34:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

susie8642

铜虫 (小有名气)
自己顶一下，太困惑了。。。
回复此楼
2楼2011-03-29 08:20:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★ ★
susie8642(金币+1):谢谢参与
zzy870720z(金币+1): 谢谢交流 2011-03-30 08:48:37
susie8642(金币+1): 2011-03-30 18:55:23
引用回帖:
Originally posted by susie8642 at 2011-03-29 08:20:40:
自己顶一下，太困惑了。。。

看不见你的具体过程，可以把具体过程截图发上来。。。。
赞一下(2人) 回复此楼
3楼2011-03-30 02:05:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by zyj8119 at 2011-03-30 02:05:32:
看不见你的具体过程，可以把具体过程截图发上来。。。。

没有什么过程，就是给晶格不断施加静水压，然后得到enthalpy，然后结果差很远。。。求教啊！
赞一下 回复此楼
4楼2011-03-30 18:53:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★
susie8642(金币+1):谢谢参与
不断施加静水压，enthalpy是怎么得到的呢？
赞一下(1人) 回复此楼
5楼2011-03-30 21:03:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by bingmou at 2011-03-30 21:03:42:
不断施加静水压，enthalpy是怎么得到的呢？

是在.castep文件中得到的
赞一下 回复此楼
6楼2011-03-30 22:33:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by susie8642 at 2011-03-30 22:33:32:
是在.castep文件中得到的

不过不知道这样得到的对不对啊
赞一下 回复此楼
7楼2011-03-30 22:54:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seek269

铁虫 (正式写手)
★
susie8642(金币+1):谢谢参与
susie8642(金币+1): 2011-03-31 11:37:48
有可能你方法搞错了
赞一下(1人) 回复此楼
8楼2011-03-31 11:33:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
我也觉得很有可能是这样的，但是不知道正确的方法到底是什么呢，请指教啊，好迷茫啊
赞一下 回复此楼
9楼2011-03-31 11:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★ ★
susie8642(金币+1):谢谢参与
zzy870720z(金币+1): 谢谢指点 2011-03-31 17:20:00
susie8642(金币+1): 非常感谢 2011-03-31 17:26:41
你把Si的c和beta相在一系列压强下的enthalpy都算出来，作图，看看有没有焦点，如果有交点就说明有相变
赞一下(2人) 回复此楼
10楼2011-03-31 15:58:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
非常感谢，我就是这样做的，但是交点在20GPa，不知道是哪里出了问题呢？
赞一下 回复此楼
11楼2011-03-31 17:26:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by fengshiquan at 2011-03-31 15:58:11:
你把Si的c和beta相在一系列压强下的enthalpy都算出来，作图，看看有没有焦点，如果有交点就说明有相变

麻烦你再详细说一下，看看我哪里错了
赞一下 回复此楼
12楼2011-03-31 17:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
zzy870720z(金币+1): 谢谢指教 2011-04-01 08:25:52
引用回帖:
Originally posted by susie8642 at 2011-03-31 17:27:05:
麻烦你再详细说一下，看看我哪里错了

就是做两相的焓值图，看焓值图有没有交点，如果有，交点处对应的压强即为相变压强
赞一下(1人) 回复此楼
14楼2011-03-31 19:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by fengshiquan at 2011-03-31 19:14:13:
就是做两相的焓值图，看焓值图有没有交点，如果有，交点处对应的压强即为相变压强

非常感谢！这个焓值是不是就是。castep输出的那个enthalpy值呢？相变过程需要对晶胞进行什么设置吗？
赞一下 回复此楼
15楼2011-03-31 20:50:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
susie8642(金币+5): 2011-04-01 08:13:51
zzy870720z(金币+1): 谢谢补充 2011-04-01 08:26:04
就是castep输出的那个enthalpy值，相变过程不需要对晶胞进行什么特殊设置
赞一下(1人) 回复此楼
16楼2011-03-31 21:57:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by xk6891 at 2011-03-31 17:38:08:
个人感觉这个值还是较为正常，加静水压的算法准确度不高感觉，不如改变晶格常数，试一试

非常感谢！可以讲一讲怎样通过改变晶格常数进行测试的方法吗？
赞一下 回复此楼
17楼2011-04-01 08:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by fengshiquan at 2011-03-31 21:57:32:
就是castep输出的那个enthalpy值，相变过程不需要对晶胞进行什么特殊设置

非常感谢！那也就是我的测试方式没有什么问题了，那为什么我的结果会差那么远呢？我看还有用cell volume进行测试的，如图所示，但是在加静水压下，材料的volume，castep并没有给出，这个方式请问您知道吗？是怎样具体操作的呢？
赞一下 回复此楼
18楼2011-04-01 08:18:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
susie8642(金币+2): 2011-04-01 09:38:42
引用回帖:
Originally posted by susie8642 at 2011-04-01 08:18:31:
非常感谢！那也就是我的测试方式没有什么问题了，那为什么我的结果会差那么远呢？我看还有用cell volume进行测试的，如图所示，但是在加静水压下，材料的volume，castep并没有给出，这个方式请问您知道吗？是怎 ...

在加静水压下，材料的volume，castep应该会给出，看有没有相变可以看一下有没有体积突变，这是相变的特征，有突变说明可能是相变点，然后再进一步判断
赞一下 回复此楼
20楼2011-04-01 09:17:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by fengshiquan at 2011-04-01 09:17:43:
在加静水压下，材料的volume，castep应该会给出，看有没有相变可以看一下有没有体积突变，这是相变的特征，有突变说明可能是相变点，然后再进一步判断

非常感谢！由于加静水压后，在计算前castep会提示将结构转化为p1，这样在.castep文件中就没有volume这一项的输出了，我想请问一下这时从哪里可以得到volume值呢？
赞一下 回复此楼
21楼2011-04-01 09:40:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
gzqdyouxia(金币+1): 谢谢提示 2011-04-01 10:14:49
susie8642(金币+1): 2011-04-01 10:48:50
我在加静水压后，在.castep文件中就有volume这一项的输出了，不知道你的是什么情况
赞一下(1人) 回复此楼
22楼2011-04-01 10:03:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
★
gzqdyouxia(金币+1): 谢谢交流 2011-04-01 10:44:52
susie8642(金币+1): 2011-04-01 10:50:29
引用回帖:
Originally posted by susie8642 at 2011-04-01 08:13:22:
非常感谢！可以讲一讲怎样通过改变晶格常数进行测试的方法吗？

加压相当于改变晶格常数的大小，这样的话，可以减小晶格常数达到加压的目的，但是这个过程中压力的变化是如何的，就是减少0.1埃相当于加多少压，还需要想办法得到
赞一下(1人) 回复此楼
23楼2011-04-01 10:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by fengshiquan at 2011-04-01 10:03:51:
我在加静水压后，在.castep文件中就有volume这一项的输出了，不知道你的是什么情况

对c-Si结构是这样的,但是我还尝试了beta-Si结构，这个时候就没有volume输出了，不知道是什么原因，你可以也试一下
赞一下 回复此楼
24楼2011-04-01 10:50:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by xk6891 at 2011-04-01 10:34:12:
加压相当于改变晶格常数的大小，这样的话，可以减小晶格常数达到加压的目的，但是这个过程中压力的变化是如何的，就是减少0.1埃相当于加多少压，还需要想办法得到

非常感谢！这个方法就与改变volume是近似的吧？不知这种加压相变问题最经典的方法是哪个呢？
赞一下 回复此楼
25楼2011-04-01 10:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
susie8642(金币+1): 2011-04-01 11:33:10
susie8642(金币+1): 2011-04-01 11:35:38
引用回帖:
Originally posted by susie8642 at 2011-04-01 10:51:20:
非常感谢！这个方法就与改变volume是近似的吧？不知这种加压相变问题最经典的方法是哪个呢？

看到的和听到的结合周围人作的， 改变volume据说是最好的，看一下理论把，还是
赞一下 回复此楼
26楼2011-04-01 11:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by xk6891 at 2011-04-01 11:29:22:
看到的和听到的结合周围人作的， 改变volume据说是最好的，看一下理论把，还是

非常感谢！改变volume就是您之前提及的改变晶格常数的方式是吧？想确认一下，呵呵
赞一下 回复此楼
27楼2011-04-01 11:36:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
susie8642(金币+1): 3 2011-04-01 11:59:37
引用回帖:
Originally posted by susie8642 at 2011-04-01 11:36:22:
非常感谢！改变volume就是您之前提及的改变晶格常数的方式是吧？想确认一下，呵呵

恩，是
回复此楼
28楼2011-04-01 11:49:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by xk6891 at 2011-04-01 11:49:29:
恩，是

经过计算，还是不太对，再次请教！按照改变晶格常数方法，建立了volume与相变的关系，见图。发现在volume为29时发生相变。这时，在volume与pressure的关系图中去找，发现相变压强在2或30。感觉不太对，使用老方法用pressure和enthalpy的关系，发现相变点在21左右。实际相变点在12.5GPa左右，怎么都不对，很苦恼。求解答。。。







[ Last edited by susie8642 on 2011-4-7 at 16:17 ]
赞一下 回复此楼
29楼2011-04-07 16:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
引用回帖:
Originally posted by susie8642 at 2011-04-07 16:13:26:
经过计算，还是不太对，再次请教！按照改变晶格常数方法，建立了volume与相变的关系，见图。发现在volume为29时发生相变。这时，在volume与pressure的关系图中去找，发现相变压强在2或30。感觉不太对，使用老 ...

这个就不知道了，呵呵
赞一下 回复此楼
30楼2011-04-08 12:23:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by xk6891 at 2011-04-08 12:23:16:
这个就不知道了，呵呵

那么请问我现在方法是不是已经正确了呢？我又把点取得密了一些，发现比较接近正确值了。
赞一下 回复此楼
31楼2011-04-08 13:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
★
franch(金币+1): 谢谢回帖交流 2011-04-08 17:26:50
引用回帖:
Originally posted by susie8642 at 2011-04-08 13:36:53:
那么请问我现在方法是不是已经正确了呢？我又把点取得密了一些，发现比较接近正确值了。

“我又把点取得密了一些，发现比较接近正确值了”如果这样的话，是不是可以提高精度，进一步试一试
赞一下(1人) 回复此楼
32楼2011-04-08 14:34:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by xk6891 at 2011-04-08 14:34:55:
我实在不敢妄加指导，只是个人经验

请问静水压怎么加呢？谢谢！
赞一下(1人) 回复此楼
34楼2011-04-08 21:14:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by xk6891 at 2011-04-08 14:34:55:
我实在不敢妄加指导，只是个人经验

已经非常感谢您了！受益匪浅！只是不知道这样的方式是否能够应用于其他晶胞的相变中，而不是正好凑数凑上了，所以想问一下。
赞一下 回复此楼
35楼2011-04-09 08:13:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by frank_zhan at 2011-04-08 21:14:53:
请问静水压怎么加呢？谢谢！

使用castep几何优化中自带的施压选项进行的，不知是否正确？
赞一下 回复此楼
36楼2011-04-09 08:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xk6891

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by susie8642 at 2011-04-09 08:13:01:
已经非常感谢您了！受益匪浅！只是不知道这样的方式是否能够应用于其他晶胞的相变中，而不是正好凑数凑上了，所以想问一下。

没帮上什么忙，不谢，呵呵。多交流。
赞一下(1人) 回复此楼
37楼2011-04-09 09:31:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by susie8642 at 2011-04-09 08:13:41:
使用castep几何优化中自带的施压选项进行的，不知是否正确？

谢谢，我看见了。
再请问一下这位老兄知道怎么给计算体系加上电场吗？
赞一下(1人) 回复此楼
38楼2011-04-12 09:49:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

susie8642

铜虫 (小有名气)
引用回帖:
Originally posted by frank_zhan at 2011-04-12 09:49:33:
谢谢，我看见了。
再请问一下这位老兄知道怎么给计算体系加上电场吗？

这个还真不知道，不好意思哈
赞一下 回复此楼
39楼2011-04-12 11:26:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
加压之后是不是就几何优化就可以了？
赞一下 回复此楼
40楼2012-03-15 23:07:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
456913楼: Originally posted by susie8642 at 2011-04-08 13:36:53
那么请问我现在方法是不是已经正确了呢？我又把点取得密了一些，发现比较接近正确值了。

请问“我又把点取得密了一些，发现比较接近正确值了”，您这时是用MS自带的施压选项进行的加压，还是通过改变晶格常数加压？

谢谢！
赞一下 回复此楼
41楼2012-07-18 10:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
tangyanxmc19楼
2011-04-01 08:28   回复  
susie8642(金币+1):谢谢参与
祝福
相关版块跳转 我要订阅楼主 susie8642 的主题更新
411/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭