D_e   = 1.6          # eV
alpha = 3.028        # Å^−1
r_e   = 1.411        # Å
dt    = 0.1          # s
k_B   = 8.617333e-5  # eV/K
m     = 18.998403    # amu