24小时热门版块排行榜    

查看: 1303  |  回复: 10

鲤鱼飞

铁虫 (正式写手)

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

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


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

请问各位大虫虫们,在做P-V时如何设置的啊??
回复此楼
CAEP
已阅   回复此楼   关注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
记得以前读的文献中提到相变是通过比较两相能量随压力变化来确定,交叉处为相变,能量低为稳定新象。 可能是结构优化时只在原对称性附近优化, 而不是所有空间群搜索, 不容易发现体积部连续的一级相变。个人意见,仅供参考。
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
这种加压的做法,是不行的。
很简单的想想: 两个不同的相,如同两个相的原胞里面的原子数可能不同,你再怎么通过这种方式的加压也看到相变。
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 ]
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之类的软件。
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
10楼2010-09-16 09:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 鲤鱼飞 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题,想知道的进来 +3 且听虎啸 2026-08-07 5/250 2026-08-07 16:50 by wdc28
[基金申请] 关于filecode +3 布布和一二 2026-08-07 4/200 2026-08-07 16:34 by 医学老男孩
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 基金中了 +11 laoda193707 2026-08-06 11/550 2026-08-07 10:54 by TT和NN
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见