批量得到多环的NICS

123qwertybobo

如果认真看过上篇帖子http://bbs.keinsci.com/thread-22561-1-1.html，那么参考下就会写批量gjf的脚本（和本sh的目的无关）。这次相对会繁琐一点，但原理还是简单。前提是懂得在linux下使用Mulitiwfn求NICS。过程就是将"100\n24\n中心坐标\n环的序号\n相应的屏蔽张量”输入给Multiwfn，就能得到最终的NICS1zz。另外，上次可以同时做多个环的NICS，所以本shell必须能满足该功能。假设有2个lab文件，分别是Plane1.lab和Plane2.lab，那么会得到4个点的屏蔽张量NMR1.tmp2~NMR4.tmp2，两个质心坐标MS1.tmp3和MS2.tmp3，以及最终的NICS所用的Bq1-1.tmp4、Bq1-2.tmp4、Bq2-1.tmp4和Bq2-2.tmp4。之所以又去算质心是方便我验证数据。使用：因为我习惯是利用gjf得到同名文件夹，再把gjf批量复制进去再提交得到含Bq的log文件。所以Find_NICS.sh就是利用gjf依次进入同名文件夹处理log后，再依次把结果保存到起始目录的NICS.txt中。若有不同的习惯，可相应改之。另附删除tmp文件和lab文件可自行删除的语句，在Find_NICS.sh的末尾。在时间允许下，多问和多查是没有坏处的（虽然我的水平也很有限）。当然不允许的话，用就完了。

下回预告：批量求两点距离的sh。