The combination of intramolecular charge constraints with the neglect of molecular electronegativity allows us to maintain the limit that infinitely separated nucleobases will obtain a ground state density equal to n 0 (r). On one hand, this procedure avoids a known problem with fluctuating charge models where infinitely separated molecules can experience long range charge transfer when their electronegativities are different. On the other, all intermolecular charge transfer is also not allowed to occur when two nucleobases approach each other at reasonable interaction distances. This does not mean that the energetic effects of charge transfer are completely ignored. Instead, we model this contribution to the total energy by observing that when large amounts of intermolecular polarization occurs, there should also be an accompanied lowering of energy by sharing electron orbitals of the other molecule. From our DFT calculations (as described in the section 3), we find a strong linear correlation between these two contributions (see supporting information). Additionally, we find that for hydrogen-bonded complexes, the sharing of such orbitals gives the appearance of simply enhancing the polarization of two fragments. To this end, we introduce a phenomenological scaling factor to the polarization energy, s pol-ct . Combining the frozen energy terms in equations 7, 22, and 23, the first and second order kinetic-XC terms in equations 18 -21, and our polarization scaling to the total energy in equation 6, the total energy instead reads