¿Es Orca rompiendo erróneamente bonos? ¿Por qué puede sólo optimizar la molécula con COPTO activado?

Question

Me gustaría pedir ayuda en el uso de Orca (v. 4.0.0) para la realización de una simple optimización.

En primer lugar tengo esta estructura (principalmente optimizado con MMFF94s campo de fuerza de Avogadro):

Por eso, he tratado de llevar una simple optimización mediante el uso de Orca y sólo como una prueba he utilizado la siguiente entrada:

! PAL8 PBE def2-SVP def2/JK RIJONX TightOpt TightSCF CPCM(DMSO) KeepDens NormalPrint Grid4 FinalGrid5
%method
    Z_MaxIter 150
end
%scf 
    MaxIter 150 
end

*xyz 0 1
C       4.026085     -1.212051      0.000031
C      -3.313454      0.000008     -0.000131
C       5.413748     -1.205853     -0.000365
C      -4.026091      1.212051     -0.000333
C       6.111884      0.000013     -0.000433
C      -5.413752      1.205847     -0.000315
C       5.413733      1.205871     -0.000018
C      -6.111885     -0.000021     -0.000102
C       4.026071      1.212051      0.000379
C      -5.413728     -1.205875      0.000032
C       3.313452     -0.000005      0.000350
C      -4.026066     -1.212050      0.000016
C       1.893686     -0.000008      0.000523
C       0.678508      0.000000      0.000397
C      -0.678520      0.000009      0.000187
C      -1.893672      0.000016     -0.000043
H       3.483003     -2.148319      0.000143
H       5.951276     -2.145874     -0.000768
H      -3.483012      2.148322     -0.000572
H       7.194627      0.000021     -0.001168
H      -5.951286      2.145866     -0.000334
H       5.951249      2.145899     -0.000149
H      -7.194627     -0.000032      0.000311
H       3.482975      2.148312      0.000762
H      -5.951240     -2.145905      0.000285
H      -3.482967     -2.148309      0.000048

Así, después de correr, vi que de alguna manera Orca fue la "separación" de todos los bonos y llegar a esto:

Sin embargo, mi trabajo sigue sin detectar ningún problema y que sólo termina porque golpea el máximo permitido de iteraciones.

También probé a cambiar el conjunto de base, pero sin éxito. También, he jugado mucho con los funcionales, que ha cambiado entre LDA y GGA sin éxito, semiempírico, no han mostrado ningún buen resultado.

La única cosa que podría utilizar para resolver este problema fue por ajuste de la COPT en lugar de OPTAR sólo, que los cambios Orca a coordenadas cartesianas (pero ya estaba usando coordenadas cartesianas?).

Me pregunto por qué lo que está sucediendo. Me estoy perdiendo algo? Lo que está mal con mi aporte a la causa de Orca para "dividir" todos los átomos separados? Y por qué Orca sólo puede optimizar de esta molécula con COPT activado?

Answer 1

Parece ser un error con la optimización en coordenadas internas. Como ya he mencionado en los comentarios, ORCA no podía convergen en la relajación de la geometría paso a la hora de la optimización se realiza en coordenadas internas. El cambio a la optimización en coordenadas Cartesianas (coordsys cartesian) o viejo interna (coordsys redundant_old) en la %geom sección de la entrada, o el uso de la ! COPT simple entrada de la OP ha hecho, se soluciona el problema.

P. S. OP podría informar de como un error en la ORCA foro para que los desarrolladores pueden tomar un vistazo más de cerca a este caso.