Author Topic: What should I do if the B-matrix is singular?  (Read 4592 times)

willows9803

  • Newbie
  • *
  • Posts: 9
  • Karma: +0/-0
What should I do if the B-matrix is singular?
« on: February 08, 2008, 03:23:26 pm »
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

  • Global Moderator
  • Sr. Member
  • *****
  • Posts: 258
  • Karma: +0/-0
    • Hattig's Group at the RUB
Re: What should I do if the B-matrix is singular?
« Reply #1 on: February 12, 2008, 09:55:35 pm »
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

  • Newbie
  • *
  • Posts: 9
  • Karma: +0/-0
Re: What should I do if the B-matrix is singular?
« Reply #2 on: February 19, 2008, 02:39:29 pm »
Yes, first cartesian, then internal, it works!
Thanks a lot, christof! :D
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