Jump to content

Rep:Mod:ROBFINN

From ChemWiki

NH3 Molecule

Summary

Calculation Method RB3LYP
Basic set 6-31G(d,p)
Final Energy / au -56.55776873
RSM Gradient 0.00000485
Point Group C3V
N-H bond distance / Å 1.01798
H-N-H bond angle / ° 105.741

Item Table

         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
Predicted change in Energy=-5.986260D-10
 Optimization completed.
    -- Stationary point found.

Jmol Image

NH3 Molecule

The optimisation file is liked to here

Frequency Analysis

wavenumber / cm-1 intensity / arbitrary units image
1090 145
1694 14
1694 14
3461 1
3590 0
3590 0

From the 3N-6 rule, six modes would be expected. Mode 2 and 3, and 5 and 6 are degenerate. Modes 1, 2 and 3 are bending vibrations. Modes 4, 5 and 6 are stretching vibrations. Mode 4 is highly symmetrical. Mode 1 is the umbrella mode. You would expect to see two bands in an experimental spectrum of gaseous ammonia. Mode four is symmetrical, so won't produce a dipole moment so won't be seen. Modes 2&3 and 5&6 are degenerate so will appear as two peaks. Mode 5&6 have an intensity of zero so won't be seen.

Charge Analysis

The charge on the nitrogen is -1.125 and the charge on the hydrogen is +0.375. This is due to nitrogen being more electronegative than hydrogen, so nitrogen has a higher tendency to attract the electron density from the covalent bond, making the nitrogen more negative.

H2 Molecule

Summary

Calculation Method RB3LYP
Basic set 6-31G(d,p)
Final Energy / au -1.17853936
RSM Gradient 0.00000017
Point Group Dinfh
H-H bond distance / Å 0.74279

Item Table

         Item               Value     Threshold  Converged?
 Maximum Force            0.000000     0.000450     YES
 RMS     Force            0.000000     0.000300     YES
 Maximum Displacement     0.000000     0.001800     YES
 RMS     Displacement     0.000001     0.001200     YES
 Predicted change in Energy=-1.164080D-13
 Optimization completed.
    -- Stationary point found.

Jmol Image

H2 Molecule

The optimisation file is liked to here

Frequency Analysis

wavenumber / cm-1 intensity / arbitrary units
4466 0

Linear molecules use the rule 3N-5. Therefore, we can see there should only be one vibrational mode.

Charge Analysis

The charge on both hydrogens is zero, as they have the same electronegativity.

N2 Molecule

Summary

Calculation Method RB3LYP
Basic set 6-31G(d,p)
Final Energy / au -109.52412868
RSM Gradient 0.00000060
Point Group Dinfh
N-N bond distance / Å 1.10550

Item Table

         Item               Value     Threshold  Converged?
 Maximum Force            0.000001     0.000450     YES
 RMS     Force            0.000001     0.000300     YES
 Maximum Displacement     0.000000     0.001800     YES
 RMS     Displacement     0.000000     0.001200     YES
 Predicted change in Energy=-3.400991D-13
 Optimization completed.
    -- Stationary point found.

Jmol Image

N2 Molecule

The optimisation file is liked to here

Frequency Analysis

wavenumber / cm-1 intensity / arbitrary units
2457 0

Linear molecules use the rule 3N-5. Therefore, we can see there should only be one vibrational mode.

Charge Analysis

The charge on both nitrogens is zero, as they have the same electronegativity.

trans-bis(1-(diethylphosphino)-N-((diethylphosphino)methyl)-N-(2,6-difluorobenzyl)methanamine)-bis(dinitrogen)-chromium [AQEZIH]

[Jonathan D. Egbert, Molly O'Hagan, Eric S. Wiedner, R. Morris Bullock, Nicholas A. Piro, W. Scott Kassel, Michael T. Mock, Chemical Communications, 2016, 52, 9343, DOI: 10.1039/C6CC03449G]

N-N bond distance is equal to 1.127 Å. The crystal structure has a longer N-N bond length when compared to the computational distance. This implies that the bond is weaker. This could be due to the central chromium being electron withdrawing, which makes the bond weaker and consequently longer. The crystal structure is also a solid, so more electron withdrawing groups could interact with the N-N bond, spreading electron density over more atoms, making the bond weaker and longer.

Haber-Bosch Process

Molecule Energy / au
NH3 -56.55777
N2 -109.52413
H2 -1.17854

N2 + 3H2 → 2NH3

ΔE = [2xE(NH3)] - [E(N2) + 3xE(H2)]

  = [-113.11554] - [-109.52413 - 3.53562]
  = -0.05579 au
  = -146.5 kJ/mol


The reaction is exothermic, therefore the product (ammonia) is more stable as energy is released during the reaction.

CH4 Molecule

Summary

Calculation Method RB3LYP
Basic set 6-31G(d,p)
Final Energy / au -40.52401404
RSM Gradient 0.00003263
Point Group Td
C-H bond distance / Å 1.09197
H-C-H bond angle / ° 109.471

Item Table

         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
 Predicted change in Energy=-2.256043D-08
 Optimization completed.
    -- Stationary point found.

Jmol Image

CH4 Molecule

The optimisation file is liked to here

Frequency Analysis

wavenumber / cm-1 intensity / arbitrary units
1356 14
1356 14
1356 14
1579 0
1579 0
3046 0
3162 25
3162 25
3162 25

From the 3N-6 rule, nine modes would be expected. Mode 1&2&3, and 4&5, and 7&8&9 are degenerate. Modes 1, 2, 3, 4 and 5 are bending vibrations. Modes 6, 7, 8 and 9 are stretching vibrations. Mode 6 is highly symmetrical. You would see two peaks for methane. Modes 1&2&3 and 7&8&9 are degenerate so will produce two peaks. The rest of the modes have zero intensity so won't be seen.

Charge Analysis

The charge on the carbon is -0.930 and the charge on the hydrogen is +0.233. This is due to carbon being more electronegative than hydrogen, so carbon has a higher tendency to attract the electron density from the covalent bond, making the carbon more negative.

Molecular Orbitals

HOMO

These three molecular orbitals are degenerate with an energy of -0.38831 au. These are the highest occupied molecular orbital in the molecule, and are each occupied by two electrons with opposite spins. These are bonding orbitals and form σ bonds. Each MO has contribution from 1s atomic orbital from hydrogen and a 2p atomic orbital from carbon, which also forms an antibonding σ MO. This does contribute to the bond order of the molecule.

LUMO

This is the lowest unoccupied molecular orbital in the molecule, with an energy of 0.11824 au. It is unoccupied. The 2s AO of carbon and the 1s AO of hydrogen combind to form this MO and a σ MO. It is an antibonding orbital. This doesn't contribute to the bond order of the molecule.

Carbon 1s

This MO is contributed to by the 1s AO in carbon, with an energy of -10.16707 au. Due to it's highly negative energy, it doesn't contribute to the bonding.

This id the first σ molecular orbital with an energy of -0.69041 au. This is occupied with two electrons with opposite spins and is a bonding orbital. The 2s AO from carbon and the 1s AO from hydrogen contribute to this. This does contribute to the bond order of the molecule.

***

These are three degenerate σ* molecular orbitals. They are unoccupied and each have an energy of 0.17677 au. A 2p AO from carbon and a 1s AO from hydrogen combind to form each MO. This doesn't contribute to the bond order of the molecule.

Independent Cl2 Molecule

Summary

Calculation Method RB3LYP
Basic set 6-31G(d,p)
Final Energy / au -920.34987886
RSM Gradient 0.00002511
Point Group Dinfh
Cl-Cl bond distance / Å 2.04174

Item Table

         Item               Value     Threshold  Converged?
 Maximum Force            0.000043     0.000450     YES
 RMS     Force            0.000043     0.000300     YES
 Maximum Displacement     0.000121     0.001800     YES
 RMS     Displacement     0.000172     0.001200     YES
 Predicted change in Energy=-5.277255D-09
 Optimization completed.
    -- Stationary point found.

Jmol Image

Cl2 Molecule

The optimisation file is liked to here

Frequency Analysis

wavenumber / cm-1 intensity / arbitrary units
520 0

Linear molecules use the rule 3N-5. Therefore, we can see there should only be one vibrational mode.

Charge Analysis

The charge on both chlorines is zero, as they have the same electronegativity.

Marking

Note: All grades and comments are provisional and subject to change until your grades are officially returned via blackboard. Please do not contact anyone about anything to do with the marking of this lab until you have received your grade from blackboard.

Wiki structure and presentation 1/1

Is your wiki page clear and easy to follow, with consistent formatting?

YES, overall a very good wiki well done!

Do you effectively use tables, figures and subheadings to communicate your work?

YES

NH3 1/1

Have you completed the calculation and given a link to the file?

YES

Have you included summary and item tables in your wiki?

YES

Have you included a 3d jmol file or an image of the finished structure?

YES

Have you included the bond lengths and angles asked for?

YES

Have you included the “display vibrations” table?

YES

Have you added a table to your wiki listing the wavenumber and intensity of each vibration?

YES

Did you do the optional extra of adding images of the vibrations?

YES

Have you included answers to the questions about vibrations and charges in the lab script?

YES

N2 and H2 0.5/0.5

Have you completed the calculations and included all relevant information? (summary, item table, structural information, jmol image, vibrations and charges)

YES

Crystal structure comparison 0.5/0.5

Have you included a link to a structure from the CCDC that includes a coordinated N2 or H2 molecule?

YES

Have you compared your optimised bond distance to the crystal structure bond distance?

YES

Haber-Bosch reaction energy calculation 1/1

Have you correctly calculated the energies asked for? ΔE=2*E(NH3)-[E(N2)+3*E(H2)]

YES

Have you reported your answers to the correct number of decimal places?

YES

Do your energies have the correct +/- sign?

YES

Have you answered the question, Identify which is more stable the gaseous reactants or the ammonia product?

YES

Your choice of small molecule 4/5

Have you completed the calculation and included all relevant information?

YES

Have you added information about MOs and charges on atoms?

YES, and overall your explanations are correct, clear and well presented, well done! The aspect you could have improved upon would be to explain that the LUMO and the 1 sigma MOs are bonding/antibonding counterparts. And the same for the HOMOs and the 7,8,9 sigma star MOs.

Independence 1/1

If you have finished everything else and have spare time in the lab you could:

Check one of your results against the literature, or

Do an extra calculation on another small molecule, or

Do some deeper analysis on your results so far

You calculated Cl2 independently well done!