| 查看: 825 | 回复: 4 | ||
[求助]
molpro 的SAPT 中出错
|
|
这是输入文件 basis=aug-cc-pvtz ca=2101.2 cb=2102.2 dist=[2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.1,3.2,3.3] ang do i=1,#dist r=dist(i) hf edm=energy dummy,Ne {hf; save,$ca} ema=energy sapt;monomerA dummy,Ar {hf; start,atdens; save,$cb} emb=energy sapt;monomerB sapt,SAPT_LEVEL=3;intermol,ca=$ca,cb=$cb eint_hf=(edm-ema-emb)*1000. mH delta_hf=eint_hf-e1pol-e1ex-e2ind-e2exind elst(i)=E1pol; exch(i)=E1ex ind(i)=E2ind; exind(i)=E2exind disp(i)=E2disp; exdisp(i)=E2exdisp etot(i)=E12tot;eint_sapt(i)=eint_hf+e2disp+e2exdisp data,truncate,$ca enddo {table,dist,elst,exch,ind,exind,disp,exdisp,etot,eint_sapt ftyp,d,d,d,d,d,d,d,d,d,d plot} 错误信息 i=1 算完后，算i=2 时，出现 ?DUMP RECORD 2101.2 CORRESPONDS TO DIFFERENT BASIS OR MOLECULE |
回复此楼» 猜你喜欢
国家级人才团体课题组招收2026届博士
已经有1人回复
-
气相色谱点不着火
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有231人回复
-
化工安全有关的复习提纲
已经有0人回复
-
国家级人才课题组招收2026年入学博士
已经有2人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
散金
已经有35人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
已经有13人回复
fmklz
木虫 (正式写手)
搞电器工业化的化学打工仔
- 应助: 2 (幼儿园)
- 金币: 2874.3
- 散金: 50
- 红花: 4
- 帖子: 352
- 在线: 383.4小时
- 虫号: 770620
- 注册: 2009-05-14
- 性别: GG
- 专业: 能源化工
2楼2012-11-07 11:20:01
3楼2012-11-07 19:25:59
fmklz
木虫 (正式写手)
搞电器工业化的化学打工仔
- 应助: 2 (幼儿园)
- 金币: 2874.3
- 散金: 50
- 红花: 4
- 帖子: 352
- 在线: 383.4小时
- 虫号: 770620
- 注册: 2009-05-14
- 性别: GG
- 专业: 能源化工
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
qishiyy: 金币+20, ★★★★★最佳答案 2012-11-10 14:48:11
qishiyy: 金币+20, ★★★★★最佳答案 2012-11-10 14:48:11
|
我没有完全明白你的意思。简单说两句吧 1.你要是想达到这个目的呢？-----那只设R为变量即可，其他为常数，就可以了。 2.个人觉得实际情况，R较远时，构型变化不大，但较近时，构型肯定有些变化，在固定r值时，优化其他参数，来观察其他是否变化，若变化，便说明单体构型要变。 对于2情况，你要想保持hf波涵连续性，用个for循环。假设每一个R对应其他参数（固定R优化后）有r2，r3，alpha1，alpha2等等。 R=[,,,,,,,,] distance1=[,,,,,,,] distance2=[,,,,,,,,] angle1=[,,,,,,,,] angle2=[,,,,,,,] !这下面几个数组与相应的R保持一一对应 i=0; do ith=1,#R i=i+1 r1(i)=R(ith) r2(i)=distance1(ith) r3(i)=distance2(ith) alpha1(i)=angle1(ith) alpha2(i)=angle2(ith) ....... enddo 这样就把不同r的移动写入一个文件，再用start，save波涵，来继承上一个构型的波涵。 |
4楼2012-11-07 21:41:52
5楼2012-11-10 14:50:01