vibration frequency of H2

[jia... at vt.edu]

Dear all,

I optimized the H2 system with tighter convergence criteria
    EPS_SCF 1E-8
    MAX_FORCE 5E-5

Then I use following script to calculate the vibration frequency:

&VIBRATIONAL_ANALYSIS
   INTENSITIES
   NPROC_REP 16
   DX 0.001
   &MODE_SELECTIVE
     &INVOLVED_ATOMS
       INVOLVED_ATOMS 1 2
     &END INVOLVED_ATOMS
   &END MODE_SELECTIVE
   &PRINT
     &PROGRAM_RUN_INFO ON
     &END
   &END
 &END

I get the vibration is 429.830  1221.671  2595.744  5189.575


But when I delete the &MODE_SELECTIVE session using 
&VIBRATIONAL_ANALYSIS
   INTENSITIES
   NPROC_REP 16
   DX 0.001
   &PRINT
     &PROGRAM_RUN_INFO ON
     &END
   &END
 &END
I get different vibration 22.484153  22.484154  3556.728838

How does the MODE_selective affects the vibration?

I have attached the input and output files for geometry optimization and 
vibration frequency analysis.

Sorry to bother you and thanks for your help! 
Jiamin Wang

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20170304/ba322e2b/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: H2_test.zip
Type: application/zip
Size: 107472 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20170304/ba322e2b/attachment.zip>