¿Orca está rompiendo incorrectamente los enlaces? ¿Por qué solo puede optimizar la molécula con COPT activado?

Question

Me gustaría pedir ayuda para usar Orca (v. 4.0.0) para llevar a cabo una optimización simple.

Primero, tengo esta estructura (optimizada principalmente usando el campo de fuerza MMFF94s en Avogadro):

Así que intenté realizar una optimización simple usando Orca y solo como prueba utilicé 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í que después de ejecutarlo, vi que de alguna manera Orca estaba "dividiendo" todos los enlaces y obteniendo esto:

Sin embargo, mi trabajo sigue adelante sin detectar ningún problema y simplemente termina porque alcanza el número máximo de iteraciones permitidas.

También intenté cambiar el conjunto de bases, pero sin éxito. Además, he jugado mucho con funcionales, cambiando entre LDA y GGA sin éxito, los semiempíricos tampoco han mostrado buenos resultados.

Lo único que pude usar para resolver este problema fue configurar COPT en lugar de OPT solamente, lo que cambia Orca a coordenadas cartesianas (¿pero ya estaba utilizando coordenadas cartesianas?).

Solo me pregunto por qué está sucediendo esto. ¿Estoy olvidando algo? ¿Qué hay de malo en mi entrada que hace que Orca "divida" todos los átomos? ¿Y por qué Orca solo puede optimizar esta molécula con COPT activado?

Answer 1

Parece ser un error con la optimización en coordenadas internas. Como mencioné en los comentarios, ORCA no pudo converger en el paso de relajación de la geometría cuando la optimización se realiza en coordenadas internas. Cambiar a la optimización en coordenadas cartesianas (coordsys cartesian) o en las antiguas coordenadas internas (coordsys redundant_old) en la sección %geom de la entrada, o usar la entrada simple ! COPT como ha hecho OP, resuelve el problema.

PD: OP podría informarlo como un error en el foro de ORCA para que los desarrolladores puedan analizar más de cerca este caso.