Hi NWChem users,
I am trying to calculate the vibration frequency of a configuration where adatom moving near the optimized position on CNT surface. I first optimized the whole stucture (221 atoms), and then restart the job with:
restart cnt_5_10p1vib
set gen_hess:actlist 221 # 221st atom is the adatom
set hessian:compress .true.
vib
reuse
temp 1 300
animate
end
task dft freq numerical
I have two questions:
1. Even thought I put those two set directives in the restart job (as suggested in http://www.nwchemsw.org/index.php/Special:AWCforum/st/id518/printing_hessian_with_greater...), I am still getting all atom frequencies with a 663x663 matrix. Is there any trick here to freeze the atoms in CNT and actually have only adatom moving?
2. I also tried to restart again but this time I changed the temperature to 2000 K. However, the output frequencies of normal modes are exactly the same with those at 300 K. The only difference in the output file is in this block:
Temperature = 300.00K
frequency scaling parameter = 1.0000
ZeroPoint correction to Energy = 1852.825 kcal/mol ( 2.952665 au)
Thermal correction to Energy = 1856.996 kcal/mol ( 2.959312 au)
Thermal correction to Enthalpy = 1857.592 kcal/mol ( 2.960261 au)
Total Entropy = 58.653 cal/molK
 Translational = 49.102 cal/molK (mol. weight =2322.2348)
 Rotational = 0.000 cal/molK (symmetry # = 1)
 Vibrational = 9.551 cal/molK
Cv (constant volume heat capacity) = 54.633 cal/molK
 Translational = 2.979 cal/molK
 Rotational = 1.986 cal/molK
 Vibrational = 49.668 cal/molK
Temperature = 2000.00K
frequency scaling parameter = 1.0000
ZeroPoint correction to Energy = 1852.825 kcal/mol ( 2.952665 au)
Thermal correction to Energy = 3162.591 kcal/mol ( 5.039911 au)
Thermal correction to Enthalpy = 3166.564 kcal/mol ( 5.046242 au)
Total Entropy = 1191.330 cal/molK
 Translational = 58.522 cal/molK (mol. weight =2322.2348)
 Rotational = 0.000 cal/molK (symmetry # = 1)
 Vibrational = 1132.808 cal/molK
Cv (constant volume heat capacity) = 1095.424 cal/molK
 Translational = 2.979 cal/molK
 Rotational = 1.986 cal/molK
 Vibrational = 1090.458 cal/molK
Why are those frequencies the same at different temperature? Is there anything wrong in my procedure?
Thank you so much!
Longtao
