The high-velocity molecular jet driven by Class 0 protostar IRAS 04166+2706 exhibits a unique saw-tooth velocity pattern. It consists of a series of well-aligned symmetric knots with similar averaged speeds, whose speeds at peaks of emission decreases roughly linearly away from the origin. Recent ALMA observations of knots R6 and B6 reveal kinematic behavior with expansion velocity increasing linearly from the axis to the edge. This pattern can be formed by a spherically expanding wind with axial density concentration. In this picture, the diverging velocity profile naturally possesses an increasing expansion velocity away from the axis, resulting in a tooth-like feature on the position-velocity diagram through projection. Such geometric picture predicts a correspondence between the slopes of the teeth and the outflow inclination angles, and the same inclination angle of 52$^\circ$ of the IRAS 04166+2706 can generally explain the whole pattern. Aided by numerical simulations in the framework of unified wind model by Shang et al. (2006), the observed velocity pattern can indeed be generated. A proper geometrical distribution of the jet and wind material is essential to the reconstruction the ejection history of the system.

We explore kinematics and morphologies of molecular outflows driven by young protostars using magnetohydrodynamic simulations in the context of the unified wind model of Shang et al. The model explains the observed high-velocity jet and low-velocity shell features. In this work we investigate how these characteristics are affected by the underlying temperature and magnetic field strength. We study the problem of a warm wind running into a cold ambient toroid by using a tracer field that keeps track of the wind material. While an isothermal equation of state is adopted, the effective temperature is determined locally based on the wind mass fraction. In the unified wind model, the density of the wind is cylindrically stratified and highly concentrated toward the outflow axis. Our simulations show that for a sufficiently magnetized wind, the jet identity can be well maintained even at high temperatures. However, for a high temperature wind with low magnetization, the thermal pressure of the wind gas can drive material away from the axis, making the jet less collimated as it propagates. We also study the role of the poloidal magnetic field of the toroid. It is shown that the wind-ambient interface becomes more resistant to corrugation when the poloidal field is present, and the poloidal field that bunches up within the toroid prevents the swept-up material from being compressed into a thin layer. This suggests that the ambient poloidal field may play a role in producing a smoother and thicker swept-up shell structure in the molecular outflow.

The ejection of matter in the close vicinity of a young stellar object is investigated, treating the accretion disk as a gravitationally bound reservoir of matter. By solving the resistive MHD equations in 2D axisymmetry using our version of the Zeus-3D code with newly implemented resistivity, we study the effect of magnetic diffusivity in the magnetospheric accretion-ejection mechanism. Physical resistivity was included in the whole computational domain so that reconnection is enabled by the physical as well as the numerical resistivity. We show, for the first time, that quasi-stationary fast ejecta of matter, which we call {\em micro-ejections}, of small mass and angular momentum fluxes, can be launched from a purely resistive magnetosphere. They are produced by a combination of pressure gradient and magnetic forces, in presence of ongoing magnetic reconnection along the boundary layer between the star and the disk, where a current sheet is formed. Mass flux of micro-ejection increases with increasing magnetic field strength and stellar rotation rate, and is not dependent on the disk to corona density ratio and amount of resistivity.

Stellar magnetosphere and accretion disk interact, and a result should be outflow launched from the innermost vicinity of a protostellar object. We simulated physical conditions in this region by resistive MHD simulations. Outflows resembling the observed ones do not happen in the closest vicinity, except for quasi-stationary funnel flows onto the star, but could occur at few tens of stellar radii above the star. Numerical simulations we performed using our resistive version of ZEUS-3D code, ZEUS347.