Author Topic: reduced mass of diatomic molecule in aoforce output  (Read 288 times)

tonyst

  • Newbie
  • *
  • Posts: 3
  • Karma: +0/-0
reduced mass of diatomic molecule in aoforce output
« on: May 11, 2020, 02:04:30 pm »
Dear Turbomole enthusiasts,

an analytical force constant calculation of IrN yielded the following output:

Code: [Select]
       mode               1        2        3        4        5        6

     frequency           0.00     0.00     0.00     0.00     0.00  1239.80

     symmetry                                                         a1 

        IR                -        -        -        -        -       YES
|dDIP/dQ|   (a.u.)     0.0000   0.0000   0.0000   0.0000   0.0000   0.0056
intensity (km/mol)       0.00     0.00     0.00     0.00     0.00    56.16
intensity (  %   )       0.00     0.00     0.00     0.00     0.00   100.00
 
       RAMAN              -        -        -        -        -       YES
 
  1   ir          x   0.00000  0.00000  1.00000  0.00000  0.00000  0.00000
                  y   0.00000  0.00000  0.00000  0.00000  1.00000  0.00000
                  z   0.00000  0.70711  0.00000  0.00000  0.00000 -0.07268
  2   n           x   0.00000  0.00000  0.00000  1.00000  0.00000  0.00000
                  y   1.00000  0.00000  0.00000  0.00000  0.00000  0.00000
                  z   0.00000  0.70711  0.00000  0.00000  0.00000  0.99736

reduced mass(g/mol)    14.007  103.113  192.220   14.007  192.220   14.948

The reduced mass is calculated to be 14.948, but should it not be 192.220*14.007/(192.220+14.007)=13.056? What am I missing?

Thank you for your help

/Edit: and for N2 it is 14?

Code: [Select]
       mode               1        2        3        4        5        6

     frequency           0.00     0.00     0.00     0.00     0.00  2348.41

     symmetry                                                         a1g

        IR                -        -        -        -        -       NO
|dDIP/dQ|   (a.u.)     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
intensity (km/mol)       0.00     0.00     0.00     0.00     0.00     0.00
intensity (  %   )       0.00     0.00     0.00     0.00     0.00     0.00
 
       RAMAN              -        -        -        -        -       YES
 
  1   n           x   0.00000  0.00000  0.00000  1.00000  0.00000  0.00000
                  y   0.00000  1.00000  0.00000  0.00000  0.00000  0.00000
                  z   0.70711  0.00000  0.00000  0.00000  0.00000  0.70711
  2   n           x   0.00000  0.00000  1.00000  0.00000  0.00000  0.00000
                  y   0.00000  0.00000  0.00000  0.00000  1.00000  0.00000
                  z   0.70711  0.00000  0.00000  0.00000  0.00000 -0.70711

reduced mass(g/mol)    14.003   14.003   14.003   14.003   14.003   14.003

/Edit2:

Taking the wavenumber and force constant from that calculation I derive a reduced mass of 13.055474:

Code: [Select]
You have: (0.7594293 hartree / atomiclength^2)/(1239.80 cm^-1 * c * 2 * pi)^2*N_A
You want: g/mol
* 13.055474

I am really curious what that reduced mass in the aoforce output means.

« Last Edit: May 13, 2020, 10:44:34 am by tonyst »