| 查看: 925 | 回复: 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 ] |
回复此楼» 猜你喜欢
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生
已经有0人回复
-
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有272人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
求助Cu2+1O的CIF文件(PDF: 05-0667)
已经有1人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
KAUST（阿卜杜拉国王科技大学）MXene 器件方向博士后招聘
已经有0人回复
-
沙特阿拉伯阿卜杜拉国王科技大学（KAUST）电池方向博士后招聘
已经有0人回复
-
福州大学新能源材料与工程研究院招收2026年入学博士
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样的审稿意见怎样回复？
已经有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