Inorganic Computational Lab:gab116
EX3 section
BH3
B3LYP/6-31G(d,p) level
Item Value Threshold Converged? Maximum Force 0.000018 0.000450 YES RMS Force 0.000009 0.000300 YES Maximum Displacement 0.000070 0.001800 YES RMS Displacement 0.000035 0.001200 YES
Frequency analysis log file GAB116_BH3_FREQ.LOG
Low frequencies --- -10.3498 -3.4492 -1.2454 -0.0056 0.4779 3.2165 Low frequencies --- 1162.9519 1213.1527 1213.1554
optimised BH3 molecule |
Vibrational spectrum for BH3
wavenumber (cm-1) | Intensity (arbitrary units) | symmetry | IR active? | type |
1163 | 93 | A2" | yes | out-of-plane bend |
1213 | 14 | E' | very slight | bend |
1213 | 14 | E' | very slight | bend |
2582 | 0 | A1' | no | symmetric stretch |
2716 | 126 | E' | yes | asymmetric stretch |
2716 | 126 | E' | yes | asymmetric stretch |
There are fewer vibrational peaks than there are vibrations as not all vibrations are IR active, e.g. at 2582 cm-1. In addition, there are two pairs of vibrations that occur at the same wavenumber and hence overlap (at 2716 cm-1 and 1213 cm-1). As a result, there are three peaks in the spectrum even though there are six vibrations.
MO Diagram for BH3
Source: http://www.huntresearchgroup.org.uk/teaching/teaching_comp_lab_year2a/Tut_MO_diagram_BH3.pdf
Are there any significant differences between the real and LCAO MOs?
The real MOs are more diffuse than the LCAO MOs.
What does this say about the accuracy and usefulness of qualitative MO theory?
While qualitative MO theory is useful in terms of accurately predicting the MO's location, basic shape and relative AO contribution, it is not able to predict the volume occupied by the real MOs.
Smf115 (talk) 11:49, 28 May 2018 (BST)Clear inclusion of the calculated MOs in the diagram and nice evaluation of the useful nature of qualitative MO theory. To improve, there could have been more discussion on the similarities and differences of the calculated and qualitative MOs, for example, noticing differences in AO contirbutions particularly in the 2e' and 3a1' MOs.
NH3
B3LYP/6-31G(d,p) level
Item Value Threshold Converged? Maximum Force 0.000013 0.000450 YES RMS Force 0.000006 0.000300 YES Maximum Displacement 0.000039 0.001800 YES RMS Displacement 0.000013 0.001200 YES
Frequency analysis log file GAB116_NH3_FREQ.LOG
Low frequencies --- -8.4661 -8.4184 -0.0028 0.0337 0.1931 26.4322 Low frequencies --- 1089.7605 1694.1862 1694.1866
optimised NH3 molecule |
NH3BH3
B3LYP/6-31G(d,p) level
Item Value Threshold Converged? Maximum Force 0.000114 0.000450 YES RMS Force 0.000063 0.000300 YES Maximum Displacement 0.000621 0.001800 YES RMS Displacement 0.000355 0.001200 YES
Frequency analysis log file GAB116_NH3BH3_FREQ.LOG
Low frequencies --- -0.0617 -0.0457 -0.0066 21.6788 21.6848 40.5422 Low frequencies --- 266.0173 632.3610 640.1362
optimised NH3BH3 molecule |
B-N Dative Bond Energy in NH3BH3
NH3 + BH3 -> H3B:NH3
E(NH3)= -56.55777 au to 5 decimal places
E(BH3)= -26.61532 au to 5 d.p.
E(NH3BH3)= -83.22469 au to 5 d.p.
ΔE=E(NH3BH3)-[E(NH3)+E(BH3)]= -0.05160 au to 5 d.p. = -135 kJ/mol (to the nearest 1 kJ/mol)
(1.00000 au = 2625.5002 kJ/mol, source: http://www.colby.edu/chemistry/PChem/Hartree.html)
Based on your energy calculation is the B-N dative bond weak, medium or strong? What comparison have you made to come to this conclusion?
The B-N dative bond is medium strength, compared to the strong C-C bond in the similarly structured ethane and the weak intermolecular Hydrogen bond in water.
C-C bond energy= 348 kJ/mol, (source: P. Atkins and J. de Paula, Physical Chemistry, 10th Edition, 2014, OUP, 986). Hydrogen bond energy in water= 22 kJ/mol, (source: P. Atkins and J. de Paula, Physical Chemistry, 10th Edition, 2014, OUP, 674).
Smf115 (talk) 12:13, 28 May 2018 (BST)Correct calculation and nice consideration given to the accuracy of the reported energy values. Good comparisons made and referenced literature values.
BBr3
B3LYP/6-31G(d,p) LANL2DZ level
Item Value Threshold Converged? Maximum Force 0.000008 0.000450 YES RMS Force 0.000004 0.000300 YES Maximum Displacement 0.000036 0.001800 YES RMS Displacement 0.000018 0.001200 YES
Frequency analysis log file Gab116_BBr3_FREQ_scanserver.LOG
Low frequencies --- -0.0137 -0.0064 -0.0046 2.4315 2.4315 4.8421 Low frequencies --- 155.9631 155.9651 267.7052
optimised BBr3 molecule |
Project section: Aromaticity
Benzene
B3LYP/6-31G(d,p) level
Item Value Threshold Converged? Maximum Force 0.000197 0.000450 YES RMS Force 0.000085 0.000300 YES Maximum Displacement 0.000780 0.001800 YES RMS Displacement 0.000333 0.001200 YES
Frequency analysis log file GAB116_BENZENE_FREQ.LOG
Low frequencies --- -3.5606 -3.5606 -0.0088 -0.0043 -0.0043 10.0905 Low frequencies --- 413.9582 413.9582 621.1416
optimised Benzene molecule |
Borazine
B3LYP/6-31G(d,p) level
Item Value Threshold Converged? Maximum Force 0.000218 0.000450 YES RMS Force 0.000069 0.000300 YES Maximum Displacement 0.000332 0.001800 YES RMS Displacement 0.000106 0.001200 YES
Frequency analysis log file GAB116_BORAZINE_FREQ.LOG
Low frequencies --- -12.5959 -12.3825 -8.9170 -0.0099 0.0381 0.0783 Low frequencies --- 289.1143 289.1234 403.9108
optimised Borazine molecule |
blue = nitrogen
pink = boron
NBO Charge Analysis
Smf115 (talk) 11:40, 1 June 2018 (BST)Good charge analysis with the same colour range clearly used to highlight the charge distributions across both molecules. Further analysis could have considered the symmetry or the overall charges of the molecule.
Molecular Orbital Comparison
The Concept of Aromaticity
Aromatic systems generally have the following properties (source: https://onlinelibrary.wiley.com/doi/epdf/10.1002/chem.200700250):
1. They posses a resonance stability that mean their bond energies are more stable than their alkene equivalents and have bond lengths that are between single and double bonds.
2. When placed in an external magnetic field, a pi electron ring current is induced.
3. Aromatics are involved in reactions where the pi structure is conserved.
Simple aromatic molecules obey the empirical Hückel's rule (source: J. Clayden, N. Greeves and S. Warren, Organic Chemistry, 2nd Edition, 2012, OUP, 161); for planar compounds with a cyclic, contiguous array of p-orbitals perpendicular to the plane of the rings, there must be 4n+2 p electrons, where n is a natural number. For example, benzene is sp2 hybridized and so there are 6 pz orbitals perpendicular to the plane and 6 p electrons involved in the aromaticity and hence n=1 and Hückel's rule is obeyed.
This rule leads to a simple description of aromaticity as the delocalised pi overlap of pz orbitals above and below the plane, as is shown in the image of MO 17 below.
In fact, bonding in aromatics involves a wide variety of AOs and overlaps, as can be seen from the MO comparison above. Not all aromatic MOs involve only pi orbitals and many can involve sigma interactions, including MO 7 above.
More complex interpretations of aromaticity can include sigma-only aromaticity. For example, there are main group hydride clusters that exhibit no pi aromaticity and only sigma aromaticity. σ-aromaticity and σ-antiaromaticity in saturated inorganic rings is examined in the following paper: Z.-H. Li, D. Moran, K.-N. Fan and P. von Ragué Schleyer, J. Phys. Chem. A, 2005, 109 (16), 3711–3716 (https://pubs.acs.org/doi/abs/10.1021/jp048541o).
Smf115 (talk) 11:38, 1 June 2018 (BST)Good understanding of some of the key concepts of aromaticity and reference to the MOs visualised to illustrate sigma aromaticity. To improve, the answer could have been more developed in places with a greater discussion of the properties stated, the quantitative approaches to aromaticity or around where the classical picture breaks down for example.
Smf115 (talk) 11:38, 1 June 2018 (BST)Overall , a good wiki report with a solid first section in particular.