The analysis of intramolecular atomic energies and charges, obtained from QTAIM, strongly and unambiguously suggests that the most stable HL + form of a linear triamine H 2 N(CH 2 ) 2 NH(CH 2 ) 3 NH 2 (L) has the secondary nitrogen atom protonated, contradicting earlier reports. This is supported by the G(aq) value (at the RX3LYP/6-311++G(d,p) level of theory in solvent, CPCM/UAKS) which is lower by about 4 kcal mol À1 (equivalent to three log units of a protonation constant log K H ) when compared with remaining HL + tautomers. The most stable H 2 L 2+ form of L found from QTAIM has two terminal nitrogen atoms protonated, which is in agreement with the existing views. A new rule governing the protonation sequence of linear polyamines is proposed stating that at thermodynamic equilibrium the protonation sequence in a solution results in the formation of tautomers which: (i) minimizes the differences in the atomic energies and net atomic charges of both terminal atoms, N and C, in all partially protonated forms, HL + and H 2 L 2+ , and (ii) have a symmetrical distribution of the atomic energies and net atomic charges on the terminal fragments -CH 2 -NH 2 of a linear molecule. Several conditions, involving the atomic energies and charges of the terminal N-and carbon atoms, are defined and incorporated in a protocol which results in a definite prediction of the protonation sequence of linear aliphatic triamines.