Figure 5.1:

Morse potential for H_2 and HF.

Code for Figure 5.1

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
# Converted from morse.m - Morse potential curves
import numpy as np
from misc import save_ascii

r = np.arange(0, 5.01, 0.01)

D, B, re = 109.5, 1.94, 0.742
V1 = D * (np.exp(-2*B*(r-re)) - 2*np.exp(-B*(r-re)))

D, B, re = 140.5, 2.22, 0.917
V2 = D * (np.exp(-2*B*(r-re)) - 2*np.exp(-B*(r-re)))

temp = np.column_stack([r, V1, V2])
save_ascii('morse.dat', temp)