Direct variational calculations, where the Born-Oppenheimer (BO) approximation is not assumed, are performed for all 26 bound rovibrational states corresponding to the lowest rotational excitation (i.e. the N = 1 states) of the tritium molecule (T-2). The non-BO energies are used to determine the ortho-para isomerization energies. All-particle explicitly correlated Gaussian basis functions are employed in the calculations and over 11 000 Gaussians independently generated for each state are used. The exponential parameters of the Gaussians are optimized with the aid of analytically calculated energy gradient determined with respect to these parameters. The non-BO wave functions are used to calculate expectation values of the inter-particle distances and the triton-triton correlation functions. (C) 2015 Elsevier B.V. All rights reserved.