计算化学公社

标题: 求助orca和gaussian联用算ts自定义溶剂部分语句 [打印本页]

作者
Author: thunder55 时间: 2023-4-2 11:02
提示: 该帖被管理员或版主屏蔽

作者
Author: sobereva 时间: 2023-4-2 14:40
光添加一行cpcmset一点用处也没有，这又不会写入到自动产生的ORCA输入文件里。一行行阅读orca.sh，搞清楚脚本产生ORCA输入文件的原理，在恰当位置插入溶剂模型设置

作者
Author: thunder55 时间: 2023-4-2 16:32
谢谢抽空社长回复！orca.sh文件的内容如下：

#Set the number of parallel cores and maximum memory utilized by each core (MB)
nprocs=20
maxcore=1000
#Set calculation level
level="r2SCAN-3c"
#Set parameters for other aspects
numset="tightscf noautostart"
cpcmset="epsilon=32.17 refrac=1.47"
#ORCA executable path
orcapath="..."

read atoms derivs charge spin < $2

if [ $derivs == "2" ] ; then
task="engrad freq"
elif [ $derivs == "1" ] ; then
task="engrad"
fi

#Create ORCA input file
echo "Generating mol.inp"
cat > mol.inp <<EOF
! $level $numset $task BOHRS nopop
%pal nprocs $nprocs end
%maxcore $maxcore
* xyz $charge $spin
$(sed -n 2,$(($atoms+1))p < $2 | cut -c 1-72)
*
EOF

#Run ORCA
echo "Running ORCA..."
$orcapath mol.inp > mol.out
echo "ORCA running finished!"

echo "Extracting data from ORCA outputs via extderi"
./extorca $3 $atoms $derivs

mv mol.gbw tmp.gbw -f
rm mol.* mol_* -f
mv tmp.gbw mol.gbw -f

orca创建输入文件应该是这部分吧：
#Create ORCA input file
echo "Generating mol.inp"
cat > mol.inp <<EOF
! $level $numset $task BOHRS nopop
%pal nprocs $nprocs end
%maxcore $maxcore
* xyz $charge $spin
$(sed -n 2,$(($atoms+1))p < $2 | cut -c 1-72)
*
EOF

是把! $level $numset $task BOHRS nopop这个语句添加为这样吗？! $level $cpcmset $numset $task BOHRS nopop
或者该在哪个恰当位置加入溶剂模型设置呢，怎么添加设置？请社长指导一下，谢谢！

作者
Author: sobereva 时间: 2023-4-3 10:01

thunder55 发表于 2023-4-2 16:32
谢谢抽空社长回复！orca.sh文件的内容如下：

#Set the number of parallel cores and maximum memory ut ...

搞懂ORCA输入文件写法，再了解点shell编程的基础知识，自然就知道了。还不清楚就实际测试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)