| 查看: 998 | 回复: 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 |
回复此楼» 猜你喜欢
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
分子动力学模拟合作需求
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有211人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有1人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
fmklz
木虫 (正式写手)
搞电器工业化的化学打工仔
- 应助: 2 (幼儿园)
- 金币: 2874.3
- 散金: 50
- 红花: 4
- 帖子: 352
- 在线: 383.7小时
- 虫号: 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.7小时
- 虫号: 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
投票:
30