From NWChem
You are viewing a single post from the thread title above
-
Bert Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
|
|
Forum Vet
Threads 5
Posts 598
|
|
9:43:32 AM PDT - Tue, Sep 6th 2011 |
|
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 pmHi 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.
|
|
|
AWC's:
2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC