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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助成功】castep Si高压相变分析
411/1返回列表
查看: 4180  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 085600材料与化工 326 求调剂 +5 热爱生活ing 2026-03-09 5/250 2026-03-14 02:39 by JourneyLucky
[考研] 化学工程321分求调剂（南京工业，浙江工业） +3 大米饭！ 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭