| 查看: 1770 | 回复: 4 | ||||
dongwen176木虫 (小有名气)
|
[求助]
MOLCAS中采用CASPT2计算CO2二价阳离子1a2态不收敛怎么办
|
|
各位大叔，我最近在计算CO2二价阳离子相关的态，拉势能曲线时，1a2态的第二个根的曲线十分靠近第一个根，固定键长大于2.7angstrom的点都不收敛，原因就是2.7以后两个根的能量十分的接近；对于这类情况，不知道该怎么修改输入文件来解决不收敛的问题呢？ 以下是输出文件片断： Final state energy(ies): ------------------------ root number 1 E = -186.50735064 a.u. root number 2 E = -186.50421657 a.u. |
回复此楼» 猜你喜欢
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有271人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 海基科技2012年工艺过程计算模拟软件gPROMS技术培训
已经有5人回复
- 妈妈得了间质性肺炎，这两天持续发烧，怎么办？？急！
已经有14人回复
- 怎么样禁止绕开队列提交任务？
已经有21人回复
- VASP计算论文一篇（含response，支持信息和正文修改痕迹）(校稿已经上传)
已经有48人回复
- 三维凝固模拟的时候，怎么用matlab画成三维体啊
已经有14人回复
- MPI的安装(管理员和普通用户均可)
已经有34人回复
- 难道每次计算前都咬进行收敛性测试吗？
已经有14人回复
- vasp4.6收敛的结构在vasp5.2 计算不收敛
已经有7人回复
- 纯理论计算投稿的问题
已经有21人回复
- 收敛精度
已经有7人回复
- 关于vasp收敛和优化
已经有24人回复
- Ni2+和Co2+是否也像Fe2+一样易被氧化？
已经有6人回复
- 求助DFT计算中CuCl2在γ-Al2O3 表面的吸附结构？
已经有5人回复
- 【讨论】dlpoly能量不收敛
已经有3人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13358
- 散金: 27307
- 红花: 302
- 帖子: 4684
- 在线: 929小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★
dongwen176(金币+20): 解释的很详细 2011-10-09 20:31:42
小红豆(金币+5): 谢谢 2011-10-09 21:34:06
dongwen176(金币+20): 解释的很详细 2011-10-09 20:31:42
小红豆(金币+5): 谢谢 2011-10-09 21:34:06
|
避免根的跳跃，需要做态平均计算。 首先选几个典型结构（键长：大、中、小；键角：大、中、小），在每个结构i，分别作N_i个根的casscf计算，N_i要取得稍大一些，可能要反复多试几次。数一下casscf级别的激发能（估算的绝热激发能；不要看垂直激发能）有多少个位于研究范围，假设是n_i。所以N_i一定要大于等于n_i，由于N_i=n_i的情况一般无法事先知道，所以至少要测试到N_i=n_i+1。 在正式的casscf + mrci或caspt2计算中，进行平均的电子态个数应当是n_i的最大值。 对于有对称性的体系，若研究全部低电子态，以上测试需要对每个不可约表示分别做。 |
2楼2011-10-09 13:58:16
dongwen176
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1831.3
- 红花: 1
- 帖子: 179
- 在线: 195.6小时
- 虫号: 91819
- 注册: 2005-09-07
- 性别: GG
- 专业: 理论和计算化学
3楼2011-10-09 20:10:54
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13358
- 散金: 27307
- 红花: 302
- 帖子: 4684
- 在线: 929小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
|
对不起，上面用“绝热激发能”的说法不合适。应该是激发态在典型结构的单点能，相对于基态平衡结构的能量。只需要casscf级别的一系列单点计算，不需要对激发态做结构优化。 典型结构就是有代表性的，比如CO2++，如果只关心中心对称的线形结构，就改变键长。可以选在（某个特定激发态的，或者基态的）平衡键长Re附近，小键长可以选在2/3*Re埃附近，大键长可以选在3/2*Re至2*Re之间。如果研究势能面和离解过程，还要加上一个更大的大键长（如8.0埃）。 如果只关心某个激发态，可以把它的权重定为1（如果是简并态，各个分量的权重也都是1），不感兴趣的态，可以把权重都设为一个小值，通常是0.1-0.5之间，只要不影响收敛，不造成根跳跃，当然是越小越好。 如果同时研究多个电子态，那么所有的权重一般都是1（但3d金属、Rydberg态的计算是例外）。 此外还有一种动态优化权重的方法，用得不多。 |
4楼2011-10-09 23:58:30
dongwen176
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1831.3
- 红花: 1
- 帖子: 179
- 在线: 195.6小时
- 虫号: 91819
- 注册: 2005-09-07
- 性别: GG
- 专业: 理论和计算化学
5楼2011-11-01 17:18:23
30