计算化学公社

标题: 求助:COLUMBUS的reordermo报错问题 [打印本页]

作者
Author:
jmliu    时间: 2023-2-11 12:11
标题: 求助:COLUMBUS的reordermo报错问题
新手最近正在学习Columbus做计算,目前已经成功计算MCSCF的单点能。但是体系需要进行换轨道。

根据tutorial,应该是使用自带程序reordermo。但是由于手册中的案例无需换轨道便没有给出具体案例。Usually, the active orbitals obtained in the first MCSCF calculation do not correspond to the active space that we have in mind.
When that happens, it is necessary to reorder the molecular orbital file and run the MCSCF calculation again.
COLUMBUS distribution contains an auxiliary program called reordermo.x to do this task.
This program may be accessed through the 7) Utilities option in colinp.




我计算的体系是C1对称性,74个电子,在cas(10,8)下进行计算。需要对32,33轨道进行交换。
进入程序Utilities的reordermo后提示:
Enter MOs to reorder as triples
   <DRT> <initial position> <final position>
Please enter a complete permutation for every DRT
根据提示输入了 1 32 33

出现报错:Error occurred in reordermo.x  errno=256, check runc.error

输出文件reordermols报错信息如下:
===== Program reordermo =====
------------------------------------------------------------------------
&input
mocoef='MOCOEFS/mocoef_mc.sp'
reordered='MOCOEFS/mocoef.new'
newmo(32,1)=33,
&end
------------------------------------------------------------------------
implicit mapping error newmo(33, 1)= 0
bummer (warning):incomplete permutation for drt # 1
mapping error:
               the new mo number on the right hand side appears
               several times for one DRT.
implicit mapping error:
               there was no complete permutation defined, hence
               collision with the implicit mapping, i.e. newmo(j,i)=j
bummer (fatal):input errors!2



新手对COLUMBUS的使用还很陌生,恳请大佬多多指教,谢谢!!




作者
Author:
beefly    时间: 2023-2-11 12:50
和molpro、molcas的mcscf比起来,columbus的mcscf没有任何优势,不仅速度慢一两个数量级,而且态平均的个数不能太多
作者
Author:
jmliu    时间: 2023-2-11 13:45
beefly 发表于 2023-2-11 12:50
和molpro、molcas的mcscf比起来,columbus的mcscf没有任何优势,不仅速度慢一两个数量级,而且态平均的个数 ...

是我用NEWTON-X做Trajectory surface hopping的动力学模拟。这个算cas和mrci用columbus比较方便。
作者
Author:
kkkjjjlll    时间: 2025-11-14 09:17
columbus调换轨道需要重复两次,比如调换调换轨道20和21,需要输入 1 20 21,以及 1 21 20,两个轨道顺序需要调换一下再输入一次




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3