MRD:DY2917
Molecular Reaction Dynamics
H+H2 System
The experiment is analysed under the condition:
distance | momentum |
---|---|
rab =0.74 | pab =0 |
rbc =2.30 | pbc =-2.7 |

Transition state Definition
Define the following factors: L=length of the linearly placed hydrogen= rab +rbc S=difference between rab and rbc = rab -rbc
On the potential energy surface, the transition state is the saddle point, which is defined by
∂U/∂L=0 ; ∂U/∂L2>0
∂U/∂S=0 ; ∂U/∂S2<0

From Figure 1., the transition state must lie on the line which moves through the lowest energy surface shown on the graph, which in this case is the reaction path, the Bold Black line.
As for a Local Minimum:
∂U/∂L=0 ; ∂U/∂L2>0
∂U/∂S=0 ; ∂U/∂S2>0
From Figure 2., there are several turning points. However, the transition state is the point at the middle where S=0.This is because it is a symmetric reaction, which means that at the transition state rab =rbc . Any other turning points are known as local minimum points. Reason for the reaction pathway to vibrate along the potential surface is due to vibration energy. Good explanation for finding the position of your transition state. Though, your second partial derivatives are missing a squared sign. Sf3014 (talk) 13:00, 14 May 2019 (BST)
Transition State Position

Using the same parameters as above, a Internuclear Distances VS Time graph has been plotted where the transition state is shown as the intercept of the blue and orange line.
Therefore, A B and C are equally separated.
The parameters are then set to:
distance | momentum |
---|---|
rab =1.20 | pab =0 |
rbc =1.20 | pbc =0 |
These parameters are set because the saddle point must satisfy rab =rbc .
The energy surface diagram is then analysed by MEP and the data are exported as Excel.
Since TS has the lowest potential energy,searching for the lowest V would give a corresponding distance:
rab =rbc =rts =0.908 Å
The rab and rbc distances at the transition state can also be found on you internuclear vs time graph at the point you stated. You could also show the internuclear distance vs time graph for parameters at the transition state, to back up your statement. Sf3014 (talk) 13:18, 14 May 2019 (BST)
MEP and Dynamic Trajectory

distance | momentum |
---|---|
rab =0.908 | pab =0 |
rbc =0.918 | pbc =0 |
Both MEP and Dynamic calculation methods have been used to generate the Internuclear Distance VS Time plot.

As can be seen from Fig 3. and Fig 4., the distance of B-C fluctuates after t=0.4 min for dynamic calculation while it is a straight line calculated by MEP.
Dynamic calculation has included the motion of the molecule, which means the fluctuation of the blue line is caused by the actual vibration of H2 molecule. MEP, however, only shows the lowest energy pathway.
Another difference shown on the graph is the gradient of the orange and green line. Shown on Fig 3. The gradient is decreasing with respect to time, this is because as the distance between H2 and H atom increases, the repulsion between them decreases, and therefore, from MEP calculation, it would show that the rate of separation is decreasing. Good comparison between dynamic and MEP calculations. Although, there is some confusion in your reasoning, ie you stated that the B-C distance is gaining vibrational energy when instead it is A-B which is the blue line shown in the graph. Also, how do you know the time is in minutes? I would think t is faster than seconds. Sf3014 (talk) 13:36, 14 May 2019 (BST)


By comparing Fig 5. and Fig 6. it can be clearly seen that Dynamic calculation is based on the conservation of energy, while MEP only consider stating the minimum energy pathway. Good observation but explanation is needed, you should compare your statement with a written description for the behaviour of your groups. Sf3014 (talk) 13:38, 14 May 2019 (BST)
Reactive and unreactive trajectory
pAB | pBC | Etot | Reactive? | Description of the dynamics |
---|---|---|---|---|
-1.25 | -2.5 | -99.0184 | Yes | Vibrating H2 molecule approaches HC atom and successfully interact, HBC molecule vibrates away |
-1.5 | -2.0 | -100.456 | No | Vibrating H2 molecule approaches HC atom but does not have enough energy to overcome the transition state, the reaction does't proceed and HAB molecule vibrates back |
-1.5 | -2.5 | -98.9549 | Yes | Vibrating H2 molecule approaches HC atom and successfully interact, however, at a faster rate |
-2.5 | -5.0 | -84.9557 | No | Vibrating H2 molecule approaches HC atom, reactants have crossed the reaction barrier, transition state is reached where HB atom 'bounces' between HA and HC, however, HAB molecule is reformed and vibrates back |
-2.5 | -5.2 | -84.4160 | Yes | Vibrating H2 molecule approaches HC atom, reactants have crossed the reaction barrier, transition state is reached where HB atom 'bounces' between HA and HC, finally, HBC molecule is formed and vibrates away. |
From these figures, it can be concluded that
1)for reaction to occur, molecules must posses energy higher than the reaction barrier.
2)for a reaction, higher energy the molecule possess, faster the rate of reaction
3)even if the molecules have energy higher than the transition state, the reaction still might not occur.





Well done but if you comment on reaction time you should back it up with a graph that shows the time. Sf3014 (talk) 13:55, 14 May 2019 (BST)
Transition State Theory
Assumption made:
Two basic assumptions:
1.Nuclear motion is stationary compared to electronic motion due to the large difference in mass (equivalent to the Born-Oppenheimer approximation)
2.Reactant molecules are distributed according to the Boltzmann distribution.
Three unique assumptions:
3.If the reactants passes the transition state in the direction of forming product, they can not turn back and reform reactants.
4.In the transition state, the motion of molecules in the direction of reaction coordinates can be separated from the other motions, it is not effected by the other motions and can be treated classically as Translation.
5.Assumption 2 still holds even in the absence of an equilibrium between reactants and products.
As shown on Fig 10., even if the reactants have passed through the transition state, it still reforms back to the reactant, this is against assumption 3.
Where did you get this information from? You should add a reference. Sf3014 (talk) 14:00, 14 May 2019 (BST)
F - H - H system
PES inspections
Reaction type


Reaction | Reaction type | Transition state position |
---|---|---|
F + H2 → HF + H | Exothermic | rFH = 1.813 , rHH = 0.744 |
H + HF → F + H2 | Endothermic | rHH = 1.813 , rHF = 0.744 |

For F + H2 → HF + H reaction, the H-H bond is broken while a new H-F bond is formed.
Since H-F (565 Kj/moL) bond has a larger bond energy than H-H (432 Kj/moL), it suggests that this reaction is exothermic
As for the H + HF → F + H2 reaction, it is the opposite, H-F bond is broken and H-H bond is formed and it is deduced to be endothermic. In graph explanation, state that the energy is a potential energy. Its good practice to define AB and BC distances for your reaction (AB is HF and BC is HH), so your explanation is clear when comparing to the graph. You have no statement explaining fig. 12 and 13, the question asked you to refer to Hammond’s postulate. Also, where did you get these bond energies from? You need to reference.. Sf3014 (talk) 18:20, 14 May 2019 (BST)
Activation energy
The activation energy for this reaction F + H2 → HF + H is: E=0.20 Kcal/moL
The activation energy for this reaction H + HF → F + H2 is: E=30 Kcal/moL


How did you get to these numbers for your activation energies? Some explanation is needed. Sf3014 (talk) 18:21, 14 May 2019 (BST)
Reaction Dynamics
Energy conservation


distance | momentum |
---|---|
rFH = 2.0 | pFH = -0.5 |
rHH = 0.8 | pHH = -1 |

As shown on Fig 18., the potential energy decreases as the kinetic energy increases, the increase in KE corresponds to the heat generated during this exothermic reaction.
Some description needed for fig.17, which shows the transfer of energy between vibrational and translational energies. You forgot to include how this process can be confirmed experimentally. Sf3014 (talk) 18:22, 14 May 2019 (BST)
Polanyi's empirical rules
F + H2 = FH + H is an exothermic reaction, the transition state resembles the reactants and therefore, has a early barrier.
FH + H = F + H2 is an endothermic reaction, the transition state resembles the products and hence, possesses a late barrier.
Polanyi's rules states that:
1) Translation energy is most effective for crossing early barriers
2) Vibration energy is most effective for crossing late barriers
Therefore, it can be predicted that F + H2 = FH + H can occur under Low Vibration High Translation and FH + H = F + H2 under High Vibration Low Translation. This corresponds to the experimental data shown above.
In what way does this correspond to the data above? Give examples. Also, reference is needed for Polanyi’s rules. Sf3014 (talk) 18:23, 14 May 2019 (BST)