Three-dimensional simulations of disk accretion to an inclined dipole.: I.: Magnetospheric flows at different Θ

We present results of fully three-dimensional MHD simulations of disk accretion to a slowly rotating magnetized star with its dipole moment inclined at an angle Theta to the rotation axis of the disk (which is aligned with the spin axis of the star). The main goal was to investigate the pattern of magnetospheric flow and the disk-star interaction for a variety of inclination angles Theta. We observed that at Theta = 0degrees, the disk stops at magnetospheric radius r(m), and matter flows to the star through axisymmetric funnel flows, as observed in earlier axisymmetric simulations. However, when the dipole moment of the star is inclined, then the flow becomes nonaxisymmetric. The nonaxisymmetry becomes notable at very small inclination angles, Theta similar to 2degrees-5degrees. The pattern of magnetospheric flow is different at different Theta. For relatively small angles, Theta less than or similar to 30degrees, the densest matter flows to the star mostly in two streams, which follow paths to the closest magnetic pole. The streams typically corotate with the star, but they may precess about the star for Theta less than or similar to 10degrees. At intermediate angles, 30degrees less than or similar to Theta less than or similar to 60degrees, the streams become more complicated and often split into four streams. For even larger angles, Theta greater than or similar to 60degrees, matter accretes in two streams, but their geometry is different from that of the streams at small Theta. Magnetic braking changes the structure of the inner regions of the disk. It creates a region of lower density ("gap") for r(m) less than or similar to r less than or similar to 4r(m). A ring of higher density forms at r similar to r(m) for Theta less than or similar to 30degrees. For r less than or similar to (2-3)r(m), the azimuthal velocities are sub-Keplerian. The inner region of the disk at r similar to r(m) is warped. The warping is due to the tendency of matter to corotate with an inclined magnetosphere. The normal of the inner warped part of the disk is close to the magnetic axis of the dipole. The accreting matter brings positive angular momentum to the ( slowly rotating) star, tending to spin it up. The corresponding torque NZ depends only weakly on Theta. The angular momentum flux to the star near the star's surface is transported predominantly by the magnetic field; the matter component contributes similar to1% of the total flux. The torques N-X and N-Y are also calculated, and these may give a slow precession of the symmetry axis of the star. The angle Theta was fixed in simulations because the timescale of its evolution is much longer than that of the simulations. Results of simulations are important for understanding the nature of classical T Tauri stars, cataclysmic variables, and X-ray pulsars. These stars often show complicated spectral and photometric variability patterns, which may be connected with the structure of magnetospheric flows. The magnetospheric structure of stars with different Theta can give different variability patterns in observed light curves. This can provide information about inclination angles Theta in different stars. A notable result of the present simulations is the formation of multiple streams in the accretion flows near the star for intermediate inclination angles. This may give short-scale, quasi-periodic variability in the light curves of some stars.