basis set error

  • ccummins
  • Topic Author
  • Offline
  • New Member
  • New Member
More
7 months 2 weeks ago #1386 by ccummins
basis set error was created by ccummins
I had a problem trying to run a calculation on SnCl3 anion.  I used the same input as found in the following test input file:

MINP_SnCl2_def2-SVP_def2-ECP-28_DF-HF_geomopt

I just changed the structure to match SnCl3 and made the charge=-1.  I am getting the following error: 
 Basis set cc-pvtz for atom 1 (Sn) is not found!
 
 Fatal error in exec integ.
 Program will stop.

This is perplexing, as the MINP uses basis=def2-SVP.

For completeness, here is my MINP file:

# DF-HF geometry optimization for SnCl2 with the def2-SVP basis set and the def2-ECP-28 pseudopotential for Sn
basis=def2-SVP
calc=df-hf
mem=3000MB
gopt=full
charge=-1
orbloco=ibo

geom=xyz
4
model.mopin
Sn        0.000000000000        0.000000000000        0.000000000000
Cl        2.390000000000        0.000000000000        0.000000000000
Cl       -0.415020000000        0.000000000000       -2.353690000000
Cl       -0.415020000000       -2.038360000000        1.176850000000

(Note, I was having a similar problem trying to carry out a calculation using a Mo atom.  there must be something wrong with how I am setting the basis, but I have not been able to solve it yet.)

Thanks for any assistance with this!


 

Please Log in or Create an account to join the conversation.

  • kallay
  • Offline
  • Administrator
  • Administrator
  • Mihaly Kallay
More
7 months 2 weeks ago #1387 by kallay
Replied by kallay on topic basis set error
The cc-pVTZ basis set is needed for the IBO localization since the minimal basis required for IBOs is constructed from this basis set. If you remove orbloco=ibo, it will run smoothly.

Best regards,
Mihaly Kallay

Please Log in or Create an account to join the conversation.

  • ccummins
  • Topic Author
  • Offline
  • New Member
  • New Member
More
7 months 2 weeks ago #1388 by ccummins
Replied by ccummins on topic basis set error
Does this mean I cannot construct IBOs for heavy elements such as Sn? On the basis set exchange, I do not see cc-pVTZ available for Sn, although cc-pVTZ-PP is available there. Thanks!!

Please Log in or Create an account to join the conversation.

  • kallay
  • Offline
  • Administrator
  • Administrator
  • Mihaly Kallay
More
7 months 2 weeks ago #1389 by kallay
Replied by kallay on topic basis set error
I have not yet used IBOs for heavy element. Please check the literature what the recommendation is concerning heavy elements, and you can add the required basis sets by hand. Probably cc-pVTZ-PP is fine.

Best regards,
Mihaly Kallay

Please Log in or Create an account to join the conversation.

  • ccummins
  • Topic Author
  • Offline
  • New Member
  • New Member
More
7 months 1 week ago #1391 by ccummins
Replied by ccummins on topic basis set error
Thanks. How would I tell mrcc to look for cc-pVTZ-PP instead of cc-pVTZ when running orbloco=ibo? I don't see a keyword for that.

I have also a theory question: when running calc=LNO-CCSD(T) and orbloco=ibo or orbloco=pm, are the resulting localized orbitals in the MOLDEN file constructed with inclusion of electron correlation? I understand that mrcc may be generating the orbitals in the following way: density matrix -> cholesky guess -> ibo or pm. Is this correct?

Thanks very much.
The following user(s) said Thank You: kaisheng

Please Log in or Create an account to join the conversation.

  • kallay
  • Offline
  • Administrator
  • Administrator
  • Mihaly Kallay
More
7 months 1 week ago #1392 by kallay
Replied by kallay on topic basis set error
For the moment, there no simple way to tell the program to use the PP basis sets. The corresponding code that generates the minimal basis set should be rewritten.
The IBOs are localized HF orbitals, correlation is not included. Yes, Cholesky localized orbitals are used as initial guess, but the HF density is used for that purpose.

Best regards,
Mihaly Kallay

Please Log in or Create an account to join the conversation.

Time to create page: 0.104 seconds
Powered by Kunena Forum