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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】压力下几何结构的变化问题。
111/2返回列表12下一页
查看: 1075  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

鲤鱼飞

铁虫 (正式写手)

[交流] 【求助】压力下几何结构的变化问题。

如题。
我构建金属BCC单胞转原胞后，在不同压力下求P-V的关系。但直到压力达到5000GPa，a从2.62一直压缩到了1.56了，但夹角仍然为109.47°，这是由于BCC的对称性决定的角度么？国外报道在150GPa会发生BCC- 菱方的转变，也就是角度变化了。
但我在constraint下面，没办法在BCC的情况下去掉夹角的固定呐。如下图。


当我将BCC转为P1后，倒是可以都不固定了。但在初始优化时夹角都有变化了...加压力到150Gpa，夹角也有变化。我就不知道是压力引起的，还是自身优化引起的。

请问各位大虫虫们，在做P-V时如何设置的啊？？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
CAEP
1楼 2010-09-14 09:36:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)
  不要沉了。
回复此楼
CAEP
2楼2010-09-14 19:10:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
★
ym23(金币+1):欢迎常来交流 2010-09-14 22:30:46
鲤鱼飞(金币+1): 2010-09-14 23:09:35
记得以前读的文献中提到相变是通过比较两相能量随压力变化来确定，交叉处为相变，能量低为稳定新象。 可能是结构优化时只在原对称性附近优化， 而不是所有空间群搜索， 不容易发现体积部连续的一级相变。个人意见，仅供参考。
赞一下(1人) 回复此楼
3楼2010-09-14 21:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)
谢谢了。
回复此楼
CAEP
4楼2010-09-14 23:09:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
aylayl08(金币+2):谢谢侯博士回复 2010-09-15 15:51:41
鲤鱼飞(金币+1): 2010-09-15 22:06:43
这种加压的做法，是不行的。
很简单的想想： 两个不同的相，如同两个相的原胞里面的原子数可能不同，你再怎么通过这种方式的加压也看到相变。
赞一下(1人) 回复此楼
5楼2010-09-15 08:13:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★
aylayl08(金币+1):谢谢 2010-09-15 15:51:54
鲤鱼飞(金币+1): 2010-09-15 22:06:39
引用回帖:
Originally posted by babaleo at 2010-09-14 21:43:06:
记得以前读的文献中提到相变是通过比较两相能量随压力变化来确定，交叉处为相变，能量低为稳定新象。 可能是结构优化时只在原对称性附近优化， 而不是所有空间群搜索， 不容易发现体积部连续的一级相变。个人意见 ...

这一楼将的是一种方式。

另外一种方式，建立一个大的超原胞，做NPT的MD。

[ Last edited by valenhou001 on 2010-9-15 at 12:28 ]
赞一下(1人) 回复此楼
6楼2010-09-15 08:14:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟
引用回帖:
Originally posted by valenhou001 at 2010-09-15 08:14:26:


这是一种方式。

另外一种方式，建立一个大的超原胞，做NPT的MD。

侯老师您好，麻烦您讲授下这种方法的大致过程和INCAR的必要设置
赞一下 回复此楼
7楼2010-09-15 08:21:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★
aylayl08(金币+1):谢谢 2010-09-15 15:52:04
鲤鱼飞(金币+1): 2010-09-15 22:06:30
VASP里面没有实现NPT的分子动力学方法。实现不了我说的第二种方式来研究固体在高压下的相变。采用像PWSCF, CPMD之类的软件。
赞一下(1人) 回复此楼
8楼2010-09-15 12:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)
引用回帖:
Originally posted by valenhou001 at 2010-09-15 12:29:22:
VASP里面没有实现NPT的分子动力学方法。实现不了我说的第二种方式来研究固体在高压下的相变。采用像PWSCF, CPMD之类的软件。

那您的意思是在MS里面用MD的方法做，选择NVT吧，是一种体积守恒和温度不变的相变？ 如果是压力致的相变，我想不能在NPT等温等压下做吧？ 您的意思是否是先设置压力为150GPa或200GPa等可能的相变点，然后做NPT？？
赞一下 回复此楼
CAEP
9楼2010-09-15 22:11:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
zzy870720z(金币+2):谢谢侯老师指点 2010-09-16 15:48:13
鲤鱼飞(金币+1): 2010-09-16 18:08:24
晶体结构发生相变时，体积也必然会变的。你看MS里面的castep是否实现了NPT的系综的模拟。

建议你先找找相关资料读读，对NVT和NPT系综的模拟有点了解。

你可以参考这篇文章是怎么做的：此文对上面提到的两种方法都用到了来研究ZnSe的高压相变情况。
2009 J. Phys.: Condens. Matter 21 125403
The structural phase transition of ZnSe under hydrostatic and nonhydrostatic compressions: an ab initio molecular dynamics study
http://iopscience.iop.org/0953-8984/21/12/125403
赞一下(1人) 回复此楼
10楼2010-09-16 09:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 鲤鱼飞 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭