The analysis of Li I 6708A line was performed for 6 rotational phases distributed over the whole rotational period (9.5 years) of HD166473. The magnetic field model was constructed based on the polarimetric measurements from (Mathys et al. 20007). For each observed phase the modulus of the magnetic field was also estimated from simulation of the Fe II 6147A, 6149A and Pr III 6706.7A line profiles taking into account Zeeman magnetic splitting. The lithium abundance in each phase was obtained from fitting the observed Li I 6708A profile with the synthetic one calculated assuming Paschen-Back splitting and estimated magnetic field characteristics from Pr III 6706.7A line profile.

The magnetic star HD 21699 possesses a unique magnetic field structure where the magnetic dipole is displaced from the centre by 0.4 +/- 0.1 of the stellar radius (perpendicularly to the magnetic axis), as a result, the magnetic poles are situated close to one another on the stellar surface with an angular separation of 55$^o$ and not 180$^o$ as seen in the case of a centred dipole. Respectively, the two magnetic poles form a large "magnetic spot". High-resolution spectra were obtained allowing He I and Si II abundance variations to be studied as a function of rotational phase. The results show that the helium abundance is concentrated in one hemisphere of the star, near the magnetic poles and it is comparatively weaker in another hemisphere, where magnetic field lines are horizontal with respect to the stellar surface. At the same time, the silicon abundance is greatest between longitudes of 180 - 320$^o$, the same place where the helium abundance is the weakest. These abundance variations (with rotational phase) support predictions made by the theory of atomic diffusion in the presence of a magnetic field. Simultaneously, these result support the possibility of the formation of unusual structures in stellar magnetic fields. Analysis of vertical stratification of the silicon and helium abundances shows that the boundaries of an abundance jump (in the two step model) are similar for each element; $\tau_{5000}$ = 0.8-1.2 for helium and 0.5-1.3 for silicon. The elemental abundances in the layers of effective formation of selected absorption lines for various phases are also correlated with the excitation energies of low transition levels: abundances are enhanced for higher excitation energy and higher optical depth within the applied model atmosphere.

The surface structure of the magnetic field in stars is one part of the information about the surface inhomogeneities like brightness, movement or chemical composition, which is contained in the integral radiation flux but distorted beyond recognition by losses of information on the topographical arrangement over the surface because of convolution processes and partial invisibility. Therefore, the complicated processes, beginning with the representation of the map of the surface field distribution and ending with the resulting phase-curves of the integral magnetic field, are simulated by a computer programme. We use this programme as a tool for the analysis of the magnetic field structure out of the really observed data in phase diagrams by variation of the parameters and by fitting the calculated curves to the observed ones.

We investigate the radial velocity and the magnetic field of the CP star alpha^2 CVn. The observed variation of the magnetic field is compared with that of our model. We search for a relation between the magnetic field and the distribution of the chemical elements. The period in the radial velocities is constant over a time interval of about 100 years.

We search for a relation between the published distributions of different elements and the calculated magnetic field structure, following from a dipole-quadrupole configuration, of the CP2 star CU Vir. The highest concentration of individual chemical elements on the stellar surface coincides obviously with the regions of the highest values of the magnetic field strength.