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