We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We present the first sub-arcsecond images of 49 Ceti in the [C I]3P1-3P0 emission and the 614 micron dust continuum emission observed with ALMA, as well as that in the CO(J=3-2) emission prepared by using the ALMA archival data. The spatial distribution of the 614 micron dust continuum emission is found to have a broad-ring structure with a radius of about 100 au around the central star. A substantial amount of gas is also associated with 49 Ceti. The [C I] emission map shows two peaks inside the dust ring, and its overall extent is comparable to that of the dust continuum emission and the CO emission. We find that the [C I]/CO(J=3-2) intensity ratio significantly varies along the major axis. The ratio takes the minimum value of 1.8 around the dust peak position, and increases inwards and outwards. The enhanced ratio around the central star (~3) likely originates from the stellar UV radiation, while that in the outer disk (~10) from the interstellar UV radiation. Such complex distributions of the [C I] and CO(J=3-2) emission will be a key to understand the origin of the gas in 49 Ceti, and will also provide a stringent constraint on physical and chemical models of gaseous debris disks.

Protostellar flares are rapid magnetic energy release events associated with formation of hot plasma in protostars. In the previous models of protostellar flares, the interaction between a protostellar magnetosphere with the surrounding disk plays crucial roles in building-up and releasing the magnetic energy. However, it remains unclear if protostars indeed have magnetospheres because vigorous disk accretion and strong disk magnetic fields in the protostellar phase may destroy the magnetosphere. Considering this possibility, we investigate the energy accumulation and release processes in the absence of a magnetosphere using a three-dimensional magnetohydrodynamic simulation. Our simulation reveals that protostellar flares are repeatedly produced even in such a case. Unlike in the magnetospheric models, the protostar accumulates magnetic energy by acquiring large-scale magnetic fields from the disk by accretion. Protostellar flares occur when a portion of the large-scale magnetic fields are removed from the protostar as a result of magnetic reconnection. Protostellar flares in the simulation are consistent with observations; the released magnetic energy (up to $\sim 3\times 10^{38}$ erg) is large enough to drive observed flares, and the flares produce hot ejecta. The expelled magnetic fields enhance accretion, and the energy build-up and release processes are repeated as a result. The magnetic flux removal via reconnection leads to redistribution of magnetic fields in the inner disk. We therefore consider that protostellar flares will play an important role in the evolution of the disk magnetic fields in the vicinity of protostars.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We investigate the formation of molecular clouds from atomic gas by using three-dimensional magnetohydrodynamic simulations, including non-equilibrium chemical reactions and heating/cooling processes. We consider super-Alfv\'enic head-on colliding flows of atomic gas possessing the two-phase structure that consists of HI clouds and surrounding warm diffuse gas. We examine how the formation of molecular clouds depends on the angle $\theta$ between the upstream flow and the mean magnetic field. We find that there is a critical angle $\theta_\mathrm{cr}$ above which the shock-amplified magnetic field controls the post-shock gas dynamics. If the atomic gas is compressed almost along the mean magnetic field ($\theta\ll\theta_\mathrm{cr}$), super-Alfv\'enic anisotropic turbulence is maintained by the accretion of the highly inhomogeneous upstream atomic gas. As a result, a greatly extended turbulence-dominated post-shock layer is generated. Around $\theta\sim \theta_\mathrm{cr}$, the shock-amplified magnetic field weakens the post-shock turbulence, leading to a dense post-shock layer. For $\theta\gg \theta_\mathrm{cr}$, the strong magnetic pressure suppresses the formation of cold dense clouds. Efficient molecular cloud formation is expected if $\theta$ is less than a few times $\theta_\mathrm{cr}$. Developing an analytic model and performing a parameter survey, we obtain an analytic formula for the critical angle as a function of the mean density, collision speed, and field strength of the upstream atomic gas. The critical angle is found to be less than $\sim 15^\circ$ as long as the field strength is larger than $1~\mu$G, indicating that the probability of occurrence of compression with $\theta<\theta_\mathrm{cr}$ is limited if shock waves come from various directions.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present the results of a global, three-dimensional magnetohydrodynamics simulation of an accretion disk with a rotating, weakly magnetized central star. The disk is threaded by a weak, large-scale poloidal magnetic field, and the central star has no strong stellar magnetosphere initially. Our simulation investigates the structure of the accretion flows from a turbulent accretion disk onto the star. The simulation reveals that fast accretion onto the star at high latitudes occurs even without a stellar magnetosphere. We find that the failed disk wind becomes the fast, high-latitude accretion as a result of angular momentum exchange mediated by magnetic fields well above the disk, where the Lorentz force that decelerates the rotational motion of gas can be comparable to the centrifugal force. Unlike the classical magnetospheric accretion scenario, fast accretion streams are not guided by magnetic fields of the stellar magnetosphere. Nevertheless, the accretion velocity reaches the free-fall velocity at the stellar surface due to the efficient angular momentum loss at a distant place from the star. This study provides a possible explanation why Herbig Ae/Be stars whose magnetic fields are generally not strong enough to form magnetospheres also show indications of fast accretion. A magnetically driven jet is not formed from the disk in our model. The differential rotation cannot generate sufficiently strong magnetic fields for the jet acceleration because the Parker instability interrupts the field amplification.

The effect of misalignment between the magnetic field $\magB$ and the angular momentum $\Jang$ 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 $\alpha$ (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 $\theta=0^\circ$ ($\Jang \parallel \magB$) being always greater than that with $\theta=90^\circ$ ($\Jang \perp \magB$), independently of $\alpha$, meaning that circumstellar disks form more easily form in a core with $\theta=0^\circ$. The ideal MHD simulations, in contrast, show the the central angular momentum of the core with $\theta=90^\circ$ being greater than with $\theta=0^\circ$ for small $\alpha$, and is smaller for large $\alpha$. 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 $\theta=90^\circ$ ) or slowly (for $\theta=0^\circ$) 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).

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

We perform three-dimensional radiation non-ideal magnetohydrodynamics simulations and investigate the impact of the Hall effect on the angular momentum evolution in the collapsing cloud cores in which the magnetic field $\magB$ and angular momentum $\angJ$ are misaligned with each other. We find that the Hall effect notably changes the magnetic torques in the pseudo-disk, and strengthens and weakens the magnetic braking in cores with an acute and obtuse relative angles between $\magB$ and $\angJ$, respectively. This suggests that the bimodal evolution of the disk size may occur in early disk evolutionary phase even if $\magB$ and $\angJ$ are randomly distributed. We show that a counter-rotating envelope form in the upper envelope of the pseudo-disk in cloud cores with obtuse relative angles. We also find that a counter-rotating region forms at the midplane of the pseudo-disk in cloud cores with acute relative angles. The former and latter types of counter-rotating envelopes may be associated with the YSOs with a large ($r\sim100$ AU) and small ($r\lesssim10$ AU) disks, respectively.

Recent observations suggest that intensive molecular cloud collision can trigger massive star/cluster formation. The most important physical process caused by the collision is a shock compression. In this paper, the influence of a shock wave on the evolution of a molecular cloud is studied numerically by using isothermal magnetohydrodynamics (MHD) simulations with the effect of self-gravity. Adaptive-mesh-refinement and sink particle techniques are used to follow long-time evolution of the shocked cloud. We find that the shock compression of turbulent inhomogeneous molecular cloud creates massive filaments, which lie perpendicularly to the background magnetic field as we have pointed out in a previous paper. The massive filament shows global collapse along the filament, which feeds a sink particle located at the collapse center. We observe high accretion rate dot{M}_acc > 10^{-4} M_sun/yr that is high enough to allow the formation of even O-type stars. The most massive sink particle achieves M>50 M_sun in a few times 10^5 yr after the onset of the filament collapse.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We have detected [C I] 3P1-3P0 emissions in the gaseous debris disks of 49 Ceti and Beta Pictoris with the 10 m telescope of the Atacama Submillimeter Telescope Experiment, which is the first detection of such emissions. The line profiles of [C I] are found to resemble those of CO(J=3-2) observed with the same telescope and the Atacama Large Millimeter/submillimeter Array. This result suggests that atomic carbon (C) coexists with CO in the debris disks, and is likely formed by the photodissociation of CO. Assuming an optically thin [C I] emission with the excitation temperature ranging from 30 to 100 K, the column density of C is evaluated to be (2.2+-0.2)x10^17 and (2.5+-0.7)x10^16 cm^-2 for 49 Ceti and Beta Pictoris, respectively. The C/CO column density ratio is thus derived to be 54+-19 and 69+-42 for 49 Ceti and Beta Pictoris, respectively. These ratios are higher than those of molecular clouds and diffuse clouds by an order of magnitude. The unusually high ratios of C to CO are likely attributed to a lack of H2 molecules needed to reproduce CO molecules efficiently from C. This result implies a small number of H2 molecules in the gas disk; i.e., there is an appreciable contribution of secondary gas from dust grains.

In this study, we investigate the dispersive properties of smoothed particle magnetohydrodynamics (SPM) in a strongly magnetized medium by using linear analysis. In modern SPM, a correction term proportional to the divergence of the magnetic fields is subtracted from the equation of motion to avoid a numerical instability arising in a strongly magnetized medium. From the linear analysis, it is found that SPM with the correction term suffer from significant dispersive errors, especially for slow waves propagating along magnetic fields. The phase velocity for all wave numbers is significantly larger than the exact solution and has a peak at a finite wavenumber. These excessively large dispersive errors occur because magnetic fields contribute an unphysical repulsive force along magnetic fields. The dispersive errors cannot be reduced, even with a larger smoothing length and smoother kernel functions such as the Gaussian or quintic spline kernels. We perform the linear analysis for this problem and find that the dispersive errors can be removed completely while keeping SPM stable if the correction term is reduced by half. These findings are confirmed by several simple numerical experiments.

The formation process of circumstellar disks is still controversial because of the interplay of complex physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we investigate the effect of the Hall current term on the formation of the circumstellar disk using three- dimensional simulations. In our simulations, all non-ideal effects as well as the radiation transfer are considered. The size of the disk is significantly affected by a simple difference in the inherent properties of the prestellar core, namely whether the rotation vector and the magnetic field are parallel or anti- parallel. In the former case, only a very small disk (< 1 AU) is formed. On the other hand, in the latter case, a massive and large (> 20 AU) disk is formed in the early phase of protostar formation. Since the parallel and anti-parallel properties do not readily change, we expect that the parallel and anti-parallel properties are also important in the subsequent disk evolution and the difference between the two cases is maintained or enhanced. This result suggests that the disk size distribution of the Class 0 young stellar objects is bimodal. Thus, the disk evolution can be categorized into two cases and we may call the parallel and anti-parallel systems as Ortho-disk and Para-disk, respectively. We also show that the anti-rotating envelopes against the disk-rotation appear with a size of 200 AU. We predict that the anti-rotating envelope will be found in the future observations.

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 $\beta$ of the centre of the first core becomes large, $\beta>10^4$. 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, $\beta\sim 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 $\mu=4$). The disc has a radius of $r \sim 1$ AU at the protostar formation epoch. We confirm that the disc is rotationally supported. We also show that the disc is massive ($Q\sim 1$) and that gravitational instability may play an important role in the subsequent disc evolution.

We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two-fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells or bubbles on a galactic scale. First we estimate the ensemble-averaged growth rate of molecular clouds over a timescale larger than a million years. Next we perform radiation hydrodynamics simulations to evaluate the destruction rate of magnetized molecular clouds by the stellar FUV radiation. We also investigate the resultant star formation efficiency within a cloud which amounts to a low value (a few percent) if we adopt the power-law exponent -2.5 for the mass distribution of stars in the cloud. We finally describe the time evolution of the mass function of molecular clouds over a long timescale (>1Myr) and discuss the steady state exponent of the power-law slope in various environments.

In this paper, the nonlinear evolution of a bistable interstellar medium is investigated using two-dimensional simulations with a realistic cooling rate, thermal conduction, and physical viscosity. The calculations are performed using periodic boundary conditions without any external dynamical forcing. As the initial condition, a spatially uniform unstable gas under thermal equilibrium is considered. At the initial stage, the unstable gas quickly segregates into two phases, or cold neutral medium (CNM) and warm neutral medium (WNM). Then, self-sustained turbulence with velocity dispersion of $0.1-0.2\;\mathrm{km\;s^{-1}}$ is observed in which the CNM moves around in the WNM. We find that the interfacial medium (IFM) between the CNM and WNM plays an important role in sustaining the turbulence. The self-sustaining mechanism can be divided into two steps. First, thermal conduction drives fast flows streaming into concave CNM surfaces towards the WNM. The kinetic energy of the fast flows in the IFM is incorporated into that of the CNM through the phase transition. Second, turbulence inside the CNM deforms interfaces and forms other concave CNM surfaces, leading to fast flows in the IFM. This drives the first step again and a cycle is established by which turbulent motions are self-sustained.

In this paper, we present an explicit scheme for Ohmic dissipation with smoothed particle magneto-hydrodynamics (SPMHD). We propose a SPH discretization of Ohmic dissipation and solve Ohmic dissipation part of induction equation with the supertime-stepping method (STS) which allows us to take a longer time-step than Courant-Friedrich-Levy stability condition. Our scheme is second-order accurate in space and first-order accurate in time. Our numerical experiments show that optimal choice of the parameters of STS for Ohmic dissipation of SPMHD is {\nu}sts ~ 0.01 and Nsts ~ 5.

In this paper, we investigate the structure of condensation fronts from warm diffuse gas to cold neutral medium (CNM) under the plane parallel geometry. The solutions have two parameters, the pressure of the CNM and the mass flux across the transition front, and their ranges are much wider than previously thought. First, we consider the pressure range where the three phases, the CNM, the unstable phase, and the warm neutral medium, can coexist in the pressure equilibrium. In a wide range of the mass flux, we find solutions connecting the CNM and the unstable phase. Moreover, we find solutions in larger pressure range where there is only one thermal equilibrium state or the CNM. These solutions can be realized in shock-compressed regions that are promising sites of molecular cloud formation. We also find remarkable properties in our solutions. Heat conduction becomes less important with increasing mass flux, and the thickness of the transition layer is characterized by the cooling length instead of the Field length.

In this paper, we develop a new method for magnetohydrodynamics (MHD) using smoothed particle hydrodynamics (SPH). To describe MHD shocks accurately, the Godunov method is applied to SPH instead of artificial dissipation terms. In the interaction between particles, we solve a nonlinear Riemann problem with magnetic pressure for compressive waves and apply the method of characteristics for Alfv{\'e}n waves. An extensive series of MHD test calculations is performed. In all test calculations, we compare the results of our SPH code with those of a finite-volume method with an approximate Riemann solver, and confirm excellent agreement.

We study the escape of cosmic-ray (CR) protons accelerated at a supernova remnant (SNR) by numerically solving a diffusion-convection equation from the vicinity of the shock front to the region far away from the front. We consider the amplifications of Alfven waves generated by the escaping CR particles and their effects on CR escape into interstellar medium (ISM). We find that the amplification of the waves significantly delays the escape of the particles even far away from the shock front (on a scale of the SNR). This means that the energy spectrum of CR particles measured through gamma-ray observations at molecular clouds around SNRs is seriously affected by the particle scattering by the waves.

We perform a linear perturbation analysis of expanding shells driven by expansions of HII regions. The ambient gas is assumed to be uniform. As an unperturbed state, we develop a semi-analytic method for deriving the time evolution of the density profile across the thickness. It is found that the time evolution of the density profile can be divided into three evolutionary phases, deceleration-dominated, intermediate, and self-gravity-dominated phases. The density peak moves relatively from the shock front to the contact discontinuity as the shell expands. We perform a linear analysis taking into account the asymmetric density profile obtained by the semi-analytic method, and imposing the boundary conditions for the shock front and the contact discontinuity while the evolutionary effect of the shell is neglected. It is found that the growth rate is enhanced compared with the previous studies based on the thin-shell approximation. This is due to the boundary effect of the contact discontinuity and asymmetric density profile that were not taken into account in previous works.

We investigate the gravitational fragmentation of expanding shells driven by HII regions using the three-dimensional Lagrangian simulation codes based on the Riemann solver, called Godunov smoothed particle hydrodynamics. The ambient gas is assumed to be uniform. In order to attain high resolution to resolve the geometrically thin dense shell, we calculate not the whole but a part of the shell. We find that perturbations begin to grow earlier than the prediction of the linear analysis under the thin-shell approximation. The wavenumber of the most unstable mode is larger than that in the thin-shell linear analysis. The development of the gravitational instability is accompanied by the significant deformation of the contact discontinuity. These results are consistent with a linear analysis presented by Iwasaki et al. (2011) that have taken into account the density profile across the thickness and approximate shock and contact discontinuity boundary conditions. We derive useful analytic formulae for the fragment scale and the epoch when the gravitational instability begins to grow.

In this paper, the stability of a dynamically condensing radiative gas layer is investigated by linear analysis. Our own time-dependent, self-similar solutions describing a dynamical condensing radiative gas layer are used as an unperturbed state. We consider perturbations that are both perpendicular and parallel to the direction of condensation. The transverse wave number of the perturbation is defined by $k$. For $k=0$, it is found that the condensing gas layer is unstable. However, the growth rate is too low to become nonlinear during dynamical condensation. For $k\ne0$, in general, perturbation equations for constant wave number cannot be reduced to an eigenvalue problem due to the unsteady unperturbed state. Therefore, direct numerical integration of the perturbation equations is performed. For comparison, an eigenvalue problem neglecting the time evolution of the unperturbed state is also solved and both results agree well. The gas layer is unstable for all wave numbers, and the growth rate depends a little on wave number. The behaviour of the perturbation is specified by $kL_\mathrm{cool}$ at the centre, where the cooling length, $L_\mathrm{cool}$, represents the length that a sound wave can travel during the cooling time. For $kL_\mathrm{cool}\gg1$, the perturbation grows isobarically. For $kL_\mathrm{cool}\ll1$, the perturbation grows because each part has a different collapse time without interaction. Since the growth rate is sufficiently high, it is not long before the perturbations become nonlinear during the dynamical condensation. Therefore, according to the linear analysis, the cooling layer is expected to split into fragments with various scales.

A new self-similar solution describing the dynamical condensation of a radiative gas is investigated under a plane-parallel geometry. The dynamical condensation is caused by thermal instability. The solution is applicable to generic flow with a net cooling rate per unit volume and time $\propto \rho^2 T^\alpha$, where $\rho$, $T$ and $\alpha$ are density, temperature and a free parameter, respectively. Given $\alpha$, a family of self-similar solutions with one parameter $\eta$ is found in which the central density and pressure evolve as follows: $\rho(x=0,t)\propto (t_\mathrm{c}-t)^{-\eta/(2-\alpha)}$ and $P(x=0,t)\propto (t_\mathrm{c}-t)^{(1-\eta)/(1-\alpha)}$, where $t_\mathrm{c}$ is an epoch when the central density becomes infinite. For $\eta\sim 0$, the solution describes the isochoric mode, whereas for $\eta\sim1$, the solution describes the isobaric mode. The self-similar solutions exist in the range between the two limits; that is, for $0<\eta<1$. No self-similar solution is found for $\alpha>1$. We compare the obtained self-similar solutions with the results of one-dimensional hydrodynamical simulations. In a converging flow, the results of the numerical simulations agree well with the self-similar solutions in the high-density limit. Our self-similar solutions are applicable to the formation of interstellar clouds (HI cloud and molecular cloud) by thermal instability.

In this paper we investigate gravitational instability of shocked gas layers using linear analysis. An unperturbed state is a self-gravitating isothermal layer which grows with time by the accretion of gas through shock fronts due to a cloud-cloud collision. Since the unperturbed state is not static, and cannot be described by a self-similar solution, we numerically solved the perturbation equations and directly integrated them over time. We took account of the distribution of physical quantities across the thickness. Linearized Rankine-Hugoniot relations were imposed at shock fronts as boundary conditions. The following results are found from our unsteady linear analysis: the perturbation initially evolves in oscillatory mode, and begins to grow at a certain epoch. The wavenumber of the fastest growing mode is given by $k=2\sqrt{2\pi G\rho_\mathrm{E} {\cal M\mit}}/c_\mathrm{s}$, where $\rho_\mathrm{E}, c_\mathrm{s}$ and $\cal M\mit$ are the density of parent clouds, the sound velocity and the Mach number of the collision velocity, respectively. For this mode, the transition epoch from oscillatory to growing mode is given by $t_g = 1.2/\sqrt{2\pi G\rho_\mathrm{E} {\cal M\mit}}$. The epoch at which the fastest growing mode becomes non-linear is given by $2.4\delta_0^{-0.1}/\sqrt{2\pi G \rho_\mathrm{E}{\cal M\mit}}$, where $\delta_0$ is the initial amplitude of the perturbation of the column density. As an application of our linear analysis, we investigate criteria for collision-induced fragmentation. Collision-induced fragmentation will occur only when parent clouds are cold, or $\alpha_0=5c_\mathrm{s}^2 R/2G M < 1$, where $R$ and $M$ are the radius and the mass of parent clouds, respectively.