Hello everyone,
I am working with short peptides (~400-500 atoms). When I am doing a minimization using jobex in TurboMole 7.7, usually in ired, sometimes I get this error during statpt: "geoired.f: lowest eigenvalue is too small". I was unable to find any information or what this means.
I am wondering why this is and if there is a way to fix it. Things we have tried and somewhat worked:
- small displacements of atoms before the start of the optimization (within 0.1 A)
- $redund_inp keyword in control set to "metric 2" (this works in majority of cases, but also occasionally fails)
- starting optimization in cartesian coordinates (without ired) for a few steps (tested 25) and then continuing in ired.
Despite all of these workarounds, which mostly work, statpt still occasionally fails in some cases (~1%) and I'm not sure whether this can be fixed or how to fix it. Can someone please elaborate what this error means and what I can do with it?
I'm attaching the control and coord files, as well as statpt output. cood.non_displaced.inp are the original coordinates of the structure, which was pre-optimized using xTB-GFN2 method. Upon visual inspection, I didn't notice any glaring issues with the structure itself. The structure is supposed to have certain restraints (to keep certain atoms together; a tetrahedron created by atoms 52, 118, 400 and 447), which might play some role, but in other similar calculations this didn't cause issues.
Thank you for your time and responses,
jmkormanik