Rotation Entropy for non-Abelian point group sym and CC-methods.

From NWChem

Viewed 814 times, With a total of 3 Posts
Jump to: navigation, search

Gets Around
Threads 15
Posts 51
Hi.

Rotation entropy of molecule depends on its pont group and symmetry number thereof, as described http://cccbdb.nist.gov/thermo.asp (eq 19), but when using CC-methods only abelian point groups is allowed.

In the case of highly symmetric compounds such as cubane or dodecahedrane incorrect values of rotation entropy are obtained.

Gets Around
Threads 15
Posts 51
reusing existing hessian with proper symmetry group solve the problem.

Gets Around
Threads 11
Posts 98
Hi Vladimir,

I am glad that you were able to solve this problem. Could you post an example input file demonstrating your solution? I will want to perform similar calculations in the future.

Thanks,
Matt

Gets Around
Threads 15
Posts 51
Quote:Mernst Mar 19th 1:34 pm
Hi Vladimir,

I am glad that you were able to solve this problem. Could you post an example input file demonstrating your solution? I will want to perform similar calculations in the future.

Thanks,
Matt

Hi Matt, here is an example.
As described in http://www.nwchem-sw.org/index.php/Vibration I reuse the previously computed hessian by specifying "reuse" in the "freq" module input block.
title "N8-cubane CCSD/cc-pVDZ hessian"
memory stack 100 mb heap 100 mb global 100 mb
geometry
 symmetry Oh
 N                     0.76104280     0.76104280     0.76104278
 N                    -0.76104280    -0.76104280     0.76104278
 N                    -0.76104280     0.76104280    -0.76104278
 N                     0.76104280    -0.76104280    -0.76104278
 N                    -0.76104280    -0.76104280    -0.76104278
 N                     0.76104280     0.76104280    -0.76104278
 N                     0.76104280    -0.76104280     0.76104278
 N                    -0.76104280     0.76104280     0.76104278
end
freq
 reuse N8.hess
end
task tce freq


It is important that the order of the atoms during the calculation of hessian and its subsequent reusing were exactly the same. Check it out in the output file.


Forum >> NWChem's corner >> NWChem functionality



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