Source author record

Shu-ichiro Inutsuka

Shu-ichiro Inutsuka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.EP astro-ph.GA astro-ph.HE astro-ph physics.plasm-ph astro-ph.IM physics.flu-dyn astro-ph.CO nucl-th hep-ph nlin.PS nucl-ex physics.comp-ph

Catalog footprint

What is connected

82works

14topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Effects of magnetic field orientations in dense cores on gas kinematics in protostellar envelopes

Theoretically, misalignment between the magnetic field and rotational axis in a dense core is considered to be dynamically important in the star formation process, however, extent of this influence remains observationally unclear. For a sample of 32 Class 0 and I protostars in the Perseus Molecular Cloud, we analyzed gas motions using C$^{18}$O data from the SMA MASSES survey and the magnetic field structures using 850 $μ$m polarimetric data from the JCMT BISTRO-1 survey and archive. We do not find any significant correlation between the velocity gradients in the C$^{18}$O emission in the protostellar envelopes at a 1,000 au scale and the misalignment between the outflows and magnetic field orientations in the dense cores at a 4,000 au scale, and there is also no correlation between the velocity gradients and the angular dispersions of the magnetic fields. However, a significant dependence on the misalignment angles emerges after we normalize the rotational motion by the infalling motion, where the ratios increase from $\lesssim1$ to $\gtrsim1$ with increasing misalignment angles. This suggests that the misalignment could prompt angular momentum transportation to the envelope scale but is not a dominant factor in determining the envelope rotation, and other parameters, like mass accretion in protostellar sources, also play an important role. These results remain valid after taking into account projection effects. The comparison between our estimated angular momentum in the protostellar envelopes and the sizes of the known protostellar disks suggests that significant angular momentum is likely lost between radii of $\sim$1,000-100 au in protostellar envelopes.

preprint2022arXiv

Velocity structure of the 50 pc-long NGC 6334 filamentary cloud: Hints of multiple compressions and their impact on the cloud properties?

[Abridged] The interstellar medium is observed to be organised in filamentary structures, as well as neutral (HI) and ionized (HII) bubbles. The expanding nature of these bubbles makes them shape their surroundings and possibly play a role in the formation and evolution of interstellar filaments. We present APEX $^{13}$CO and C$^{18}$O(2-1) observations of the NGC 6334 molecular cloud. We investigate the gas velocity structure along and across the 50 pc-long cloud and towards 75 identified velocity-coherent-filaments (VCFs). We measure a wealth of velocity gradients along the VCFs. We derive the column density and velocity power spectra of the VCFs. These power spectra are well represented with power laws showing similar slopes for both quantities (with a mean of about -2), albeit some differ by up to a factor of two. The position velocity diagrams perpendicular to three VCFs show the V-shaped velocity pattern, corresponding to a bent structure in velocity space with the filament at the tip of the V surrounded by an extended structure connected to it with a velocity gradient. This velocity structure is qualitatively similar to that resulting from numerical simulations of filament formation from large-scale compression from propagating shock fronts. In addition, the radial profiles perpendicular to these VCFs hint to small-scale internal impacts from neighbouring HII bubbles. The observed opposite curvature in velocity space towards the VCFs points to various origins of large-scale external compressions from propagating HI bubbles. This suggests the plausible importance of multiple HI compressions, separated in space and time, in the formation and evolution of molecular clouds and their star formation history. These latter atomic compressions due to past and distant star formation events are complemented by the impact of HII bubbles from present time and local star formation activity.

preprint2021arXiv

The Effects of Cosmic-Ray Diffusion and Radiative Cooling on the Galactic Wind from the Milky Way

The effects of cosmic-ray diffusion and radiative cooling on the structure of the Galactic wind are studied using a steady state approximation. It is known that realistic cooling processes suppress the wind from launching. The effects of cosmic-ray diffusion are also supposed to be unfavorable for launching the wind. Both of these effects have not been studied simultaneously in a steady-state approximation of the wind. We find 327,254 solutions of the steady-state Galactic wind and confirm that: the effect of cosmic-ray pressure depends on the Alfv{é}n Mach number, the mass flux carried by the wind does not depend on the cosmic-ray pressure directly (but depends on the thermal pressure), and the typical conditions found in the Galaxy may correspond to the wind solution that provides metal polluted matters at a height of $\sim300$~kpc from the disk.

preprint2020arXiv

A low-velocity bipolar outflow from a deeply embedded object in Taurus revealed by the Atacama Compact Array

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 Atacama Large Millimeter/submillimeter Array$-$Atacama Compact Array (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) emitters in our sample. All of 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.

preprint2020arXiv

Distribution and kinematics of 26Al in the Galactic disc

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.

preprint2020arXiv

FRagmentation and Evolution of dense cores Judged by ALMA (FREJA). I (Overview). Inner $\sim$1000 au structures of prestellar/protostellar cores in Taurus

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 Atacama Large Millimeter/submillimeter Array-Atacama Compact Array (ALMA-ACA) 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 occur in a region smaller than 0.1 pc, which may determine the final core mass associated with individual protostar formation before starting the dynamical collapse of the core with central density of $\sim$(0.3-1) $\times$ 10$^6$ cm$^{-3}$.

preprint2020arXiv

From diffuse gas to dense molecular cloud cores

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.

preprint2020arXiv

Secular Gravitational Instability of Drifting Dust in Protoplanetary Disks: Formation of Dusty Rings without Significant Gas Substructures

Secular gravitational instability (GI) is one of the promising mechanisms for creating annular substructures and planetesimals in protoplanetary disks. We perform numerical simulations of the secular GI in a radially extended disk with inward drifting dust grains. The results show that, even in the presence of the dust diffusion, the dust rings form via the secular GI while the dust grains are moving inward, and the dust surface density increases by a factor of ten. Once the secular GI develops into a nonlinear regime, the total mass of the resultant rings can be a significant fraction of the dust disk mass. In this way, a large amount of drifting dust grains can be collected in the dusty rings and stored for planetesimal formation. In contrast to the emergence of remarkable dust substructures, the secular GI does not create significant gas substructures. This result indicates that observations of a gas density profile near the disk midplane enable us to distinguish the mechanisms for creating the annular substructures in the observed disks. The resultant rings start decaying once they enter the inner region stable to the secular GI. Since the ring-gap contrast smoothly decreases, it seems possible that the rings are observed even in the stable region. We also discuss the likely outcome of the non-linear growth and indicate the possibility that a significantly developed region of the secular GI may appear as a gap-like substructure in dust continuum emission since dust growth into larger solid bodies and planetesimal formation reduce the total emissivity.

preprint2020arXiv

The molecular cloud lifecycle

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.

preprint2020arXiv

The Origin of the Stellar Mass Distribution and Multiplicity

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.

preprint2019arXiv

High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions

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.

preprint2016arXiv

A revised condition for self-gravitational fragmentation of protoplanetary disks

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.

preprint2016arXiv

An Origin of Multiple Ring Structure and Hidden Planets in HL Tau: A Unified Picture by Secular Gravitational Instability

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 $α$ 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.

preprint2016arXiv

Dust Dynamics in Protoplanetary Disk Winds Driven by Magneto-Rotational Turbulence: A Mechanism for Floating Dust Grains with Characteristic Sizes

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 $μ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.

preprint2016arXiv

Formation of HI Clouds in Shock-compressed Interstellar Medium: Physical Origin of Angular Correlation Between Filamentary Structure and Magnetic Field

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.

preprint2016arXiv

Revealing a detailed mass distribution of a high-density core MC27/L1521F in Taurus with ALMA

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.

preprint2016arXiv

Stochastic Particle Acceleration in Turbulence Generated by the Magnetorotational Instability

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.

preprint2016arXiv

The Origin of the Rotation Profiles in Star Forming Clouds

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.

preprint2015arXiv

An Extension of Godunov SPH: Application to Negative Pressure Media

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.

preprint2015arXiv

An origin of arc structures deeply embedded in dense molecular cloud cores

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.

preprint2015arXiv

Atmospheric Escape by Magnetically Driven Wind from Gaseous Planets II --Effects of Magnetic Diffusion--

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.

preprint2015arXiv

Development of a Method for the Observation of Lightning in Protoplanetary Disks Using Ion Lines

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π$ in azimuth and $25 \mathrm{au}<r<50 \mathrm{au}$ is $1200σ$- to $4000σ$-detectable. The lower limits of the radii of $5σ$-detectable LMG clumps are between 1.6 au and 5.3 au, depending on the models.

preprint2015arXiv

Diffusion of Cosmic Rays in a Multiphase Interstellar Medium Swept-Up by a Supernova Remnant Blast Wave

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.

preprint2015arXiv

Effects of Ohmic and ambipolar diffusion on the formation and evolution of the first cores, protostars and circumstellar discs

We investigate the formation and evolution of a first core, protostar, and circumstellar disc with a three-dimensional non-ideal (including both Ohmic and ambipolar diffusion) radiation magnetohydrodynamics simulation. We found that the magnetic flux is largely removed by magnetic diffusion in the first core phase and that the plasma $β$ of the centre of the first core becomes large, $β>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, $β\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 $μ=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.

preprint2015arXiv

Formation of terrestrial planets in disks evolving via disk winds and implications for the origin of the solar system's terrestrial planets

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.

preprint2015arXiv

On the Radius Anomaly of Hot Jupiters: Reexamination of the Possibility and Impact of Layered Convection

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.

preprint2015arXiv

Significant Gas-to-Dust Ratio Asymmetry and Variation in the Disk of HD 142527 and the Indication of Gas Depletion

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.

preprint2015arXiv

The Formation and Destruction of Molecular Clouds and Galactic Star Formation

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.

preprint2014arXiv

ALMA observations of a high-density core in Taurus: dynamical gas interaction at the possible site of a multiple star formation

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.

preprint2014arXiv

Atmospheric escape by magnetically driven wind from gaseous planets

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.

preprint2014arXiv

Effects of radiation transfer on the structure of self-gravitating disks, their fragmentation and evolution of the fragments

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 $Ω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.

preprint2014arXiv

MHD Simulations of Global Accretion Disks with Vertical Magnetic Fields

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

preprint2014arXiv

N-body Simulations of Terrestrial Planet Formation under the Influence of a Hot Jupiter

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.

preprint2014arXiv

On the Viability of the Magnetorotational Instability in Circumplanetary Disks

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.

preprint2014arXiv

Self-Sustained Turbulence without Dynamical Forcing: A Two-Dimensional Study of a Bistable Interstellar Medium

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.

preprint2014arXiv

The Nonlinear Ohm's Law: Plasma Heating by Strong Electric Fields and its Effects on the Ionization Balance in Protoplanetary Disks

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.

preprint2014arXiv

Two-Component Secular Gravitational Instability in a Protoplanetary Disk: A Possible Mechanism for Creating Ring-Like Structures

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.

preprint2013arXiv

A new scheme of causal viscous hydrodynamics for relativistic heavy-ion collisions: A Riemann solver for quark-gluon plasma

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.

preprint2013arXiv

A Semi-Analytical Description for the Formation and Gravitational Evolution of Protoplanetary Disks

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.

preprint2013arXiv

An explicit scheme for ohmic dissipation with smoothed particle magneto-hydrodynamics

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 νsts ~ 0.01 and Nsts ~ 5.

preprint2013arXiv

Conditions for Circumstellar Disk Formation: Effects of Initial Conditions and Sink Treatment

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.

preprint2013arXiv

Crowding Out of Giants by Dwarfs: An Origin for the Lack of Companion Planets in Hot Jupiter Systems

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.

preprint2013arXiv

From Filamentary Networks to Dense Cores in Molecular Clouds: Toward a New Paradigm for Star Formation

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.

preprint2013arXiv

Interdependence of Electric Discharge and The Magnetorotational Instability in Protoplanetary Disks

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.

preprint2013arXiv

Local Enhancement of Surface Density in the Protoplanetary Ring Surrounding HD 142527

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.

preprint2012arXiv

Enhanced Dissipation Rate of Magnetic Field in Striped Pulsar Winds by the Effect of Turbulence

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 σ-problem in the Crab pulsar wind nebulae.

preprint2012arXiv

Formation of Turbulent and Magnetized Molecular Clouds via Accretion Flows of HI Clouds

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.

preprint2012arXiv

Structure of dynamical condensation fronts in the interstellar medium

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.

preprint2012arXiv

Toward Understanding the Origin of Turbulence in Molecular Clouds: Small Scale Structures as Units of Dynamical Multi-Phase Interstellar Medium

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

preprint2011arXiv

A Fast and Accurate Calculation Scheme for Ionization Degrees in Protoplanetary and Circumplanetary Disks with Charged Dust Grains

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.

preprint2011arXiv

A Fast Numerical Scheme for Causal Relativistic Hydrodynamics with Dissipation

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.

preprint2011arXiv

Effect of Magnetic Braking on the Circumstellar Disk Formation in a Strongly Magnetized Cloud

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.

preprint2011arXiv

Gravitational Fragmentation of Expanding Shells. I. Linear Analysis

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.

preprint2011arXiv

Gravitational Fragmentation of Expanding Shells. II. Three-dimensional Simulations

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.

preprint2011arXiv

Recurrent Planet Formation and Intermittent Protostellar Outflows Induced by Episodic Mass Accretion

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.

preprint2011arXiv

Secular Gravitational Instability of a Dust Layer in Shear Turbulence

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.

preprint2011arXiv

Smoothed particle magnetohydrodynamics with a Riemann solver and the method of characteristics

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{é}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.

preprint2011arXiv

Toward Understanding the Cosmic-ray Acceleration at Young Supernova Remnants Interacting with Interstellar Clouds: Possible Applications to RX J1713.7-3946

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 $π^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.

preprint2010arXiv

Analysis of the Relaxation Process using Non-Relativistic Kinetic Equation

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.

preprint2010arXiv

Emergence of Protoplanetary Disks and Successive Formation of Gaseous Planets by Gravitational Instability

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.

preprint2010arXiv

Formation Process of the Circumstellar Disk: Long-term Simulations in the Main Accretion Phase of Star Formation

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.

preprint2010arXiv

Gas Accretion onto a Protoplanet and Formation of a Gas Giant Planet

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.

preprint2010arXiv

N-Body Simulation of Planetesimal Formation through Gravitational Instability of a Dust Layer in Laminar Gas Disk

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.

preprint2010arXiv

Protoplanetary Disk Winds by Magnetorotational Instability : Formation of an Inner Hole and a Crucial Assist for Planet Formation

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.

preprint2010arXiv

The relativistic kinetic dispersion relation: Comparison of the relativistic Bhatnagar-Gross-Krook model and Grad's 14-moment expansion

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.

preprint2010arXiv

Two-Dimensional Study of the Propagation of Planetary Wake and the Indication to Gap Opening in an Inviscid Protoplanetary Disk

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.

preprint2010arXiv

Two-Step Acceleration Model of Cosmic Rays at Middle-Aged SNR

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.

preprint2009arXiv

Binary Formation with Different Metallicities: Dependence on Initial Conditions

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.

preprint2009arXiv

First Direct Simulation of Brown Dwarf Formation in a Compact Cloud Core

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.

preprint2009arXiv

Kelvin-Helmholtz instabilities with Godunov SPH

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.

preprint2009arXiv

N-Body Simulation of Planetesimal Formation through Gravitational Instability and Coagulation. II. Accretion Model

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.

preprint2009arXiv

Orbital Evolution of a Particle Interacting with a Single Planet in a Protoplanetary Disk

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.

preprint2009arXiv

Protostellar Jet and Outflow in the Collapsing Cloud Core

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.

preprint2009arXiv

The Circumbinary Outflow: A Protostellar Outflow Driven by a Circumbinary Disk

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.

preprint2009arXiv

Turbulence and Magnetic Field Amplification in Supernova Remnants: Interactions Between A Strong Shock Wave and Multi-Phase Interstellar Medium

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.

preprint2009arXiv

Two-Fluid MHD Simulations of Converging Hi Flows in the Interstellar Medium. II: Are Molecular Clouds Generated Directly from Warm Neutral Medium?

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.

preprint2008arXiv

Two-Fluid MHD Simulations of Converging HI Flows in the Interstellar Medium. I: Methodology and Basic Results

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.

preprint2007arXiv

Dynamical Formation of the Dark Molecular Hydrogen Clouds around Diffuse HII Regions

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)

preprint2006arXiv

Dynamical Expansion of Ionization and Dissociation Front around a Massive Star : A Starburst Mechanism

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 ε\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.

preprint2006arXiv

Outflows driven by Giant Protoplanets

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.

preprint2006arXiv

Solar Winds Driven by Nonlinear Low-Frequency Alfven Waves from the Photosphere : Parametric Study for Fast/Slow Winds and Disappearance of Solar Winds

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

Shu-ichiro Inutsuka

What is connected

Connect this record

See the researcher in context

Building this map preview

82 published item(s)

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

Effects of magnetic field orientations in dense cores on gas kinematics in protostellar envelopes

Velocity structure of the 50 pc-long NGC 6334 filamentary cloud: Hints of multiple compressions and their impact on the cloud properties?

The Effects of Cosmic-Ray Diffusion and Radiative Cooling on the Galactic Wind from the Milky Way

A low-velocity bipolar outflow from a deeply embedded object in Taurus revealed by the Atacama Compact Array

Distribution and kinematics of 26Al in the Galactic disc

FRagmentation and Evolution of dense cores Judged by ALMA (FREJA). I (Overview). Inner $\sim$1000 au structures of prestellar/protostellar cores in Taurus

From diffuse gas to dense molecular cloud cores

Secular Gravitational Instability of Drifting Dust in Protoplanetary Disks: Formation of Dusty Rings without Significant Gas Substructures

The molecular cloud lifecycle

The Origin of the Stellar Mass Distribution and Multiplicity

High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions

A revised condition for self-gravitational fragmentation of protoplanetary disks

An Origin of Multiple Ring Structure and Hidden Planets in HL Tau: A Unified Picture by Secular Gravitational Instability

Dust Dynamics in Protoplanetary Disk Winds Driven by Magneto-Rotational Turbulence: A Mechanism for Floating Dust Grains with Characteristic Sizes

Formation of HI Clouds in Shock-compressed Interstellar Medium: Physical Origin of Angular Correlation Between Filamentary Structure and Magnetic Field

Revealing a detailed mass distribution of a high-density core MC27/L1521F in Taurus with ALMA

Stochastic Particle Acceleration in Turbulence Generated by the Magnetorotational Instability

The Origin of the Rotation Profiles in Star Forming Clouds

An Extension of Godunov SPH: Application to Negative Pressure Media

An origin of arc structures deeply embedded in dense molecular cloud cores

Atmospheric Escape by Magnetically Driven Wind from Gaseous Planets II --Effects of Magnetic Diffusion--

Development of a Method for the Observation of Lightning in Protoplanetary Disks Using Ion Lines

Diffusion of Cosmic Rays in a Multiphase Interstellar Medium Swept-Up by a Supernova Remnant Blast Wave

Effects of Ohmic and ambipolar diffusion on the formation and evolution of the first cores, protostars and circumstellar discs

Formation of terrestrial planets in disks evolving via disk winds and implications for the origin of the solar system's terrestrial planets

On the Radius Anomaly of Hot Jupiters: Reexamination of the Possibility and Impact of Layered Convection

Significant Gas-to-Dust Ratio Asymmetry and Variation in the Disk of HD 142527 and the Indication of Gas Depletion

The Formation and Destruction of Molecular Clouds and Galactic Star Formation

ALMA observations of a high-density core in Taurus: dynamical gas interaction at the possible site of a multiple star formation

Atmospheric escape by magnetically driven wind from gaseous planets

Effects of radiation transfer on the structure of self-gravitating disks, their fragmentation and evolution of the fragments

MHD Simulations of Global Accretion Disks with Vertical Magnetic Fields

N-body Simulations of Terrestrial Planet Formation under the Influence of a Hot Jupiter

On the Viability of the Magnetorotational Instability in Circumplanetary Disks

Self-Sustained Turbulence without Dynamical Forcing: A Two-Dimensional Study of a Bistable Interstellar Medium

The Nonlinear Ohm's Law: Plasma Heating by Strong Electric Fields and its Effects on the Ionization Balance in Protoplanetary Disks

Two-Component Secular Gravitational Instability in a Protoplanetary Disk: A Possible Mechanism for Creating Ring-Like Structures

A new scheme of causal viscous hydrodynamics for relativistic heavy-ion collisions: A Riemann solver for quark-gluon plasma

A Semi-Analytical Description for the Formation and Gravitational Evolution of Protoplanetary Disks

An explicit scheme for ohmic dissipation with smoothed particle magneto-hydrodynamics

Conditions for Circumstellar Disk Formation: Effects of Initial Conditions and Sink Treatment

Crowding Out of Giants by Dwarfs: An Origin for the Lack of Companion Planets in Hot Jupiter Systems

From Filamentary Networks to Dense Cores in Molecular Clouds: Toward a New Paradigm for Star Formation

Interdependence of Electric Discharge and The Magnetorotational Instability in Protoplanetary Disks

Local Enhancement of Surface Density in the Protoplanetary Ring Surrounding HD 142527

Enhanced Dissipation Rate of Magnetic Field in Striped Pulsar Winds by the Effect of Turbulence

Formation of Turbulent and Magnetized Molecular Clouds via Accretion Flows of HI Clouds

Structure of dynamical condensation fronts in the interstellar medium

Toward Understanding the Origin of Turbulence in Molecular Clouds: Small Scale Structures as Units of Dynamical Multi-Phase Interstellar Medium

A Fast and Accurate Calculation Scheme for Ionization Degrees in Protoplanetary and Circumplanetary Disks with Charged Dust Grains

A Fast Numerical Scheme for Causal Relativistic Hydrodynamics with Dissipation

Effect of Magnetic Braking on the Circumstellar Disk Formation in a Strongly Magnetized Cloud

Gravitational Fragmentation of Expanding Shells. I. Linear Analysis

Gravitational Fragmentation of Expanding Shells. II. Three-dimensional Simulations

Recurrent Planet Formation and Intermittent Protostellar Outflows Induced by Episodic Mass Accretion

Secular Gravitational Instability of a Dust Layer in Shear Turbulence

Smoothed particle magnetohydrodynamics with a Riemann solver and the method of characteristics

Toward Understanding the Cosmic-ray Acceleration at Young Supernova Remnants Interacting with Interstellar Clouds: Possible Applications to RX J1713.7-3946

Analysis of the Relaxation Process using Non-Relativistic Kinetic Equation

Emergence of Protoplanetary Disks and Successive Formation of Gaseous Planets by Gravitational Instability

Formation Process of the Circumstellar Disk: Long-term Simulations in the Main Accretion Phase of Star Formation

Gas Accretion onto a Protoplanet and Formation of a Gas Giant Planet

N-Body Simulation of Planetesimal Formation through Gravitational Instability of a Dust Layer in Laminar Gas Disk

Protoplanetary Disk Winds by Magnetorotational Instability : Formation of an Inner Hole and a Crucial Assist for Planet Formation

The relativistic kinetic dispersion relation: Comparison of the relativistic Bhatnagar-Gross-Krook model and Grad's 14-moment expansion

Two-Dimensional Study of the Propagation of Planetary Wake and the Indication to Gap Opening in an Inviscid Protoplanetary Disk

Two-Step Acceleration Model of Cosmic Rays at Middle-Aged SNR

Binary Formation with Different Metallicities: Dependence on Initial Conditions

First Direct Simulation of Brown Dwarf Formation in a Compact Cloud Core

Kelvin-Helmholtz instabilities with Godunov SPH

N-Body Simulation of Planetesimal Formation through Gravitational Instability and Coagulation. II. Accretion Model

Orbital Evolution of a Particle Interacting with a Single Planet in a Protoplanetary Disk

Protostellar Jet and Outflow in the Collapsing Cloud Core