Abstract
We describe a numerical scheme for magnetohydrodynamics simulations of dust–gas mixture by extending smoothed particle magnetohydrodynamics. We employ the single-species particle approach to ...describe dust–gas mixture with several modifications from the previous studies. We assume that the charged and neutral dust can be treated as single-fluid, that the electromagnetic force acts on the gas, and that that acting on the charged dust is negligible. The validity of these assumptions in the context of protostar formation is not obvious and is extensively evaluated. By investigating the electromagnetic force and electric current with terminal velocity approximation, it is found that as the dust size increases, the contribution of dust to them becomes smaller and negligible. We conclude that our assumption that the electromagnetic force on the dusts is negligible is valid for the dust size with
a
d
≳ 10
μ
m. On the other hand, they do not produce the numerical artifact for the dust
a
d
≲ 10
μ
m in the envelope and disk, where the perfect coupling between gas and dust is realized. However, we also found that our assumptions may break down in outflow (or under an environment with very strong magnetic field and low density) for the dust
a
d
≲ 10
μ
m. We conclude that our assumptions are valid in almost all cases where macroscopic dust dynamics is important in the context of protostar formation. We conduct numerical tests of dusty waves, dusty magnetohydrodynamics shocks, and gravitational collapse of magnetized cloud cores with our simulation code. The results show that our numerical scheme well reproduces the dust dynamics in the magnetized medium.
The formation and early evolution of low-mass young stellar objects (YSOs) are investigated using three-dimensional non-ideal magnetohydrodynamics simulations. We investigate the evolution of YSOs up ...to after protostar formation, at which protostellar mass reaches . We particularly focus on the impact of the dust model on the evolution. We found that a circumstellar disk is formed in all simulations, regardless of the dust model. Disk size is approximately 10 au at the protostar formation epoch, and it increases to several tens of au at after protostar formation. The disk mass is comparable to the central protostellar mass, and gravitational instability develops. In simulations with small dust sizes, the warp of the pseudodisk develops after protostar formation. The warp strengthens magnetic braking in the disk and decreases disk size. Ion-neutral drift can occur in the infalling envelope when the typical dust size is and the protostar (plus disk) mass is . The outflow activity is anticorrelated to the dust size, and the strong outflow appears with small dust grains.
The effect of misalignment between the magnetic field and the angular momentum of molecular cloud cores on the angular momentum evolution during the gravitational collapse is investigated by ideal ...and non-ideal MHD simulations. For the non-ideal effect, we consider the ohmic and ambipolar diffusion. Previous studies that considered the misalignment reported qualitatively contradicting results. Magnetic braking was reported as being either strengthened or weakened by misalignment in different studies. We conducted simulations of cloud core collapse by varying the stability parameter (the ratio of the thermal to gravitational energy of the core) with and without including magnetic diffusion. The non-ideal MHD simulations show the central angular momentum of the core, with θ = 0° ( ) being always greater than that with θ = 90° ( ), independently of , meaning that circumstellar disks form more easily in a core with θ = 0°. The ideal MHD simulations, in contrast, show the central angular momentum of the core with θ = 90° being greater than with θ = 0° for small and smaller for large . Inspection of the angular momentum evolution of the fluid elements reveals three mechanisms contributing to the evolution of the angular momentum: (i) magnetic braking in the isothermal collapse phase, (ii) selective accretion of the rapidly (for θ = 90°) or slowly (for θ = 0°) rotating fluid elements to the central region, and (iii) magnetic braking in the first core and the disk. The difference between the ideal and non-ideal simulations arises from the different efficiencies of (iii).
Fragmentation of protoplanetary discs due to gravitational instabilities is a candidate of a formation mechanism of binary stars, brown dwarfs, and gaseous giant planets. The condition for the ...fragmentation has been thought that the disc cooling time-scale is comparable to its dynamical time-scale. However, some numerical simulations suggest that the fragmentation does not occur even if the cooling time is small enough, or the fragmentation can occur even when the cooling is inefficient. To reveal a realistic condition for fragmentation of self-gravitating discs, we perform two-dimensional numerical simulations that take into account the effect of the irradiation of the central star and radiation cooling of the disc, and precisely investigate the structure of the spiral arms formed in the protoplanetary discs. We show that the Toomre Q parameter in the spiral arms is an essential parameter for fragmentation. The spiral arms fragment only when Q < 0.6 in the spiral arms. We have further confirmed that this fragmentation condition observed in the numerical simulations can be obtained from the linear stability analysis for the self-gravitating spiral arms. These results indicate that the process of fragmentation of protoplanetary discs is divided into two stages: formation of the spiral arms in the discs; and fragmentation of the spiral arm. Our work reduces the condition for the fragmentation of the protoplanetary discs to the condition of the formation of the spiral arm that satisfies Q < 0.6.
We investigate the formation and evolution of a first core, protostar, and circumstellar disc with a three-dimensional non-ideal (including both Ohmic and ambipolar diffusion) radiation ...magnetohydrodynamics simulation. We found that the magnetic flux is largely removed by magnetic diffusion in the first-core phase and that the plasma β of the centre of the first core becomes large, β > 104. Thus, proper treatment of first-core phase is crucial in investigating the formation of protostar and disc. On the other hand, in an ideal simulation, β ∼ 10 at the centre of the first core. The simulations with magnetic diffusion show that the circumstellar disc forms at almost the same time of protostar formation even with a relatively strong initial magnetic field (the value for the initial mass-to-flux ratio of the cloud core relative to the critical value is μ = 4). The disc has a radius of r ∼ 1 AU at the protostar formation epoch. We confirm that the disc is rotationally supported. We also show that the disc is massive (Q ∼ 1) and that gravitational instability may play an important role in the subsequent disc evolution.
By performing a one-dimensional magnetohydrodynamic simulation with radiative cooling and thermal conduction, we show that the coronal heating and the fast solar wind acceleration in the coronal ...holes are natural consequences of the footpoint fluctuations of the magnetic fields at the photosphere. We initially set up a static open flux tube with a temperature of 10 super(4) K rooted at the photosphere. We impose transverse photospheric motions corresponding to the granulations with a velocity < dv sub( )> = 0.7 km s super(-1) and a period between 20 s and 30 minutes, which generate outgoing Alfven waves. We self-consistently treat these waves and the plasma heating. After attenuation in the chromosphere by 85% of the initial energy flux, the outgoing AlfVen waves enter the corona and contribute to the heating and acceleration of the plasma mainly by the nonlinear generation of the compressive waves and shocks. Our result clearly shows that the initially cool and static atmosphere is naturally heated up to 10 super(6) K and accelerated to 800 km s super(-1).
We investigate the structure of self-gravitating discs, their fragmentation and the evolution of the fragments (the clumps) using both an analytic approach and three-dimensional radiation ...hydrodynamics simulations starting from molecular cores. The simulations show that non-local radiative transfer determines the disc temperature. We find the disc structure is well described by an analytical model of a quasi-steady self-gravitating disc with radial radiative transfer. Because the radiative process is not local and radiation from the interstellar medium cannot be ignored, the local radiative cooling is not balanced with the viscous heating in a massive disc around a low-mass star. In our simulations, there are cases in which the disc does not fragment even though it satisfies the fragmentation criterion based on disc cooling time (...). This indicates that, at least, the criterion is not a sufficient condition for fragmentation. We determine the parameter range for the host cloud core in which disc fragmentation occurs. In addition, we show that the temperature evolution of the centre of the clump is close to that of typical first cores, and that the minimum initial mass of clumps is about a few Jupiter masses. (ProQuest: ... denotes formulae/symbols omitted.)
More than a half of the asteroids in the main belt have irregular shapes with ratios of the minor to major axis lengths of less than 0.6. One of the mechanisms that create such shapes is collisions ...between asteroids. The relationship between the shapes of collisional outcomes and impact conditions such as impact velocities may provide information on the collisional environments and its evolutionary stages when those asteroids are created. In this study, we perform numerical simulations of collisional destruction of asteroids with radii 50 km and subsequent gravitational reaccumulation using smoothed-particle hydrodynamics for elastic dynamics with self-gravity, a model of rock fractures, and a model of friction in completely damaged rock. We systematically vary the impact velocity from 50 to 400 m s−1 and the impact angle from 5° to 45°. We investigate shapes of the largest remnants resulting from collisional simulations. As a result, various shapes (bilobed, spherical, flat, elongated, and hemispherical shapes) are formed through equal-mass and low-velocity (50−400 m s−1) impacts. We clarify a range of the impact angle and velocity to form each shape. Our results indicate that irregular shapes, especially flat shapes, of asteroids with diameters larger than 80 km are likely to be formed through similar-mass and low-velocity impacts, which are likely to occur in the primordial environment prior to the formation of Jupiter.
Smoothed particle hydrodynamics is reformulated in terms of the convolution of the original hydrodynamics equations, and the new evolution equations for the particles are derived. The same evolution ...equation of motion is also derived using a new action principle. The force acting on each particle is determined by solving the Riemann problem. The use of the Riemann solver strengthens the method, making it accurate for the study of phenomena with strong shocks. The prescription for the variable smoothing length is shown. These techniques are implemented in strict conservation form. The results of a few test problems are also shown.
Context.
Despite recent observational and theoretical advances in mapping the magnetic fields associated with molecular clouds, their three-dimensional (3D) morphology remains unresolved. ...Multi-wavelength and multi-scale observations will allow us to paint a comprehensive picture of the magnetic fields of these star-forming regions.
Aims.
We reconstructed the 3D magnetic field morphology associated with the Perseus molecular cloud and compared it with predictions of cloud-formation models. These cloud-formation models predict a bending of magnetic fields associated with filamentary molecular clouds. We compared the orientation and direction of this field bending with our 3D magnetic-field view of the Perseus cloud.
Methods.
We used previous line-of-sight and plane-of-sky magnetic field observations as well as Galactic magnetic field models to reconstruct the complete 3D magnetic field vectors and morphology associated with the Perseus cloud.
Results.
We approximated the 3D magnetic field morphology of the cloud as a concave arc that points in the decreasing longitude direction in the plane of the sky (from our point of view). This field morphology preserves a memory of the Galactic magnetic field. In order to compare this morphology to cloud-formation model predictions, we assume that the cloud retains a memory of its most recent interaction. After incorporating velocity observations, we find that the line-of-sight magnetic field observations are consistent with predictions of shock-cloud-interaction models.
Conclusions.
To our knowledge, this is the first time that the 3D magnetic fields of a molecular cloud have been reconstructed. We find the 3D magnetic field morphology of the Perseus cloud to be consistent with the predictions of the shock-cloud-interaction model that describes the formation mechanism of filamentary molecular clouds.