What should I do if the B-matrix is singular?

[willows9803]

I tried to run dscf on a big molecule with 52 atoms. However, the relax found the B-matrix was singular (the internal coordinates were linearly dependent) and stopped abnormally. Is there anything I can do to eliminate this singularity in the Define script?

Thanks a lot!

[christof.haettig]

Hello,

it's quite unusual that define is not able to construct useful redundant internal coordinated with the 'ired' option. But it can happen for molecules with a very complicated structure (several rings and cages, possibly mixed with intermolecular connections are a nightmare). Than that easiest way out is to do the geometry optimization in cartesian coordinates. This is not very efficient, but sometimes the only solution. After some geometry cycles in cartesian coordinates one can try again to construct redundant internal coordinates.

christof

[willows9803]

Yes, first cartesian, then internal, it works!
Thanks a lot, christof!

Hello,

it's quite unusual that define is not able to construct useful redundant internal coordinated with the 'ired' option. But it can happen for molecules with a very complicated structure (several rings and cages, possibly mixed with intermolecular connections are a nightmare). Than that easiest way out is to do the geometry optimization in cartesian coordinates. This is not very efficient, but sometimes the only solution. After some geometry cycles in cartesian coordinates one can try again to construct redundant internal coordinates.

christof