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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请教大家一个关于压力、温度模拟的问题
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1355  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

future_wl

木虫 (著名写手)

[交流] 请教大家一个关于压力、温度模拟的问题

如果一个晶体各个轴向上原子之间结合力的强弱不一样，那么施加温度后，会出现有的键伸缩容易，有的键伸缩困难；如果施加不同的压力会不会也出现类似的情况：键强的会变化的慢，键弱的变化的快。
因为软件很难模拟温度的影响，可否用压力的模拟来代替，来算出具体化学键的强弱。
希望大家多多指导！
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-09-25 23:44:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
future_wl: 金币+10, 感谢你细致的讲解！ 2012-09-26 08:53:44
WDD880227: 金币+1, 多谢耐心指教 2012-09-26 16:39:44
不知楼主打算怎样用加温或加压后的模拟来考察化学键的强度？难道要开振动频率？可是振动频率与键的强度之间并无直接的对应关系啊。要用化学键断裂极限温度，或者错位极限压力？

我的理解是，化学键的强度是用键能来衡量，即把断键看成反应，来看反应热和反应能垒。通常对于有机分子晶体，可以认为固体中的键能和气相中的单分子的键能基本相等，所以可以直接对气相单分子计算键断势能面，从而得到键断裂前后的键能。因为晶体堆积能比反应能小1-2个数量级，它只能使化学键在势能面中的平衡位置稍稍移动，基本还是在基态势能面的底部，除非极高温高压。实际上高温高压影响的，只是原子在势能面上的访问范围，而固体中的化学键的势能面基本还是和气相中的一致，顶多加个介电常数校正和固体堆积能对化学反应能的校正。

上面说的是有机分子晶体，对于周期固体，特别是金属，无机盐、配位化合物，与有机分子晶体中的化学键能的研究完全不同，需要用另外的量来衡量。什么量呢？

对于金属和金属氧化物固体，可以用空穴形成能来估算。当然，移走一个原子，形成一个空穴，可能涉及断裂多根配位键。因此这个量只能间接反映某个化学键的强弱。

表面原子的解离能。表面反应更可能对应实际的键解离反应。

还有更复杂的，假如实验上已经知道某种电子激发导致某个键断裂或某个反应发生，此时可以通过电子激发态计算来计算对应的激发能，从而真实的了解固体中的化学键的强度。不过，目前第一原理固体的激发态计算还不成熟，误差比较大。

不知道这些是不是楼主正在考虑的。一起讨论。
赞一下 回复此楼
3楼2012-09-26 08:37:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
3楼: Originally posted by ChemiAndy at 2012-09-26 08:37:34
不知楼主打算怎样用加温或加压后的模拟来考察化学键的强度？难道要开振动频率？可是振动频率与键的强度之间并无直接的对应关系啊。要用化学键断裂极限温度，或者错位极限压力？

我的理解是，化学键的强度是用键能 ...

你好，我是打算用晶格振动里面的力常数来描述化学键的强弱，不知道可不可行。
vasp可以计算力常数矩阵，我想的是：如果体系存在两种化学键，较强的那个对应的两个原子之间的力常数数值应该相对较大（对于体系基态性质的计算证明了这一点）。
现在的问题是：实验上观测到随温度的升高，较弱的化学键膨胀的越来越快，而较强的化学键膨胀的越来越慢。目前软件很难模拟温度的变化，我就想能不能通过施加压力来观察力常数的变化。
请指教!
赞一下 回复此楼
4楼2012-09-26 09:42:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
如果所比较的化学键是同一类型的，比如，都是Si-O键，应该可以；但是如果是两种不同类型的键，则力常数没有可比性，因为力常数还与质量有关。

加温看化学键强度类似于微扰法，即把温度T看作扰动，然后看键长的变化；如果将键长变化看作温度T的线性函数，则为
Δr=kΔT
然后比较k的大小。

现在你们要用P取代T，然后看：
Δr=kΔP
没有什么不可以，就是一个扰动而已。不过扰动要尽可能小，使响应保持在线性区，否则，这个展开可能需要到高级。
赞一下 回复此楼
5楼2012-09-27 15:20:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
不如看海2楼
2012-09-26 00:16   回复  
相关版块跳转 我要订阅楼主 future_wl 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学0703-一志愿211-338分求调剂 +7 vants 2026-04-05 7/350 2026-04-05 18:17 by cql1109
[考研] 0855求调剂材料 +10 红桃灼灼 2026-04-04 10/500 2026-04-05 18:08 by 猪会飞
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 一志愿华中农业大学0710（A）初试329分 求调剂 +4 一名26考研生 2026-04-04 4/200 2026-04-05 10:01 by barlinike
[考研] 081700，311，求调剂 +12 冬十三 2026-04-04 13/650 2026-04-05 08:43 by qlm5820
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 0835学硕299求调剂 08大类可接受 +5 useryy 2026-04-03 5/250 2026-04-04 20:07 by 蓝云思雨
[考研] 。 +5 雾与海 2026-04-02 6/300 2026-04-04 19:53 by 蓝云思雨
[考研] 考研调剂 +4 zybz冲冲冲 2026-04-03 6/300 2026-04-04 13:08 by zybz冲冲冲
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[考研] 266分，求材料相关专业调剂 +13 哇呼哼呼哼 2026-03-30 15/750 2026-04-03 15:24 by arrow8852
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 085600，320分求调剂 +6 大馋小子 2026-04-02 6/300 2026-04-02 21:54 by dongzh2009
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 085602化工求调剂（331分） +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
信息提示
关闭
请填处理意见
关闭