| 查看: 1087 | 回复: 5 | ||
[求助]
怎样外推到CBS极限
|
|
用CCSD(T)+AVDZ,AVTZ,AVQZ,AV5Z分别算了一个双原子体系100个点的单点能，怎样进行基组外推? 如果用mixed exponential/Gaussian function : E(x)=Ecbs+b*exp(-(x-1))+c*exp(-(x-1)^2) 进行拟合应该怎样操作？ 1）是将对应相同键长的avdz ，t，q，5z的四个数据单独拟合出b和c？ 这种方法我试过，每个组合拟合的b和c是不一样的，不会是将所有100个键长的点都一个个拟合得到的Ecbs吧？ 2）还是将所有400个点一起拟合？ 应该怎样操作。 [ Last edited by kent1022 on 2012-5-10 at 21:35 ] |
回复此楼» 猜你喜欢
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
关于合成MIL-125的求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有224人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样的审稿意见怎样回复？
已经有7人回复
- 怎样判断DEAE-纤维素柱上样超载？急求！！
已经有8人回复
- 12%乙酸溶液怎样配制
已经有4人回复
- 关于动力学极限电流密度
已经有3人回复
- 求助：使用B-methyl CBS 不对称还原酮为醇，CBS使用溶液的还是固体的好？
已经有4人回复
- 我收到的DS2019表和邀请函
已经有8人回复
- 大家认为英语的好差对论文的录用起到什么程度的影响
已经有35人回复
- 基金填写完无法打印，有遇到这种情况的么
已经有10人回复
- 【求助】怎样将羧基与二氯亚砜反应生成酰氯？(很急)
已经有41人回复
mchen10
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1508.2
- 散金: 41
- 红花: 6
- 帖子: 425
- 在线: 159.5小时
- 虫号: 695973
- 注册: 2009-02-04
- 性别: GG
- 专业: 理论和计算化学
2楼2012-05-11 00:52:15
3楼2012-05-11 08:57:56
mchen10
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1508.2
- 散金: 41
- 红花: 6
- 帖子: 425
- 在线: 159.5小时
- 虫号: 695973
- 注册: 2009-02-04
- 性别: GG
- 专业: 理论和计算化学
4楼2012-05-11 23:20:39
|
试了一下，确实可以用批处理搞定。 不过对于离解能De和平衡核间距Re应该怎样确定？ 看到过有文章把De也用上面的mixed exponential/Gaussian function公式拟合得到，不过没发现有文章对Re也是这样拟合得到。 另外，可以通过某一势能函数拟合各CBS单点能得到Re和De 。 我不清楚这两种情况的结果分别是怎样的，于是把这两种情况分别试了一下，得到的两组Re和De的值还是有差别的。 现在的问题是是（1）将Re通过外推得到，在拟合单点能时固定，从而得到De ？ 还是应该（2）将Re和De均通过CBS的单点能用某种势能函数拟合得到？ 比如，我之前分别用AVDZ,AVTZ,AVQZ,AV5Z扫描的单点能都是把Re固定在各自优化的Re结果上用M-S函数拟合出De |
5楼2012-05-13 21:15:31
mchen10
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1508.2
- 散金: 41
- 红花: 6
- 帖子: 425
- 在线: 159.5小时
- 虫号: 695973
- 注册: 2009-02-04
- 性别: GG
- 专业: 理论和计算化学
6楼2012-05-14 01:06:12
20