Dielectric Properties at MP2 and CCSD(T) level of theory

From NWChem

You are viewing a single post from the thread title above
Jump to: navigation, search

Click here for full thread
  • Bert Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Profile
    Send PM
Forum Vet
Threads 5
Posts 598
Sven,

This was my mistake. Response properties in TCE are only available for CCSD, CCSD(T), and CCSDT.

As to why the CCSD(T) is further away from experiment then MP2, this might be geometry dependent. Have you tried optimizing at the CCSD(T) level and to get the dipole moment there.

Bert


Quote: Sep 3rd 6:50 pm
Hi Bert!

It seems that the problem is that for mbpt2 the code is not able to do response properties or calculations for the dipole moment. I think you mentioned in your first answer in the forum that response properties are available for MBPT and CC calculations. So I am not sure if I have a mistake in my input, if the calculation is available in general. The input file for the job is:

start LiF

memory total 2500 mb

geometry
 symmetry c1
Li 0 0 0.00
F 0 0 1.56
end

basis spherical
 Li library cc-pvtz
F library cc-pvtz
end

tce
freeze core
MBPT2
io ga
2eorb
dipole
tilesize 16
end

set tce:lineresp T
set tce:afreq 0.000
set tce:respaxis T T T

mp2
 freeze core
end

task mp2 optimize
task mp2 freq numerical

task tce energy

tce
freeze core
ccsd(t)
io ga
2eorb
dipole
tilesize 16
end

task tce energy

The MP2 part is running fine and is giving the SCF and SCF+MP2 dipole moments. If you perform mbpt2 calculations in tce it does not seem to perform any calculations of the dipole moment or polarizability. Only the SCF part is given in the output (agrees with the results of the MP2 routine). On the other hand similar settings produce the desired output for CCSD(T) calculations (by the way I do not understand why the CCSD(T) results for the dipole moment or worse than the MP2 results.).

If you wish I could send you see output file.

Best wishes and thanks in advance

Sven

P.S.: Tried to email you this file but I am not sure if it worked, since a error occurred every time I tried to email. Consequently I am posting it here as well.


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


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