For exploring the preferred site for hydrogen bond formation, theoretical calculations have been performed for a number of six-member, nonaromatic rings allowing for alternative protonation on the ring nitrogens. Gas-phase protonation studies for test molecules indicate that the B3LYP/aug-cc-pvtz and QCISD(T)(CBS) calculations approach the experimental values within about 1 kcal/mol with considerable improvement for relative enthalpies and free energies. Relative free energies I calculated at the IEF-PCM/B3LYP/aug-cc-pvtz level predict favorable protonation on the tertiary rather than on the secondary nitrogen both in aqueous solution and in a dichloromethane solvent for saturated rings. Protonation on a nitrogen atom next to a C=C bond is disfavored due to a large increase in internal energy. Monte Carlo simulations considering a counterion and Ewald summation for the long-range electrostatic effects for a 0.1 molar model system predict Delta G(solv)/MC values generally less negative than from the IEF-PCM calculations. These results make the protonation on the tertiary nitrogen even more favored. The solute-solvent pair-energy distribution depends sensitively on the applied model. In conclusion, the freely moving anion has been considered as the most: relevant model with overall neutrality for the system and applying the least restrictions.