User Tools

Site Tools


vibrational_modes

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
vibrational_modes [2011/02/18 13:13]
127.0.0.1 external edit
vibrational_modes [2014/09/12 15:41] (current)
Line 1: Line 1:
 ===== Vibrational modes ===== ===== Vibrational modes =====
  
-In this section we describe how to calculatethe ​vibrational modes of molecules using dynamical matrix. Vibrational modes are obtained from eigen values ​of the dynamical matrix. The elements of the matrix are defined as second ​derivastives ​of energy with respect to geometry. The first derivatives with respect to geometry are forces.The second derivatives are obtained from the different of two forces and displacement. At first calculate optimization structure of the molecule and then move with atom in first direction and do SCF calculation,​ after that move with the same atom but in opposite direction and do again SCF calculation. Then evaluate the matrix element as different of this two forces devided by displacement. This procedure repeat for each direction and each atom of molecule. In order to calculate vibrational modes the dynamical matrix is mass weighted. Vibrational modes are calculated from eigen values ​of dynamical matrix.+In this section we describe how to calculate the vibrational modes of molecules using dynamical matrix. Vibrational modes are obtained from eigenvalues ​of the dynamical matrix. The elements of the matrix are defined as second ​derivatives ​of energy with respect to geometry. The first derivatives with respect to geometry are forces.The second derivatives are obtained from the different of two forces and displacement. At first calculate optimization structure of the molecule and then move with atom in first direction and do SCF calculation,​ after that move with the same atom but in opposite direction and do again SCF calculation. Then evaluate the matrix element as different of this two forces devided by displacement. This procedure repeat for each direction and each atom of molecule. In order to calculate vibrational modes the dynamical matrix is mass weighted. Vibrational modes are calculated from eigenvalues ​of dynamical matrix.
  
 Dynamical matrix calculation set the keywords ''​idynmat''​ in section ''&​OPTION&''​ in input file ''​fireball in''​ Dynamical matrix calculation set the keywords ''​idynmat''​ in section ''&​OPTION&''​ in input file ''​fireball in''​
Line 20: Line 20:
 </​code>​ </​code>​
  
-Parametters ​of dynamical matrix calculation are set in file ''​dyn.optional''​ placed in work directory.+Parameters ​of dynamical matrix calculation are set in file ''​dyn.optional''​ placed in work directory.
  
-=== dyn.optional:​ ===+Calculation of all vibrational modes 
 +=== dyn.optional: ​ === 
 + 
 +<​code>​ 
 +0.03              !elementary displacement 
 +1 1 1             ​!dimension vector 
 +dyn.dat ​          ​!output file   
 +.false. ​           !flag on when only some of atoms are moved, otherwise move all atoms 
 +                      
 +</​code>​ 
 +Calculation of vibrational modes only for selected atoms 
 +=== dyn.optional:  ​===
  
 <​code>​ <​code>​
Line 33: Line 44:
 4                        ​ 4                        ​
 </​code>​ </​code>​
 +
  
 The vibrational modes are writen into file ''​out''​. The vibrational modes are writen into file ''​out''​.
  
 {{:​benzen:​benzen_modes.png|vibrational modes}} {{:​benzen:​benzen_modes.png|vibrational modes}}
vibrational_modes.1298031180.txt.gz · Last modified: 2013/06/21 13:01 (external edit)