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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 【求助】castep中怎样设置温度
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 738  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

doudouding

木虫 (小有名气)

[交流] 【求助】castep中怎样设置温度

想计算体系在不同温度下的弹性常数，没找到在哪里能设置温度，请高手指点

[ Last edited by freshgirl on 2009-6-25 at 15:36 ]
回复此楼

» 猜你喜欢
1楼 2009-01-05 23:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

心宁

金虫 (正式写手)

哪里能设置温度?

★ ★
fegg7502(金币+2,VIP+0):thank you very much!
请问大家，有谁在castep模块，改变温度进行计算的，怎么改变温度呢？
Supercondcor(274056850) 21:17:48
castep不能改变温度进行计算
桔迷(183331864) 21:19:11
当然可以 castep 和 dmol的分子动力学
桔迷(183331864) 21:19:41
不过做好死机的准备 计算量相当大
sixhair (30473721) 21:20:12
  就看你服务器强不强 不过不一定要计算玩的

桔迷(183331864) 21:21:01
嗯
娃娃(35861885) 21:21:27
那怎么设置啊，请教一下
桔迷(183331864) 21:23:31
不知道　你看看ｃａｓｔｅｐ　的ｓｅｔｕｐ选项中有个ｄｙｎａｍｉｃｓ
桔迷(183331864) 21:23:49
算恒定温度用NVT系综
娃娃(35861885) 21:26:09
sixhair ：知道怎么改变温度吗
娃娃(35861885) 21:26:38
桔迷：如果要改变温度要选什么系综a
桔迷(183331864) 21:27:12
ＮＶＴ
桔迷(183331864) 21:27:32
自己改变呗
娃娃(35861885) 21:28:06
只要选定这个ＮＶＴ系综了，温度就直接改变就行了吗？
桔迷(183331864) 21:28:34
自己改温度　
桔迷(183331864) 21:28:56
改完了算啊　体系结构稳定就行了
桔迷(183331864) 21:29:34
要做退火的话那要自己编程了。。我也不太懂
娃娃(35861885) 21:32:10
算体系的话，是不是做好先优化一下结构，然后，算性质的时候set up 里面Task选定Dynamics，properties里面算则要计算的性质就行，这样子对吗？
桔迷(183331864) 21:33:20
还有步长　分子动力学很麻烦的　是动态过程你要等体系结构稳定为止

娃娃(35861885) 21:34:47
可以选择迭代的次数少一些吗？
桔迷(183331864) 21:35:54
这不是迭代！！这是真实的运动轨迹计算，你去看一下分子动力学的书吧
娃娃(35861885) 21:36:23
奥
桔迷(183331864) 21:37:22
计算精度的设定需要一定的经验，这个我帮不了你　当然高精度应该没问题　只是耗时多些
娃娃(35861885) 21:37:45
奥，好的，谢谢啊
赞一下(10人) 回复此楼
志存高远须积厚流广－坚持不懈定斗转星移！
2楼2009-01-05 23:19:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

doudouding

木虫 (小有名气)
明白了，就是说castep中没法用第一性原理的方法计算不同温度下体系的性质。

再问一个问题，实验上在不同温度下测量得到了材料的弹性常数，在用castep做计算时计算得到的结果应该和什么温度下的弹性常数做比较呢？也就说，castep程序中默认的温度是多少？
赞一下 回复此楼
3楼2009-01-06 08:18:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sss3012

银虫 (正式写手)
★
fegg7502(金币+1,VIP+0):thank you very much!
第一性原理计算得到的都是0K下的结果
赞一下(10人) 回复此楼
4楼2009-01-06 22:53:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuner

铁虫 (正式写手)
★
fegg7502(金币+1,VIP+0):thank you very much!
0K,这是第一原理的本质决定的哈
赞一下(10人) 回复此楼
5楼2009-01-06 22:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

心宁

金虫 (正式写手)

希望能得到有这方面计算经验的专家指点

★ ★ ★
fegg7502(金币+3,VIP+0):thank you very much!
4L 的结论我总是听到 不知道是不是真的 ？

Thermodynamic calculations
In principle, if the molecular geometry, energy, and vibrational frequencies are known, the thermodynamics properties (enthalpy, entropy, free energy) of the molecule can be calculated (Stull and Prophet, 1971). The results can be displayed using the Thermodynamic properties dialog, or viewed in the VAMP output file.

There is, however, a fundamental inconsistency inherent in the parameterization of semiempirical methods to reproduce heats of formation. The energy calculated by VAMP is the internal energy of a hypothetical motionless (Born-Oppenheimer) state. To relate this energy to heats of formation at 298 K, an atom-based scheme is used, assuming that the energy difference between the Born-Oppenheimer state and the molecule at 298 K can be treated in an additive fashion. This means that the zero-point energy and the energy required to warm the molecule to 298 K are assumed to be identical for isomers. This is clearly not the case and can lead to errors of up to 5 kcal mol-1.

A further problem arises in the calculation of the thermodynamic properties. As the methods are parameterized to reproduce the heat of formation at 298 K, this temperature must be used as the reference point for thermodynamics calculations. The calculation, however, gives an energy for a motionless molecule. This leads to a basic inconsistency. Thermodynamics calculations yield entropy, enthalpy, and free energy changes between the Born-Oppenheimer state and the temperature in question. The heat of formation scale must, however, be adjusted so that the value at 298 K is the same as that given in the normal calculation.

The result is that, in effect, we have to work with a fictitious Born-Oppenheimer energy in thermodynamics calculations. The best procedure is simply to ignore the "heat of formation" and to calculate thermodynamic quantities for reactions based on the Born-Oppenheimer energy and the calculated entropy and enthalpy changes, as for ab initio calculations.

Using the thermodynamic data
Once the enthalpy, entropy, and free energy of the molecule are calculated as a function of temperature, the results can be used to compute the heat of formation at any temperature. The heat of formation at 298 K is written to the VAMP output file and can be adjusted for different temperatures as follows:


Display the results using the Thermodynamic properties dialog or view them in the VAMP output file.
For the temperature of interest, say 100 K, compute the temperature-dependent contribution, enthalpy minus temperature times entropy (E - TS).
Compute E - TS at 298 K. If you did not compute the thermodynamics at exactly 298 K, pick the temperature closest to it (e.g., 300 K) or estimate it by interpolation.
Compute the difference in E - TS at the two temperatures and adjust the heat of formation at 298 K by that amount. For example, say that E - TS at 100 K is 0.5 kcal mol-1 and the value at 298 K is 2.2 kcal mol-1. You would compute the heat of formation at 100 K by subtracting 1.7 kcal mol-1 (2.2 - 0.5) from the heat of formation at 298 K given in the VAMP output file.

Further information
赞一下(10人) 回复此楼
志存高远须积厚流广－坚持不懈定斗转星移！
6楼2009-01-07 09:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

veryman

木虫 (正式写手)

$_$
★ ★
wuchenwf(金币+2,VIP+0):xiexie
引用回帖:
Originally posted by kuner at 2009-1-6 22:57:
0K,这是第一原理的本质决定的哈

第一原理目前看来只能算基态，就是0K下的状态，这是在KS方程的近似时的一个假设条件。
不够就算是基态，我们也可以得到很多信息了。
赞一下(10人) 回复此楼
7楼2009-01-07 16:31:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 doudouding 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +6 汪！？！ 2026-04-08 6/300 2026-04-08 21:38 by 朱云虎202
[考研] 320分人工智能调剂 +9 振—TZ 2026-04-03 10/500 2026-04-08 19:56 by 振—TZ
[考研] 生物学308分求调剂（一志愿华东师大）做过分子实验 +6 相信必会光芒万� 2026-04-07 7/350 2026-04-08 16:49 by tjzhao
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 一志愿211电子信息347求调剂 +3 554916 2026-04-03 3/150 2026-04-07 23:22 by 如若时光倒流
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 材料工程302分求调剂 +13 zyx上岸！ 2026-04-04 13/650 2026-04-07 11:14 by 诗与自由
[考研] 292求调剂 +4 lilllllxccc 2026-04-05 5/250 2026-04-07 09:29 by 纺大杨老师
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 求调剂，一志愿郑州大学材料与化工专硕，英二数二342分，求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 302分 085601求调剂推荐 +11 zyx上岸！ 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 309求调剂 +6 刘刘刘1231 2026-04-02 7/350 2026-04-04 13:41 by liucky
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 346求调剂 +5 郑诚乐 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
信息提示
关闭
请填处理意见
关闭