With the help of 3D MHD simulations we investigate the collapse and fragmentation of rotating turbulent prestellar core embedded into turbulent medium. The numerical code is based on a high resolution Godunov-type finite-difference scheme. Initial turbulence is represented by the ensemble of Alfven waves with power law spectrum. Our computations show that under realistic parameters two bound fragments can appear when the density increases at $100-10^3$ times. The distance between the fragments is about $0.1$ of the initial core radius and their orbital period is comparable to the initial free fall time of the core. These results can explain the origin of binary stars with separation $0.001-0.01$ pc in the Galaxy field.

A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gasdynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined during the calculations of the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the radiative transfer in continuum and in molecular lines. We presented the results for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud starting from a quasi-static state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. We have shown that the influence of the magnetic field and rotation can lead to a redistribution of angular momentum in the cloud and the formation of a characteristic rotational velocity structure. As a result, the distribution of the velocity centroid of the molecular lines can acquire an hourglass shape. We plan to use the developed program package together with a model for the chemical evolution to interpret and model observed starless and protostellar cores.