Jump to content

Rep:Mod:YYL0508

From ChemWiki

Loh Yong Yao
CID: 00605583

BH3 Optimisation

3-21G Optimisation

Link to BH3 3-21G Optimisation Log File

The process of the optimisation can be observed from the minimisation of total energy and gradient as seen in the graphs below:




The summary of the optimisation is as follows:
File type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 3-21G
Final Energy:-26.46226429 a.u.
Gradient: 0.00008851 a.u.
Dipole Moment: 0.00 D
Point Group: D3h
Job Time: 17.0 s

         Item               Value     Threshold  Converged?
 Maximum Force            0.000220     0.000450     YES
 RMS     Force            0.000106     0.000300     YES
 Maximum Displacement     0.000940     0.001800     YES
 RMS     Displacement     0.000447     0.001200     YES
 Predicted change in Energy=-1.672479D-07
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.1948         -DE/DX =   -0.0002              !
 ! R2    R(1,3)                  1.1944         -DE/DX =   -0.0001              !
 ! R3    R(1,4)                  1.1947         -DE/DX =   -0.0002              !
 ! A1    A(2,1,3)              119.986          -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              120.0157         -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              119.9983         -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------


6-31g(d,p) Optimisation

Link to BH3 6-31g(d,p) Optimisation Log File
The summary of the optimisation is as follows:
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Energy: -26.61532372 a.u.
Gradient: 0.00001746 a.u.
Dipole Moment: 0.00 D
Point Group: D3h
Job Time: 6.0 s

Item               Value     Threshold  Converged?
 Maximum Force            0.000019     0.000450     YES
 RMS     Force            0.000013     0.000300     YES
 Maximum Displacement     0.000090     0.001800     YES
 RMS     Displacement     0.000047     0.001200     YES
 Predicted change in Energy=-2.991560D-09
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.1923         -DE/DX =    0.0                 !
 ! R2    R(1,3)                  1.1923         -DE/DX =    0.0                 !
 ! R3    R(1,4)                  1.1924         -DE/DX =    0.0                 !
 ! A1    A(2,1,3)              120.0126         -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              119.9883         -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              119.9991         -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------

Summary

Optimised BH3 Bond Distance: 1.19 Å (comparable to lit value = 1.19 Å)[1]
Optimised H-B-H Bond Angle: 120.0 ° (consistent with D3h Point Group)
Energy for 3-21G Optimised Structure: -26.46226429 a.u.
Energy for 6-31G(d,p) Optimised Structure: -26.61532372 a.u.

GaBr3 Optimisation

LanL2DZ Optimisation (calculated on HPC)

Link to GaBr3 LanL2DZ Optimisation Log File
The summary of the optimisation is as follows:
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: LANL2DZ
Energy: -41.70082783 a.u.
Gradient: 0.00000016 a.u.
Dipole Moment: 0.00 D
Point Group: D3h
Job Time: 13.9 s

Item               Value     Threshold  Converged?
 Maximum Force            0.000000     0.000450     YES
 RMS     Force            0.000000     0.000300     YES
 Maximum Displacement     0.000003     0.001800     YES
 RMS     Displacement     0.000002     0.001200     YES
 Predicted change in Energy=-1.282685D-12
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  2.3502         -DE/DX =    0.0                 !
 ! R2    R(1,3)                  2.3502         -DE/DX =    0.0                 !
 ! R3    R(1,4)                  2.3502         -DE/DX =    0.0                 !
 ! A1    A(2,1,3)              120.0            -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              120.0            -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              120.0            -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------

Optimised GaBr3 Bond Distance: 2.35 Å (comparable to lit value = 2.3525 Å) [1]
Optimised Br-Ga-Br Bond Angle: 120.0 ° (consistent with D3h Point Group)
Energy for LanL2DZ Optimised Structure: -41.70082783 a.u.

BBr3 Optimisation

6-31G(d,p) + LanL2DZ Optimisation

Link to BBr3 6-31G(d,p) + LanL2DZ Optimisation Log File
The summary of the optimisation is as follows:
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: GEN
Energy: -64.43644917 a.u.
Gradient: 0.00001229 a.u.
Dipole Moment: 0.00 D
Point Group: D3h
Job Time: 20.7 s

 Item               Value     Threshold  Converged?
 Maximum Force            0.000035     0.000450     YES
 RMS     Force            0.000020     0.000300     YES
 Maximum Displacement     0.000207     0.001800     YES
 RMS     Displacement     0.000125     0.001200     YES
 Predicted change in Energy=-5.432275D-09
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.934          -DE/DX =    0.0                 !
 ! R2    R(1,3)                  1.934          -DE/DX =    0.0                 !
 ! R3    R(1,4)                  1.9339         -DE/DX =    0.0                 !
 ! A1    A(2,1,3)              119.9878         -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              120.0008         -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              120.0113         -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------


Optimised B-Br Bond Distance: 1.93 Å (comparable to lit value = 1.893)[1]
Optimised Br-B-Br Bond Angle: 120.0 ° (consistent with D3h Point Group)
Energy for LanL2DZ Optimised Structure: -64.43644917 a.u.

Analysis of Results

Comparison between BH3, GaBr3 and BBr3

Bond Length Comparison
Molecule Bond Length (Å)
BH3
1.19
GaBr3
2.35
BBr3
1.93


H and Br are similar in that they both require one more electron to attain a full valence shell. Comparing H and Br as ligands with B as the central atom, the bond length of central atom to ligand increases from 1.19 Å to 1.93 Å. This is probably due to the fact that B is in the 2nd row of the periodic table and has much better orbital overlap with H than with Br, which is in the 4th row. This is because there is a better match in orbital size between atoms in the same row of the periodic table. Also, the orbitals are more diffuse down the group and overlap less effectively. Also, Br is a much larger atom as compared to H, hence, the internuclear distance between B and Br is larger than that between B and H. Therefore, the B-H bond is shorter than the B-Br bond.

Comparing B and Ga as the central atom, they are both in group 3, but Ga is 2 rows below B, and is in the same row as the ligand Br. Hence, Ga is a much larger atom as compared to B, hence the Ga-Br bond length, which is the internuclear distance, is larger as compared to the internuclear distance between a small B atom and Br. Also, orbitals are more diffused down the group, hence, there is less effective overlap and the bond is weaker and longer. This outweighs the fact that there is better orbital size match between Ga and Br.

A bond can loosely be defined as the overlap of atomic orbitals which allows for the delocalisation of electrons between two nuclei, lowering the energy of the electrons. In some structures, this overlap is very small, possibly due to the distance between the central atom and the ligand and hence, stabilisation is very small. However, there is still a bond between the atoms, albeit a weak one. GaussView draws bonds based on a distance critera, so the fact that GaussView hasn't drawn bonds doesn't mean they are not there. Rather, it just means that the bond is longer than the pre-set distance criteria for a covalent bond.

This is shown in the 7 different strutures obtained when BH3 was optimised using GaussView. Initially, the H atoms are too far apart from the central B atom, and no bond is drawn between the B and H atoms. However, as the optimisation went on, a bond was eventually formed in the later structures as the H atoms were moved closer to the central B atoms and satisfied the distance criteria in GaussView. This also coincides with an increase in orbital overlap and stabilisation of the electrons.

Further Analysis of BH3

Frequency Calculations for BH3

Link to Log File for Frequency Calculations of BH3
The summary of the frequency calculation is as follows:
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Energy: -26.61532364 a.u.
Gradient: 0.00000003 a.u.
Dipole Moment: 0.00 D
Point Group: D3h
Job Time: 6.0 s

Low frequencies ---   -9.3212   -9.3057   -0.0697    0.0001    0.5380    2.6444
Low frequencies --- 1162.9905 1213.1497 1213.1499


 Item               Value     Threshold  Converged?
 Maximum Force            0.000000     0.000002     YES
 RMS     Force            0.000000     0.000001     YES
 Maximum Displacement     0.000000     0.000006     YES
 RMS     Displacement     0.000000     0.000004     YES
 Predicted change in Energy=-1.929123D-14
 Optimization completed.
    -- Stationary point found.

Vibrations

Vibrations of BH3
No.
Form of Vibration
Frequency (cm-1)
Intensity
Symmetry D3h Point Group
1

All 3 H atoms moving up and down perpendicular to the plane of the molecule.
The H atoms move in the same direction at the same time. This is a bending mode.
1163
93
A2"
2

All 3 H atoms moving left and right parallel to the plane of the molecule.
Only 2 H atoms have synchronised movement. The other H atom moves in an opposite direction. This is a bending mode.
1213
14
E'
3

2 H atoms moving towards each other in the plane of the molecule. This is a bending mode.
1213
14
E'
4

All 3 H atoms move in and out in a concerted fashion in the plane of the molecule.
B atom is stationary. This is a symmetric stretching mode.
2583
0
Totally symmetric, A1'
5

2 H atoms move in and out in an alternating fashion in the plane of the molecule. This is an asymmetric stretching mode.
2716
126
E'
6

3 H atoms move in and out in the plane of the molecule.
However, only two of them are synchronized. This is an asymmetric stretching mode.
2716
126
E'


IR Spectrum



There are 6 vibrational modes, however, one of them has 0 intensity, while there are two pairs of degenerate vibrational modes. The mode with 0 intensity is the symmetrical stretch mode, and it has no intensity due to the fact that the stretching does not change the dipole moment of the molecule. There must be a change in dipole moment in order for the vibrational mode to show up on the IR spectrum. Hence, only 3 vibrational modes are observed in the IR spectrum.

Further Analysis of GaBr3

Frequency calculations for GaBr3

Link to Log File for Frequency calculation for GaBr3
Link to Dspace
The summary of the frequency calculation is as follows:
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: LanL2DZ
Energy: -41.70082783 a.u.
Gradient: 0.00000011 a.u.
Dipole Moment: 0.00 D
Point Group: D3h
Job Time: 11.3 s

Low frequencies ---   -0.5252   -0.5247   -0.0024   -0.0010    0.0235    1.2010
Low frequencies ---   76.3744   76.3753   99.6982

Lowest real normal mode is 76.3744 cm-1

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

IR Spectrum

Comparison between the vibrations of BH3 and GaBr3

Comparison between the vibrations of BH3 and GaBr3
Symmetry D3h Point Group
Description
BH3 Frequency (cm-1)
GaBr3 Frequency (cm-1)
A2"
All 3 substitutent atoms are moving up and down perpendicular to the plane of the molecule.
The atoms move in the same direction at the same time.
This is a bending mode.
1162
100
E'
All 3 substituent atoms are moving left and right parallel to the plane of the molecule.
The atoms move in the same direction at the same time.
This is a bending mode.
1213
76
E'
2 substituent atoms moving towards each other in the plane of the molecule.
This is a bending mode.
1213
76
A1'
All 3 substituent atoms move in and out in a concerted fashion in the plane of the molecule.
This is a symmetric stretching mode.
2583
197
E'
2 substituent atoms move in and out in an alternating fashion in the plane of the molecule.
This is an asymmetric stretching mode.
2716
316
E'
All 3 substituent atoms move in and out in the plane of the molecule.
However, only two of them are synchronized.
This is an asymmetric stretching mode.
2716
316


The large difference in values of the frequencies for BH3 compared to GaBr3 indicate a large difference in bond strength. Comparing the frequencies and the vibrational modes of the molecules, there has been a reordering of modes for GaBr3 relative to BH3. The spectra for both molecules are similar in that both have 3 signals, and this is because both have a vibrational mode that has 0 intensity as well as two pairs of degenerate modes.
A2" and E' modes are associated with wagging and scissoring respectively, hence they are of lower energies, while A1' and the other E' modes are associated with symmetrical stretching and asymmetrical stretching respectively, and stretching modes of vibrations are of higher energies.

The frequency analysis is essentially the second derivative of the potential energy surface, while the optimisation when the first derivative of the potential energy surface is zero (i.e. a turning point). Hence the same method and basis set must be used for both the optimisation and frequency analysis to ensure that the frequency analysis is conducted for the same potential energy surface that optimisation was carried out on.

The purpose of carrying out a frequency calculation is twofold. Firstly, it provides predicted IR and Raman modes for comparison with experimental values. Secondly and more importantly, it tells us if the optimisation that we have performed coincides with a minimum point on the potential energy surface. If all the frequencies are positive, then the optimisation is that of the minimum energy. If all the frequencies are negative, then the optimisation is that of a transition state.

The first 6 low frequencies represent the translational and rotational frequencies of the molecule, while the last 3 frequencies represent the vibrational frequencies of the molecule. This is the origin of the 3N - 6 formula for the number of vibrational frequencies of a non-linear molecule.

MO analysis

BH3

Link to Log File
Dspace Link


There are no real significant differences between the real and LCAO MOs, hence, qualitative MO theory can be considered to be reasonably accurate and really useful.

NBO Analysis

6-31G(d,p) Optimisation

Link to NH3 6-31G (d,p) Optimisation Log File
The summary of the optimisation is as follows:
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Energy: -56.55776854 a.u.
Gradient: 0.00003886 a.u.
Dipole Moment: 1.85 D
Point Group: C3v
Job Time: 11.0 s

Item               Value     Threshold  Converged?
 Maximum Force            0.000077     0.000450     YES
 RMS     Force            0.000044     0.000300     YES
 Maximum Displacement     0.000174     0.001800     YES
 RMS     Displacement     0.000112     0.001200     YES
 Predicted change in Energy=-1.601487D-08
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.018          -DE/DX =   -0.0001              !
 ! R2    R(1,3)                  1.018          -DE/DX =   -0.0001              !
 ! R3    R(1,4)                  1.0181         -DE/DX =   -0.0001              !
 ! A1    A(2,1,3)              105.7546         -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              105.7438         -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              105.743          -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)           -111.8739         -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------


Frequency

Link to NH3 Frequency Log File
Link to Dspace
The summary of the frequency calculation is as follows:
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Energy: -56.55776872 a.u.
Gradient: 0.00000029 a.u.
Dipole Moment: 1.85 D
Point Group: D3h
Job Time: 20.7 s

Low frequencies ---   -7.8411   -6.1113   -3.4648   -0.0019   -0.0018   -0.0008
Low frequencies --- 1089.3523 1693.9243 1693.9288


 Item               Value     Threshold  Converged?
 Maximum Force            0.000000     0.000002     YES
 RMS     Force            0.000000     0.000001     YES
 Maximum Displacement     0.000003     0.000006     YES
 RMS     Displacement     0.000001     0.000004     YES
 Predicted change in Energy=-1.142027D-12
 Optimization completed.
    -- Stationary point found.
 

MO

Link to NH3 MO Log File
DSpace Link to population analysis

NBO


Charge Range: -1.125 to 1.125
Charge on N: -1.125
Charge on each H: +0.375

Association Energies

Optimisation

Link to NH3BH3 6-31G(d,p) Optimisation Log File
The summary of the optimisation is as follows:
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Energy: -83.22468911 a.u.
Gradient: 0.00000009 a.u.
Dipole Moment: 5.56 D
Point Group: C1
Job Time: 49.0 s

Item               Value     Threshold  Converged?
 Maximum Force            0.000000     0.000002     YES
 RMS     Force            0.000000     0.000001     YES
 Maximum Displacement     0.000000     0.000006     YES
 RMS     Displacement     0.000000     0.000004     YES
 Predicted change in Energy=-1.536215D-14
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,8)                  1.0185         -DE/DX =    0.0                 !
 ! R2    R(2,8)                  1.0185         -DE/DX =    0.0                 !
 ! R3    R(3,8)                  1.0185         -DE/DX =    0.0                 !
 ! R4    R(4,7)                  1.2098         -DE/DX =    0.0                 !
 ! R5    R(5,7)                  1.2098         -DE/DX =    0.0                 !
 ! R6    R(6,7)                  1.2098         -DE/DX =    0.0                 !
 ! R7    R(7,8)                  1.6677         -DE/DX =    0.0                 !
 ! A1    A(4,7,5)              113.8738         -DE/DX =    0.0                 !
 ! A2    A(4,7,6)              113.8738         -DE/DX =    0.0                 !
 ! A3    A(4,7,8)              104.5974         -DE/DX =    0.0                 !
 ! A4    A(5,7,6)              113.8738         -DE/DX =    0.0                 !
 ! A5    A(5,7,8)              104.5975         -DE/DX =    0.0                 !
 ! A6    A(6,7,8)              104.5975         -DE/DX =    0.0                 !
 ! A7    A(1,8,2)              107.876          -DE/DX =    0.0                 !
 ! A8    A(1,8,3)              107.876          -DE/DX =    0.0                 !
 ! A9    A(1,8,7)              111.0228         -DE/DX =    0.0                 !
 ! A10   A(2,8,3)              107.8758         -DE/DX =    0.0                 !
 ! A11   A(2,8,7)              111.023          -DE/DX =    0.0                 !
 ! A12   A(3,8,7)              111.023          -DE/DX =    0.0                 !
 ! D1    D(4,7,8,1)           -180.0            -DE/DX =    0.0                 !
 ! D2    D(4,7,8,2)            -60.0            -DE/DX =    0.0                 !
 ! D3    D(4,7,8,3)             60.0            -DE/DX =    0.0                 !
 ! D4    D(5,7,8,1)            -60.0            -DE/DX =    0.0                 !
 ! D5    D(5,7,8,2)             60.0            -DE/DX =    0.0                 !
 ! D6    D(5,7,8,3)            180.0            -DE/DX =    0.0                 !
 ! D7    D(6,7,8,1)             60.0            -DE/DX =    0.0                 !
 ! D8    D(6,7,8,2)           -180.0            -DE/DX =    0.0                 !
 ! D9    D(6,7,8,3)            -60.0            -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------

Frequency

Link to NH3BH3 Frequency Log File
Dspace Link
The summary of the frequency calculation is as follows:
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Energy: -83.22468911 a.u.
Gradient: 0.00000023 a.u.
Dipole Moment: 5.56 D
Point Group: C1
Job Time: 104.8 s

Low frequencies ---   -2.9083   -1.3024   -0.0008    0.0007    0.0010    3.6217
Low frequencies ---  263.3601  632.9766  638.4443


Item               Value     Threshold  Converged?
 Maximum Force            0.000000     0.000002     YES
 RMS     Force            0.000000     0.000001     YES
 Maximum Displacement     0.000001     0.000006     YES
 RMS     Displacement     0.000001     0.000004     YES
 Predicted change in Energy=-8.207670D-13
 Optimization completed.
    -- Stationary point found.

Calculation

The dissociation energy can be calculated as follows:
E(NH3) = -56.55776854 a.u.
E(BH3) = -26.61532372 a.u.
E(NH3BH3) = -83.22468911 a.u
ΔE = E(NH3BH3) - [E(NH3) + E(BH3)] = -82.22468911 - (-56.55776854 + -26.61532372) = -0.05159685 a.u. = -135.5 kJ mol-1
This gives the association energy, hence the dissociation energy is +135.5 kJ mol-1
This is a reasonable value for as typical dissociation energies are in the range of 100 to 500 kJ mol-1.
Additionally, I would consider the association between NH3 and BH3 to be a weak one as it is a dative covalent bond, hence, a relatively low value of +135.5 is reflective of a relatively low dissociation energy when compared to typical values, which indicates that the bond is weak.

References

<references> [1]