blob: 4036a404fde0d5fe9f54d298a2780f4f2ed6d9ab (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
//function for calculating the various vibrational energy levels of CO molecule
function[energy]=F(V)
energy=((V+.5)*h*v)/e;
endfunction
//function for converting eV to cm-1
function[energy]=G(V)
energy=((V+.5)*h*v*8065)/e;
endfunction
//variable initialization
m1=12; //molar mass of C atom (amu)
m2=16; //molar mass of O atom (amu)
u=1.68*10^-27; //atomic mass unit (kg)
k=1870 //force constant of CO molecule (N/m)
h=6.6*10^-34; //Plank's constant (joule second)
e=1.602*10^-19; //charge of electron (Coulomb)
//calculation of energy levels
mu=((m1*m2)/(m1+m2))*u; //reduced mass of CO molecule (kg)
v=(1/(2*%pi))*sqrt(k/mu); //frequency of vibration of CO molecule (Hz)
e1=F(0); //energy of 1st level of CO molecule (eV)
E1=G(0); //energy of 1st level of CO molecule (cm-1)
e2=F(1); //energy of 2nd level of CO molecule (eV)
E2=G(1); //energy of 2nd level of CO molecule (cm-1)
e3=F(2); //energy of 3rd level of CO molecule (eV)
E3=G(2); //energy of 3rd level of CO molecule (cm-1)
printf("\nE = %.3f eV, %.3f eV, %.3f eV..........................\n = %.1f cm-1, %.1f cm-1, %.1f cm-1.................",e1,e2,e3,E1,E2,E3);
|
