Magnetic moments of charged and neutral mesons are calculated with the use of the relativistic Hamiltonian derived from the path integral form of the $q_1\bar q_2$ Green's function. The magnetic moments are shown to be expressed via the average quark energies which are defined by the fundamental quantities: the string tension $\sigma$, the current quark masses, and the strong coupling constant $\alpha_s$. Resulting values for vector, axial, and tensor light and $K$ mesons agree well with all available lattice data.

Theory of strong decays defines in addition to decay widths, also the channel coupling and the mass shifts of the levels above the decay thresholds. In the standard decay models of the 3P0 type the decay vertex is taken to be a phenomenological constant "gamma" and such a choice leads to large mass shifts of all meson levels due to real and virtual decays, the latter giving a divergent contribution. Here we show that taking the microscopic details of decay vertex into account, one obtains new string width coefficient, which strongly suppresses virtual decay contribution. In addition for a realistic space structure of the decay vertex of highly excited states, the decay matrix elements appear to be strongly different from those, where the constant "gamma" is used. From our analysis also follows that so-called flattening potential can imitate the effects of intermediate decay channels.

Radiative decays of X(3872) with $J^{PC}=1^{++}$ are studied in the coupled-channel approach, where the $c\bar c$ states are described by relativistic string Hamiltonian, while for the decay channels $DD^*$ a string breaking mechanism is used. Within this method a sharp peak and correct mass shift of the $2 {}^3P_1$ charmonium state just to the $D^0D^{*0}$ threshold was already obtained for a prescribed channel coupling to the $DD^*$ decay channels. For the same value of coupling the normalized wave function (w.f.) of X(3872) acquires admixture of the $1 {}^3P_1$ component with the w.f. fraction $c_1=0.153 (\theta=8.8^\circ$), which increases the transition rate $\Gamma(X(3872)\rightarrow J/\psi\gamma)$ up to 50-70 keV, making the ratio $R=\frac{\mathcal{B}(X(3872)\rightarrow \psi^{\prime}\gamma)}{\mathcal{B}(X(3872)\rightarrow J/\psi \gamma)}=0.8\pm 0.20 (th)$ significantly smaller, as compared to $R\simeq 5$ for X(3872) as a purely $2 {}^3P_1$ state.

The masses of $\psi((n+1) {}^3S_1)$ and $\psi(n {}^3D_1)$ are calculated using the relativistic string Hamiltonian with "linear+gluon-exchange" potential. They occur in the range 4.5-5.8 GeV, in particular, $M(3D)=4.54$ GeV, $M(5S)=4.79$ GeV, $M(4D)=4.85$ GeV are calculated with accuracy $\sim 50$ MeV. Linear Regge trajectories: $M^2(nS)=M^2(\psi(4.42))+ 2.91$ GeV$^2 (n-4)$ $(n\geq 4) $ and $M^2(nD)=(4.54^2+ 2.88 (n-3))$ GeV$^2$ ($n\geq 3$) are obtained only for higher charmonium excitations. They have slopes two times larger than those of light mesons and give good description of calculated masses. These masses are compared with enhancements in some recent $e^+e^-$ experiments.

The vector and pseudoscalar decay constants are calculated in the framework of the Field Correlator Method. Di-electron widths: $\Gamma_{ee}(J/\psi)=5.41$ keV, $\Gamma_{ee}(\psi'(3686))=2.47$ keV, $\Gamma_{ee}(\psi''(3770))=0.248$ keV, in good agreement with experiment, are obtained with the same coupling, $\alpha_s=0.165$, in QCD radiative corrections. We show that the larger $\alpha_s=0.191\pm 0.004$ is needed to reach agreement with experiment for $\Gamma_{\gamma\gamma}(\eta_c)=7.22$ keV, $\Gamma_{\gamma\gamma} (\chi(^3P_0))=3.3$ keV, $\Gamma_{\gamma\gamma}(\chi(^3P_2))= 0.54$ keV, and also for $\Gamma(J/\psi\to 3g)=59.5$ keV, $\Gamma(J/\psi\to \gamma 2g)=5.7$ keV. Meanwhile even larger $\alpha_s=0.238$ gives rise to good description of $\Gamma(\psi'\to 3g)=52.7$ keV, $\Gamma(\psi'\to \gamma 2g)= 3.5$ keV, and provides correct ratio of the branching fractions: $\frac{\mathcal{B}(J/\psi\to light hadrons)}{\mathcal{B}(\psi'\to light hadrons)}=0.24.$

The di-electron widths of $\psi(4040)$, $\psi(4160)$, and $\psi(4415)$, and their ratios are shown to be in good agreement with experiment, if in all cases the $S-D$ mixing with a large mixing angle $\theta\approx 34^\circ$ is taken. Arguments are presented why continuum states give small contributions to the wave functions at the origin. We find that the Y(4360) resonance, considered as a pure $3 {}^3D_1$ state, would have very small di-electron width, $\Gamma_{ee}(Y(4360))=0.060$ keV. On the contrary, for large mixing between the $4 {}^3S_1$ and $3 {}^3D_1$ states with the mixing angle $\theta=34.8^\circ$, $\Gamma_{ee}(\psi(4415))=0.57$ keV coincides with the experimental number, while a second physical resonance, probably Y(4360), has also a rather large $\Gamma_{ee} (Y(\sim 4400))=0.61$ keV. For the higher resonance Y(4660), considered as a pure $5 {}^3S_1$ state, we predict the di-electron width $\Gamma_{ee}(Y(4660))=0.70$ keV, but it becomes significantly smaller, namely 0.31 keV, if the mixing angle between the $5 {}^3S_1$ and $4 {}^3D_1$ states $\theta=34^\circ$. The mass and di-electron width of the $6 {}^3S_1$ charmonium state are calculated.

The hyperfine splittings in heavy quarkonia are studied using new experimental data on the di-electron widths. The smearing of the spin-spin interaction is taken into account, while the radius of smearing is fixed by the known $J/\psi-\eta_c(1S)$ and $\psi(2S)-\eta'_c(2S)$ splittings and appears to be small, $r_{ss} \approx 0.06$ fm. Nevertheless, even with such a small radius an essential suppression of the hyperfine splittings ($\sim 50%)$ is observed in bottomonium. For the $nS b\bar b$ states $(n=1,2,...6)$ the values we predict (in MeV) are 28, 12, 10, 6, 6, and 3, respectively. In single-channel approximation for the $3S$ and $4S$ charmonium states the splittings 16(2) MeV and 12(4) MeV are obtained.