MP2 Correlation energy too big

From NWChem

Viewed 2095 times, With a total of 2 Posts
Jump to: navigation, search

Gets Around
Threads 24
Posts 86
I run RI-MP2 calculation of cubane molecule (C8H8) in a different bases.
cc-pCVTZ, cc-pCVQZ looks good, but with cc-pCV5Z basis I have got unrelailable MP2 energy.

my input file is (molecular geometry is almost equilibrium):
title "cubane cc-pV5Z MP2-RI energy"

geometry
 symmetry Oh
 H                     1.41063149     1.41063149     1.41063149
 C                     0.78371973     0.78371973     0.78371973
end

basis
 H library cc-pV5Z
 C library cc-pCV5Z
end

basis "ri-mp2 basis"
 H library cc-pV5Z-ri
 C library cc-pV5Z-ri
end

task rimp2


output energy is:
          -------------------------------------------
          SCF energy                -307.507255162945
          Correlation energy       -4369.270156917300
          Total MP2 energy         -4676.777412080244
          -------------------------------------------
Edited On 10:24:52 PM PDT - Fri, May 3rd 2013 by Vladimir

Gets Around
Threads 24
Posts 86
Use of spherical functions can do the trick!

title "cubane cc-pV5Z MP2-RI energy"

geometry
symmetry Oh
H 1.41063149 1.41063149 1.41063149
C 0.78371973 0.78371973 0.78371973
end

basis spherical
H library cc-pV5Z
C library cc-pCV5Z
end

basis "ri-mp2 basis"
H library cc-pV5Z-ri
C library cc-pV5Z-ri
end

scf
direct
end

task rimp2


output energy is:
-------------------------------------------
SCF energy                -307.507287526028
Correlation energy          -1.815754093360
Total MP2 energy          -309.323041619388
-------------------------------------------


many thanx to Endlesslake
Edited On 8:40:43 AM PDT - Mon, May 6th 2013 by Vladimir

Gets Around
Threads 3
Posts 85
Dunning basis sets are designed to be used with spherical functions so the cartesian functions cause very bad linear dependency problems in some cases. Spherical functions aren't a trick; they are the right way to run this calculation.


Forum >> NWChem's corner >> Running NWChem



Who's here now Members 0 Guests 1 Bots/Crawler 0


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC