| 查看: 850 | 回复: 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 年 CSC/BBSRC phd 机器学习 催化剂设计方向
已经有0人回复
-
法国里昂大学招生CSC 资助博士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有61人回复
-
法国里昂大学接收CSC 资助访问学者
已经有0人回复
-
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）
已经有1人回复
-
帮导师招2026CSC博士（巴塞罗那自治大学UAB-CSC博士项目）
已经有25人回复
-
有没有LIMS客户的资源，求推荐
已经有0人回复
-
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）
已经有21人回复
-
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）
已经有20人回复
-
深圳理工大学招聘AI4Science与分子模拟方向科研人员
已经有0人回复
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