Jump to content

Yx6015IntroToMolecularModelling2

From ChemWiki

NH3 Molecule

Information
NH3 Ammonia
Calculation Details
Method Basis Set Point Group
RB3LYP 6-31G(d,p) C3v
Calculation Results
E(RB3LYP) /a.u. RMS Gradient /a.u. Dipole Moment /D
-56.55776873 0.00000323 1.8465
Click here for the log file.
         Item               Value     Threshold  Converged?
 Maximum Force            0.000004     0.000450     YES
 RMS     Force            0.000004     0.000300     YES
 Maximum Displacement     0.000072     0.001800     YES
 RMS     Displacement     0.000035     0.001200     YES
Vibration and Charge Window
  • N-H Bond Length: 1.01792 Angstrom
  • H-N-H Bond Angle: 105.741°
  • Modes from the 3N-6 Rule: 6
  • Degenerate Modes: 2, 3 and 5, 6
  • Bending Vibration Modes: 1, 2 and 3
  • Stretching Vibration Modes: 4, 5 and 6
  • Highly Symmetric Mode: 4
  • "Umbrella" Mode: 1
  • Number of Experimental Spectrum Bands of Gaseous Ammonia: 4

N2 and H2 Molecule

Information Information
N2 Nitrogen H2 Hydrogen
Calculation Details
Method Basis Set Point Group Method Basis Set Point Group
RB3LYP 6-31G(d,p) D∞h RB3LYP 6-31G(d,p) D∞h
Calculation Results Calculation Results
E(RB3LYP) /a.u. RMS Gradient /a.u. Dipole Moment /D E(RB3LYP) /a.u. RMS Gradient /a.u. Dipole Moment /D
-109.52412868 0.00000060 0 -1.17853936 0.00000017 0
Click here for the log file. Click here for the log file.
 N2      Item               Value     Threshold  Converged?     H2	Item               Value     Threshold  Converged?
 Maximum Force            0.000001     0.000450     YES 	Maximum Force            0.000000     0.000450     YES
 RMS     Force            0.000001     0.000300     YES		RMS     Force            0.000000     0.000300     YES
 Maximum Displacement     0.000000     0.001800     YES		Maximum Displacement     0.000000     0.001800     YES
 RMS     Displacement     0.000000     0.001200     YES		RMS     Displacement     0.000001     0.001200     YES

Reactivity

We can determine the energy for the reaction of N2 + 3H2 -> 2NH3.

  • E(NH3)= -56.55776873 a.u.
  • 2*E(NH3)= -113.11553746 a.u.
  • E(N2)= -109.52412868 a.u.
  • E(H2)= -1.17853936 a.u.
  • 3*E(H2)= -3.53561808 a.u.
  • ΔE=2*E(NH3)-[E(N2)+3*E(H2)]= -0.05579070 a.u.= -146.47 kJ/mol

The products are more stable. The energy change is very different from experimental enthalpy change (ΔH = −92.4 kJ mol−1). This shows the inadequacy of using simple DFT methods to calculate thermodynamic data.

CH4Molecule

Information
CH4 Methane
Calculation Details
Method Basis Set Point Group
RB3LYP 6-31G(d,p) Td
Calculation Results
E(RB3LYP) /a.u. RMS Gradient /a.u. Dipole Moment /D
-40.52401404 0.00003263 0
Click here for the log file.
         Item               Value     Threshold  Converged?
 Maximum Force            0.000063     0.000450     YES
 RMS     Force            0.000034     0.000300     YES
 Maximum Displacement     0.000179     0.001800     YES
 RMS     Displacement     0.000095     0.001200     YES
Vibration properties
Mode Frequency /cm-1 Intensity Mode Frequency /cm-1 Intensity Mode Frequency /cm-1 Intensity
1 1356.20 14.1008 4 1578.58 0 7 3162.33 25.3343
2 1356.20 14.1008 5 1578.58 0 8 3162.33 25.3343
3 1356.20 14.1008 6 3046.46 0 9 3162.33 25.3343
This is a bending motion of the entire molecule, with medium intensity. 3 degenerate modes due to 3 different orientation of bending. These are 2 degenerate symmetric bends and 1 symmetric stretching. There is no change in dipole moment thus resulting in it having 0 intensity and not absorbing any infrared radiation. These are due to 3 different stretching of the molecule with degenerate energy levels. As stretching is more energetic than bending, it occurs at higher energy level with higher intensity.
Charge Distribution
Unsurprisingly, as the center carbon atom is more electronegative than the hydrogen atom, more electron is distributed around the center carbon.
Molecular Orbitals
MO10(E=+0.52915) This is a high energy unoccupied anti-bonding MO, and is similar to MO11 and MO12, with different spatial orientation.
MO9 (E=+0.17677) This is a high energy unoccupied anti-bonding MO, and is similar to MO8 and MO7, with different spatial orientation. Notably, this shows less nodal region as compared to MO10, possibly resulting in its lowering in energy.
MO6 (E=+0.11824) This is a anti-bonding LUMO. This is visually similar to MO2, with the main difference being the size of the orbital (i.e. MO2 is smaller). While no nodal region is shown, the large size of the orbital contributed to its high energy.
MO5 (E=-0.38831) This is a bonding HOMO, and is similar to MO4 and MO3, with different spatial orientation. The nodal region is noticeably smaller, resulting in its low energy.
MO1 (E=-10.16707) This is a very low energy occupied bonding MO and is entirely around the center carbon atom. This is possibly entirely due to the 1s orbital of the carbon atom.

CH3+ and CH3- Molecules

Information Information
CH3+ Carbocation CH3- Carbanion
Calculation Details
Method Basis Set Point Group Method Basis Set Point Group
RB3LYP 6-31G(d,p) D3h RB3LYP 6-31G(d,p) C3v
Calculation Results Calculation Results
E(RB3LYP) /a.u. RMS Gradient /a.u. Dipole Moment /D E(RB3LYP) /a.u. RMS Gradient /a.u. Dipole Moment /D
-39.48468012 0.00010036 0 -39.79602794 0.00010504 1.7169
Click here for the log file. Click here for the log file.
CH3+    Item               Value     Threshold  Converged?     CH3-    Item               Value     Threshold  Converged?
Maximum Force            0.000201     0.000450     YES 	       Maximum Force            0.000173     0.000450     YES
RMS     Force            0.000131     0.000300     YES         RMS     Force            0.000114     0.000300     YES
Maximum Displacement     0.000569     0.001800     YES         Maximum Displacement     0.000510     0.001800     YES
RMS     Displacement     0.000373     0.001200     YES         RMS     Displacement     0.000372     0.001200     YES


Electron Density Comparison
CH3+ Description CH3-
The images are produced as electrostatic potential mapped onto a surface of electron density on the same scale, with blue being electron deficient and red being electron rich. Carbocation is one of the most important intermediate in organic chemistry and reacts as a electrophile. This can be clearly seen by the electron distribution. On the other hand, carbanion has a electron - rich lone pair, and it is the lone pair that provides its nucleophilicity and reactivity.