
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 socalled flattening potential can
imitate the effects of intermediate decay channels.

Radiative decays of X(3872) with $J^{PC}=1^{++}$ are studied in the
coupledchannel 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 5070 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+gluonexchange"
potential. They occur in the range 4.55.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 (n4)$
$(n\geq 4) $ and $M^2(nD)=(4.54^2+ 2.88 (n3))$ 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. Dielectron 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 dielectron 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 $SD$ 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 dielectron 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 dielectron 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 dielectron width of the $6 {}^3S_1$ charmonium
state are calculated.

The hyperfine splittings in heavy quarkonia are studied using new
experimental data on the dielectron widths. The smearing of the spinspin
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 singlechannel
approximation for the $3S$ and $4S$ charmonium states the splittings 16(2) MeV
and 12(4) MeV are obtained.