| 查看: 1029 | 回复: 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 |
回复此楼» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有160人回复
-
祈福青基
已经有34人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
牛来！米来！面来！
已经有8人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
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
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
5楼2012-11-10 14:50:01
