Physical Chemistry 3: — Chemical Kinetics — - Christian-Albrechts ...

Physical Chemistry 3: — Chemical Kinetics — - Christian-Albrechts ... Physical Chemistry 3: — Chemical Kinetics — - Christian-Albrechts ...

21.10.2014 Views

6.3 Trajectory Calculations 143 6.3 Trajectory Calculations The “exact” dynamics of the molecules on a potential energy hypersurface can be followed by classical trajectory calculations. I Phase Space: 6-dimensional space spanned by the spatial coordinates ( )and the conjugated momenta ( ): ↔ (6.12) where = = · (6.13) I I Hamilton function: = + = 2 + () (6.14) 2 Note that depends only on and depends only on . Equations of motion: • for a 1 particle system: (1) = = = − () y · = − = − since contains only not . (2) y y • for an atomic system: = + = (6.15) (6.16) = 2 2 + (6.17) = 1 × = 1 µ × · 3X =1 · = = · (6.18) (6.19) 2 2 + ( 1 2 3−6 ) (6.20) · = − (6.21) · = (6.22) Thus we obtain a set of coupled differential equations which can be solved using the Runge-Kutta or Gear methods.

6.3 Trajectory Calculations 144 • Thermal rate constants ( ) follow by suitable averaging over the initial conditions ( vibrational states, rotational states, angles, ...). • Major problems: Zero-point energy? Tunneling? Surface crossings? I Figure 6.6: Trajectory calculations. I Figure 6.7: Trajectory calculations.

6.3 Trajectory Calculations 143<br />

6.3 Trajectory Calculations<br />

The “exact” dynamics of the molecules on a potential energy hypersurface can be<br />

followed by classical trajectory calculations.<br />

I Phase Space: 6-dimensional space spanned by the spatial coordinates ( )and<br />

the conjugated momenta ( ):<br />

↔ (6.12)<br />

where<br />

= <br />

<br />

<br />

= <br />

·<br />

(6.13)<br />

I<br />

I<br />

Hamilton function:<br />

= + = 2<br />

+ () (6.14)<br />

2<br />

Note that depends only on and depends only on .<br />

Equations of motion:<br />

• for a 1 particle system:<br />

(1)<br />

= = <br />

<br />

= −<br />

()<br />

<br />

y<br />

·<br />

= − <br />

= − <br />

since contains only not .<br />

(2)<br />

y<br />

y<br />

• for an atomic system:<br />

= + =<br />

(6.15)<br />

(6.16)<br />

= 2<br />

2 + (6.17)<br />

<br />

<br />

= 1 × = 1 µ<br />

× <br />

·<br />

3X<br />

=1<br />

·<br />

= <br />

<br />

= ·<br />

(6.18)<br />

(6.19)<br />

<br />

2 2 + ( 1 2 3−6 ) (6.20)<br />

·<br />

= − <br />

(6.21)<br />

<br />

·<br />

= <br />

(6.22)<br />

<br />

Thus we obtain a set of coupled differential equations which can be solved using<br />

the Runge-Kutta or Gear methods.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!