| 查看: 2801 | 回复: 11 | |||
[交流]
【求助】用Gaussian计算分子的熵时，不同分子的熵的计算基准是一样的吗？谢谢
|
| 请教大虾，用Gaussian计算分子的熵时，不同分子（如甲烷、乙烷等）的熵的计算基准是一样的吗？是不是默认0K时所有分子的熵均为零呢？谢谢啦 |
回复此楼» 猜你喜欢
深圳理工大学招聘AI4Science与分子模拟方向科研人员
已经有0人回复
-
全奖博士 英国利物浦大学 × 台湾清华大学 联合培养
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有77人回复
-
寻博导
已经有13人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有9人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问能用MS软件进行热膨胀系数、自由焓、熵等热力学参数的计算吗？
已经有20人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 怎么用gaussian计算计算双原子分子势能曲线？
已经有15人回复
- 关于熵变的计算，有人知道怎么计算大的体系的熵变吗？比如将近一千个残基的体系？
已经有11人回复
- Gaussian对分子构型优化加溶剂化与否对最后结果有影响吗？
已经有8人回复
- gaussian03计算分子化学硬度
已经有6人回复
- 【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算
已经有9人回复
- 【求助】绝对熵值的计算
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
鲍红丽课题组 研究生招生启事
+1/477
-
南京大学智能驱动与感知材料实验室诚招博士/硕士/科研助理/访问学生
+1/281
-
专业技术开发及第三方检测
+1/87
-
物理学 调剂
+1/86
-
浙江农林大学招收调剂
+1/85
-
材料/化学相关专业2026级学术/专业型硕士研究生
+1/82
-
厦门大学电池工况表征博士后招聘（固态电池 / 层状氧化物正极 / 原位-工况表征方向）
+1/78
-
青岛大学泰山学者课题组招2026年申请考核制博士
+1/34
-
中科院生态环境研究中心国重实验室招聘客座研究生1-2名
+1/33
-
0854电子信息调剂，集成电路，芯片方向，闽南师范大学光电芯片研发实验室
+1/33
-
中科院化学所 宋延林 课题组招聘合成化学方向博士后（开展打印合成化学方向研究）
+1/32
-
2026年工科硕士调剂-上海大学全国重点实验室团队-材料数据挖掘方向-研究生3-5人
+1/32
-
【沈阳农业大学机械专硕调剂】
+2/30
-
南京医科大学-膜蛋白结构与功能课题组-招收2026级博士研究生
+1/16
-
上海科技大学物质学院招收有机化学方向2026级工程博士
+1/8
-
发paper
+1/5
-
杭州北航国新院徐梁教授招聘博士（控制/机器视觉/机器人导航操作方向）
+1/2
-
东北农业大学（211）水利工程招生
+1/1
-
澳大利亚科学院院士Prof.Gordon Wallace, University of Wollongong 全奖博士机会
+1/1
-
2026年浙江清华长三角研究院智慧能源中心招聘博士后-电气方向
+1/1
2楼2010-12-15 22:50:10
|
若所有分子在绝对零度时熵为零的话，我计算了一下乙烷在在接近绝对零度时（0.001K，1.0 atm ）时的熵，却出现了负值，这是为什么呢? 部分输出结果如下，详细输出文件见附件http://d.namipan.com/d/be753ea92 ... 39c4feef9f628850100 ------------------- - Thermochemistry - ------------------- Temperature 0.001 Kelvin. Pressure 1.00000 Atm. Atom 1 has atomic number 6 and mass 12.00000 Atom 2 has atomic number 1 and mass 1.00783 Atom 3 has atomic number 1 and mass 1.00783 Atom 4 has atomic number 1 and mass 1.00783 Atom 5 has atomic number 6 and mass 12.00000 Atom 6 has atomic number 1 and mass 1.00783 Atom 7 has atomic number 1 and mass 1.00783 Atom 8 has atomic number 1 and mass 1.00783 Molecular mass: 30.04695 amu. Principal axes and moments of inertia in atomic units: 1 2 3 EIGENVALUES -- 22.50754 90.72976 90.73040 X 1.00000 0.00000 0.00000 Y 0.00000 1.00000 0.00000 Z 0.00000 0.00000 1.00000 This molecule is an asymmetric top. Rotational symmetry number 1. Rotational temperatures (Kelvin) 3.84821 0.95464 0.95463 Rotational constants (GHZ): 80.18384 19.89139 19.89125 Zero-point vibrational energy 197533.4 (Joules/Mol) 47.21161 (Kcal/Mol) Vibrational temperatures: 449.61 1197.14 1197.72 1453.03 1777.62 (Kelvin) 1778.20 2062.27 2092.39 2203.85 2204.29 2211.96 2212.65 4384.44 4386.18 4458.40 4458.62 4493.50 4493.68 Zero-point correction= 0.075236 (Hartree/Particle) Thermal correction to Energy= 0.075237 Thermal correction to Enthalpy= 0.075237 Thermal correction to Gibbs Free Energy= 0.075237 Sum of electronic and zero-point Energies= -79.755183 Sum of electronic and thermal Energies= -79.755183 Sum of electronic and thermal Enthalpies= -79.755183 Sum of electronic and thermal Free Energies= -79.755183 E (Thermal) CV S KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin Total 47.212 5.962 -44.209 Electronic 0.000 0.000 0.000 Translational 0.000 2.981 -26.490 Rotational 0.000 2.981 -17.719 Vibrational 47.212 0.000 0.000 |
3楼2010-12-15 23:14:08
4楼2010-12-16 09:58:40
5楼2010-12-16 13:04:10
|
Based on my understanding, computational programms have numerical error and limitations. Sometimes when a program takes a "BAD" variable (e.g. too closed to 0), it could not function normally. Gaussian might do the same, its entropy calculation might be good for most of cases (like 1K-9999K probabally), but it is not suitable for a very small T like 0.001, which won't affect its usage much.. I was talking to my self.. again.. |
6楼2010-12-16 14:37:45
7楼2010-12-16 17:01:31
★ ★ ★
freesea(金币+1): 2010-12-16 23:22:32
小红豆(金币+3):鼓励交流 2010-12-17 10:01:32
freesea(金币+1): 2010-12-16 23:22:32
小红豆(金币+3):鼓励交流 2010-12-17 10:01:32
|
【商榷】即使对理想气体，在温度很低时Gaussian的计算结果也是错误的。因为Gaussian对于平动熵和转动熵的计算中使用了配分函数，参见 http://www.gaussian.com/g_whitepap/thermo.htm 而有关的配分函数的表达式则采用了经典近似，这种近似在温度很低时并不成立。兹以平动为例，其经典近似条件为 (V/N)^(1/3) >> h/(3mkT)^(1/2) [参见：朱文浩，顾毓沁编著，《统计物理学基础》，清华大学出版社，192页] 式中V为体积，N为分子数。显然，当T趋于0K时，上述近似条件失效。因此有关的公式不能用于温度很低的情况，否则会出现许多奇怪的事。 按理说，在此种情况下Gaussian应该给出警告，但遗憾的是它并没有给出这种警告。因此，如果把Gaussian完全当做黑箱来用是很危险的。 |
8楼2010-12-16 21:08:21
9楼2010-12-16 23:39:22
10楼2010-12-17 08:37:53
11楼2010-12-17 08:43:19
12楼2010-12-17 09:25:11
