请教如何用orca来重复这一段使用dmol3的计算过程

牧生

同事请我帮他用主流的程序复现一下该段文字描述的计算过程，并计算出ΔE1 。作为新手，我虽然看得出来这一段写的极为辣鸡，但我主要用gmx，偶尔用用orca，我也并不会用orca来复现和计算这个ΔE1。  同事描述ΔE1为Reaction heat，N=2，且数值为–1.81 ev。

All spin-polarized DFT calculations are performed using the DMol3 code . All calculations are performed using the generalized gradient approximation of the B3LYP functional with Grimme method for DFT-D correction, the solvation model is COSMO using water as solvent (dielectric constant: 78.54), and the effective core potentials (ECPs) with double-numeric quality basis set with polarization functions (basis set is DNP). The Fermi coating method of 0.005 Ha is used to improve the computational performance. The convergence limits for energy change, maximum force, and maximum displacement are 2×10−5 Ha, 0.004 HaÅ−1, and 0.005 Å, respectively.

sobereva

GaussView画出反应物和产物结构，然后按下文创建几何优化输入文件（建议B97-3c级别），然后对得到的结构按下文创建wB97M-V/def2-TZVP级别的输入文件算单点，最后求差就完了。注意开启溶剂模型，优化期间建议用CPCM，单点期间建议用SMD
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490（http://bbs.keinsci.com/thread-13560-1-1.html）
（ORCA没法严格重复Dmol3的级别，如今的ORCA 5不再明确支持过时的COSMO，而且Dmol3用的数值原子轨道基组在ORCA里用不了。上面我说的计算级别远远比此文里的好）

我不知道C6F13是个什么东西。注意体系如果是高度柔性的，是讲究人的话应当用molclus（http://www.keinsci.com/research/molclus.html）先做构象搜索，至少产物这边理应做构象搜索。

在使用B97-3c优化之前，可以先让Molclus调用xtb在GFN-xTB级别下做构象的初筛。如果特别不讲究、懒得做构象搜索，就直接用xtb在GFN-xTB下对体系做个预优化，然后再B97-3c优化。