24小时热门版块排行榜    

查看: 1341  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 鲤鱼飞 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-30 01:58 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:49 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 01:14 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见