Herschel observations have recently revealed that interstellar molecular clouds consist of many filaments. Polarization observations in optical and infrared wavelengths indicate that the magnetic ...field often runs perpendicular to the filament. In this article, we study the magnetohydrostatic configuration of isothermal gas in which the thermal pressure and the Lorentz force are balanced against the self-gravity, and the magnetic field is globally perpendicular to the axis of the filament. The model is controlled by three parameters: center-to-surface density ratio (rho sub(c)/rho sub(s)), plasma beta of surrounding interstellar gas ( beta sub(0)), and the radius of the hypothetical parent cloud normalized by the scale-height (R' sub(0)), although there remains freedom in how the mass is distributed against the magnetic flux (mass loading). In the case where R' sub(0) is small enough, the magnetic field plays a role in confining the gas. However, the magnetic field generally has the effect of supporting the cloud. There is a maximum line-mass (mass per unit length) above which the cloud is not supported against gravity. Compared with the maximum line-mass of a nonmagnetized cloud (2c super(2) s/G, where c sub(s) and G represent, respectively, the isothermal sound speed and the gravitational constant), that of the magnetized filament is larger than the nonmagnetized one. The maximum line-mass is numerically obtained as lambda sub(max) Asymptotically = to 0.24Phi sub(cl)/G super(1 /2) + 1.66c super(2) sub(s)/G, where Phi sub(cl) represents one half of the magnetic flux threading the filament per unit length. The maximum mass of the filamentary cloud is shown to be significantly affected by the magnetic field when the magnetic flux per unit length exceeds Phi sub(cl) gap 3 pc mu G(c sub(s)/1 90 m s super(-1)) super(2).
In this study, we investigate the line emissions from cold molecular gas based on our previous "radiation-driven fountain model," which reliably explains the spectral energy distribution of the ...nearest type 2 Seyfert galaxy, the Circinus galaxy. Using a snapshot of the best-fit radiation-hydrodynamic model for the central pc, in which non-equilibrium X-ray-dominated region chemistry is solved, we conduct post-processed non-local thermodynamic equilibrium radiation transfer simulations for the CO lines. We obtain a spectral line energy distribution with a peak around , and its distribution suggests that the lines are not thermalized. However, for a given line of sight, the optical depth distribution is highly non-uniform between and . The CO-to-H2 conversion factor ( ), which can be directly obtained from the results and is not a constant, depends strongly on the integrated intensity and differs from the fiducial value for local objects. exhibits a large dispersion of more than one order of magnitude, reflecting the non-uniform internal structure of a "torus." In addition, we found that the physical conditions differ between grid cells on a scale of a few parsecs along the observed lines of sight; therefore, a specific observed line ratio does not necessarily represent a single physical state of the interstellar medium.
Abstract
Filamentary structures are recognized as a fundamental component of interstellar molecular clouds in observations made by the Herschel satellite. These filaments, especially massive ...filaments, often extend in a direction perpendicular to the interstellar magnetic field. Furthermore, the filaments sometimes have an apparently negative temperature gradient—that is, their temperatures decrease toward their centers. In this paper, we study the magnetohydrostatic equilibrium state of negative-indexed polytropic gas with the magnetic field running perpendicular to the axis of the filament. The model is controlled by four parameters: center-to-surface density ratio (
ρ
c
/
ρ
s
), plasma
β
of the surrounding gas, radius of the parent cloud
R
0
′
normalized by the scale height, and the polytropic index
N
. The steepness of the temperature gradient is represented by
N
. We found that the envelope of the column density profile becomes shallow when the temperature gradient is large. This reconciles the inconsistency between the observed profiles and those expected from the isothermal models. We compared the maximum line mass (mass per unit length), above which there is no equilibrium, with that of the isothermal nonmagnetized filament. We obtained an empirical formula to express the maximum line mass of a magnetized polytropic filament as
λ
m
a
x
≃
λ
0
,
m
a
x
(
N
)
/
M
⊙
p
c
−
1
2
+
5.9
1.0
+
1.2
/
N
1
/
2
Φ
c
l
/
1
μ
G
p
c
2
1
/
2
M
⊙
pc
−1
, where
λ
0
,
max
(
N
)
represents the maximum line mass of the nonmagnetized filament and Φ
cl
indicates one-half of the magnetic flux threading the filament per unit length. Although the negative-indexed polytrope makes the maximum line mass decrease, compared with that of the isothermal model, a magnetic field threading the filament increases the line mass.
The mechanisms causing millimeter-wave polarization in protoplanetary disks are under debate. To disentangle the polarization mechanisms, we observe the protoplanetary disk around HL Tau at 3.1 mm ...with the Atacama Large Millimeter/submillimeter Array (ALMA), which had the polarization detected with CARMA at 1.3 mm. We successfully detect the ring-like azimuthal polarized emission at 3.1 mm. This indicates that dust grains are aligned with the major axis being in the azimuthal direction, which is consistent with the theory of radiative alignment of elongated dust grains, where the major axis of dust grains is perpendicular to the radiation flux. Furthermore, the morphology of the polarization vectors at 3.1 mm is completely different from those at 1.3 mm. We interpret the polarization at 3.1 mm to be dominated by the grain alignment with the radiative flux producing azimuthal polarization vectors, while the self-scattering dominates at 1.3 mm and produces the polarization vectors parallel to the minor axis of the disk. By modeling the total polarization fraction with a single grain population model, the maximum grain size is constrained to be , which is smaller than the previous predictions based on the spectral index between ALMA at 3 mm and the Very Large Array at 7 mm.
The observational expectation of polarization measurements of thermal dust radiation is investigated to find information on molecular outflows based on magnetohydrodynamical (MHD) and ...radiation-transfer simulations. There are two major proposed models for driving of molecular outflows: (1) molecular gas is accelerated by magnetic pressure gradient or magnetocentrifugal wind mechanism before the magnetic field and molecular gas are decoupled, (2) the linear momentum of a highly collimated jet is transferred to the ambient molecular gas. In order to distinguish between these two models, it is crucial to observe the configuration of the magnetic field. An observation of a toroidal magnetic field would be strong evidence that the first model is appropriate. We calculated the polarization distribution of thermal dust radiation due to the alignment of dust grains along the magnetic field using molecular outflow data obtained from two-dimensional axisymmetric MHD simulations. An asymmetric distribution around the
$z$
-axis is characteristic for magnetic fields composed of both poloidal and toroidal components. We found that the outflow has a low polarization degree compared with the envelope and that the envelope and outflow have different polarization directions (B-vector); i.e., the magnetic field within the envelope is parallel to the global magnetic field lines while the magnetic field of the outflow is perpendicular to it. We, then, demonstrated that the point-symmetric (rather than axisymmetric) distributions of low polarization regions indicate that molecular outflows are likely to be magnetically driven. Observations of this polarization distribution with tools such as ALMA would confirm the origin of the molecular outflow.
Abstract
Filamentary molecular clouds are regarded as the place where newborn stars form. In particular, a hub region, a place where it appears as if several filaments are colliding, often indicates ...active star formation. To understand the star formation in filament structures, we investigate the collisions between two filaments using two-dimensional magnetohydrodynamical simulations. As a model of filaments, we assume that the filaments are in magnetohydrostatic equilibrium under a global magnetic field perpendicular to the filament axis. We set two identical filaments with an infinite length and make them collide with a zero-impact parameter (head-on). When the two filaments collide while sharing the same magnetic flux, we found two types of evolution after a merged filament is formed: runaway radial collapse and stable oscillation with a finite amplitude. The condition for the radial collapse is independent of the collision velocity and is given by the total line mass of the two filaments exceeding the magnetically critical line mass for which no magnetohydrostatic solution exists. The radial collapse proceeds in a self-similar manner, resulting in a unique distribution irrespective of the various initial line masses of the filament, as the collapse progresses. When the total line mass is less massive than the magnetically critical line mass, the merged filament oscillates, and the density distribution is well-fitted by a magnetohydrostatic equilibrium solution. The condition necessary for the radial collapse is also applicable to the collision whose direction is perpendicular to the global magnetic field.
We examine the linear stability of an isothermal filamentary cloud permeated by a perpendicular magnetic field. Our model cloud is assumed to be supported by gas pressure against self-gravity in the ...unperturbed state. For simplicity, the density distribution is assumed to be symmetric around the axis. Also for simplicity, the initial magnetic field is assumed to be uniform, and turbulence is not taken into account. The perturbation equation is formulated to be an eigenvalue problem. The growth rate is obtained as a function of the wavenumber for fragmentation along the axis and the magnetic field strength. The growth rate depends critically on the outer boundary. If the displacement vanishes in regions very far from the cloud axis (fixed boundary), cloud fragmentation is suppressed by a moderate magnetic field, which means the plasma beta is below 1.67 on the cloud axis. If the displacement is constant along the magnetic field in regions very far from the cloud, the cloud is unstable even when the magnetic field is infinitely strong. The cloud is deformed by circulation in the plane perpendicular to the magnetic field. The unstable mode is not likely to induce dynamical collapse, since it is excited even when the whole cloud is magnetically subcritical. For both boundary conditions, the magnetic field increases the wavelength of the most unstable mode. We find that the magnetic force suppresses compression perpendicular to the magnetic field especially in regions of low density.
We report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation. We take into account many physical processes ...required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail. We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.