MRD:01535442

Molecular Reaction Dynamics: Applications to Triatomic systems

In this report, we will be investigating the reaction dynamics of two triatomic systems, H-H-H and F-H-H. This includes investigation of their transition states, reaction coordinates and potential energy surfaces, and how these affect the outcome of chemical reactions.

Exercise 1: H + H₂ System

Figure 1: Internuclear Distance vs Time Graph for H + H₂ system, showing estimate for r_ts

On a potential energy surface diagram, how is the transition state mathematically defined?

∂V(r_i)/∂r_i=0 defines the point in the potential enrgy surface diagram where to gradient is zero, being defined as the maximum on the minimum energy curve.

How can the transition state be identified, and how can it be distinguished from a local minimum of the potential energy surface?

The transitions state can be identified by mapping trajectories near the supposed transition state, and observe whether the line goes towards the products or reactants. It can be distinguished from a local minimum by looking at the second derivative, which will be positive when the function is a local minimum, and negative if it is a maximum point.

Report your best estimate of the transition state position (r_ts) and explain your reasoning illustrating it with a “Internuclear Distances vs Time” plot for a relevant trajectory

The best estimate for r_ts (the transition state position) is 90.8 pm. If you take a look at the internuclear distance vs time graph (Figure 1), it shows that across time, there is no oscillatory behaviour in the triatomic system that alters the distance from the r_ts, therefore 90.8 pm is the closest estimate for this position. Note that distance A-B and B-C are the same, so line A-B in the graph is behind line B-C.

Figure 2: Contour plot of system displaced from transtition state. Calculation type: Dynamics

Comment on how the mep and the trajectory you just calculated differ

The minimum energy path or mep takes a path through the potential surface where it does not move up or down the contour, rather it takes the the path of lowest energy througout (see Figure 2). The trajectory calculated by dynamics (Figure 3) follows the same overall path as the mep, though it's much more wavy. This is due to the vibration of the diatomic molecule. Conservation of energy is present in both as the dynamic trajectory oscillates between the same energy contours.

Figure 3: Contour plot of system displaced from transtion state. Calculation type: MEP

Complete the table below by adding the total energy, whether the trajectory is reactive or unreactive, and provide a plot of the trajectory and a small description for what happens along the trajectory.

p1/ g.mol-1.pm.fs-1	p2/ g.mol-1.pm.fs-1	Etot	Reactive?	Description of the dynamics	Illustration of the trajectory
-2.56	-5.1	-414.280	Yes
-3.1	-4.1	-420.077	No
-3.1	-5.1	-413.977	Yes
-5.1	-10.1	-357.277
-5.1	-10.6	-349.477