As seen in Figure , which demonstrates the final energy profiles E int after accounting for E pol and E ct for some hydrogen bonded complexes, the fit resembles the QM energy profiles very well, especially for hydrogen bonding within the plane of the nucleobases. Recalling that the frozen energy E frz is a sum of three terms, the quality of the fit for both E frz and E int also demonstrates the fit for E Pauli and E disp . Frozen energies near the equilibrium points are well represented by our model, and energies away from this position that deviate from the QM values at least have the sought after slope of the profile. We should note that constant energy shifts in ŷ or ẑ, such as with G:U or C:C in Figure are due to shifts in the energy from the orthogonal x direction. Should the energy be improved upon at the equilibrium point, this shift would disappear. Separation Distance (Ang.)