Conformer Search (GOAT)#
构象搜索(GOAT)

For flexible molecules, it is important to identify the low energy conformational space to provide the most accurate model of the investigated system. The global optimizer algorithm (GOAT) within ORCA can help to automatically find the global minimum structure and the ensemble around it. GOAT can be envoked by the simple input keyword GOAT:
对于柔性分子,识别低能构象空间至关重要,以提供所研究系统的最准确模型。ORCA 中的全局优化算法(GOAT)能够自动寻找全局最小结构及其周围的集合。通过简单的输入关键词 GOAT 即可调用 GOAT。

!GOAT XTB

*XYZFILE 0 1 inp.xyz

Note 注释

As a GOAT run involves a large number of geometry optimizations, it is recommended to parallelize the calculation, e.g. via the PAL4 keyword. Further, specifically for large systems, the usage of fast semi-empirical methods like GFN2-xTB speeds up the calculation. Nevertheless, GOAT can be used with any method in ORCA that is capable of optimizing geometries.
由于 GOAT 运行涉及大量几何优化,建议通过 PAL4 关键字等方法进行并行计算。此外,对于大型系统,使用 GFN2-xTB 等快速半经验方法可加速计算。尽管如此,GOAT 仍可与 ORCA 中任何具备优化几何能力的方法结合使用。

Example: Diclofenac# 示例:双氯芬酸 #

As an example, we will use to find the global minimum conformation of Diclofenac starting from the PubChem structure.
作为一个示例,我们将从 PubChem 结构出发,寻找双氯芬酸的全局最小构象。

../_images/diclofenac.png

Figure: Diclofenac 图:双氯芬酸#

The input looks like:
输入看起来像:

!GOAT XTB

*XYZFILE 0 1 diclofenac.xyz

After the successful GOAT run, the global minimum structure is stored in the basename.globalminimum.xyz file and the full structure ensemble in the basename.finalensemble.xyz file. In this case, GOAT found 17 unique conformers. The corresponding energies and weights are reported in the output file:
在成功完成 GOAT 运行后,全局最小结构存储在 basename.globalminimum.xyz 文件中,而完整结构集合则存储在 basename.finalensemble.xyz 文件中。在此情况下,GOAT 发现了 17 个独特的构象。相应的能量和权重在输出文件中报告:

                     ------------------------------
                                           _
                                          | |
                          __ _  ___   __ _| |_
                         / _` |/ _ \ / _` | __|
                        | (_| | (_) | (_| | |_
                         \__, |\___/ \__,_|\__|
                         __/ |
                        |___/

                      A Global Optimizer Algorithm
                     ------------------------------
[...]
                Global minimum found!
                Writing structure to diclofenac.globalminimum.xyz

                # Final ensemble info #
                Conformer     Energy     Degen.   % total   % cumul.
                              (kcal/mol)
                ------------------------------------------------------
                        0     0.000          1      75.54      75.54
                        1     0.976          1      14.56      90.10
                        2     1.991          1       2.62      92.72
                        3     2.028          1       2.46      95.18
                        4     2.413          1       1.29      96.47
                        5     2.572          1       0.98      97.45
                        6     2.621          1       0.91      98.36
                        7     2.799          1       0.67      99.03
                        8     3.339          1       0.27      99.30
                        9     3.420          1       0.24      99.53
                       10     3.781          1       0.13      99.66
                       11     3.827          1       0.12      99.78
                       12     4.027          1       0.08      99.86
                       13     4.373          1       0.05      99.91
                       14     4.614          1       0.03      99.94
                       15     4.624          1       0.03      99.97
                       16     4.752          1       0.02     100.00
                       17     5.951          1       0.00     100.00

                Conformers below 3 kcal/mol: 8
                Lowest energy conformer    : -55.771201 Eh
                Sconf at 298.15 K          :  1.83 cal/(molK)
                Gconf at 298.15 K          : -0.17 kcal/mol

                Writing final ensemble to diclofenac.finalensemble.xyz

We see, that the conformer space of Diclofenac is dominated by the two conformers 1 and 2, with their weights summing up to 90.1 %. The 10 lowest conformers look like:
我们观察到,双氯芬酸的构象空间主要由构象 1 和 2 主导,它们的权重总和达到 90.1%。前 10 个最低能量构象如下所示:

../_images/diclofenac.gif

Figure: The 10 lowest conformers of Diclofenac found by GOAT.
图:GOAT 发现的 Diclofenac 的 10 个最低能构象。
#

Structures# 结构

Diclofenac from PubChem 3033
来自 PubChem 3033 的双氯芬酸
30

Cl      1.95410      1.15000      -2.50780
Cl      1.13770     -1.63920       2.11360
 O     -3.26200     -2.92840      -1.06470
 O     -2.79060     -1.91080       0.90920
 N      0.26790     -0.20510      -0.39900
 C     -2.06400      0.51390      -0.37690
 C     -0.73130      0.71780      -0.01920
 C     -2.47610     -0.68300      -1.17030
 C      1.65710     -0.24820      -0.17950
 C     -3.03820      1.43500       0.00810
 C     -0.37280      1.84290       0.72340
 C     -2.67970      2.56000       0.75060
 C     -1.34700      2.76400       1.10830
 C      2.53530      0.34770      -1.09180
 C      2.17400     -0.88650       0.95340
 C     -2.84800     -1.87490      -0.31230
 C      3.91240      0.30580      -0.87390
 C      3.55110     -0.92850       1.17130
 C      4.42030     -0.33240       0.25760
 H     -1.70860     -0.97920      -1.89300
 H     -3.36140     -0.42660      -1.76760
 H     -0.08610     -1.11460      -0.67800
 H     -4.08120      1.28850      -0.26040
 H      0.65690      2.02780       1.01670
 H     -3.43820      3.27690       1.05110
 H     -1.06830      3.63990       1.68680
 H      4.60370      0.76540      -1.57580
 H      3.96350     -1.42150       2.04800
 H      5.49250     -0.36510       0.42740
 H     -3.50250     -3.70110      -0.51020