The ground bound states in the H2+, D2+, HD+ and
HT+ ions are determined to a high accuracy by using the
multi-box strategy developed (Frolov A M 2001 Phys. Rev.
E 64 036704) for high precision variational calculations
in three-body systems. In this approach the ground state
wavefunctions are represented as the variational sums of generalized
exponents written in the relative inter-particle coordinates
r32,r31 and r21. Our present variational energies
and bound state properties are determined without making use of
the Born-Oppenheimer approximation. The computed expectation
values for various properties in the H2+, D2+, HD+
and HT+ ions are significantly more accurate than values
known for these systems from earlier computations. The observed
agreement between the computed and predicted cusp values is very
good for electron-nuclear cusps. In contrast with this, the
computed nuclear-nuclear cusps and nuclear-nuclear
delta functions differ significantly from the predicted values.