Mass spectrum of heavy tetraquarks in variational approach

Within the framework of the quark model and the variational method, the bound states of four heavy quarks (tetraquarks) are investigated. The basis variational wave functions are chosen in the Gaussian form. The matrix elements of the Hamiltonian are calculated analytically. Numerical calculation of the energy levels of tetraquarks is performed using a program in the Matlab system. Hyperfine structure of the spectrum is calculated. To increase the accuracy of the calculation, relativistic corrections are taken into account.