26Al is a short-lived radioactive isotope thought to be injected into the interstellar medium (ISM) by massive stellar winds and supernovae. However, all-sky maps of 26Al emission show a distribution with a much larger scale height and faster rotation speed than either massive stars or the cold ISM. We investigate the origin of this discrepancy using an N-body+hydrodynamics simulation of a Milky-Way-like galaxy, self-consistently including self-gravity, star formation, stellar feedback, and 26Al production. We find no evidence that the Milky Way's spiral structure explains the 26Al anomaly. Stars and the 26Al bubbles they produce form along spiral arms, but, because our simulation produces material arms that arise spontaneously rather than propagating arms forced by an external potential, star formation occurs at arm centres rather than leading edges. As a result, we find a scale height and rotation speed for 26Al similar to that of the cold ISM. However, we also show that a synthetic 26Al emission map produced for a possible Solar position at the edge of a large 26Al bubble recovers many of the major qualitative features of the observed 26Al sky. This suggests that the observed anomalous 26Al distribution is the product of foreground emission from the 26Al produced by a nearby, recent supernova.

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.

The formation and early evolution of low mass young stellar objects (YSOs) are investigated using three-dimensional non-ideal magneto-hydrodynamics simulations. We investigate the evolution of YSOs up to ~ 10^4 yr after protostar formation, at which protostellar mass reaches ~ 0.1 M_\odot . 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 ~ 10^4 yr after protostar formation. Disk mass is comparable to central protostellar mass and gravitational instability develops. In the simulations with small dust size, the warp of the pseudodisk develops ~ 10^4 yr after protostar formation. The warp strengthens magnetic braking in the disk and decreases disk size. Ion-neutral drift can occur in the infalling envelope under the conditions that the typical dust size is a \gtrsim 0.2{\mu}m and the protostar (plus disk) mass is M \gtrsim 0.1 M_\odot. The outflow activity is anti-correlated to the dust size and the strong outflow appears with small dust grains.

The first hydrostatic core, the first quasi-hydrostatic object formed during the star formation process, is still the observational missing link between the prestellar and protostellar phases, mainly due to its short lifetime. Although we have not established a clear method to identify this rare object, recent theoretical studies predict that the first core has millimeter continuum emission and low-velocity outflow with a wide opening angle. An extensive continuum/outflow survey toward a large number of $"$starless$"$ cores in nearby star-forming regions works as a pathfinder. We observed 32 prestellar cores in Taurus with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in 1.3 mm continuum and molecular lines using the ALMA-ACA stand-alone mode. Among the targets, MC35-mm centered at one of the densest $"$starless$"$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitter in our sample. All the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

We have performed survey-type observations in 1 mm continuum and molecular lines toward dense cores (32 prestellar + 7 protostellar) with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in the Taurus molecular clouds using the ALMA-ACA (Atacama Compact Array) stand-alone mode with an angular resolution of 6.$''$5 ($\sim$900 au). The primary purpose of this study is to investigate the innermost part of dense cores toward understanding the initial condition of star formation. In the protostellar cores, contributions from protostellar disks dominate the observed continuum flux with a range of 35-90\% except for the very low-luminosity object. For the prestellar cores, we have successfully confirmed continuum emission from dense gas with a density of $\gtrsim$3 $\times$10$^5$ cm$^{-3}$ toward approximately one-third of the targets. Thanks to the lower spatial frequency coverage with the ACA-7 m array, the detection rate is significantly higher than that of the previous surveys, which have 0 or 1 continuum detected sources among large number of starless samples using the ALMA Main array. The statistical counting method tells us that the lifetime of the prestellar cores until protostar formation therein approaches the free-fall time as the density increases. Among the prestellar cores, at least two targets have possible internal substructures, which are detected in continuum emission with the size scale of $\sim$1000 au if we consider the molecular line (C$^{18}$O and N$_2$D$^{+}$) distributions. These results suggest that small scale fragmentation/coalescence processes within a $\lesssim$0.1 pc scale region occur to determine the final core mass leading to protostars before the dynamical collapse whose central density range is $\sim$(0.3-1) $\times$10$^6$ cm$^{-3}$.

In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems.

Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation.

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

We investigate the roles of magnetically driven disk wind (MDW) and thermally driven photoevaporative wind (PEW) in the long-time evolution of protoplanetary disks. We start simulations from the early phase in which the disk mass is $0.118\,{\mathrm{M}_{\odot}}$ around a $1\,{\mathrm{M}_{\odot}}$ star and track the evolution until the disk is completely dispersed. We incorporate the mass loss by PEW and the mass loss and magnetic braking (wind torque) by MDW, in addition to the viscous accretion, viscous heating, and stellar irradiation. We find that MDW and PEW respectively have different roles: magnetically driven wind ejects materials from an inner disk in the early phase, whereas photoevaporation has a dominant role in the late phase in the outer ($\gtrsim1\,$au) disk. The disk lifetime, which depends on the combination of MDW, PEW, and viscous accretion, shows a large variation of $\sim1$-$20\,$Myr; the gas is dispersed mainly by the MDW and the PEW in the cases with a low viscosity and the lifetime is sensitive to the mass-loss rate and torque of the MDW, whereas the lifetime is insensitive to these parameters when the viscosity is high. Even in disks with very weak turbulence, the cooperation of MDW and PEW enables the disk dispersal within a few Myr.

The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disruptions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no flat remnants with the ratio of minor to major axis lengths c/a < 0.5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain significantly flat asteroids. We suggest that these flat asteroids may be interlopers or formed through low-velocity collisions between member asteroids.

RX J0046.5$-$7308 is a shell-type supernova remnant (SNR) in the Small Magellanic Cloud (SMC). We carried out new $^{12}$CO($J$ = 1-0, 3-2) observations toward the SNR using Mopra and the Atacama Submillimeter Telescope Experiment. We found eight molecular clouds (A-H) along the X-ray shell of the SNR. The typical cloud size and mass are $\sim$10-15 pc and $\sim$1000-3000 $M_{\odot}$, respectively. The X-ray shell is slightly deformed and has the brightest peak in the southwestern shell where two molecular clouds A and B are located. The four molecular clouds A, B, F, and G have high intensity ratios of $^{12}$CO($J$ = 3-2) / $^{12}$CO($J$ = 1-0) $> 1.2$, which are not attributable to any identified internal infrared sources or high-mass stars. The HI cavity and its expanding motion are found toward the SNR, which are likely created by strong stellar winds from a massive progenitor. We suggest that the molecular clouds A-D, F, and G and HI clouds within the wind-blown cavity at $V_\mathrm{LSR} = 117.1$-122.5 km s$^{-1}$ are to be associated with the SNR. The X-ray spectroscopy reveals the dynamical age of $26000^{+1000}_{-2000}$ yr and the progenitor mass of $\gtrsim 30$ $M_{\odot}$, which is also consistent with the proposed scenario. We determine physical conditions of the giant molecular cloud LIRS 36A using the large velocity gradient analysis with archival datasets of the Atacama Large Millimeter/submillimeter Array; the kinematic temperature is $72^{+50}_{-37}$ K and the number density of molecular hydrogen is $1500^{+600}_{-300}$ cm$^{-3}$. The next generation of $\gamma$-ray observations will allow us to study the pion-decay $\gamma$-rays from the molecular clouds in the SMC SNR.

Various instabilities have been proposed as a promising mechanism to accumulate dust. Moreover, some of them are expected to lead to the multiple-ring structure formation and the planetesimal formation in protoplanetary disks. In a turbulent gaseous disk, the growth of the instabilities and the dust accumulation are quenched by turbulent diffusion of dust grains. The diffusion process has been often modeled by a diffusion term in the continuity equation for the dust density. The dust diffusion model, however, does not guarantee the angular momentum conservation in a disk. In this study, we first formulate equations that describe the dust diffusion and also conserve the total angular momentum of a disk. Second, we perform the linear perturbation analysis on the secular gravitational instability (GI) using the equations. The results show that the secular GI is a monotonically growing mode, contrary to the result of previous analyses that found it overstable. We find that the overstability is caused by the non-conservation of the angular momentum. Third, we find a new axisymmetric instability due to the combination of the dust-gas friction and the turbulent gas viscosity, which we refer to as two-component viscous gravitational instability (TVGI). The most unstable wavelength of TVGI is comparable to or smaller than the gas scale height. TVGI accumulates dust grains efficiently, which indicates that TVGI is a promising mechanism for the formation of multiple-ring-like structures and planetesimals. Finally, we examine the validity of the ring formation via the secular GI and TVGI in the HL Tau disk and find both instabilities can create multiple rings whose width is about 10 au at orbital radii larger than 50 au.

The angular momentum of a molecular cloud core plays a key role in star formation, since it is directly related to the outflow and the jet emanating from the new-born star and it eventually results in the formation of the protoplanetary disk. However, the origin of the core rotation and its time evolution are not well understood. Recent observations reveal that molecular clouds exhibit a ubiquity of filamentary structures and that star forming cores are associated with the densest filaments. Since these results suggest that dense cores form primarily in filaments, the mechanism of core formation from filament fragmentation should explain the distribution of the angular momentum of these cores. In this paper we analyze the relation between velocity fluctuations along the filament close to equilibrium and the angular momentum of the cores formed along its crest. We first find that an isotropic velocity fluctuation that follows the three-dimensional Kolmogorov spectrum does not reproduce the observed angular momentum of molecular cloud cores. We then identify the need for a large power at small scales and study the effect of three power spectrum models. We show that the one-dimensional Kolmogorov power spectrum with a slope -5/3 and an anisotropic model with reasonable parameters are compatible with the observations. Our results stress the importance of more detailed and systematic observations of both the velocity structure along filaments and the angular momentum distribution of molecular cloud cores to determine the validity of the mechanism of core formation from filamentary molecular clouds.

The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

The magnetohydrodynamics (MHD) of protoplanetary disks are strongly subject to the non-ideal MHD effects arising from the low ionization fraction of the disk gas. A strong electric field induced by gas motions can heat ionized gas particles and can thereby affect the ionization balance in the disks. Our previous studies revealed that in dusty protoplanetary disks, the Ohmic conductivity decreases with increasing electric field strength until the electrical breakdown of the disk gas occurs. In this study, we extend our previous work to more general cases where both electric and magnetic fields affect the motion of plasma particles, allowing us to study the impacts of plasma heating on all non-ideal MHD effects: Ohmic, Hall, and ambipolar diffusion. We find that the upper limit on the electric current we previously derived applies even in the presence of magnetic fields. Although the Hall and ambipolar resistivities can either increase or decrease with electric field strength depending on the abundance of charged dust grains, the Ohmic resistivity always increases with electric field strength. An order-of-magnitude estimate suggests that a large-scale electric current generated by gas motions in the inner part of protoplanetary disks could exceed the upper limit. This implies that MHD motions of the inner disk, such as the motion driven by the Hall-shear instability, could either get suppressed or trigger electrical breakdown (lightning discharge). This may have important implications for gas accretion and chondrule formation in the inner part of protoplanetary disks.

We present a detailed morphological study of TeV gamma rays, synchrotron radiation, and interstellar gas in the young Type Ia supernova remnant (SNR) RCW 86. We find that the interstellar atomic gas shows good spatial correlation with the gamma rays, indicating that the TeV gamma rays from RCW 86 are likely to be dominantly of hadronic origin. In contrast, the spatial correlation between the interstellar molecular cloud and the TeV gamma rays is poor in the southeastern shell of the SNR. We argue that this poor correlation can be attributed to the low-energy cosmic rays (~ 1 TeV) not penetrating into the dense molecular cloud due to an enhancement of the turbulent magnetic field around the dense cloud of ~ 10-100 $\mu$G. We also find that the southwestern shell, which is bright in both synchrotron X-ray and radio continuum radiation, shows a significant gamma-ray excess compared with the interstellar proton column density, suggesting that leptonic gamma rays via inverse Compton scattering possibly contributes along with hadronic gamma rays. The total cosmic-ray energies of the young TeV gamma-ray SNRs$-$RX J1713.7$-$3946, Vela Jr, HESS J1731$-$347, and RCW 86$-$are roughly similar, which indicates that cosmic rays can be accelerated in both the core-collapse and Type Ia supernovae. The total energy of cosmic rays derived using the gas density, ~ 10$^{48}$-10$^{49}$ erg, gives a safe lower limit due mainly to the low filling factor of interstellar gas within the shell.

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.

Light curve observations of a recently discovered interstellar object 1I/`Oumuamua suggest that this object has an extremely elongated shape with the axis ratio 0.3 or smaller. Planetesimal collisions can produce irregular shapes including elongated shapes. In this paper, we suggest that the extremely elongated shape of 1I/`Oumuamua may be the result of such an impact. To find detailed impact conditions to form the extremely elongated objects, we conduct numerical simulations of planetesimal collisions using Smoothed Particle Hydrodynamics method for elastic dynamics with self-gravity and interparticle friction. Impacts into strengthless target planetesimals with radius 50 m are conducted with various ratios of impactor mass to target mass q, friction angles phi_d, impact velocities v_imp, and impact angles theta_imp. We find that impacts with q \geq 0.5, phi_d \geq 40 degrees, v_imp \leq 40 degrees, and theta_imp \leq 30 degrees produce remnants with the ratio of intermediate to major axis length less than 0.3. This impact condition suggests that the parent protoplanetary disk in the planetesimal collision stage was weakly turbulent (alpha < 10^{-4} for the inner disk) and composed of planetesimals smaller than ~ 7 km to ensure small impact velocity.

Recent observations of the nearby Galactic molecular clouds indicate that the dense gas in molecular clouds have quasi-universal properties on star formation, and observational studies of extra galaxies have shown a galactic-scale correlation between the star formation rate (SFR) and surface density of molecular gas. To reach a comprehensive understanding of both properties, it is important to quantify the fractional mass of the dense gas in molecular clouds f_DG. In particular, for the Milky Way (MW), there are no previous studies resolving the f_DG disk over a scale of several kpc. In this study, the f_DG was measured over 5kpc in the first quadrant of the MW, based on the CO J=1-0 data in l=10-50 deg obtained as part of the FOREST Unbiased Galactic Plane Imaging Survey with the Nobeyama 45-m Telescope (FUGIN) project. The total molecular mass was measured using 12CO, and the dense gas mass was estimated using C18O. The fractional masses including f_DG in the region within ~30% of the distances to the tangential points of the Galactic rotation (e.g., the Galactic Bar, Far-3kpc Arm, Norma Arm, Scutum Arm, Sagittarius Arm, and inter-arm regions) were measured. As a result, an averaged f_DG of 2.9^{+2.6}_{-2.6} % was obtained for the entirety of the target region. This low value suggests that dense gas formation is the primary factor of inefficient star formation in galaxies. It was also found that the f_DG shows large variations depending on the structures in the MW disk. The f_DG in the Galactic arms were estimated to be ~4-5%, while those in the bar and inter-arm regions were as small as ~0.1-0.4%. These results indicate that the formation/destruction processes of the dense gas and their timescales are different for different regions in the MW, leading to the differences in SFRs.

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 carried out new $^{12}$CO($J$ = 1-0, 3-2) observations of a N63A supernova remnant (SNR) from the LMC using ALMA and ASTE. We find three giant molecular clouds toward the northeast, east, and near the center of the SNR. Using the ALMA data, we spatially resolved clumpy molecular clouds embedded within the optical nebulae in both the shock-ionized and photoionized lobes discovered by previous H$\alpha$ and [S II] observations. The total mass of the molecular clouds is $\sim$$800$ $M_{\odot}$ for the shock-ionized region and $\sim$$1700$ $M_{\odot}$ for the photoionized region. Spatially resolved X-ray spectroscopy reveals that the absorbing column densities toward the molecular clouds are $\sim$$1.5$-$6.0\times10^{21}$ cm$^{-2}$, which are $\sim$$1.5$-$15$ times less than the averaged interstellar proton column densities for each region. This means that the X-rays are produced not only behind the molecular clouds, but also in front of them. We conclude that the dense molecular clouds have been completely engulfed by the shock waves, but have still survived erosion owing to their high-density and short interacting time. The X-ray spectrum toward the gas clumps is well explained by an absorbed power-law or high-temperature plasma models in addition to the thermal plasma components, implying that the shock-cloud interaction is efficiently working for both the cases through the shock ionization and magnetic field amplification. If the hadronic gamma-ray is dominant in the GeV band, the total energy of cosmic-ray protons is calculated to be $\sim$$0.3$-$1.4\times10^{49}$ erg with the estimated ISM proton density of $\sim$$190\pm90$ cm$^{-3}$, containing both the shock-ionized gas and neutral atomic hydrogen.

We review the role that magnetic field may have on the formation and evolution of molecular clouds. After a brief presentation and main assumptions leading to ideal MHD equations, their most important correction, namely the ion-neutral drift is described. The nature of the multi-phase interstellar medium (ISM) and the thermal processes that allows this gas to become denser are presented. Then we discuss our current knowledge of compressible magnetized turbulence, thought to play a fundamental role in the ISM. We also describe what is known regarding the correlation between the magnetic and the density fields. Then the influence that magnetic field may have on the interstellar filaments and the molecular clouds is discussed, notably the role it may have on the prestellar dense cores as well as regarding the formation of stellar clusters. Finally we briefly review its possible effects on the formation of molecular clouds themselves. We argue that given the magnetic intensities that have been measured, it is likely that magnetic field is i) responsible of reducing the star formation rate in dense molecular cloud gas by a factor of a few, ii) strongly shaping the interstellar gas by generating a lot of filaments and reducing the numbers of clumps, cores and stars, although its exact influence remains to be better understood. % by a factor on the order of at least 2. Moreover at small scales, magnetic braking is likely a dominant process that strongly modifies the outcome of the star formation process. Finally, we stress that by inducing the formation of more massive stars, magnetic field could possibly enhance the impact of stellar feedback.

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).

Uranus has a tilted rotation axis, which is supposed to be caused by a giant impact. In general, an impact event also changes the internal compositional distribution and drives mass ejection from the planet, which may provide the origin of satellites. Previous studies of the impact simulation of Uranus investigated the resultant angular momentum and the ejected mass distribution. However, the effect of changing the initial condition of the thermal and compositional structure is not studied. In this paper, we perform hydrodynamics simulations for the impact events of Uranus-size ice giants composed of a water core surrounded by a hydrogen envelope using two variant methods of the smoothed particle hydrodynamics. We find that the higher entropy target loses its envelope more efficiently than the low entropy target. However, the higher entropy target gains more angular momentum than the lower entropy target since the higher entropy target has more expanded envelope. We discuss the efficiency of angular momentum transport and the amount of the ejected mass and find a simple analytical model to roughly reproduce the outcomes of numerical simulations. We suggest the range of possible initial conditions for the giant impact on proto-Uranus that reproduces the present rotation tilt of Uranus and sufficiently provides the total angular momentum of the satellite system that can be created from the fragments from the giant impact.

We perform radiation hydrodynamical simulations of protostellar collapse in spherical symmetry, with a special focus on very low-mass objects, i.e. brown dwarfs and sub-brown dwarfs. The inclusion of a realistic equation of state that includes the effect of hydrogen dissociation allows for a modeling of the complete process from the beginning of the collapse until the formation of the protostar. We solve the frequency-dependent radiative transfer equation without any diffusion approximation, using realistic dust and gas opacities. Our results show that the properties of the protostar are essentially independent of the initial conditions, which had previously only been confirmed for higher mass ranges. For very low mass initial conditions, however, we find that the first core phase of the collapse shows some significant differences in the time evolution, with the first core lifetime increasing dramatically because of the reduced accretion rate from the surrounding envelope. We consider the observational implications of this. We also investigate the opposite case of a collapse without any first core phase, which may occur for very unstable initial conditions. In the appendix, we describe a severe numerical problem that causes an unphysical expansion after the formation of the protostar, which may affect other attempts at similar calculations of self-gravitational collapse. We explain the origin of the unphysical behavior and present a solution that can be used in similar investigations.

Pairs of azimuthal intensity decrements at near symmetric locations have been seen in a number of protoplanetary disks. They are most commonly interpreted as the two shadows cast by a highly misaligned inner disk. Direct evidence of such an inner disk, however, remain largely illusive, except in rare cases. In 2012, a pair of such shadows were discovered in scattered light observations of the near face-on disk around 2MASS J16042165-2130284, a transitional object with a cavity $\sim$60 AU in radius. The star itself is a `dipper', with quasi-periodic dimming events on its light curve, commonly hypothesized as caused by extinctions by transiting dusty structures in the inner disk. Here, we report the detection of a gas disk inside the cavity using ALMA observations with $\sim0$\farcs2 angular resolution. A twisted butterfly pattern is found in the moment 1 map of CO (3-2) emission line towards the center, which is the key signature of a high misalignment between the inner and outer disks. In addition, the counterparts of the shadows are seen in both dust continuum emission and gas emission maps, consistent with these regions being cooler than their surroundings. Our findings strongly support the hypothesized misaligned-inner-disk origin of the shadows in the J1604-2130 disk. Finally, the inclination of inner disk would be close to -45 $^{\circ}$ in contrast with 45 $^{\circ}$; it is possible that its internal asymmetric structures cause the variations on the light curve of the host star.

More than a half of asteroids in the main belt have irregular shapes with the ratios of the minor to major axis lengths less than 0.6. One of the mechanisms to create such shapes is collisions between asteroids. The relationship between 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 fracture of rock, and a model of friction of completely damaged rock. We systematically vary the impact velocity from 50 m/s to 400 m/s and the impact angle from 5 degrees to 45 degrees. 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) 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.

Recent numerical simulations of magnetized accretion discs show that the radial-azimuthal component of the stress tensor due to the magnetorotational instability (MRI) is well represented by a power-law function of the gas pressure rather than a linear relation which has been used in most of the accretion disc studies. The exponent of this power-law function which depends on the net flux of the imposed magnetic field is reported in the range between zero and unity. However, the physical consequences of this power-law stress-pressure relation within the framework of the standard disc model have not been explored so far. In this study, the structure of an accretion disc with a power-law stress-pressure relation is studied using analytical solutions in the steady-state and time-dependent cases. The derived solutions are applicable to different accreting systems, and as an illustrative example, we explore structure of protoplanetary discs using these solutions. We show that the slopes of the radial surface density and temperature distributions become steeper with decreasing the stress exponent. However, if the disc opacity is dominated by icy grains and value of the stress exponent is less than about $0.5$, the surface density and temperature profiles become so steep that make them unreliable. We also obtain analytical solutions for the protoplanetary discs which are irradiated by the host star. Using these solutions, we find that the effect of the irradiation becomes more significant with decreasing the stress exponent.

We present Nobeyama 45m telescope C18O, 13CO, and 12CO(1-0) mapping observations towards an interstellar filament in the Taurus molecular cloud. We investigate the gas velocity structure along the filament and in its surrounding parent cloud. The filament is detected in the optically thin C18O emission as a single velocity component, ~1pc long, ~0.06pc wide structure. The C18O emission traces dust column densities larger than ~5x10^21 cm-2. The line-of-sight (LOS) velocity fluctuates along the filament crest with an average amplitude of ~0.2 km/s. The 13CO and 12CO integrated intensity maps show spatially extended emission around the elongated filament. We identify three extended structures with LOS velocities redshifted and blueshifted with respect to the average velocity of the filament identified in C18O. Based on combined analyses of velocity integrated channel maps and intensity variations of the optically thick 12CO spectra on and off the filament, we propose a 3-dimensional structure of the cloud surrounding the filament. We further suggest a multi-interaction scenario where sheet-like extended structures interact, in space and time, with the filament and are responsible for its compression and/or disruption, playing an important role in the star formation history of the filament. We also identify, towards the same field, a very faint filament showing a velocity field compatible with the filament formation process proposed by Inoue et al. (2017), where a filament is formed due to convergence of a flow of matter generated by the bending of the ambient magnetic field structure induced by an interstellar shock compression.

We report ALMA Cycle 3 observations in CO isotopes toward a dense core, MC27/L1521F in Taurus, which is considered to be at an early stage of multiple star formation in a turbulent environment. Although most of the high-density parts of this core are considered to be as cold as $\sim$10 K, high-angular resolution ($\sim$20 au) observations in $^{12}$CO ($J$ = 3--2) revealed complex warm ($>$15--60 K) filamentary/clumpy structures with the sizes from a few tens of au to $\sim$1,000 au. The interferometric observations of $^{13}$CO and C$^{18}$O show that the densest part with arc-like morphologies associated with the previously identified protostar and condensations are slightly redshifted from the systemic velocity of the core. We suggest that the warm CO clouds may be consequences of shock heating induced by interactions among the different density/velocity components that originated from the turbulent motions in the core. However, such a small-scale and fast turbulent motion does not correspond to a simple extension of the line-width-size relation (i.e., Larson'{}s law), and thus the actual origin remains to be studied. The high-angular resolution CO observations are expected to be essential in detecting small-scale turbulent motions in dense cores and to investigate protostar formation therein.

N49 (LHA 120-N49) is a bright X-ray supernova remnant (SNR) in the Large Magellanic Cloud. We present new $^{12}$CO($J$ = 1-0, 3-2), HI, and 1.4 GHz radio-continuum observations of the SNR N49 using Mopra, ASTE, ALMA, and ATCA. We have newly identified three HI clouds using ATCA with an angular resolution of ~20": one associated with the SNR and the others located in front of the SNR. Both the CO and HI clouds in the velocity range from 280-291 km s$^{-1}$ are spatially correlated with both the soft X-rays (0.2-1.2 keV) and the hard X-rays (2.0-7.0 keV) of N49 on a ~10 pc scale. CO 3-2/1-0 intensity ratios indicate higher values of the CO cloud toward the SNR shell with an angular resolution of ~45", and thus a strong interaction was suggested. Using the ALMA, we have spatially resolved CO clumps embedded within or along the southeastern rim of N49 with an angular resolution of ~3''. Three of the CO clumps are rim-brightened on a 0.7-2 pc scale in both hard X-rays and the radio continuum$:$ this provides further evidence for dynamical interactions between the CO clumps and the SNR shock wave. The enhancement of the radio synchrotron radiation can be understood in terms of magnetic-field amplification around the CO clumps via a shock-cloud interaction. We also present a possible scenario in which the recombining plasma that dominates the hard X-rays from N49 was formed via thermal conduction between the SNR shock waves and the cold$/$dense molecular clumps.

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.

Recent radio observations towards nearby galaxies started to map the whole disk and to identify giant molecular clouds (GMCs) even in the regions between galactic spiral structures. Observed variations of GMC mass functions in different galactic environment indicates that massive GMCs preferentially reside along galactic spiral structures whereas inter-arm regions have many small GMCs. Based on the phase transition dynamics from magnetized warm neutral medium to molecular clouds, Kobayashi et al. 2017 proposes a semi-analytical evolutionary description for GMC mass functions including cloud-cloud collision (CCC) process. Their results show that CCC is less dominant in shaping the mass function of GMCs compared with the accretion of dense HI gas driven by the propagation of supersonic shock waves. However, their formulation does not take into account the possible enhancement of star formation by CCC. Radio observations within the Milky Way indicate the importance of CCC for the formation of star clusters and massive stars. In this article, we reformulate the time evolution equation largely modified from Kobayashi et al. 2017 so that we additionally compute star formation subsequently taking place in CCC clouds. Our results suggest that, although CCC events between smaller clouds outnumber the ones between massive GMCs, CCC-driven star formation is mostly driven by massive GMCs > 10^5.5 Msun (where Msun is the solar mass). The resultant cumulative CCC-driven star formation may amount to a few 10 percent of the total star formation in the Milky Way and nearby galaxies.

We develop a new numerical scheme for solving the radiative transfer equation in a spherically symmetric system. This scheme does not rely on any kind of diffusion approximation and it is accurate for optically thin, thick, and intermediate systems. In the limit of a homogeneously distributed extinction coefficient, our method is very accurate and exceptionally fast. We combine this fast method with a slower, but more generally applicable method to describe realistic problems. We perform various test calculations including a simplified protostellar collapse simulation. We also discuss possible future improvements.

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.

Gas giant planets may form early-on during the evolution of protostellar discs, while these are relatively massive. We study how Jupiter-mass planet-seeds (termed protoplanets) evolve in massive, but gravitationally stable (Q>1.5), discs using radiative hydrodynamic simulations. We find that the protoplanet initially migrates inwards rapidly, until it opens up a gap in the disc. Thereafter, it either continues to migrate inwards on a much longer timescale or starts migrating outwards. Outward migration occurs when the protoplanet resides within a gap with gravitationally unstable edges, as a high fraction of the accreted gas is high angular momentum gas from outside the protoplanet's orbit. The effect of radiative heating from the protoplanet is critical in determining the direction of the migration and the eccentricity of the protoplanet. Gap opening is facilitated by efficient cooling that may not be captured by the commonly used beta-cooling approximation. The protoplanet initially accretes at a high rate (1e-3Mj/yr), and its accretion luminosity could be a few tenths of the host star's luminosity, making the protoplanet easily observable (albeit only for a short time). Due to the high gas accretion rate, the protoplanet generally grows above the deuterium-burning mass-limit. Protoplanet radiative feedback reduces its mass growth so that its final mass is near the brown dwarf-planet boundary. The fate of a young planet-seed is diverse and could vary from a gas giant planet on a circular orbit at a few AU from the central star to a brown dwarf on an eccentric, wide orbit.

We present ALMA 0.87 mm continuum, HCO+ J=4--3 emission line, and CO J=3--2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two component transitional disk system with a radial dust gap of 0."2 +/- 0."05, an azimuthal gap in the HCO+ J=4--3 moment zero map, as well as two bridge-like features in the gas data. Interestingly these features in the gas disk have no analogue in the dust disk making them of particular interest. We modeled the dust disk using the Monte Carlo radiative transfer code HOCHUNK3D (Whitney et al. 2013) using a two disk components. We find that there is a radial gap that extends from 15-60 au in all grain sizes which differs from previous work.

We perform non-linear simulation of secular gravitational instability (GI) in protoplanetary disks that has been proposed as a mechanism of the planetesimal formation and the multiple ring formation. Since the timescale of the growth of the secular GI is much longer than the Keplerian rotation period, we develop a new numerical scheme for a long term calculation utilizing the concept of symplectic integrator. With our new scheme, we first investigate the non-linear development of the secular GI in a disk without a pressure gradient in the initial state. We find that the surface density of dust increases by more than a factor of one hundred while that of gas does not increase even by a factor of two, which results in the formation of dust-dominated rings. A line mass of the dust ring tends to be very close to the critical line mass of a self-gravitating isothermal filament. Our results indicate that the non-linear growth of the secular GI provides a powerful mechanism to concentrate the dust. We also find that the dust ring formed via the non-linear growth of the secular GI migrates inward with a low velocity, which is driven by the self-gravity of the ring. We give a semi-analytical expression for the inward migration speed of the dusty ring.

We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.

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.

Magnetorotational instability (MRI) has a potential to generate the vigorous turbulence in protoplanetary disks, although its turbulence strength and accretion stress remains debatable because of the uncertainty of MRI with low ionization fraction. We focus on the heating of electrons by strong electric fields which amplifies nonideal magnetohydrodynamic effects. The heated electrons frequently collide with and stick to dust grains, which in turn decreases the ionization fraction and is expected to weaken the turbulent motion driven by MRI. In order to quantitatively investigate the nonlinear evolution of MRI including the electron heating, we perform magnetohydrodynamical simulation with the unstratified shearing box. We introduce a simple analytic resistivity model depending on the current density by mimicking resistivity given by the calculation of ionization. Our simulation confirms that the electron heating suppresses magnetic turbulence when the electron heating occurs with low current density. We find a clear correlation between magnetic stress and its current density, which means that the magnetic stress is proportional to the squared current density. When the turbulent motion is completely suppressed, laminar accretion flow is caused by ordered magnetic field. We give an analytical description of the laminar state by using a solution of linear perturbation equations with resistivity. We also propose a formula that successfully predicts the accretion stress in the presence of the electron heating.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

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 report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within the dust ring. These spirals are trailing and appear to be about 4 times brighter than their surrounding medium. Their kinematics is consistent with Keplerian rotation at an inclination of 23 degree. The apparent two-arm-spiral pattern is best explained by tidal disturbances created by an unseen companion located at 60--80 au, with dust confined in the pressure bumps created outside this companion orbit. An additional companion at r of 30 au, coinciding with the peak CO brightness and a large pitch angle of the spiral, would help to explain the overall emptiness of the cavity. Alternative mechanisms to excite the spirals are discussed. The origin of the large pitch angle detected here remain puzzling.

We investigate the formation of circumstellar disks and outflows subsequent to the collapse of molecular cloud cores with the magnetic field and turbulence. Numerical simulations are performed by using an adaptive mesh refinement to follow the evolution up to $\sim 1000$~yr after the formation of a protostar. In the simulations, circumstellar disks are formed around the protostars; those in magnetized models are considerably smaller than those in nonmagnetized models, but their size increases with time. The models with stronger magnetic field tends to produce smaller disks. During evolution in the magnetized models, the mass ratios of a disk to a protostar is approximately constant at $\sim 1-10$\%. The circumstellar disks are aligned according to their angular momentum, and the outflows accelerate along the magnetic field on the $10-100$~au scale; this produces a disk that is misaligned with the outflow. The outflows are classified into two types: a magneto-centrifugal wind and a spiral flow. In the latter, because of the geometry, the axis of rotation is misaligned with the magnetic field. The magnetic field has an internal structure in the cloud cores, which also causes misalignment between the outflows and the magnetic field on the scale of the cloud core. The distribution of the angular momentum vectors in a core also has a non-monotonic internal structure. This should create a time-dependent accretion of angular momenta onto the circumstellar disk. Therefore, the circumstellar disks are expected to change their orientation as well as their sizes in the long-term evolutions.

Recent radio observations show that the giant molecular cloud (GMC) mass functions noticeably vary across galactic disks. High-resolution magnetohydrodynamics simulations show that multiple episodes of compression are required for creating a molecular cloud in the magnetized interstellar medium. In this article, we formulate the evolution equation for the GMC mass function to reproduce the observed profiles, for which multiple compression are driven by the network of expanding shells due to HII regions and supernova remnants. We introduce the cloud-cloud collision (CCC) terms in the evolution equation in contrast to the previous work (Inutsuka et al. 2015). The computed time evolution suggests that the GMC mass function slope is governed by the ratio of GMC formation timescale to its dispersal timescale, and that the CCC effect is limited only in the massive-end of the mass function. In addition, we identify a gas resurrection channel that allows the gas dispersed by massive stars to regenerate GMC populations or to accrete onto the pre-existing GMCs. Our results show that almost all of the dispersed gas contribute to the mass growth of pre-existing GMCs in arm regions whereas less than 60 per cent in inter-arm regions. Our results also predict that GMC mass functions have a single power-law exponent in the mass range < 10^5.5 Msun (where Msun represents the solar mass), which is well characterized by GMC self-growth and dispersal timescales. Measurement of the GMC mass function slope provides a powerful method to constrain those GMC timescales and the gas resurrecting factor in various environment across galactic disks.

Godunov Smoothed Particle Hydrodynamics (Godunov SPH) method is a computational fluid dynamics method that utilizes a Riemann solver and achieves the second-order accuracy in space. In this paper, we extend the Godunov SPH method to elastic dynamics by incorporating deviatoric stress tensor that represents the stress for shear deformation or anisotropic compression. Analogously to the formulation of the original Godunov SPH method, we formulate the equation of motion, the equation of energy, and the time evolution equation of deviatoric stress tensor so that the resulting discretized system achieves the second-order accuracy in space. The standard SPH method tends to suffer from the tensile instability that results in unphysical clustering of particles especially in tension-dominated region. We find that the tensile instability can be suppressed by selecting appropriate interpolation for density distribution in the equation of motion for the Godunov SPH method even in the case of elastic dynamics. Several test calculations for elastic dynamics are performed, and the accuracy and versatility of the present method are shown.

Angular momentum distribution and its redistribution are of crucial importance in formation and evolution of circumstellar disks. Many molecular line observations toward young stellar objects indicate that radial distributions of the specific angular momentum $j$ are more or less universal. In small scales, typically $R\lesssim$ 100 AU, the specific angular momenta distribute like $j\propto r^{1/2}$, indicating existence of rotationally supported disk. In outer regions, $R\gtrsim$ 5000 AU, $j$ increases as the radius increases and the slope is steeper than unity, which is supposed to reflect the original angular momentum distributions in the maternal molecular clouds. And lastly there is a connecting region, 100 AU $\lesssim R \lesssim$ 5000 AU, in which $j$-distribution looks almost flat. While this is often interpreted as a consequence of conservation of the specific angular momentum, it actually is insufficient and requires a stronger condition that the initial distribution of $j$ must be spatially uniform. However, this requirement is unrealistic and inconsistent with observations. In this work, we propose a simple alternative explanation; the flat $j$ profile is produced by strong prolongation due to the large velocity gradient in the accreting flow no matter what the initial $j$ distribution is. We provide a simple analytic model for gravitational collapse of molecular clouds. This model can be used to estimate ages of protostars based solely on the observed rotation profile. We demonstrate its validity in comparison with hydrodynamic simulations, and apply the model to young stellar objects such as L1527 IRS, TMC-1A and B335.

Recent ALMA observation has revealed multiple ring structures formed in a protoplanetary disk around HL Tau. Prior to the ALMA observation of HL Tau, theoretical analysis of secular gravitational instability (GI) described a possible formation of multiple ring structures with separations of 13 AU around a radius of 100 AU in protoplanetary disks under certain conditions. In this article, we reanalyze the viability of secular GI by adopting the physical values inferred from the observations. We derive the radial distributions of the most unstable wavelength and the growth timescale of secular GI and verify that secular GI can form the ring structures observed in HL Tau. When a turbulent viscosity coefficient $\alpha$ remains small in inner region of the disk, secular GI grows in the whole disk. Thus, the formation of planetary mass objects should occur first in the inner region as a result of gravitational fragmentation after the non-linear growth of secular GI. In this case, resulting objects are expected to create the gaps at $r$ ~ 10 AU and ~ 30 AU. As a result, all ring structures in HL Tau can be created by secular GI. If this scenario is realized in HL Tau, the outer region corresponds to the earlier growth phase of the most unstable mode of secular GI, and the inner region corresponds to the outcome of the non-linear growth of secular GI. Therefore, this interpretation suggests that we are possibly witnessing both the beginning and end of planet formation in HL Tau.

Recent observations of neutral Galactic interstellar medium showed that filamentary structures of HI clouds are aligned with the interstellar magnetic field. Many interesting applications are proposed based on the alignment such as measurement of magnetic field strength through the Chandrasekhar-Fermi method and removal of polarized foreground dust emissions for the detection of inflationary polarized emission in the cosmic microwave background radiation. However, the physical origin of the alignment remains to be explained. To understand the alignment mechanism, we examine formation of HI clouds triggered by shock compression of diffuse warm neutral medium using three-dimensional magnetohydrodynamic simulations with the effects of optically thin cooling and heating. We show that the shock-compressed diffuse interstellar medium of density n~1 cm^-3 evolves into HI clouds with typical density n~50 cm^-3 via thermal instability driven by cooling, which is consistent with previous studies. We apply a machine vision transformation developed by Clark et al.(2014) to the resulting column density structures obtained by the simulations in order to measure angle correlation between filamentary structures of HI clouds and magnetic field. We find that the orientation of HI filaments depends on the environmental turbulent velocity field, particularly on the strength of shear strain in the direction of the magnetic field, which is controlled by the angle between the shock propagation direction and upstream magnetic field. When the strain along the magnetic field is weak, filamentary components of HI clouds basically lie perpendicular to the magnetic field. However, the filaments have come to align with the magnetic field, if we enhance the turbulent strain along the magnetic field or if we set turbulence in the preshock medium.

Context: Debris disks are important observational clues for understanding planetary-system formation process. In particular, faint warm debris disks may be related to late planet formation near 1 AU. A systematic search of faint warm debris disks is necessary to reveal terrestrial planet formation. Aims: Faint warm debris disks show excess emission that peaks at mid-IR wavelengths. Thus we explore debris disks using the AKARI mid-IR all-sky point source catalog (PSC), a product of the second generation unbiased IR all-sky survey. Methods : We investigate IR excess emission for 678 isolated main-sequence stars for which there are 18 micron detections in the AKARI mid-IR all-sky catalog by comparing their fluxes with the predicted fluxes of the photospheres based on optical to near-IR fluxes and model spectra. The near-IR fluxes are first taken from the 2MASS PSC. However, 286 stars with Ks<4.5 in our sample have large flux errors in the 2MASS photometry due to saturation. Thus we have measured accurate J, H, and Ks band fluxes, applying neutral density (ND) filters for Simultaneous InfraRed Imager for Unbiased Survey (SIRIUS) on IRSF, the \phi 1.4 m near-IR telescope in South Africa, and improved the flux accuracy from 14% to 1.8% on average. Results: We identified 53 debris-disk candidates including eight new detections from our sample of 678 main-sequence stars. The detection rate of debris disks for this work is ~8%, which is comparable with those in previous works by Spitzer and Herschel. Conclusion: The importance of this study is the detection of faint warm debris disks around nearby field stars. At least nine objects have a large amount of dust for their ages, which cannot be explained by the conventional steady-state collisional cascade model.

We present the results of ALMA observations of dust continuum emission and molecular rotational lines toward a dense core MC27 (aka L1521F) in Taurus, which is considered to be at a very early stage of star formation. The detailed column density distribution on size scales from a few tens AU to ~10,000 AU scale are revealed by combining the ALMA (12 m array + 7 m array) data with the published/unpublished single-dish data. The high angular resolution observations at 0.87 mm with a synthesized beam size of ~0."74 x 0."32 reveal that a protostellar source, MMS-1, is not spatially resolved and lacks associated gas emission, while a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part of MC27/L1521F (r $\lesssim$ 3000 AU) is N(H2) ~r$^{-0.4}$, clearly flatter than that of the outer part, ~r$^{-1.0}$. The complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.

Fragmentation of protoplanetary disks 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 disk cooling timescale is comparable to its dynamical timescale. 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 disks, we perform two-dimensional numerical simulations that take into account the effect of the irradiation of the central star and radiation cooling of the disk, and precisely investigate the structure of the spiral arms formed in the protoplanetary disks. 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 disks is divided into two stages: formation of the spiral arms in the disks; and fragmentation of the spiral arm. Our work reduces the condition for the fragmentation of the protoplanetary disks to the condition of the formation of the spiral arm that satisfies Q < 0.6.

We investigate stochastic particle acceleration in accretion flows. It is believed that the magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation without back reaction to the field. The results show that the CRs randomly gain or lose their energies through the interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: The diffusion coefficient in direction of the unperturbed flow is about twenty times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic, and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration efficiently works for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction suffer unphysical runaway acceleration by the Lorentz transformation, which needs to be taken with special care.

We investigate the dynamics of dust grains with various sizes in protoplanetary disk winds driven by magnetorotational turbulence, by simulating the time evolution of the dust grain distribution in the vertical direction. Small dust grains, which are well coupled to the gas, are dragged upward with the upflowing gas, while large grains remain near the midplane of a disk. Intermediate--size grains float at several scale heights from the midplane in time-averated force balance between the downward gravity and the upward gas drag. For the minimum mass solar nebula at 1 AU, dust grains with size of 20 -- 40 $\mu m$ float at 5-10 scale heights from the midplane. Considering the dependence on the distance from the central star, smaller-size grains remain only in an outer region of the disk, while larger-size grains are distributed in a broader region. We also discuss the implication of our result to the observation of dusty material around young stellar objects.

Detailed observations of the solar system planets reveal a wide variety of local atmospheric conditions. Astronomical observations have revealed a variety of extrasolar planets none of which resembles any of the solar system planets in full. Instead, the most massive amongst the extrasolar planets, the gas giants, appear very similar to the class of (young) Brown Dwarfs which are amongst the oldest objects in the universe. Despite of this diversity, solar system planets, extrasolar planets and Brown Dwarfs have broadly similar global temperatures between 300K and 2500K. In consequence, clouds of different chemical species form in their atmospheres. While the details of these clouds differ, the fundamental physical processes are the same. Further to this, all these objects were observed to produce radio and X-ray emission. While both kinds of radiation are well studied on Earth and to a lesser extent on the solar system planets, the occurrence of emission that potentially originate from accelerated electrons on Brown Dwarfs, extrasolar planets and protoplanetary disks is not well understood yet. This paper offers an interdisciplinary view on electrification processes and their feedback on their hosting environment in meteorology, volcanology, planetology and research on extrasolar planets and planet formation.

The modification of Smoothed Particle Hydrodynamics (SPH) method with Riemann Solver is called Godunov SPH. We further extend the Godunov SPH to the description of a medium with negative pressure. Under certain circumstances, the SPH method shows an unphysical instability that results in particle clustering. This instability is called the tensile instability. The tensile instability occurs in positive pressure regions in a regular fluid if a very large number of neighbor particles are used with certain shapes of kernel functions, and it is significant in negative pressure regions that emerge in stretched elastic bodies. We must suppress the tensile instability in SPH for calculations of elastic bodies. In this study, we develop a new technique to remove the tensile instability by extending the Godunov SPH method and conducting a linear stability analysis of the equation of motion for the extended method. We find that the tensile instability can be suppressed by choosing an appropriate order of interpolation in the equation of motion of the Godunov SPH method. We also derive an analytic solution for a Riemann solver for a simple equation of state of an elastic body, and construct a Godunov SPH method for the equation of state that allows negative pressure.

Observations have revealed that a significant number of hot Jupiters have anomalously large radii. Layered convection induced by compositional inhomogeneity has been proposed to account for the radius anomaly of hot Jupiters. To reexamine the impact of the compositional inhomogeneity, we perform an evolutionary calculation by determining convection regime at each evolutionary time step according to the criteria from linear analyses. It is shown that the impact is limited in the case of the monotonic gradient of heavy element abundance. The layered convection is absent for the first 1 Gyr from the formation of hot Jupiters and instead overturning convection develops. The super-adiabaticity of the temperature gradient is limited by the neutrally stable state for the Ledoux stability criterion. The effect of the increased mass of heavy elements essentially compensates the effect of the delayed contraction on the planetary radius caused by compositional inhomogeneity. In addition, even in the case where the layered convection is artificially imposed, this mechanism requires extremely thin layers (~ 10^1-10^3 cm) to account for the observed radius anomaly. The long-term stability of such thin layers remains to be studied. Therefore, if the criteria adopted in this paper are adequate, it might be difficult to explain the inflated radii of hot Jupiters by monotonic gradient of heavy element abundance alone.

In this paper, we propose observational methods for detecting lightning in protoplanetary disks. We do so by calculating the critical electric field strength in the lightning matrix gas (LMG), the parts of the disk where the electric field is strong enough to cause lightning. That electric field accelerates multiple positive ion species to characteristic terminal velocities. In this paper, we present three distinct discharge models, with corresponding critical electric fields. We simulate the position-velocity diagrams and the integrated emission maps for the models. We calculate the measure of sensitivity values for detection of the models, and for distinguishing between the models. At the distance of TW-Hya (54pc), LMG that occupies $2\pi$ in azimuth and $25 \mathrm{au}<r<50 \mathrm{au}$ is $1200\sigma$- to $4000\sigma$-detectable. The lower limits of the radii of $5\sigma$-detectable LMG clumps are between 1.6 au and 5.3 au, depending on the models.

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

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 roles of Alfvenic waves in the weakly-ionized atmosphere of hot Jupiters by carrying out non-ideal magnetohydrodynamic (MHD) simulations with Ohmic diffusion in one-dimensional magnetic flux tubes. Turbulence at the surface excites Alfven waves and they propagate upward to drive hot (~ 10^4 K) outflows. The magnetic diffusion plays an important role in the dissipation of the Alfvenic waves in the weakly ionized atmosphere of hot Jupiters. The mass-loss rate of the spontaneously driven planetary wind is considerably reduced, in comparison with that obtained from ideal MHD simulations because the Alfvenic waves are severely damped at low altitudes in the atmosphere, whereas the wave heating is still important in the heating of the upper atmosphere. Dependence on the surface temperature, planetary radius, and velocity dispersion at the surface is also investigated. We find an inversion phenomenon of the transmitted wave energy flux; the energy flux carried by Alfven waves in the upper atmosphere has a nonmonotonic correlation with the input energy flux from the surface in a certain range of the surface temperature because the resistivity is determined by the global physical properties of the atmosphere in a complicated manner. We also point out that the heating and mass loss are expected only in limited zones if the open magnetic field is confined in the limited regions.

We report detections of two candidate distant submillimeter galaxies (SMGs), MM J154506.4$-$344318 and MM J154132.7$-$350320, which are discovered in the AzTEC/ASTE 1.1 mm survey toward the Lupus-I star-forming region. The two objects have 1.1 mm flux densities of 43.9 and 27.1 mJy, and have Herschel/SPIRE counterparts as well. The Submillimeter Array counterpart to the former SMG is identified at 890 $\mu$m and 1.3 mm. Photometric redshift estimates using all available data from the mid-infrared to the radio suggest that the redshifts of the two SMGs are $z_{\rm photo} \simeq$ 4-5 and 3, respectively. Near-infrared objects are found very close to the SMGs and they are consistent with low-$z$ ellipticals, suggesting that the high apparent luminosities can be attributed to gravitational magnification. The cumulative number counts at $S_{\rm 1.1mm} \ge 25$ mJy, combined with other two 1.1-mm brightest sources, are $0.70 ^{+0.56}_{-0.34}$ deg$^{-2}$, which is consistent with a model prediction that accounts for flux magnification due to strong gravitational lensing. Unexpectedly, a $z > 3$ SMG and a Galactic dense starless core (e.g., a first hydrostatic core) could be similar in the mid-infrared to millimeter spectral energy distributions and spatial structures at least at $\gtrsim 1"$. This indicates that it is necessary to distinguish the two possibilities by means of broad band photometry from the optical to centimeter and spectroscopy to determine the redshift, when a compact object is identified toward Galactic star-forming regions.

Supernova remnants (SNRs) are one of the most energetic astrophysical events and are thought to be the dominant source of Galactic cosmic rays (CRs). A recent report on observations from the Fermi satellite has shown a signature of pion decay in the gamma-ray spectra of SNRs. This provides strong evidence that high-energy protons are accelerated in SNRs. The actual gamma-ray emission from pion decay should depend on the diffusion of CRs in the interstellar medium. In order to quantitatively analyse the diffusion of high-energy CRs from acceleration sites, we have performed test particle numerical simulations of CR protons using a three-dimensional magnetohydrodynamics (MHD) simulation of an interstellar medium swept-up by a blast wave. We analyse the diffusion of CRs at a length scale of order a few pc in our simulated SNR, and find the diffusion of CRs is precisely described by a Bohm diffusion, which is required for efficient acceleration at least for particles with energies above 30 TeV for a realistic interstellar medium. Although we find the possibility of a superdiffusive process (travel distance proportional to t^0.75) in our simulations, its effect on CR diffusion at the length scale of the turbulence in the SNR is limited.

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.

Recent three-dimensional magnetohydrodynamical simulations have identified a disk wind by which gas materials are lost from the surface of a protoplanetary disk, which can significantly alter the evolution of the inner disk and the formation of terrestrial planets. A simultaneous description of the realistic evolution of the gaseous and solid components in a disk may provide a clue for solving the problem of the mass concentration of the terrestrial planets in the solar system. We simulate the formation of terrestrial planets from planetary embryos in a disk that evolves via magnetorotational instability and a disk wind. The aim is to examine the effects of a disk wind on the orbital evolution and final configuration of planetary systems. We perform N-body simulations of sixty 0.1 Earth-mass embryos in an evolving disk. The evolution of the gas surface density of the disk is tracked by solving a one-dimensional diffusion equation with a sink term that accounts for the disk wind. We find that even in the case of a weak disk wind, the radial slope of the gas surface density of the inner disk becomes shallower, which slows or halts the type I migration of embryos. If the effect of the disk wind is strong, the disk profile is significantly altered (e.g., positive surface density gradient, inside-out evacuation), leading to outward migration of embryos inside ~ 1 AU. Disk winds play an essential role in terrestrial planet formation inside a few AU by changing the disk profile. In addition, embryos can undergo convergent migration to ~ 1 AU in certainly probable conditions. In such a case, the characteristic features of the solar system's terrestrial planets (e.g., mass concentration around 1 AU, late giant impact) may be reproduced.

We report preliminary results of large-scale distribution toward the Magellanic supernova remnant N132D using Mopra and Chandra archival datasets. We identified a cavity-like CO structure along the X-ray shell toward the southern half of it. The total mass of associating molecular gas is $\sim10^4 M_\odot$, which is smaller than the previous study by an order of magnitude. Further observations using ALMA, ASTE, and Mopra will reveal the detailed spatial structures and its physical conditions.

We investigated the formation of arc-like structures in the infalling envelope around protostars, motivated by the recent Atacama Large Millimeter/Submillimeter Array (ALMA) observations of the high-density molecular cloud core, MC27/L1521F. We performed self-gravitational hydrodynamical numerical simulations with an adaptive mesh refinement code. A filamentary cloud with a 0.1~pc width fragments into cloud cores because of perturbations due to weak turbulence. The cloud core undergoes gravitational collapse to form multiple protostars, and gravitational torque from the orbiting protostars produces arc structures extending up to a 1000~AU scale. As well as on a spatial extent, the velocity ranges of the arc structures, $\sim0.5\,\mathrm{km\,s}^{-1}$, are in agreement with the ALMA observations. We also found that circumstellar disks are often misaligned in triple system. The misalignment is caused by the tidal interaction between the protostars when they undergo close encounters because of a highly eccentric orbit of the tight binary pair.

The ionization state of the gas plays a key role in the MHD of protoplanetary disks. However, the ionization state can depend on the gas dynamics, because electric fields induced by MHD turbulence can heat up plasmas and thereby affect the ionization balance. To study this nonlinear feedback, we construct an ionization model that includes plasma heating by electric fields and impact ionization by heated electrons, as well as charging of dust grains. We show that when plasma sticking onto grains is the dominant recombination process, the electron abundance in the gas decreases with increasing electric field strength. This is a natural consequence of electron-grain collisions whose frequency increases with electron's random velocity. The decreasing electron abundance may lead to a self-regulation of MHD turbulence. In some cases, not only the electron abundance but also the electric current decreases with increasing field strength in a certain field range. The resulting N-shaped current--field relation violates the fundamental assumption of the non-relativistic MHD that the electric field is uniquely determined by the current density. At even higher field strengths, impact ionization causes an abrupt increase of the electric current as expected by previous studies. We find that this discharge current is multi-valued (i.e., the current--field relation is S-shaped) under some circumstances, and that the intermediate branch is unstable. The N/S-shaped current--field relations may yield hysteresis in the evolution of MHD turbulence in some parts of protoplanetary disks.

We have carried out a spectral analysis of the Suzaku X-ray data in the 0.4-12 keV range toward the shell-type very-high-energy {\gamma}-ray supernova remnant RX J1713.7-3946. The aims of this analysis are to estimate detailed X-rays spectral properties at a high angular resolution up to 2 arcmin, and to compare them with the interstellar gas. The X-ray spectrum is non-thermal and used to calculate absorbing column density, photon index, and absorption-corrected X-ray flux. The photon index varies significantly from 2.1 to 2.9. It is shown that the X-ray intensity is well correlated with the photon index, especially in the west region, with a correlation coefficient of 0.81. The X-ray intensity tends to increase with the averaged interstellar gas density while the dispersion is relatively large. The hardest spectra having the photon index less than 2.4 are found outside of the central 10 arcmin of the SNR, from the north to the southeast (~430 arcmin^2) and from the southwest to the northwest (~150 arcmin^2). The former region shows low interstellar gas density, while the latter high interstellar gas density. We present discussion for possible scenarios which explain the distribution of the photon index and its relationship with the interstellar gas.

We investigate structure of self-gravitating disks, their fragmentation and evolution of the fragments (the clumps) using both analytic approach and three-dimensional radiation hydrodynamics simulations starting from molecular cores. The simulations show that non-local radiative transfer determines disk temperature. We find the disk structure is well described by an analytical model of quasi-steady self-gravitating disk with radial radiative transfer. Because the radiative process is not local and radiation from the interstellar medium cannot be ignored, the local radiative cooling would not be balanced with the viscous heating in a massive disk around a low mass star. In our simulations, there are cases in which the disk does not fragment even though it satisfies the fragmentation criterion based on disk cooling time ($Q \sim 1$ and $\Omega t_{\rm cool}\sim 1$). This indicates that at least the criterion is not sufficient condition for fragmentation. We determine the parameter range for the host cloud core in which disk fragmentation occurs. In addition, we show that the temperature evolution of the center of the clump is close to that of typical first cores and the minimum initial mass of clumps to be about a few Jupiter mass.

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections ($\kappa$ And b, two $\sim$60 M$_{\rm J}$ brown dwarf companions in the Pleiades, PZ Tel B, and CD$-$35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of $\sim$30--100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at $\sim$5 M$_{\rm Jup}$, with a single power law distribution. We find that $p(M, a) \propto M^{-0.65\pm0.60} a^{-0.85\pm0.39}$ (1$\sigma$ errors) provides an adequate fit to our data, with 1.0--3.1\% (68\% confidence) of stars hosting 5--70 $M_{\rm Jup}$ companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs.

The instability in protoplanetary disks due to gas-dust friction and self-gravity of gas and dust is investigated by linear analysis. In the case where the dust to gas ratio is enhanced and turbulence is week, the instability grows, even in gravitationally stable disks, on a timescale of order $10^{4- 5}$yr at a radius of order 100AU. If we ignore the dynamical feedback from dust grains in the gas equation of motion, the instability reduces to the so-called "secular gravitational instability", which was investigated previously as an instability of dust in a fixed background gas flow. In this work, we solve the equations of motion for both gas and dust consistently and find that long-wavelength perturbations are stable, in contrast to the secular gravitational instability in the simplified treatment. This may indicate that we should not neglect small terms in equation of motion if the growth rate is small. The instability is expected to form ring structures in protoplanetary disks. The width of the ring formed at a radius of 100 AU is a few tens of AU. Therefore, the instability is a candidate for the formation mechanism of observed ring-like structures in disks. Another aspect of the instability is the accumulation of dust grains, and hence the instability may play an important role in the formation of planetesimals, rocky protoplanets, and cores of gas giants located at radii $\sim$100 AU. If these objects survive the dispersal of the gaseous component of the disk, they may be the origin of debris disks.

We calculate the mass loss driven by MHD waves from hot Jupiters by using MHD simulations in one-dimensional flux tubes. If a gaseous planet has magnetic field, MHD waves are excited by turbulence at the surface, dissipate at the upper atmosphere, and drive gas outflows. Our calculation shows that mass loss rates are comparable to the observed mass loss rates of hot Jupiters, therefore it is suggested that gas flow driven by MHD waves can plays an important role in the mass loss from gaseous planets. The mass loss rate varies dramatically with radius and mass of a planet: a gaseous planet with a small mass but with inflated radius produces very large mass loss rate. We also derive an analytical expression for the dependence of mass loss rate on planet radius and mass that is in good agreement with the numerical calculation. The mass loss rate also depends on the amplitude of velocity dispersion at the surface of a planet. Thus we expect to infer the condition of the surface and the internal structure of a gaseous planet from future observations of mass loss rate from various exoplanets.

Starless dense cores eventually collapse dynamically, which forms protostars inside them, and the physical properties of the cores determine the nature of the forming protostars. We report ALMA observations of dust continuum emission and molecular rotational lines toward MC27 or L1521F, which is considered to be very close to the first protostellar core phase. We found a few starless high-density cores, one of which has a very high density of ~10^7 cm^-3, within a region of several hundred AU around a very low-luminosity protostar detected by Spitzer. A very compact bipolar outflow with a dynamical timescale of a few hundred years was found toward the protostar. The molecular line observation shows several cores with an arc-like structure, possibly due to the dynamical gas interaction. These complex structures revealed in the present observations suggest that the initial condition of star formation is highly dynamical in nature, which is considered to be a key factor in understanding fundamental issues of star formation such as the formation of multiple stars and the origin of the initial mass function of stars.

We investigate the formation of multiple-planet systems in the presence of a hot Jupiter using extended N-body simulations that are performed simultaneously with semi-analytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of a hot Jupiter. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal-mass; the largest planet constitutes more than 50 percent of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper (Ogihara et al. 2013), we have found a new physical mechanism of induced migration of the hot Jupiter, which is called a crowding-out. If the hot Jupiter opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the hot Jupiter remains. We also discuss angular momentum transfer between the planets and disk.

We examine whether the magnetorotational instability (MRI) can serve as a mechanism of angular momentum transport in circumplanetary disks. For the MRI to operate the ionization degree must be sufficiently high and the magnetic pressure must be sufficiently lower than the gas pressure. We calculate the spatial distribution of the ionization degree and search for the MRI-active region where the two criteria are met. We find that there can be thin active layers at the disk surface depending on the model parameters, however, we find hardly any region which can sustain well-developed MRI turbulence; when the magnetic field is enhanced by MRI turbulence at the disk surface layer, a magnetically dominated atmosphere encroaches on a lower altitude and a region of well-developed MRI turbulence becomes smaller. We conclude that if there are no angular momentum transfer mechanisms other than MRI in gravitationally stable circumplanetary disks, gas is likely to pile up until disks become gravitationally unstable, and massive disks may survive for a long time.

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.

(Abridged) We report results of three dimensional MHD simulations of global accretion disks threaded with weak vertical magnetic fields. We perform the simulations in the spherical coordinates with different temperature profiles and accordingly different rotation profiles. In the cases with a spatially constant temperature, because the rotation frequency is vertically constant in the equilibrium condition, general properties of the turbulence are quantitatively similar to those obtained in local shearing box simulations. On the other hand, in the cases with a radially variable temperature profile, the vertical differential rotation, which is inevitable in the equilibrium condition, winds up the magnetic field lines, in addition to the usual radial differential rotation. As a result, the coherent wound magnetic fields contribute to the Maxwell stress in the surface regions. We obtain nondimensional density and velocity fluctuations ~0.1-0.2 at the midplane. The azimuthal power spectra of the magnetic fields show shallow slopes, ~m^0- m^{-1}, than those of velocity and density. We observe intermittent and structured disk winds driven by the Poynting flux associated with the MHD turbulence. The Poynting flux originating from magnetic tension is injected from the regions above a scale height towards both the midplane and the surfaces. Related to this, sound waves are directed to the midplane from the surface regions. The mass accretion mainly occurs near the surfaces and the gas near the midplane slowly moves outward, which causes large-scale meridional circulations. The vertical magnetic fields are also dragged inward in the surface regions, while they stochastically move outward and inward around the midplane. We also discuss an observational implication of induced spiral structure in the simulated turbulent disks to protoplanetary disks.

Recent studies of the nearest star-forming clouds of the Galaxy at submillimeter wavelengths with the Herschel Space Observatory have provided us with unprecedented images of the initial and boundary conditions of the star formation process. The Herschel results emphasize the role of interstellar filaments in the star formation process and connect remarkably well with nearly a decade's worth of numerical simulations and theory that have consistently shown that the ISM should be highly filamentary on all scales and star formation is intimately related to self-gravitating filaments. In this review, we trace how the apparent complexity of cloud structure and star formation is governed by relatively simple universal processes - from filamentary clumps to galactic scales. We emphasize two crucial and complementary aspects: (i) the key observational results obtained with Herschel over the past three years, along with relevant new results obtained from the ground on the kinematics of interstellar structures, and (ii) the key existing theoretical models and the many numerical simulations of interstellar cloud structure and star formation. We then synthesize a comprehensive physical picture that arises from the confrontation of these observations and simulations.

We investigate formation of close-in terrestrial planets from planetary embryos under the influence of a hot Jupiter (HJ) using gravitational N-body simulations that include gravitational interactions between the gas disk and the terrestrial planet (e.g., type I migration). Our simulations show that several terrestrial planets efficiently form outside the orbit of the HJ, making a chain of planets, and all of them gravitationally interact directly or indirectly with the HJ through resonance, which leads to inward migration of the HJ. We call this mechanism of induced migration of the HJ as "crowding out." The HJ is eventually lost by collision with the central star, and only several terrestrial planets remain. We also find that the efficiency of the crowding-out effect depends on model parameters; for example, the heavier the disk is, the more efficient the crowding out is. When planet formation occurs in a massive disk, the HJ can be lost to the central star and is never observed. On the other hand, for a less massive disk, the HJ and terrestrial planets can coexist; however, the companion planets can be below the detection limit of current observations. In both cases, systems with the HJ and terrestrial planets have little chance for detection. Therefore, our model naturally explains the lack of companion planets in HJ systems regardless of the disk mass. In effect, our model provide a theoretical prediction for future observations; additional planets can be discovered just outside the HJ, and their masses should generally be small.

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

RX J1713.7-3946 is the most remarkable very-high-energy \gamma-ray supernova remnant which emits synchrotron X-rays without thermal features. We made a comparative study of CO, HI and X-rays in order to better understand the relationship between the X-rays, and the molecular and atomic gas. The results indicate that the X-rays are enhanced around the CO and HI clumps on a pc scale but are decreased inside the clumps on a 0.1 pc scale. Magnetohydrodynamic numerical simulations of the shock interaction with molecular and atomic gas indicate that the interaction between the shock waves and the clumps excite turbulence which amplifies the magnetic field around the clumps (Inoue et al. 2012). We suggest that the amplified magnetic field around the CO and HI clumps enhances the synchrotron X-rays and possibly the acceleration of cosmic-ray electrons.

In this article, we present a state-of-the-art algorithm for solving the relativistic viscous hydrodynamics equation with the QCD equation of state. The numerical method is based on the second-order Godunov method and has less numerical dissipation, which is crucial in describing of quark-gluon plasma in high-energy heavy-ion collisions. We apply the algorithm to several numerical test problems such as sound wave propagation, shock tube and blast wave problems. In sound wave propagation, the intrinsic numerical viscosity is measured and its explicit expression is shown, which is the second-order of spatial resolution both in the presence and absence of physical viscosity. The expression of the numerical viscosity can be used to determine the maximum cell size in order to accurately measure the effect of physical viscosity in the numerical simulation.

We report the results of our three-dimensional radiation hydrodynamics simulation of collapsing unmagnetized molecular cloud cores. We investigate the formation and evolution of the circumstellar disk and the clumps formed by disk fragmentation. Our simulation shows that disk fragmentation occurs in the early phase of circumstellar disk evolution and many clumps form. The clump can be represented by a polytrope sphere of index $n \sim 3$ and $n \gtrsim 4$ at central temperature $T_c \lesssim100$ K and $T_c \gtrsim 100$ K, respectively. We demonstrate, numerically and theoretically, that the maximum mass of the clump, beyond which it inevitably collapses, is $\sim 0.03 M_{\odot}$. The entropy of the clump increases during its evolution, implying that evolution is chiefly determined by mass accretion from the disk rather than by radiative cooling. Although most of the clumps rapidly migrate inward and finally fall onto the protostar, a few clumps remain in the disk. The central density and temperature of the surviving clump rapidly increase and the clump undergoes a second collapse within 1000 - 2000 years after its formation. In our simulation, three second cores of masses $0.2\msun$, $0.15\msun$, and $0.06\msun$ formed. These are protostars or brown dwarfs rather than protoplanets. For the clumps to survive as planetary-mass objects, the rapid mass accretion should be prevented by some mechanisms.

The formation of a circumstellar disk in collapsing cloud cores is investigated with three-dimensional magnetohydrodynamic simulations. We prepare four types of initial cloud having different density profiles and calculate their evolution with or without a sink. To investigate the effect of magnetic dissipation on disk formation, the Ohmic dissipation is considered in some models. Calculations show that disk formation is very sensitive to both the initial cloud configuration and the sink treatment. The disk size considerably differs in clouds with different density profiles even when the initial clouds have almost the same mass-to-flux ratio. Only a very small disk (\sim 10 AU in size) appears in clouds with a uniform density profile, whereas a large disk (\sim 100 AU in size) forms in clouds with a Bonnor-Ebert density profile. In addition, a large sink accretion radius numerically impedes disk formation during the main accretion phase and tends to foster the misleading notion that disk formation is completely suppressed by magnetic braking. The protostellar outflow is also greatly affected by the sink properties. A sink accretion radius of \lesssim 1 AU and sink threshold density of \gtrsim 10^ 13 cm^-3 are necessary for investigating disk formation during the main accretion phase.

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.

We study how the magnetorotational instability (MRI) in protoplanetary disks is affected by the electric discharge caused by the electric field in the resistive MHD. We have performed three-dimensional shearing box simulations with various values of plasma beta and electrical breakdown models. We find the self-sustainment of the MRI in spite of the high resistivity. The instability gives rise to the large electric field that causes the electrical breakdown, and the breakdown maintains the high ionization degree required for the instability. The condition for this self-sustained MRI is set by the balance between the energy supply from the shearing motion and the energy consumed by the Ohmic dissipation. We apply the condition to various disk models and study where the active, self-sustained, and dead zones of MRI are. In the fiducial minimum-mass solar nebula (MMSN) model, the newly-found sustained zone occupies only the limited volume of the disk. In the late-phase gas-depleted disk models, however, the sustained zone occupies larger volume of the disk.

We investigate the formation process of self-gravitating protoplanetary disks in unmagnetized molecular clouds. The angular momentum is redistributed by the action of gravitational torques in the massive disk during its early formation. We develop a simplified one-dimensional accretion disk model that takes into account the infall of gas from the envelope onto the disk and the transfer of angular momentum in the disk with an effective viscosity. First we evaluate the gas accretion rate from the cloud core onto the disk by approximately estimating the effects of gas pressure and gravity acting on the cloud core. We formulate the effective viscosity as a function of the Toomre Q parameter that measures the local gravitational stability of the rotating thin disk. We use a function for viscosity that changes sensitively with Q when the disk is gravitationally unstable. We find a strong self-regulation mechanism in the disk evolution. During the formation stage of protoplanetary disks, the evolution of the surface density does not depend on the other details of the modeling of effective viscosity, such as the prefactor of the viscosity coefficient. Next, to verify our model, we compare the time evolution of the disk calculated with our formulation with that of three-dimensional hydrodynamical simulations. The structures of the resultant disks from the one-dimensional accretion disk model agree well with those of the three-dimensional simulations. Our model is a useful tool for the further modeling of chemistry, radiative transfer, and planet formation in protoplanetary disks.

Using 3D MHD simulation with the effects of radiative cooling/heating, chemical reactions, and thermal conduction, we investigate the formation of molecular cloud in the ISM. We consider the formation of molecular cloud by accretion of the HI clouds as suggested in recent observations. The simulation shows that the initial HI medium is compressed and piled up behind the shock waves induced by the accretion flows. Since the initial medium is highly inhomogeneous as a consequence of the thermal instability, the formed molecular cloud becomes very turbulent owing to the development of the Richtmyer-Meshkov instability. The structure of the post shock region is composed of dense cold gas (T<100 K) and diffuse warm gas (T>1,000 K), which are spatially well mixed owing to the turbulence. Because the energy source of the turbulence is the accretion flows, the turbulence is highly anisotropic biased toward the direction of accretion flows. The kinetic energy of the turbulence dominates the thermal, magnetic, and gravitational energies in the total 10 Myr evolution. However, the kinetic energy measured by using the CO-fraction-weighted density is comparable to the other energies. This suggests that the true kinetic energy of turbulence in molecular cloud as a hole can be much larger than the kinetic energy of turbulence estimated by line-width of molecular emissions. The clumps in the molecular cloud show statistically homogeneous evolution as follows: The typical plasma beta of the clumps is roughly constant <beta>~ 0.4. The size-velocity dispersion relation show dv ~ 1.5 km s^{-1} (l/1 pc)^{0.5}, irrespective of the density. The clumps evolve toward magnetically supercritical, gravitationally unstable cores. The clumps seem to evolve into cores that satisfy the condition for fragmentation into binary. These statistical properties may provide the initial condition of star formation.

In this letter we report on turbulent acceleration of the dissipation of magnetic field in the postshock re- gion of a Poynting flux-dominated flow, such as the Crab pulsar wind nebula. We have performed two- dimensional resistive relativistic magnetohydrodynamics simulations of subsonic turbulence driven by the Richtmyer-Meshkov instability at the shock fronts of the Poynting flux-dominated flows in pulsar winds. We find that turbulence stretches current sheets which substantially enhances the dissipation of magnetic field, and that most of the initial magnetic field energy is dissipated within a few eddy-turnover times. We also develop a simple analytical model for turbulent dissipation of magnetic field that agrees well with our simulations. The analytical model indicates that the dissipation rate does not depend on resistivity even in the small resistivity limit. Our findings can possibly alleviate the {\sigma}-problem in the Crab pulsar wind nebulae.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

In order to investigate the origin of the interstellar turbulence, detailed observations in the CO J=1--0 and 3--2 lines have been carried out in an interacting region of a molecular cloud with an HII region. As a result, several 1,000 to 10,000 AU scale cloudlets with small velocity dispersion are detected, whose systemic velocities have a relatively large scatter of a few km/s. It is suggested that the cloud is composed of small-scale dense and cold structures and their overlapping effect makes it appear to be a turbulent entity as a whole. This picture strongly supports the two-phase model of turbulent medium driven by thermal instability proposed previously. On the surface of the present cloud, the turbulence is likely to be driven by thermal instability following ionization shock compression and UV irradiation. Those small scale structures with line width of ~ 0.6 km/s have a relatively high CO line ratio of J=3--2 to 1--0, 1 < R(3-2/1-0) < 2. The large velocity gradient analysis implies that the 0.6 km/s width component cloudlets have an average density of 10^{3-4} cm^{-3}, which is relatively high at cloud edges, but their masses are only < 0.05 M_{sun}.

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.

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

We perform a linear stability analysis of a dust layer in a turbulent gas disk. Youdin (2011) investigated the secular gravitational instability of a dust layer using hydrodynamic equations with a turbulent diffusion term. We obtain essentially the same result independently of Youdin (2011). In the present analysis, we restrict the area of interest to small dust particles, while investigating the secular gravitational instability in a more rigorous manner. We discuss the time evolution of the dust surface density distribution using a stochastic model and derive the advection-diffusion equation. The validity of the analysis by Youdin (2011) is confirmed in the strong drag limit. We demonstrate quantitatively that the finite thickness of a dust layer weakens the secular gravitational instability and that the density-dependent diffusion coefficient changes the growth rate. We apply the obtained results to the turbulence driven by the shear instability and find that the secular gravitational instability is faster than the radial drift when the gas density is three times as large as that in the minimum-mass disk model. If the dust particles are larger than chondrules, the secular gravitational instability grows within the lifetime of a protoplanetary disk.

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.

RX J1713.7$-$3946 is the most remarkable TeV $\gamma$-ray SNR which emits $\gamma$-rays in the highest energy range. We made a new combined analysis of CO and \ion{H}{1} in the SNR and derived the total protons in the interstellar medium (ISM). We have found that the inclusion of the \ion{H}{1} gas provides a significantly better spatial match between the TeV $\gamma$-rays and ISM protons than the H$_2$ gas alone. In particular, the southeastern rim of the $\gamma$-ray shell has a counterpart only in the \ion{H}{1}. The finding shows that the ISM proton distribution is consistent with the hadronic scenario that comic ray (CR) protons react with ISM protons to produce the $\gamma$-rays. This provides another step forward for the hadronic origin of the $\gamma$-rays by offering one of the necessary conditions missing in the previous hadronic interpretations. We argue that the highly inhomogeneous distribution of the ISM protons is crucial in the origin of the $\gamma$-rays. Most of the neutral gas was likely swept up by the stellar wind of an OB star prior to the SNe to form a low-density cavity and a swept-up dense wall. The cavity explains the low-density site where the diffusive shock acceleration of charged particles takes place with suppressed thermal X-rays, whereas the CR protons can reach the target protons in the wall to produce the $\gamma$-rays. The present finding allows us to estimate the total CR proton energy to be $\sim 10^{48}$ ergs, 0.1% of the total energy of a SNe.

Using three-dimensional magnetohydrodynamics simulations, we investigate general properties of a blast wave shock interacting with interstellar clouds. The shock-cloud interaction generates a turbulent shell through the vorticity generations. In the turbulent shell, the magnetic field is amplified as a result of turbulent dynamo action. In the case of a young supernova remnant, the corresponding strength of the magnetic field is approximately 1 mG. The propagation speed of the shock wave is significantly stalled in the clouds because of the high density, while the shock maintains a high velocity in the diffuse surrounding. In addition, when the shock wave hits the clouds, reflection shock waves are generated that propagate back into the shocked shell. From these simulation results, many observational characteristics of a young SNR RX J1713.7-3946 that is suggested to be interacting with molecular clouds, can be explained as follows: The reflection shocks can accelerate particles in the turbulent downstream region where the magnetic field strength reaches 1mG, which causes short-time variability of synchrotron X-rays. Since the shock velocity is stalled locally in the clouds, the temperature in the shocked cloud is suppressed below 1 keV. Thus, thermal X-ray line emissions would be faint even if the SNR is interacting with molecular clouds. We also find that the photon index of the $\pi^0$-decay gamma rays generated by cosmic-ray protons can be 1.5, because the penetration depth of high-energy particles into the clumpy clouds depends on their energy, which is consistent with the recent gamma-ray observation by Fermi space telescope.

We develop a fast and accurate calculation method for ionization degrees in protoplanetary and circumplanetary disks including dust grains. We apply our method to calculate the ionization degree of circumplanetary disks. It is important to understand the structure and evolution of protoplanetary/circumplanetary disks since they are thought to be the sites of planet/satellite formation. The turbulence that causes gas accretion is supposed to be driven by magnetorotational instability (MRI) that occurs only when the ionization degree is high enough for magnetic field to be coupled to gas. We calculate the ionization degrees in circumplanetary disks to estimate the sizes of MRI-inactive regions. We properly include the effect of dust grains because they efficiently capture charged particles and make ionization degree lower. Inclusion of dust grains complicates the reaction equations and requires expensive computation. In order to accelerate the calculation of ionization reactions, we develop a semianalytic method based on the charge distribution model proposed previously. This method enables us to study the ionization state of disks for a wide range of model parameters. For a previous model of circum-Jovian disk, we find that an MRI-inactive region covers almost all regions even without dust grains. This suggests that the gas accretion rates in circumplanetary disks are much smaller than previously thought.

In this paper, we develop a stable and fast numerical scheme for relativistic dissipative hydrodynamics based on Israel-Stewart theory. Israel-Stewart theory is a stable and causal description of dissipation in relativistic hydrodynamics although it includes relaxation process with the timescale for collision of constituent particles, which introduces stiff equations and makes practical numerical calculation difficult. In our new scheme, we use Strang's splitting method, and use the piecewise exact solutions for solving the extremely short timescale problem. In addition, since we split the calculations into inviscid step and dissipative step, Riemann solver can be used for obtaining numerical flux for the inviscid step. The use of Riemann solver enables us to capture shocks very accurately. Simple numerical examples are shown. The present scheme can be applied to various high energy phenomena of astrophysics and nuclear physics.

Using resistive magnetohydrodynamics simulation, we investigate circumstellar disk formation in a strongly magnetized cloud. As the initial state, an isolated cloud core embedded in a low-density interstellar medium with a uniform magnetic field is adopted. The cloud evolution is calculated until almost all gas inside the initial cloud falls onto either the circumstellar disk or a protostar, and a part of the gas is ejected into the interstellar medium by the protostellar outflow driven by the circumstellar disk. In the early main accretion phase, the disk size is limited to \sim 10 AU because the angular momentum of the circumstellar disk is effectively transferred by both magnetic braking and the protostellar outflow. In the later main accretion phase, however, the circumstellar disk grows rapidly and exceeds 100 AU by the end of the main accretion phase. This rapid growth of the circumstellar disk is caused by the depletion of the infalling envelope, while magnetic braking is effective when the infalling envelope is more massive than the circumstellar disk. The infalling envelope cannot brake the circumstellar disk when the latter is more massive than the former. In addition, the protostellar outflow weakens and disappears in the later main accretion phase, because the outflow is powered by gas accretion onto the circumstellar disk. Although the circumstellar disk formed in a magnetized cloud is considerably smaller than that in an unmagnetized cloud, a circumstellar disk exceeding 100 AU can form even in a strongly magnetized cloud.

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.

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 report high-resolution 1.6 $\micron$ polarized intensity ($PI$) images of the circumstellar disk around the Herbig Ae star AB Aur at a radial distance of 22 AU ($0."15$) up to 554 AU (3.$"$85), which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarimetry. We revealed complicated and asymmetrical structures in the inner part ($\lesssim$140 AU) of the disk, while confirming the previously reported outer ($r$ $\gtrsim$200 AU) spiral structure. We have imaged a double ring structure at $\sim$40 and $\sim$100 AU and a ring-like gap between the two. We found a significant discrepancy of inclination angles between two rings, which may indicate that the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is $\sim$45 AU or less) within two rings as well as three prominent $PI$ peaks at $\sim$40 AU. The observed structures, including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, provide essential information for understanding the formation mechanism of recently detected wide-orbit ($r$ $>$20 AU) planets.

The formation and evolution of a circumstellar disk in magnetized cloud cores is investigated from prestellar core stage until sim 10^4 yr after protostar formation. In the circumstellar disk, fragmentation first occurs due to gravitational instability in a magnetically inactive region, and substellar-mass objects appear. The substellar-mass objects lose their orbital angular momenta by gravitational interaction with the massive circumstellar disk and finally fall onto the protostar. After this fall, the circumstellar disk increases its mass by mass accretion and again induces fragmentation. The formation and falling of substellar-mass objects are repeated in the circumstellar disk until the end of the main accretion phase. In this process, the mass of fragments remain small, because the circumstellar disk loses its mass by fragmentation and subsequent falling of fragments before it becomes very massive. In addition, when fragments orbit near the protostar, they disturb the inner disk region and promote mass accretion onto the protostar. The orbital motion of substellar-mass objects clearly synchronizes with the time variation of the accretion luminosity of the protostar. Moreover, as the objects fall, the protostar shows a strong brightening for a short duration. The intermittent protostellar outflows are also driven by the circumstellar disk whose magnetic field lines are highly tangled owing to the orbital motion of fragments. The time-variable protostellar luminosity and intermittent outflows may be a clue for detecting planetary-mass objects in the circumstellar disk.

Recent gamma-ray observations of middle-aged supernova remnants revealed a mysterious broken power-law spectrum. Using three-dimensional magnetohydrodynamics simulations, we show that the interaction between a supernova blast wave and interstellar clouds formed by thermal instability generates multiple reflected shocks. The typical Mach numbers of the reflected shocks are shown to be M ~ 2 depending on the density contrast between the diffuse intercloud gas and clouds. These secondary shocks can further energize cosmic-ray particles originally accelerated at the blast-wave shock. This "two-step" acceleration scenario reproduces the observed gamma-ray spectrum and predicts the high-energy spectral index ranging approximately from 3 to 4.

We analyze the physical processes of gap formation in an inviscid protoplanetary disk with an embedded protoplanet using two-dimensional local shearing-sheet model. Spiral density wave launched by the planet shocks and the angular momentum carried by the wave is transferred to the background flow. The exchange of the angular momentum can affect the mass flux in the vicinity of the planet to form an underdense region, or gap, around the planetary orbit. We first perform weakly non-linear analyses to show that the specific vorticity formed by shock dissipation of density wave can be a source of mass flux in the vicinity of the planet, and that the gap can be opened even for low-mass planets unless the migration of the planet is substantial. We then perform high resolution numerical simulations to check analytic consideration. By comparing the gap opening timescale and type I migration timescale, we propose a criterion for the formation of underdense region around the planetary orbit that is qualitatively different from previous studies. The minimum mass required for the planet to form a dip is twice as small as previous studies if we incorporate the standard values of type I migration timescale, but it can be much smaller if there is a location in the disk where type I migration is halted.

RX J1713.7-3946 is one of the TeV {\gamma}-ray supernova remnants (SNRs) emitting synchrotron X rays. The SNR is associated with molecular gas located at ~1 kpc. We made new molecular observations toward the dense cloud cores, peaks A, C and D, in the SNR in the 12CO(J=2-1) and 13CO(J=2-1) transitions at angular resolution of 90". The most intense core in 13CO, peak C, was also mapped in the 12CO(J=4-3) transition at angular resolution of 38". Peak C shows strong signs of active star formation including bipolar outflow and a far-infrared protostellar source and has a steep gradient with a r^{-2.2$\pm$0.4} variation in the average density within radius r. Peak C and the other dense cloud cores are rim-brightened in synchrotron X rays, suggesting that the dense cloud cores are embedded within or on the outer boundary of the SNR shell. This confirms the earlier suggestion that the X rays are physically associated with the molecular gas (Fukui et al. 2003). We present a scenario where the densest molecular core, peak C, survived against the blast wave and is now embedded within the SNR. Numerical simulations of the shock-cloud interaction indicate that a dense clump can indeed survive shock erosion, since shock propagation speed is stalled in the dense clump. Additionally, the shock-cloud interaction induces turbulence and magnetic field amplification around the dense clump that may facilitate particle acceleration in the lower-density inter-clump space leading to the enhanced synchrotron X rays around dense cores.

We use resistive magnetohydrodynamical simulations with the nested grid technique to study the formation of protoplanetary disks around protostars from molecular cloud cores that provide the realistic environments for planet formation. We find that gaseous planetary-mass objects are formed in the early evolutionary phase by gravitational instability in regions that are decoupled from the magnetic field and surrounded by the injection points of the magnetohydrodynamical outflows during the formation phase of protoplanetary disks. Magnetic decoupling enables massive disks to form and these are subject to gravitational instability, even at ~ 10 AU. The frequent formation of planetary-mass objects in the disk suggests the possibility of constructing a hybrid planet formation scenario, where the rocky planets form later under the influence of the giant planets in the protoplanetary disk.

By constructing a global model based on 3D local magnetohydrodynamical (MHD) simulations, we show that the disk wind driven by magnetorotational instability (MRI) plays a significant role in the dispersal of the gas component of proto-planetary disks. Because the mass loss time scale by the MRI-driven disk winds is proportional to the local Keplerian rotation period, a gas disk dynamically evaporates from the inner region with possibly creating a gradually expanding inner hole, while a sizable amount of the gas remains in the outer region. The disk wind is highly time-dependent with quasi-periodicity of several times Keplerian rotation period at each radius, which will be observed as time-variability of protostar-protoplanetary disk systems. These features persistently hold even if a dead zone exists because the disk winds are driven from the surface regions where ionizing cosmic rays and high energy photons can penetrate. Moreover, the predicted inside-out clearing significantly suppresses the infall of boulders to a central star and the Type I migration of proto-planets which are favorable for the formation and survival of planets.

We investigate the formation process of planetesimals from the dust layer by the gravitational instability in the gas disk using local $N$-body simulations. The gas is modeled as a background laminar flow. We study the formation process of planetesimals and its dependence on the strength of the gas drag. Our simulation results show that the formation process is divided into three stages qualitatively: the formation of wake-like density structures, the creation of planetesimal seeds, and their collisional growth. The linear analysis of the dissipative gravitational instability shows that the dust layer is secularly unstable although Toomre's $Q$ value is larger than unity. However, in the initial stage, the growth time of the gravitational instability is longer than that of the dust sedimentation and the decrease in the velocity dispersion. Thus, the velocity dispersion decreases and the disk shrinks vertically. As the velocity dispersion becomes sufficiently small, the gravitational instability finally becomes dominant. Then wake-like density structures are formed by the gravitational instability. These structures fragment into planetesimal seeds. The seeds grow rapidly owing to mutual collisions.

In this paper, we study the Cauchy problem of the linearized kinetic equations for the models of Marle and Anderson-Witting, and compare these dispersion relations with the 14-moment theory. First, we propose a modification of the Marle model to improve the resultant transport coefficients in accord with those obtained by the full Boltzmann equation. Using the modified Marle model and Anderson-Witting model, we calculate dispersion relations that are kinetically correct within the validity of the BGK approximation. The 14-moment theory that includes the time derivative of dissipation currents has causal structure, in contrast to the acausal first-order Chapman-Enskog approximation. However, the dispersion relation of the 14-moment theory does not accurately describe the result of the kinetic equation. Thus, our calculation indicates that keeping these second-order terms does not simply correspond to improving the physical description of the relativistic hydrodynamics.

We study the linearized kinetic equation of relaxation model which was proposed by Bhatnagar, Gross and Krook (also called BGK model) and solve the dispersion relation. Using the solution of the dispersion relation, we analyze the relaxation of the macroscopic mode and kinetic mode. Since BGK model is not based on the expansion in the mean free path in contrast to the Chapman-Enskog expansion, the solution can describe accurate relaxation of initial disturbance with any wavelength. This non-relativistic analysis gives suggestions for our next work of relativistic analysis of relaxation.

We investigate gas accretion onto a protoplanet, by considering the thermal effect of gas in three-dimensional hydrodynamical simulations, in which the wide region from a protoplanetary gas disk to a Jovian radius planet is resolved using the nested-grid method. We estimate the mass accretion rate and growth timescale of gas giant planets. The mass accretion rate increases with protoplanet mass for M_p < M_cri, while it becomes saturated or decreases for M_p > M_cri, where M_cri = 0.036 M_Jup (a_p/1AU)^0.75, and M_Jup and a_p are the Jovian mass and the orbital radius, respectively. The growth timescale of a gas giant planet or the timescale of the gas accretion onto the protoplanet is about 10^5 yr, that is two orders of magnitude shorter than the growth timescale of the solid core. The thermal effects barely affect the mass accretion rate because the gravitational energy dominates the thermal energy around the protoplanet. The mass accretion rate obtained in our local simulations agrees quantitatively well with those obtained in global simulations with coarser spatial resolution. The mass accretion rate is mainly determined by the protoplanet mass and the property of the protoplanetary disk. We find that the mass accretion rate is correctly calculated when the Hill or Bondi radius is sufficiently resolved. Using the oligarchic growth of protoplanets, we discuss the formation timescale of gas giant planets.

The formation and evolution of the circumstellar disk in unmagnetized molecular clouds is investigated using three-dimensional hydrodynamic simulations from the prestellar core until the end of the main accretion phase. In collapsing clouds, the first (adiabatic) core with a size of ~10AU forms prior to the formation of the protostar. At its formation, the first core has a thick disk-like structure, and is mainly supported by the thermal pressure. After the protostar formation, it decreases the thickness gradually, and becomes supported by the centrifugal force. We found that the first core is a precursor of the circumstellar disk. This indicates that the circumstellar disk is formed before the protostar formation with a size of ~10AU, which means that no protoplanetary disk smaller than <10AU exists. Reflecting the thermodynamics of the collapsing gas, at the protostar formation epoch, the circumstellar disk has a mass of ~0.01-0.1 solar mass, while the protostar has a mass of ~10^-3 solar mass. Thus, just after the protostar formation, the circumstellar disk is about 10-100 times more massive than the protostar. Even in the main accretion phase that lasts for ~10^5yr, the circumstellar disk mass dominates the protostellar mass. Such a massive disk is unstable to gravitational instability, and tends to show fragmentation. Our calculations indicate that the planet or brown-dwarf mass object may form in the circumstellar disk in the main accretion phase. In addition, the mass accretion rate onto the protostar shows strong time variability that is caused by the perturbation of proto-planets and/or the spiral arms in the circumstellar disk. Such variability provides a useful signature for detecting the planet-sized companion in the circumstellar disk around very young protostars.

Numerical simulations for the non-linear development of Kelvin-Helmholtz instability in two different density layers have been performed with the particle-based method (Godunov SPH) developed by Inutsuka (2002). The Godunov SPH can describe the Kelvin-Helmholtz instability even with a high density contrast, while the standard SPH shows the absence of the instability across a density gradient (Agertz et al. 2007). The interaction of a dense blob with a hot ambient medium has been performed also. The Godunov SPH describes the formation and evolution of the fingers due to the combinations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities. The blob test result coincides well with the results of the grid-based codes. An inaccurate handling of a density gradient in the standard SPH has been pointed out as the direct reason of the absence of the instabilities. An unphysical force happens at the density gradient even in a pressure equilibrium, and repulses particles from the initial density discontinuity. Therefore, the initial perturbation damps, and a gap forms at the discontinuity. The unphysical force has been studied in terms of the consistency of a numerical scheme. Contrary to the standard SPH, the momentum equation of the Godunov SPH doesnt use the particle approximation, and has been derived from the kernel convolution or a new Lagrangian function. The new Lagrangian function used in the Godunov SPH is more analogous to the real Lagrangian function for continuum. The momentum equation of the Godunov SPH has much better linear consistency, so the unphysical force is greatly reduced compared to the standard SPH in a high density contrast.

The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local $N$-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of non-axisymmetric structures, the creation of planetesimal seeds, and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to $L_y^{3/2}$, where $L_y$ is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of $L_x$, where $L_x$ is the size of the computational domain in the radial direction if $L_x$ is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.

Formation of interstellar clouds as a consequence of thermal instability is studied using two-dimensional two-fluid magnetohydrodynamic simulations. We consider the situation of converging, supersonic flows of warm neutral medium in the interstellar medium that generate a shocked slab of thermally unstable gas in which clouds form. We found, as speculated in paper I, that in the shocked slab magnetic pressure dominates thermal pressure and the thermal instability grows in the isochorically cooling, thermally unstable slab that leads formation of HI clouds whose number density is typically n < 100 cm^-3, even if the angle between magnetic field and converging flows is small. We also found that even if there is a large dispersion of magnetic field, evolution of the shocked slab is essentially determined by the angle between the mean magnetic field and converging flows. Thus, the direct formation of molecular clouds by piling up warm neutral medium does not seem a typical molecular cloud formation process, unless the direction of supersonic converging flows is biased to the orientation of mean magnetic field by some mechanism. However, when the angle is small, the HI shell generated as a result of converging flows is massive and possibly evolves into molecular clouds, provided gas in the massive HI shell is piled up again along the magnetic field line. We expect that another subsequent shock wave can pile up again the gas of the massive shell and produce a larger cloud. We thus emphasize the importance of multiple episodes of converging flows, as a typical formation process of molecular clouds.

We examine emission from a young protostellar object (YPO) with three-dimensional ideal MHD simulations and three-dimensional non-local thermodynamic equilibrium (non-LTE) line transfer calculations, and show the first results. To calculate the emission field, we employed a snapshot result of an MHD simulation having young bipolar outflows and a dense protostellar disk (a young circumstellar disk) embedded in an infalling envelope. Synthesized line emission of two molecular species (CO and SiO) show that subthermally excited SiO lines as a high density tracer can provide a better probe of the complex velocity field of a YPO, compared to fully thermalized CO lines. In a YPO at the earliest stage when the outflows are still embedded in the collapsing envelope, infall, rotation and outflow motions have similar speeds. We find that the combined velocity field of these components introduces a great complexity in the line emissions through varying optical thickness and emissivity, such as asymmetric double-horn profiles. We show that the rotation of the outflows, one of the features that characterizes an outflow driven by magneto-centrifugal forces, appears clearly in velocity channel maps and intensity-weighted mean velocity (first moment of velocity) maps. The somewhat irregular morphology of the line emission at this youngest stage is dissimilar to a more evolved object such as young Class 0. High angular resolution observation by e.g., the Atacama Large Millimeter/submillimeter Array (ALMA) telescope can reveal these features.

Brown dwarf formation and star formation efficiency are studied using a nested grid simulation that covers five orders of magnitude in spatial scale (10^4 - 0.1AU). Starting with a rotating magnetized compact cloud with a mass of 0.22 M_sun, we follow the cloud evolution until the end of main accretion phase. Outflow of about 5 km/s emerges about 100 yr before the protostar formation and does not disappear until the end of the calculation. The mass accretion rate declines from 10^-6 M_sun/yr to 10^-8 - 10^-12 M_sun/yr in a short time (about 10^4 yr) after the protostar formation. This is because (1) a large fraction of mass is ejected from the host cloud by the protostellar outflow and (2) the gas escapes from the host cloud by the thermal pressure. At the end of the calculation, 74% (167 M_Jup) of the total mass (225 M_Jup) is outflowing from the protostar, in which 34% (77 M_Jup) of the total mass is ejected by the protostellar outflow with supersonic velocity and 40% (90 M_Jup) escapes with subsonic velocity. On the other hand, 20% (45 M_Jup) is converted into the protostar and 6% (13 M_Jup) remains as the circumstellar disk. Thus, the star formation efficiency is epsilon = 0.2. The resultant protostellar mass is in the mass range of brown dwarfs. Our results indicate that brown dwarfs can be formed in compact cores in the same manner as hydrogen-burning stars, and the magnetic field and protostellar outflow are essential in determining the star formation efficiency and stellar mass.

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower-metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for < 10^-4 Z_sun. With initial mass < 0.1 M_sun, if fragments are ejected in embryo from the host clouds by multi-body interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.

The protostellar outflows have indispensable role in the formation of single stars, because they carry off the excess angular momentum from the centre of the shrinking gas cloud, and permits further collapse to form a star. On the other hand, a significant fraction of stars is supposed to be born as binaries with circumbinary disk that are frequently observed. Here, we investigate the evolution of a magnetized rotating cloud using three-dimensional resistive MHD nested-grid code, and show that the outflow is driven by the circumbinary disk and has an important role even in the binary formation. After the adiabatic core formation in the collapsing cloud core, the magnetic flux is significantly removed from the centre of the cloud by the Ohmic dissipation. Since this removal makes the magnetic braking ineffective, the adiabatic core continuously acquires the angular momentum to induce fragmentation and subsequent binary formation. The magnetic field accumulates in the circumbinary disk where the removal and accretion of magnetic field are balanced, and finally drives circumbinary outflow. This result explains the spectacular morphology of some specific young stellar objects such as L1551 IRS5. We can infer that most of the bipolar molecular outflows observed by low density tracers (i.e., CO) would correspond to circumbinary or circum-multiple outflows found in this report, since most of the young stellar objects are supposed to be binaries or multiples.

We investigate the driving mechanism of outflows and jets in star formation process using resistive MHD nested grid simulations. We found two distinct flows in the collapsing cloud core: Low-velocity outflows (sim 5 km/s) with a wide opening angle, driven from the first adiabatic core, and high-velocity jets (sim 50 km/s) with good collimation, driven from the protostar. High-velocity jets are enclosed by low-velocity outflow. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration. Therefore, the low-velocity outflow from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, high-velocity jet from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity outflows correspond to the observed molecular outflows, while high-velocity jets correspond to the observed optical jets. We suggest that the protostellar outflow and the jet are driven by different cores (the first adiabatic core and protostar), rather than that the outflow being entrained by the jet.

We investigate the effects of viscosity on disk-planet interaction and discuss how type I migration of planets is modified. We have performed a linear calculation using shearing-sheet approximation and obtained the detailed, high resolution density structure around the planet embedded in a viscous disk with a wide range of viscous coefficients. We use a time-dependent formalism that is useful in investigating the effects of various physical processes on disk-planet interaction. We find that the density structure in the vicinity of the planet is modified and the main contribution to the torque comes from this region, in contrast to inviscid case. Although it is not possible to derive total torque acting on the planet within the shearing-sheet approximation, the one-sided torque can be very different from the inviscid case, depending on the Reynolds number. This effect has been neglected so far but our results indicate that the interaction between a viscous disk and a planet can be qualitatively different from an inviscid case and the details of the density structure in the vicinity of the planet is critically important.

Instabilities of the dust layer in a protoplanetary disk are investigated. It is known that the streaming instability develops and dust density concentration occurs in a situation where the initial dust density is uniform. This work considers the effect of initial dust density gradient vertical to the midplane. Dust and gas are treated as different fluids. Pressure of dust fluid is assumed to be zero. The gas friction time is assumed to be constant. Axisymmetric two-dimensional numerical simulation was performed using the spectral method. We found that an instability develops with a growth rate on the order of the Keplerian angular velocity even if the gas friction time multiplied by the Keplerian angular velocity is as small as 0.001. This instability is powered by two sources: (1) the vertical shear of the azimuthal velocity, and (2) the relative motion of dust and gas coupled with the dust density fluctuation due to advection. This instability diffuses dust by turbulent advection and the maximum dust density decreases. This means that the dust concentration by the streaming instability which is seen in the case of a uniform initial dust density becomes ineffective as dust density gradient increases by the dust settling toward the midplane.

The radioactive isotope $^{60}$Fe ($T_{1/2} = 1.5 $ Myr) was present in the early solar system. It is unlikely that it was injected directly into the nascent solar system by a single, nearby supernova. It is proposed instead that it was inherited during the molecular cloud stage from several supernovae belonging to previous episodes of star formation. The expected abundance of $^{60}$Fe in star forming regions is estimated taking into account the stochasticity of the star-forming process, and it is showed that many molecular clouds are expected to contain $^{60}$Fe (and possibly $^{26}$Al [$T_{1/2} = 0.74 $ Myr]) at a level compatible with that of the nascent solar system. Therefore, no special explanation is needed to account for our solar system's formation.

We investigate the motion of a particle around a low mass planet embedded in a non-turbulent gaseous disk. We take into account the effect of the gas structure that is modified by the gravitational interaction between the planet. We derive an analytic formula that describes the change of the semi-major axis of the particle due to the encounter with the planet using local approximation in distant encounter regime. Our final formula includes the effects of steady, axisymmetric radial gas flow, the global gas pressure gradient in the disk, planet gravity, and the structure of the gas flow modified by the planet's gravity. We compare the analytic results with numerical calculations, and indicate that our formula well describes the secular evolution of the dust particles' semi-major axes well, especially for small particles with large drag coefficient. We discuss the conditions for dust gap opening around a low mass planet and radial distribution of dust particles. Our formula may provide a useful tool for calculating radial distribution of particles in a disk around the planet.

We examine MHD simulations of the propagation of a strong shock wave through the interstellar two-phase medium composed of small-scale cloudlets and diffuse warm neutral medium in two-dimensional geometry. The pre-shock two-phase medium is provided as a natural consequence of the thermal instability that is expected to be ubiquitous in the interstellar medium. We show that the shock-compressed shell becomes turbulent owing to the preshock density inhomogeneity and magnetic field amplification takes place in the shell. The maximum field strength is determined by the condition that plasma beta ~ 1, which gives the field strength on the order of 1 mG in the case of shock velocity ~ 1,000 km/s. The strongly magnetized region shows filamentary and knot-like structures in two-dimensional simulations. The spatial scale of the regions with magnetic field of 1 mG in our simulation is roughly 0.05 pc which is comparable to the spatial scale of the X-ray hot spots recently discovered in supernova remnants where the magnetic field strength is indicated to be amplified up to the order of 1 mG. This result may also suggest that the turbulent region with locally strong magnetic field is expected to be spread out in the region with frequent supernova explosions, such as in the Galactic center and starburst galaxies.

By performing local three-dimensional MHD simulations of stratified accretion disks, we investigate disk winds driven by MHD turbulence. Initially given weak vertical magnetic fields are effectively amplified by magnetorotational instability and winding due to differential rotation. Large scale channel flows develop most effectively at 1.5 - 2 times the scale heights where the magnetic pressure is comparable to but slightly smaller than the gas pressure. The breakup of these channel flows drives structured disk winds by transporting the Poynting flux to the gas. These features are universally observed in the simulations of various initial fields. This disk wind process should play an essential role in the dynamical evaporation of proto-planetary disks. The breakup of channel flows also excites the momentum fluxes associated with Alfvenic and (magneto-)sonic waves toward the mid-plane, which possibly contribute to the sedimentation of small dust grains in protoplanetary disks.

We simulate the acceleration processes of collisionless particles in a shock structure with magnetohydrodynamical (MHD) fluctuations. The electromagnetic field is represented as a sum of MHD shock solution ($\Mag_0, \Ele_0$) and torsional Alfven modes spectra ($\delta \Mag, \delta \Ele $). We represent fluctuation modes in logarithmic wavenumber space. Since the electromagnetic fields are represented analytically, our simulations can easily cover as large as eight orders of magnitude in resonant frequency, and do not suffer from spatial limitations of box size or grid spacing. We deterministically calculate the particle trajectories under the Lorenz force for time interval of up to ten years, with a time step of $\sim 0.5 \sec$. This is sufficient to resolve Larmor frequencies without a stochastic treatment. Simulations show that the efficiency of the first order Fermi acceleration can be parametrized by the fluctuation amplitude $\eta \equiv < \delta B^2 > ^{\frac 1 2} {B_0}^{-1}$ . Convergence of the numerical results is shown by increasing the number of wave modes in Fourier space while fixing $\eta$. Efficiency of the first order Fermi acceleration has a maximum at $ \eta \simeq 10^1$. The acceleration rate depends on the angle between the shock normal and $\Mag_0$, and is highest when the angle is zero. Our method will provide a convenient tool for comparing collisionless turbulence theories with, for example, observations of bipolar structure of super nova remnants (SNRs) and shell-like synchrotron-radiating structure.

We develop an unconditionally stable numerical method for solving the coupling between two fluids (frictional forces/heatings, ionization, and recombination), and investigate the dynamical condensation process of thermally unstable gas that is provided by the shock waves in a weakly ionized and magnetized interstellar medium by using two-dimensional two-fluid magnetohydrodynamical simulations. If we neglect the effect of magnetic field, it is known that condensation driven by thermal instability can generate high density clouds whose physical condition corresponds to molecular clouds (precursor of molecular clouds). In this paper, we study the effect of magnetic field on the evolution of supersonic converging HI flows and focus on the case in which the orientation of magnetic field to converging flows is orthogonal. We show that the magnetic pressure gradient parallel to the flows prevents the formation of high density and high column density clouds, but instead generates fragmented, filamentary HI clouds. With this restricted geometry, magnetic field drastically diminishes the opportunity of fast molecular cloud formation directly from the warm neutral medium, in contrast to the case without magnetic field.

Jet driving and fragmentation process in collapsing primordial cloud are studied using three-dimensional MHD nested grid simulations. Starting from a rotating magnetized spherical cloud with the number density of n=10^3 cm^-3, we follow the evolution of the cloud up to the stellar density n=10^22 cm^-3. We calculate 36 models parameterizing the initial magnetic and rotational energies (\gamma_0, \beta_0). In the collapsing primordial clouds, the cloud evolutions are characterized by the ratio of the initial rotational to magnetic energy, \gamma_0/\beta_0. The Lorentz force significantly affects the cloud evolution when \gamma_0 > \beta_0, while the centrifugal force is more dominant than the Lorentz force when \beta_0 > \gamma_0. When the cloud rotates rapidly with angular velocity of \Omega_0 > 10^-17 (n/10^3 cm^-3)^2/3 s^-1 and \beta_0 > \gamma_0, fragmentation occurs before the protostar is formed, but no jet appears after the protostar formation. On the other hand, a strong jet appears after the protostar formation without fragmentation when the initial cloud has the magnetic field of B_0 > 10^-9 (n/10^3 cm^-3)^2/3 G and \gamma_0 > \beta_0. Our results indicate that proto-Population III stars frequently show fragmentation and protostellar jet. Population III stars are therefore born as binary or multiple stellar systems, and they can drive strong jets, which disturb the interstellar medium significantly, as well as in the present-day star formation, and thus they may induce the formation of next generation stars.

We study the effect of poloidal magnetic field on type I planetary migration by linear perturbation analysis in the shearing-sheet approximation and the analytic results are compared with numerical calculations. In contrast to the unmagnetized case, the basic equations that describe the wake due to the planet in the disk allow magnetic resonances at which density perturbation diverges. In order to simplify the problem, we consider the case without magneto-rotational instability. We perform two sets of analyses: two-dimensional and three-dimensional. In two-dimensional analysis, we find the generalization of the torque formula previously known in unmagnetized case. In three-dimensional calculations, we focus on the disk with very strong magnetic field and derive a new analytic formula for the torque exerted on the planet. We find that when Alfven velocity is much larger than sound speed, two-dimensional torque is suppressed and three-dimensional modes dominate, in contrast to the unmagnetized case.

Improving our understanding of the initial conditions and earliest stages of star formation is crucial to gain insight into the origin of stellar masses, multiple systems, and protoplanetary disks. We review the properties of low-mass dense cores as derived from recent millimeter/submillimeter observations of nearby molecular clouds and discuss them in the context of various contemporary scenarios for cloud core formation and evolution. None of the extreme scenarios can explain all observations. Pure laminar ambipolar diffusion has relatively long growth times for typical ionization levels and has difficulty satisfying core lifetime constraints. Purely hydrodynamic pictures have trouble accounting for the inefficiency of core formation and the detailed velocity structure of individual cores. A possible favorable scenario is a mixed model involving gravitational fragmentation of turbulent molecular clouds close to magnetic criticality. The evolution of the magnetic field and angular momentum in individual cloud cores after the onset of gravitational collapse is also discussed. In particular, we stress the importance of radiation-magnetohydrodynamical processes and resistive MHD effects during the protostellar phase. We also emphasize the role of the formation of the short-lived first (protostellar) core in providing a chance for sub-fragmentation into binary systems and triggering MHD outflows. Future submillimeter facilities such as Herschel and ALMA will soon provide major new observational constraints in this field. On the theoretical side, an important challenge for the future will be to link the formation of molecular clouds and prestellar cores in a coherent picture.

We investigate the accretion of angular momentum onto a protoplanet system using three-dimensional hydrodynamical simulations. We consider a local region around a protoplanet in a protoplanetary disk with sufficient spatial resolution. We describe the structure of the gas flow onto and around the protoplanet in detail. We find that the gas flows onto the protoplanet system in the vertical direction crossing the shock front near the Hill radius of the protoplanet, which is qualitatively different from the picture established by two-dimensional simulations. The specific angular momentum of the gas accreted by the protoplanet system increases with the protoplanet mass. At Jovian orbit, when the protoplanet mass M_p is M_p < 1 M_J, where M_J is Jovian mass, the specific angular momentum increases as j \propto M_p. On the other hand, it increases as j \propto M_p^2/3 when the protoplanet mass is M_p > 1 M_J. The stronger dependence of the specific angular momentum on the protoplanet mass for M_p < 1 M_J is due to thermal pressure of the gas. The estimated total angular momentum of a system of a gas giant planet and a circumplanetary disk is two-orders of magnitude larger than those of the present gas giant planets in the solar system. A large fraction of the total angular momentum contributes to the formation of the circumplanetary disk. We also discuss the satellite formation from the circumplanetary disk.

The fragmentation process of primordial-gas cores during prestellar collapse is studied using three-dimensional nested-grid hydrodynamics. Starting from the initial central number density of n \sim10^3 cm^-3, we follow the evolution of rotating spherical cores up to the stellar density n \simeq 10^{22} cm^-3. An initial condition of the cores is specified by three parameters: the ratios of the rotation and thermal energies to the gravitational energy (\beta_0, and \alpha_0, respectively), and the amplitude of the bar-mode density perturbation (A_\phi). Cores with rotation \beta_0 > 10^{-6} are found to fragment during the collapse. The fragmentation condition hardly depends on either the initial thermal energy \alpha_0 or amplitude of bar-mode perturbation A_\phi. Since the critical rotation parameter for fragmentation is lower than that expected in the first star formation, binaries or multiples are also common for the first stars.

Fragmentation and binary formation processes are studied using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n=10^4 cm^-3) to the stellar core (n \simeq 10^22 cm^-3). We calculated 147 models with different initial magnetic, rotational, and thermal energies, and the amplitudes of the non-axisymmetric perturbation. In a collapsing cloud, fragmentation is mainly controlled by the initial ratio of the rotational to the magnetic energy, regardless of the initial thermal energy and amplitude of the non-axisymmetric perturbation. When the clouds have large rotational energies in relation to magnetic energies, fragmentation occurs in the low-density evolution phase (10^12 cm^-3 < n < 10^15 cm^-3) with separations of 3-300 AU. Fragments that appeared in this phase are expected to evolve into wide binary systems. On the other hand, fragmentation does not occur in the low-density evolution phase, when initial clouds have large magnetic energies in relation to the rotational energies. In these clouds, fragmentation only occurs in the high-density evolution phase (n > 10^17 cm^-3) after the clouds experience significant reduction of the magnetic field owing to Ohmic dissipation in the period of 10^12 cm^-3 < n < 10^15 cm^-3. Fragments appearing in this phase have separations of < 0.3 AU, and are expected to evolve into close binary systems. As a result, we found two typical fragmentation epochs, which cause different stellar separations. Although these typical separations are disturbed in the subsequent gas accretion phase, we might be able to observe two peaks of binary separations in extremely young stellar groups.

The evolution of the magnetic field and angular momentum in the collapsing cloud core is studied using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n=10^4 cm^-3) to the stellar core (n \simeq 10^22 cm^-3). The magnetic field strengths at the center of the clouds converge to a certain value as the clouds collapse, when the clouds have the same angular momenta but different strengths of the magnetic fields at the initial state. For 10^12 cm^-3 < n < 10^16 cm^-3, Ohmic dissipation considerably removes the magnetic field from the collapsing cloud core, and the magnetic field lines, which are strongly twisted for n <10^12 cm^-3, is de-collimated. The magnetic field lines are twisted and amplified again for nc > 10^16 cm^-3, because the magnetic field is recoupled with the warm gas. Finally, protostars at their formation epoch have 0.1-1kG of the magnetic fields, which are comparable to observations. The magnetic field strength of protostar slightly depends on the angular momentum of the host cloud. The protostar formed from the slowly rotating cloud core has a stronger magnetic field. The evolution of the angular momentum is closely related to the evolution of the magnetic field. The angular momentum in the collapsing cloud is removed by the magnetic effect. The formed protostars have 0.1-2 days of the rotation period at their formation epoch, which are slightly shorter than the observation. This indicates that the further removal mechanism of the angular momentum such as interaction between the protostar and disk, wind gas or jet is important in further evolution of the protostar.

The evolutionary conditions for the dissipative continuous magnetohydrodynamic (MHD) shocks are studied. We modify Hada's approach in the stability analysis of the MHD shock waves. The matching conditions between perturbed shock structure and asymptotic wave modes shows that all types of the MHD shocks, including the intermediate shocks, are evolutionary and perturbed solutions are uniquely defined. We also adopt our formalism to the MHD shocks in the system with resistivity without viscosity, which is often used in numerical simulation, and show that all types of shocks that are found in the system satisfy the evolutionary condition and perturbed solutions are uniquely defined. These results suggest that the intermediate shocks may appear in reality.

We investigate the dynamical condensation process in a magnetized thermally bistable medium. We perform one-dimensional two fluid numerical simulations that describe the neutral and ionized components in the interstellar medium with purely transverse magnetic fields. We find that the clouds that are formed as a consequence of the thermal instability always have a magnetic field strength on the order of a few microgauss irrespective of the initial strength. This shows good agreement with the measurements of the magnetic field strength in diffuse clouds. We show analytically that the final magnetic field strength in the clouds is determined by the balance between ambipolar diffusion and the accumulation of the magnetic field due to condensation.

The driving mechanism of jets and outflows in star formation process is studied using resistive MHD nested grid simulations. We calculated cloud evolution from the molecular cloud core to the stellar core. In the collapsing cloud core, we found two distinct flows: Low-velocity flows (sim 5 km/s) with a wide opening angle, driven from the adiabatic core, and high-velocity flows (sim 30 km/s) with good collimation, driven from the protostar. High-velocity flows are enclosed by low-velocity flows after protostar formation. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration; therefore, flows from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, flows from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity flows correspond to the observed molecular outflows, while high-velocity flows correspond to the observed optical jets. We suggest that the outflow and the jet are driven by different cores, rather than that the outflow being entrained by the jet.

We examine the triggering process of molecular cloud formation around diffuse HII regions. We calculate the time evolution of the shell as well as of the HII region in a two-phase neutral medium, solving the UV and FUV radiative transfer, the thermal and chemical processes in the time-dependent hydrodynamics code. In the cold neutral medium, the ambient gas is swept up in the cold (T \sim 30-40K) and dense (n \sim 10^3 cm^-3) shell around the HII region. In the shell, H_2 molecules are formed from the swept-up HI gas, but CO molecules are hardly formed. The reformation of H_2 molecules is more efficient with a higher-mass central star. The physical and chemical properties of gas in the shell are just intermediate between those of the neutral medium and molecular clouds observed by the CO emission. The dense shell with cold HI/H_2 gas easily becomes gravitationally unstable, and breaks up into small clouds. The cooling layer just behind the shock front also suffers from thermal instability, and will fragment into cloudlets with some translational motions. We suggest that the predicted cold ``dark'' HI/H_2 gas should be detected as the HI self-absorption (HISA) feature. We have sought such features in recent observational data, and found shell-like HISA features around the giant HII regions, W4 and W5. The shell-like HISA feature shows good spatial correlation with dust emission, but poor correlation with CO emission. Our quantitative analysis shows that the HISA cloud can be as cold as T \sim a few x 10K. (abridged)

We examine observational characteristics of multi-phase turbulent flows in the diffuse interstellar medium (ISM) using a synthetic radiation field of atomic and molecular lines. We consider the multi-phase ISM which is formed by thermal instability under the irradiation of UV photons with moderate visual extinction $A_V\sim 1$. Radiation field maps of C$^{+}$, C$^0$, and CO line emissions were generated by calculating the non-local thermodynamic equilibrium (nonLTE) level populations from the results of high resolution hydrodynamic simulations of diffuse ISM models. By analyzing synthetic radiation field of carbon lines of [\ion{C}{2}] 158 $\mu$m, [\ion{C}{1}] $^3P_2-^3P_1$ (809 GHz), $^3P_1-^3P_0$ (492 GHz), and CO rotational transitions, we found a high ratio between the lines of high- and low-excitation energies in the diffuse multi-phase interstellar medium. This shows that simultaneous observations of the lines of warm- and cold-gas tracers will be useful in examining the thermal structure, and hence the origin of diffuse interstellar clouds.

We performed N-body simulations of a dust layer without a gas component and examined the formation process of planetesimals. We found that the formation process of planetesimals can be divided into three stages: the formation of non-axisymmetric wake-like structures, the creation of aggregates, and the collisional growth of the aggregates. Finally, a few large aggregates and many small aggregates are formed. The mass of the largest aggregate is larger than the mass predicted by the linear perturbation theory. We examined the dependence of system parameters on the planetesimal formation. We found that the mass of the largest aggregates increase as the size of the computational domain increases. However the ratio of the aggregate mass to the total mass $M_\mathrm{aggr}/M_\mathrm{total}$ is almost constant $0.8-0.9$. The mass of the largest aggregate increases with the optical depth and the Hill radius of particles.

We analyze the structure and stability of the transition layer (or front) that connects the cold neutral medium and warm neutral medium in the plane-parallel geometry. Such fronts appear in recent numerical simulations of a thermally bistable interstellar medium. The front becomes an evaporation or condensation front depending on the surrounding pressure. The stability analysis is performed in both long- and short-wavelength approximations. We find that the plane-parallel evaporation front is unstable under corrugational deformations, whereas the condensation front seems to be stable. The instability is analogous to the Darrieus-Landau instability in combustion front. The growth rate of the instability is proportional to the speed of the evaporation flow and the corrugation wavenumber for modes with wavelength much longer than the thickness of the front, and it is suppressed at scales approximately equal to the thickness of the front. The timescale of the instability is smaller than the cooling timescale of the warm neutral medium ($\sim 1$ Myr), and can be as small as the cooling timescale of the cold neutral medium ($\sim 0.01-0.1$ Myr). Thus, this instability should be one of the processes for driving the interstellar turbulence.

We investigate outflows driven by a giant protoplanet using three-dimensional MHD nested grid simulations. We consider a local region around the protoplanet in the protoplanetary disk, and calculate three models: (a) unmagnetized disk model, (b) magnetized disk model having magnetic field azimuthally parallel to the disk, and (c) magnetic field perpendicular to the disk. Outflows with velocities, at least, 10 km/s are driven by the protoplanets in both magnetized disk models, while outflow does not appear in unmagnetized disk model. Tube-like outflows along the azimuthal direction of the protoplanetary disk appear in model with magnetic field being parallel to the disk. In this model, the magnetically dominated regions (i.e., density gap) are clearly contrasted from other regions and spiral waves appear near the protoplanet. On the other hand, in model with magnetic field being perpendicular to the disk, outflows are driven by a protoplanet with cone-like structure just as seen in the outflow driven by a protostar. Magnetic field lines are strongly twisted near the protoplanet and the outflows have well-collimated structures in this model.These outflows can be landmarks for searching exo-protoplanets in their formation stages. Our results indicate that the accretion rate onto the protoplanet tend to have a larger value than that expected from previous hydrodynamical calculations, since a fraction of the angular momentum of circum-planetary disk is removed by outflows, enhanced non-axisymmetric patterns caused by magnetic field, and magnetic braking. Possible implications for observation are also briefly discussed.

We quantitatively examine the significance of star formation triggered in the swept-up shell around an expanding HII region. If the swept-up molecular gas is sufficiently massive, new OB stars massive enough to repeat the triggering process will form in the shell. We determine the lower limit (M_thr) for the mass of the star that sweeps up the molecular gas, where at least one new star with mass M_* > M_thr forms after the shell fragmentation. To calculate the threshold stellar mass, M_thr, we examine how massive molecular shells can form around various central stars, by performing detailed numerical radiation hydrodynamics calculations. The mass of the photodissociated gas is generally larger than the mass of the photoionized gas. However, the swept-up molecular mass exceeds the photodissociated mass with a higher-mass star of M_* > 20 Msun. The accumulated molecular mass generally increases with the stellar mass, and amounts to 10^{4-5} Msun for M_* > 20 Msun with an ambient density of n \sim 100/cc. The threshold stellar mass is M_thr \sim 18 Msun with the star-formation efficiency of \epsilon \sim 0.1 and n \sim 100/cc. We examine the generality of this mode of run-away triggering for different sets of parameters, and found that M_thr \sim 15-20 Msun in various situations. If the ambient density is too high or the star-formation efficiency is too low, the triggering is not run-away, but a single event.

We show that the coronal heating and the fast solar wind acceleration in the coronal holes are natural consequence of the footpoint fluctuations of the magnetic fields at the photosphere, by performing one-dimensional magnetohydrodynamical simulation with radiative cooling and thermal conduction. We initially set up a static open flux tube with temperature 10^4K rooted at the photosphere. We impose transverse photospheric motions corresponding to the granulations with velocity <dv> = 0.7km/s and period between 20 seconds 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 initial cool and static atmosphere is naturally heated up to 10^6K and accelerated to 800km/s.

We show that the coronal heating and the acceleration of the fast solar wind in the coronal holes are natural consequence of the footpoint fluctuations of the magnetic fields at the photosphere by one-dimensional, time-dependent, and nonlinear magnetohydrodynamical simulation with radiative cooling and thermal conduction. We impose low-frequency (<0.05Hz) transverse photospheric motions, corresponding to the granulations, with velocity <dv> = 0.7$km/s. In spite of the attenuation in the chromosphere by the reflection, the sufficient energy of the generated outgoing Alfven waves transmit into the corona to heat and accelerate of the plasma by nonlinear dissipation. Our result clearly shows that the initial cool (10^4K) and static atmosphere is naturally heated up to 10^6K and accelerated to 800km/s, and explain recent SoHO observations and Interplanetary Scintillation measurements.

(abridged) We investigate how the properties of the corona and solar wind in the open coronal holes depend on the properties of the magnetic fields and their footpoint motions at the surface, by perfoming 1D MHD simulations from the photosphere to 0.3 or 0.1AU. We impose low-frequency (<0.05Hz) transverse fluctuations of the field lines at the photosphere with various amplitude, spectrum, and polarization in the open flux tubes with different photospheric field strength, B, and super-radial expansion of the cross section, f_max. We find that a transonic solar wind is the universal consequence. The atmosphere is also stably heated up to >10^6K by the dissipation of the Alfven waves through compressive-wave generation and wave reflection in the case of the sufficient wave input with photospheric amplitude, <dv> > 0.7km/s. The density, and accordingly the mass flux, of solar winds show a quite sensitive dependence on <dv> because of an unstable aspect of the heating by the nonlinear Alfven waves. A case with <dv>=0.4km/s gives ~50 times smaller mass flux than the fiducial case for the fast wind with <dv>=0.7km/s; solar wind almost disappears only if <dv> becomes half. We also find that the solar wind speed has a positive correlation with B/f_max, which is consistent with recent observations. We finally show that both fast and slow solar winds can be explained by the single process, the dissipation of the low-frequency Alfven waves, with different sets of <dv> and B/f_max. Our simulations naturally explain the observed (i) anticorrelation of the solar wind speed and the coronal temperature and (ii) larger amplitude of the Alfvenic fluctuations in the fast winds. In Appendix, we also explain our implementation of the outgoing boundary condition of the MHD waves with some numerical tests.

Supersonic turbulent motions are the remarkable properties of interstellar medium. Previous numerical simulations have demonstrated that the thermal instability in a shock-compressed layer produces the supersonic turbulent motion that does not decay. In this paper we focus on two- and three-dimensional numerical simulations of the non-linear development of simple thermal instability incorporating physical viscosity but without any external forcing, in order to isolate the effects of various processes responsible for the long-lasting turbulent motion. As the initial condition for our simulations, we set up spatially uniform gas with thermally unstable temperature in a box with periodic boundaries. After the linear growth stage of the thermal instability, two-phase medium forms where cold clumps are embedded in warm medium, and turbulent fluid flow clearly visible as translational motions of the cold clumps does not decay in a viscous dissipation timescale. The amplitude of the turbulent velocity increases when we reduce the Prandtl number that is the non-dimensional ratio of kinetic viscosity to thermal conduction: the saturation amplitude does not change when we increase the viscosity and thermal conduction coefficients simultaneously in order to keep the Prandtl number. This shows that the thermal conduction plays an important role in maintaining turbulent motions against viscous dissipation. The amplitude also increases when we increase the ratio of the computational domain length $L$ to the cooling length $\lambda_{\rm c}$ that is defined by the product of the cooling time and the sound speed, as long as $L \la 100 \lambda_{\rm c}$.

The protostellar jets driven by the formation of the first stars are studied by using MHD nested grid simulations. Starting from a slowly rotating spherical cloud of 5.1 times 10^4 Msun permeated by a uniform magnetic field, we follow the evolution from the central number density n = 10^3 cm^-3 to n simeq 10^23 cm^-3. We calculate four models that differ in initial magnetic field strengths and angular velocities. In all models, protostars of simeq 10^-3 Msun are formed at n simeq 10^22 cm^-3 in accordance with one-dimensional calculations. By this epoch, the magnetic flux density is amplified by 10 orders of magnitude from the initial value. Consequently, the formed protostar possesses the magnetic field of \sim 10^6 G that is much larger than the flux density of the present counterparts, reflecting the fact that the dissipation of a magnetic field is ineffective in primordial gas clouds. If the initial magnetic field B > 10^-9 (n/10^3 cm^-3)^2/3 G, the protostellar jet is launched and its velocities reaches \sim 70 km s^-1 by the time the protostellar mass becomes (4-6) times 10^-3 Msun, and a fraction (3-10%) of the accreting matter is blown off from the central region. Owing to the interaction of these ejecta with surrounding matter, expanding bow shocks are created at both heads of the jet. If this jet continues to sweep out the surrounding gas that otherwise accretes onto the central star or circumstellar disk, the final mass of the first star can be substantially reduced. In addition, dense post-shock regions behind the bow shocks are expected to promote the chemical reactions (formation of H_2 and HD), and this provides possible environments for subsequent low-mass star formation in the early universe.

Recent high resolution numerical simulations have suggested that the interstellar atomic hydrogen clouds have a complex two-phase structure. Since molecular clouds form through the contraction of HI gas, the question arises as to whether this structure is maintained in the molecular phase or not. Here we investigate whether the warm neutral atomic hydrogen (WNM) can exist in molecular clouds. We calculate how far a piece of WNM which is not heated by the UV photons could penetrate into the cloud, and find that in the absence of any heating it is unlikely that large fraction of WNM survives inside high pressure molecular clouds. We then consider two possible heating mechanisms, namely dissipation of turbulent energy and dissipation of MHD waves propagating in the WNM inside the cloud. We find that the second one is sufficient to allow the existence of WNM inside a molecular cloud of size $\simeq$ 1 pc having pressure equal to $\simeq 10 \times P_{\rm ISM}$. This result suggests the possibility that channels of magnetised WMN may provide efficient energy injection for sustaining internal turbulence which otherwise decays in a crossing time.

We analyze the stability of dust layer in protoplanetary disk to understand the effect of the relative motion between gas and dust. The previous analyses not including the effect of relative motion between gas and dust show that the shear-induced turbulence may prevent the dust grains from settling sufficiently to be gravitationally unstable. We determine the growth rate of Kelvin-Helmholtz instability in wide range of parameter space, and propose a possible path toward the planetesimal formation through the gravitational instability. We expect the density of dust layer becomes $\rho_\mathrm{d}/\rho_\mathrm{g} \sim 200$ if the dust grains can grow up to 5m.

We analyze the dynamical expansion of the HII region, photodissociation region, and the swept-up shell, solving the UV- and FUV-radiative transfer, the thermal and chemical processes in the time-dependent hydrodynamics code. Following our previous paper, we investigate the time evolutions with various ambient number densities and central stars. Our calculations show that basic evolution is qualitatively similar among our models with different parameters. The molecular gas is finally accumulated in the shell, and the gravitational fragmentation of the shell is generally expected. The quantitative differences among models are well understood with analytic scaling relations. The detailed physical and chemical structure of the shell is mainly determined by the incident FUV flux and the column density of the shell, which also follow the scaling relations. The time of shell-fragmentation, and the mass of the gathered molecular gas are sensitive tothe ambient number density. In the case of the lower number density, the shell-fragmentation occurs over a longer timescale, and the accumulated molecular gas is more massive. The variations with different central stars are more moderate. The time of the shell-fragmentation differs by a factor of several with the various stars of M_* = 12-101 M_sun. According to our numerical results, we conclude that the expanding HII region should be an efficient trigger for star formation in molecular clouds if the mass of the ambient molecular material is large enough.

Developing a theory of low-mass star formation ($\sim 0.1$ to 3~M$_{\odot}$) remains one of the most elusive and important goals of theoretical astrophysics. The star-formation process is the outcome of the complex dynamics of interstellar gas involving non-linear interactions of turbulence, gravity, magnetic field and radiation. The evolution of protostellar condensations, from the moment they are assembled by turbulent flows to the time they reach stellar densities, spans an enormous range of scales, resulting in a major computational challenge for simulations. Since the previous Protostars and Planets conference, dramatic advances in the development of new numerical algorithmic techniques have been successfully implemented on large scale parallel supercomputers. Among such techniques, Adaptive Mesh Refinement and Smooth Particle Hydrodynamics have provided frameworks to simulate the process of low-mass star formation with a very large dynamic range. It is now feasible to explore the turbulent fragmentation of molecular clouds and the gravitational collapse of cores into stars self-consistently within the same calculation. The increased sophistication of these powerful methods comes with substantial caveats associated with the use of the techniques and the interpretation of the numerical results. In this review, we examine what has been accomplished in the field and present a critique of both numerical methods and scientific results. We stress that computational simulations should obey the available observational constraints and demonstrate numerical convergence. Failing this, results of numerical simulations do not advance our understanding of low-mass star formation.

The stellar core formation and high speed jets driven by the formed core are studied by using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n = 10^6 cm^-3, r_c = 4.6 times 10^4 AU) to the stellar core (n \simeq 10^23 cm^-3, r_c \simeq 1 solar radius). We resolve cloud structure over 7 orders of magnitude in spatial extent and over 17 orders of magnitude in density contrast. For comparison, we calculate two models: resistive and ideal MHD models. Both models have the same initial condition, but the former includes dissipation process of magnetic field while the latter does not. The magnetic fluxes in resistive MHD model are extracted from the first core during 10^12 cm^-3 < n < 10^16 cm^-3 by Ohmic dissipation. Magnetic flux density of the formed stellar core (n \simeq 10^20 cm^-3) in resistive MHD model is two orders of magnitude smaller than that in ideal MHD model. Since magnetic braking is less effective in resistive MHD model, rapidly rotating stellar core (the second core) is formed. After stellar core formation, the magnetic field of the core is largely amplified both by magneto-rotational instability and the shearing motion between the stellar core and ambient medium. As a consequence, high speed (simeq 45 km,s^-1) jets are driven by the second core, which results in strong mass ejection. A cocoon-like structure around the second core also forms with clear bow shocks.

We estimate the evaporation timescale for spherical HI clouds consisting of the cold neutral medium surrounded by the warm neutral medium. We focus on clouds smaller than 1pc, which corresponds to tiny HI clouds recently discovered by Braun & Kanekar and Stanimirovi{\'c} & Heiles. By performing one-dimensional spherically symmetric numerical simulations of the two-phase interstellar medium (ISM), we derive the timescales as a function of the cloud size and of pressure of the ambient warm medium. We find that the evaporation timescale of the clouds of 0.01 pc is about 1Myr with standard ISM pressure, $p/k_{B}\sim 10^{3.5}$ K cm$^{-3}$, and for clouds larger than about 0.1 pc it depends strongly on the pressure. In high pressure cases, there exists a critical radius for clouds growing as a function of pressure, but the minimum critical size is $\sim$ 0.03 pc for a standard environment. If tiny HI clouds exist ubiquitously, our analysis suggests two implications: tiny HI clouds are formed continuously with the timescale of 1Myr, or the ambient pressure around the clouds is much higher than the standard ISM pressure. We also find that the results agree well with those obtained by assuming quasi-steady state evolution. The cloud-size dependence of the timescale is well explained by an analytic approximate formula derived by Nagashima, Koyama & Inutsuka. We also compare it with the evaporation rate given by McKee & Cowie.

We analyze the ionization state of the magnetohydrodynamically turbulent protoplanetary disks and propose a new mechanism of sustaining ionization. First, we show that in the quasi-steady state of turbulence driven by magnetorotational instability in a typical protoplanetary disk with dust grains, the amount of energy dissipation should be sufficient for providing the ionization energy that is required for activating magnetorotational instability. Second, we show that in the disk with dust grains the energetic electrons that compose electric currents in weakly ionized gas can provide collisional ionization, depending on the actual saturation state of magnetorotational turbulence. On the other hand, we show that in the protoplanetary disks with the reduced effect of dust grains, the turbulent motion can homogenize the ionization degree, leading to the activation of magnetorotational instability even in the absence of other ionization processes. The results in this Letter indicate that most of the regions in protoplanetary disks remain magnetically active, and we thus require a change in the theoretical modeling of planet formation.

We analyze the evaporation and condensation of spherical and cylindrical HI clouds of the cold neutral medium surrounded by the warm neutral medium. Because the interstellar medium including those two phases is well described as a thermally bistable fluid, it is useful to apply pattern formation theories to the dynamics of the interface between the two phases. Assuming isobaric evolution of fluids and a simple cubic form of the heat-loss function, we show the curvature effects of the interface. We find that approximate solutions for spherical clouds are in good agreement with numerically obtained solutions. We extend our analysis to general curved fronts taking into account the curvature effects explicitly. We find that the curvature effects always stabilise curved interfaces under assumptions such as isobaric evolution we adopt in this Letter.

We analyze the dynamical expansion of the HII region and outer photodissociation region (PDR) around a massive star by solving the UV and FUV radiation transfer and the thermal and chemical processes in a time-dependent hydrodynamics code. We focus on the physical structure of the shell swept up by the shock front (SF) preceding the ionization front (IF). After the IF reaches the initial Stromgren radius, the SF emerges in front of the IF and the geometrically thin shell bounded with the IF and the SF is formed. The gas density inside the shell is about 10-100 times as high as the ambient gas density. Initially the dissociation fronts expands faster than IF and the PDR is formed outside the HII region. Thereafter the IF and SF gradually overtakes the proceeding dissociation fronts (DFs), and eventually DFs are taken in the shell. The chemical composition within the shell is initially atomic, but hydrogen and carbon monoxide molecules are gradually formed. This is partly because the IF and SF overtake DFs and SF enters the molecular region, and partly because the reformation timescales of the molecules become shorter than the dynamical timescale. The gas shell becomes dominated by the molecular gas by the time of gravitational fragmentation, which agrees with some recent observations. A simple estimation of star formation rate in the shell can provide a significant star formation rate in our galaxy.

We present the dynamics of a thermally bistable medium using one-dimensional numerical calculations, including cooling, heating, thermal conduction, and physical viscosity.We set up a two-phase medium from a thermally unstable one-phase medium and follow its long-term evolution. To clarify the role of thermal conduction, we compare the results of the two models, with and without thermal conduction. The calculations show that the thermal conduction helps to generate the kinetic energy of translational motions of the clouds. Next, we focus on spatial resolution because we have to resolve the Field length $\lambda_{\rm F}$, which is the characteristic length scale of the thermal conduction. The results show convergence only when thermal conduction is included and a large enough cell number is used. We find it necessary to maintain a cell size of less than $\lambda_{\rm F}/3$ to achieve a converged motion in the two-phase medium. We refer to the constraint that $\lambda_{\rm F}/3$ be resolved as the ``Field condition''. The inclusion of thermal conduction to satisfy the Field condition is crucial to numerical studies of thermal instability (TI) and may be important for studies of the turbulent interstellar two-phase medium: the calculations of TI without thermal conduction may be susceptible to contamination by artificial phenomena that do not converge with increasing resolutions.

The saturation level of the magnetorotational instability (MRI) is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted and the vertical component of gravity is ignored, so that the evolution of the MRI is followed in a small local part of the disk. We focus on the dependence of the saturation level of the stress on the gas pressure, which is a key assumption in the standard alpha disk model. From our numerical experiments it is found that there is a weak power-law relation between the saturation level of the Maxwell stress and the gas pressure in the nonlinear regime; the higher the gas pressure, the larger the stress. Although the power-law index depends slightly on the initial field geometry, the relationship between stress and gas pressure is independent of the initial field strength, and is unaffected by Ohmic dissipation if the magnetic Reynolds number is at least 10. The relationship is the same in adiabatic calculations, where pressure increases over time, and nearly-isothermal calculations, where pressure varies little with time. Our numerical results are qualitatively consistent with an idea that the saturation level of the MRI is determined by a balance between the growth of the MRI and the dissipation of the field through reconnection. The quantitative interpretation of the pressure-stress relation, however, may require advances in the theoretical understanding of non-steady magnetic reconnection.

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.

We have investigated the evolution and distribution of molecules in collapsing prestellar cores via numerical chemical models, adopting the Larson-Penston solution and its delayed analogues to study collapse. Molecular abundances and distributions in a collapsing core are determined by the balance among the dynamical, chemical and adsorption time scales. When the central density n_H of a prestellar core with the Larson-Penston flow rises to 3 10^6 cm^{-3}, the CCS and CO column densities are calculated to show central holes of radius 7000 AU and 4000 AU, respectively, while the column density of N2H+ is centrally peaked. These predictions are consistent with observations of L1544. If the dynamical time scale of the core is larger than that of the Larson-Penston solution owing to magnetic fields, rotation, or turbulence, the column densities of CO and CCS are smaller, and their holes are larger than in the Larson-Penston core with the same central gas density. On the other hand, N2H+ and NH3 are more abundant in the more slowly collapsing core. Therefore, molecular distributions can probe the collapse time scale of prestellar cores. Deuterium fractionation has also been studied via numerical calculations. The deuterium fraction in molecules increases as a core evolves and molecular depletion onto grains proceeds. When the central density of the core is n_H=3 10^6 cm^{-3}, the ratio DCO+/HCO+ at the center is in the range 0.06-0.27, depending on the collapse time scale and adsorption energy; this range is in reasonable agreement with the observed value in L1544.

The propagation of a shock wave into an interstellar medium is investigated by two-dimensional numerical hydrodynamic calculation with cooling, heating and thermal conduction. We present results of the high-resolution two-dimensional calculations to follow the fragmentation due to the thermal instability in a shock-compressed layer. We find that geometrically thin cooling layer behind the shock front fragments into small cloudlets. The cloudlets have supersonic velocity dispersion in the warm neutral medium in which the fragments are embedded as cold condensations. The fragments tend to coalesce and become larger clouds.

We analytically investigate the formation of an HII region in the accreting envelope of a newborn star. Special care is taken to examine the role of ionizing radiation force. This effect modifies velocity and density distributions and thereby affects the expansion of the HII region. As a result, the upper limit of the stellar mass imposed by the growth of an HII region around a forming star is increased by a larger factor than the previous estimate. In particular, for a star forming out of metal-free gas, this mechanism does not impose a firm upper limit on its mass.

The nonlinear evolution and the saturation mechanism of the magnetorotational instability (MRI) are investigated using three-dimensional resistive MHD simulations. A local shearing box is used for our numerical analysis and the simulations begin with a purely vertical magnetic field. We find that the magnetic stress in the nonlinear stage of the MRI is strongly fluctuating. The time evolution shows the quasi-periodic recurrence of spike-shape variations typically for a few orbits which correspond to the rapid amplification of the magnetic field by the nonlinear growth of a two-channel solution followed by the decay through magnetic reconnections. The increase rate of the total energy in the shearing box system is analytically related to the volume-averaged torque in the system. We find that at the saturated state this energy gain of the system is balanced with the increase of the thermal energy mostly due to the joule heating. The spike-shape time evolution is a general feature of the nonlinear evolution of the MRI in the disks threaded by vertical fields and can be seen if the effective magnetic Reynolds number is larger than about unity.

We investigate the propagation of a shock wave into a warm neutral medium and cold neutral medium by one-dimensional hydrodynamic calculations with detailed treatment of thermal and chemical processes. Our main result shows that thermal instability inside the shock-compressed layer produces a geometrically thin, dense layer in which a large amount of hydrogen molecules is formed. Linear stability analysis suggests that this thermally-collapsed layer will fragment into small molecular cloudlets. We expect that frequent compression due to supernova explosions, stellar winds, spiral density waves, etc., in the galaxy make the interstellar medium occupied by these small molecular cloudlets.