| 查看: 824 | 回复: 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 |
回复此楼» 猜你喜欢
PbS量子点如何确定浓度
已经有11人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有223人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
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